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

CRISPR/Cas9-mediated disruption of CjACOS5 confers no-pollen formation on sugi trees (Cryptomeria japonica D. Don)

Sugi (Cryptomeria japonica D. Don) is an economically important coniferous tree in Japan. However, abundant sugi pollen grains are dispersed and transported by the wind each spring and cause a severe pollen allergy syndrome (Japanese cedar pollinosis). The use of pollen-free sugi that cannot produce pollen has been thought as a countermeasure to Japanese cedar pollinosis. The sugi CjACOS5 gene is an ortholog of Arabidopsis ACOS5 and rice OsACOS12, which encode an acyl-CoA synthetase that is involved in the synthesis of sporopollenin in pollen walls. To generate pollen-free sugi, we mutated CjACOS5 using the CRISPR/Cas9 system. As a result of sugi transformation mediated by Agrobacterium tumefaciens harboring the CjACOS5-targeted CRISPR/Cas9 vector, 1 bp-deleted homo biallelic mutant lines were obtained. Chimeric mutant lines harboring both mutant and wild-type CjACOS5 genes were also generated. The homo biallelic mutant lines had no-pollen in male strobili, whereas chimeric mutant lines had male strobili with or without pollen grains. Our results suggest that CjACOS5 is essential for the production of pollen in sugi and that its disruption is useful for the generation of pollen-free sugi. In addition to conventional transgenic technology, genome editing technology, including CRISPR/Cas9, can confer new traits on sugi.

Construction of a reference transcriptome for the analysis of male sterility in sugi (Cryptomeria japonica D. Don) focusing on MALE STERILITY 1 (MS1)

Sugi (Cryptomeria japonica D. Don) is an important conifer used for afforestation in Japan. As the genome of this species is 11 Gbps, it is too large to assemble within a short timeframe. Transcriptomics is one approach that can address this deficiency. Here we designed a workflow consisting of three stages to de novo assemble transcriptome using Oases and Trinity. The three transcriptomic stage used were independent assembly, automatic and semi-manual integration, and refinement by filtering out potential contamination. We identified a set of 49,795 cDNA and an equal number of translated proteins. According to the benchmark set by BUSCO, 87.01% of cDNAs identified were complete genes, and 78.47% were complete and single-copy genes. Compared to other full-length cDNA resources collected by Sanger and PacBio sequencers, the extent of the coverage in our dataset was the highest, indicating that these data can be safely used for further studies. When two tissue-specific libraries were compared, there were significant expression differences between male strobili and leaf and bark sets. Moreover, subtle expression difference between male-fertile and sterile libraries were detected. Orthologous genes from other model plants and conifer species were identified. We demonstrated that our transcriptome assembly output (CJ3006NRE) can serve as a reference transcriptome for future functional genomics and evolutionary biology studies.

The Transcriptome of Cunninghamia lanceolata male/female cone reveal the association between MIKC MADS-box genes and reproductive organs development

BACKGROUND

Cunninghamia lanceolata (Chinese fir), a member of the conifer family Cupressaceae, is one of the most popular cultivated trees for wood production in China. Continuous research is being performed to improve C. lanceolata breeding values. Given the high rate of seed abortion (one of the reasons being the failure of ovule and pollen development) in C. lanceolata, the proper formation of female/male cones could theoretically increase the number of offspring in future generations. MIKC MADS-box genes are well-known for their roles in the flower/cone development and comprise the typical/atypical floral development model for both angiosperms and gymnosperms.

RESULTS

We performed a transcriptomic analysis to find genes differentially expressed between female and male cones at a single, carefully determined developmental stage, focusing on the MIKC MADS-box genes. We finally obtained 47 unique MIKC MADS-box genes from C. lanceolata and divided these genes into separate branches. 27 out of the 47 MIKC MADS-box genes showed differential expression between female and male cones, and most of them were not expressed in leaves. Out of these 27 genes, most B-class genes (AP3/PI) were up-regulated in the male cone, while TM8 genes were up-regulated in the female cone. Then, with no obvious overall preference for AG (class C + D) genes in female/male cones, it seems likely that these genes are involved in the development of both cones. Finally, a small number of genes such as GGM7, SVP, AGL15, that were specifically expressed in female/male cones, making them candidate genes for sex-specific cone development.

CONCLUSIONS

Our study identified a number of MIKC MADS-box genes showing differential expression between female and male cones in C. lanceolata, illustrating a potential link of these genes with C. lanceolata cone development. On the basis of this, we postulated a possible cone development model for C. lanceolata. The gene expression library showing differential expression between female and male cones shown here, can be used to discover unknown regulatory networks related to sex-specific cone development in the future.

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.

Scanning RNA-Seq and RAD-Seq approach to develop SNP markers closely linked to MALE STERILITY 1 (MS1) in Cryptomeria japonica D. Don

Cryptomeria japonica is a major forestry tree species in Japan. Male sterility of the species is caused by a recessive gene, which shows dysfunction of pollen development and results in no dispersed pollen. Because the pollen of C. japonica induces pollinosis, breeding of pollen-free C. japonica is desired. In this study, single nucleotide polymorphism (SNP) markers located at 1.78 and 0.58 cM to a male sterility locus (MS1) were identified from an analysis of RNA-Seq and RAD-Seq, respectively. SNPs closely linked to MS1 were first scanned by a method similar to MutMap, where a type of index was calculated to measure the strength of the linkage between a marker sequence and MS1. Linkage analysis of selected SNP markers confirmed a higher efficiency of the current method to construct a partial map around MS1. Allele-specific PCR primer pair for the most closely linked SNP with MS1 was developed as a codominant marker, and visualization of the PCR products on an agarose gel enabled rapid screening of male sterile C. japonica. The allele-specific primers developed in this study would be useful for establishing the selection of male sterile C. japonica.

Fine mapping of the male-sterile genes (MS1, MS2, MS3, and MS4) and development of SNP markers for marker-assisted selection in Japanese cedar (Cryptomeria japonica D. Don)

Pollinosis caused by Japanese cedar (Cryptomeria japonica) is a widespread social problem in Japan. To date, 23 male-sterile C. japonica trees have been selected in Japan to address pollinosis, from which four male-sterility loci (MS1, MS2, MS3, and MS4) have been identified from test crossing results. For efficient breeding of male-sterile C. japonica trees, more male-sterile individuals and individuals heterozygous for male-sterile genes are required. Therefore, we aimed to develop DNA markers for marker-assisted selection of four types of male-sterile genes from populations without a family structure. First, for four families exhibiting segregation of each male-sterile locus (MS1, MS2, MS3, and MS4), genome-wide single-nucleotide polymorphism and insertion/deletion (indel) genotyping was performed using the Axiom myDesign Targeted Genotyping Array method. Four high-density linkage maps for mapping the MS1, MS2, MS3, and MS4 families were constructed, which included 4923, 1722, 1896, and 2247 markers, respectively. In these maps, 15, 4, 2, and 2 markers were located 0.0, 3.3, 1.1, and 0.0 cM from the MS1, MS2, MS3, and MS4 loci, respectively. Second, for the markers located 0.0 cM from a male-sterile locus (i.e., MS1 and MS4), to clarify the most tightly linked markers, we calculated the prediction rate of male-sterile gene genotypes from marker genotypes for 78 trees. The markers with the highest prediction rates were AX-174127446 (0.95) for MS1 and AX-174121522 (1.00) for MS4. The AX-174121522 marker was considered to be suitable for selecting trees homozygous or heterozygous for the MS4 gene from plus-trees without a pollination test, which requires a large amount of time and effort. The nearest markers to the male-sterile loci found in this study may facilitate the isolation of male-sterile genes in C. japonica via combination with the draft genomic sequence that is currently being collated.

Identification of novel putative causative genes and genetic marker for male sterility in Japanese cedar (Cryptomeria japonica D.Don)

Background

Japanese cedar (Cryptomeria japonica) is an important tree for Japanese forestry. Male-sterile marker development in Japanese cedar would facilitate selection of male-sterile plus trees, addressing the widespread social problem of pollinosis and facilitating the identification of heterozygotes, which are useful for breeding.

Results

This study used next-generation sequencing for single-nucleotide polymorphism discovery in libraries constructed from several organs, including male-sterile and male-fertile strobili. The single-nucleotide polymorphisms obtained were used to construct a high-density linkage map, which enabled identification of a locus on linkage group 9 strongly correlated with male-sterile trait. Expressed sequence tags corresponding to 11 marker loci from 5 isotigs were associated with this locus within 33.4-34.5 cM. These marker loci explained 100% of the phenotypic variation. Several homologs of these sequences are associated with male sterility in rice or Arabidopsis, including a pre-mRNA splicing factor, a DEAD-box protein, a glycosyl hydrolase, and a galactosyltransferase. These proteins are thus candidates for the causal male-sterile gene at the ms-1 locus. After we used a SNaPshot assay to develop markers for marker-assisted selection (MAS), we tested F₂ progeny between male-sterile and wild-type plus trees to validate the markers and extrapolated the testing to a larger plus-tree population. We found that two developed from one of the candidates for the causal gene were suitable for MAS.

Conclusions

More than half of the ESTs and SNPs we collected were new, enlarging the genomic basis for genetic research on Japanese cedar. We developed two SNP markers aimed at MAS that distinguished individuals carrying the male-sterile trait with 100% accuracy, as well as individuals heterozygous at the male-sterile locus, even outside the mapping population. These markers should enable practical MAS for conifer breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4581-5) contains supplementary material, which is available to authorized users.

Comprehensive analysis of small RNAs expressed in developing male strobili of Cryptomeria japonica

Deep sequencing of small RNAs (sRNAs) in developing male strobili of second-generation offspring originating from a nuclear genic male sterile tree of Cryptomeria japonica were performed to characterize sRNA populations in the male strobili at early pollen developmental stages. Comparing to sequences of microRNA (miRNA) families of plant species and sRNAs expressed in the reproductive organs of representative vascular plants, 37 conserved miRNA families were detected, of which eight were ubiquitously expressed in the reproductive organs of land plant species. In contrast, miR1083 was common in male reproductive organs of gymnosperm species but absent in angiosperm species. In addition to conserved miRNAs, 199 novel miRNAs candidates were predicted. The expression patterns of the obtained sRNAs were further investigated to detect the differentially expressed (DE) sRNAs between genic male sterile and fertile individuals. A total of 969 DE sRNAs were obtained and only three known miRNA families were included among them. These results suggest that both conserved and species-specific sRNAs contribute to the development of male strobili in C. japonica.

Phytoplasma-conserved phyllogen proteins induce phyllody across the Plantae by degrading floral MADS domain proteins

ABCE-class MADS domain transcription factors (MTFs) are key regulators of floral organ development in angiosperms. Aberrant expression of these genes can result in abnormal floral traits such as phyllody. Phyllogen is a virulence factor conserved in phytoplasmas, plant pathogenic bacteria of the class Mollicutes. It triggers phyllody in Arabidopsis thaliana by inducing degradation of A- and E-class MTFs. However, it is still unknown whether phyllogen can induce phyllody in plants other than A. thaliana, although phytoplasma-associated phyllody symptoms are observed in a broad range of angiosperms. In this study, phyllogen was shown to cause phyllody phenotypes in several eudicot species belonging to three different families. Moreover, phyllogen can interact with MTFs of not only angiosperm species including eudicots and monocots but also gymnosperms and a fern, and induce their degradation. These results suggest that phyllogen induces phyllody in angiosperms and inhibits MTF function in diverse plant species.

Determination of male strobilus developmental stages by cytological and gene expression analyses in Japanese cedar (Cryptomeria japonica)

The molecular mechanisms that control male strobilus development in conifers are largely unknown because the developmental stages and related genes have not yet been characterized. The determination of male strobilus developmental stages will contribute to genetic research and reproductive biology in conifers. Our objectives in this study were to determine the developmental stages of male strobili by cytological and transcriptome analysis, and to determine the stages at which aberrant morphology is observed in a male-sterile mutant of Cryptomeria japonica D. Don to better understand the molecular mechanisms that control male strobilus and pollen development. Male strobilus development was observed for 8 months, from initiation to pollen dispersal. A set of 19,209 expressed sequence tags (ESTs) collected from a male reproductive library and a pollen library was used for microarray analysis. We divided male strobilus development into 10 stages by cytological and transcriptome analysis. Eight clusters (7324 ESTs) exhibited major changes in transcriptome profiles during male strobili and pollen development in C. japonica Two clusters showed a gradual increase and decline in transcript abundance, respectively, while the other six clusters exhibited stage-specific changes. The stages at which the male sterility trait of Sosyun was expressed were identified using information on male strobilus and pollen developmental stages and gene expression profiles. Aberrant morphology was observed cytologically at Stage 6 (microspore stage), and differences in expression patterns compared with wild type were observed at Stage 4 (tetrad stage).

Prevalent Exon-Intron Structural Changes in the APETALA1/FRUITFULL, SEPALLATA, AGAMOUS-LIKE6, and FLOWERING LOCUS C MADS-Box Gene Subfamilies Provide New Insights into Their Evolution

AP1/FUL, SEP, AGL6, and FLC subfamily genes play important roles in flower development. The phylogenetic relationships among them, however, have been controversial, which impedes our understanding of the origin and functional divergence of these genes. One possible reason for the controversy may be the problems caused by changes in the exon-intron structure of genes, which, according to recent studies, may generate non-homologous sites and hamper the homology-based sequence alignment. In this study, we first performed exon-by-exon alignments of these and three outgroup subfamilies (SOC1, AG, and STK). Phylogenetic trees reconstructed based on these matrices show improved resolution and better congruence with species phylogeny. In the context of these phylogenies, we traced evolutionary changes of exon-intron structures in each subfamily. We found that structural changes have occurred frequently following gene duplication and speciation events. Notably, exons 7 and 8 (if present) suffered more structural changes than others. With the knowledge of exon-intron structural changes, we generated more reasonable alignments containing all the focal subfamilies. The resulting trees showed that the SEP subfamily is sister to the monophyletic group formed by AP1/FUL and FLC subfamily genes and that the AGL6 subfamily forms a sister group to the three abovementioned subfamilies. Based on this topology, we inferred the evolutionary history of exon-intron structural changes among different subfamilies. Particularly, we found that the eighth exon originated before the divergence of AP1/FUL, FLC, SEP, and AGL6 subfamilies and degenerated in the ancestral FLC-like gene. These results provide new insights into the origin and evolution of the AP1/FUL, FLC, SEP, and AGL6 subfamilies.

Evolutionary rate variation in two conifer species, Taxodium distichum (L.) Rich. var. distichum (baldcypress) and Cryptomeria japonica (Thunb. ex L.f.) D. Don (Sugi, Japanese cedar)

With the advance of sequencing technologies, large-scale data of expressed sequence tags and full-length cDNA sequences have been reported for several conifer species. Comparative analyses of evolutionary rates among diverse taxa provide insights into taxon-specific molecular evolutionary features and into the origin of variation in evolutionary rates within genomes and between species. Here, we estimated evolutionary rates in two conifer species, Taxodium distichum and Cryptomeria japonica, to illuminate the molecular evolutionary features of these species, using hundreds of genes and employing Chamaecyparis obtusa as an outgroup. Our results show that the mutation rates based on synonymous substitution rates (dS) of T. distichum and C. japonica are approximately 0.67 × 10(-9) and 0.59 × 10(-9)/site/year, respectively, which are 15-25 times lower than those of annual angiosperms. We found a significant positive correlation between dS and GC3. This implies that a local mutation bias, such as context dependency of the mutation bias, exists within the genomes of T. distichum and C. japonica, and/or that selection acts on synonymous sites in these species. In addition, the means of the ratios of synonymous to nonsynonymous substitution rate in the two species are almost the same, suggesting that the average intensity of functional constraint is constant between the lineages. Finally, we tested the possibility of positive selection based on the site model, and detected one candidate gene for positive selection.

MADS goes genomic in conifers: towards determining the ancestral set of MADS-box genes in seed plants

BACKGROUND AND AIMS

MADS-box genes comprise a gene family coding for transcription factors. This gene family expanded greatly during land plant evolution such that the number of MADS-box genes ranges from one or two in green algae to around 100 in angiosperms. Given the crucial functions of MADS-box genes for nearly all aspects of plant development, the expansion of this gene family probably contributed to the increasing complexity of plants. However, the expansion of MADS-box genes during one important step of land plant evolution, namely the origin of seed plants, remains poorly understood due to the previous lack of whole-genome data for gymnosperms.

METHODS

The newly available genome sequences of Picea abies, Picea glauca and Pinus taeda were used to identify the complete set of MADS-box genes in these conifers. In addition, MADS-box genes were identified in the growing number of transcriptomes available for gymnosperms. With these datasets, phylogenies were constructed to determine the ancestral set of MADS-box genes of seed plants and to infer the ancestral functions of these genes.

KEY RESULTS

Type I MADS-box genes are under-represented in gymnosperms and only a minimum of two Type I MADS-box genes have been present in the most recent common ancestor (MRCA) of seed plants. In contrast, a large number of Type II MADS-box genes were found in gymnosperms. The MRCA of extant seed plants probably possessed at least 11-14 Type II MADS-box genes. In gymnosperms two duplications of Type II MADS-box genes were found, such that the MRCA of extant gymnosperms had at least 14-16 Type II MADS-box genes.

CONCLUSIONS

The implied ancestral set of MADS-box genes for seed plants shows simplicity for Type I MADS-box genes and remarkable complexity for Type II MADS-box genes in terms of phylogeny and putative functions. The analysis of transcriptome data reveals that gymnosperm MADS-box genes are expressed in a great variety of tissues, indicating diverse roles of MADS-box genes for the development of gymnosperms. This study is the first that provides a comprehensive overview of MADS-box genes in conifers and thus will provide a framework for future work on MADS-box genes in seed plants.

De novo assembly of maritime pine transcriptome: implications for forest breeding and biotechnology

Maritime pine (Pinus pinasterAit.) 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.

Transcriptome sequencing and profiling of expressed genes in cambial zone and differentiating xylem of Japanese cedar (Cryptomeria japonica)

BACKGROUND

Forest trees have ecological and economic importance, and Japanese cedar has highly valued wood attributes. Thus, studies of molecular aspects of wood formation offer practical information that may be used for screening and forward genetics approaches to improving wood quality.

RESULTS

After identifying expressed sequence tags in Japanese cedar tissue undergoing xylogenesis, we designed a custom cDNA microarray to compare expression of highly regulated genes throughout a growing season. This led to identification of candidate genes involved both in wood formation and later cessation of growth and dormancy. Based on homology to orthologous protein groups, the genes were assigned to functional classes. A high proportion of sequences fell into functional classes related to posttranscriptional modification and signal transduction, while transcription factors and genes involved in the metabolism of sugars, cell-wall synthesis and lignification, and cold hardiness were among other classes of genes identified as having a potential role in xylem formation and seasonal wood formation.

CONCLUSIONS

We obtained 55,051 unique sequences by next-generation sequencing of a cDNA library prepared from cambial meristem and derivative cells. Previous studies on conifers have identified unique sequences expressed in developing xylem, but this is the first comprehensive study utilizing a collection of expressed sequence tags for expression studies related to xylem formation in Japanese cedar, which belongs to a different lineage than the Pinaceae. Our characterization of these sequences should allow comparative studies of genome evolution and functional genetics of wood species.

A novel approach of preventing Japanese cedar pollen dispersal that is the cause of Japanese cedar pollinosis (JCP) using pollen-specific fungal infection

In Japan, Japanese cedar pollen dispersal is one of the major causes of pollinosis. Sydowia japonica is an ascomycetous fungus that grows exclusively on the male strobili of Japanese cedar, suggesting a possible mechanism for controlling pollen dispersal. To evaluate this possibility, eleven isolates of S. japonica were collected from around Japan and used as an inoculum to male strobili of Japanese cedar. The treatment demonstrated that the fungus infected only the pollen and prevented pollen dispersal. The fungus did not cause any additional symptoms to other parts of Japanese cedar, such as needles, stems, and buds. All S. japonica isolates collected around Japan could serve to control pollen dispersal. Periodic observation of the fungal pathogenesis with stereomicroscope and scanning electron microscope showed that hyphal fragments and conidia of S. japonica germinated on the surface of male strobili, and the germ tube entered pollen sacs through opening microsporophylls. Within the pollen sacs, the hyphae penetrated pollen gradually, such that all pollen was infected by the fungus by approximately one month before the pollen dispersal season. The infected pollen was destroyed due to the fungal infection and was never released. Our data suggests a novel approach of preventing pollen dispersal using pollen-specific fungal infection.

An EST dataset for Metasequoia glyptostroboides buds: the first EST resource for molecular genomics studies in Metasequoia

The "living fossil" Metasequoia glyptostroboides Hu et Cheng, commonly known as dawn redwood or Chinese redwood, is the only living species in the genus and is valued for its essential oil and crude extracts that have great potential for anti-fungal activity. Despite its paleontological significance and economical value as a rare relict species, genomic resources of Metasequoia are very limited. In order to gain insight into the molecular mechanisms behind the formation of reproductive buds and the transition from vegetative phase to reproductive phase in Metasequoia, we performed sequencing of expressed sequence tags from Metasequoia vegetative buds and female buds. By using the 454 pyrosequencing technology, a total of 1,571,764 high-quality reads were generated, among which 733,128 were from vegetative buds and 775,636 were from female buds. These EST reads were clustered and assembled into 114,124 putative unique transcripts (PUTs) with an average length of 536 bp. The 97,565 PUTs that were at least 100 bp in length were functionally annotated by a similarity search against public databases and assigned with Gene Ontology (GO) terms. A total of 59 known floral gene families and 190 isotigs involved in hormone regulation were captured in the dataset. Furthermore, a set of PUTs differentially expressed in vegetative and reproductive buds, as well as SSR motifs and high confidence SNPs, were identified. This is the first large-scale expressed sequence tags ever generated in Metasequoia and the first evidence for floral genes in this critically endangered deciduous conifer species.

Towards decoding the conifer giga-genome

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world's conifer forests.

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.

A white spruce gene catalog for conifer genome analyses

Several angiosperm plant genomes, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), poplar (Populus trichocarpa), and grapevine (Vitis vinifera), have been sequenced, but the lack of reference genomes in gymnosperm phyla reduces our understanding of plant evolution and restricts the potential impacts of genomics research. A gene catalog was developed for the conifer tree Picea glauca (white spruce) through large-scale expressed sequence tag sequencing and full-length cDNA sequencing to facilitate genome characterizations, comparative genomics, and gene mapping. The resource incorporates new and publicly available sequences into 27,720 cDNA clusters, 23,589 of which are represented by full-length insert cDNAs. Expressed sequence tags, mate-pair cDNA clone analysis, and custom sequencing were integrated through an iterative process to improve the accuracy of clustering outcomes. The entire catalog spans 30 Mb of unique transcribed sequence. We estimated that the P. glauca nuclear genome contains up to 32,520 transcribed genes owing to incomplete, partially sequenced, and unsampled transcripts and that its transcriptome could span up to 47 Mb. These estimates are in the same range as the Arabidopsis and rice transcriptomes. Next-generation methods confirmed and enhanced the catalog by providing deeper coverage for rare transcripts, by extending many incomplete clusters, and by augmenting the overall transcriptome coverage to 38 Mb of unique sequence. Genomic sample sequencing at 8.5% of the 19.8-Gb P. glauca genome identified 1,495 clusters representing highly repeated sequences among the cDNA clusters. With a conifer transcriptome in full view, functional and protein domain annotations clearly highlighted the divergences between conifers and angiosperms, likely reflecting their respective evolutionary paths.

Gene discovery from Jatropha curcas by sequencing of ESTs from normalized and full-length enriched cDNA library from developing seeds

BACKGROUND

Jatropha curcas L. is promoted as an important non-edible biodiesel crop worldwide. Jatropha oil, which is a triacylglycerol, can be directly blended with petro-diesel or transesterified with methanol and used as biodiesel. Genetic improvement in jatropha is needed to increase the seed yield, oil content, drought and pest resistance, and to modify oil composition so that it becomes a technically and economically preferred source for biodiesel production. However, genetic improvement efforts in jatropha could not take advantage of genetic engineering methods due to lack of cloned genes from this species. To overcome this hurdle, the current gene discovery project was initiated with an objective of isolating as many functional genes as possible from J. curcas by large scale sequencing of expressed sequence tags (ESTs).

RESULTS

A normalized and full-length enriched cDNA library was constructed from developing seeds of J. curcas. The cDNA library contained about 1 × 10(6) clones and average insert size of the clones was 2.1 kb. Totally 12,084 ESTs were sequenced to average high quality read length of 576 bp. Contig analysis revealed 2258 contigs and 4751 singletons. Contig size ranged from 2-23 and there were 7333 ESTs in the contigs. This resulted in 7009 unigenes which were annotated by BLASTX. It showed 3982 unigenes with significant similarity to known genes and 2836 unigenes with significant similarity to genes of unknown, hypothetical and putative proteins. The remaining 191 unigenes which did not show similarity with any genes in the public database may encode for unique genes. Functional classification revealed unigenes related to broad range of cellular, molecular and biological functions. Among the 7009 unigenes, 6233 unigenes were identified to be potential full-length genes.

CONCLUSIONS

The high quality normalized cDNA library was constructed from developing seeds of J. curcas for the first time and 7009 unigenes coding for diverse biological functions including oil biosynthesis were identified. These genes will serve as invaluable genetic resource for crop improvement in jatropha to make it an ideal and profitable crop for biodiesel production.

Molecular interactions of orthologues of floral homeotic proteins from the gymnosperm Gnetum gnemon provide a clue to the evolutionary origin of 'floral quartets'

Several lines of evidence suggest that the identity of floral organs in angiosperms is specified by multimeric transcription factor complexes composed of MADS-domain proteins. These bind to specific cis-regulatory elements ('CArG-boxes') of their target genes involving DNA-loop formation, thus constituting 'floral quartets'. Gymnosperms, angiosperms' closest relatives, contain orthologues of floral homeotic genes, but when and how the interactions constituting floral quartets were established during evolution has remained unknown. We have comprehensively studied the dimerization and DNA-binding of several classes of MADS-domain proteins from the gymnosperm Gnetum gnemon. Determination of protein-protein and protein-DNA interactions by yeast two-hybrid, in vitro pull-down and electrophoretic mobility shift assays revealed complex patterns of homo- and heterodimerization among orthologues of floral homeotic class B, class C and class E proteins and B(sister) proteins. Using DNase I footprint assays we demonstrate that both orthologues of class B with C proteins, and orthologues of class C proteins alone, but not orthologues of class B proteins alone can loop DNA in floral quartet-like complexes. This is in contrast to class B and class C proteins from angiosperms, which require other factors such as class E floral homeotic proteins to 'glue' them together in multimeric complexes. Our findings suggest that the evolutionary origin of floral quartet formation is based on the interaction of different DNA-bound homodimers, does not depend on class E proteins, and predates the origin of angiosperms.

Genomics and bioinformatics resources for crop improvement

Recent remarkable innovations in platforms for omics-based research and application development provide crucial resources to promote research in model and applied plant species. A combinatorial approach using multiple omics platforms and integration of their outcomes is now an effective strategy for clarifying molecular systems integral to improving plant productivity. Furthermore, promotion of comparative genomics among model and applied plants allows us to grasp the biological properties of each species and to accelerate gene discovery and functional analyses of genes. Bioinformatics platforms and their associated databases are also essential for the effective design of approaches making the best use of genomic resources, including resource integration. We review recent advances in research platforms and resources in plant omics together with related databases and advances in technology.

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 naked and the dead: the ABCs of gymnosperm reproduction and the origin of the angiosperm flower

20 years after establishment of the ABC model many of the molecular mechanisms underlying development of the angiosperm flower are relatively well understood. Central players in the gene regulatory network controlling flower development are SQUA-like, DEF/GLO-like, AG-like and AGL6/SEP1-like MIKC-type MADS-domain transcription factors. These provide class A, class B, class C and the more recently defined class E floral homeotic functions, respectively. There is evidence that the floral homeotic proteins recognize the DNA of target genes in an organ-specific way as multimeric protein complexes, thus constituting 'floral quartets'. In contrast to the detailed insights into flower development, how the flower originated during evolution has remained enigmatic. However, while orthologues of all classes of floral homeotic genes appear to be absent from all non-seed plants, DEF/GLO-like, AG-like, and AGL6-like genes have been found in diverse extant gymnosperms, the closest relatives of the angiosperms. While SQUA-like and SEP1-like MADS-box genes appear to be absent from extant gymnosperms, reconstruction of MADS-box gene phylogeny surprisingly suggests that the most recent common ancestor of gymnosperms and angiosperms possessed representatives of both genes, but that these have been lost in the lineage that led to extant gymnosperms. Expression studies and genetic complementation experiments indicate that both angiosperm and gymnosperm AG-like and DEF/GLO-like genes have conserved functions in the specification of reproductive organs and in distinguishing male from female organs, respectively. Based on these findings novel models about the molecular basis of flower origin, involving changes in the expression patterns of DEF/GLO-like or AGL6/SEP1/SQUA-like genes in reproductive structures, were developed. While in angiosperms SEP1-like proteins play an important role in floral quartet formation, preliminary evidence suggests that gymnosperm DEF/GLO-like and AG-like proteins alone can already form floral quartet-like complexes, further corroborating the view that the formation of floral quartet-like complexes predated flower origin during evolution.

Functional genomics using RIKEN Arabidopsis thaliana full-length cDNAs

Full-length cDNAs are essential for the correct annotation of genomic sequences as well as for the functional analysis of genes and their products. We have isolated about 240,000 RIKEN Arabidopsis full-length (RAFL) cDNA clones. These clones were clustered into about 17,000 non-redundant cDNA groups, i.e., about 60% of all Arabidopsis predicted genes. The sequence information of the RAFL cDNAs is useful for promoter analysis, and for the correct annotation of predicted transcriptional units and gene products. We prepared cDNA microarrays containing independent full-length cDNA groups and studied the expression profiles of genes under various stress- and hormone-treatment conditions, and in various mutants and transgenic plants. These expression profiling studies have shown the expression levels of many genes as a detailed snapshot describing the state of a biological system in planta under various conditions. We have applied RAFL cDNAs to the functional analysis of proteins using the full-length cDNA over-expressing (FOX) gene hunting system and the wheat germ cell-free protein synthesis system. The RAFL cDNA collection was also used for determination of the domain structure of proteins by NMR. In this review, we summarize the present state and perspectives of functional genomics using RAFL cDNAs.

Phenome analysis in plant species using loss-of-function and gain-of-function mutants

Analysis of genetic mutations is one of the most effective ways to investigate gene function. We now have methods that allow for mass production of mutant lines and cells in a variety of model species. Recently, large numbers of mutant lines have been generated by both 'loss-of-function' and 'gain-of-function' techniques. In parallel, phenotypic information covering various mutant resources has been acquired and released in web-based databases. As a result, significant progress in comprehensive phenotype analysis is being made through the use of these tools. Arabidopsis and rice are two major model plant species in which genome sequencing projects have been completed. Arabidopsis is the most widely used experimental plant, with a large number of mutant resources and several examples of systematic phenotype analysis. Rice is a major crop species and is used as a model plant, with an increasing number of mutant resources. Other plant species are also being employed in functional genetics research. In this review, the present status of mutant resources for large-scale studies of gene function in plant research and the current perspective on using loss-of-function and gain-of-function mutants in phenome research will be discussed.

Genomics of growth traits in forest trees

Growth traits in trees are fundamental components of adaptation in a forest ecosystem and of productivity in planted forests. A number of processes determine tree growth, which are controlled by genetic and epigenetic factors that respond dynamically to environmental signals throughout centuries. Advances in genomics have allowed an increased comprehension of the complex mechanisms of tree growth and adaptation. Yet, the application of genomics to improving forest productivity and sustainability still entails capturing a large proportion of the total genetic variation controlling the component traits. Nonetheless, genetics and genomics are unifying disciplines that will serve well to dissect the variables and mechanisms of tree growth and development.