Quantifying the effect of human practices on S. cerevisiae vineyard metapopulation diversity.

Saccharomyces cerevisiae is the main actor of wine fermentation but at present, still little is known about the factors impacting its distribution in the vineyards. In this study, 23 vineyards and 7 cellars were sampled over 2 consecutive years in the Bordeaux and Bergerac regions. The impact of geography and farming system and the relation between grape and vat populations were evaluated using a collection of 1374 S. cerevisiae merlot grape isolates and 289 vat isolates analyzed at 17 microsatellites loci. A very high genetic diversity of S. cerevisiae strains was obtained from grape samples, higher in conventional farming system than in organic one. The geographic appellation and the wine estate significantly impact the S. cerevisiae population structure, whereas the type of farming system has a weak global effect. When comparing cellar and vineyard populations, we evidenced the tight connection between the two compartments, based on the high proportion of grape isolates (25%) related to the commercial starters used in the cellar and on the estimation of bidirectional geneflows between the vineyard and the cellar compartments.

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.

CReasPy-Cloning: A Method for Simultaneous Cloning and Engineering of Megabase-Sized Genomes in Yeast Using the CRISPR-Cas9 System.

Over the past decade, a new strategy was developed to bypass the difficulties to genetically engineer some microbial species by transferring (or "cloning") their genome into another organism that is amenable to efficient genetic modifications and therefore acts as a living workbench. As such, the yeast Saccharomyces cerevisiae has been used to clone and engineer genomes from viruses, bacteria, and algae. The cloning step requires the insertion of yeast genetic elements in the genome of interest, in order to drive its replication and maintenance as an artificial chromosome in the host cell. Current methods used to introduce these genetic elements are still unsatisfactory, due either to their random nature (transposon) or the requirement for unique restriction sites at specific positions (TAR cloning). Here we describe the CReasPy-cloning, a new method that combines both the ability of Cas9 to cleave DNA at a user-specified locus and the yeast's highly efficient homologous recombination to simultaneously clone and engineer a bacterial chromosome in yeast. Using the 0.816 Mbp genome of Mycoplasma pneumoniae as a proof of concept, we demonstrate that our method can be used to introduce the yeast genetic element at any location in the bacterial chromosome while simultaneously deleting various genes or group of genes. We also show that CReasPy-cloning can be used to edit up to three independent genomic loci at the same time with an efficiency high enough to warrant the screening of a small (<50) number of clones, allowing for significantly shortened genome engineering cycle times.

Correction for Rosa et al., "Draft Genome Sequence of the Candida zemplinina (syn., Starmerella bacillaris) Type Strain CBS 9494".

[This corrects the article DOI: 10.1128/MRA.00872-18.].

Draft Genome Sequence of the Candida zemplinina (syn., Starmerella bacillaris) Type Strain CBS 9494 [corrected].

Starmerella bacillaris is an ascomycetous yeast ubiquitously present in grapes and fermenting grape musts. In this report, we present the draft genome sequence of the S. bacillaris type strain CBS 9494, isolated from sweet botrytized wines, which will contribute to the study of this genetically heterogeneous wine yeast species.

MiR-10a and HOXB4 are overexpressed in atypical myeloproliferative neoplasms.

BACKGROUND: Atypical Myeloproliferative Neoplasms (aMPN) share characteristics of MPN and Myelodysplastic Syndromes. Although abnormalities in cytokine signaling are common in MPN, the pathophysiology of atypical MPN still remains elusive. Since deregulation of microRNAs is involved in the biology of various cancers, we studied the miRNome of aMPN patients. METHODS: MiRNome and mutations in epigenetic regulator genes ASXL1, TET2, DNMT3A, EZH2 and IDH1/2 were explored in aMPN patients. Epigenetic regulation of miR-10a and HOXB4 expression was investigated by treating hematopoietic cell lines with 5-aza-2'deoxycytidine, valproic acid and retinoic acid. Functional effects of miR-10a overexpression on cell proliferation, differentiation and self-renewal were studied by transducing CD34+ cells with lentiviral vectors encoding the pri-miR-10a precursor. RESULTS: MiR-10a was identified as the most significantly up-regulated microRNA in aMPN. MiR-10a expression correlated with that of HOXB4, sitting in the same genomic locus. The transcription of these two genes was increased by DNA demethylation and histone acetylation, both necessary for optimal expression induction by retinoic acid. Moreover, miR-10a and HOXB4 overexpression seemed associated with DNMT3A mutation in hematological malignancies. However, overexpression of miR-10a had no effect on proliferation, differentiation or self-renewal of normal hematopoietic progenitors. CONCLUSIONS: MiR-10a and HOXB4 are overexpressed in aMPN. This overexpression seems to be the result of abnormalities in epigenetic regulation mechanisms. Our data suggest that miR-10a could represent a simple marker of transcription at this genomic locus including HOXB4, widely recognized as involved in stem cell expansion.

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.

Brettanomyces bruxellensis population survey reveals a diploid-triploid complex structured according to substrate of isolation and geographical distribution.

Brettanomyces bruxellensis is a unicellular fungus of increasing industrial and scientific interest over the past 15 years. Previous studies revealed high genotypic diversity amongst B. bruxellensis strains as well as strain-dependent phenotypic characteristics. Genomic assemblies revealed that some strains harbour triploid genomes and based upon prior genotyping it was inferred that a triploid population was widely dispersed across Australian wine regions. We performed an intraspecific diversity genotypic survey of 1488 B. bruxellensis isolates from 29 countries, 5 continents and 9 different fermentation niches. Using microsatellite analysis in combination with different statistical approaches, we demonstrate that the studied population is structured according to ploidy level, substrate of isolation and geographical origin of the strains, underlying the relative importance of each factor. We found that geographical origin has a different contribution to the population structure according to the substrate of origin, suggesting an anthropic influence on the spatial biodiversity of this microorganism of industrial interest. The observed clustering was correlated to variable stress response, as strains from different groups displayed variation in tolerance to the wine preservative sulfur dioxide (SO2). The potential contribution of the triploid state for adaptation to industrial fermentations and dissemination of the species B. bruxellensis is discussed.

Modulation of oncogenic miRNA biogenesis using functionalized polyamines.

MicroRNAs are key factors in the regulation of gene expression and their deregulation has been directly linked to various pathologies such as cancer. The use of small molecules to tackle the overexpression of oncogenic miRNAs has proved its efficacy and holds the promise for therapeutic applications. Here we describe the screening of a 640-compound library and the identification of polyamine derivatives interfering with in vitro Dicer-mediated processing of the oncogenic miR-372 precursor (pre-miR-372). The most active inhibitor is a spermine-amidine conjugate that binds to the pre-miR-372 with a KD of 0.15 µM, and inhibits its in vitro processing with a IC50 of 1.06 µM. The inhibition of miR-372 biogenesis was confirmed in gastric cancer cells overexpressing miR-372 and a specific inhibition of proliferation through de-repression of the tumor suppressor LATS2 protein, a miR-372 target, was observed. This compound modifies the expression of a small set of miRNAs and its selective biological activity has been confirmed in patient-derived ex vivo cultures of gastric carcinoma. Polyamine derivatives are promising starting materials for future studies about the inhibition of oncogenic miRNAs and, to the best of our knowledge, this is the first report about the application of functionalized polyamines as miRNAs interfering agents.

Biogeography of Oenococcus oeni Reveals Distinctive but Nonspecific Populations in Wine-Producing Regions.

Understanding the mechanisms behind the typicity of regional wines inevitably brings attention to microorganisms associated with their production. Oenococcus oeni is the main bacterial species involved in wine and cider making. It develops after the yeast-driven alcoholic fermentation and performs the malolactic fermentation, which improves the taste and aromatic complexity of most wines. Here, we have evaluated the diversity and specificity of O. oeni strains in six regions. A total of 235 wines and ciders were collected during spontaneous malolactic fermentations and used to isolate 3,212 bacterial colonies. They were typed by multilocus variable analysis, which disclosed a total of 514 O. oeni strains. Their phylogenetic relationships were evaluated by a second typing method based on single nucleotide polymorphism (SNP) analysis. Taken together, the results indicate that each region holds a high diversity of strains that constitute a unique population. However, strains present in each region belong to diverse phylogenetic groups, and the same groups can be detected in different regions, indicating that strains are not genetically adapted to regions. In contrast, greater strain identity was seen for cider, white wine, or red wine of Burgundy, suggesting that genetic adaptation to these products occurred. IMPORTANCE: This study reports the isolation, genotyping, and geographic distribution analysis of the largest collection of O. oeni strains performed to date. It reveals that there is very high diversity of strains in each region, the majority of them being detected in a single region. The study also reports the development of an SNP genotyping method that is useful for analyzing the distribution of O. oeni phylogroups. The results show that strains are not genetically adapted to regions but to specific types of wines. They reveal new phylogroups of strains, particularly two phylogroups associated with white wines and red wines of Burgundy. Taken together, the results shed light on the diversity and specificity of wild strains of O. oeni, which is crucial for understanding their real contribution to the unique properties of wines.

Cellar-Associated Saccharomyces cerevisiae Population Structure Revealed High-Level Diversity and Perennial Persistence at Sauternes Wine Estates.

UNLABELLED: Three wine estates (designated A, B, and C) were sampled in Sauternes, a typical appellation of the Bordeaux wine area producing sweet white wine. From those wine estates, 551 yeast strains were collected between 2012 and 2014, added to 102 older strains from 1992 to 2011 from wine estate C. All the strains were analyzed through 15 microsatellite markers, resulting in 503 unique Saccharomyces cerevisiae genotypes, revealing high genetic diversity and a low presence of commercial yeast starters. Population analysis performed using Fst genetic distance or ancestry profiles revealed that the two closest wine estates, B and C, which have juxtaposed vineyard plots and common seasonal staff, share more related isolates with each other than with wine estate A, indicating exchange between estates. The characterization of isolates collected 23 years ago at wine estate C in relation to recent isolates obtained at wine estate B revealed the long-term persistence of isolates. Last, during the 2014 harvest period, a temporal succession of ancestral subpopulations related to the different batches associated with the selective picking of noble rotted grapes was highlighted. IMPORTANCE: High genetic diversity of S. cerevisiae isolates from spontaneous fermentation on wine estates in the Sauternes appellation of Bordeaux was revealed. Only 7% of all Sauternes strains were considered genetically related to specific commercial strains. The long-term persistence (over 20 years) of S. cerevisiae profiles on a given wine estate is highlighted.

Microsatellite analysis of Saccharomyces uvarum diversity.

Considered as a sister species of Saccharomyces cerevisiae, S. uvarum is, to a lesser extent, an interesting species for fundamental and applied research studies. Despite its potential interest as a new gene pool for fermenting agents, the intraspecific molecular genetic diversity of this species is still poorly investigated. In this study, we report the use of nine microsatellite markers to describe S. uvarum genetic diversity and population structure among 108 isolates from various geographical and substrate origins (wine, cider and natural sources). Our combined microsatellite markers set allowed differentiating 89 genotypes. In contrast to S. cerevisiae genetic diversity, wild and human origin isolates were intertwined. A total of 75% of strains were proven to be homozygotes and estimated heterozygosity suggests a selfing rate above 0.95 for the different population tested here. From this point of view, the S. uvarum life cycle appears to be more closely related to S. paradoxus or S. cerevisiae of natural resources than S. cerevisiae wine isolates. Population structure could not be correlated to distinct geographic or technological origins, suggesting lower differentiation that may result from a large exchange between human and natural populations mediated by insects or human activities.

Decoding the oak genome: public release of sequence data, assembly, annotation and publication strategies.

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole-genome shotgun (WGS) approach, without the use of costly and time-consuming methods, such as fosmid or BAC clone-based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS-FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired-end and mate-pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired-end reads and contaminants detected, resulting in a total of 17,910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.

The yeast Starmerella bacillaris (synonym Candida zemplinina) shows high genetic diversity in winemaking environments.

The yeast Candida zemplinina (Starmerella bacillaris) is frequently isolated from grape and wine environments. Its enological use in mixed fermentation with Saccharomyces cerevisiae has been extensively investigated these last few years, and several interesting features including low ethanol production, fructophily, glycerol and other metabolites production, have been described. In addition, molecular tools allowing the characterization of yeast populations have been developed, both at the inter- and intraspecific levels. However, most of these fingerprinting methods are not compatible with population genetics or ecological studies. In this work, we developed 10 microsatellite markers for the C. zemplinina species that were used for the genotyping of 163 strains from nature or various enological regions (28 vineyards/wineries from seven countries). We show that the genetic diversity of C. zemplinina is shaped by geographical localization. Populations isolated from winemaking environments are quite diverse at the genetic level: neither clonal-like behaviour nor specific genetic signature were associated with the different vineyards/wineries. Altogether, these results suggest that C. zemplinina is not under selective pressure in winemaking environments.

Hanseniaspora uvarum from Winemaking Environments Show Spatial and Temporal Genetic Clustering.

Hanseniaspora uvarum is one of the most abundant yeast species found on grapes and in grape must, at least before the onset of alcoholic fermentation (AF) which is usually performed by Saccharomyces species. The aim of this study was to characterize the genetic and phenotypic variability within the H. uvarum species. One hundred and fifteen strains isolated from winemaking environments in different geographical origins were analyzed using 11 microsatellite markers and a subset of 47 strains were analyzed by AFLP. H. uvarum isolates clustered mainly on the basis of their geographical localization as revealed by microsatellites. In addition, a strong clustering based on year of isolation was evidenced, indicating that the genetic diversity of H. uvarum isolates was related to both spatial and temporal variations. Conversely, clustering analysis based on AFLP data provided a different picture with groups showing no particular characteristics, but provided higher strain discrimination. This result indicated that AFLP approaches are inadequate to establish the genetic relationship between individuals, but allowed good strain discrimination. At the phenotypic level, several extracellular enzymatic activities of enological relevance (pectinase, chitinase, protease, ß-glucosidase) were measured but showed low diversity. The impact of environmental factors of enological interest (temperature, anaerobia, and copper addition) on growth was also assessed and showed poor variation. Altogether, this work provided both new analytical tool (microsatellites) and new insights into the genetic and phenotypic diversity of H. uvarum, a yeast species that has previously been identified as a potential candidate for co-inoculation in grape must, but whose intraspecific variability had never been fully assessed.

Development of microsatellite markers for the rapid and reliable genotyping of Brettanomyces bruxellensis at strain level.

Although many yeasts are useful for food production and beverage, some species may cause spoilage with important economic loss. This is the case of Dekkera/Brettanomyces bruxellensis, a contaminant species that is mainly associated with fermented beverages (wine, beer, cider and traditional drinks). To better control Brettanomyces spoilage, rapid and reliable genotyping methods are necessary to determine the origins of the spoilage, to assess the effectiveness of preventive treatments and to develop new control strategies. Despite several previously published typing methods, ranging from classical molecular methods (RAPD, AFLP, REA-PFGE, mtDNA restriction analysis) to more engineered technologies (infrared spectroscopy), there is still a lack of a rapid, reliable and universal genotyping approach. In this work, we developed eight polymorphic microsatellites markers for the Brettanomyces/Dekkera bruxellensis species. Microsatellite typing was applied to the genetic analysis of wine and beer isolates from Europe, Australia and South Africa. Our results suggest that B. bruxellensis is a highly disseminated species, with some strains isolated from different continents being closely related at the genetic level. We also focused on strains isolated from two Bordeaux wineries on different substrates (grapes, red wines) and for different vintages (over half a century). We showed that all B. bruxellensis strains within a cellar are strongly related at the genetic level, suggesting that one clonal population may cause spoilage over decades. The microsatellite tool now paves the way for future population genetics research of the B. bruxellensis species.

Winemaking and bioprocesses strongly shaped the genetic diversity of the ubiquitous yeast Torulaspora delbrueckii.

The yeast Torulaspora delbrueckii is associated with several human activities including oenology, bakery, distillery, dairy industry, etc. In addition to its biotechnological applications, T. delbrueckii is frequently isolated in natural environments (plant, soil, insect). T. delbrueckii is thus a remarkable ubiquitous yeast species with both wild and anthropic habitats, and appears to be a perfect yeast model to search for evidence of human domestication. For that purpose, we developed eight microsatellite markers that were used for the genotyping of 110 strains from various substrates and geographical origins. Microsatellite analysis showed four genetic clusters: two groups contained most nature strains from Old World and Americas respectively, and two clusters were associated with winemaking and other bioprocesses. Analysis of molecular variance (AMOVA) confirmed that human activities significantly shaped the genetic variability of T. delbrueckii species. Natural isolates are differentiated on the basis of geographical localisation, as expected for wild population. The domestication of T. delbrueckii probably dates back to the Roman Empire for winemaking (∼ 1900 years ago), and to the Neolithic era for bioprocesses (∼ 4000 years ago). Microsatellite analysis also provided valuable data regarding the life-cycle of the species, suggesting a mostly diploid homothallic life. In addition to population genetics and ecological studies, the microsatellite tool will be particularly useful for further biotechnological development of T. delbrueckii strains for winemaking and other bioprocesses.

High-density linkage mapping in a pine tree reveals a genomic region associated with inbreeding depression and provides clues to the extent and distribution of meiotic recombination.

BACKGROUND: The availability of a large expressed sequence tags (EST) resource and recent advances in high-throughput genotyping technology have made it possible to develop highly multiplexed SNP arrays for multi-objective genetic applications, including the construction of meiotic maps. Such approaches are particularly useful in species with a large genome size, precluding the use of whole-genome shotgun assembly with current technologies. RESULTS: In this study, a 12 k-SNP genotyping array was developed for maritime pine from an extensive EST resource assembled into a unigene set. The offspring of three-generation outbred and inbred mapping pedigrees were then genotyped. The inbred pedigree consisted of a classical F2 population resulting from the selfing of a single inter-provenance (Landes x Corsica) hybrid tree, whereas the outbred pedigree (G2) resulted from a controlled cross of two intra-provenance (Landes x Landes) hybrid trees. This resulted in the generation of three linkage maps based on SNP markers: one from the parental genotype of the F2 population (1,131 markers in 1,708 centimorgan (cM)), and one for each parent of the G2 population (1,015 and 1,110 markers in 1,447 and 1,425 cM for the female and male parents, respectively). A comparison of segregation patterns in the progeny obtained from the two types of mating (inbreeding and outbreeding) led to the identification of a chromosomal region carrying an embryo viability locus with a semi-lethal allele. Following selfing and segregation, zygote mortality resulted in a deficit of Corsican homozygous genotypes in the F2 population. This dataset was also used to study the extent and distribution of meiotic recombination along the length of the chromosomes and the effect of sex and/or genetic background on recombination. The genetic background of trees in which meiotic recombination occurred was found to have a significant effect on the frequency of recombination. Furthermore, only a small proportion of the recombination hot- and cold-spots were common to all three genotypes, suggesting that the spatial pattern of recombination was genetically variable. CONCLUSION: This study led to the development of classical genomic tools for this ecologically and economically important species. It also identified a chromosomal region bearing a semi-lethal recessive allele and demonstrated the genetic variability of recombination rate over the genome.

Dynamics of bacterial communities in cockles (Cerastoderma edule) with respect to trematode parasite (Bucephalus minimus) infestation.

The bacterial communities associated with the cockle (Cerastoderma edule) were investigated at the individual level through a 10-month monitoring programme. Temporal changes and those changes associated with a common parasite of the cockle, Bucephalus minimus, were investigated by monthly sampling of individuals, selected based on their shell length (cohort monitoring). Cockle bacterial community abundance (CBCA) and diversity (CBCD) were estimated by epifluorescence microscopy counts and automated ribosomal intergenic spacer analysis, respectively. CBCA showed a temporal pattern peaking at 30 × 10(6) cells per gram of cockle flesh and intervalval liquid in October and a significant 1.8-fold increase linked with B. minimus occurrence. CBCD was characterized by 112 ± 26 intergenic transcribed spacer (ITS) per individual and showed a relative homology between individuals (52 ± 6%, Jaccard similarity) in spite of more than 30% of rare ITS. Consistent with an undisturbed evolution of the condition index of the studied cohort individuals as an estimate of their physiological state, neither temporal nor parasite-induced change in CBCA has been related to marked changes in CBCD.

Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: oak.

BACKGROUND: The Fagaceae family comprises about 1,000 woody species worldwide. About half belong to the Quercus family. These oaks are often a source of raw material for biomass wood and fiber. Pedunculate and sessile oaks, are among the most important deciduous forest tree species in Europe. Despite their ecological and economical importance, very few genomic resources have yet been generated for these species. Here, we describe the development of an EST catalogue that will support ecosystem genomics studies, where geneticists, ecophysiologists, molecular biologists and ecologists join their efforts for understanding, monitoring and predicting functional genetic diversity. RESULTS: We generated 145,827 sequence reads from 20 cDNA libraries using the Sanger method. Unexploitable chromatograms and quality checking lead us to eliminate 19,941 sequences. Finally a total of 125,925 ESTs were retained from 111,361 cDNA clones. Pyrosequencing was also conducted for 14 libraries, generating 1,948,579 reads, from which 370,566 sequences (19.0%) were eliminated, resulting in 1,578,192 sequences. Following clustering and assembly using TGICL pipeline, 1,704,117 EST sequences collapsed into 69,154 tentative contigs and 153,517 singletons, providing 222,671 non-redundant sequences (including alternative transcripts). We also assembled the sequences using MIRA and PartiGene software and compared the three unigene sets. Gene ontology annotation was then assigned to 29,303 unigene elements. Blast search against the SWISS-PROT database revealed putative homologs for 32,810 (14.7%) unigene elements, but more extensive search with Pfam, Refseq_protein, Refseq_RNA and eight gene indices revealed homology for 67.4% of them. The EST catalogue was examined for putative homologs of candidate genes involved in bud phenology, cuticle formation, phenylpropanoids biosynthesis and cell wall formation. Our results suggest a good coverage of genes involved in these traits. Comparative orthologous sequences (COS) with other plant gene models were identified and allow to unravel the oak paleo-history. Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were searched, resulting in 52,834 SSRs and 36,411 SNPs. All of these are available through the Oak Contig Browser http://genotoul-contigbrowser.toulouse.inra.fr:9092/Quercus_robur/index.html. CONCLUSIONS: This genomic resource provides a unique tool to discover genes of interest, study the oak transcriptome, and develop new markers to investigate functional diversity in natural populations.