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

A single-nucleotide substitution of CjTKPR1 determines pollen production in the gymnosperm plant Cryptomeria japonica

Pollinosis, also known as pollen allergy or hay fever, is a global problem caused by pollen produced by various plant species. The wind-pollinated Japanese cedar (Cryptomeria japonica) is the largest contributor to severe pollinosis in Japan, where increasing proportions of people have been affected in recent decades. The MALE STERILITY 4 (MS4) locus of Japanese cedar controls pollen production, and its homozygous mutants (ms4/ms4) show abnormal pollen development after the tetrad stage and produce no mature pollen. In this study, we narrowed down the MS4 locus by fine mapping in Japanese cedar and found TETRAKETIDE α-PYRONE REDUCTASE 1 (TKPR1) gene in this region. Transformation experiments using Arabidopsis thaliana showed that single-nucleotide substitution ("T" to "C" at 244-nt position) of CjTKPR1 determines pollen production. Broad conservation of TKPR1 beyond plant division could lead to the creation of pollen-free plants not only for Japanese cedar but also for broader plant species.

CRISPR/Cas9-mediated targeted mutagenesis in Japanese cedar (Cryptomeria japonica D. Don)

Cryptomeria japonica (Japanese cedar or sugi) is one of the most important coniferous tree species in Japan and breeding programs for this species have been launched since 1950s. Genome editing technology can be used to shorten the breeding period. In this study, we performed targeted mutagenesis using the CRISPR/Cas9 system in C. japonica. First, the CRISPR/Cas9 system was tested using green fluorescent protein (GFP)-expressing transgenic embryogenic tissue lines. Knock-out efficiency of GFP ranged from 3.1 to 41.4% depending on U6 promoters and target sequences. The GFP knock-out region was mottled in many lines, indicating genome editing in individual cells. However, in 101 of 102 mutated individuals (> 99%) from 6 GFP knock-out lines, embryos had a single mutation pattern. Next, we knocked out the endogenous C. japonica magnesium chelatase subunit I (CjChlI) gene using two guide RNA targets. Green, pale green, and albino phenotypes were obtained in the gene-edited cell lines. Sequence analysis revealed random deletions, insertions, and replacements in the target region. Thus, targeted mutagenesis using the CRISPR/Cas9 system can be used to modify the C. japonica genome.

Selection and Optimization of Reference Genes for MicroRNA Expression Normalization by qRT-PCR in Chinese Cedar (Cryptomeria fortunei) under Multiple Stresses

MicroRNA (miRNA) expression analysis is very important for investigating its functions. To date, no research on reference genes (RGs) for miRNAs in gymnosperms, including Cryptomeria fortunei, has been reported. Here, ten miRNAs (i.e., pab-miR159a, cln-miR162, cas-miR166d, pab-miR395b, ppt-miR894, cln-miR6725, novel1, novel6, novel14 and novel16) and three common RGs (U6, 5S and 18S) were selected as candidate RGs. qRT-PCR was used to analyse their expressions in C. fortunei under various experimental conditions, including multiple stresses (cold, heat, drought, salt, abscisic acid and gibberellin) and in various tissues (roots, stems, tender needles, cones and seeds). Four algorithms (delta Ct, geNorm, NormFinder and BestKeeper) were employed to assess the stability of candidate RG expression; the geometric mean and RefFinder program were used to comprehensively evaluate RG stability. According to the results, novel16, cln-miR6725, novel1 and U6 were the most stable RGs for studying C. fortunei miRNA expression. In addition, the expression of three target miRNAs (aly-miR164c-5p, aly-miR168a-5p and smo-miR396) was examined to verify that the selected RGs are suitable for miRNA expression normalisation. This study may aid further investigations of miRNA expression/function in the response of C. fortunei to abiotic stress and provides an important basis for the standardisation of miRNA expression in other gymnosperm species.

Factors Affecting the Number of Pollen Grains per Male Strobilus in Japanese Cedar (Cryptomeria japonica)

Japanese cedar (Cryptomeria japonica) is the most important timber species in Japan; however, its pollen is the primary cause of pollinosis in Japan. The total number of pollen grains produced by a single tree is determined by the number of male strobili (male flowers) and the number of pollen grains per male strobilus. While the number of male strobili is a visible and well-investigated trait, little is known about the number of pollen grains per male strobilus. We hypothesized that genetic and environmental factors affect the pollen number per male strobilus and explored the factors that affect pollen production and genetic variation among clones. We counted pollen numbers of 523 male strobili from 26 clones using a cell counter method that we recently developed. Piecewise Structural Equation Modeling (pSEM) revealed that the pollen number is mostly affected by genetic variation, male strobilus weight, and pollen size. Although we collected samples from locations with different environmental conditions, statistical modeling succeeded in predicting pollen numbers for different clones sampled from branches facing different directions. Comparison of predicted pollen numbers revealed that they varied >3-fold among the 26 clones. The determination of the factors affecting pollen number and a precise evaluation of genetic variation will contribute to breeding strategies to counter pollinosis. Furthermore, the combination of our efficient counting method and statistical modeling will provide a powerful tool not only for Japanese cedar but also for other plant species.

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.

An improved pollen number counting method using a cell counter and mesh columns

BACKGROUND

The determination of pollen number is important in evolutionary, agricultural, and medical studies. Tree species of the Cupressaceae family cause serious pollinosis worldwide. Although Japanese cedar (Cryptomeria japonica) is the most important forestry species in Japan, it is also the biggest cause of pollinosis in the country. Japanese cedar trees have been selected for growth speed and superior morphological traits and then cloned. These clones may vary in their pollen production, but there has been little research on how many pollen grains are produced by a single male strobilus (flower). A recently reported method for counting pollen number with a cell counter was applicable to Arabidopsis species and wheat, but was not suitable for Japanese cedar because the strobilus does not open with heating (e.g. 60 °C, overnight).

RESULTS

Here, we report an improved pollen counting method for Japanese cedar using a precise and rapid cell counter in combination with home-made mesh columns. The male strobilus was gently crushed using a pestle. Large and small debris were then removed using 100- and 20-μm mesh columns, respectively. We successfully detected pollen sizes and numbers that differed between two clones using this method.

CONCLUSIONS

This improved method is not only suitable for counting pollen from Japanese cedar, but could also be applied to other species of the Cupressaceae family with hard scale tissue covering the pollen. Moreover, this method could be applied to a broader range of plant species, such as wheat, because there is no need to wait for anthesis and debris can be removed efficiently.

Transcriptome Analysis in Male Strobilus Induction by Gibberellin Treatment in Cryptomeria japonica D. Don

The plant hormone gibberellin (GA) is known to regulate elongating growth, seed germination, and the initiation of flower bud formation, and it has been postulated that GAs originally had functions in reproductive processes. Studies on the mechanism of induction of flowering by GA have been performed in Arabidopsis and other model plants. In coniferous trees, reproductive organ induction by GAs is known to occur, but there are few reports on the molecular mechanism in this system. To clarify the gene expression dynamics of the GA induction of the male strobilus in Cryptomeria japonica, we performed comprehensive gene expression analysis using a microarray. A GA-treated group and a nontreated group were allowed to set, and individual trees were sampled over a 6-week time course. A total of 881 genes exhibiting changed expression was identified. In the GA-treated group, genes related to ‘stress response’ and to ‘cell wall’ were initially enriched, and genes related to ‘transcription’ and ‘transcription factor activity’ were enriched at later stages. This analysis also clarified the dynamics of the expression of genes related to GA signaling transduction following GA treatment, permitting us to compare and contrast with the expression dynamics of genes implicated in signal transduction responses to other plant hormones. These results suggested that various plant hormones have complex influences on the male strobilus induction. Additionally, principal component analysis (PCA) using expression patterns of the genes that exhibited sequence similarity with flower bud or floral organ formation-related genes of Arabidopsis was performed. PCA suggested that gene expression leading to male strobilus formation in C. japonica became conspicuous within one week of GA treatment. Together, these findings help to clarify the evolution of the mechanism of induction of reproductive organs by GA.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

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.

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.