SummaryMaritime pine (Pinus pinaster Ait.) is a widely distributed conifer species in Southwestern Europe and one of the most advanced models for conifer research. In the current work, comprehensive characterization of the maritime pine transcriptome was performed using a combination of two different next-generation sequencing platforms, 454 and Illumina. De novo assembly of the transcriptome provided a catalogue of 26 020 unique transcripts in maritime pine trees and a collection of 9641 full-length cDNAs. Quality of the transcriptome assembly was validated by RT-PCR amplification of selected transcripts for structural and regulatory genes. Transcription factors and enzyme-encoding transcripts were annotated. Furthermore, the available sequencing data permitted the identification of polymorphisms and the establishment of robust single nucleotide polymorphism (SNP) and simple-sequence repeat (SSR) databases for genotyping applications and integration of translational genomics in maritime pine breeding programmes. All our data are freely available at SustainpineDB, the P. pinaster expressional database. Results reported here on the maritime pine transcriptome represent a valuable resource for future basic and applied studies on this ecological and economically important pine species.

BackgroundThe 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.ResultsIn 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.ConclusionThis 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.

Corynespora cassiicola is an Ascomycetes fungus with a broad host range and diverse life styles. Mostly known as a necrotrophic plant pathogen, it has also been associated with rare cases of human infection. In the rubber tree, this fungus causes the Corynespora leaf fall (CLF) disease, which increasingly affects natural rubber production in Asia and Africa. It has also been found as an endophyte in South American rubber plantations where no CLF outbreak has yet occurred. The C. cassiicola species is genetically highly diverse, but no clear relationship has been evidenced between phylogenetic lineage and pathogenicity. Cassiicolin, a small glycosylated secreted protein effector, is thought to be involved in the necrotrophic interaction with the rubber tree but some virulent C. cassiicola isolates do not have a cassiicolin gene. This study set out to identify other putative effectors involved in CLF. The genome of a highly virulent C. cassiicola isolate from the rubber tree (CCP) was sequenced and assembled. In silico prediction revealed 2870 putative effectors, comprising CAZymes, lipases, peptidases, secreted proteins and enzymes associated with secondary metabolism. Comparison with the genomes of 44 other fungal species, focusing on effector content, revealed a striking proximity with phylogenetically unrelated species (Colletotrichum acutatum, Colletotrichum gloesporioides, Fusarium oxysporum, nectria hematococca, and Botrosphaeria dothidea) sharing life style plasticity and broad host range. Candidate effectors involved in the compatible interaction with the rubber tree were identified by transcriptomic analysis. Differentially expressed genes included 92 putative effectors, among which cassiicolin and two other secreted singleton proteins. Finally, the genomes of 35 C. cassiicola isolates representing the genetic diversity of the species were sequenced and assembled, and putative effectors identified. At the intraspecific level, effector-based classification was found to be highly consistent with the phylogenomic trees. Identification of lineage-specific effectors is a key step toward understanding C. cassiicola virulence and host specialization mechanisms.

A novel category of major intrinsic proteins which share weak similarities with previously identified aquaporin subfamilies was recently identified in land plants, and named X (for unrecognized) intrinsic proteins (XIPs). Because XIPs are still ranked as uncharacterized proteins, their further molecular characterization is required. Herein, a systematic fine-scale analysis of XIP sequences found in flowering plant databases revealed that XIPs are found in at least five groups. The phylogenetic relationship of these five groups with the phylogenetic organization of angiosperms revealed an original pattern of evolution for the XIP subfamily through distinct angiosperm taxon-specific clades. Of all flowering plant having XIPs, the genus Populus encompasses the broadest panel and the highest polymorphism of XIP isoforms, with nine PtXIP sequences distributed within three XIP groups. Comprehensive PtXIP gene expression patterns showed that only two isoforms (PtXIP2;1 and PtXIP3;2) were transcribed in vegetative tissues. However, their patterns are contrasted, PtXIP2;1 was ubiquitously accumulated whereas PtXIP3;2 was predominantly detected in wood and to a lesser extent in roots. Furthermore, only PtXIP2;1 exhibited a differential expression in leaves and stems of drought-, salicylic acid-, or wounding-challenged plants. Unexpectedly, the PtXIPs displayed different abilities to alter water transport upon expression in Xenopus laevis oocytes. PtXIP2;1 and PtXIP3;3 transported water while other PtXIPs did not.

A global DNA methylation and proteomics approach was used to investigate somatic embryo maturation in hybrid larch. Each developmental step during somatic embryogenesis was associated with a distinct and significantly different global DNA methylation level: from 45.8% mC for undifferentiated somatic embryos (1-week proliferation) to 61.5% mC for immature somatic embryos (1-week maturation), while maturation was associated with a decrease in DNA methylation to 53.4% for mature cotyledonary somatic embryos (8-weeks maturation). The presence of 5-azacytidine (hypo-methylating agent) or hydroxyurea (hyper-methylating agent) in the maturation medium altered the global DNA methylation status of the embryogenic cultures, and significantly reduced both their relative growth rate and embryogenic potential, suggesting an important role for DNA methylation in embryogenesis. Maturation was also assessed by examining changes in the total protein profile. Storage proteins, identified as legumin- and vicilin-like, appeared at the precotyledonary stage. In the proteomic study, total soluble proteins were extracted from embryos after 1 and 8 weeks of maturation, and separated by two-dimensional gel electrophoresis. There were 147 spots which showed significant differences between the stages of maturation; they were found to be involved mainly in primary metabolism and the stabilization of the resulting metabolites. This indicated that the somatic embryo was still metabolically active at 8 weeks of maturation. This is the first report of analyses of global DNA methylation (including the effects of hyper- and hypo-treatments) and proteome during somatic embryogenesis in hybrid larch, and thus provides novel insights into maturation of conifer somatic embryos.