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

Climate-Associated Genetic Variation and Projected Genetic Offsets for Cryptomeria japonica D. Don Under Future Climate Scenarios

Revealing the spatial distribution of adaptive genetic variation is both a challenging and crucial task in evolutionary ecology, essential for understanding local adaptation within species, and in management, for predicting species responses to future climate change. This understanding is particularly important for long-lived tree species, which may not be able to migrate quickly enough to adapt to rapid climate changes and may need to rely on their standing genetic variation. In this study, we focused on Cryptomeria japonica, a major component of Japan's temperate forests and an important forestry species adapted to the humid environment of monsoon Asia. We extracted climate-associated genetic variation from the entire genome and evaluated its distribution and vulnerability under future climate scenarios using spatial modeling techniques. We analyzed 31,676 high-quality SNPs from 249 individuals across 22 natural populations of C. japonica, covering its entire distribution range. We identified 239 candidate climate-associated SNPs and found winter temperature, summer precipitation, and winter precipitation as the most significant factors explaining the genetic variation in these SNPs. The climate-associated genetic variation deviated from non-associated (neutral) genetic variation in the opposite (the Sea of Japan and Pacific Ocean) sides of Japanese archipelago, suggesting natural selection of different climate conditions in these regions. Difference in estimated allele frequency at the climate-associated loci (genetic offset) between the present and future (2090 in the SSP5-8.5 scenario) climate conditions was predicted to be larger in three areas (not only southwestern Japan but also coastal area on the Sea of Japan side and inland area on the Pacific Ocean side in northeastern Japan). This prediction implies the discrepancy between standing genetic variation at the present and that adaptive to the future climate in these areas, which underscores the necessity for proactive management to adjust the adaptive genetic variation.

A chromosome-level genome assembly of a model conifer plant, the Japanese cedar, Cryptomeria japonica D. Don

BACKGROUND

The Japanese cedar (Cryptomeria japonica D. Don) is one of the most important Japanese forest trees, occupying approximately 44% of artificial forests and planted in East Asia, the Azores Archipelago, and certain islands in the Indian Ocean. Although the huge genome of the species (ca. 9 Gbp) with abundant repeat elements may have represented an obstacle for genetic analysis, this species is easily propagated by cutting, flowered by gibberellic acid, transformed by Agrobacterium, and edited by CRISPR/Cas9. These characteristics of C. japonica recommend it as a model conifer species for which reference genome sequences are necessary.

RESULTS

Herein, we report the first chromosome-level assembly of C. japonica (2n = 22) using third-generation selfed progeny (estimated homozygosity rate = 0.96). Young leaf tissue was used to extract high molecular weight DNA (> 50 kb) for HiFi PacBio long-read sequencing and to construct an Hi-C/Omni-C library for Illumina short-read sequencing. The 29× and 26× genome coverage of HiFi and Illumina reads, respectively, for de novo assembly yielded 2,651 contigs (9.1 Gbp, N50 contig size 12.0 Mbp). Hi-C analysis mapped 97% of the nucleotides on 11 chromosomes. The assembly was verified through comparison with a consensus linkage map comprising 7,781 markers. BUSCO analysis identified ∼ 91% conserved genes.

CONCLUSIONS

Annotations of genes and comparisons of repeat elements with other Cupressaceae and Pinaceae species provide a fundamental resource for conifer research.

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

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.

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.

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.

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.

High-density genetic map construction and quantitative trait loci identification for growth traits in (Taxodium distichum var. distichum × T. mucronatum) × T. mucronatum

BACKGROUND

'Zhongshanshan' is the general designation for the superior interspecific hybrid clones of Taxodium species, which is widely grown for economic and ecological purposes in southern China. Growth is the priority objective in 'Zhongshanshan' tree improvement. A high-density linkage map is vital to efficiently identify key quantitative trait loci (QTLs) that affect growth.

RESULTS

In total, 403.16 Gb of data, containing 2016,336 paired-end reads, was obtained after preprocessing. The average sequencing depth was 28.49 in T. distichum var. distichum, 25.18 in T. mucronatum, and 11.12 in each progeny. In total, 524,662 high-quality SLAFs were detected, of which 249,619 were polymorphic, and 6166 of the polymorphic markers met the requirements for use in constructing a genetic map. The final map harbored 6156 SLAF markers on 11 linkage groups, and was 1137.86 cM in length, with an average distance of 0.18 cM between adjacent markers. Separate QTL analyses of traits in different years by CIM detected 7 QTLs. While combining multiple-year data, 13 QTLs were detected by ICIM. 5 QTLs were repeatedly detected by the two methods, and among them, 3 significant QTLs (q6-2, q4-2 and q2-1) were detected in at least two traits. Bioinformatic analysis discoveried a gene annotated as a leucine-rich repeat receptor-like kinase gene within q4-2.

CONCLUSIONS

This map is the most saturated one constructed in a Taxodiaceae species to date, and would provide useful information for future comparative mapping, genome assembly, and marker-assisted selection.

Dissecting Key Adaptation Traits in the Polyploid Perennial Medicago sativa Using GBS-SNP Mapping

Understanding key adaptation traits is crucial to developing new cultivars with broad adaptations. The main objective of this research is to understand the genetic basis of winter hardiness (WH) and fall dormancy (FD) in alfalfa and the association between the two traits. QTL analysis was conducted in a pseudo-testcross F1 population developed from two cultivars contrasting in FD (3010 with FD = 2 and CW 1010 with FD = 10). The mapping population was evaluated in three replications at two locations (Watkinsville and Blairsville, GA). FD levels showed low to moderate correlations with WH (0.22-0.57). Assessing dormancy in winter is more reliable than in the fall in southern regions with warm winters. The mapping population was genotyped using Genotyping-by-sequencing (GBS). Single dose allele SNPs (SDA) were used for constructing linkage maps. The parental map (CW 1010) consisted of 32 linkage groups spanning 2127.5 cM with 1377 markers and an average marker density of 1.5 cM/SNP. The maternal map (3010) had 32 linkage groups spanning 2788.4 cM with 1837 SDA SNPs with an average marker density of 1.5 cM/SNP. Forty-five significant (P < 0.05) QTLs for FD and 35 QTLs for WH were detected on both male and female linkage maps. More than 75% (22/28) of the dormancy QTL detected from the 3010 parent did not share genomic regions with WH QTLs and more than 70% (12/17) dormancy QTLs detected from CW 1010 parent were localized in different genomic regions than WH QTLs. These results suggest that the two traits have independent inheritance and therefore can be improved separately in breeding programs.

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.

High density SNP and DArT-based genetic linkage maps of two closely related oil palm populations

Oil palm (Elaeis guineensis Jacq.) is an outbreeding perennial tree crop with long breeding cycles, typically 12 years. Molecular marker technologies can greatly improve the breeding efficiency of oil palm. This study reports the first use of the DArTseq platform to genotype two closely related self-pollinated oil palm populations, namely AA0768 and AA0769 with 48 and 58 progeny respectively. Genetic maps were constructed using the DArT and SNP markers generated in combination with anchor SSR markers. Both maps consisted of 16 major independent linkage groups (2n = 2× = 32) with 1399 and 1466 mapped markers for the AA0768 and AA0769 populations, respectively, including the morphological trait "shell-thickness" (Sh). The map lengths were 1873.7 and 1720.6 cM with an average marker density of 1.34 and 1.17 cM, respectively. The integrated map was 1803.1 cM long with 2066 mapped markers and average marker density of 0.87 cM. A total of 82% of the DArTseq marker sequence tags identified a single site in the published genome sequence, suggesting preferential targeting of gene-rich regions by DArTseq markers. Map integration of higher density focused around the Sh region identified closely linked markers to the Sh, with D.15322 marker 0.24 cM away from the morphological trait and 5071 bp from the transcriptional start of the published SHELL gene. Identification of the Sh marker demonstrates the robustness of using the DArTseq platform to generate high density genetic maps of oil palm with good genome coverage. Both genetic maps and integrated maps will be useful for quantitative trait loci analysis of important yield traits as well as potentially assisting the anchoring of genetic maps to genomic sequences.

New statistical methods for estimation of recombination fractions in F2 population

Background

Dominant markers in an F₂ population or a hybrid population have much less linkage information in repulsion phase than in coupling phase. Linkage analysis produces two separate complementary marker linkage maps that have little use in disease association analysis and breeding. There is a need to develop efficient statistical methods and computational algorithms to construct or merge a complete linkage dominant marker maps. The key for doing so is to efficiently estimate recombination fractions between dominant markers in repulsion phases.

Result

We proposed an expectation least square (ELS) algorithm and binomial analysis of three-point gametes (BAT) for estimating gamete frequencies from F₂ dominant and codominant marker data, respectively. The results obtained from simulated and real genotype datasets showed that the ELS algorithm was able to accurately estimate frequencies of gametes and outperformed the EM algorithm in estimating recombination fractions between dominant loci and recovering true linkage maps of 6 dominant loci in coupling and unknown linkage phases. Our BAT method also had smaller variances in estimation of two-point recombination fractions than the EM algorithm.

Conclusion

ELS is a powerful method for accurate estimation of gamete frequencies in dominant three-locus system in an F₂ population and BAT is a computationally efficient and fast method for estimating frequencies of three-point codominant gametes.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1804-8) contains supplementary material, which is available to authorized users.

Genetic linkage map construction and QTL mapping of seedling height, basal diameter and crown width of Taxodium 'Zhongshanshan 302' × T. mucronatum

Taxodium is a genus renowned for its fast growth, good form and tolerance of flooding, salt, alkalinity, disease and strong winds. In this study, a genetic linkage map was constructed using sequence-related amplified polymorphism (SRAP) and simple sequence repeat (SSR) markers based on an F₁ population containing 148 individuals generated from a cross between T. ‘Zhongshanshan 302’ and T. mucronatum. The map has a total length of 976.5 cM, with a mean distance of 7.0 cM between markers, and contains 34 linkage groups with 179 markers (171 SRAPs and 8 SSRs). Quantitative trait loci (QTLs) affecting growth traits, such as seedling height, basal diameter and crown width, were detected based on the constructed linkage map. Four significant QTLs were identified, three of which, namely qtSH-1 for seedling height, qtBD-1 for basal diameter and qtCW-1 for crown width, were located at 2.659 cM of LG7 with logarithm odds values of 3.72, 3.49 and 3.93, respectively, and explained 24.9, 27.0 and 21.7 % of the total variation of the three grown traits, respectively. Another QTL for crown width (qtCW-2) was detected at 1.0 cM on LG13, with a logarithm of odds value of 3.15, and explained 31.7 % of the total variation of crown width. This is the first report on the construction of a genetic linkage map and QTL analysis in Taxodium, laying the groundwork for the construction of a high-density genetic map and QTL mapping in the genus Taxodium.

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).

Comparative Genomics Analyses Reveal Extensive Chromosome Colinearity and Novel Quantitative Trait Loci in Eucalyptus

Dense genetic maps, along with quantitative trait loci (QTLs) detected on such maps, are powerful tools for genomics and molecular breeding studies. In the important woody genus Eucalyptus, the recent release of E. grandis genome sequence allows for sequence-based genomic comparison and searching for positional candidate genes within QTL regions. Here, dense genetic maps were constructed for E. urophylla and E. tereticornis using genomic simple sequence repeats (SSR), expressed sequence tag (EST) derived SSR, EST-derived cleaved amplified polymorphic sequence (EST-CAPS), and diversity arrays technology (DArT) markers. The E. urophylla and E. tereticornis maps comprised 700 and 585 markers across 11 linkage groups, totaling at 1,208.2 and 1,241.4 cM in length, respectively. Extensive synteny and colinearity were observed as compared to three earlier DArT-based eucalypt maps (two maps with E. grandis × E. urophylla and one map of E. globulus) and with the E. grandis genome sequence. Fifty-three QTLs for growth (10-56 months of age) and wood density (56 months) were identified in 22 discrete regions on both maps, in which only one colocalizaiton was found between growth and wood density. Novel QTLs were revealed as compared with those previously detected on DArT-based maps for similar ages in Eucalyptus. Eleven to 585 positional candidate genes were obained for a 56-month-old QTL through aligning QTL confidence interval with the E. grandis genome. These results will assist in comparative genomics studies, targeted gene characterization, and marker-assisted selection in Eucalyptus and the related taxa.

Evidence of intense chromosomal shuffling during conifer evolution

Although recent advances have been gained on genome evolution in angiosperm lineages, virtually nothing is known about karyotype evolution in the other group of seed plants, the gymnosperms. Here, we used high-density gene-based linkage mapping to compare the karyotype structure of two families of conifers (the most abundant group of gymnosperms) separated around 290 Ma: Pinaceae and Cupressaceae. We propose for the first time a model based on the fusion of 20 ancestral chromosomal blocks that may have shaped the modern karyotpes of Pinaceae (with n = 12) and Cupressaceae (with n = 11). The considerable difference in modern genome organization between these two lineages contrasts strongly with the remarkable level of synteny already reported within the Pinaceae. It also suggests a convergent evolutionary mechanism of chromosomal block shuffling that has shaped the genomes of the spermatophytes.

Genetic differentiation and evolutionary adaptation in Cryptomeria japonica

Local adaptation of plant species is a central issue for survival during global climate change, especially for long-lived forest trees, with their lengthy regeneration time and spatially limited gene flow. Identification of loci and/or genomic regions associated with local adaptation is necessary for knowledge of both evolution and molecular breeding for climate change. Cryptomeria japonica is an important species for forestry in Japan; it has a broad natural distribution and can survive in a range of different environments. The genetic structure of 14 natural populations of this species was investigated using 3930 SNP markers. Populations on the Pacific Ocean side of Japan are clearly different from those on the Japan Sea side, as discussed in previous studies. Structure analysis and population network trees show that peripheral populations, including the most northerly and southerly ones, have unique features. We found that the genetic differentiation coefficient is low, F_ST = 0.05, although it must account for the presence of important genes associated with adaptation to specific environments. In total, 208 outlier loci were detected, of which 43 were associated with environmental variables. Four clumped regions of outlier loci were detected in the genome by linkage analysis. Linkage disequilibrium (LD) was quite high in these clumps of outlier loci, which were found in linkage groups (LGs) 2, 7, 10, and 11, especially between populations of two varieties, and when interchromosomal LD was also detected. The LG7 region is characteristic of the Yakushima population, which is a large, isolated, peripheral population occupying a specific environment resulting from isolation combined with volcanic activity in the region. The detected LD may provide strong evidence for selection between varieties.

A Picea abies linkage map based on SNP markers identifies QTLs for four aspects of resistance to Heterobasidion parviporum infection

A consensus linkage map of Picea abies, an economically important conifer, was constructed based on the segregation of 686 SNP markers in a F1 progeny population consisting of 247 individuals. The total length of 1889.2 cM covered 96.5% of the estimated genome length and comprised 12 large linkage groups, corresponding to the number of haploid P. abies chromosomes. The sizes of the groups (from 5.9 to 9.9% of the total map length) correlated well with previous estimates of chromosome sizes (from 5.8 to 10.8% of total genome size). Any locus in the genome has a 97% probability to be within 10 cM from a mapped marker, which makes the map suited for QTL mapping. Infecting the progeny trees with the root rot pathogen Heterobasidion parviporum allowed for mapping of four different resistance traits: lesion length at the inoculation site, fungal spread within the sapwood, exclusion of the pathogen from the host after initial infection, and ability to prevent the infection from establishing at all. These four traits were associated with two, four, four and three QTL regions respectively of which none overlapped between the traits. Each QTL explained between 4.6 and 10.1% of the respective traits phenotypic variation. Although the QTL regions contain many more genes than the ones represented by the SNP markers, at least four markers within the confidence intervals originated from genes with known function in conifer defence; a leucoanthocyanidine reductase, which has previously been shown to upregulate during H. parviporum infection, and three intermediates of the lignification process; a hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase, a 4-coumarate CoA ligase, and a R2R3-MYB transcription factor.

Genetic control of functional traits related to photosynthesis and water use efficiency in Pinus pinaster Ait. drought response: integration of genome annotation, allele association and QTL detection for candidate gene identification

BACKGROUND

Understanding molecular mechanisms that control photosynthesis and water use efficiency in response to drought is crucial for plant species from dry areas. This study aimed to identify QTL for these traits in a Mediterranean conifer and tested their stability under drought.

RESULTS

High density linkage maps for Pinus pinaster were used in the detection of QTL for photosynthesis and water use efficiency at three water irrigation regimes. A total of 28 significant and 27 suggestive QTL were found. QTL detected for photochemical traits accounted for the higher percentage of phenotypic variance. Functional annotation of genes within the QTL suggested 58 candidate genes for the analyzed traits. Allele association analysis in selected candidate genes showed three SNPs located in a MYB transcription factor that were significantly associated with efficiency of energy capture by open PSII reaction centers and specific leaf area.

CONCLUSIONS

The integration of QTL mapping of functional traits, genome annotation and allele association yielded several candidate genes involved with molecular control of photosynthesis and water use efficiency in response to drought in a conifer species. The results obtained highlight the importance of maintaining the integrity of the photochemical machinery in P. pinaster drought response.

Use of EST-SSR markers for evaluating genetic diversity and fingerprinting celery (Apium graveolens L.) cultivars

Celery (Apium graveolens L.) is one of the most economically important vegetables worldwide, but genetic and genomic resources supporting celery molecular breeding are quite limited, thus few studies on celery have been conducted so far. In this study we made use of simple sequence repeat (SSR) markers generated from previous celery transcriptome sequencing and attempted to detect the genetic diversity and relationships of commonly used celery accessions and explore the efficiency of the primers used for cultivars identification. Analysis of molecular variance (AMOVA) of Apium graveolens L. var. dulce showed that approximately 43% of genetic diversity was within accessions, 45% among accessions, and 22% among horticultural types. The neighbor-joining tree generated by unweighted pair group method with arithmetic mean (UPGMA), and population structure analysis, as well as principal components analysis (PCA), separated the cultivars into clusters corresponding to the geographical areas where they originated. Genetic distance analysis suggested that genetic variation within Apium graveolens was quite limited. Genotypic diversity showed any combinations of 55 genic SSRs were able to distinguish the genotypes of all 30 accessions.

Demonstration of genome-wide association studies for identifying markers for wood property and male strobili traits in Cryptomeria japonica

Genome-wide association studies (GWAS) are an alternative to bi-parental QTL mapping in long-lived perennials. In the present study, we examined the potential of GWAS in conifers using 367 unrelated plus trees of Cryptomeria japonica D. Don, which is the most widely planted and commercially important tree species in Japan, and tried to detect significant associations between wood property traits and quantity of male strobili on the one hand, and 1,032 single nucleotide polymorphisms (SNPs) assigned to 1,032 genes on the other. Association analysis was performed with the mixed linear model taking into account kinship relationships and subpopulation structure. In total, 6 SNPs were found to have significant associations with the variations in phenotype. These SNPs were not associated with the positions of known genes and QTLs that have been reported to date, thus they may identify novel QTLs. These 6 SNPs were all found in sequences showing similarities with known genes, although further analysis is required to dissect the ways in which they affect wood property traits and abundance of male strobili. These presumptive QTL loci provide opportunities for improvement of C. japonica, based on a marker approach. The results suggest that GWAS has potential for use in future breeding programs in C. japonica.

High-density linkage mapping in a pine tree reveals a genomic region associated with inbreeding depression and provides clues to the extent and distribution of meiotic recombination

BACKGROUND

The availability of a large expressed sequence tags (EST) resource and recent advances in high-throughput genotyping technology have made it possible to develop highly multiplexed SNP arrays for multi-objective genetic applications, including the construction of meiotic maps. Such approaches are particularly useful in species with a large genome size, precluding the use of whole-genome shotgun assembly with current technologies.

RESULTS

In this study, a 12 k-SNP genotyping array was developed for maritime pine from an extensive EST resource assembled into a unigene set. The offspring of three-generation outbred and inbred mapping pedigrees were then genotyped. The inbred pedigree consisted of a classical F2 population resulting from the selfing of a single inter-provenance (Landes x Corsica) hybrid tree, whereas the outbred pedigree (G2) resulted from a controlled cross of two intra-provenance (Landes x Landes) hybrid trees. This resulted in the generation of three linkage maps based on SNP markers: one from the parental genotype of the F2 population (1,131 markers in 1,708 centimorgan (cM)), and one for each parent of the G2 population (1,015 and 1,110 markers in 1,447 and 1,425 cM for the female and male parents, respectively). A comparison of segregation patterns in the progeny obtained from the two types of mating (inbreeding and outbreeding) led to the identification of a chromosomal region carrying an embryo viability locus with a semi-lethal allele. Following selfing and segregation, zygote mortality resulted in a deficit of Corsican homozygous genotypes in the F2 population. This dataset was also used to study the extent and distribution of meiotic recombination along the length of the chromosomes and the effect of sex and/or genetic background on recombination. The genetic background of trees in which meiotic recombination occurred was found to have a significant effect on the frequency of recombination. Furthermore, only a small proportion of the recombination hot- and cold-spots were common to all three genotypes, suggesting that the spatial pattern of recombination was genetically variable.

CONCLUSION

This study led to the development of classical genomic tools for this ecologically and economically important species. It also identified a chromosomal region bearing a semi-lethal recessive allele and demonstrated the genetic variability of recombination rate over the genome.

Development of high-throughput SNP-based genotyping in Acacia auriculiformis x A. mangium hybrids using short-read transcriptome data

Background

Next Generation Sequencing has provided comprehensive, affordable and high-throughput DNA sequences for Single Nucleotide Polymorphism (SNP) discovery in Acacia auriculiformis and Acacia mangium. Like other non-model species, SNP detection and genotyping in Acacia are challenging due to lack of genome sequences. The main objective of this study is to develop the first high-throughput SNP genotyping assay for linkage map construction of A. auriculiformis x A. mangium hybrids.

Results

We identified a total of 37,786 putative SNPs by aligning short read transcriptome data from four parents of two Acacia hybrid mapping populations using Bowtie against 7,839 de novo transcriptome contigs. Given a set of 10 validated SNPs from two lignin genes, our in silico SNP detection approach is highly accurate (100%) compared to the traditional in vitro approach (44%). Further validation of 96 SNPs using Illumina GoldenGate Assay gave an overall assay success rate of 89.6% and conversion rate of 37.5%. We explored possible factors lowering assay success rate by predicting exon-intron boundaries and paralogous genes of Acacia contigs using Medicago truncatula genome as reference. This assessment revealed that presence of exon-intron boundary is the main cause (50%) of assay failure. Subsequent SNPs filtering and improved assay design resulted in assay success and conversion rate of 92.4% and 57.4%, respectively based on 768 SNPs genotyping. Analysis of clustering patterns revealed that 27.6% of the assays were not reproducible and flanking sequence might play a role in determining cluster compression. In addition, we identified a total of 258 and 319 polymorphic SNPs in A. auriculiformis and A. mangium natural germplasms, respectively.

Conclusion

We have successfully discovered a large number of SNP markers in A. auriculiformis x A. mangium hybrids using next generation transcriptome sequencing. By using a reference genome from the most closely related species, we converted most SNPs to successful assays. We also demonstrated that Illumina GoldenGate genotyping together with manual clustering can provide high quality genotypes for a non-model species like Acacia. These SNPs markers are not only important for linkage map construction, but will be very useful for hybrid discrimination and genetic diversity assessment of natural germplasms in the future.

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.

Annotated genetic linkage maps of Pinus pinaster Ait. from a Central Spain population using microsatellite and gene based markers

BACKGROUND

Pinus pinaster Ait. is a major resin producing species in Spain. Genetic linkage mapping can facilitate marker-assisted selection (MAS) through the identification of Quantitative Trait Loci and selection of allelic variants of interest in breeding populations. In this study, we report annotated genetic linkage maps for two individuals (C14 and C15) belonging to a breeding program aiming to increase resin production. We use different types of DNA markers, including last-generation molecular markers.

RESULTS

We obtained 13 and 14 linkage groups for C14 and C15 maps, respectively. A total of 211 and 215 markers were positioned on each map and estimated genome length was between 1,870 and 2,166 cM respectively, which represents near 65% of genome coverage. Comparative mapping with previously developed genetic linkage maps for P. pinaster based on about 60 common markers enabled aligning linkage groups to this reference map. The comparison of our annotated linkage maps and linkage maps reporting QTL information revealed 11 annotated SNPs in candidate genes that co-localized with previously reported QTLs for wood properties and water use efficiency.

CONCLUSIONS

This study provides genetic linkage maps from a Spanish population that shows high levels of genetic divergence with French populations from which segregating progenies have been previously mapped. These genetic maps will be of interest to construct a reliable consensus linkage map for the species. The importance of developing functional genetic linkage maps is highlighted, especially when working with breeding populations for its future application in MAS for traits of interest.

Genome scanning for detecting adaptive genes along environmental gradients in the Japanese conifer, Cryptomeria japonica

Local adaptation is important in evolutionary processes and speciation. We used multiple tests to identify several candidate genes that may be involved in local adaptation from 1026 loci in 14 natural populations of Cryptomeria japonica, the most economically important forestry tree in Japan. We also studied the relationships between genotypes and environmental variables to obtain information on the selective pressures acting on individual populations. Outlier loci were mapped onto a linkage map, and the positions of loci associated with specific environmental variables are considered. The outlier loci were not randomly distributed on the linkage map; linkage group 11 was identified as a genomic island of divergence. Three loci in this region were also associated with environmental variables such as mean annual temperature, daily maximum temperature, maximum snow depth, and so on. Outlier loci identified with high significance levels will be essential for conservation purposes and for future work on molecular breeding.