An improved assembly of the pearl millet reference genome using Oxford Nanopore long reads and optical mapping

Pearl millet (Pennisetum glaucum (L.)) R. Br. syn. Cenchrus americanus (L.) Morrone) is an important crop in South Asia and sub-Saharan Africa which contributes to ensuring food security. Its genome has an estimated size of 1.76 Gb and displays a high level of repetitiveness above 80%. A first assembly was previously obtained for the Tift 23D2B1-P1-P5 cultivar genotype using short-read sequencing technologies. This assembly is, however, incomplete and fragmented with around 200 Mb unplaced on chromosomes. We report here an improved quality assembly of the pearl millet Tift 23D2B1-P1-P5 cultivar genotype obtained with an approach combining Oxford Nanopore long reads and Bionano Genomics optical maps. This strategy allowed us to add around 200 Mb at the chromosome-level assembly. Moreover we strongly improved continuity in the order of the contigs and scaffolds within the chromosomes, particularly in the centromeric regions. Notably, we added more than 100 Mb around the centromeric region on chromosome 7. This new assembly also displayed a higher gene completeness with a complete BUSCO score of 98.4% using the Poales database. This more complete and higher quality assembly of the Tift 23D2B1-P1-P5 genotype now available to the community will help in the development of research on the role of structural variants, and more broadly in genomics studies and the breeding of pearl millet.

Genomic footprints of selection in early-and late-flowering pearl millet landraces

Pearl millet is among the top three-cereal production in one of the most climate vulnerable regions, sub-Saharan Africa. Its Sahelian origin makes it adapted to grow in poor sandy soils under low soil water regimes. Pearl millet is thus considered today as one of the most interesting crops to face the global warming. Flowering time, a trait highly correlated with latitude, is one of the key traits that could be modulated to face future global changes. West African pearl millet landraces, can be grouped into early- (EF) and late-flowering (LF) varieties, each flowering group playing a specific role in the functioning and resilience of Sahelian smallholders. The aim of this study was thus to detect genes linked to flowering but also linked to relevant traits within each flowering group. We thus investigated genomic and phenotypic diversity in 109 pearl millet landrace accessions, i.e., 66 early-flowering and 43 late-flowering, grown in the groundnut basin, the first area of rainfed agriculture in Senegal dominated by dry cereals (millet, maize, and sorghum) and legumes (groundnuts, cowpeas). We were able to confirm the role of PhyC gene in pearl millet flowering and identify several other genes that appear to be as much as important, such as FSR12 and HAC1. HAC1 and two other genes appear to be part of QTLs previously identified and deserve further investigation. At the same time, we were able to highlight a several genes and variants that could contribute to the improvement of pearl millet yield, especially since their impact was demonstrated across flowering cycles.

Unveiling biogeographical patterns of the ichthyofauna in the Tuichi basin, a biodiversity hotspot in the Bolivian Amazon, using environmental DNA

To date, more than 2400 valid fish species have been recorded in the Amazon basin. However, some regions remain poorly documented. This is the case in the Beni basin and in particular in one of its main sub-basins, the Tuichi, an Andean foothills rivers flowing through the Madidi National Park in the Bolivian Amazonia. The knowledge of its ichthyological diversity is, however, essential for the management and protection of aquatic ecosystems, which are threatened by the development of infrastructures (dams, factories and cities), mining and deforestation. Environmental DNA (eDNA) has been relatively little used so far in the Amazon basin. We sampled eDNA from water in 34 sites in lakes and rivers in the Beni basin including 22 sites in the Tuichi sub-basin, during the dry season. To assess the biogeographical patterns of the amazonian ichthyofauna, we implemented a metabarcoding approach using two pairs of specific primers designed and developed in our laboratory to amplify two partially overlapping CO1 fragments, one of 185bp and another of 285bp. We detected 252 fish taxa (207 at species level) among which 57 are newly identified for the Beni watershed. Species compositions are significantly different between lakes and rivers but also between rivers according to their hydrographic rank and altitude. Furthermore, the diversity patterns are related to the different hydro-ecoregions through which the Tuichi flows. The eDNA approach makes it possible to identify and complete the inventory of the ichthyofauna in this still poorly documented Amazon basin. However, taxonomic identification remains constrained by the lack of reference barcodes in public databases and does not allow the assignment of all OTUs. Our results can be taken into account in conservation and management strategies and could serve as a baseline for future studies, including on other Andean tributaries.

Species-level ichthyoplankton dynamics for 97 fishes in two major river basins of the Amazon using quantitative metabarcoding

The Amazon basin holds the world's largest freshwater fish diversity. Information on the intensity and timing of reproductive ecology of Amazonian fish is scant. We use a metabarcoding method by capture using a single probe to quantify species-level ichthyoplankton dynamics. We sampled the Marañón and the Ucayali rivers in Peru monthly for 2 years. We identified 97 species that spawned mainly during the flood start, the flood end or the receding periods, although some species had spawning activity in more than one period. This information was new for 40 of the species in the Amazon basin and 80 species in Peru. Most species ceased spawning for a month during a strong hydrological anomaly in January 2016, demonstrating the rapidity with which they react to environmental modifications during the breeding season. We also document another unreported event in the Amazon basin, the inverse phenology of species belonging to one genus (Triportheus). Overall larval flow in the Marañón was more than twice that of the Ucayali, including for most commercial species (between two and 20 times higher), whereas the Ucayali accounts for ~80% of the fisheries landings in the region. Our results are discussed in the light of the main anthropogenic threats to fishes, hydropower dam construction and the Hidrovía Amazónica, and should serve as a pre-impact baseline.

Transferability, development of simple sequence repeat (SSR) markers and application to the analysis of genetic diversity and population structure of the African fan palm (Borassus aethiopum Mart.) in Benin

Background

In Sub-Saharan Africa, Borassus aethiopum Mart. (African fan palm) is an important non-timber forest product-providing palm that faces multiple anthropogenic threats to its genetic diversity. However, this species is so far under-studied, which prevents its sustainable development as a resource. The present work is a first attempt at characterizing the genetic diversity and population structure of B. aethiopum across nine collection sites spanning the three climatic regions of Benin, West Africa, through the use of microsatellite markers.

Results

During a first phase we relied on the reported transferability of primers developed in other palm species. We find that, in disagreement with previously published results, only 22.5% of the markers tested enable amplification of B. aethiopum DNA and polymorphism detection is very low.

In a second phase, we generated a B. aethiopum-specific genomic dataset through high-throughput sequencing and used it for the de novo detection of microsatellite loci. Among the primer pairs targeting these, 11 detected polymorphisms and were further used for analyzing genetic diversity. Across the nine sites, expected heterozygosity (He) ranges from 0.263 to 0.451 with an overall average of 0.354, showing a low genetic diversity. Analysis of molecular variance (AMOVA) shows that within-site variation accounts for 53% of the genetic variation. Accordingly, the low number of migrants and positive values of the fixation index (F) in sites from both the Central (Sudano-Guinean) and the Southern (Guinean) climatic regions suggest limited gene flow between sites. The global correlation between genetic and geographic distances is weak; however, our clustering analyses indicate that B. aethiopum palms from Savè (Center) are genetically more similar to those from the North than to samples from other Central sites.

Conclusions

In the light of our results, we discuss the use of inter-species transfer vs. de novo development of microsatellite markers in genetic diversity analyses targeting under-studied species, and suggest future applications for our molecular resources. We propose that, while prominent short-range pollen and seed dispersal in Benin explain most of our results, gene flux between the Central and Northern regions, as a result of animal and/or human migrations, might underlie the Savè discrepancy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00955-y.

Pearl millet genomic vulnerability to climate change in West Africa highlights the need for regional collaboration

Climate change is already affecting agro-ecosystems and threatening food security by reducing crop productivity and increasing harvest uncertainty. Mobilizing crop diversity could be an efficient way to mitigate its impact. We test this hypothesis in pearl millet, a nutritious staple cereal cultivated in arid and low-fertility soils in sub-Saharan Africa. We analyze the genomic diversity of 173 landraces collected in West Africa together with an extensive climate dataset composed of metrics of agronomic importance. Mapping the pearl millet genomic vulnerability at the 2050 horizon based on the current genomic-climate relationships, we identify the northern edge of the current areas of cultivation of both early and late flowering varieties as being the most vulnerable to climate change. We predict that the most vulnerable areas will benefit from using landraces that already grow in equivalent climate conditions today. However, such seed-exchange scenarios will require long distance and trans-frontier assisted migrations. Leveraging genetic diversity as a climate mitigation strategy in West Africa will thus require regional collaboration.

Replacement of local crops with alternative varieties adapted to future conditions may improve food security under climate change. Here the authors apply landscape genomics and ensemble climate modelling to pearl millet in West Africa, supporting the potential of transfrontier assisted seed exchange.

Aquaporins are main contributors to root hydraulic conductivity in pearl millet [Pennisetum glaucum (L) R. Br.]

Pearl millet is a key cereal for food security in arid and semi-arid regions but its yield is increasingly threatened by water stress. Physiological mechanisms relating to conservation of soil water or increased water use efficiency can alleviate that stress. Aquaporins (AQP) are water channels that mediate root water transport, thereby influencing plant hydraulics, transpiration and soil water conservation. However, AQP remain largely uncharacterized in pearl millet. Here, we studied AQP function in root water transport in two pearl millet lines contrasting for water use efficiency (WUE). We observed that these lines also contrasted for root hydraulic conductivity (Lpr) and AQP contribution to Lpr. The line with lower WUE showed significantly higher AQP contribution to Lpr. To investigate AQP isoforms contributing to Lpr, we developed genomic approaches to first identify the entire AQP family in pearl millet and secondly, characterize the plasma membrane intrinsic proteins (PIP) gene expression profile. We identified and annotated 33 AQP genes in pearl millet, among which ten encoded PIP isoforms. PgPIP1-3 and PgPIP1-4 were significantly more expressed in the line showing lower WUE, higher Lpr and higher AQP contribution to Lpr. Overall, our study suggests that the PIP1 AQP family are the main regulators of Lpr in pearl millet and may possibly be associated with mechanisms associated to whole plant water use. This study paves the way for further investigations on AQP functions in pearl millet hydraulics and adaptation to environmental stresses.

Abandonment of pearl millet cropping and homogenization of its diversity over a 40 year period in Senegal

Cultivated diversity is considered an insurance against major climatic variability. However, since the 1980s, several studies have shown that climate variability and agricultural changes may already have locally eroded crop genetic diversity. We studied pearl millet diversity in Senegal through a comparison of pearl millet landraces collected 40 years apart. We found that more than 20% of villages visited in 1976 had stopped growing pearl millet. Despite this, its overall genetic diversity has been maintained but differentiation between early- and late-flowering accessions has been reduced. We also found stronger crop-to-wild gene flow than wild-to-crop gene flow and that wild-to-crop gene flow was weaker in 2016 than in 1976. In conclusion, our results highlight genetic homogenization in Senegal. This homogenization within cultivated pearl millet and between wild and cultivated forms is a key factor in genetic erosion and it is often overlooked. Improved assessment and conservation strategies are needed to promote and conserve both wild and cultivated pearl millet diversity.

Microsatellite markers development for Indonesian nutmeg (Myristica fragrans Houtt.) and transferability to other Myristicaceae spp

Myristica fragrans (Myristicaceae) is a tropical evergreen tree that yields the two famous spices: nutmeg and mace. Despite its socio-economic importance, the spatial distribution of its genetic diversity is barely documented. In this aim, 48 nuclear microsatellite markers were isolated of which 14 were polymorphic in M. fragrans. Number of alleles per locus ranged from 2 to 6. The level of observed heterozygosity ranged from 0.038 to 0.929 across loci. Transferability of these microsatellites in other Myristica species (M. fatua, M. argentea, and M. crassipes) and Myristicaceae species (Horsfieldia palauensis) was tested and successful. These new microsatellites will be useful for future investigation on genetic diversity and population structure of M. fragrans and phylogenetically-related species.

A single amino acid substitution (H451Y) in Leishmania calcium-dependent kinase SCAMK confers high tolerance and resistance to antimony

Background

For almost a century, antimonials have remained the first-line drugs for the treatment of leishmaniasis. However, little is known about their mode of action and clinical resistance mechanisms.

Objectives

We have previously shown that Leishmania nicotinamidase (PNC1) is an essential enzyme for parasite NAD+ homeostasis and virulence in vivo. Here, we found that parasites lacking the pnc1 gene (Δpnc1) are hypersusceptible to the active form of antimony (SbIII) and used these mutant parasites to better understand antimony’s mode of action and the mechanisms leading to resistance.

Methods

SbIII-resistant WT and Δpnc1 parasites were selected in vitro by a stepwise selection method. NAD(H)/NADP(H) dosages and quantitative RT–PCR experiments were performed to explain the susceptibility differences observed between strains. WGS and a marker-free CRISPR/Cas9 base-editing approach were used to identify and validate the role of a new resistance mutation.

Results

NAD+-depleted Δpnc1 parasites were highly susceptible to SbIII and this phenotype could be rescued by NAD+ precursor or trypanothione precursor supplementation. Δpnc1 parasites could become resistant to SbIII by an unknown mechanism. WGS revealed a unique amino acid substitution (H451Y) in an EF-hand domain of an orphan calcium-dependent kinase, recently named SCAMK. When introduced into a WT reference strain by base editing, the H451Y mutation allowed Leishmania parasites to survive at extreme concentrations of SbIII, potentiating the rapid emergence of resistant parasites.

Conclusions

These results establish that Leishmania SCAMK is a new central hub of antimony’s mode of action and resistance development, and uncover the importance of drug tolerance mutations in the evolution of parasite drug resistance.

Yam genomics supports West Africa as a major cradle of crop domestication

While there has been progress in our understanding of the origin and history of agriculture in sub-Saharan Africa, a unified perspective is still lacking on where and how major crops were domesticated in the region. Here, we investigated the domestication of African yam (Dioscorea rotundata), a key crop in early African agriculture. Using whole-genome resequencing and statistical models, we show that cultivated yam was domesticated from a forest species. We infer that the expansion of African yam agriculture started in the Niger River basin. This result, alongside with the origins of African rice and pearl millet, supports the hypothesis that the vicinity of the Niger River was a major cradle of African agriculture.

Long-fragment targeted capture for long-read sequencing of plastomes

PREMISE

Third-generation sequencing methods generate significantly longer reads than those produced using alternative sequencing methods. This provides increased possibilities for the study of biodiversity, phylogeography, and population genetics. We developed a protocol for in-solution enrichment hybridization capture of long DNA fragments applicable to complete plastid genomes.

METHODS AND RESULTS

The protocol uses cost-effective in-house probes developed via long-range PCR and was used in six non-model monocot species (Poaceae: African rice, pearl millet, fonio; and three palm species). DNA was extracted from fresh and silica gel-dried leaves. Our protocol successfully captured long-read plastome fragments (3151 bp median on average), with an enrichment rate ranging from 15% to 98%. DNA extracted from silica gel-dried leaves led to low-quality plastome assemblies when compared to DNA extracted from fresh tissue.

CONCLUSIONS

Our protocol could also be generalized to capture long sequences from specific nuclear fragments.

Metabarcoding by capture using a single COI probe (MCSP) to identify and quantify fish species in ichthyoplankton swarms

The ability to determine the composition and relative frequencies of fish species in large ichthyoplankton swarms could have extremely important ecological applications However, this task is currently hampered by methodological limitations. We proposed a new method for Amazonian species based on hybridization capture of the COI gene DNA from a distant species (Danio rerio), absent from our study area (the Amazon basin). The COI sequence of this species is approximately equidistant from all COI of Amazonian species available. By using this sequence as probe we successfully facilitated the simultaneous identification of fish larvae belonging to the order Siluriformes and to the Characiformes represented in our ichthyoplankton samples. Species relative frequencies, estimated by the number of reads, showed almost perfect correlations with true frequencies estimated by a Sanger approach, allowing the development of a quantitative approach. We also proposed a further improvement to a previous protocol, which enables lowering the sequencing effort by 40 times. This new Metabarcoding by Capture using a Single Probe (MCSP) methodology could have important implications for ecology, fisheries management and conservation in fish biodiversity hotspots worldwide. Our approach could easily be extended to other plant and animal taxa.

Identification of a Hypervirulent Pathotype of Rice yellow mottle virus: A Threat to Genetic Resistance Deployment in West-Central Africa

Rice yellow mottle virus (RYMV) causes high losses to rice production in Africa. Several sources of varietal high resistance are available but the emergence of virulent pathotypes that are able to overcome one or two resistance alleles can sometimes occur. Both resistance spectra and viral adaptability have to be taken into account to develop sustainable rice breeding strategies against RYMV. In this study, we extended previous resistance spectrum analyses by testing the rymv1-4 and rymv1-5 alleles that are carried by the rice accessions Tog5438 and Tog5674, respectively, against isolates that are representative of RYMV genetic and pathogenic diversity. Our study revealed a hypervirulent pathotype, named thereafter pathotype T', that is able to overcome all known sources of high resistance. This pathotype, which is spatially localized in West-Central Africa, appears to be more abundant than previously suspected. To better understand the adaptive processes of pathotype T', molecular determinants of resistance breakdown were identified via Sanger sequencing and validated through directed mutagenesis of an infectious clone. These analyses confirmed the key role of convergent nonsynonymous substitutions in the central part of the viral genome-linked protein to overcome RYMV1-mediated resistance. In addition, deep-sequencing analyses revealed that resistance breakdown does not always coincide with fixed mutations. Actually, virulence mutations that are present in a small proportion of the virus population can be sufficient for resistance breakdown. Considering the spatial distribution of RYMV strains in Africa and their ability to overcome the RYMV resistance genes and alleles, we established a resistance-breaking risk map to optimize strategies for the deployment of sustainable and resistant rice lines in Africa.

No Excess of Cis-Regulatory Variation Associated with Intraspecific Selection in Wild Pearl Millet (Cenchrus americanus)

Several studies suggest that cis-regulatory mutations are the favorite target of evolutionary changes, one reason being that cis-regulatory mutations might have fewer deleterious pleiotropic effects than protein-coding mutations. A review of the process also suggests that this bias towards adaptive cis-regulatory variation might be less pronounced at the intraspecific level compared with the interspecific level. In this study, we assessed the contribution of cis-regulatory variation to adaptation at the intraspecific level using populations of wild pearl millet (Cenchrus americanus ssp. monodii) sampled along an environmental gradient in Niger. From RNA sequencing of hybrids to assess allele-specific expression, we identified genes with cis-regulatory divergence between two parental accessions collected in contrasted environmental conditions. This revealed that ∼15% of transcribed genes showed cis-regulatory variation. Intersecting the gene set exhibiting cis-regulatory variation with the gene set identified as targets of selection revealed no excess of cis-acting mutations among the selected genes. We additionally found no excess of cis-regulatory variation among genes associated with adaptive traits. As our approach relied on methods identifying mainly genes submitted to strong selection pressure or with high phenotypic effect, the contribution of cis-regulatory changes to soft selection or polygenic adaptive traits remains to be tested. However our results favor the hypothesis that enrichment of adaptive cis-regulatory divergence builds up over time. For short evolutionary time-scales, cis-acting mutations are not predominantly involved in adaptive evolution associated with strong selective signal.

Human management and hybridization shape treegourd fruits in the Brazilian Amazon Basin

Local people's perceptions of cultivated and wild agrobiodiversity, as well as their management of hybridization are still understudied in Amazonia. Here we analyze domesticated treegourd (Crescentia cujete), whose versatile fruits have technological, symbolic, and medicinal uses. A wild relative (C. amazonica) of the cultivated species grows spontaneously in Amazonian flooded forests. We demonstrated, using whole chloroplast sequences and nuclear microsatellites, that the two species are strongly differentiated. Nonetheless, they hybridize readily throughout Amazonia and the proportions of admixture correlate with fruit size variation of cultivated trees. New morphotypes arise from hybridization, which are recognized by people and named as local varieties. Small hybrid fruits are used to make the important symbolic rattle (maracá), suggesting that management of hybrid trees is an ancient human practice in Amazonia. Effective conservation of Amazonian agrobiodiversity needs to incorporate this interaction between wild and cultivated populations that is managed by smallholder families. Beyond treegourd, our study clearly shows that hybridization plays an important role in tree crop phenotypic diversification and that the integration of molecular analyses and farmers’ perceptions of diversity help disentangle crop domestication history.

Myosin XI is associated with fitness and adaptation to aridity in wild pearl millet

Phenotypic changes in plants can be observed along many environmental gradients and are determined by both environmental and genetic factors. The identification of alleles associated with phenotypic variations is a rapidly developing area of research. We studied the genetic basis of phenotypic variations in 11 populations of wild pearl millet (Pennisetum glaucum) on two North-South aridity gradients, one in Niger and one in Mali. Most of the 11 phenotypic traits assessed in a common garden experiment varied between the populations studied. Moreover, the size of the inflorescence, the number of flowers and aboveground dry mass co-varied positively with a decrease in rainfall. To decipher the genetic basis of these phenotypes, we used an association mapping strategy with a mixed model. We found two SNPs on the same myosin XI contig significantly associated with variations in the average number of flowers. Both the allele frequency of the two SNPs and the average number of flowers co-varied with the rainfall gradient on the two gradients. Interestingly, this gene was also a target of selection during domestication. The Myosin XI gene is thus a good candidate for fitness-related adaptation in wild populations.

Mutations in Rice yellow mottle virus Polyprotein P2a Involved in RYMV2 Gene Resistance Breakdown

Rice yellow mottle virus (RYMV) is one of the major diseases of rice in Africa. The high resistance of the Oryza glaberrima Tog7291 accession involves a null allele of the RYMV2 gene, whose ortholog in Arabidopsis, CPR5, is a transmembrane nucleoporin involved in effector-triggered immunity. To optimize field deployment of the RYMV2 gene and improve its durability, which is often a weak point in varietal resistance, we analyzed its efficiency toward RYMV isolates representing the genetic diversity of the virus and the molecular basis of resistance breakdown. Tog7291 resistance efficiency was highly variable depending on the isolate used, with infection rates ranging from 0 to 98% of plants. Back-inoculation experiments indicated that infection cases were not due to an incomplete resistance phenotype but to the emergence of resistance-breaking (RB) variants. Interestingly, the capacity of the virus to overcome Tog7291 resistance is associated with a polymorphism at amino-acid 49 of the VPg protein which also affects capacity to overcome the previously studied RYMV1 resistance gene. This polymorphism appeared to be a main determinant of the emergence of RB variants. It acts independently of the resistance gene and rather reflects inter-species adaptation with potential consequences for the durability of resistance. RB mutations were identified by full-length or partial sequencing of the RYMV genome in infected Tog7291 plants and were validated by directed mutagenesis of an infectious viral clone. We found that Tog7291 resistance breakdown involved mutations in the putative membrane anchor domain of the polyprotein P2a. Although the precise effect of these mutations on rice/RYMV interaction is still unknown, our results offer a new perspective for the understanding of RYMV2 mediated resistance mechanisms. Interestingly, in the susceptible IR64 variety, RB variants showed low infectivity and frequent reversion to the wild-type genotype, suggesting that Tog7291 resistance breakdown is associated with a major loss of viral fitness in normally susceptible O. sativa varieties. Despite the high frequency of resistance breakdown in controlled conditions, this loss of fitness is an encouraging element with regards to RYMV2 resistance durability.

Genome scan reveals selection acting on genes linked to stress response in wild pearl millet

Uncovering genomic regions involved in adaption is a major goal in evolutionary biology. High-throughput sequencing now makes it possible to tackle this challenge in nonmodel species. Yet, despite the increasing number of methods targeted to specifically detect genomic footprints of selection, the complex demography of natural populations often causes high rates of false positive in gene discoveries. The aim of this study was to identify climate adaptations in wild pearl millet populations, Cenchrus americanus ssp. monodii. We focused on two climate gradients, one in Mali and one in Niger. We used a two-step strategy to limit false-positive outliers. First, we considered gradients as biological replicates and performed RNA sequencing of four populations at the extremities. We combined four methods-three based on differentiation among populations and one based on diversity patterns within populations-to identify outlier SNPs from a set of 87 218 high-quality SNPs. Among 11 155 contigs of pearl millet reference transcriptome, 540 exhibited selection signals as evidenced by at least one of the four methods. In a second step, we genotyped 762 samples in 11 additional populations distributed along the gradients using SNPs from the detected contigs and random SNPs as control. We further assessed selection on this large data set using a differentiation-based method and a method based on correlations with environmental variables based. Four contigs displayed consistent signatures between the four extreme and 11 additional populations, two of which were linked to abiotic and biotic stress responses.

Chloroplast sequence of treegourd (Crescentia cujete, Bignoniaceae) to study phylogeography and domestication

PREMISE OF THE STUDY

Crescentia cujete (Bignoniaceae) fruit rinds are traditionally used for storage vessels and handicrafts. We assembled its chloroplast genome and identified single-nucleotide polymorphisms (SNPs).

METHODS AND RESULTS

Using a genome skimming approach, the whole chloroplast of C. cujete was assembled using 3,106,928 sequence reads of 150 bp. The chloroplast is 154,662 bp in length, structurally divided into a large single copy region (84,788 bp), a small single copy region (18,299 bp), and two inverted repeat regions (51,575 bp) with 88 genes annotated. By resequencing the whole chloroplast, we identified 66 SNPs in C. cujete (N = 30) and 68 SNPs in C. amazonica (N = 6). Nucleotide diversity was estimated at 1.1 × 10^-3 and 3.5 × 10^-3 for C. cujete and C. amazonica, respectively.

CONCLUSIONS

This broadened C. cujete genetic toolkit will be important to study the origin, domestication, diversity, and phylogeography of treegourds in the Neotropics.

Spatial and Temporal Variation in Selection of Genes Associated with Pearl Millet Varietal Quantitative Traits In situ

Ongoing global climate changes imply new challenges for agriculture. Whether plants and crops can adapt to such rapid changes is still a widely debated question. We previously showed adaptation in the form of earlier flowering in pearl millet at the scale of a whole country over three decades. However, this analysis did not deal with variability of year to year selection. To understand and possibly manage plant and crop adaptation, we need more knowledge of how selection acts in situ. Is selection gradual, abrupt, and does it vary in space and over time? In the present study, we tracked the evolution of allele frequency in two genes associated with pearl millet phenotypic variation in situ. We sampled 17 populations of cultivated pearl millet over a period of 2 years. We tracked changes in allele frequencies in these populations by genotyping more than seven thousand individuals. We demonstrate that several allele frequencies changes are compatible with selection, by correcting allele frequency changes associated with genetic drift. We found marked variation in allele frequencies from year to year, suggesting a variable selection effect in space and over time. We estimated the strength of selection associated with variations in allele frequency. Our results suggest that the polymorphism maintained at the genes we studied is partially explained by the spatial and temporal variability of selection. In response to environmental changes, traditional pearl millet varieties could rapidly adapt thanks to this available functional variability.

Phylogeography of the genus Podococcus (Palmae/Arecaceae) in Central African rain forests: Climate stability predicts unique genetic diversity

The tropical rain forests of Central Africa contain high levels of species diversity. Paleovegetation or biodiversity patterns suggested successive contraction/expansion phases on this rain forest cover during the last glacial maximum (LGM). Consequently, the hypothesis of the existence of refugia e.g. habitat stability that harbored populations during adverse climatic periods has been proposed. Understory species are tightly associated to forest cover and consequently are ideal markers of forest dynamics. Here, we used two central African rain forest understory species of the palm genus, Podococcus, to assess the role of past climate variation on their distribution and genetic diversity. Species distribution modeling in the present and at the LGM was used to estimate areas of climatic stability. Genetic diversity and phylogeography were estimated by sequencing near complete plastomes for over 120 individuals. Areas of climatic stability were mainly located in mountainous areas like the Monts de Cristal and Monts Doudou in Gabon, but also lowland coastal forests in southeast Cameroon and northeast Gabon. Genetic diversity analyses shows a clear North-South structure of genetic diversity within one species. This divide was estimated to have originated some 500,000years ago. We show that, in Central Africa, high and unique genetic diversity is strongly correlated with inferred areas of climatic stability since the LGM. Our results further highlight the importance of coastal lowland rain forests in Central Africa as harboring not only high species diversity but also important high levels of unique genetic diversity. In the context of strong human pressure on coastal land use and destruction, such unique diversity hotspots need to be considered in future conservation planning.

Genotyping-by-Sequencing SNP Identification for Crops without a Reference Genome: Using Transcriptome Based Mapping as an Alternative Strategy

Next-generation sequencing opens the way for genomic studies of diversity even for non-model crops and animals. Genome reduction techniques are becoming progressively more popular as they allow a fraction of the genome to be sequenced for multiple individuals and/or populations. These techniques are an efficient way to explore genome diversity in non-model crops and animals for which no reference genome is available. Genome reduction techniques emerged with the development of specific pipelines such as UNEAK (Universal Network Enabled Analysis Kit) and Stacks. However, even for non-model crops and animals, transcriptomes are easier to obtain, thereby making it possible to directly map reads. We investigate the direct use of transcriptome as an alternative strategy. Our specific objective was to compare SNPs obtained from the UNEAK pipeline as well as SNPs obtained by directly mapping genotyping-by-sequencing reads on a transcriptome. We assessed the feasibility of both SNP datasets, UNEAK and transcriptome mapping, to investigate the diversity of 91 samples of wild pearl millet sampled across its distribution area. Both approaches produced several tens of thousands of single nucleotide variants, but differed in the way the variants were identified, leading to differences in the frequency spectrum associated with marked differences in the assessment of diversity. Difference in the frequency spectrum significantly biased a large set of diversity analyses as well as detection of selection approaches. However, whatever the approach, we found very similar inference of genetic structure, with three major genetic groups from West, Central, and East Africa. For non-model crops, using transcriptome data as a reference is thus a particularly promising way to obtain a more thorough analysis of datasets generated using genome reduction techniques.

Association studies including genotype by environment interactions: prospects and limits

Background

Association mapping studies offer great promise to identify polymorphisms associated with phenotypes and for understanding the genetic basis of quantitative trait variation. To date, almost all association mapping studies based on structured plant populations examined the main effects of genetic factors on the trait but did not deal with interactions between genetic factors and environment. In this paper, we propose a methodological prospect of mixed linear models to analyze genotype by environment interaction effects using association mapping designs. First, we simulated datasets to assess the power of linear mixed models to detect interaction effects. This simulation was based on two association panels composed of 90 inbreds (pearl millet) and 277 inbreds (maize).

Results

Based on the simulation approach, we reported the impact of effect size, environmental variation, allele frequency, trait heritability, and sample size on the power to detect the main effects of genetic loci and diverse effect of interactions implying these loci. Interaction effects specified in the model included SNP by environment interaction, ancestry by environment interaction, SNP by ancestry interaction and three way interactions. The method was finally used on real datasets from field experiments conducted on the two considered panels. We showed two types of interactions effects contributing to genotype by environment interactions in maize: SNP by environment interaction and ancestry by environment interaction. This last interaction suggests differential response at the population level in function of the environment.

Conclusions

Our results suggested the suitability of mixed models for the detection of diverse interaction effects. The need of samples larger than that commonly used in current plant association studies is strongly emphasized to ensure rigorous model selection and powerful interaction assessment. The use of ancestry interaction component brought valuable information complementary to other available approaches.

Association mapping, patterns of linkage disequilibrium and selection in the vicinity of the PHYTOCHROME C gene in pearl millet

Linkage analysis confirmed the association in the region of PHYC in pearl millet. The comparison of genes found in this region suggests that PHYC is the best candidate. Major efforts are currently underway to dissect the phenotype-genotype relationship in plants and animals using existing populations. This method exploits historical recombinations accumulated in these populations. However, linkage disequilibrium sometimes extends over a relatively long distance, particularly in genomic regions containing polymorphisms that have been targets for selection. In this case, many genes in the region could be statistically associated with the trait shaped by the selected polymorphism. Statistical analyses could help in identifying the best candidate genes into such a region where an association is found. In a previous study, we proposed that a fragment of the PHYTOCHROME C gene (PHYC) is associated with flowering time and morphological variations in pearl millet. In the present study, we first performed linkage analyses using three pearl millet F2 families to confirm the presence of a QTL in the vicinity of PHYC. We then analyzed a wider genomic region of ~100 kb around PHYC to pinpoint the gene that best explains the association with the trait in this region. A panel of 90 pearl millet inbred lines was used to assess the association. We used a Markov chain Monte Carlo approach to compare 75 markers distributed along this 100-kb region. We found the best candidate markers on the PHYC gene. Signatures of selection in this region were assessed in an independent data set and pointed to the same gene. These results foster confidence in the likely role of PHYC in phenotypic variation and encourage the development of functional studies.

Evolutionary history of pearl millet (Pennisetum glaucum [L.] R. Br.) and selection on flowering genes since its domestication

The plant domestication process is associated with considerable modifications of plant phenotype. The identification of the genetic basis of this adaptation is of great interest for evolutionary biology. One of the methods used to identify such genes is the detection of signatures of selection. However, domestication is generally associated with major demographic effects. It is therefore crucial to disentangle the effects of demography and selection on diversity. In this study, we investigated selection in a flowering time pathway during domestication of pearl millet. We first used a random set of 20 genes to model pearl millet domestication using approximate Bayesian computation. This analysis showed that a model with exponential growth and wild-cultivated gene flow was well supported by our data set. Under this model, the domestication date of pearl millet is estimated at around 4,800 years ago. We assessed selection in 15 pearl millet DNA sequences homologous to flowering time genes and showed that these genes underwent selection more frequently than expected. We highlighted significant signatures of selection in six pearl millet flowering time genes associated with domestication or improvement of pearl millet. Moreover, higher deviations from neutrality were found for circadian clock-associated genes. Our study provides new insights into the domestication process of pearl millet and shows that a category of genes of the flowering pathway were preferentially selected during pearl millet domestication.

Genetic basis of pearl millet adaptation along an environmental gradient investigated by a combination of genome scan and association mapping

Identifying the molecular bases of adaptation is a key issue in evolutionary biology. Genome scan is an efficient approach for identifying important molecular variation involved in adaptation. Association mapping also offers an opportunity to gain insight into genotype-phenotype relationships. Using these two approaches coupled with environmental data should help to come up with a refined picture of the evolutionary process underlying adaptation. In this study, we first conducted a selection scan analysis on a transcription factor gene family. We focused on the MADS-box gene family, a gene family which plays a crucial role in vegetative and flower development. Twenty-one pearl millet populations were sampled along an environmental gradient in West Africa. We identified one gene, i.e. PgMADS11, using Bayesian analysis to detect selection signatures. Polymorphism at this gene was also associated with flowering time variation in an association mapping framework. Finally, we found that PgMADS11 allele frequencies were closely associated with annual rainfall. Overall, we determined an efficient way to detect functional polymorphisms associated with climate variation in non-model plants by combining genome scan and association mapping. These results should help monitor the impact of recent climatic changes on plant adaptation.

Spatio-temporal dynamics of genetic diversity in Sorghum bicolor in Niger

The dynamics of crop genetic diversity need to be assessed to draw up monitoring and conservation priorities. However, few surveys have been conducted in centres of diversity. Sub-Saharan Africa is the centre of origin of sorghum. Most Sahel countries have been faced with major human, environmental and social changes in recent decades, which are suspected to cause genetic erosion. Sorghum is the second staple cereal in Niger, a centre of diversity for this crop. Niger was submitted to recurrent drought period and to major social changes during these last decades. We report here on a spatio-temporal analysis of sorghum genetic diversity, conducted in 71 villages covering the rainfall gradient and range of agro-ecological conditions in Niger's agricultural areas. We used 28 microsatellite markers and applied spatial and genetic clustering methods to investigate change in genetic diversity over a 26-year period (1976-2003). Global genetic differentiation between the two collections was very low (F (st) = 0.0025). Most of the spatial clusters presented no major differentiation, as measured by F (st), and showed stability or an increase in allelic richness, except for two of them located in eastern Niger. The genetic clusters identified by Bayesian analysis did not show a major change between the two collections in the distribution of accessions between them or in their spatial location. These results suggest that farmers' management has globally preserved sorghum genetic diversity in Niger.

Phylogeny and origin of pearl millet (Pennisetum glaucum [L.] R. Br) as revealed by microsatellite loci

During the last 12,000 years, different cultures around the world have domesticated cereal crops. Several studies investigated the evolutionary history and domestication of cereals such as wheat in the Middle East, rice in Asia or maize in America. The domestication process in Africa has led to the emergence of important cereal crops like pearl millet in Sahelian Africa. In this study, we used 27 microsatellite loci to analyze 84 wild accessions and 355 cultivated accessions originating from the whole pearl millet distribution area in Africa and Asia. We found significantly higher diversity in the wild pearl millet group. The cultivated pearl millet sample possessed 81% of the alleles and 83% of the genetic diversity of the wild pearl millet sample. Using Bayesian approaches, we identified intermediate genotypes between the cultivated and wild groups. We then analyzed the phylogenetic relationship among accessions not showing introgression and found that a monophyletic origin of cultivated pearl millet in West Africa is the most likely scenario supported by our data set.

Niger-wide assessment of in situ sorghum genetic diversity with microsatellite markers

Understanding the geographical, environmental and social patterns of genetic diversity on different spatial scales is key to the sustainable in situ management of genetic resources. However, few surveys have been conducted on crop genetic diversity using exhaustive in situ germplasm collections on a country scale and such data are missing for sorghum in sub-Saharan Africa, its centre of origin. We report here a genetic analysis of 484 sorghum varieties collected in 79 villages evenly distributed across Niger, using 28 microsatellite markers. We found a high level of SSR diversity in Niger. Diversity varied between eastern and western Niger, and allelic richness was lower in the eastern part of the country. Genetic differentiation between botanical races was the first structuring factor (Fst = 0.19), but the geographical distribution and the ethnic group to which farmers belonged were also significantly associated with genetic diversity partitioning. Gene pools are poorly differentiated among climatic zones. The geographical situation of Niger, where typical western African (guinea), central African (caudatum) and eastern Sahelian African (durra) sorghum races converge, explained the high observed genetic diversity and was responsible for the interactions among the ethnic, geographical and botanical structure revealed in our study. After correcting for the structure of botanical races, spatial correlation of genetic diversity was still detected within 100 km, which may hint at limited seed exchanges between farmers. Sorghum domestication history, in relation to the spatial organisation of human societies, is therefore key information for sorghum in situ conservation programs in sub-Saharan Africa.

Role of seed flow on the pattern and dynamics of pearl millet (Pennisetum glaucum [L.] R. Br.) genetic diversity assessed by AFLP markers: a study in south-western Niger

We studied the regional genetic diversity and seed exchange dynamics of pearl millet landraces in south-western Niger. The genetic study was based on AFLP markers. We found significant genetic differentiation between landraces in different geographical areas of south-western Niger. However, the degree of differentiation was low insofar as only 1.9% of the total molecular diversity was due to regional differentiation, suggesting a relatively high gene flow. Anthropologic studies on farming practices have suggested that seed exchanges between farmers on a large geographical scale probably make a considerable contribution to this result. In order to test this hypothesis, the effects of seed exchange on the genetic diversity of landraces was analyzed on seed samples from two distant villages in contrasting areas of south-western Niger. Seeds imported by farmers into the southern village of Sina Koara did not differ significantly from locally grown landraces. By contrast, in the northern village of Alzou, several samples were genetically different from locally grown landraces and closer to southern accessions. These data suggest that the seed flow is preferentially from south to north, i.e. from an area with more favorable rainfall conditions. The potential consequences for the genetic diversity and adaptation of northern pearl millet landraces are discussed.

Genetic diversity and gene flow among pearl millet crop/weed complex: a case study

Weedy plants with intermediate (domesticated x wild) phenotypes occur in most pearl millet fields in West Africa, even in the absence of wild populations. They are usually found, in high numbers, both inside and outside of drills. Questions pertaining to the evolutionary dynamics of diversity within the pearl millet complex (domesticated-weedy-wild forms) were addressed in this study. The diversity of the different components of this complex sampled in two pearl millet fields in two villages of southwestern Niger was assessed at both molecular (AFLP) and morphological levels. Results show that, in both fields, weedy plants found outside of drills are morphologically distinct from weedy plants found inside drills, despite their close similarity at AFLP markers. The data suggest some introgression from the wild to the weedy population but nevertheless that the gene flow between the parapatric wild and domesticated populations is very low. This challenges the traditional view that regular hybridization between domesticated and wild pearl millets explains the abundance of these weedy plants despite farmers' seed selection. The level of genetic differentiation between fields from the two villages was low when considering domesticated and weedy plants. This could be explained by high gene flow resulting from substantial seed exchanges between farmers. The fact that it is very difficult for farmers to keep their own selected seeds, and the consequent substantial seed exchanges between them, is probably the main factor accounting for the maintenance and dispersal of weedy pearl millets in the region, even in areas where no wild forms have been observed.

Haplotype variation of cpDNA in the agamic grass complex Pennisetum section Brevivalvula (Poaceae)

The Brevivalvula section of the grass polyploid complex Pennisetum shows various reproductive systems, apomixis being the most widespread. Haplotype variation of chloroplast DNA was studied in the six morphological taxa (species) of this section by using RFLP analysis in 54 plants corresponding to 14 elementary taxa, each characterized on the basis of morphology and ploidy level. Two additional species, Pennisetum glaucum and P. purpureum, which belong to another section of the same genus, were analysed for comparison. In Brevivalvula, chloroplast DNA size was estimated to range between 130 and 133 kb. Thirteen of 15 distinct haplotypes identified in the study were specific to the Brevivalvula section. They were unequally distributed among the morphotypes, the ploidy levels and sampling sites. Within the Brevivalvula section, plants of P. setosum, which are perennial and reproduce vegetatively or by agamospermy, possessed a single specific haplotype. This species differed clearly from the five other morphological species, which are known to be annual, to show either sexual or agamospermic reproduction and which shared most of the 12 other haplotypes observed in the section, suggesting the occurrence of multiple hybridization events between the taxa. Chloroplast DNA variation was highly geographically structured, suggesting low seed dispersal between sites, whereas the substantial haplotype diversity observed in the sites may indicate that agamic reproduction is responsible for the maintenance of distinct genetically isolated clones. Haplotype classification using Wagner's parsimony suggested the occurrence of bidirectional gene flow between the diploids and the polyploids, as reported already in other related apomictic complexes.

Chloroplast DNA extraction from herbaceous and woody plants for direct restriction fragment length polymorphism analysis

The technique described here is a fast and simple method of extracting chloroplast DNA (cpDNA). It overcomes the need for differential centrifugation using density gradients. The leaves do not have to be kept in the dark and lyophilized before extraction, but lyophilization is still possible. The chloroplasts are specifically lysed in a cell extract of leaves, using a non-ionic detergent. After isolation by centrifugation, the cpDNA is purified by the combined action of proteolytic enzymes and detergents, followed by the elimination of proteins using a mixture of chloroform and isoamyl alcohol. This method provided good quality restriction profiles for all species analyzed.

Expression of a self-incompatibility gene in a self-compatible line of Brassica oleracea

In cruciferous plants, self-pollination is prevented by the action of genes situated at the self-incompatibility locus or S-locus. The self-incompatibility reaction is associated with expression of stigma glycoproteins encoded by the S-locus glycoprotein (SLG) gene. Only a few cases of self-compatible plants derived from self-incompatible lines in the crucifer Brassica have been reported. In these cases, self-compatibility was generally ascribed to the action of single genes unlinked to the S-locus. In contrast, we report here a line of Brassica oleracea var acephala with a self-compatible phenotype linked to the S-locus. By means of both biochemical and immunochemical analyses, we showed that this self-compatible (Sc) line nonetheless possesses stigmatic SLGs (SLG-Sc) that are expressed with a similar spatial and temporal pattern to that described for the SLGs of self-incompatible Brassica plants. Moreover, the SLG-Sc products segregate with the self-compatibility phenotype in F2 progeny, suggesting that changes at the S-locus may be responsible for the occurrence of the self-compatibility character. A cDNA clone encoding the SLG-Sc product was isolated, and the deduced amino acid sequence showed this glycoprotein to be highly homologous to the pollen recessive S2 allele glycoprotein. Hence, self-compatibility in this Brassica Sc line correlates with the expression of a pollen recessive-like S allele in the stigma.