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

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.

Changes in annual transcriptome dynamics of a clone of Japanese cedar (Cryptomeria japonica D. Don) planted under different climate conditions

Environmental responses are critical for plant growth and survival under different climate conditions. To elucidate the underlying biological mechanisms of environmental responses in Japanese cedar (Cryptomeria japonica D. Don), the annual transcriptome dynamics of common clonal trees (Godai1) planted at three different climate sites (Yamagata, Ibaraki, and Kumamoto Prefectures) were analyzed using microarrays. Both principal component analysis (PCA) and hierarchical clustering of the microarray data indicated the transition to dormant transcriptome status occurred earlier and the transition to active growth status later in the colder region. Interestingly, PCA also indicated that the transcriptomes of trees grown under three different conditions were similar during the growth period (June to September), whereas the transcriptomes differed between sites during the dormant period (January to March). In between-site comparisons, analyses of the annual expression profiles of genes for sites 'Yamagata vs. Kumamoto', 'Yamagata vs. Ibaraki', and 'Ibaraki vs. Kumamoto' identified 1,473, 1,137, and 925 targets exhibiting significantly different expression patterns, respectively. The total of 2,505 targets that exhibited significantly different expression patterns in all three comparisons may play important roles in enabling cuttings to adapt to local environmental conditions. Partial least-squares regression analysis and Pearson correlation coefficient analysis revealed that air temperature and day length were the dominant factors controlling the expression levels of these targets. GO and Pfam enrichment analyses indicated that these targets include genes that may contribute to environmental adaptation, such as genes related to stress and abiotic stimulus responses. This study provided fundamental information regarding transcripts that may play an important role in adaptation to environmental conditions at different planting sites.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica)

BACKGROUND

Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated.

RESULTS

Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species.

CONCLUSIONS

The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

Detection of candidate gene LsACOS5 and development of InDel marker for male sterility by ddRAD-seq and resequencing analysis in lettuce

A new breeding method of F1 hybrid using male sterility would open an exciting frontier in lettuce breeding, a self-pollinating crop. Male sterility is a crucial trait in F1 hybrid breeding. It is essential to map the causative gene for using male sterility. The ms-S, male-sterile (MS) gene of 'CGN17397', was mapped to linkage group (LG) 8 by ddRAD-seq and narrowed down between two markers using two F2 populations. This region spans approximately 10.16 Mb, where 94 genes were annotated according to the lettuce reference genome sequence (version8 from 'Salinas'). The whole-genome sequencing of the MS lines 'CGN17397-MS' and male-fertile (MF) lines 'CGN17397-MF' revealed that only one gene differed in the area of Lsat_1_v5_gn_8_148221.1, a homolog of acyl-CoA synthetase5 (ACOS5), and was deleted in the MS lines. It was reported that ACOS5 was needed for pollen wall formation and that the null mutants of ACOS5 were entirely male sterility in some plants. Thus, I concluded that Lsat_1_v5_gn_8_148221.1 designated as LsACOS5 was a biologically plausible candidate gene for the ms-S locus. By using the structural polymorphism of LsACOS5, an InDel marker was developed to select the MS trait. The results obtained here provide valuable information for the genic male-sterility in lettuce.

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.

Toward genomic selection in Pinus taeda: Integrating resources to support array design in a complex conifer genome

PREMISE

An informatics approach was used for the construction of an Axiom genotyping array from heterogeneous, high-throughput sequence data to assess the complex genome of loblolly pine (Pinus taeda).

METHODS

High-throughput sequence data, sourced from exome capture and whole genome reduced-representation approaches from 2698 trees across five sequence populations, were analyzed with the improved genome assembly and annotation for the loblolly pine. A variant detection, filtering, and probe design pipeline was developed to detect true variants across and within populations. From 8.27 million variants, a total of 642,275 were evaluated and 423,695 of those were screened across a range-wide population.

RESULTS

The final informatics and screening approach delivered an Axiom array representing 46,439 high-confidence variants to the forest tree breeding and genetics community. Based on the annotated reference genome, 34% were located in or directly upstream or downstream of genic regions.

DISCUSSION

The Pita50K array represents a genome-wide resource developed from sequence data for an economically important conifer, loblolly pine. It uniquely integrates independent projects that assessed trees sampled across the native range. The challenges associated with the large and repetitive genome are addressed in the development of this resource.

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.

Identification and genetic diversity analysis of a male-sterile gene (MS1) in Japanese cedar (Cryptomeria japonica D. Don)

Identifying causative genes for a target trait in conifer reproduction is challenging for species lacking whole-genome sequences. In this study, we searched for the male-sterility gene (MS1) in Cryptomeria japonica, aiming to promote marker-assisted selection (MAS) of male-sterile C. japonica to reduce the pollinosis caused by pollen dispersal from artificial C. japonica forests in Japan. We searched for mRNA sequences expressed in male strobili and found the gene CJt020762, coding for a lipid transfer protein containing a 4-bp deletion specific to male-sterile individuals. We also found a 30-bp deletion by sequencing the entire gene of another individual with the ms1. All nine breeding materials with the allele ms1 had either a 4-bp or 30-bp deletion in gene CJt020762, both of which are expected to result in faulty gene transcription and function. Furthermore, the 30-bp deletion was detected from three of five individuals in the Ishinomaki natural forest. From our findings, CJt020762 was considered to be the causative gene of MS1. Thus, by performing MAS using two deletion mutations as a DNA marker, it will be possible to find novel breeding materials of C. japonica with the allele ms1 adapted to the unique environment of each region of the Japanese archipelago.

Genotype and transcriptome effects on somatic embryogenesis in Cryptomeria japonica

Somatic embryogenesis (SE), which is in vitro regeneration of plant bodies from somatic cells, represents a useful means of clonal propagation and genetic engineering of forest trees. While protocols to obtain calluses and induce regeneration in somatic embryos have been reported for many tree species, the knowledge of molecular mechanisms of SE development is still insufficient to achieve an efficient supply of somatic embryos required for the industrial application. Cryptomeria japonica, a conifer species widely used for plantation forestry in Japan, is one of the tree species waiting for a secure SE protocol; the probability of normal embryo development appears to depend on genotype. To discriminate the embryogenic potential of embryonal masses (EMs) and efficiently obtain normal somatic embryos of C. japonica, we investigated the effects of genotype and transcriptome on the variation in embryogenic potential. Using an induction experiment with 12 EMs each from six genotypes, we showed that embryogenic potential differs between/within genotypes. Comparisons of gene expression profiles among EMs with different embryogenic potentials revealed that 742 differently expressed genes were mainly associated with pattern forming and metabolism. Thus, we suggest that not only genotype but also gene expression profiles can determine success in SE development. Consistent with previous findings for other conifer species, genes encoding leafy cotyledon, wuschel, germin-like proteins, and glutathione-S-transferases are likely to be involved in SE development in C. japonica and indeed highly expressed in EMs with high-embryogenic potential; therefore, these proteins represent candidate markers for distinguishing embryogenic potential.

Alternative splicing and its regulatory role in woody plants

Alternative splicing (AS) is an important post-transcriptional process to enhance proteome diversity in eukaryotic organisms. In plants, numerous reports have primarily focused on AS analysis in model plant species or herbaceous plants, leading to a notable lack of research on AS in woody plants. More importantly, emerging evidence indicates that many important traits, including wood formation and stress resistance, in woody plants are controlled by AS. In this review article, we summarize the current progress of all kinds of AS studies in different tree species at various stages of development and in response to various stresses, revealing the significant role played by AS in woody plants, as well as the similar properties and differential regulation within their herbaceous counterparts. Furthermore, we propose several potential strategies to facilitate the functional characterization of splicing factors in woody plants and evaluate a general pipeline for the systematic characterization of splicing isoforms in a complex AS regulatory network. The utilization of genetic studies and high-throughput omics integration approaches to analyze AS genes and splicing factors is likely to further advance our understanding of AS modulation in woody plants.

A 3K Axiom SNP array from a transcriptome-wide SNP resource sheds new light on the genetic diversity and structure of the iconic subtropical conifer tree Araucaria angustifolia (Bert.) Kuntze

High-throughput SNP genotyping has become a precondition to move to higher precision and wider genome coverage genetic analysis of natural and breeding populations of non-model species. We developed a 44,318 annotated SNP catalog for Araucaria angustifolia, a grandiose subtropical conifer tree, one of the only two native Brazilian gymnosperms, critically endangered due to its valuable wood and seeds. Following transcriptome assembly and annotation, SNPs were discovered from RNA-seq and pooled RAD-seq data. From the SNP catalog, an Axiom® SNP array with 3,038 validated SNPs was developed and used to provide a comprehensive look at the genetic diversity and structure of 15 populations across the natural range of the species. RNA-seq was a far superior source of SNPs when compared to RAD-seq in terms of conversion rate to polymorphic markers on the array, likely due to the more efficient complexity reduction of the huge conifer genome. By matching microsatellite and SNP data on the same set of A. angustifolia individuals, we show that SNPs reflect more precisely the actual genome-wide patterns of genetic diversity and structure, challenging previous microsatellite-based assessments. Moreover, SNPs corroborated the known major north-south genetic cline, but allowed a more accurate attribution to regional versus among-population differentiation, indicating the potential to select ancestry-informative markers. The availability of a public, user-friendly 3K SNP array for A. angustifolia and a catalog of 44,318 SNPs predicted to provide ~29,000 informative SNPs across ~20,000 loci across the genome, will allow tackling still unsettled questions on its evolutionary history, toward a more comprehensive picture of the origin, past dynamics and future trend of the species' genetic resources. Additionally, but not less importantly, the SNP array described, unlocks the potential to adopt genomic prediction methods to accelerate the still very timid efforts of systematic tree breeding of A. angustifolia.

SNP Genotyping with Target Amplicon Sequencing Using a Multiplexed Primer Panel and Its Application to Genomic Prediction in Japanese Cedar, Cryptomeria japonica (L.f.) D.Don

Along with progress in sequencing technology and accumulating knowledge of genome and gene sequences, molecular breeding techniques have been developed for predicting the genetic potential of individual genotypes and for selecting superior individuals. For Japanese cedar (Cryptomeria japonica (L.f.) D.Don), which is the most common coniferous species in Japanese forestry, we constructed a custom primer panel for target amplicon sequencing in order to simultaneously determine 3034 informative single nucleotide polymorphisms (SNPs). We performed primary evaluation of the custom primer panel with actual sequencing and in silico PCR. Genotyped SNPs had a distribution over almost the entire region of the C. japonica linkage map and verified the high reproducibility of genotype calls compared to SNPs obtained by genotyping arrays. Genotyping was performed for 576 individuals of the F1 population, and genomic prediction models were constructed for growth and wood property-related traits using the genotypes. Amplicon sequencing with the custom primer panel enables efficient obtaining genotype data in order to perform genomic prediction, manage clones, and advance forest tree breeding.

Somatic Embryogenesis and Plant Regeneration from Sugi (Japanese Cedar, Cryptomeria japonica D. Don, Cupressaceae) Seed Families by Marker Assisted Selection for the Male Sterility Allele ms1

One of the possible countermeasures for pollinosis caused by sugi (Cryptomeria japonica), a serious public health problem in Japan, is the use of male sterile plants (MSPs; pollen-free plants). However, the production efficiencies of MSPs raised by conventional methods are extremely poor, time consuming, and resulting in a high seedling cost. Here, we report the development of a novel technique for efficient production of MSPs, which combines marker-assisted selection (MAS) and somatic embryogenesis (SE). SE from four full sib seed families of sugi, carrying the male sterility gene MS1, was initiated using megagametophyte explants that originated from four seed collections taken at one-week intervals during the month of July 2017. Embryogenic cell lines (ECLs) were achieved in all families, with initiation rates varying from 0.6% to 59%. Somatic embryos were produced from genetic marker-selected male sterile ECLs on medium containing maltose, abscisic acid (ABA), polyethylene glycol (PEG), and activated charcoal (AC). Subsequently, high frequencies of germination and plant conversion (≥76%) were obtained on plant growth regulator-free medium. Regenerated plantlets were acclimatized successfully, and the initial growth of male sterile somatic plants was monitored in the field.

Marker-Assisted Selection of Trees with MALE STERILITY 1 in Cryptomeria japonica D. Don

The practical use of marker-assisted selection (MAS) is limited in conifers because of the difficulty with developing markers due to a rapid decrease in linkage disequilibrium, the limited genomic information available, and the diverse genetic backgrounds among the breeding material collections. First, in this study, two families were produced by artificial crossing between two male-sterile trees, ‘Shindai11’ and ‘Shindai12’, and a plus tree, ‘Suzu-2’ (Ms1/ms1) (S11-S and S12-S families, respectively). The segregation ratio between the male-sterile and male-fertile trees did not deviate significantly from the expected 1:1 ratio in either family. These results clearly suggested that the male-sterile gene of ‘Shindai11’ and ‘Shindai12’ is MALE STERILITY 1 (MS1). Since it is difficult to understand the relative positions of each marker, due to the lack of a linkage map which all the closely linked markers previously reported are mapped on, we constructed a partial linkage map of the region encompassing MS1 using the S11-S and S12-S families. For the S11-S and S12-S families, 19 and 18 markers were mapped onto the partial linkage maps of the MS1 region, respectively. There was collinearity (conserved gene order) between the two partial linkage maps. Two markers (CJt020762_ms1-1 and reCj19250_2335) were mapped to the same position as the MS1 locus on both maps. Of these markers, we used CJt020762 for the MAS in this study. According to the MAS results for 650 trees from six prefectures of Japan (603 trees from breeding materials and 47 trees from the Ishinomaki natural population), five trees in Niigata Prefecture and one tree in Yamagata Prefecture had heterozygous ms1-1, and three trees in Miyagi Prefecture had heterozygous ms1-2. The results obtained in this study suggested that ms1-1 and ms1-2 have different geographical distributions. Since MAS can be used effectively to reduce the labor and time required for selection of trees with a male-sterile gene, the research should help ensure that the quantity of breeding materials will increase to assist future tree-breeding efforts.

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.

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.

Effects of day length- and temperature-regulated genes on annual transcriptome dynamics in Japanese cedar (Cryptomeria japonica D. Don), a gymnosperm indeterminate species

Seasonal phenomena in plants are primarily affected by day length and temperature. The shoot transcriptomes of trees grown in the field and a controlled-environment chamber were compared to characterize genes that control annual rhythms and the effects of day length- and temperature-regulated genes in the gymnosperm Japanese cedar (Cryptomeria japonica D. Don), which exhibits seasonally indeterminate growth. Annual transcriptome dynamics were clearly demonstrated by principal component analysis using microarray data obtained under field-grown conditions. Analysis of microarray data from trees grown in a controlled chamber identified 2,314 targets exhibiting significantly different expression patterns under short-day (SD) and long-day conditions, and 2,045 targets exhibited significantly different expression patterns at 15°C (LT; low temperature) versus 25°C. Interestingly, although growth was suppressed under both SD and LT conditions, approximately 80% of the SD- and LT-regulated targets differed, suggesting that each factor plays a unique role in the annual cycle. The top 1,000 up-regulated targets in the growth/dormant period in the field coincided with more than 50% of the SD- and LT-regulated targets, and gene co-expression network analysis of the annual transcriptome indicated a close relationship between the SD- and LT-regulated targets. These results indicate that the respective effects of day length and temperature interact to control annual transcriptome dynamics. Well-known upstream genes of signaling pathways responsive to environmental conditions, such as the core clock (LHY/CjLHYb and CCA1/CjLHYa) and PEBP family (MFT) genes, exhibited unique expression patterns in Japanese cedar compared with previous reports in other species, suggesting that these genes control differences in seasonal regulation mechanisms between species. The results of this study provide new insights into seasonal regulation of transcription in Japanese cedar.

Construction of genetic linkage map and identification of a novel major locus for resistance to pine wood nematode in Japanese black pine (Pinus thunbergii)

BACKGROUND

Pine wilt disease (PWD), which is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is currently the greatest threat to pine forests in Europe and East Asian countries including Japan. Constructing a detailed linkage map of DNA markers and identifying PWD resistance genes/loci lead to improved resistance in Pinus thunbergii, as well as other Pinus species that are also susceptible to PWD.

RESULTS

A total F₁ mapping population of 188 individuals derived from a cross between the PWD-resistant P. thunbergii varieties 'Tanabe 54' (resistant rank 2 to PWD) and 'Tosashimizu 63' (resistant rank 4 to PWD) was inoculated with PWN, and was evaluated for disease symptoms. To perform linkage analysis for PWN resistance, a set of three maps was constructed; two parental maps generated using the integrated two-way pseudo-testcross method, and a consensus map with population-type cross-pollination. The linkage map of 'Tanabe 54' consisted of 167 loci, and covered 14 linkage groups (LGs), with a total genetic distance of 1214.6 cM. The linkage map of 'Tosashimizu 63' consisted of 252 loci, and covered 14 LGs, with a total genetic distance of 1422.1 cM. The integrated consensus map comprised 12 LGs with the basic chromosome number of P. thunbergii, and a total genetic distance of 1403.6 cM. Results from quantitative trait loci (QTL) analysis using phenotype data and linkage maps indicated that PWN resistance is controlled by a single dominant allele, which was derived from the 'Tanabe 54' female parent. This major QTL was located on linkage group 3 and was designated PWD1 for PINE WILT DISEASE 1.

CONCLUSIONS

The PWD1 locus is a major resistance QTL located on the Pinus consensus LG03 that acts in a dominant manner to confer pine wood nematode resistance. Information from the present study will be useful for P. thunbergii breeding programs to improve resistance to PWD, and also to help identify susceptibility genes in Pinus species.

Genetic Mapping of Prince Rupprecht's Larch (Larix principis-rupprechtii Mayr) by Specific-Locus Amplified Fragment Sequencing

A high-density genetic linkage map is essential for plant genetics and genomics research. However, due to the deficiency of genomic data and high-quality molecular markers, no genetic map has been published for Prince Rupprecht’s larch (Larix principis-rupprechtii Mayr), a conifer species with high ecological and commercial value in northern China. In this study, 145 F1 progeny individuals from an intraspecific cross between two elite clones of L. principis-rupprechtii and their parents were employed to construct the first genetic map in this important tree species using specific-locus amplified fragment sequencing (SLAF-seq). After preprocessing, the procedure yielded 300.20 Gb of raw data containing 1501.22 M pair-end reads. A total of 324,352 SNP markers were detected and 122,785 of them were polymorphic, with a polymorphism rate of 37.86%. Ultimately, 6099 SNPs were organized into a genetic map containing 12 linkage groups, consistent with the haploid chromosome number of larch and most other species in the Pinaceae family. The linkage map spanned 2415.58 cM and covered 99.6% of the L. principis-rupprechtii genome with an average of 0.4 cM between adjacent markers. To the best of our knowledge, this map is the first reference map for L. principis-rupprechtii, as well as the densest one obtained in larch species thus far. The genome-wide SNPs and the high-resolution genetic map will provide a foundation for future quantitative trait loci mapping, map-based cloning, marker-assisted selection, comparative genomics, and genome sequence assembly for larch trees.

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 Specific Nuclear Genetic Loci and Genes That Interact With the Mitochondrial Genome and Contribute to Fecundity in Caenorhabditis elegans

Previous studies have found that fecundity is a multigenic trait regulated, in part, by mitochondrial-nuclear (mit-n) genetic interactions. However, the identification of specific nuclear genetic loci or genes interacting with the mitochondrial genome and contributing to the quantitative trait fecundity is an unsolved issue. Here, a panel of recombinant inbred advanced intercrossed lines (RIAILs), established from a cross between the N2 and CB4856 strains of C. elegans, were used to characterize the underlying genetic basis of mit-n genetic interactions related to fecundity. Sixty-seven single nucleotide polymorphisms (SNPs) were identified by association mapping to be linked with fecundity among 115 SNPs linked to mitotype. This indicated significant epistatic effects between nuclear and mitochondria genetics on fecundity. In addition, two specific nuclear genetic loci interacting with the mitochondrial genome and contributing to fecundity were identified. A significant reduction in fecundity was observed in the RIAILs that carried CB4856 mitochondria and a N2 genotype at locus 1 or a CB4856 genotype at locus 2 relative to the wild-type strains. Then, a hybrid strain (CNC10) was established, which was bred as homoplasmic for the CB4856 mtDNA genome and N2 genotype at locus 1 in the CB4856 nuclear background. The mean fecundity of CNC10 was half the fecundity of the control strain. Several functional characteristics of the mitochondria in CNC10 were also influenced by mit-n interactions. Overall, experimental evidence was presented that specific nuclear genetic loci or genes have interactions with the mitochondrial genome and are associated with fecundity. In total, 18 genes were identified using integrative approaches to have interactions with the mitochondrial genome and to contribute to fecundity.

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.

Transcriptome dynamics of rooting zone and aboveground parts of cuttings during adventitious root formation in Cryptomeria japonica D. Don

BACKGROUND

Adventitious root formation is an essential physiological process for successful propagation of cuttings in various plant species. Because coniferous species are highly heterozygous, propagation of cuttings is of great practical use in breeding. Although various factors influence adventitious root formation, little is known of the associated regulatory mechanisms. Whereas adventitious roots generally form from the base of cuttings, this process is accompanied by physiological changes in leaves, which supply assimilates and metabolites. Herein, we present microarray analyses of transcriptome dynamics during adventitious root formation in whole cuttings in the coniferous species, Cryptomeria japonica.

RESULTS

Temporal patterns of gene expression were determined in the base, the middle, and needles of cuttings at eight time points during adventitious root formation. Global gene expression at the base had diverged from that in the middle by 3-h post-insertion, and changed little in the subsequent 3-days post-insertion, and global gene expression in needles altered characteristically at 3- and 6-weeks post-insertion. In Gene Ontology enrichment analysis of major gene clusters based on hierarchical clustering, the expression profiles of genes related to carbohydrates, plant hormones, and other categories indicated multiple biological changes that were involved in adventitious root formation.

CONCLUSIONS

The present comprehensive transcriptome analyses indicate major transcriptional turning and contribute to the understanding of the biological processes and molecular factors that influence adventitious root formation in C. japonica.

Potential of Genome-Wide Studies in Unrelated Plus Trees of a Coniferous Species, Cryptomeria japonica (Japanese Cedar)

A genome-wide association study (GWAS) was conducted on more than 30,000 single nucleotide polymorphisms (SNPs) in unrelated first-generation plus tree genotypes from three populations of Japanese cedar Cryptomeria japonica D. Don with genomic prediction for traits of growth, wood properties and male fecundity. Among the assessed populations, genetic characteristics including the extent of linkage disequilibrium (LD) and genetic structure differed and these differences are considered to be due to differences in genetic background. Through population-independent GWAS, several significant SNPs found close to the regions associated with each of these traits and shared in common across the populations were identified. The accuracies of genomic predictions were dependent on the traits and populations and reflected the genetic architecture of traits and genetic characteristics. Prediction accuracies using SNPs selected based on GWAS results were similar to those using all SNPs for several combinations of traits and populations. We discussed the application of genome-wide studies for C. japonica improvement.