Cost-effective enrichment hybridization capture of chloroplast genomes at deep multiplexing levels for population genetics and phylogeography studies

Development and characterization of nuclear microsatellite markers for the African walnut Coula edulis Baill (Coulaceae)

PREMISE OF THE STUDY

Coula edulis Baill (Coulaceae) is a common tree species in the Guineo-Congolian forests producing an edible fruit known as African walnut, which is an important food and income resource for rural populations. However, the species suffers from a deficit of natural regeneration. We developed here nuclear microsatellite markers for C. edulis to be able to study the genetic structure of its natural populations and gene flow.

METHODS AND RESULTS

A genomic library was obtained using the Illumina platform, and 21 polymorphic microsatellite loci were developed. The polymorphic microsatellites displayed eight to 22 alleles per locus (average: 14.2), with a mean expected heterozygosity ranging from 0.33 to 0.72 in five populations from Central and West Africa.

CONCLUSIONS

The high polymorphism of the nuclear microsatellite markers developed makes them useful to investigate gene flow and the organization of genetic diversity in C. edulis, and to assess whether particular genetic resources require conservation efforts.

Complete chloroplast genome of two nutmeg species, Myristica argentea Warb. 1891 and Myristica fatua Houtt. 1774 (Myristicaceae)

Myristica argentea Warb. 1891 and M. fatua Houtt. 1774 are two South-East Asian food tree species. They are harvested from the wild or cultivated for local uses as a condiment (nutmeg and mace), medicine, and source of wood. In this study, we reconstructed the complete chloroplast (cp) genomes of these two species from whole genome sequencing data using the Illumina NovaSeq platform. The genome sizes of M. argentea and M. fatua were respectively 155,871 base pairs (bp) and 155,898 bp, including 126 genes and an overall GC content of 39.20% in both species. Our study provides useful resources for future evolutionary research and diversity analysis of Myristica species.

Rapid evolution of an RNA virus to escape recognition by a rice nucleotide-binding and leucine-rich repeat domain immune receptor

Viral diseases are a major limitation for crop production, and their control is crucial for sustainable food supply. We investigated by a combination of functional genetics and experimental evolution the resistance of rice to the rice yellow mottle virus (RYMV), which is among the most devastating rice pathogens in Africa, and the mechanisms underlying the extremely fast adaptation of the virus to its host. We found that the RYMV3 gene that protects rice against the virus codes for a nucleotide-binding and leucine-rich repeat domain immune receptor (NLRs) from the Mla-like clade of NLRs. RYMV3 detects the virus by forming a recognition complex with the viral coat protein (CP). The virus escapes efficiently from detection by mutations in its CP, some of which interfere with the formation of the recognition complex. This study establishes that NLRs also confer in monocotyledonous plants immunity to viruses, and reveals an unexpected functional diversity for NLRs of the Mla clade that were previously only known as fungal disease resistance proteins. In addition, it provides precise insight into the mechanisms by which viruses adapt to plant immunity and gives important knowledge for the development of sustainable resistance against viral diseases of cereals.

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.

Development and characterization of microsatellite markers in the African timber tree species Cylicodiscus gabunensis (Fabaceae)

BACKGROUND

Cylicodiscus gabunensis (Fabaceae) or 'Okan' is a Central African multipurpose timber species that is used for heavy construction and traditional medicine. Despite being currently heavily exploited, the species shows a low population density and a natural regeneration deficit in dense forest.

METHODS AND RESULTS

We aimed to characterize polymorphic microsatellite markers that can be used to study patterns of genetic structure and gene flow (mating system, pollen and seed dispersal and inbreeding depression) and ultimately, help to develop sustainable forest management practices. We developed 24 polymorphic markers that can be amplified in three PCR multiplexes that were tested in 647 individuals of C. gabunensis from three populations. The number of alleles ranged from 4 to 27 and the average observed and expected heterozygosity across loci and populations were Ho = 0.585 (0.081-0.936) and He = 0.510 (0.076-0.914), respectively.

CONCLUSIONS

This set of markers is a useful tool for exploring intra-specific diversity, genetic structure and gene flow of C. gabunensis.

Bactris gasipaes Kunth var. gasipaes complete plastome and phylogenetic analysis

Bactris gasipaes var. gasipaes (Arecaceae, Palmae) is an economically and socially important plant species for populations across tropical South and Central America. It has been domesticated from its wild variety, B. gasipaes var. chichagui, since pre-Columbian times. In this study, we sequenced the plastome of the cultivated variety, B. gasipaes Kunth var. gasipaes and compared it with the published plastome of the wild variety. The chloroplast sequence obtained was 156,580 bp. The cultivated chloroplast sequence was conserved compared to the wild type sequence with 99.8% of nucleotide identity. We did, however, identify multiple Single Nucleotide Variants (SNVs), insertions, microsatellites and a resolved region of missing nucleotides. A SNV in one of the core barcode markers (matK) was detected between the wild and cultivated accessions. Phylogenetic analysis was carried out across the Arecaceae family and compared to previous reports, resulting in an identical topology. This study is a step forward in understanding the genome evolution of this species.

Sequence Capture of Mitochondrial Genome with PCR-Generated Baits Provides New Insights into the Biogeography of the Genus Abies Mill

Mitochondrial DNA (mtDNA), being maternally inherited in plants of the family Pinaceae, is an important source of phylogeographic information. However, its use is hindered by a low mutation rate and frequent structure rearrangements. In the present study, we tested the method of genomic libraries enrichment with mtDNA via the sequence capture method yielding mtDNA data which were further used to reconstruct the phylogenetic tree of the genus Abies. The baits for hybrid capture were obtained by long-range PCR using primers designed on the basis of the assembly of Abies sibirica Ledeb. mitochondrial genome. Mitochondrial genomes of Picea sitchensis (Bong.) Carr., Larix sibirica Ledeb., and Keteleeria davidiana (Bertrand) Beissn. were used as an outgroup. The resulting phylogenetic tree consists of two sister branches, including the Eurasian and American species, respectively, with some exceptions. The subclade of A. sachalinensis (F. Schmidt) Mast. and A. veitchii Lindl. (Japan and Sakhalin islands) occupies a basal position in the branch of American firs, probably due to the complex history of fir migrations from North America to Eurasia. The tree has high support for majority of clades. For species represented by more than one sample an intraspecific variability was found which is suitable to design mtDNA markers for phylogeographic and population studies.

Phylogenomic relationships and historical biogeography in the South American vegetable ivory palms (Phytelepheae)

The palm tribe Phytelepheae form a clade of three genera and eight species whose phylogenetic relationships and historical biogeography are not fully understood. Based on morphological similarities and phylogenetic relatedness, it has been suggested that Phytelephas seemannii and Phytelephas schottii are synonyms of Phytelephas macrocarpa, implying the existence of only six species within the Phytelepheae. In addition, uncertainty in their phylogenetic relationships in turn results in blurred biogeographic history. We inferred the phylogenomic relationships in the Phytelepheae by target-capturing 176 nuclear genes and estimated divergence times by using four fossils for time calibration. We lastly explored the biogeographic history of the tribe by inferring its ancestral range evolution. Our phylogenomic trees showed that P. seemannii and P. schottii are not closely related with P. macrocarpa, and therefore, support the existence of eight species in the Phytelepheae. The ancestor of the tribe was widely-distributed in the Chocó, Magdalena, and Amazonia during the Miocene at 19.25 Ma. Early diversification in Phytelephas at 5.27 Ma could have occurred by trans-Andean vicariance after the western Andes uplifted rapidly at ∼ 10 Ma. Our results show the utility of phylogenomic approaches to shed light on species relationships and their biogeographic history.

Characterizing Microbiomes via Sequencing of Marker Loci: Techniques To Improve Throughput, Account for Cross-Contamination, and Reduce Cost

New approaches to characterizing microbiomes via high-throughput sequencing provide impressive gains in efficiency and cost reduction compared to approaches that were standard just a few years ago. However, the speed of method development has been such that staying abreast of the latest technological advances is challenging. Moreover, shifting laboratory protocols to include new methods can be expensive and time consuming. To facilitate adoption of new techniques, we provide a guide and review of recent advances that are relevant for single-locus sequence-based study of microbiomes—from extraction to library preparation—including a primer regarding the use of liquid-handling automation in small-scale academic settings. Additionally, we describe several amendments to published techniques to improve throughput, track contamination, and reduce cost. Notably, we suggest adding synthetic DNA molecules to each sample during nucleic acid extraction, thus providing a method of documenting incidences of cross-contamination. We also describe a dual-indexing scheme for Illumina sequencers that allows multiplexing of many thousands of samples with minimal PhiX input. Collectively, the techniques that we describe demonstrate that laboratory technology need not impose strict limitations on the scale of molecular microbial ecology studies.

IMPORTANCE New methods to characterize microbiomes reduce technology-imposed limitations to study design, but many new approaches have not been widely adopted. Here, we present techniques to increase throughput and reduce contamination alongside a thorough review of current best practices.

Genetic Diversity of Rice stripe necrosis virus and New Insights into Evolution of the Genus Benyvirus

The rice stripe necrosis virus (RSNV) has been reported to infect rice in several countries in Africa and South America, but limited genomic data are currently publicly available. Here, eleven RSNV genomes were entirely sequenced, including the first corpus of RSNV genomes of African isolates. The genetic variability was differently distributed along the two genomic segments. The segment RNA1, within which clusters of polymorphisms were identified, showed a higher nucleotidic variability than did the beet necrotic yellow vein virus (BNYVV) RNA1 segment. The diversity patterns of both viruses were similar in the RNA2 segment, except for an in-frame insertion of 243 nucleotides located in the RSNV tgbp1 gene. Recombination events were detected into RNA1 and RNA2 segments, in particular in the two most divergent RSNV isolates from Colombia and Sierra Leone. In contrast to BNYVV, the RSNV molecular diversity had a geographical structure with two main RSNV lineages distributed in America and in Africa. Our data on the genetic diversity of RSNV revealed unexpected differences with BNYVV suggesting a complex evolutionary history of the genus Benyvirus.

Green giant-a tiny chloroplast genome with mighty power to produce high-value proteins: history and phylogeny

Free-living cyanobacteria were entrapped by eukaryotic cells ~2 billion years ago, ultimately giving rise to chloroplasts. After a century of debate, the presence of chloroplast DNA was demonstrated in the 1960s. The first chloroplast genomes were sequenced in the 1980s, followed by ~100 vegetable, fruit, cereal, beverage, oil and starch/sugar crop chloroplast genomes in the past three decades. Foreign genes were expressed in isolated chloroplasts or intact plant cells in the late 1980s and stably integrated into chloroplast genomes, with typically maternal inheritance shown in the 1990s. Since then, chloroplast genomes conferred the highest reported levels of tolerance or resistance to biotic or abiotic stress. Although launching products with agronomic traits in important crops using this concept has been elusive, commercial products developed include enzymes used in everyday life from processing fruit juice, to enhancing water absorption of cotton fibre or removal of stains as laundry detergents and in dye removal in the textile industry. Plastid genome sequences have revealed the framework of green plant phylogeny as well as the intricate history of plastid genome transfer events to other eukaryotes. Discordant historical signals among plastid genes suggest possible variable constraints across the plastome and further understanding and mitigation of these constraints may yield new opportunities for bioengineering. In this review, we trace the evolutionary history of chloroplasts, status of autonomy and recent advances in products developed for everyday use or those advanced to the clinic, including treatment of COVID-19 patients and SARS-CoV-2 vaccine.

Utility of the Mitochondrial Genome in Plant Taxonomic Studies

Size, structure, and sequence content lability of plant mitochondrial genome (mtDNA) across species has sharply limited its use in taxonomic studies. Historically, mtDNA variation has been first investigated with RFLPs, while the development of universal primers then allowed studying sequence polymorphisms within short genomic regions (<3 kb). The recent advent of NGS technologies now offers new opportunities by greatly facilitating the assembly of longer mtDNA regions, and even full mitogenomes. Phylogenetic works aiming at comparing signals from different genomic compartments (i.e., nucleus, chloroplast, and mitochondria) have been developed on a few plant lineages, and have been shown especially relevant in groups with contrasted inheritance of organelle genomes. This chapter first reviews the main characteristics of mtDNA and the application offered in taxonomic studies. It then presents tips for best sequencing protocol based on NGS data to be routinely used in mtDNA-based phylogenetic studies.

Genomic Profiling: The Strengths and Limitations of Chloroplast Genome-Based Plant Variety Authentication

Genomic profiling is a suitable tool for variety authentication and has applications in both operational quality and regulatory raw material control. It can be used to differentiate species or varieties and to identify admixtures as well as field contaminants. To establish a molecular profile, reliable and very accurate sequence data are required. As a result of the influence of the pollinator plant, nuclear genome-based authentication is in most cases not suitable for a direct application on the fruit. Sequences must be used that come exclusively from the localized mother plant. Parts of the fruit of maternal origin, e.g., components derived from the blossom, are suitable as a basis for this. Alternatively, DNA from cell organelles that are maternally inherited, such as mitochondria or chloroplasts, can be used. The latter will be discussed in this review in closer detail. Although individual gene segments on the chloroplast genome are already used for species differentiation in barcoding studies on plants, little is known about the usefulness of the entire chloroplast genome for intraspecies differentiation in general and for differentiation between modern varieties in particular. Results from the literature as well as from our own work suggest that chloroplast genome sequences are indeed very well-suited for the differentiation of old varieties. On the other hand, they are less or not suitable for the genetic differentiation of modern cultivars, because they are often too closely related.

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.

Sociogenetic Organization of the Red Honey Ant (Melophorus bagoti)

Kin selection and inclusive fitness are thought to be key factors explaining the reproductive altruism displayed by workers in eusocial insect species. However, when a colony's queen has mated with <2 males, workers may increase their fitness by producing their own male offspring. Conversely, when the queen has mated with ≥2 males, workers are expected to increase their inclusive fitness by eschewing the production of their sons and preventing other workers from reproducing as well. Here, we investigated sociogenetic structure and worker reproduction in the red honey ant, Melophorus bagoti. Morphometric analyses revealed that workers belong to one of two distinct subcastes: they are either majors or minors. Using DNA microsatellite markers, we showed that all the colonies had a single, multiple-mated queen and that there was no relationship between worker patriline and worker subcaste. Furthermore, we found that workers were producing males in the presence of the queen, which contrasts with the predictions of inclusive fitness theory. Although our results are based on a small sample, they can serve as the foundation for future research examining worker reproduction in M. bagoti.

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.

Applied Barcoding: The Practicalities of DNA Testing for Herbals

DNA barcoding is a widely accepted technique for the identification of plant materials, and its application to the authentication of commercial medicinal plants has attracted significant attention. The incorporation of DNA-based technologies into the quality testing protocols of international pharmacopoeias represents a step-change in status, requiring the establishment of standardized, reliable and reproducible methods. The process by which this can be achieved for any herbal medicine is described, using Hypericum perforatum L. (St John's Wort) and potential adulterant Hypericum species as a case study. A range of practical issues are considered including quality control of DNA sequences from public repositories and the construction of individual curated databases, choice of DNA barcode region(s) and the identification of informative polymorphic nucleotide sequences. A decision tree informs the structure of the manuscript and provides a template to guide the development of future DNA barcode tests for herbals.

Adaptation Evolution and Phylogenetic Analyses of Species in Chinese Allium Section Pallasia and Related Species Based on Complete Chloroplast Genome Sequences

The section Pallasia is one of the components of the genus Allium subgenus Allium (Amaryllidaceae), and species relationship in this section is still not resolved very well, which hinders further evolutionary and adaptive studies. Here, the complete chloroplast genomes of five sect. Pallasia species were reported, and a comparative analysis was performed with other three related Allium species. The genome size of the eight species ranged from 151,672 bp to 153,339 bp in length, GC content changed from 36.7% to 36.8%, and 130 genes (except Allium pallasii), 37 tRNA, and 8 rRNA were identified in each genome. By analyzing the IR/LSC and IR/SSC boundary, A. pallasii exhibited differences compared with other seven species. Phylogenetic analysis achieved high supports in each branch, seven of the eight Allium species cluster into a group, and A. pallasii exhibit a close relationship with A. obliquum. Higher pairwise Ka/Ks ratios were found in A. schoenoprasoides compared to A. caeruleum and A. macrostemon while a lower value of Ka/Ks ratios was detected between A. caeruleum and A. macrostemon. This study will be a great contribution to the future phylogenetic and adaptive research in Allium.

Comparative Plastid Genomics of Neotropical Bulbophyllum (Orchidaceae; Epidendroideae)

Pantropical Bulbophyllum, with ∼2,200 species, is one of the largest genera in Orchidaceae. Although phylogenetics and taxonomy of the ∼60 American species in the genus are generally well understood, some species complexes need more study to clearly delimit their component species and provide information about their evolutionary history. Previous research has suggested that the plastid genome includes phylogenetic markers capable of providing resolution at low taxonomic levels, and thus it could be an effective tool if these divergent regions can be identified. In this study, we sequenced the complete plastid genome of eight Bulbophyllum species, representing five of six Neotropical taxonomic sections. All plastomes conserve the typical quadripartite structure, and, although the general structure of plastid genomes is conserved, differences in ndh-gene composition and total length were detected. Total length was determined by contraction and expansion of the small single-copy region, a result of an independent loss of the seven ndh genes. Selection analyses indicated that protein-coding genes were generally well conserved, but in four genes, we identified 95 putative sites under positive selection. Furthermore, a total of 54 polymorphic simple sequence repeats were identified, for which we developed amplification primers. In addition, we propose 10 regions with potential to improve phylogenetic analyses of Neotropical Bulbophyllum species.

Plastid genomes of the North American Rhus integrifolia-ovata complex and phylogenomic implications of inverted repeat structural evolution in Rhus L

Plastid genomes (plastomes) represent rich sources of information for phylogenomics, from higher-level studies to below the species level. The genus Rhus (sumac) has received a significant amount of study from phylogenetic and biogeographic perspectives, but genomic studies in this genus are lacking. Rhus integrifolia and R. ovata are two shrubby species of high ecological importance in the southwestern USA and Mexico, where they occupy coastal scrub and chaparral habitats. They hybridize frequently, representing a fascinating system in which to investigate the opposing effects of hybridization and divergent selection, yet are poorly characterized from a genomic perspective. In this study, complete plastid genomes were sequenced for one accession of R. integrifolia and one each of R. ovata from California and Arizona. Sequence variation among these three accessions was characterized, and PCR primers potentially useful in phylogeographic studies were designed. Phylogenomic analyses were conducted based on a robustly supported phylogenetic framework based on 52 complete plastomes across the order Sapindales. Repeat content, rather than the size of the inverted repeat, had a stronger relative association with total plastome length across Sapindales when analyzed with phylogenetic least squares regression. Variation at the inverted repeat boundary within Rhus was striking, resulting in major shifts and independent gene losses. Specifically, rps19 was lost independently in the R. integrifolia-ovata complex and in R. chinensis, with a further loss of rps22 and a major contraction of the inverted repeat in two accessions of the latter. Rhus represents a promising novel system to study plastome structural variation of photosynthetic angiosperms at and below the species level.

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.

Plastome phylogeography in two African rain forest legume trees reveals that Dahomey Gap populations originate from the Cameroon volcanic line

Paleo-environmental data show that the distribution of African rain forests was affected by Quaternary climate changes. In particular, the Dahomey Gap (DG) - a 200 km wide savanna corridor currently separating the West African and Central African rain forest blocks and containing relict rain forest fragments - was forested during the mid-Holocene and possibly during previous interglacial periods, whereas it was dominated by open vegetation (savanna) during glacial periods. Genetic signatures of past population fragmentation and demographic changes have been found in some African forest plant species using nuclear markers, but such events appear not to have been synchronous or shared across species. To better understand the colonization history of the DG by rain forest trees through seed dispersal, the plastid genomes of two widespread African forest legume trees, Anthonotha macrophylla and Distemonanthus benthamianus, were sequenced in 47 individuals for each species, providing unprecedented phylogenetic resolution of their maternal lineages (857 and 115 SNPs, respectively). Both species exhibit distinct lineages separating three regions: 1. Upper Guinea (UG, i.e. the West African forest block), 2. the area ranging from the DG to the Cameroon volcanic line (CVL), and 3. Lower Guinea (LG, the western part of the Central African forest block) where three lineages co-occur. In both species, the DG populations (including southern Nigeria west of Cross River) exhibit much lower genetic diversity than UG and LG populations, and their plastid lineages originate from the CVL, confirming the role of the CVL as an ancient forest refuge. Despite the similar phylogeographic structures displayed by A. macrophylla and D. benthamianus, molecular dating indicates very contrasting ages of lineage divergence (UG diverged from LG since c. 7 Ma and 0.7 Ma, respectively) and DG colonization (probably following the Mid Pleistocene Transition and the Last Glacial Maximum, respectively). The stability of forest refuge areas and repeated similar forest shrinking/expanding events during successive glacial periods might explain why similar phylogeographic patterns can be generated over contrasting timescales.

Strategies for reducing per-sample costs in target capture sequencing for phylogenomics and population genomics in plants

The reduced cost of high-throughput sequencing and the development of gene sets with wide phylogenetic applicability has led to the rise of sequence capture methods as a plausible platform for both phylogenomics and population genomics in plants. An important consideration in large targeted sequencing projects is the per-sample cost, which can be inflated when using off-the-shelf kits or reagents not purchased in bulk. Here, we discuss methods to reduce per-sample costs in high-throughput targeted sequencing projects. We review the minimal equipment and consumable requirements for targeted sequencing while comparing several alternatives to reduce bulk costs in DNA extraction, library preparation, target enrichment, and sequencing. We consider how each of the workflow alterations may be affected by DNA quality (e.g., fresh vs. herbarium tissue), genome size, and the phylogenetic scale of the project. We provide a cost calculator for researchers considering targeted sequencing to use when designing projects, and identify challenges for future development of low-cost sequencing in non-model plant systems.

Can we use it? On the utility of de novo and reference-based assembly of Nanopore data for plant plastome sequencing

The chloroplast genome harbors plenty of valuable information for phylogenetic research. Illumina short-read data is generally used for de novo assembly of whole plastomes. PacBio or Oxford Nanopore long reads are additionally employed in hybrid approaches to enable assembly across the highly similar inverted repeats of a chloroplast genome. Unlike for PacBio, plastome assemblies based solely on Nanopore reads are rarely found, due to their high error rate and non-random error profile. However, the actual quality decline connected to their use has rarely been quantified. Furthermore, no study has employed reference-based assembly using Nanopore reads, which is common with Illumina data. Using Leucanthemum Mill. as an example, we compared the sequence quality of seven chloroplast genome assemblies of the same species, using combinations of two sequencing platforms and three analysis pipelines. In addition, we assessed the factors which might influence Nanopore assembly quality during sequence generation and bioinformatic processing. The consensus sequence derived from de novo assembly of Nanopore data had a sequence identity of 99.59% compared to Illumina short-read de novo assembly. Most of the errors detected were indels (81.5%), and a large majority of them is part of homopolymer regions. The quality of reference-based assembly is heavily dependent upon the choice of a close-enough reference. When using a reference with 0.83% sequence divergence from the studied species, mapping of Nanopore reads results in a consensus comparable to that from Nanopore de novo assembly, and of only slightly inferior quality compared to a reference-based assembly with Illumina data. For optimal de novo assembly of Nanopore data, appropriate filtering of contaminants and chimeric sequences, as well as employing moderate read coverage, is essential. Based on these results, we conclude that Nanopore long reads are a suitable alternative to Illumina short reads in plastome phylogenomics. Few errors remain in the finalized assembly, which can be easily masked in phylogenetic analyses without loss in analytical accuracy. The easily applicable and cost-effective technology might warrant more attention by researchers dealing with plant chloroplast genomes.

New microsatellite markers for Dacryodes edulis (Burseraceae), an indigenous fruit tree species from Central Africa

Microsatellites were designed and characterized in the African fruit tree species Dacryodes edulis (Burseraceae). The fruits are commercialized throughout Central Africa and the species is present in forested environments as well as cultivated systems. The high variability of these markers makes them suitable to investigate the structure of genetic diversity in this important food tree species from Central Africa. From a genomic library obtained by next-generation sequencing, 21 new polymorphic microsatellite loci were developed. Tested on 95 individuals from four populations coming from three countries of the Congo Basin, the microsatellites displayed two to 20 alleles (mean 7.5; expected heterozygosity 0.003 to 0.937, mean 0.666). The transferability of microsatellites was effective for four other Dacryodes species (D. buettneri, D. igaganga, D. osika, D. pubescens). This set of newly developed microsatellite markers will be useful for assessing the genetic diversity and differentiation as well as gene flow patterns of D. edulis in tropical forests from Central Africa.

Skim-Sequencing Reveals the Likely Origin of the Enigmatic Endangered Sunflower Helianthus schweinitzii

Resolving the origin of endangered taxa is an essential component of conservation. This information can be used to guide efforts of bolstering genetic diversity, and also enables species recovery and future evolutionary studies. Here, we used low-coverage whole genome sequencing to clarify the origin of Helianthus schweinitzii, an endangered tetraploid sunflower that is endemic to the Piedmont Plateau in the eastern United States. We surveyed four accessions representing four populations of H. schweinitzii and 38 accessions of six purported parental species. Using de novo approaches, we assembled 87,004 bp of the chloroplast genome and 6770 bp of the nuclear 35S rDNA. Phylogenetic reconstructions based on the chloroplast genome revealed no reciprocal monophyly of taxa. In contrast, nuclear rDNA data strongly supported the currently accepted sections of the genus Helianthus. Information from combined cpDNA and rDNA provided evidence that H. schweinitzii is likely an allo-tetraploid that formed as a result of hybridization between the diploids Helianthus giganteus and Helianthus microcephalus.

Advances in Molecular Genetics and Genomics of African Rice (Oryza glaberrima Steud)

African rice (Oryza glaberrima) has a pool of genes for resistance to diverse biotic and abiotic stresses, making it an important genetic resource for rice improvement. African rice has potential for breeding for climate resilience and adapting rice cultivation to climate change. Over the last decade, there have been tremendous technological and analytical advances in genomics that have dramatically altered the landscape of rice research. Here we review the remarkable advances in knowledge that have been witnessed in the last few years in the area of genetics and genomics of African rice. Advances in cheap DNA sequencing technologies have fuelled development of numerous genomic and transcriptomic resources. Genomics has been pivotal in elucidating the genetic architecture of important traits thereby providing a basis for unlocking important trait variation. Whole genome re-sequencing studies have provided great insights on the domestication process, though key studies continue giving conflicting conclusions and theories. However, the genomic resources of African rice appear to be under-utilized as there seems to be little evidence that these vast resources are being productively exploited for example in practical rice improvement programmes. Challenges in deploying African rice genetic resources in rice improvement and the genomics efforts made in addressing them are highlighted.

Development of high-resolution DNA barcodes for Dioscorea species discrimination and phylogenetic analysis

The genus Dioscorea is widely distributed in tropical and subtropical regions, and is economically important in terms of food supply and pharmaceutical applications. However, DNA barcodes are relatively unsuccessful in discriminating between Dioscorea species, with the highest discrimination rate (23.26%) derived from matK sequences. In this study, we compared genic and intergenic regions of three Dioscorea chloroplast genomes and found that the density of SNPs and indels in intergenic sites was about twice and seven times higher than that of SNPs and indels in the genic regions, respectively. A total of 52 primer pairs covering highly variable regions were designed and seven pairs of primers had 80%-100% PCR success rate. PCR amplicons of 73 Dioscorea individuals and assembled sequences of 47 Dioscorea SRAs were used for estimating intraspecific and interspecific divergence for the seven loci: The rpoB-trnC locus had the highest interspecific divergence. Automatic barcoding gap discovery (ABGD), Poisson tree processes (PTP), and generalized mixed Yule coalescence (GMYC) analysis were applied for species delimitation based on the seven loci and successfully identified the majority of species, except for species in the Enantiophyllum section. Phylogenetic analysis of 51 Dioscorea individuals (28 species) showed that most individuals belonging to the same species tended to cluster in the same group. Our results suggest that the variable loci derived from comparative analysis of plastid genome sequences could be good DNA barcode candidates for taxonomic analysis and species delimitation.

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.

Evolution of Oryza chloroplast genomes promoted adaptation to diverse ecological habitats

The course, tempo and mode of chloroplast genome evolution remain largely unknown, resulting in limited knowledge about how plant plastome gene and genome evolve during the process of recent plant speciation. Here, we report the complete plastomes of 22 closely related Oryza species in chronologically ordered stages and generate the first precise map of genomic structural variation, to our knowledge. The occurrence rapidity was estimated on average to be ~7 insertions and ~15 deletions per Myr. Relatively fewer deletions than insertions result in an increased repeat density that causes the observed growth of Oryza chloroplast genome sizes. Genome-wide scanning identified 14 positively selected genes that are relevant to photosynthesis system, eight of which were found independently in shade-tolerant or sun-loving rice species. psaA seemed positively selected in both shade-tolerant and sun-loving rice species. The results show that adaptive evolution of chloroplast genes makes rice species adapt to diverse ecological habitats related to sunlight preferences.

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.

Sequence capture using AFLP-generated baits: A cost-effective method for high-throughput phylogenetic and phylogeographic analysis

Target sequence capture is an efficient technique to enrich specific genomic regions for high-throughput sequencing in ecological and evolutionary studies. In recent years, many sequence capture approaches have been proposed, but most of them rely on commercial synthetic baits which make the experiment expensive. Here, we present a novel sequence capture approach called AFLP-based genome sequence capture (AFLP Capture). This method uses the AFLP (amplified fragment length polymorphism) technique to generate homemade capture baits without the need for prior genome information, thus is applicable to any organisms. In this approach, biotinylated AFLP fragments representing a random fraction of the genome are used as baits to capture the homologous fragments from genomic shotgun sequencing libraries. In a trial study, by using AFLP Capture, we successfully obtained 511 orthologous loci (>700,000 bp in total length) from 11 Odorrana species and more than 100,000 single nucleotide polymorphisms (SNPs) in four analyzed individuals of an Odorrana species. This result shows that our method can be used to address questions of various evolutionary depths (from interspecies level to intraspecies level). We also discuss the flexibility in bait preparation and how the sequencing data are analyzed. In summary, AFLP Capture is a rapid and flexible tool and can significantly reduce the experimental cost for phylogenetic studies that require analyzing genome-scale data (hundreds or thousands of loci).

Phylogenomics of the Major Tropical Plant Family Annonaceae Using Targeted Enrichment of Nuclear Genes

Targeted enrichment and sequencing of hundreds of nuclear loci for phylogenetic reconstruction is becoming an important tool for plant systematics and evolution. Annonaceae is a major pantropical plant family with 110 genera and ca. 2,450 species, occurring across all major and minor tropical forests of the world. Baits were designed by sequencing the transcriptomes of five species from two of the largest Annonaceae subfamilies. Orthologous loci were identified. The resulting baiting kit was used to reconstruct phylogenetic relationships at two different levels using concatenated and gene tree approaches: a family wide Annonaceae analysis sampling 65 genera and a species level analysis of tribe Piptostigmateae sampling 29 species with multiple individuals per species. DNA extraction was undertaken mainly on silicagel dried leaves, with two samples from herbarium dried leaves. Our kit targets 469 exons (364,653 bp of sequence data), successfully capturing sequences from across Annonaceae. Silicagel dried and herbarium DNA worked equally well. We present for the first time a nuclear gene-based phylogenetic tree at the generic level based on 317 supercontigs. Results mainly confirm previous chloroplast based studies. However, several new relationships are found and discussed. We show significant differences in branch lengths between the two large subfamilies Annonoideae and Malmeoideae. A new tribe, Annickieae, is erected containing a single African genus Annickia. We also reconstructed a well-resolved species-level phylogenetic tree of the Piptostigmteae tribe. Our baiting kit is useful for reconstructing well-supported phylogenetic relationships within Annonaceae at different taxonomic levels. The nuclear genome is mainly concordant with plastome information with a few exceptions. Moreover, we find that substitution rate heterogeneity between the two subfamilies is also found within the nuclear compartment, and not just plastomes and ribosomal DNA as previously shown. Our results have implications for understanding the biogeography, molecular dating and evolution of Annonaceae.

Complete chloroplast genome sequences of Dioscorea: Characterization, genomic resources, and phylogenetic analyses

Dioscorea L., the largest genus of the family Dioscoreaceae with over 600 species, is not only an important food but also a medicinal plant. The identification and classification of Dioscorea L. is a rather difficult task. In this study, we sequenced five Dioscorea chloroplast genomes, and analyzed with four other chloroplast genomes of Dioscorea species from GenBank. The Dioscorea chloroplast genomes displayed the typical quadripartite structure of angiosperms, which consisted of a pair of inverted repeats separated by a large single-copy region, and a small single-copy region. The location and distribution of repeat sequences and microsatellites were determined, and the rapidly evolving chloroplast genome regions (trnK-trnQ, trnS-trnG, trnC-petN, trnE-trnT, petG-trnW-trnP, ndhF, trnL-rpl32, and ycf1) were detected. Phylogenetic relationships of Dioscorea inferred from chloroplast genomes obtained high support even in shortest internodes. Thus, chloroplast genome sequences provide potential molecular markers and genomic resources for phylogeny and species identification.

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.

Reassembling haplotypes in a mixture of pooled amplicons when the relative concentrations are known: A proof-of-concept study on the efficient design of next-generation sequencing strategies

Next-generation sequencing can be costly and labour intensive. Usually, the sequencing cost per sample is reduced by pooling amplified DNA = amplicons) derived from different individuals on the same sequencing lane. Barcodes unique to each amplicon permit short-read sequences to be assigned appropriately. However, the cost of the library preparation increases with the number of barcodes used. We propose an alternative to barcoding: by using different known proportions of individually-derived amplicons in a pooled sample, each is characterised a priori by an expected depth of coverage. We have developed a Hidden Markov Model that uses these expected proportions to reconstruct the input sequences. We apply this method to pools of mitochondrial DNA amplicons extracted from kangaroo meat, genus Macropus. Our experiments indicate that the sequence coverage can be efficiently used to index the short-reads and that we can reassemble the input haplotypes when secondary factors impacting the coverage are controlled. We therefore demonstrate that, by combining our approach with standard barcoding, the cost of the library preparation is reduced to a third.

A 250 plastome phylogeny of the grass family (Poaceae): topological support under different data partitions

The systematics of grasses has advanced through applications of plastome phylogenomics, although studies have been largely limited to subfamilies or other subgroups of Poaceae. Here we present a plastome phylogenomic analysis of 250 complete plastomes (179 genera) sampled from 44 of the 52 tribes of Poaceae. Plastome sequences were determined from high throughput sequencing libraries and the assemblies represent over 28.7 Mbases of sequence data. Phylogenetic signal was characterized in 14 partitions, including (1) complete plastomes; (2) protein coding regions; (3) noncoding regions; and (4) three loci commonly used in single and multi-gene studies of grasses. Each of the four main partitions was further refined, alternatively including or excluding positively selected codons and also the gaps introduced by the alignment. All 76 protein coding plastome loci were found to be predominantly under purifying selection, but specific codons were found to be under positive selection in 65 loci. The loci that have been widely used in multi-gene phylogenetic studies had among the highest proportions of positively selected codons, suggesting caution in the interpretation of these earlier results. Plastome phylogenomic analyses confirmed the backbone topology for Poaceae with maximum bootstrap support (BP). Among the 14 analyses, 82 clades out of 309 resolved were maximally supported in all trees. Analyses of newly sequenced plastomes were in agreement with current classifications. Five of seven partitions in which alignment gaps were removed retrieved Panicoideae as sister to the remaining PACMAD subfamilies. Alternative topologies were recovered in trees from partitions that included alignment gaps. This suggests that ambiguities in aligning these uncertain regions might introduce a false signal. Resolution of these and other critical branch points in the phylogeny of Poaceae will help to better understand the selective forces that drove the radiation of the BOP and PACMAD clades comprising more than 99.9% of grass diversity.

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.

Evolution in the Amphi-Atlantic tropical genus Guibourtia (Fabaceae, Detarioideae), combining NGS phylogeny and morphology

Tropical rain forests support a remarkable diversity of tree species, questioning how and when this diversity arose. The genus Guibourtia (Fabaceae, Detarioideae), characterized by two South American and 13 African tree species growing in various tropical biomes, is an interesting model to address the role of biogeographic processes and adaptation to contrasted environments on species diversification. Combining whole plastid genome sequencing and morphological characters analysis, we studied the timing of speciation and diversification processes in Guibourtia through molecular dating and ancestral habitats reconstruction. All species except G. demeusei and G. copallifera appear monophyletic. Dispersal from Africa to America across the Atlantic Ocean is the most plausible hypothesis to explain the occurrence of Neotropical Guibourtia species, which diverged ca. 11.8 Ma from their closest African relatives. The diversification of the three main clades of African Guibourtia is concomitant to Miocene global climate changes, highlighting pre-Quaternary speciation events. These clades differ by their reproductive characters, which validates the three subgenera previously described: Pseudocopaiva, Guibourtia and Gorskia. Within most monophyletic species, plastid lineages start diverging from each other during the Pliocene or early Pleistocene, suggesting that these species already arose during this period. The multiple transitions between rain forests and dry forests/savannahs inferred here through the plastid phylogeny in each Guibourtia subgenus address thus new questions about the role of phylogenetic relationships in shaping ecological niche and morphological similarity among taxa.

The complete chloroplast genome of Primulina and two novel strategies for development of high polymorphic loci for population genetic and phylogenetic studies

Background

Primulina Hance is an emerging model for studying evolutionary divergence, adaptation and speciation of the karst flora. However, phylogenetic relationships within the genus have not been resolved due to low variation detected in the cpDNA regions. Chloroplast genomes can provide important information for phylogenetic and population genetic studies. Recent advances in next-generation sequencing (NGS) techniques greatly facilitate sequencing whole chloroplast genomes for multiple individuals. Consequently, novel strategies for development of highly polymorphic loci for population genetic and phylogenetic studies based on NGS data are needed.

Methods

For development of high polymorphic loci for population genetic and phylogenetic studies, two novel strategies are proposed here. The first protocol develops lineage-specific highly variable markers from the true high variation regions (Con_Seas) across whole cp genomes, instead of traditional noncoding regions. The pipeline has been integrated into a single perl script, and named "Con_Sea_Identification_and_PIC_Calculation". The second method assembles chloroplast fragments (poTs) and sub-super-marker (CpContigs) through our "SACRing" pipeline. This approach can fundamentally alter the strategies used in phylogenetic and population genetic studies based on cp markers, facilitating a transition from traditional Sanger sequencing to RAD-Seq. Both of these scripts are available at https://github.com/scbgfengchao/.

Results

Three complete Primulina chloroplast genomes were assembled from genome survey data, and then two novel strategies were developed to yield highly polymorphic markers. For experimental evaluation of the first protocol, a set of Primulina species were used for PCR amplification. The results showed that these newly developed markers are more variable than traditional ones, and seem to be a better choice for phylogenetic and population studies in Primulina. The second method was also successfully applied in population genetic studies of 21 individuals from three natural populations of Primulina.

Conclusions

These two novel strategies may provide a pathway for similar research in other non-model species. The newly developed high polymorphic loci in this study will promote further the phylogenetic and population genetic studies in Primulina and other genera of the family Gesneriaceae.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1067-z) contains supplementary material, which is available to authorized users.

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.

Molecular basis of African yam domestication: analyses of selection point to root development, starch biosynthesis, and photosynthesis related genes

BACKGROUND

After cereals, root and tuber crops are the main source of starch in the human diet. Starch biosynthesis was certainly a significant target for selection during the domestication of these crops. But domestication of these root and tubers crops is also associated with gigantism of storage organs and changes of habitat.

RESULTS

We studied here, the molecular basis of domestication in African yam, Dioscorea rotundata. The genomic diversity in the cultivated species is roughly 30% less important than its wild relatives. Two percent of all the genes studied showed evidences of selection. Two genes associated with the earliest stages of starch biosynthesis and storage, the sucrose synthase 4 and the sucrose-phosphate synthase 1 showed evidence of selection. An adventitious root development gene, a SCARECROW-LIKE gene was also selected during yam domestication. Significant selection for genes associated with photosynthesis and phototropism were associated with wild to cultivated change of habitat. If the wild species grow as vines in the shade of their tree tutors, cultivated yam grows in full light in open fields.

CONCLUSIONS

Major rewiring of aerial development and adaptation for efficient photosynthesis in full light characterized yam domestication.