Coconut genome assembly enables evolutionary analysis of palms and highlights signaling pathways involved in salt tolerance

Coconut (Cocos nucifera) is the emblematic palm of tropical coastal areas all around the globe. It provides vital resources to millions of farmers. In an effort to better understand its evolutionary history and to develop genomic tools for its improvement, a sequence draft was recently released. Here, we present a dense linkage map (8402 SNPs) aiming to assemble the large genome of coconut (2.42 Gbp, 2n = 32) into 16 pseudomolecules. As a result, 47% of the sequences (representing 77% of the genes) were assigned to 16 linkage groups and ordered. We observed segregation distortion in chromosome Cn15, which is a signature of strong selection among pollen grains, favouring the maternal allele. Comparing our results with the genome of the oil palm Elaeis guineensis allowed us to identify major events in the evolutionary history of palms. We find that coconut underwent a massive transposable element invasion in the last million years, which could be related to the fluctuations of sea level during the glaciations at Pleistocene that would have triggered a population bottleneck. Finally, to better understand the facultative halophyte trait of coconut, we conducted an RNA-seq experiment on leaves to identify key players of signaling pathways involved in salt stress response. Altogether, our findings represent a valuable resource for the coconut breeding community.

The genome draft of coconut (Cocos nucifera)

Coconut palm (Cocos nucifera,2n = 32), a member of genus Cocos and family Arecaceae (Palmaceae), is an important tropical fruit and oil crop. Currently, coconut palm is cultivated in 93 countries, including Central and South America, East and West Africa, Southeast Asia and the Pacific Islands, with a total growth area of more than 12 million hectares [1]. Coconut palm is generally classified into 2 main categories: “Tall” (flowering 8–10 years after planting) and “Dwarf” (flowering 4–6 years after planting), based on morphological characteristics and breeding habits. This Palmae species has a long growth period before reproductive years, which hinders conventional breeding progress. In spite of initial successes, improvements made by conventional breeding have been very slow. In the present study, we obtained de novo sequences of the Cocos nucifera genome: a major genomic resource that could be used to facilitate molecular breeding in Cocos nucifera and accelerate the breeding process in this important crop. A total of 419.67 gigabases (Gb) of raw reads were generated by the Illumina HiSeq 2000 platform using a series of paired-end and mate-pair libraries, covering the predicted Cocos nucifera genome length (2.42 Gb, variety “Hainan Tall”) to an estimated ×173.32 read depth. A total scaffold length of 2.20 Gb was generated (N50 = 418 Kb), representing 90.91% of the genome. The coconut genome was predicted to harbor 28 039 protein-coding genes, which is less than in Phoenix dactylifera (PDK30: 28 889), Phoenix dactylifera (DPV01: 41 660), and Elaeis guineensis (EG5: 34 802). BUSCO evaluation demonstrated that the obtained scaffold sequences covered 90.8% of the coconut genome and that the genome annotation was 74.1% complete. Genome annotation results revealed that 72.75% of the coconut genome consisted of transposable elements, of which long-terminal repeat retrotransposons elements (LTRs) accounted for the largest proportion (92.23%). Comparative analysis of the antiporter gene family and ion channel gene families between C. nucifera and Arabidopsis thaliana indicated that significant gene expansion may have occurred in the coconut involving Na⁺/H⁺ antiporter, carnitine/acylcarnitine translocase, potassium-dependent sodium-calcium exchanger, and potassium channel genes. Despite its agronomic importance, C. nucifera is still under-studied. In this report, we present a draft genome of C. nucifera and provide genomic information that will facilitate future functional genomics and molecular-assisted breeding in this crop species.

Improving transcriptome de novo assembly by using a reference genome of a related species: Translational genomics from oil palm to coconut

The palms are a family of tropical origin and one of the main constituents of the ecosystems of these regions around the world. The two main species of palm represent different challenges: coconut (Cocos nucifera L.) is a source of multiple goods and services in tropical communities, while oil palm (Elaeis guineensis Jacq) is the main protagonist of the oil market. In this study, we present a workflow that exploits the comparative genomics between a target species (coconut) and a reference species (oil palm) to improve the transcriptomic data, providing a proteome useful to answer functional or evolutionary questions. This workflow reduces redundancy and fragmentation, two inherent problems of transcriptomic data, while preserving the functional representation of the target species. Our approach was validated in Arabidopsis thaliana using Arabidopsis lyrata and Capsella rubella as references species. This analysis showed the high sensitivity and specificity of our strategy, relatively independent of the reference proteome. The workflow increased the length of proteins products in A. thaliana by 13%, allowing, often, to recover 100% of the protein sequence length. In addition redundancy was reduced by a factor greater than 3. In coconut, the approach generated 29,366 proteins, 1,246 of these proteins deriving from new contigs obtained with the BRANCH software. The coconut proteome presented a functional profile similar to that observed in rice and an important number of metabolic pathways related to secondary metabolism. The new sequences found with BRANCH software were enriched in functions related to biotic stress. Our strategy can be used as a complementary step to de novo transcriptome assembly to get a representative proteome of a target species. The results of the current analysis are available on the website PalmComparomics (http://palm-comparomics.southgreen.fr/).

Ploidy and domestication are associated with genome size variation in Palms

PREMISE OF THE STUDY

The genome size of a species (C-value) is associated with growth, development and adaptation to environmental changes. Angiosperm C-values range 1200-fold and frequently vary within species, although little is known about the impacts of domestication on genome size. Genome size variation among related species of palms is of evolutionary significance because changes characterize clades and may be associated with polyploidy, transposon amplifications, deletions, or rearrangements. Further knowledge of genome size will provide crucial information needed for planning of whole genome sequencing and accurate annotations. We studied the genome size of Cocos nucifera and its variation among cultivars, and compared it to values for related palms from the Attaleinae subtribe.

METHODS

Flow cytometric analysis of isolated nuclei from young palm leaves was used to estimate genome sizes of 23 coconut cultivars (Talls, Dwarfs, and hybrids) worldwide and 17 Cocoseae species. Ancestral genome size was reconstructed on a maximum likelihood phylogeny of Attaleinae from seven WRKY loci.

KEY RESULTS

The coconut genome is large-averaging 5.966 pg-and shows intraspecific variation associated with domestication. Variation among Tall coconuts was significantly greater than among Dwarfs. Attaleinae genomes showed moderate size variation across genera, except polyploids Jubaeopsis caffra, Voanioala gerardii, Beccariophoenix alfredii, and Allagoptera caudescens, which had larger genomes.

CONCLUSIONS

Our results contribute to the understanding of the relationship between domestication and genome size in long-lived tree crops and provide a basis for whole-genome sequencing of the coconut and other domesticated plants. Polyploidy evolved independently in two clades within Attaleinae.

Zika virus infection complicated by Guillain-Barre syndrome--case report, French Polynesia, December 2013

Zika fever, considered as an emerging disease of arboviral origin, because of its expanding geographic area, is known as a benign infection usually presenting as an influenza-like illness with cutaneous rash. So far, Zika virus infection has never led to hospitalisation. We describe the first case of Guillain-Barré syndrome (GBS) occurring immediately after a Zika virus infection, during the current Zika and type 1 and 3 dengue fever co-epidemics in French Polynesia.

The presence of coconut in southern Panama in pre-Columbian times: clearing up the confusion

BACKGROUND

The pre-Columbian presence of coconut on the Pacific coast of Panama is attested by a number of independent written accounts. However, recent papers question their accuracy and conclude that coconut was introduced to the region by the Spaniards after their conquests. Scope In order to examine the value of such claims, an extensive search was conducted of the relevant historical accounts of coconut in America and in the Orient.

KEY RESULTS

The Spanish chronicler Oviedo (1478-1557) is found to have effectively used fruit and seed size to distinguish coconut from other palms. In addition, it is shown that he has been inaccurately faulted with incorrectly representing a cluster of coconuts. The original drawing, a cluster of a native Bactris, was in the marginalia and was only assigned to coconut after Oviedo's death. Finally, the location is identified of a coastal Panamanian site described by Pedro Mártir de Anglería and where tidal dispersal of coconuts was observed.

CONCLUSIONS

This previously overlooked evidence confirms the pre-historical presence of coconut in Panama. Genetic data indicate that it must have been brought there directly or indirectly from the Philippines. But when, where and by whom remains a subject of research. Further molecular marker studies, computer simulation of natural drift and archaeological research could contribute to this research.

Independent origins of cultivated coconut (Cocos nucifera L.) in the old world tropics

As a portable source of food, water, fuel, and construction materials, the coconut (Cocos nucifera L.) played a fundamental role in human migrations and the development of civilization across the humid tropics. Here we investigated the coconut's domestication history and its population genetic structure as it relates to human dispersal patterns. A sample of 1,322 coconut accessions, representing the geographical and phenotypic diversity of the species, was examined using ten microsatellite loci. Bayesian analyses reveal two highly genetically differentiated subpopulations that correspond to the Pacific and Indo-Atlantic oceanic basins. This pattern suggests independent origins of coconut cultivation in these two world regions, with persistent population structure on a global scale despite long-term human cultivation and dispersal. Pacific coconuts show additional genetic substructure corresponding to phenotypic and geographical subgroups; moreover, the traits that are most clearly associated with selection under human cultivation (dwarf habit, self-pollination, and “niu vai” fruit morphology) arose only in the Pacific. Coconuts that show evidence of genetic admixture between the Pacific and Indo-Atlantic groups occur primarily in the southwestern Indian Ocean. This pattern is consistent with human introductions of Pacific coconuts along the ancient Austronesian trade route connecting Madagascar to Southeast Asia. Admixture in coastal east Africa may also reflect later historic Arab trading along the Indian Ocean coastline. We propose two geographical origins of coconut cultivation: island Southeast Asia and southern margins of the Indian subcontinent.

Population structures of Brazilian tall coconut (Cocos nucifera L.) by microsatellite markers

Coconut palms of the Tall group were introduced to Brazil from the Cape Verde Islands in 1553. The present study sought to evaluate the genetic diversity among and within Brazilian Tall coconut populations. Samples were collected of 195 trees from 10 populations. Genetic diversity was accessed by investigating 13 simple sequence repeats (SSR) loci. This provided a total of 68 alleles, ranging from 2 to 13 alleles per locus, with an average of 5.23. The mean values of gene diversity (H_e ) and observed heterozygosity (H_o ) were 0.459 and 0.443, respectively. The genetic differentiation among populations was estimated at θ^P=0.1600and the estimated apparent outcrossing rate was t_a = 0.92. Estimates of genetic distances between the populations varied from 0.034 to 0.390. Genetic distance and the corresponding clustering analysis indicate the formation of two groups. The first consists of the Baía Formosa, Georgino Avelino, and São José do Mipibu populations and the second consists of the Japoatã, Pacatuba, and Praia do Forte populations. The correlation matrix between genetic and geographic distances was positive and significant at a 1% probability. Taken together, our results suggest a spatial structuring of the genetic variability among the populations. Geographically closer populations exhibited greater similarities.

Identification of molecular markers associated with mite resistance in coconut (Cocos nucifera L.)

Coconut mite ( Aceria guerreronis ‘Keifer’) has become a major threat to Indian coconut ( Coçcos nucifera L.) cultivators and the processing industry. Chemical and biological control measures have proved to be costly, ineffective, and ecologically undesirable. Planting mite-resistant coconut cultivars is the most effective method of preventing yield loss and should form a major component of any integrated pest management stratagem. Coconut genotypes, and mite-resistant and -susceptible accessions were collected from different parts of South India. Thirty-two simple sequence repeat (SSR) and 7 RAPD primers were used for molecular analyses. In single-marker analysis, 9 SSR and 4 RAPD markers associated with mite resistance were identified. In stepwise multiple regression analysis of SSRs, a combination of 6 markers showed 100% association with mite infestation. Stepwise multiple regression analysis for RAPD data revealed that a combination of 3 markers accounted for 83.86% of mite resistance in the selected materials. Combined stepwise multiple regression analysis of RAPD and SSR data showed that a combination of 5 markers explained 100% of the association with mite resistance in coconut. Markers associated with mite resistance are important in coconut breeding programs and will facilitate the selection of mite-resistant plants at an early stage as well as mother plants for breeding programs.

Comparison of the anaesthetic requirement with target-controlled infusion of propofol to insert the laryngeal tube vs. the laryngeal mask

BACKGROUND AND OBJECTIVE

The target effect-site concentration of propofol to insert a laryngeal mask airway was recently reported as almost 5 microg mL(-1). The present study aimed to determine the target effect-site concentration with target-controlled infusion of propofol to place classical larnygeal mask airway or current laryngeal tube in adult patients.

METHODS

We included 40 patients scheduled for short gynaecological and radiological procedures under general anaesthesia in a randomized, double-blind manner using the Dixon's up-and-down statistical method. Monitoring included standard cardiorespiratory monitors, and bispectral index monitoring was used for all patients. Anaesthesia was conducted with a target-controlled infusion system: Diprifusor. The initial target plasma concentration of propofol was 5 microg mL(-1), and was changed stepwise by 0.5 microg mL(-1) increments according to Dixon's up-and-down method. Criteria for acceptable insertion were: Muzi's score < or = 2, and mean arterial blood pressure, heart rate or bispectral index variation <20% the baseline values.

RESULTS

Target effect-site concentration of propofol required to insert laryngeal tube was 6.3 +/- 0.3 microg mL(-1) with Dixon method and ED50 was 6.1 microg mL(-1) (5.9-6.4) with logistic regression method. In the case of larnygeal mask airway they were 7.3 +/- 0.2 microg mL(-1) (Dixon method) and 7.3 microg mL(-1) (7.1-7.5; with logistic regression) respectively (P < 0.05). ED95 (logistic regression) was 6.8 microg mL(-1) (5.9-7.6) for laryngeal tube and 7.7 microg mL(-1) (7.3-8.0) for larnygeal mask airway (P < 0.05). Haemodynamic incidents were 55% in the larnygeal mask airway group vs. 30% in the laryngeal tube group (P < 0.05).

CONCLUSIONS

The target effect-site concentration of propofol for insertion of laryngeal tube was lower than for larnygeal mask airway (P < 0.05), with a consequent reduction of the propofol induced haemodynamic side-effects.

QTL analysis of fruit components in the progeny of a Rennell Island Tall coconut (Cocos nucifera L.) individual

We investigated the genetic factors controlling fruit components in coconut by performing QTL analyses for fruit component weights and ratios in a segregating progeny of a Rennell Island Tall genotype. The underlying linkage map of this population was already established in a previous study, as well as QTL analyses for fruit production, which were used to complement our results. The addition of 53 new markers (mainly SSRs) led to minor amendments in the map. A total of 52 putative QTLs were identified for the 11 traits under study. Thirty-four of them were grouped in six small clusters, which probably correspond to single pleiotropic genes. Some additional QTLs located apart from these clusters also had relatively large effects on the individual traits. The QTLs for fruit component weight, endosperm humidity and fruit production were found at different locations in the genome, suggesting that efficient marker-assisted selection for yield can be achieved by selecting QTLs for the individual components. The detected QTLs descend from a genotype belonging to the "Pacific" coconut group. Based on the known molecular and phenotypic differences between "Pacific" and "Indo-Atlantic" coconuts, we suggest that a large fraction of coconut genetic diversity is still to be investigated by studying populations derived from crosses between these groups.

GENECLASS2: a software for genetic assignment and first-generation migrant detection

GENECLASS2 is a software that computes various genetic assignment criteria to assign or exclude reference populations as the origin of diploid or haploid individuals, as well as of groups of individuals, on the basis of multilocus genotype data. In addition to traditional assignment aims, the program allows the specific task of first-generation migrant detection. It includes several Monte Carlo resampling algorithms that compute for each individual its probability of belonging to each reference population or to be a resident (i.e., not a first-generation migrant) in the population where it was sampled. A user-friendly interface facilitates the treatment of large datasets.

Analytical bayesian approach for assigning individuals to populations

We propose a general formulation of the Bayesian method for assigning individuals to a population among a predetermined set of reference populations using molecular marker information. Compared to previously published methods, ours allows us to consider different types of prior information about allele frequencies by using a Dirichlet prior probability distribution. It also makes it possible to assign a set of individuals assumed to belong to the same population with increased accuracy using their pooled genotype data. The efficiency of the method is illustrated by application to a group of closely related coconut populations. An interesting feature of the Bayesian procedure is the way it handles imprecise information. With a poor or even incomplete dataset, assignment is still be possible and gives valid results: poor data quality is reflected in an ambiguous result rather than in a false conclusion.

Construction of a linkage map of the Rennell Island Tall coconut type (Cocos nuciferaL.) and QTL analysis for yield characters

AFLP and SSR DNA markers were used to construct a linkage map in the coconut (Cocos nucifera L.; 2n = 32) type Rennell Island Tall (RIT). A total of 227 markers were arranged into 16 linkage groups. The total genome length corresponded to 1971 cM for the RIT map, with 5–23 markers per linkage group. QTL analysis for yield characters in two consecutive sampling periods identified nine loci. Three and two QTLs were detected for number of bunches and one and three QTLs for number of nuts. The correlation of trait values between characters and evaluation periods is partially reflected in identical QTLs. The QTLs represent characters that are important in coconut breeding. The cosegregation of markers with these QTLs provides an opportunity for marker-assisted selection in coconut breeding programmes.Key words: coconut, QTL, AFLP, SSR, marker-assisted selection (MAS).

Construction of a linkage map of the Rennell Island Tall coconut type (Cocos nucifera L.) and QTL analysis for yield characters

AFLP and SSR DNA markers were used to construct a linkage map in the coconut (Cocos nucifera L.; 2n = 32) type Rennell Island Tall (RIT). A total of 227 markers were arranged into 16 linkage groups. The total genome length corresponded to 1971 cM for the RIT map, with 5-23 markers per linkage group. QTL analysis for yield characters in two consecutive sampling periods identified nine loci. Three and two QTLs were detected for number of bunches and one and three QTLs for number of nuts. The correlation of trait values between characters and evaluation periods is partially reflected in identical QTLs. The QTLs represent characters that are important in coconut breeding. The cosegregation of markers with these QTLs provides an opportunity for marker-assisted selection in coconut breeding programmes.

Interactions of DNA with a new electron-deficient tentacle porphyrin: meso-tetrakis[2,3,5,6-tetrafluoro-4-(2-t rimethylammoniumethyl-amine)phenyl]porphy rin

A new electron-deficient tentacle porphyrin meso-tetrakis[2,3,5,6-tetrafluoro-4-(2-trimethylammoniumethylamine )phenyl]porphyrin (TthetaF4TAP) has been synthesized. The binding interactions of TthetaF4TAP with DNA polymers were studied for comparison to those of an electron-deficient tentacle porphyrin and an electron-rich tentacle porphyrin; these previously studied porphyrins bind to DNA primarily by intercalative and outside-binding modes, respectively. The three tentacle porphyrins have similar size and shape. The basicity of TthetaF4TAP indicated that it has electronic characteristics similar to those of the intercalating electron-deficient tentacle porphyrin. However, TthetaF4TAP binds to calf thymus DNA, [poly(dA-dT)]2, and [poly(dG-dC)]2 in a self-stacking, outside-binding manner under all conditions. Evidence for this binding mode included a significant hypochromicity of the Soret band, a conservative induced CD spectrum, and the absence of an increase in DNA solution viscosity. As found previously for the electron-rich porphyrin, the results suggest that combinations of closely related self-stacked forms coexist. The mix of forms depended on the DNA and the solution conditions. There are probably differences in the detailed features of the self-stacking adducts for the two types of tentacle porphyrins, especially at high R (ratio of porphyrin to DNA). At low R values, the induced CD signal of TthetaF4TAP/CT DNA resembled that of TthetaF4TAP/[poly(dA-dT)]2, suggesting that TthetaF4TAP binds preferentially at AT regions. Competitive binding experiments gave evidence that TthetaF4TAP binds preferentially to [poly(dA-dT)]2 over [poly (dG-dC)]2. Thus, despite the long, positively charged, flexible substituents on the porphyrin, the binding of TthetaF4TAP is significantly affected by base-pair composition. Similar characteristics were found previously for the electron-rich tentacle porphyrin. Thus, significant changes in electron richness have relatively minor effects on this outside binding selectivity for AT regions. TthetaF4TAP is the first porphyrin with electron deficiency and shape similar to intercalating porphyrins that does not appear to intercalate. All porphyrins reported to intercalate have had pyridinium substituents. Thus, the electronic distribution in the porphyrin ring, not just the overall electron richness, may play a role in facilitating intercalation.

Analysis of the genetic effects for several traits in oil palm (Elaeis guineensis Jacq.) populations. I. Population means

An oil palm experiment was set up in the Ivory Coast to compare the effects of crossing and selfing within two origins, Deli and La Mé, on the mean and the variability of Deli x Mé between-origin hybrids. The originality of the experiment lay in the crossing plan, which provided access to genetic parameters related to additivity, dominance and different components of epistasis. This first part covers the analysis of the components of the mean. The parents used were obtained from four palms, two from each origin. Those of La Mé origin were half-sibs. The common parent came from a wild stand in the Ivory Coast. Those of Deli origin were from two different populations bred in Southeast Asia for several generations from a narrow genetic base. These four parents gave rise to nine Deli x La Mé hybrid populations with double-cousintype links. The additive component is more important within the Deli origin than within the La Mé origin. This may be explained by the large genetic divergence between the two Deli parent palms. On the other hand, the additive(*) additive epistasis is more substantial within the La Mé origin, probably because of inbreeding. The discussion concentrates on how this information should be used when choosing parents to be crossed and tested and to produce improved populations. The crossing plan proposed can be of general use and is suitable for other species in a reciprocal recurrent selection programme.