The homomorphic self-incompatibility system in Oleaceae is controlled by a hemizygous genomic region expressing a gibberellin pathway gene

In flowering plants, outcrossing is commonly ensured by self-incompatibility (SI) systems. These can be homomorphic (typically with many different allelic specificities) or can accompany flower heteromorphism (mostly with just two specificities and corresponding floral types). The SI system of the Oleaceae family is unusual, with the long-term maintenance of only two specificities but often without flower morphology differences. To elucidate the genomic architecture and molecular basis of this SI system, we obtained chromosome-scale genome assemblies of Phillyrea angustifolia individuals and related them to a genetic map. The S-locus region proved to have a segregating 543-kb indel unique to one specificity, suggesting a hemizygous region, as observed in all distylous systems so far studied at the genomic level. Only one of the predicted genes in this indel region is found in the olive tree, Olea europaea, genome, also within a segregating indel. We describe complete association between the presence/absence of this gene and the SI types determined for individuals of seven distantly related Oleaceae species. This gene is predicted to be involved in catabolism of the gibberellic acid (GA) hormone, and experimental manipulation of GA levels in developing buds modified the male and female SI responses of the two specificities in different ways. Our results provide a unique example of a homomorphic SI system, where a single conserved gibberellin-related gene in a hemizygous indel underlies the long-term maintenance of two groups of reproductive compatibility.

Whole-genome sequencing and comparative genomics reveal candidate genes associated with quality traits in Dioscorea alata

BACKGROUND

Quality traits are essential determinants of consumer preferences. Dioscorea alata (Greater Yam), is a starchy tuber crop in tropical regions. However, a comprehensive understanding of the genetic basis underlying yam tuber quality remains elusive. To address this knowledge gap, we employed population genomics and candidate gene association approaches to unravel the genetic factors influencing the quality attributes of boiled yam.

METHODS AND RESULTS

Comparative genomics analysis of 45 plant species revealed numerous novel genes absent in the existing D. alata gene annotation. This approach, adding 48% more genes, significantly enhanced the functional annotation of three crucial metabolic pathways associated with boiled yam quality traits: pentose and glucuronate interconversions, starch and sucrose metabolism, and flavonoid biosynthesis. In addition, the whole-genome sequencing of 127 genotypes identified 27 genes under selection and 22 genes linked to texture, starch content, and color through a candidate gene association analysis. Notably, five genes involved in starch content and cell wall composition, including 1,3-beta Glucan synthase, β-amylase, and Pectin methyl esterase, were common to both approaches and their expression levels were assessed by transcriptomic data.

CONCLUSIONS

The analysis of the whole-genome of 127 genotypes of D. alata and the study of three specific pathways allowed the identification of important genes for tuber quality. Our findings provide insights into the genetic basis of yam quality traits and will help the enhancement of yam tuber quality through breeding programs.

Socio-hydrogeological survey and assessment of organic pollutants to highlight and trace back pollution fluxes threatening a coastal groundwater-dependent ecosystem

Despite being a vector of pollution towards connected ecosystems, groundwater is often underestimated or not taken into account in management frameworks. To fill this gap, we propose to add socio-economic data to hydrogeological investigations to identify past and present pollution sources linked to human activities at watershed scale in order to forecast threats towards groundwater-dependent ecosystems (GDEs). The aim of this paper is to demonstrate, by a cross-disciplinary approach, the added value of socio-hydrogeological investigations to tackle anthropogenic pollution fluxes towards a GDE and to contribute to more sustainable management of groundwater resources. A survey combining chemical compounds analysis, data compilation, land use analysis and field investigations with a questionnaire was carried out on the Biguglia lagoon plain (France). Results show a pollution with a two-fold source, both agricultural and domestic, in all water bodies of the plain. The pesticide analysis reveals the presence of 10 molecules, including domestic compounds, with concentrations exceeding European groundwater quality standards for individual pesticides, as well as pesticides already banned for twenty years. On the basis of both the field survey and the questionnaire, agricultural pollution has been identified as very local highlighting the storage capacity of the aquifer, whereas domestic pollution is diffuse over the plain and attributed to sewage network effluents and septic tanks. Domestic compounds present shorter residence time within the aquifer highlighting continuous inputs, related to consumption habits of the population. Under the Water Framework Directive (WFD), member states are required to preserve the good ecological status, quality and quantity of water bodies. However, for GDEs it is difficult to achieve the 'good status' required without considering the groundwater's pollutant storage capacity and pollution legacy. To help resolve this issue, socio-hydrogeology has proved to be an efficient tool as well as for implementing effective protection measures for Mediterranean GDEs.

The role of groundwater in CO2 production and carbon storage in Mediterranean peatlands: An isotope geochemistry approach

Peatlands are permanent wetlands recognized for ecosystem services such as biodiversity conservation and carbon storage capacity. Little information is available about their response to global change, the reason why most Earth system climate models consider a linear increase in the release of greenhouse gases (GHG), such as CO₂, with increasing temperatures. Nevertheless, numerous studies suggest that an increase in the temperature may not imply a decrease in photosynthesis and carbon storage rates if water availability is sufficient, the latter being under the control of local hydrology mechanisms. Mediterranean peatlands well illustrate this fact. Since they are groundwater-dependent, they are hydrologically resilient to the strong seasonality of hydroclimatic conditions, especially during the summer drought. In the present study, we demonstrate that, even if such peatlands release CO₂ into the atmosphere, they can maintain a carbon storage capacity. To this end, a geochemical study disentangles the origin and fate of carbon within a Mediterranean peatland at the watershed scale. Field parameters, major ions, dissolved organic and inorganic carbon content and associated δ¹³C values allow for characterizing the seasonality of hydrochemical mechanisms and carbon input from an alluvial aquifer (where rain, river, shallow, and deep groundwater flows are mixing) to the peatland. The inorganic and organic content of peat soil and δ¹³C values of total organic matter and CO₂ complete the dataset, making it possible to provide arguments in favour of lower organic matter oxidation compared to primary production. Overall, this study highlights the groundwater role in the fluxes of CO₂ at the peatland-atmosphere interface, and more broadly the need to understand the interactions between the water and carbon cycles to build better models of the future evolution of the global climate.

Exploring disease perception in Behçet's syndrome: combining a quantitative and a qualitative study based on a narrative medicine approach

BACKGROUND

Behçet Syndrome (BS) has a significant psychological and social impact on patients, caregivers and families. The present study aims at exploring disease perception in BS patients, using both a co-designed survey and the narrative medicine (NM) approach.

METHODS

An ad-hoc questionnaire was co-designed by clinicians expert in BS, BS patients and caregivers and BS adult patients were invited to answer the online questionnaires. Cluster analysis was used to analyse data from the survey and to identify groups of patients with diverse disease perception. To further explore real-life perspectives, the stories of illness of a smaller group of adult BS patients were anonymously collected online and analysed by means of text, sentiment and qualitative analysis.

RESULTS

Two hundred and seven patients answered the survey and forty-three stories were collected. The cluster analysis highlighted that accepting or not the disease has a strong impact on the daily life, on how BS patients perceive themselves and in terms of hope for the future. The stories revealed that patients often address common issues, such as the long and complex journey faced from the disease onset until the BS diagnosis, which was strongly connected to the concept of time and perceived as an exhausting period of their lives.

CONCLUSION

To our knowledge, this is the first study that addressed disease perception also applying the NM principles in BS. The current perception that BS patients have of their disease should encourage the BS scientific and patient community in joining forces in order to improve the journey of BS patients.

From genotype to phenotype: Genetic redundancy and the maintenance of an adaptive polymorphism in the context of high gene flow

A central question in evolution is how several adaptive phenotypes are maintained within a species. Theory predicts that the genetic determination of a trait, and in particular the amounts of redundancy in the mapping of genotypes to phenotypes, mediates evolutionary outcomes of phenotypic selection. In Mediterranean wild thyme, numerous discrete chemical phenotypes (chemotypes) occur in close geographic proximity. Chemotypes are defined by the predominant monoterpene produced by individual plants in their essential oil. In this study, we analyze the ecological genetics of six chemotypes nested within two well‐established chemical families (hereafter ecotypes). Ecotypes, and chemotypes within ecotypes, are spatially segregated, and their distributions track local differences in the abiotic environment. By combining population genomic, phenotypic, and environmental data from 700 individuals, we show how the genetics of ecotype determination mediates this evolutionary response. Variation in three terpene‐synthase loci explains variation in ecotype identity, with one single locus accounting for as much as 78% of this variation. Phenotypic selection combined with low segregating genotypic redundancy of ecotypes leaves a clear footprint at the genomic level: alleles associated with ecotype identity track environmental variation despite extensive gene flow. Different chemotypes within each ecotype differentially track environmental variation. Their identity is determined by multiple loci and displays a wider range of genotypic redundancy that dilutes phenotypic selection on their characteristic alleles. Our study thus provides a novel illustration of how genetic redundancy of a phenotype modulates the ability of selection to maintain adaptive differentiation. Identifying the precise genetics of the chemical polymorphism in thyme is the next crucial step for our understanding of the origin and maintenance of a polymorphism that is present in many aromatic plants.

Assessing the hydrogeological resilience of a groundwater-dependent Mediterranean peatland: Impact of global change and role of water management strategies

Mediterranean peatlands remain largely under-documented, except for detailed biological data such as fauna and flora taxa lists, and yet are increasingly threatened by water withdrawal and agriculture practices. This lack of information, particularly on their hydrogeological functioning, makes it impossible to evaluate their response to changing climate conditions. A pilot study on a representative Mediterranean peatland on the island of Corsica (France) was conducted to evaluate recharge modalities in the peatland using a coupled water-level monitoring, geochemical and isotope multi-tracing approach (electric conductivity, major ions, δ¹⁸O, δ²H, ³H, ⁸⁷Sr/⁸⁶Sr). The goal was to understand how water budgets in peatland ecosystems were maintained throughout the year, especially during the summer. Despite the remarkable stability of the peatland water level, the recharge contributions of varied water bodies through an alluvial aquifer vary significantly from one season to another. An end-member mixing analysis (EMMA) indicates that the peatland is mainly recharged by an alluvial aquifer. During fall-winter, the alluvial aquifer on which the peatland depends is recharged by the rainfall, a river, and shallow groundwater (colluvium). During spring-summer, water supply is provided mostly by a river, shallow, and deep groundwater (fractured granite). However, this specific hydrogeological functioning is not taken into account by environmental management policies making peatlands vulnerable to anthropogenic and climatic pressures. Thus, their actual status regarding water and aquatic environment management policies is discussed to provide recommendations for better consideration and preservation.

Temporal offset between precipitation and water uptake of Mediterranean pine trees varies with elevation and season

For climate models that use paleo-environment data to predict future climate change, tree-ring isotope variations are one important archive for the reconstruction of paleo-hydrological conditions. Due to the rather complicated pathway of water, starting from precipitation until its uptake by trees and the final incorporation of its components into tree-ring cellulose, a closer inspection of seasonal variations of tree water uptake is important. In this study, branch and needle samples of two pine species (Pinus pinaster and Pinus nigra subsp. laricio) and several water compartments (precipitation, creek, soil) were sampled over a two-year period and analyzed for the temporal variations of their oxygen and hydrogen stable isotope ratios (δ¹⁸O and δ²H) at five sites over an elevation gradient from sea level to around 1600 m a.s.l. on the Mediterranean island of Corsica (France). A new model was established to disentangle temporal relationships of source water uptake of trees. It uses a calculation method that incorporates the two processes mostly expected to affect source water composition: mixing of waters and evaporation. The model results showed that the temporal offset from precipitation to water uptake is not constant and varies with elevation and season. Overall, seasonal source water origin was shown to be dominated by precipitation from autumn and spring. While autumn precipitation was a more important water source for trees growing at mid- (~800-1000 m a.s.l) and high-elevation (~1600 m a.s.l.) sites, trees at coastal sites mostly took up water from late winter and spring. These findings show that predicted decreases in precipitation amounts during the wet season in the Mediterranean can have strong impacts on water availability for pine trees, especially at higher elevations.

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.

The wild grape genome sequence provides insights into the transition from dioecy to hermaphroditism during grape domestication

Background

A key step in domestication of the grapevine was the transition from separate sexes (dioecy) in wild Vitis vinifera ssp. sylvestris (V. sylvestris) to hermaphroditism in cultivated Vitis vinifera ssp. sativa (V. vinifera). It is known that V. sylvestris has an XY system and V. vinifera a modified Y haplotype (Yh) and that the sex locus is small, but it has not previously been precisely characterized.

Results

We generate a high-quality de novo reference genome for V. sylvestris, onto which we map whole-genome re-sequencing data of a cross to locate the sex locus. Assembly of the full X, Y, and Yh haplotypes of V. sylvestris and V. vinifera sex locus and examining their gene content and expression profiles during flower development in wild and cultivated accessions show that truncation and deletion of tapetum and pollen development genes on the X haplotype likely causes male sterility, while the upregulation of a Y allele of a cytokinin regulator (APRT3) may cause female sterility. The downregulation of this cytokinin regulator in the Yh haplotype may be sufficient to trigger reversal to hermaphroditism. Molecular dating of X and Y haplotypes is consistent with the sex locus being as old as the Vitis genus, but the mechanism by which recombination was suppressed remains undetermined.

Conclusions

We describe the genomic and evolutionary characterization of the sex locus of cultivated and wild grapevine, providing a coherent model of sex determination in the latter and for transition from dioecy to hermaphroditism during domestication.

Genome-wide association mapping of QTLs implied in potato virus Y population sizes in pepper: evidence for widespread resistance QTL pyramiding

In this study, we looked for genetic factors in the pepper (Capsicum annuum) germplasm that control the number of potato virus Y (PVY) particles entering the plant (i.e. effective population size at inoculation) and the PVY accumulation at the systemic level (i.e. census population size). Using genotyping-by-sequencing (GBS) in a core collection of 256 pepper accessions, we obtained 10 307 single nucleotide polymorphisms (SNPs) covering the whole genome. Genome-wide association studies (GWAS) detected seven SNPs significantly associated with the virus population size at inoculation and/or systemic level on chromosomes 4, 6, 9 and 12. Two SNPs on chromosome 4 associated with both PVY population sizes map closely to the major resistance gene pvr2 encoding the eukaryotic initiation factor 4E. No obvious candidates for resistance were identified in the confidence intervals for the other chromosomes. SNPs detected on chromosomes 6 and 12 colocalized with resistance quantitative trait loci (QTLs) previously identified with a biparental population. These results show the efficiency of GBS and GWAS in C. annuum, indicate highly consistent results between GWAS and classical QTL mapping, and suggest that resistance QTLs identified with a biparental population are representative of a much larger collection of pepper accessions. Moreover, the resistance alleles at these different loci were more frequently combined than expected by chance in the core collection, indicating widespread pyramiding of resistance QTLs and widespread combination of resistance QTLs and major effect genes. Such pyramiding may increase resistance efficiency and/or durability.

Evolutionary transcriptomics reveals the origins of olives and the genomic changes associated with their domestication

The olive (Olea europaea L. subsp. europaea) is one of the oldest and most socio-economically important cultivated perennial crop in the Mediterranean region. Yet, its origins are still under debate and the genetic bases of the phenotypic changes associated with its domestication are unknown. We generated RNA-sequencing data for 68 wild and cultivated olive trees to study the genetic diversity and structure both at the transcription and sequence levels. To localize putative genes or expression pathways targeted by artificial selection during domestication, we employed a two-step approach in which we identified differentially expressed genes and screened the transcriptome for signatures of selection. Our analyses support a major domestication event in the eastern part of the Mediterranean basin followed by dispersion towards the West and subsequent admixture with western wild olives. While we found large changes in gene expression when comparing cultivated and wild olives, we found no major signature of selection on coding variants and weak signals primarily affected transcription factors. Our results indicated that the domestication of olives resulted in only moderate genomic consequences and that the domestication syndrome is mainly related to changes in gene expression, consistent with its evolutionary history and life history traits.

A novel high-density grapevine (Vitis vinifera L.) integrated linkage map using GBS in a half-diallel population

A half-diallel population involving five elite grapevine cultivars was generated and genotyped by GBS, and highly-informative segregation data was used to construct a high-density genetic map for Vitis vinifera L. Grapevine is one of the most relevant fruit crops in the world. Deeper genetic knowledge could assist modern grapevine breeding programs to develop new wine grape varieties able to face climate change effects. To assist in the rapid identification of markers for crop yield components, grape quality traits and adaptation potential, we generated a large Vitis vinifera L. population (N = 624) by crossing five red wine cultivars in a half-diallel scheme, which was subsequently sequenced by an efficient GBS procedure. A high number of fully informative genetic variants was detected using a novel mapping approach capable of reconstructing local haplotypes from adjacent biallelic SNPs, which were subsequently used to construct the densest consensus genetic map available for the cultivated grapevine to date. This 1378.3-cM map integrates 10 bi-parental consensus maps and orders 4437 markers in 3353 unique positions on 19 chromosomes. Markers are well distributed all along the grapevine reference genome, covering up to 98.8% of its genomic sequence. Additionally, a good agreement was observed between genetic and physical orders, adding confidence in the quality of this map. Collectively, our results pave the way for future genetic studies (such as fine QTL mapping) aimed to understand the complex relationship between genotypic and phenotypic variation in the cultivated grapevine. In addition, the method used (which efficiently delivers a high number of fully informative markers) could be of interest to other outbred organisms, notably perennial fruit crops.

Pervasive hybridizations in the history of wheat relatives

Cultivated wheats are derived from an intricate history of three genomes, A, B, and D, present in both diploid and polyploid species. It was recently proposed that the D genome originated from an ancient hybridization between the A and B lineages. However, this result has been questioned, and a robust phylogeny of wheat relatives is still lacking. Using transcriptome data from all diploid species and a new methodological approach, our comprehensive phylogenomic analysis revealed that more than half of the species descend from an ancient hybridization event but with a more complex scenario involving a different parent than previously thought-Aegilops mutica, an overlooked wild species-instead of the B genome. We also detected other extensive gene flow events that could explain long-standing controversies in the classification of wheat relatives.

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.

Combinations of geoenvironmental data underline coastal aquifer anthropogenic nitrate legacy through groundwater vulnerability mapping methods

Groundwater quality is strongly dependent on land use. Past and current anthropogenic activities can lead to the diffusion of contaminants in aquifers. This diffusion can threaten the resource exploitation for decades, thereby endangering the ecological health of groundwater dependent ecosystems. Thus, groundwater stakeholders need methods for long-term management which integrate groundwater vulnerability. This study was conducted on the shallow alluvial aquifer of the groundwater-dependent Biguglia lagoon on Corsica Island, France. The aquifer is exposed to anthropogenic contamination for many decades with nitrate contamination legacy linked to agricultural activities, uncontrolled urbanization and sewage leakages. In most cases, vulnerability mapping is done in the objective of comparing groundwater situation regarding an on-going contamination process. But the question is still pending for aquifers where contamination is inherited from past practices or contaminations and where anthropogenic influences have changed through time. To propose an effective and innovative method for territorial management in Mediterranean alluvial aquifers, four index-based groundwater vulnerability mapping methods were tested and compared: two intrinsic vulnerability mapping methods (DRASTIC and SINTACS) and two specific vulnerability mapping methods (Modified-DRASTIC and SI), the latter integrating land use in the accuracy of groundwater vulnerability. Novelty is coming from the comparison between vulnerability maps and their application and validation in a hydrosystem affected by nitrate legacy-type contamination. The specific vulnerability mapping methods are more likely to represent the current pressures to which groundwater are subject. Thus, specific vulnerability methods such as the SI one revealed here very relevant to assess groundwater quality and to react retrospectively. The comparison between groundwater nitrate legacy and intrinsic groundwater vulnerability methods appeared also useful to define priority protection areas in long-term territorial management planning (EU Water Framework Direction). In this sense, the SINTACS method seems to be the more appropriate in the Mediterranean and alluvial context of this study.

Riverine carbon dioxide evasion along a high-relief watercourse derived from seasonal dynamics of the water-atmosphere gas exchange

The high-relief catchment of the Tavignanu River (Corsica Island, France) with an elevation range from sea level to 2622 m above sea level was investigated for its riverine carbon budget and stable carbon isotopes. Major riverine dissolved inorganic carbon (DIC or TCO₂) sources depended on seasons and sub-catchment lithology. In winter δ¹³C_DIC values of -2 to -7‰ (VPDB) indicated influences of atmospheric CO₂. δ¹³C_DIC values decreased gradually to values between -9 and -12‰ in July, which indicates elevated soil CO₂ contribution. An observed downstream increase in the total amount of carbon species correlated with inputs from carbonate bearing tributaries and evaporation effects in summer. Main channel partial pressure of CO₂ (pCO₂) was seasonally highly variable in the upper silicate catchment and the lower coastal plain, where summer values exceed up to six times atmospheric levels. During winter, the central section of the Tavignanu River was found to be undersaturated with respect to atmospheric CO₂ levels. The median values for main channel pCO₂ were below atmospheric levels in winter and spring and above in summer and autumn. The annual carbon flux across the air-water boundary (F_CO2) along the Tavignanu River was calculated with (0.77 ± 0.24) Gg C yr^-1, which is about seven times higher than the riverine TCO₂ transport to the sea of about 0.11 Gg C yr^-1. While large sections of the river experienced year-round atmospheric CO₂ uptake or equilibrium, the river as a whole was a small but continuous net source of carbon to the atmosphere. This underlines the important, but so far not well-constrained, contributions of smaller streams and rivers to the terrestrial carbon flux and the need of incorporating them into future global carbon cycle models.

Delayed nitrate dispersion within a coastal aquifer provides constraints on land-use evolution and nitrate contamination in the past

Identifying sources of anthropogenic pollution, and assessing the fate and residence time of pollutants in aquifers is important for the management of groundwater resources, and the ecological health of groundwater dependent ecosystems. This study investigates anthropogenic contamination in the shallow alluvial aquifer of the Marana-Casinca, hydraulically connected to the Biguglia lagoon (Corsica, France). A multi-tracer approach, combining geochemical and environmental isotopic data (δ¹⁸O-H₂O, δ²H-H₂O, ³H, δ¹⁵N-NO₃^-, δ¹⁸O-NO₃^-, δ¹¹B), and groundwater residence-time tracers (³H and CFCs) was carried out in 2016, and integrated with a study of land use evolution in the catchment during the last century. Groundwater NO₃^- concentrations, ranged between 2 mg/L and up to 30 mg/L, displaying the degradation of groundwater quality induced by anthropogenic activities (agricultural activities). Comparatively high δ¹⁵N-NO₃^- values (up to 19.7‰) in combination with δ¹¹B values that were significantly lower (between 23‰ and 26‰) than the seawater background are indicative of sewage contamination. The ongoing deterioration of groundwater quality can be attributed to the uncontrolled urbanization development all over the alluvial plain, with numerous sewage leakages from the sanitation network and private sewage systems. Integration of contaminant and water-residence time data revealed a progressive accumulation of pollutants with time in the groundwater, particularly in areas with major anthropogenic pressure and slow dynamic groundwater flow. Our approach provides time-dependent insight into nitrogen pollution in the studied aquifer over the past decades, revealing a systematic change in the dominant NO₃^- source, from agricultural to sewage contamination. Yet, today's low groundwater quality is to large parts due to legacy pollution from land-use practices several decades ago, underlining the poor self-remediating capacity of this hydrosystem. Our results can be taken as warning that groundwater pollution that happened in the recent past, or today, may have dire impacts on the quality of groundwater-dependent ecosystems in the future.

Allele-specific expression and genetic determinants of transcriptomic variations in response to mild water deficit in tomato

Characterizing the natural diversity of gene expression across environments is an important step in understanding how genotype-by-environment interactions shape phenotypes. Here, we analyzed the impact of water deficit onto gene expression levels in tomato at the genome-wide scale. We sequenced the transcriptome of growing leaves and fruit pericarps at cell expansion stage in a cherry and a large fruited accession and their F₁ hybrid grown under two watering regimes. Gene expression levels were steadily affected by the genotype and the watering regime. Whereas phenotypes showed mostly additive inheritance, ~80% of the genes displayed non-additive inheritance. By comparing allele-specific expression (ASE) in the F₁ hybrid to the allelic expression in both parental lines, respectively, 3005 genes in leaf and 2857 genes in fruit deviated from 1:1 ratio independently of the watering regime. Among these genes, ~55% were controlled by cis factors, ~25% by trans factors and ~20% by a combination of both types of factors. A total of 328 genes in leaf and 113 in fruit exhibited significant ASE-by-watering regime interaction, among which ~80% presented trans-by-watering regime interaction, suggesting a response to water deficit mediated through a majority of trans-acting loci in tomato. We cross-validated the expression levels of 274 transcripts in fruit and leaves of 124 recombinant inbred lines (RILs) and identified 163 expression quantitative trait loci (eQTLs) mostly confirming the divergences identified by ASE. Combining phenotypic and expression data, we observed a complex network of variation between genes encoding enzymes involved in the sugar metabolism.

An Analytical Method for Assessing Recharge Using Groundwater Travel Time in Dupuit-Forchheimer Aquifers

An analytical solution to calculate the recharge of unconfined aquifers with Dupuit-Forchheimer type flow conditions is proposed. This solution is derived from an existing closed-form analytical solution initially developed to determine groundwater travel time when the recharge of the aquifer is known. This existing solution has been modified to determine recharge when groundwater travel time is known. An illustration is given with a field case example for the Bonifacio aquifer of the island of Corsica (France), in the Mediterranean. In this aquifer, previously established differences in groundwater residence time between two water samples were determined from anthropogenic atmospheric gas (chlorofluorocarbons and sulfur hexafluoride) measurements. The time difference is entered into the new analytical solution to determine recharge. The calculations yield a value of average recharge that agrees with the results obtained by several other methods that were presented in previous studies to assess the recharge of the Bonifacio aquifer. Also presented in this study is a sensitivity analysis of the new analytical solution, to quantify the influence of different parameters that control recharge: hydraulic conductivity, effective porosity and the groundwater travel time. This study illustrates how geochemical data can be combined with physical models to measure recharge. Such an approach could be adopted in other homogeneous aquifers worldwide that satisfy Dupuit-Forchheimer type flow conditions.

Using insects to detect, monitor and predict the distribution of Xylella fastidiosa: a case study in Corsica

We sampled ca 2500 specimens of Philaenus spumarius (Hemiptera: Aphrophoridae) throughout Corsica without a priori knowledge on the presence of symptoms on plants. We screened 448 specimens for the presence of Xylella fastidiosa (Xf) using qPCR and a custom nested PCR. qPCR appeared versatile and under-estimated the prevalence of Xf. Nested PCR showed that Xf was present in all populations. Molecular results were validated by prediction on the distribution of Xf made from tests conducted on plants, which shows the pertinence of using vectors in risk assessment studies. Xf was detected in tenerals and adults. Thus, P. spumarius could acquire Xf from its host plant, mostly Cistus monspeliensis in Corsica, which may act as reservoir for the next season. This contrasts with other observations and suggests that management strategies may have to be adapted on a case-by-case basis. At least two genetic entities and several variants of Xf not yet identified on plants were present in the insects, which suggests ancient introductions of Xf and a probable underestimation of the current diversity of the strains present in Corsica. Interestingly 6% of the specimens carried two subspecies of Xf. Studies are required to better characterize the strains present in Corsica and to determine how the disease was introduced, spread and why no sign of a potential epidemic was detected earlier. This study shows that, when sensitive enough methods are implemented, spittlebugs (and more specifically P. spumarius for which species distribution modelling shows it could be a good sentinel for Europe) can be used to predict and better assess the exact distribution of Xf. Furthermore, Xf multiply only in their foregut and does not become circulative, which facilitates its detection.

Patterns of Genome-Wide Nucleotide Diversity in the Gynodioecious Plant Thymus vulgaris Are Compatible with Recent Sweeps of Cytoplasmic Genes

Gynodioecy is a sexual dimorphism where females coexist with hermaphrodite individuals. In most cases, this dimorphism involves the interaction of cytoplasmic male sterility (CMS) genes and nuclear restorer genes. Two scenarios can account for how these interactions maintain gynodioecy. Either CMS genes recurrently enter populations at low frequency via mutation or migration and go to fixation unimpeded (successive sweeps), or CMS genes maintain polymorphism over evolutionary time through interactions with a nuclear restorer allele (balanced polymorphism). To distinguish between these scenarios, we used transcriptome sequencing in gynodioecious Thymus vulgaris and surveyed genome-wide diversity in 18 naturally occurring individuals sampled from populations at a local geographic scale. We contrast the amount and patterns of nucleotide diversity in the nuclear and cytoplasmic genome, and find ample diversity at the nuclear level (π = 0.019 at synonymous sites) but reduced genetic diversity and an excess of rare polymorphisms in the cytoplasmic genome relative to the nuclear genome. Our finding is incompatible with the maintenance of gynodioecy via scenarios invoking long-term balancing selection, and instead suggests the recent fixation of CMS lineages in the populations studied.

Genotype-Environment Interaction Shapes the Microbial Assemblage in Grapevine's Phyllosphere and Carposphere: An NGS Approach

Plant surface or phyllosphere is the habitat of hyperdiverse microbial communities and it is always exposed to the fluctuating environmental factors, which is thought to be one of the potential drivers of microbial community structuring. Impact of grapevine genotypes in variable environmental factors (i.e., at different geographic locations) on the phyllosphere has never been studied and is the main objective of this report. Using high throughput short amplicon sequencing of 16S rRNA genes and internal transcribed spacer (ITS), we analyzed the impacts of genotypes of Vitis Vinifera (coming from three genetic pool), on the microbial (bacterial and fungal) assemblage in the phyllosphere. First, we performed the analysis of the phyllosphere microbiome while using fifteen genotypes that were chosen to maximize intra-specific diversity and grown in two Mediterranean vineyards. Then, the same analysis was performed on five commercially important varieties of Vitis vinifera that were sampled from three different French agro-climatic zones (or terroir: a combination of climate, soils, and human practices). Our study revealed that, at a particular geographic location, genotypes have an impact on microbial assemblage in the phyllosphere and carposphere of leaf and fruit (or berries), respectively, which is more prominent on the carposphere but the effect of terroir was much stronger than the genotype when the leaf phyllosphere of five grapevine varieties grown in different agro-climatic zones was compared. Impacts of the season and exterior plant organs (leaf and berries) on microbial taxa structuring in the phyllosphere was also assessed and presented in this report.

Microbial communities in natural rubber coagula during maturation: impacts on technological properties of dry natural rubber

AIM

To characterize microbial communities present in natural rubber (NR) coagula from Hevea brasiliensis latex during maturation and identify microbial taxa (bacteria and fungi) having an impact on dry NR properties.

METHODS AND RESULTS

Microbial community dynamics in NR coagula maturated under controlled conditions were compared and related with the evolution of dry NR properties. The pyrosequencing of 16S (119 837 effective reads) and 18S (131 879 effective reads) rRNA gene regions was performed on 21 samples covering different maturation times and two aeration conditions. Results showed a relatively high bacterial richness (Chao1 estimates of 200-1000) associated with significant bacterial dynamics. Lactic acid bacteria (LAB) were dominant in the first days of maturation. Then, in aerobic conditions, development of Actinobacteria represented by the family Microbacteriaceae was associated with alkalinization of the samples and a higher sensitivity of NR to thermo-oxidation as evaluated by its plasticity retention index (PRI). In anaerobiosis, the reduced development of bacteria, mostly LAB present, was associated with improved NR properties (higher initial plasticity P₀ and PRI).

CONCLUSIONS

The involvement of micro-organisms in the evolution of dry NR properties during the maturation of NR coagula was confirmed. The importance of the structure and dynamics of microbial communities is specifically highlighted.

SIGNIFICANCE AND IMPACT OF THE STUDY

Natural rubber is a key elastomer for the tyre industry and for a variety of other applications. The majority of raw NR is obtained by natural coagulation of H. brasiliensis latex under the activity of micro-organisms. An improved understanding of the microbial communities involved in the maturation of NR coagula may lead to an improvement in the production process of raw NR to provide a better consistency in NR quality.

How to optimize the precision of allele and haplotype frequency estimates using pooled-sequencing data

Sequencing pools of individuals rather than individuals separately reduces the costs of estimating allele frequencies at many loci in many populations. Theoretical and empirical studies show that sequencing pools comprising a limited number of individuals (typically fewer than 50) provides reliable allele frequency estimates, provided that the DNA pooling and DNA sequencing steps are carefully controlled. Unequal contributions of different individuals to the DNA pool and the mean and variance in sequencing depth both can affect the standard error of allele frequency estimates. To our knowledge, no study separately investigated the effect of these two factors on allele frequency estimates; so that there is currently no method to a priori estimate the relative importance of unequal individual DNA contributions independently of sequencing depth. We develop a new analytical model for allele frequency estimation that explicitly distinguishes these two effects. Our model shows that the DNA pooling variance in a pooled sequencing experiment depends solely on two factors: the number of individuals within the pool and the coefficient of variation of individual DNA contributions to the pool. We present a new method to experimentally estimate this coefficient of variation when planning a pooled sequencing design where samples are either pooled before or after DNA extraction. Using this analytical and experimental framework, we provide guidelines to optimize the design of pooled sequencing experiments. Finally, we sequence replicated pools of inbred lines of the plant Medicago truncatula and show that the predictions from our model generally hold true when estimating the frequency of known multilocus haplotypes using pooled sequencing.

Domestication rewired gene expression and nucleotide diversity patterns in tomato

Plant domestication has led to considerable phenotypic modifications from wild species to modern varieties. However, although changes in key traits have been well documented, less is known about the underlying molecular mechanisms, such as the reduction of molecular diversity or global gene co-expression patterns. In this study, we used a combination of gene expression and population genetics in wild and crop tomato to decipher the footprints of domestication. We found a set of 1729 differentially expressed genes (DEG) between the two genetic groups, belonging to 17 clusters of co-expressed DEG, suggesting that domestication affected not only individual genes but also regulatory networks. Five co-expression clusters were enriched in functional terms involving carbohydrate metabolism or epigenetic regulation of gene expression. We detected differences in nucleotide diversity between the crop and wild groups specific to DEG. Our study provides an extensive profiling of the rewiring of gene co-expression induced by the domestication syndrome in one of the main crop species.

Epistatic determinism of durum wheat resistance to the wheat spindle streak mosaic virus

KEY MESSAGE

The resistance of durum wheat to the Wheat spindle streak mosaic virus (WSSMV) is controlled by two main QTLs on chromosomes 7A and 7B, with a huge epistatic effect. Wheat spindle streak mosaic virus (WSSMV) is a major disease of durum wheat in Europe and North America. Breeding WSSMV-resistant cultivars is currently the only way to control the virus since no treatment is available. This paper reports studies of the inheritance of WSSMV resistance using two related durum wheat populations obtained by crossing two elite cultivars with a WSSMV-resistant emmer cultivar. In 2012 and 2015, 354 recombinant inbred lines (RIL) were phenotyped using visual notations, ELISA and qPCR and genotyped using locus targeted capture and sequencing. This allowed us to build a consensus genetic map of 8568 markers and identify three chromosomal regions involved in WSSMV resistance. Two major regions (located on chromosomes 7A and 7B) jointly explain, on the basis of epistatic interactions, up to 43% of the phenotypic variation. Flanking sequences of our genetic markers are provided to facilitate future marker-assisted selection of WSSMV-resistant cultivars.

Strontium isotopes as tracers of water-rocks interactions, mixing processes and residence time indicator of groundwater within the granite-carbonate coastal aquifer of Bonifacio (Corsica, France)

This study aims at identifying the water-rock interactions and mixing rates within a complex granite-carbonate coastal aquifer under high touristic pressure. Investigations have been carried out within the coastal aquifer of Bonifacio (southern Corsica, France) mainly composed of continental granitic weathering products and marine calcarenite sediments filling a granitic depression. A multi-tracer approach combining physico-chemical parameters, major ions, selected trace elements, stable isotopes of the water molecule and ⁸⁷Sr/⁸⁶Sr ratios measurements is undertaken for 20 groundwater samples during the low water period in November 2014. 5 rock samples of the sedimentary deposits and surrounding granites are also analysed. First, the water-rock interactions processes governing the groundwater mineralization are described in order to fix the hydrogeochemical background. Secondly, the flow conditions are refined through the quantification of inter aquifer levels mixing, and thirdly, the kinetics of water-rock interaction based on groundwater residence time from a previous study using CFCs and SF₆ are quantified for the two main flow lines. A regional contrast in the groundwater recharge altitude allowed the oxygene-18 to be useful combined with the ⁸⁷Sr/⁸⁶Sr ratios to differentiate the groundwater origins and to compute the mixing rates, revealing the real extension of the watershed and the availability of the resource. The results also highlight a very good correlation between the groundwater residence time and the spatial evolution of ⁸⁷Sr/⁸⁶Sr ratios, allowing water-rock interaction kinetics to be defined empirically for the two main flow lines through the calcarenites. These results demonstrate the efficiency of strontium isotopes as tracers of water-rock interaction kinetics and by extension their relevance as a proxy of groundwater residence time, fundamental parameter documenting the long term sustainability of the hydrosystem.

A large set of 26 new reference transcriptomes dedicated to comparative population genomics in crops and wild relatives

We produced a unique large data set of reference transcriptomes to obtain new knowledge about the evolution of plant genomes and crop domestication. For this purpose, we validated a RNA-Seq data assembly protocol to perform comparative population genomics. For the validation, we assessed and compared the quality of de novo Illumina short-read assemblies using data from two crops for which an annotated reference genome was available, namely grapevine and sorghum. We used the same protocol for the release of 26 new transcriptomes of crop plants and wild relatives, including still understudied crops such as yam, pearl millet and fonio. The species list has a wide taxonomic representation with the inclusion of 15 monocots and 11 eudicots. All contigs were annotated using BLAST, prot4EST and Blast2GO. A strong originality of the data set is that each crop is associated with close relative species, which will permit whole-genome comparative evolutionary studies between crops and their wild-related species. This large resource will thus serve research communities working on both crops and model organisms. All the data are available at http://arcad-bioinformatics.southgreen.fr/.

A replicated climate change field experiment reveals rapid evolutionary response in an ecologically important soil invertebrate

Whether species can respond evolutionarily to current climate change is crucial for the persistence of many species. Yet, very few studies have examined genetic responses to climate change in manipulated experiments carried out in natural field conditions. We examined the evolutionary response to climate change in a common annelid worm using a controlled replicated experiment where climatic conditions were manipulated in a natural setting. Analyzing the transcribed genome of 15 local populations, we found that about 12% of the genetic polymorphisms exhibit differences in allele frequencies associated to changes in soil temperature and soil moisture. This shows an evolutionary response to realistic climate change happening over short-time scale, and calls for incorporating evolution into models predicting future response of species to climate change. It also shows that designed climate change experiments coupled with genome sequencing offer great potential to test for the occurrence (or lack) of an evolutionary response.

Endemic insular and coastal Tunisian date palm genetic diversity

The breeding of crop species relies on the valorisation of ancestral or wild varieties to enrich the cultivated germplasm. The Tunisian date palm genetic patrimony is being threatened by diversity loss and global climate change. We have conducted a genetic study to evaluate the potential of spontaneous coastal resources to improve the currently exploited Tunisian date palm genetic pool. Eighteen microsatellite loci of Phoenix dactylifera L. were used to compare the genetic diversity of coastal accessions from Kerkennah, Djerba, Gabès and continental date palm accessions from Tozeur. A collection of 105 date palms from the four regions was analysed. This study has provided us with an extensive understanding of the local genetic diversity and its distribution. The coastal date palm genotypes exhibit a high and specific genetic diversity. These genotypes are certainly an untapped reservoir of agronomically important genes to improve cultivated germplasm in continental date palm.

Genetic structure of the date palm (Phoenix dactylifera) in the Old World reveals a strong differentiation between eastern and western populations

Since the publication of this paper, it has become apparent that an error was made in the legend to Fig. 3 and the colours referring to occidental and oriental are the wrong way round. The authors apologise for this error, and a correct version of the legend to Fig. 3 is given below.

Genetic structure of the date palm (Phoenix dactylifera) in the Old World reveals a strong differentiation between eastern and western populations

BACKGROUND AND AIMS

Date palms (Phoenix dactylifera, Arecaceae) are of great economic and ecological value to the oasis agriculture of arid and semi-arid areas. However, despite the availability of a large date palm germplasm spreading from the Atlantic shores to Southern Asia, improvement of the species is being hampered by a lack of information on global genetic diversity and population structure. In order to contribute to the varietal improvement of date palms and to provide new insights on the influence of geographic origins and human activity on the genetic structure of the date palm, this study analysed the diversity of the species.

METHODS

Genetic diversity levels and population genetic structure were investigated through the genotyping of a collection of 295 date palm accessions ranging from Mauritania to Pakistan using a set of 18 simple sequence repeat (SSR) markers and a plastid minisatellite.

KEY RESULTS

Using a Bayesian clustering approach, the date palm genotypes can be structured into two different gene pools: the first, termed the Eastern pool, consists of accessions from Asia and Djibouti, whilst the second, termed the Western pool, consists of accessions from Africa. These results confirm the existence of two ancient gene pools that have contributed to the current date palm diversity. The presence of admixed genotypes is also noted, which points at gene flows between eastern and western origins, mostly from east to west, following a human-mediated diffusion of the species.

CONCLUSIONS

This study assesses the distribution and level of genetic diversity of accessible date palm resources, provides new insights on the geographic origins and genetic history of the cultivated component of this species, and confirms the existence of at least two domestication origins. Furthermore, the strong genetic structure clearly established here is a prerequisite for any breeding programme exploiting the effective polymorphism related to each gene pool.

Assessment of DNA encapsulation, a new room-temperature DNA storage method

A new procedure for room-temperature storage of DNA was evaluated whereby DNA samples from human tissue, bacteria, and plants were stored under an anoxic and anhydrous atmosphere in small glass vials fitted in stainless-steel, laser-sealed capsules (DNAshells(®)). Samples were stored in DNAshells(®) at room temperature for various periods of time to assess any degradation and compare it to frozen control samples and those stored in GenTegra™ tubes. The study included analysis of the effect of accelerated aging by using a high temperature (76°C) at 50% relative humidity. No detectable DNA degradation was seen in samples stored in DNAshells(®) at room temperature for 18 months. Polymerase chain reaction experiments, pulsed field gel electrophoresis, and amplified fragment length polymorphism analyses also demonstrated that the protective properties of DNAshells(®) are not affected by storage under extreme conditions (76°C, 50% humidity) for 30 hours, guaranteeing 100 years without DNA sample degradation. However, after 30 hours of storage at 76°C, it was necessary to include adjustments to the process in order to avoid DNA loss. Successful protection of DNA was obtained for 1 week and even 1 month of storage at high temperature by adding trehalose, which provides a protective matrix. This study demonstrates the many advantages of using DNAshells(®) for room-temperature storage, particularly in terms of long-term stability, safety, transport, and applications for molecular biology research.

A small XY chromosomal region explains sex determination in wild dioecious V. vinifera and the reversal to hermaphroditism in domesticated grapevines

BACKGROUND

In Vitis vinifera L., domestication induced a dramatic change in flower morphology: the wild sylvestris subspecies is dioecious while hermaphroditism is largely predominant in the domesticated subsp. V. v. vinifera. The characterisation of polymorphisms in genes underlying the sex-determining chromosomal region may help clarify the history of domestication in grapevine and the evolution of sex chromosomes in plants. In the genus Vitis, sex determination is putatively controlled by one major locus with three alleles, male M, hermaphrodite H and female F, with an allelic dominance M > H > F. Previous genetic studies located the sex locus on chromosome 2. We used DNA polymorphisms of geographically diverse V. vinifera genotypes to confirm the position of this locus, to characterise the genetic diversity and traces of selection in candidate genes, and to explore the origin of hermaphroditism.

RESULTS

In V. v. sylvestris, a sex-determining region of 154.8 kb, also present in other Vitis species, spans less than 1% of chromosome 2. It displays haplotype diversity, linkage disequilibrium and differentiation that typically correspond to a small XY sex-determining region with XY males and XX females. In male alleles, traces of purifying selection were found for a trehalose phosphatase, an exostosin and a WRKY transcription factor, with strikingly low polymorphism levels between distant geographic regions. Both diversity and network analysis revealed that H alleles are more closely related to M than to F alleles.

CONCLUSIONS

Hermaphrodite alleles appear to derive from male alleles of wild grapevines, with successive recombination events allowing import of diversity from the X into the Y chromosomal region and slowing down the expansion of the region into a full heteromorphic chromosome. Our data are consistent with multiple domestication events and show traces of introgression from other Asian Vitis species into the cultivated grapevine gene pool.

Transcriptome population genomics reveals severe bottleneck and domestication cost in the African rice (Oryza glaberrima)

The African cultivated rice (Oryza glaberrima) was domesticated in West Africa 3000 years ago. Although less cultivated than the Asian rice (O. sativa), O. glaberrima landraces often display interesting adaptation to rustic environment (e.g. drought). Here, using RNA-seq technology, we were able to compare more than 12,000 transcripts between 9 O. glaberrima, 10 wild O. barthii and one O. meridionalis individuals. With a synonymous nucleotide diversity πs = 0.0006 per site, O. glaberrima appears as the least genetically diverse crop grass ever documented. Using approximate Bayesian computation, we estimated that O. glaberrima experienced a severe bottleneck during domestication. This demographic scenario almost fully accounts for the pattern of genetic diversity across O. glaberrima genome as we detected very few outliers regions where positive selection may have further impacted genetic diversity. Moreover, the large excess of derived nonsynonymous substitution that we detected suggests that the O. glaberrima population suffered from the 'cost of domestication'. In addition, we used this genome-scale data set to demonstrate that (i) O. barthii genetic diversity is positively correlated with recombination rate and negatively with gene density, (ii) expression level is negatively correlated with evolutionary constraint, and (iii) one region on chromosome 5 (position 4-6 Mb) exhibits a clear signature of introgression with a yet unidentified Oryza species. This work represents the first genome-wide survey of the African rice genetic diversity and paves the way for further comparison between the African and the Asian rice, notably regarding the genetics underlying domestication traits.

In silico mining of microsatellites in coding sequences of the date palm (Arecaceae) genome, characterization, and transferability

PREMISE OF THE STUDY

To complement existing sets of primarily dinucleotide microsatellite loci from noncoding sequences of date palm, we developed primers for tri- and hexanucleotide microsatellite loci identified within genes. Due to their conserved genomic locations, the primers should be useful in other palm taxa, and their utility was tested in seven other Phoenix species and in Chamaerops, Livistona, and Hyphaene. •

METHODS AND RESULTS

Tandem repeat motifs of 3-6 bp were searched using a simple sequence repeat (SSR)-pipeline package in coding portions of the date palm draft genome sequence. Fifteen loci produced highly consistent amplification, intraspecific polymorphisms, and stepwise mutation patterns. •

CONCLUSIONS

These microsatellite loci showed sufficient levels of variability and transferability to make them useful for population genetic, selection signature, and interspecific gene flow studies in Phoenix and other Coryphoideae genera.

Developing nuclear DNA phylogenetic markers in the angiosperm genus Leucadendron (Proteaceae): a next-generation sequencing transcriptomic approach

Despite the recent advances in generating molecular data, reconstructing species-level phylogenies for non-models groups remains a challenge. The use of a number of independent genes is required to resolve phylogenetic relationships, especially for groups displaying low polymorphism. In such cases, low-copy nuclear exons and non-coding regions, such as 3' untranslated regions (3'-UTRs) or introns, constitute a potentially interesting source of nuclear DNA variation. Here, we present a methodology meant to identify new nuclear orthologous markers using both public-nucleotide databases and transcriptomic data generated for the group of interest by using next generation sequencing technology. To identify PCR primers for a non-model group, the genus Leucadendron (Proteaceae), we adopted a framework aimed at minimizing the probability of paralogy and maximizing polymorphism. We anchored when possible the right-hand primer into the 3'-UTR and the left-hand primer into the coding region. Seven new nuclear markers emerged from this search strategy, three of those included 3'-UTRs. We further compared the phylogenetic potential between our new markers and the ribosomal internal transcribed spacer region (ITS). The sequenced 3'-UTRs yielded higher polymorphism rates than the ITS region did. We did not find strong incongruences with the phylogenetic signal contained in the ITS region and the seven new designed markers but they strongly improved the phylogeny of the genus Leucadendron. Overall, this methodology is efficient in isolating orthologous loci and is valid for any non-model group given the availability of transcriptomic data.

High genetic and epigenetic stability in Coffea arabica plants derived from embryogenic suspensions and secondary embryogenesis as revealed by AFLP, MSAP and the phenotypic variation rate

Embryogenic suspensions that involve extensive cell division are risky in respect to genome and epigenome instability. Elevated frequencies of somaclonal variation in embryogenic suspension-derived plants were reported in many species, including coffee. This problem could be overcome by using culture conditions that allow moderate cell proliferation. In view of true-to-type large-scale propagation of C. arabica hybrids, suspension protocols based on low 2,4-D concentrations and short proliferation periods were developed. As mechanisms leading to somaclonal variation are often complex, the phenotypic, genetic and epigenetic changes were jointly assessed so as to accurately evaluate the conformity of suspension-derived plants. The effects of embryogenic suspensions and secondary embryogenesis, used as proliferation systems, on the genetic conformity of somatic embryogenesis-derived plants (emblings) were assessed in two hybrids. When applied over a 6 month period, both systems ensured very low somaclonal variation rates, as observed through massive phenotypic observations in field plots (0.74% from 200,000 plant). Molecular AFLP and MSAP analyses performed on 145 three year-old emblings showed that polymorphism between mother plants and emblings was extremely low, i.e. ranges of 0-0.003% and 0.07-0.18% respectively, with no significant difference between the proliferation systems for the two hybrids. No embling was found to cumulate more than three methylation polymorphisms. No relation was established between the variant phenotype (27 variants studied) and a particular MSAP pattern. Chromosome counting showed that 7 of the 11 variant emblings analyzed were characterized by the loss of 1-3 chromosomes. This work showed that both embryogenic suspensions and secondary embryogenesis are reliable for true-to-type propagation of elite material. Molecular analyses revealed that genetic and epigenetic alterations are particularly limited during coffee somatic embryogenesis. The main change in most of the rare phenotypic variants was aneuploidy, indicating that mitotic aberrations play a major role in somaclonal variation in coffee.

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.

Male-specific DNA markers provide genetic evidence of an XY chromosome system, a recombination arrest and allow the tracing of paternal lineages in date palm

Whether sex chromosomes are differentiated is an important aspect of our knowledge of dioecious plants, such as date palm (Phoenix dactylifera). In this crop plant, the female individuals produce dates, and are thus the more valuable sex. However, there is no way to identify the sex of date palm plants before reproductive age, and the sex-determining mechanism is still unclear. To identify sex-linked microsatellite markers, we surveyed a set of 52 male and 55 female genotypes representing the geographical diversity of the species. We found three genetically linked loci that are heterozygous only in males. Male-specific alleles allowed us to identify the gender in 100% of individuals. These results confirm the existence of an XY chromosomal system with a nonrecombining XY-like region in the date palm genome. The distribution of Y haplotypes in western and eastern haplogroups allowed us to trace two male ancestral paternal lineages that account for all known Y diversity in date palm. The very low diversity associated with Y haplotypes is consistent with clonal paternal transmission of a nonrecombining male-determining region. Our results establish the date palm as a biological model with one of the most ancient sex chromosomes in flowering plants.

A composite map of expressed sequences in maize

A maize genetic map based mainly on expressed sequences has been constructed. The map incorporates data from four segregating populations. Three recombinant inbred line populations were derived from the nonreciprocal crosses between three inbred lines. A map derived from an independent F2 progeny from one of the crosses was also used. With a total of 521 genotyped individuals, accuracy in gene order is expected. Five sources of markers were used: (i) 109 loci corresponding to 69 genes of known function, (ii) 39 loci controlling protein position shifts revealed by two-dimensional electrophoresis, (iii) 8 isozyme loci, (iv) 17 loci corresponding to 14 sequenced cDNAs for which no homology was found in gene banks, and (v) 102 loci corresponding to 81 anonymous probes. As many loci were common to all maps, we tested heterogeneity between recombination fractions. The comparison of recombination fractions revealed: (i) a good correspondence between the maps derived from the same cross, (ii) few significant differences in interval distances, and (iii) global differences, which can reach 20% of the total map length. A composite map of 275 loci covering 1765 cM has been constructed. Key words : Zea mays L., RFLP, genetic map, molecular markers, proteins.

Loss of genetic diversity as a signature of apricot domestication and diffusion into the Mediterranean Basin

BACKGROUND

Domestication generally implies a loss of diversity in crop species relative to their wild ancestors because of genetic drift through bottleneck effects. Compared to native Mediterranean fruit species like olive and grape, the loss of genetic diversity is expected to be more substantial for fruit species introduced into Mediterranean areas such as apricot (Prunus armeniaca L.), which was probably primarily domesticated in China. By comparing genetic diversity among regional apricot gene pools in several Mediterranean areas, we investigated the loss of genetic diversity associated with apricot selection and diffusion into the Mediterranean Basin.

RESULTS

According to the geographic origin of apricots and using Bayesian clustering of genotypes, Mediterranean apricot (207 genotypes) was structured into three main gene pools: 'Irano-Caucasian', 'North Mediterranean Basin' and 'South Mediterranean Basin'. Among the 25 microsatellite markers used, only one displayed deviations from the frequencies expected under neutrality. Similar genetic diversity parameters were obtained within each of the three main clusters using both all SSR loci and only 24 SSR loci based on the assumption of neutrality. A significant loss of genetic diversity, as assessed by the allelic richness and private allelic richness, was revealed from the 'Irano-Caucasian' gene pool, considered as a secondary centre of diversification, to the northern and southwestern Mediterranean Basin. A substantial proportion of shared alleles was specifically detected when comparing gene pools from the 'North Mediterranean Basin' and 'South Mediterranean Basin' to the secondary centre of diversification.

CONCLUSIONS

A marked domestication bottleneck was detected with microsatellite markers in the Mediterranean apricot material, depicting a global image of two diffusion routes from the 'Irano-Caucasian' gene pool: North Mediterranean and Southwest Mediterranean. This study generated genetic insight that will be useful for management of Mediterranean apricot germplasm as well as genetic selection programs related to adaptive traits.

Deciphering the genetics of flowering time by an association study on candidate genes in bread wheat (Triticum aestivum L.)

Earliness is very important for the adaptation of wheat to environmental conditions and the achievement of high grain yield. A detailed knowledge of key genetic components of the life cycle would enable an easier control by the breeders. The objective of the study was to investigate the effect of candidate genes on flowering time. Using a collection of hexaploid wheat composed of 235 lines from diverse geographical origins, we conducted an association study for six candidate genes for flowering time and its components (vernalization sensitivity and earliness per se). The effect on the variation of earliness components of polymorphisms within the copies of each gene was tested in ANOVA models accounting for the underlying genetic structure. The collection was structured in five groups that minimized the residual covariance. Vernalization requirement and lateness tend to increase according to the mean latitude of each group. Heading date for an autumnal sowing was mainly determined by the earliness per se. Except for the Constans (CO) gene orthologous of the barley HvCO3, all gene polymorphisms had a significant impact on earliness components. The three traits used to quantify vernalization requirement were primarily associated with polymorphisms at Vrn-1 and then at Vrn-3 and Luminidependens (LD) genes. We found a good correspondence between spring/winter types and genotypes at the three homeologous copies of Vrn-1. Earliness per se was mainly explained by polymorphisms at Vrn-3 and to a lesser extent at Vrn-1, Hd-1 and Gigantea (GI) genes. Vernalization requirement and earliness as a function of geographical origin, as well as the possible role of the breeding practices in the geographical distribution of the alleles and the hypothetical adaptive value of the candidate genes, are discussed.

Multigenic phylogeny and analysis of tree incongruences in Triticeae (Poaceae)

Background

Introgressive events (e.g., hybridization, gene flow, horizontal gene transfer) and incomplete lineage sorting of ancestral polymorphisms are a challenge for phylogenetic analyses since different genes may exhibit conflicting genealogical histories. Grasses of the Triticeae tribe provide a particularly striking example of incongruence among gene trees. Previous phylogenies, mostly inferred with one gene, are in conflict for several taxon positions. Therefore, obtaining a resolved picture of relationships among genera and species of this tribe has been a challenging task. Here, we obtain the most comprehensive molecular dataset to date in Triticeae, including one chloroplastic and 26 nuclear genes. We aim to test whether it is possible to infer phylogenetic relationships in the face of (potentially) large-scale introgressive events and/or incomplete lineage sorting; to identify parts of the evolutionary history that have not evolved in a tree-like manner; and to decipher the biological causes of gene-tree conflicts in this tribe.

Results

We obtain resolved phylogenetic hypotheses using the supermatrix and Bayesian Concordance Factors (BCF) approaches despite numerous incongruences among gene trees. These phylogenies suggest the existence of 4-5 major clades within Triticeae, with Psathyrostachys and Hordeum being the deepest genera. In addition, we construct a multigenic network that highlights parts of the Triticeae history that have not evolved in a tree-like manner. Dasypyrum, Heteranthelium and genera of clade V, grouping Secale, Taeniatherum, Triticum and Aegilops, have evolved in a reticulated manner. Their relationships are thus better represented by the multigenic network than by the supermatrix or BCF trees. Noteworthy, we demonstrate that gene-tree incongruences increase with genetic distance and are greater in telomeric than centromeric genes. Together, our results suggest that recombination is the main factor decoupling gene trees from multigenic trees.

Conclusions

Our study is the first to propose a comprehensive, multigenic phylogeny of Triticeae. It clarifies several aspects of the relationships among genera and species of this tribe, and pinpoints biological groups with likely reticulate evolution. Importantly, this study extends previous results obtained in Drosophila by demonstrating that recombination can exacerbate gene-tree conflicts in phylogenetic reconstructions.

High throughput analysis of grape genetic diversity as a tool for germplasm collection management

Using 20 SSR markers well scattered across the 19 grape chromosomes, we analyzed 4,370 accessions of the INRA grape repository at Vassal, mostly cultivars of Vitis vinifera subsp. sativa (3,727), but also accessions of V. vinifera subsp. sylvestris (80), interspecific hybrids (364), and rootstocks (199). The analysis revealed 2,836 SSR single profiles: 2,323 sativa cultivars, 72 wild individuals (sylvestris), 306 interspecific hybrids, and 135 rootstocks, corresponding to 2,739 different cultivars in all. A total of 524 alleles were detected, with a mean of 26.20 alleles per locus. For the 2,323 cultivars of V. vinifera, 338 alleles were detected with a mean of 16.9 alleles per locus. The mean genetic diversity (GDI) was 0.797 and the level of heterozygosity was 0.76, with broad variation from 0.20 to 1. Interspecific hybrids and rootstocks were more heterozygous and more diverse (GDI = 0.839 and 0.865, respectively) than V. vinifera cultivars (GDI = 0.769), Vitis vinifera subsp. sylvestris being the least divergent with GDI = 0.708. Principal coordinates analysis distinguished the four groups. Slight clonal polymorphism was detected. The limit between clonal variation and cultivar polymorphism was set at four allelic differences out of 40. SSR markers were useful as a complementary tool to traditional ampelography for cultivar identification. Finally, a set of nine SSR markers was defined that was sufficient to distinguish 99.8% of the analyzed accessions. This set is suitable for routine characterization and will be valuable for germplasm management.