Comparative analysis of expressed sequence tags of conifers and angiosperms reveals sequences specifically conserved in conifers

A protocol for Agrobacterium-mediated transformation of Japanese cedar, Sugi (Cryptomeria japonica D. Don) using embryogenic tissue explants

Sugi (Cryptomeria japonica D. Don) is the most important afforestation coniferous tree in Japan. Coniferous trees normally have a long juvenile period and require a long cultivation time for breeding. Through a traditional breeding project that began in the 1950s, first generation plus trees with excellent traits were selected primarily from artificial forests and used as seedlings. Recently, the second generation plus trees obtained by crossing between plus trees have been selected. In light of this situation, the improvement of Sugi by a transgenic approach is effective in terms of shortening the breeding period compared with conventional crossing-dependent approaches. There are three key points to an efficient Agrobacterium-mediated transformation system: (1) establishment of explants with high regeneration ability, (2) optimal co-cultivation conditions for explants and Agrobacterium, and (3) efficient elimination of Agrobacterium. Here we describe a protocol for Agrobacterium-mediated transformation of Sugi that meets the above criteria using embryogenic tissues as explants isolated from immature seeds obtained by crossing.

DNA methylome of the 20-gigabase Norway spruce genome

DNA methylation plays important roles in many biological processes, such as silencing of transposable elements, imprinting, and regulating gene expression. Many studies of DNA methylation have shown its essential roles in angiosperms (flowering plants). However, few studies have examined the roles and patterns of DNA methylation in gymnosperms. Here, we present genome-wide high coverage single-base resolution methylation maps of Norway spruce (Picea abies) from both needles and somatic embryogenesis culture cells via whole genome bisulfite sequencing. On average, DNA methylation levels of CG and CHG of Norway spruce were higher than most other plants studied. CHH methylation was found at a relatively low level; however, at least one copy of most of the RNA-directed DNA methylation pathway genes was found in Norway spruce, and CHH methylation was correlated with levels of siRNAs. In comparison with needles, somatic embryogenesis culture cells that are used for clonally propagating spruce trees showed lower levels of CG and CHG methylation but higher level of CHH methylation, suggesting that like in other species, these culture cells show abnormal methylation patterns.

Transcriptome profiling and in silico analysis of somatic embryos in Japanese larch (Larix leptolepis)

Japanese larch (Larix leptolepis) is an ecologically and economically important species mainly grown in northeastern China, Japan and Europe. However, erratic flowering and poor germplasm resources caused by high embryo abortion rates have hampered breeding of Larix species. Somatic embryogenesis (SE) is an effective tool for the production of L. leptolepis with desirable characteristics, such as expression of totipotency, preparation of synthetic seeds, and genetic transformation. However, public genomic resources for this species are limited. We sequenced 591,759 raw expressed sequence tags (ESTs) from a 454 sequencing cDNA library of L. leptolepis somatic embryos, resulting in 572,403 high-quality reads. These reads were assembled into 70,927 unique sequences (UniGenes), including 32,321 contigs and 38,606 singletons. After removal of low-quality sequences, 65,115 UniGenes were annotated using the UniProtKB program. Based on their sequence similarity with known proteins, the matched 30,372 sequences from 664 species were estimated to represent approximately 19,000 unique genes. Gene ontology analysis revealed 21,324 UniGenes assigned to 51 categories. By Kyoto Encyclopedia of Genes and Genomes mapping, 25,773 transcripts were associated with 160 biochemical pathways. Further analysis screened four signal transduction pathways represented by 337 enzymes and 17 secondary metabolites. In silico analysis reveals that 207 UniESTs in Larix are homologous to MAPKs genes identified from other model plants, which may be involved in regulating SE development. This study provides an initial insight into the Larix transcriptomes of the pro-embryogenic mass and is a sound basis for future studies.

KEY MESSAGE

We constructed a large, full-length 454 sequencing cDNA library of Larix leptolepis during somatic embryogenesis. More than 590,000 sequences were obtained and a deep-coverage EST database was constructed.

A second generation framework for the analysis of microsatellites in expressed sequence tags and the development of EST-SSR markers for a conifer, Cryptomeria japonica

Background

Microsatellites or simple sequence repeats (SSRs) in expressed sequence tags (ESTs) are useful resources for genome analysis because of their abundance, functionality and polymorphism. The advent of commercial second generation sequencing machines has lead to new strategies for developing EST-SSR markers, necessitating the development of bioinformatic framework that can keep pace with the increasing quality and quantity of sequence data produced. We describe an open scheme for analyzing ESTs and developing EST-SSR markers from reads collected by Sanger sequencing and pyrosequencing of sugi (Cryptomeria japonica).

Results

We collected 141,097 sequence reads by Sanger sequencing and 1,333,444 by pyrosequencing. After trimming contaminant and low quality sequences, 118,319 Sanger and 1,201,150 pyrosequencing reads were passed to the MIRA assembler, generating 81,284 contigs that were analysed for SSRs. 4,059 SSRs were found in 3,694 (4.54%) contigs, giving an SSR frequency lower than that in seven other plant species with gene indices (5.4–21.9%). The average GC content of the SSR-containing contigs was 41.55%, compared to 40.23% for all contigs. Tri-SSRs were the most common SSRs; the most common motif was AT, which was found in 655 (46.3%) di-SSRs, followed by the AAG motif, found in 342 (25.9%) tri-SSRs. Most (72.8%) tri-SSRs were in coding regions, but 55.6% of the di-SSRs were in non-coding regions; the AT motif was most abundant in 3′ untranslated regions. Gene ontology (GO) annotations showed that six GO terms were significantly overrepresented within SSR-containing contigs. Forty–four EST-SSR markers were developed from 192 primer pairs using two pipelines: read2Marker and the newly-developed CMiB, which combines several open tools. Markers resulting from both pipelines showed no differences in PCR success rate and polymorphisms, but PCR success and polymorphism were significantly affected by the expected PCR product size and number of SSR repeats, respectively. EST-SSR markers exhibited less polymorphism than genomic SSRs.

Conclusions

We have created a new open pipeline for developing EST-SSR markers and applied it in a comprehensive analysis of EST-SSRs and EST-SSR markers in C. japonica. The results will be useful in genomic analyses of conifers and other non-model species.

The construction of a high-density linkage map for identifying SNP markers that are tightly linked to a nuclear-recessive major gene for male sterility in Cryptomeria japonica D. Don

BACKGROUND

High-density linkage maps facilitate the mapping of target genes and the construction of partial linkage maps around target loci to develop markers for marker-assisted selection (MAS). MAS is quite challenging in conifers because of their large, complex, and poorly-characterized genomes. Our goal was to construct a high-density linkage map to facilitate the identification of markers that are tightly linked to a major recessive male-sterile gene (ms1) for MAS in C. japonica, a species that is important in Japanese afforestation but which causes serious social pollinosis problems.

RESULTS

We constructed a high-density saturated genetic linkage map for C. japonica using expressed sequence-derived co-dominant single nucleotide polymorphism (SNP) markers, most of which were genotyped using the GoldenGate genotyping assay. A total of 1261 markers were assigned to 11 linkage groups with an observed map length of 1405.2 cM and a mean distance between two adjacent markers of 1.1 cM; the number of linkage groups matched the basic chromosome number in C. japonica. Using this map, we located ms1 on the 9th linkage group and constructed a partial linkage map around the ms1 locus. This enabled us to identify a marker (hrmSNP970_sf) that is closely linked to the ms1 gene, being separated from it by only 0.5 cM.

CONCLUSIONS

Using the high-density map, we located the ms1 gene on the 9th linkage group and constructed a partial linkage map around the ms1 locus. The map distance between the ms1 gene and the tightly linked marker was only 0.5 cM. The identification of markers that are tightly linked to the ms1 gene will facilitate the early selection of male-sterile trees, which should expedite C. japonica breeding programs aimed at alleviating pollinosis problems without harming productivity.

The first insight into the Taxus genome via fosmid library construction and end sequencing

Taxus mairei is a critically endangered and commercially important cultured medicinal gymnosperm in China and forms an important medicinal resource, but the research of its genome is absent. In this study, we constructed a T. mairei fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of one million clones with an average insert size of about 39 kb, amounting to 3.9 genome equivalents. Fosmid stability assays indicate that T. mairei DNA was stable during propagation in the fosmid system. End sequencing of both 5' and 3' ends of 968 individual clones generated 1,923 sequences after trimming, with an average sequence length of 839 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 560 (29.1%) significant hits (E < e(-5)). Repetitive sequences analysis revealed that 20.8% of end sequences are repetitive elements, which were composed of retroelements, DNA transposons, satellites, simple repeats, and low complexity sequences. The distribution pattern of various repeat types was found to be more similar to the gymnosperm Pinus and Picea than to the monocot and dicot. The satellites of T. mairei were significantly longer than those of P. taeda and P. glauca. The tetra-nucleotide repeats of T. mairei were much longer than those of P. glauca and P. taeda. The fosmid library and the fosmid end sequences, for the first time, will serve as a useful resource for large-scale genome sequencing, physical mapping, SSR marker development and positional cloning, and provide a better understanding of the Taxus genome.

Multilocus patterns of nucleotide polymorphism and demographic change in Taxodium distichum (Cupressaceae) in the lower Mississippi River alluvial valley

PREMISE OF THE STUDY

Studies of the geographic patterns of genetic variation can give important insights into the past population structure of species. Our study species, Taxodium distichum L. (bald-cypress), prefers riparian and wetland habitats and is widely distributed in southeastern North America and Mexico. We compared the genetic variation of T. distichum with that of its close relative, Cryptomeria japonica, which is endemic to Japan. •

METHODS

Nucleotide polymorphisms of T. distichum in the lower Mississippi River alluvial valley, USA, were examined at 10 nuclear loci. •

KEY RESULTS

The average nucleotide diversity at silent sites, π(sil), across the 10 loci in T. distichum was higher than that of C. japonica (π(sil) = 0.00732 and 0.00322, respectively). In T. distichum, Tajima's D values were each negative at 9 out of 10 loci, which suggests a recent population expansion. Maximum-likelihood and Bayesian estimations of the exponential population growth rate (g) of T. distichum populations indicated that this species had expanded approximately at the rate of 1.7-1.0 × 10(-6) per year in the past. •

CONCLUSIONS

Taxodium distichum had significantly higher nucleotide variation than C. japonica, and its patterns of polymorphism contrasted strikingly with those of the latter, which previously has been inferred to have experienced a reduction in population size.

Comparative genomics reveals conservative evolution of the xylem transcriptome in vascular plants

BACKGROUND

Wood is a valuable natural resource and a major carbon sink. Wood formation is an important developmental process in vascular plants which played a crucial role in plant evolution. Although genes involved in xylem formation have been investigated, the molecular mechanisms of xylem evolution are not well understood. We use comparative genomics to examine evolution of the xylem transcriptome to gain insights into xylem evolution.

RESULTS

The xylem transcriptome is highly conserved in conifers, but considerably divergent in angiosperms. The functional domains of genes in the xylem transcriptome are moderately to highly conserved in vascular plants, suggesting the existence of a common ancestral xylem transcriptome. Compared to the total transcriptome derived from a range of tissues, the xylem transcriptome is relatively conserved in vascular plants. Of the xylem transcriptome, cell wall genes, ancestral xylem genes, known proteins and transcription factors are relatively more conserved in vascular plants. A total of 527 putative xylem orthologs were identified, which are unevenly distributed across the Arabidopsis chromosomes with eight hot spots observed. Phylogenetic analysis revealed that evolution of the xylem transcriptome has paralleled plant evolution. We also identified 274 conifer-specific xylem unigenes, all of which are of unknown function. These xylem orthologs and conifer-specific unigenes are likely to have played a crucial role in xylem evolution.

CONCLUSIONS

Conifers have highly conserved xylem transcriptomes, while angiosperm xylem transcriptomes are relatively diversified. Vascular plants share a common ancestral xylem transcriptome. The xylem transcriptomes of vascular plants are more conserved than the total transcriptomes. Evolution of the xylem transcriptome has largely followed the trend of plant evolution.

An efficient method for developing SNP markers based on EST data combined with high resolution melting (HRM) analysis

Background

In order to identify single nucleotide polymorphisms (SNPs) efficiently from a species with a large genome, SNPs were identified from an expressed sequence tag (EST) database combined with High Resolution Melting (HRM) analysis.

Findings

A total of 574 sequence tagged sites (STSs) were generated from Cryptomeria japonica and HRM analysis was used to screen for polymorphisms in these STS markers. STSs were designed in two ways: 1) putative SNP sites were identified by comparing ESTs from specific contigs, then 226 primer pairs designed for the purpose to amplify these SNPs; 2) 348 primer pairs were randomly designed using reads from the 3' end of cDNA. HRM analysis revealed that 325 markers among eight individuals were polymorphic, and that STSs, including putative SNP sites, exhibited higher levels of polymorphism.

Conclusion

Our results indicate that the combination of SNP screening from an EST database combined with HRM analysis is a highly efficient way to develop SNP markers for expressed genes. This method will contribute to both genetic mapping and the identification of SNPs in non-model organisms.

The presence of sinapyl lignin in Ginkgo biloba cell cultures changes our views of the evolution of lignin biosynthesis

Suspension cell cultures (SCCs) from one of the oldest seed plants, Ginkgo biloba, show unpredictable alterations in the nature of the lignins, such as is the recruitment of sinapyl alcohol for lignin biosynthesis, compared with the woody tissues of the same species, which lack syringyl (S) lignins. These results show that, in this gymnosperm, the genes involved in sinapyl alcohol biosynthesis are latent and that their regulatory regions respond, by initiating gene expression, to the developmental signals and the environmental clues, which condition its in vitro culture. G. biloba SCCs not only synthesize S lignins but also their extracellular proteome contains both class III peroxidases capable of oxidizing sinapyl alcohol and enzymes involved in H2O2 production, observation which suggests that the peroxidase branch for the oxidative coupling of sinapyl alcohol units into lignins is operative. The incomplete knowledge of the G. biloba peroxidase-encoding genes led us to purify, characterize and partially sequence the peroxidase responsible for monolignol oxidation. When the major peroxidase from G. biloba SCCs (GbPrx) was purified to homogeneity, it showed absorption maxima in the visible region at 414 (Soret band), and at 543 and 570 nm, which calls to mind those shown by low-spin ferric peroxidases. However, the results also showed that the paraperoxidase-like character of GbPrx is not an obstacle for oxidizing the three monolignols compared with high-spin ferric peroxidases. Taken together, these results mean that the time at which the evolutionary gain of the segment of the route that leads to the biosynthesis of S lignins took place in seed plants needs to be revised.

A conifer genomics resource of 200,000 spruce (Picea spp.) ESTs and 6,464 high-quality, sequence-finished full-length cDNAs for Sitka spruce (Picea sitchensis)

Background

Members of the pine family (Pinaceae), especially species of spruce (Picea spp.) and pine (Pinus spp.), dominate many of the world's temperate and boreal forests. These conifer forests are of critical importance for global ecosystem stability and biodiversity. They also provide the majority of the world's wood and fiber supply and serve as a renewable resource for other industrial biomaterials. In contrast to angiosperms, functional and comparative genomics research on conifers, or other gymnosperms, is limited by the lack of a relevant reference genome sequence. Sequence-finished full-length (FL)cDNAs and large collections of expressed sequence tags (ESTs) are essential for gene discovery, functional genomics, and for future efforts of conifer genome annotation.

Results

As part of a conifer genomics program to characterize defense against insects and adaptation to local environments, and to discover genes for the production of biomaterials, we developed 20 standard, normalized or full-length enriched cDNA libraries from Sitka spruce (P. sitchensis), white spruce (P. glauca), and interior spruce (P. glauca-engelmannii complex). We sequenced and analyzed 206,875 3'- or 5'-end ESTs from these libraries, and developed a resource of 6,464 high-quality sequence-finished FLcDNAs from Sitka spruce. Clustering and assembly of 147,146 3'-end ESTs resulted in 19,941 contigs and 26,804 singletons, representing 46,745 putative unique transcripts (PUTs). The 6,464 FLcDNAs were all obtained from a single Sitka spruce genotype and represent 5,718 PUTs.

Conclusion

This paper provides detailed annotation and quality assessment of a large EST and FLcDNA resource for spruce. The 6,464 Sitka spruce FLcDNAs represent the third largest sequence-verified FLcDNA resource for any plant species, behind only rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), and the only substantial FLcDNA resource for a gymnosperm. Our emphasis on capturing FLcDNAs and ESTs from cDNA libraries representing herbivore-, wound- or elicitor-treated induced spruce tissues, along with incorporating normalization to capture rare transcripts, resulted in a rich resource for functional genomics and proteomics studies. Sequence comparisons against five plant genomes and the non-redundant GenBank protein database revealed that a substantial number of spruce transcripts have no obvious similarity to known angiosperm gene sequences. Opportunities for future applications of the sequence and clone resources for comparative and functional genomics are discussed.

Characterization of expressed sequence tags from a full-length enriched cDNA library of Cryptomeria japonica male strobili

Background

Cryptomeria japonica D. Don is one of the most commercially important conifers in Japan. However, the allergic disease caused by its pollen is a severe public health problem in Japan. Since large-scale analysis of expressed sequence tags (ESTs) in the male strobili of C. japonica should help us to clarify the overall expression of genes during the process of pollen development, we constructed a full-length enriched cDNA library that was derived from male strobili at various developmental stages.

Results

We obtained 36,011 expressed sequence tags (ESTs) from either one or both ends of 19,437 clones derived from the cDNA library of C. japonica male strobili at various developmental stages. The 19,437 cDNA clones corresponded to 10,463 transcripts. Approximately 80% of the transcripts resembled ESTs from Pinus and Picea, while approximately 75% had homologs in Arabidopsis. An analysis of homologies between ESTs from C. japonica male strobili and known pollen allergens in the Allergome Database revealed that products of 180 transcripts exhibited significant homology. Approximately 2% of the transcripts appeared to encode transcription factors. We identified twelve genes for MADS-box proteins among these transcription factors. The twelve MADS-box genes were classified as DEF/GLO/GGM13-, AG-, AGL6-, TM3- and TM8-like MIKC^C genes and type I MADS-box genes.

Conclusion

Our full-length enriched cDNA library derived from C. japonica male strobili provides information on expression of genes during the development of male reproductive organs. We provided potential allergens in C. japonica. We also provided new information about transcription factors including MADS-box genes expressed in male strobili of C. japonica. Large-scale gene discovery using full-length cDNAs is a valuable tool for studies of gymnosperm species.

Genome scan to detect genetic structure and adaptive genes of natural populations of Cryptomeria japonica

We investigated 29 natural populations of Cryptomeria japonica using 148 cleaved amplified polymorphic sequence markers to elucidate their genetic structure and identify candidate adaptive genes of this species. In accordance with the inferred evolutionary history of the species during and after the last glacial episode, the genetic diversity was higher in western populations than in northern populations. The results of phylogenetic and genetic structure analyses suggest that populations of the two main varieties of the species have clearly diverged from each other and that two of the examined loci are strongly associated with the differentiation between the two varieties. Using a coalescent simulation based on F(ST) and H(e) values, we detected five genes that had higher, and two that had lower, values than the respective 99% confidence intervals (C.I.s) that are theoretically expected intervals under a neutral infinite-island model. We also detected 13 outlier loci using a coalescent simulation based on the assumption that the 2 varieties originated from the splitting of an ancestral population. Four of these loci were detected by both methods, two of which were detected in a genetic structure analysis as loci associated with differentiation between the two varieties of the species, and are strong candidates for genes that have been subject to selection.

Genetic diversity and the genetic structure of natural populations of Chamaecyparis obtusa: implications for management and conservation

We investigated 25 natural populations of Chamaecyparis obtusa using 51 cleaved amplified polymorphic sequence (CAPS) markers, which were developed using information on sequence-tagged sites (STS) in Cryptomeria japonica. Most CAPS markers have codominant expression patterns, and are suitable for population studies because of their robustness and convenience. We estimated various genetic diversity parameters, including average heterozygosity (H(e)) and allelic richness and found that the more peripheral populations tended to have lower genetic diversity than central populations, in agreement with a previous theoretical study. The overall genetic differentiation between populations was low, but statistically significant (G(ST)=0.039), and similar to the level reported in a previous allozyme study. We attempted to detect non-neutral loci associated with local adaptation to clarify the relationship between the fixation index (F(ST)) and H(e) values for each locus and found seven candidates non-neutral loci. Phylogenetic tree analysis of the populations and Bayesian clustering analysis revealed a pattern of gradually increasing isolation of populations with increasing geographical distance. Three populations had a high degree of linkage disequilibrium, which we attribute to severe bottlenecks due to human disturbance or competition with other species during their migration from refugia after the most recent glaciation. We concluded that the small populations in western Japan and in Kanto district are more important, from a conservation perspective, than the populations in central Japan, due to their genetic divergence, relatively small sizes and restricted areas.

Structural motifs of syringyl peroxidases predate not only the gymnosperm-angiosperm divergence but also the radiation of tracheophytes

* The most distinctive variation in the monomer composition of lignins in vascular land plants is that found between the two main groups of seed plants. Thus, while gymnosperm lignins are typically composed of guaiacyl (G) units, angiosperm lignins are largely composed of similar levels of G and syringyl (S) units. * However, and contrary to what might be expected, peroxidases isolated from basal (Cycadales and Ginkgoales) and differentially evolved (Coniferales and Gnetales) gymnosperms are also able to oxidize S moieties, and this ability is independent of the presence or absence of S-type units in their lignins. * The results obtained led us to look at the protein database to search for homologies between gymnosperm peroxidases and true eudicot S-peroxidases, such as the Zinnia elegans peroxidase. * The findings showed that certain structural motifs characteristic of eudicot S-peroxidases (certain amino acid sequences and beta-sheet secondary structures) predate the gymnosperm-angiosperm divergence and the radiation of tracheophytes, since they are found not only in peroxidases from basal gymnosperms, ferns and lycopods, but also in peroxidases from the moss Physcomitrella patens (Bryopsida) and the liverwort Marchantia polymorpha (Marchantiopsida), which, as typical of bryophytes, do not have xylem tissue nor lignins.

Expressed sequence tags from loblolly pine embryos reveal similarities with angiosperm embryogenesis

The process of embryogenesis in gymnosperms differs in significant ways from the more widely studied process in angiosperms. To further our understanding of embryogenesis in gymnosperms, we have generated Expressed Sequence Tags (ESTs) from four cDNA libraries constructed from un-normalized, normalized, and subtracted RNA populations of zygotic and somatic embryos of loblolly pine (Pinus taeda L.). A total of 68,721 ESTs were generated from 68,131 cDNA clones. Following clustering and assembly, these sequences collapsed into 5,274 contigs and 6,880 singleton sequences for a total of 12,154 non-redundant sequences. Searches of a non-identical amino acid database revealed a putative homolog for 9,189 sequences, leaving 2,965 sequences with no known function. More extensive searches of additional plant sequence data sets revealed a putative homolog for all but 1,388 (11.4%) of the sequences. Using gene ontologies, a known function could be assigned for 5,495 of the 12,154 total non-redundant sequences with 13,633 associations in total assigned. When compared to approximately 72,000 sequences in a collated P. taeda transcript assembly derived from >245,000 ESTs derived from root, xylem, stem, needles, pollen cone, and shoot ESTs, 3,458 (28.5%) of the non-redundant embryo sequences were unique and thereby provide a valuable addition to development of a complete loblolly pine transcriptome. To assess similarities between angiosperm and gymnosperm embryo development, we examined our EST collection for putative homologs of angiosperm genes implicated in embryogenesis. Out of 108 angiosperm embryogenesis-related genes, homologs were present for 83 of these genes suggesting that pine contains similar genes for embryogenesis and that our RNA sampling methods were successful. We also identified sequences from the pine embryo transcriptome that have no known function and may contribute to the programming of gene expression and embryo development.