The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Low-frequency somatic mutations are heritable in tropical trees Dicorynia guianensis and Sextonia rubra

Somatic mutations potentially play a role in plant evolution, but common expectations pertaining to plant somatic mutations remain insufficiently tested. Unlike in most animals, the plant germline is assumed to be set aside late in development, leading to the expectation that plants accumulate somatic mutations along growth. Therefore, several predictions were made on the fate of somatic mutations: mutations have generally low frequency in plant tissues; mutations at high frequency have a higher chance of intergenerational transmission; branching topology of the tree dictates mutation distribution; and exposure to UV (ultraviolet) radiation increases mutagenesis. To provide insights into mutation accumulation and transmission in plants, we produced two high-quality reference genomes and a unique dataset of 60 high-coverage whole-genome sequences of two tropical tree species, Dicorynia guianensis (Fabaceae) and Sextonia rubra (Lauraceae). We identified 15,066 de novo somatic mutations in D. guianensis and 3,208 in S. rubra, surprisingly almost all found at low frequency. We demonstrate that 1) low-frequency mutations can be transmitted to the next generation; 2) mutation phylogenies deviate from the branching topology of the tree; and 3) mutation rates and mutation spectra are not demonstrably affected by differences in UV exposure. Altogether, our results suggest far more complex links between plant growth, aging, UV exposure, and mutation rates than commonly thought.

Microbial biomarkers of tree water status for next-generation biomonitoring of forest ecosystems

Next-generation biomonitoring proposes to combine machine-learning algorithms with environmental DNA data to automate the monitoring of the Earth's major ecosystems. In the present study, we searched for molecular biomarkers of tree water status to develop next-generation biomonitoring of forest ecosystems. Because phyllosphere microbial communities respond to both tree physiology and climate change, we investigated whether environmental DNA data from tree phyllosphere could be used as molecular biomarkers of tree water status in forest ecosystems. Using an amplicon sequencing approach, we analysed phyllosphere microbial communities of four tree species (Quercus ilex, Quercus robur, Pinus pinaster and Betula pendula) in a forest experiment composed of irrigated and non-irrigated plots. We used these microbial community data to train a machine-learning algorithm (Random Forest) to classify irrigated and non-irrigated trees. The Random Forest algorithm detected tree water status from phyllosphere microbial community composition with more than 90% accuracy for oak species, and more than 75% for pine and birch. Phyllosphere fungal communities were more informative than phyllosphere bacterial communities in all tree species. Seven fungal amplicon sequence variants were identified as candidates for the development of molecular biomarkers of water status in oak trees. Altogether, our results show that microbial community data from tree phyllosphere provides information on tree water status in forest ecosystems and could be included in next-generation biomonitoring programmes that would use in situ, real-time sequencing of environmental DNA to help monitor the health of European temperate forest ecosystems.

Comparative genomics of Mycoplasma feriruminatoris, a fast-growing pathogen of wild Caprinae

Mycoplasma feriruminatoris is a fast-growing Mycoplasma species isolated from wild Caprinae and first described in 2013. M. feriruminatoris isolates have been associated with arthritis, kerato conjunctivitis, pneumonia and septicemia, but were also recovered from apparently healthy animals. To better understand what defines this species, we performed a genomic survey on 14 strains collected from free-ranging or zoo-housed animals between 1987 and 2017, mostly in Europe. The average chromosome size of the M. feriruminatoris strains was 1,040±0,024 kbp, with 24 % G+C and 852±31 CDS. The core genome and pan-genome of the M. feriruminatoris species contained 628 and 1312 protein families, respectively. The M. feriruminatoris strains displayed a relatively closed pan-genome, with many features and putative virulence factors shared with species from the M. mycoides cluster, including the MIB-MIP Ig cleavage system, a repertoire of DUF285 surface proteins and a complete biosynthetic pathway for galactan. M. feriruminatoris genomes were found to be mostly syntenic, although repertoires of mobile genetic elements, including Mycoplasma Integrative and Conjugative Elements, insertion sequences, and a single plasmid varied. Phylogenetic- and gene content analyses confirmed that M. feriruminatoris was closer to the M. mycoides cluster than to the ruminant species M. yeatsii and M. putrefaciens. Ancestral genome reconstruction showed that the emergence of the M. feriruminatoris species was associated with the gain of 17 gene families, some of which encode defence enzymes and surface proteins, and the loss of 25 others, some of which are involved in sugar transport and metabolism. This comparative study suggests that the M. mycoides cluster could be extended to include M. feriruminatoris. We also find evidence that the specific organization and structure of the DnaA boxes around the oriC of M. feriruminatoris may contribute to drive the remarkable fast growth of this minimal bacterium.

Effects of azithromycin and doxycycline on the vaginal microbiota of women with urogenital Chlamydia trachomatis infection: a substudy of the Chlazidoxy randomized controlled trial

OBJECTIVES

Dysbiotic bacterial communities within the vagina are associated with Chlamydia trachomatis infection. We compared the effect of treatment with azithromycin and doxycycline on the vaginal microbiota in a cohort of women with a urogenital C. trachomatis infection randomly assigned to one of these treatments (Chlazidoxy trial).

METHODS

We analysed vaginal samples from 284 women (135 in the azithromycin group and 149 in the doxycycline group) collected at baseline and 6 weeks after treatment initiation. The vaginal microbiota was characterized using 16S rRNA gene sequencing and classified into community state types (CSTs).

RESULTS

At baseline, 75% (212/284) of the women had a high-risk microbiota (CST-III or CST-IV). A cross-sectional comparison 6 weeks after treatment showed that 15 phylotypes were differentially abundant, but this difference was not reflected at the CST (p=0.772) or diversity level (p=0.339). Between baseline and the 6-week visit, α-diversity (p=0.140) and transition probabilities between CSTs were not significantly different between the groups, and no phylotype was differentially abundant.

CONCLUSION

In women with urogenital C. trachomatis infection, the vaginal microbiota does not seem to be affected by azithromycin or doxycycline 6 weeks after treatment. Because the vaginal microbiota remains susceptible to C. trachomatis infection (with CST-III or CST-IV) following antibiotic treatment, women remain at risk of reinfection, which could originate from unprotected sexual intercourse or untreated anorectal C. trachomatis infection. This last consideration advocates for the use of doxycycline instead of azithromycin because of its higher anorectal microbiological cure rate.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

Fruit fly phylogeny imprints bacterial gut microbiota

One promising avenue for reconciling the goals of crop production and ecosystem preservation consists in the manipulation of beneficial biotic interactions, such as between insects and microbes. Insect gut microbiota can affect host fitness by contributing to development, host immunity, nutrition, or behavior. However, the determinants of gut microbiota composition and structure, including host phylogeny and host ecology, remain poorly known. Here, we used a well‐studied community of eight sympatric fruit fly species to test the contributions of fly phylogeny, fly specialization, and fly sampling environment on the composition and structure of bacterial gut microbiota. Comprising both specialists and generalists, these species belong to five genera from to two tribes of the Tephritidae family. For each fly species, one field and one laboratory samples were studied. Bacterial inventories to the genus level were produced using 16S metabarcoding with the Oxford Nanopore Technology. Sample bacterial compositions were analyzed with recent network‐based clustering techniques. Whereas gut microbiota were dominated by the Enterobacteriaceae family in all samples, microbial profiles varied across samples, mainly in relation to fly identity and sampling environment. Alpha diversity varied across samples and was higher in the Dacinae tribe than in the Ceratitinae tribe. Network analyses allowed grouping samples according to their microbial profiles. The resulting groups were very congruent with fly phylogeny, with a significant modulation of sampling environment, and with a very low impact of fly specialization. Such a strong imprint of host phylogeny in sympatric fly species, some of which share much of their host plants, suggests important control of fruit flies on their gut microbiota through vertical transmission and/or intense filtering of environmental bacteria.

Gut Microbiota and Mycobiota Evolution Is Linked to Memory Improvement after Bariatric Surgery in Obese Patients: A Pilot Study

Patients with obesity are known to exhibit gut microbiota dysbiosis and memory deficits. Bariatric surgery (BS) is currently the most efficient anti-obesity treatment and may improve both gut dysbiosis and cognition. However, no study has investigated association between changes of gut microbiota and cognitive function after BS. We prospectively evaluated 13 obese patients on anthropometric data, memory functions, and gut microbiota-mycobiota before and six months after BS. The Rey Auditory Verbal Learning Test (AVLT) and the symbol span (SS) of the Weschler Memory Scale were used to assess verbal and working memory, respectively. Fecal microbiota and mycobiota were longitudinally analyzed by 16S and ITS2 rRNA sequencing respectively. AVLT and SS scores were significantly improved after BS (AVLT scores: 9.7 ± 1.7 vs. 11.2 ± 1.9, p = 0.02, and SS scores: 9.7 ± 23.0 vs. 11.6 ± 2.9, p = 0.05). An increase in bacterial alpha-diversity, and Ruminococcaceae, Prevotella, Agaricus, Rhodotorula, Dipodascus, Malassezia, and Mucor were significantly associated with AVLT score improvement after BS, while an increase in Prevotella and a decrease in Clostridium, Akkermansia, Dipodascus and Candida were linked to SS scores improvement. We identified several changes in the microbial communities that differ according to the improvement of either the verbal or working memories, suggesting a complex gut-brain-axis that evolves after BS.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

Development of nuclear SNP markers for Mahogany (Swietenia spp.)

Swietenia species are the most valuable American tropical timbers and have been heavily overexploited for decades. The three species are listed as either vulnerable or endangered by IUCN and are included on Appendix II of CITES, yet illegal exploitation continues. Here, we used restriction associated DNA sequencing to develop a new set of 120 SNP markers for Swietenia sp., suitable for MassARRAY®iPLEX™ genotyping. These markers can be used for population genetic studies and timber tracking purposes.

Fast sequence-based microsatellite genotyping development workflow

Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20-40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.

Genetic Characterization of Cupped Oyster Resources in Europe Using Informative Single Nucleotide Polymorphism (SNP) Panels

The Pacific oyster, Crassostrea gigas, was voluntarily introduced from Japan and British Columbia into Europe in the early 1970s, mainly to replace the Portuguese oyster, Crassostrea angulata, in the French shellfish industry, following a severe disease outbreak. Since then, the two species have been in contact in southern Europe and, therefore, have the potential to exchange genes. Recent evolutionary genomic works have provided empirical evidence that C. gigas and C. angulata exhibit partial reproductive isolation. Although hybridization occurs in nature, the rate of interspecific gene flow varies across the genome, resulting in highly heterogeneous genome divergence. Taking this biological property into account is important to characterize genetic ancestry and population structure in oysters. Here, we identified a subset of ancestry-informative makers from the most differentiated regions of the genome using existing genomic resources. We developed two different panels in order to (i) easily differentiate C. gigas and C. angulata, and (ii) describe the genetic diversity and structure of the cupped oyster with a particular focus on French Atlantic populations. Our results confirm high genetic homogeneity among Pacific cupped oyster populations in France and reveal several cases of introgressions between Portuguese and Japanese oysters in France and Portugal.

Development of new SNPs loci on Quercus robur and Quercus petraea for genetic studies covering the whole species’ distribution range

We used double digest restriction site associated DNA sequencing (ddRAD) to develop new geographically informative nuclear SNP loci in Quercus robur and Quercus petraea. Genotypes derived from sequence data of 95 individuals covering the distribution range of the species were analysed to select geographically informative and polymorphic loci within Russia and Germany. We successfully screened a selected set of 119 loci on a MassARRAY® iPLEX™ platform on 190 individuals from 19 locations in Russia. The newly developed loci will be useful for genetic studies over the whole distribution range of both species.

Performance of semiconductor sequencing platform for non-invasive prenatal genetic screening for fetal aneuploidy: results from a multicenter prospective cohort study in a clinical setting

OBJECTIVE

To validate and evaluate the performance metrics of the high-throughput semiconductor sequencing platform, Ion Proton®, in non-invasive prenatal genetic screening (NIPS) for common fetal aneuploidies in a clinical setting.

METHODS

This prospective cohort study included 2505 pregnant women from eight academic genetics laboratories (695 high risk for trisomy 21 (risk ≥ 1/250) pregnancies in a validation study, and 1810 such pregnancies, without ultrasound anomalies, in a real-life NIPS clinical setting). Outcome was available for all cases in the validation cohort and for 521 in the clinical cohort. Cell-free DNA from plasma samples was sequenced using the Ion Proton sequencer, and sequencing data were analyzed using the open-access software, WISECONDOR. Performance metrics for detection of trisomies 21, 18 and 13 were calculated based on either fetal karyotype result or clinical data collected at birth. We also evaluated the failure rate and compared three methods of fetal fraction quantification (RASSF1A assay, and DEFRAG and SANEFALCON software).

RESULTS

Results from both cohorts were consistent and their gestational age was not significantly different so their data were combined to increase the sample size for analysis. Sensitivities and specificities, respectively, were as follows: for trisomy 21, 98.3% (95% CI, 93.5-99.7%) and 99.9% (95% CI, 99.4-100%); for trisomy 18, 96.7% (95% CI, 80.9-99.8%) and 100% (95% CI, 99.6-100%); and for trisomy 13, 94.1% (95% CI, 69.2-99.7%) and 100% (95% CI, 99.6-100%). Our failure rate was 1.2% initially and as low as 0.6% after retesting some of the failed samples. Fetal fraction estimation by the RASSF1A assay was consistent with DEFRAG results, and both were adequate for routine diagnosis.

CONCLUSIONS

We describe one of the largest studies evaluating Ion Proton-based NIPS and the first clinical study reporting pregnancy outcome in a large series of patients. This platform is highly efficient in detecting the three most common trisomies. Our protocol is robust and can be implemented easily in any medical genetics laboratory. Copyright © 2018 ISUOG. Published by John Wiley & Sons Ltd.

A few north Appalachian populations are the source of European black locust

The role of evolution in biological invasion studies is often overlooked. In order to evaluate the evolutionary mechanisms behind invasiveness, it is crucial to identify the source populations of the introduction. Studies in population genetics were carried out on Robinia pseudoacacia L., a North American tree which is now one of the worst invasive tree species in Europe. We realized large-scale sampling in both the invasive and native ranges: 63 populations were sampled and 818 individuals were genotyped using 113 SNPs. We identified clonal genotypes in each population and analyzed between and within range population structure, and then, we compared genetic diversity between ranges, enlarging the number of SNPs to mitigate the ascertainment bias. First, we demonstrated that European black locust was introduced from just a limited number of populations located in the Appalachian Mountains, which is in agreement with the historical documents briefly reviewed in this study. Within America, population structure reflected the effects of long-term processes, whereas in Europe it was largely impacted by human activities. Second, we showed that there is a genetic bottleneck between the ranges with a decrease in allelic richness and total number of alleles in Europe. Lastly, we found more clonality within European populations. Black locust became invasive in Europe despite being introduced from a reduced part of its native distribution. Our results suggest that human activity, such as breeding programs in Europe and the seed trade throughout the introduced range, had a major role in promoting invasion; therefore, the introduction of the missing American genetic cluster to Europe should be avoided.

The use and domestication of Theobroma cacao during the mid-Holocene in the upper Amazon

Cacao (Theobroma cacao L.) is an important economic crop, yet studies of its domestication history and early uses are limited. Traditionally, cacao is thought to have been first domesticated in Mesoamerica. However, genomic research shows that T. cacao's greatest diversity is in the upper Amazon region of northwest South America, pointing to this region as its centre of origin. Here, we report cacao use identified by three independent lines of archaeological evidence-cacao starch grains, absorbed theobromine residues and ancient DNA-dating from approximately 5,300 years ago recovered from the Santa Ana-La Florida (SALF) site in southeast Ecuador. To our knowledge, these findings constitute the earliest evidence of T. cacao use in the Americas and the first unequivocal archaeological example of its pre-Columbian use in South America. They also reveal the upper Amazon region as the oldest centre of cacao domestication yet identified.

Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting

DNA barcoding has proved difficult in a number of woody plant genera, including the ecologically important oak genus Quercus. In this study, we utilized restrictionsite-associated DNA sequencing (RAD-seq) to develop an economical single nucleotide polymorphism (SNP) DNA barcoding system that suffices to distinguish eight common, sympatric eastern North American white oak species. Two de novo clustering pipelines, PyRAD and Stacks, were used in combination with postclustering bioinformatic tools to generate a list of 291 potential SNPs, 80 of which were included in a barcoding toolkit that is easily implemented using MassARRAY mass spectrometry technology. As a proof-of-concept, we used the genotyping toolkit to infer potential hybridization between North American white oaks transplanted outside of their native range (Q. michauxii, Q. montana, Q muehlenbergii/Q. prinoides, and Q. stellata) into a horticultural collection surrounded by natural forests of locally native trees (Q. alba and Q. macrocarpa) in the living collection at The Morton Arboretum (Lisle, IL, USA). Phylogenetic and clustering analyses suggested low rates of hybridization between cultivated and native species, with the exception of one Q. michauxii mother tree, the acorns of which exhibited high admixture from either Q. alba or Q. stellata and Q. macrocarpa, and a hybrid between Q. stellata that appears to have backcrossed almost exclusively to Q. alba. Together, RAD-seq and MassARRAY technologies allow for efficient development and implementation of a multispecies barcode for one of the more challenging forest tree genera.

High-Throughput DNA sequencing of ancient wood

Reconstructing the colonization and demographic dynamics that gave rise to extant forests is essential to forecasts of forest responses to environmental changes. Classical approaches to map how population of trees changed through space and time largely rely on pollen distribution patterns, with only a limited number of studies exploiting DNA molecules preserved in wooden tree archaeological and subfossil remains. Here, we advance such analyses by applying high-throughput (HTS) DNA sequencing to wood archaeological and subfossil material for the first time, using a comprehensive sample of 167 European white oak waterlogged remains spanning a large temporal (from 550 to 9,800 years) and geographical range across Europe. The successful characterization of the endogenous DNA and exogenous microbial DNA of 140 (~83%) samples helped the identification of environmental conditions favouring long-term DNA preservation in wood remains, and started to unveil the first trends in the DNA decay process in wood material. Additionally, the maternally inherited chloroplast haplotypes of 21 samples from three periods of forest human-induced use (Neolithic, Bronze Age and Middle Ages) were found to be consistent with those of modern populations growing in the same geographic areas. Our work paves the way for further studies aiming at using ancient DNA preserved in wood to reconstruct the micro-evolutionary response of trees to climate change and human forest management.

Dealing with paralogy in RADseq data: in silico detection and single nucleotide polymorphism validation in Robinia pseudoacacia L

The RADseq technology allows researchers to efficiently develop thousands of polymorphic loci across multiple individuals with little or no prior information on the genome. However, many questions remain about the biases inherent to this technology. Notably, sequence misalignments arising from paralogy may affect the development of single nucleotide polymorphism (SNP) markers and the estimation of genetic diversity. We evaluated the impact of putative paralog loci on genetic diversity estimation during the development of SNPs from a RADseq dataset for the nonmodel tree species Robinia pseudoacacia L. We sequenced nine genotypes and analyzed the frequency of putative paralogous RAD loci as a function of both the depth of coverage and the mismatch threshold allowed between loci. Putative paralogy was detected in a very variable number of loci, from 1% to more than 20%, with the depth of coverage having a major influence on the result. Putative paralogy artificially increased the observed degree of polymorphism and resulting estimates of diversity. The choice of the depth of coverage also affected diversity estimation and SNP validation: A low threshold decreased the chances of detecting minor alleles while a high threshold increased allelic dropout. SNP validation was better for the low threshold (4×) than for the high threshold (18×) we tested. Using the strategy developed here, we were able to validate more than 80% of the SNPs tested by means of individual genotyping, resulting in a readily usable set of 330 SNPs, suitable for use in population genetics applications.

Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late-successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.

Fine-scale environmental control of hybridization in oaks

Natural hybridization is attracting much interest in modern speciation and conservation biology studies, but the underlying mechanisms remain poorly understood. In particular, it is unclear why environmental changes often increase hybridization rates. To study this question, we surveyed mating events in a mixed oak stand and developed a spatially explicit individual-based hybridization model. This model, where hybridization is frequency-dependent, pollen is nonlimiting and which allows immigrant pollen to compete with local pollen, takes into account species-specific pollen dispersal and sexual barriers to hybridization. The consequences of pollen limitation on hybridization were studied using another simple model. The results indicate that environmental changes could increase hybridization rates through two distinct mechanisms. First, by disrupting the spatial organization of communities, they should decrease the proportion of conspecific pollen available for mating, thus increasing hybridization rates. Second, by decreasing the density of conspecifics, they should increase pollen limitation and thus hybridization rates, as a consequence of chance pollination predominating over deterministic pollen competition. Altogether, our results point to a need for considering hybridization events at the appropriate level of organization and provide new insights into why hybridization rates generally increase in disturbed environments.

Current trends in microsatellite genotyping

Microsatellites have been popular molecular markers ever since their advent in the late eighties. Despite growing competition from new genotyping and sequencing techniques, the use of these versatile and cost-effective markers continues to increase, boosted by successive technical advances. First, methods for multiplexing PCR have considerably improved over the last years, thereby decreasing genotyping costs and increasing throughput. Second, next-generation sequencing technologies allow the identification of large numbers of microsatellite loci at reduced cost in non-model species. As a consequence, more stringent selection of loci is possible, thereby further enhancing multiplex quality and efficiency. However, current practices are lagging behind. By surveying recently published population genetic studies relying on simple sequence repeats, we show that more than half of the studies lack appropriate quality controls and do not make use of multiplex PCR. To make the most of the latest technical developments, we outline the need for a well-established strategy including standardized high-throughput bench protocols and specific bioinformatic tools, from primer design to allele calling.

Two highly validated multiplexes (12-plex and 8-plex) for species delimitation and parentage analysis in oaks (Quercus spp.)

Multiplex PCR is a fast and cost-effective technique allowing increased genotyping throughput of microsatellites. We developed two multiplexes for Quercus petraea and Q. robur, a 12-plex of EST-SSRs (eSSRs) and an 8-plex of genomic SSRs (gSSRs). We studied the origin of allele calling errors at the human reader and software levels. We showed that the robustness of allele identification can be improved by binning on raw peak sizes prior to genetic data analysis. We checked through simulation the power of these markers for species delimitation and hybrid detection. The resolution achieved with all 20 markers was greatly improved compared to that of previous studies based on a subset of the markers. Preliminary PCR tests suggest that these multiplexes might be useful to study other oak species as well. The strategy used for multiplex microsatellite development (from PCR conditions to the definition of allele calling rules) should be broadly applicable.

A fast and cost-effective approach to develop and map EST-SSR markers: oak as a case study

Background

Expressed Sequence Tags (ESTs) are a source of simple sequence repeats (SSRs) that can be used to develop molecular markers for genetic studies. The availability of ESTs for Quercus robur and Quercus petraea provided a unique opportunity to develop microsatellite markers to accelerate research aimed at studying adaptation of these long-lived species to their environment. As a first step toward the construction of a SSR-based linkage map of oak for quantitative trait locus (QTL) mapping, we describe the mining and survey of EST-SSRs as well as a fast and cost-effective approach (bin mapping) to assign these markers to an approximate map position. We also compared the level of polymorphism between genomic and EST-derived SSRs and address the transferability of EST-SSRs in Castanea sativa (chestnut).

Results

A catalogue of 103,000 Sanger ESTs was assembled into 28,024 unigenes from which 18.6% presented one or more SSR motifs. More than 42% of these SSRs corresponded to trinucleotides. Primer pairs were designed for 748 putative unigenes. Overall 37.7% (283) were found to amplify a single polymorphic locus in a reference full-sib pedigree of Quercus robur. The usefulness of these loci for establishing a genetic map was assessed using a bin mapping approach. Bin maps were constructed for the male and female parental tree for which framework linkage maps based on AFLP markers were available. The bin set consisting of 14 highly informative offspring selected based on the number and position of crossover sites. The female and male maps comprised 44 and 37 bins, with an average bin length of 16.5 cM and 20.99 cM, respectively. A total of 256 EST-SSRs were assigned to bins and their map position was further validated by linkage mapping. EST-SSRs were found to be less polymorphic than genomic SSRs, but their transferability rate to chestnut, a phylogenetically related species to oak, was higher.

Conclusion

We have generated a bin map for oak comprising 256 EST-SSRs. This resource constitutes a first step toward the establishment of a gene-based map for this genus that will facilitate the dissection of QTLs affecting complex traits of ecological importance.

Species relative abundance and direction of introgression in oaks

Successful hybridisation and subsequent introgression lead to the transfer of genetic material across species boundaries. In this process, species relative abundance can play a significant role. If one species is less abundant than the other, its females will receive many heterospecific gametes, increasing mate-recognition errors and thus hybridisation rate. Moreover, first-generation hybrids will also more likely mate with the more abundant species, leading to asymmetric introgression. These predictions have important fundamental consequences, especially during biological invasions or when a rare species threatened by extinction is surrounded by individuals from a related species. However, experimental tests in nature of the importance of the relative abundance of each species on hybridisation dynamics remain scarce. We assess here the impact of species relative abundance on hybridisation dynamics among four species from the European white oak species complex. A total of 2107 oak trees were genotyped at 10 microsatellite markers and Bayesian clustering methods were used to identify reference trees of each species. We then used these reference trees to simulate purebred and hybrid genotypes to determine optimal threshold for genetic assignment. With this approach, we found widespread evidence of hybridisation between all studied oak species, with high occurrence of hybrids, varying from 11% to 31% according to stand and sampling strategies. This finding suggests that hybridisation is a common phenomenon that plays a significant role in evolution of this oak species complex. In addition, we demonstrate a strong impact of species abundance on both hybridisation rate and introgression directionality.