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

The revelation of genomic breed composition using target capture sequencing: a case of Taxodium

Taxodium plants have good flood tolerance and thus were introduced into China from North America in the early 1900s. The subsequent decades of cross-breeding experiments within Taxodium have produced many new hybrid cultivars in China while also creating confusion in the genetic background of Taxodium plants. In the present study, target capture sequencing-derived SNP markers were used to reveal the genomic composition of different species and cultivars of Taxodium. The results unraveled the phylogenetic relationship within Taxodium and suggested the possibility of recent interspecific natural hybridization events. The introduced (Chinese) population is genetically similar to the native (North American) population, thus providing genetic evidence for historical introductions of Taxodium. Hybrid categories of different cultivars of Taxodium hybrid 'Zhongshanshan' were further identified, and their differences in parentage were revealed. Collectively, the findings provide evidence for understanding the genetics and hybridization of Taxodium and shed light on the future breeding and cultivation of cultivars with great ecological and economic potential.

Novel molecular markers for Taxodium breeding from the chloroplast genomes of four artificial Taxodium hybrids

Taxodium "Zhongshanshan" are a group of intraspecific Taxodium hybrids with superparental dominance and high ecological and economic value in southern China. Identifying the parentage of hybrids, especially the male parent, is critically important for genetic studies. However, the large nuclear genomes of members of the genus Taxodium pose a major challenge for the development of molecular markers. Here, we developed novel molecular markers by conducting a comparative analysis of the chloroplast genomes of four artificial Taxodium hybrids and their parents. The lengths of the whole chloroplast genome ranged from 131,942 to 132,128 bp, and the total guanine (GC) content of the chloroplast genomes ranged from 34.6% to 35.81%. A total of 120 unique genes were identified, including 83 protein-coding genes, 33 transfer RNAs, and four ribosomal RNAs. There were 69-71 simple sequence repeats were detected in the four hybrids. Phylogenetic analysis revealed that these hybrids clustered with their paternal parents. Similar findings were obtained by analysis of the GC content of protein-coding genes. Molecular markers were developed using the highly variable regions of the chloroplast genomes, and polymerase chain reaction (PCR) assays revealed that these markers were effective for identifying the male parents of these hybrids. Our findings indicate for the first time that the chloroplast genomes of Taxodium are paternally inherited. Generally, these molecular markers could facilitate breeding and genetic studies of Taxodium.

Identification of quantitative trait loci (QTLs) and candidate genes of seed Iron and zinc content in soybean [Glycine max (L.) Merr.]

BACKGROUND

Deciphering the hereditary mechanism of seed iron (Fe) and zinc (Zn) content in soybean is important and sustainable to address the "hidden hunger" that presently affects approximately 2 billion people worldwide. Therefore, in order to detect genomic regions related to soybean seed Fe and Zn content, a recombinant inbred line (RIL) population with 248 lines was assessed in four environments to detect Quantitative Trait Loci (QTLs) related to soybean seed Fe and Zn content.

RESULT

Wide variation was found in seed Fe and Zn content in four environments, and genotype, environment, and genotype × environment interactions had significant influences on both the seed Fe and Zn content. A positive correlation was observed between seed Fe content and seed Zn content, and broad-sense heritability (H2) of seed Fe and Zn content were 0.73 and 0.75, respectively. In this study, five QTLs for seed Fe content were detected with 4.57 - 32.71% of phenotypic variation explained (PVE) and logarithm of odds (LOD) scores ranging from 3.60 to 33.79. Five QTLs controlling the seed Zn content were detected, and they individually explained 3.35 to 26.48% of the phenotypic variation, with LOD scores ranging from 3.64 to 20.4. Meanwhile, 409,541 high-quality single-nucleotide variants (SNVs) and 85,102 InDels (except intergenic regions) between two bi-parental lines were identified by whole genome resequencing. A total of 12 candidate genes were reported in one major QTL for seed Fe content and two major QTLs for seed Zn content, with the help of RNA-Seq analysis, gene ontology (GO) enrichment, gene annotation, and bi-parental whole genome sequencing (WGS) data.

CONCLUSIONS

Limited studies were performed about microelement of soybean, so these results may play an important role in the biofortification of Fe and Zn and accelerate the development of marker-assisted selection (MAS) for breeding soybeans fortified with iron and zinc.

Genome Wide Association Study Identifies Candidate Genes Related to the Earlywood Tracheid Properties in Picea crassifolia Kom

Picea crassifolia Kom. is one of the timber and ecological conifers in China and its wood tracheid traits directly affect wood formation and adaptability under harsh environment. Molecular studies on P. crassifolia remain inadequate because relatively few genes have been associated with these traits. To identify markers and candidate genes that can potentially be used for genetic improvement of wood tracheid traits, we examined 106 clones of P. crassifolia, and investigated phenotypic data for 14 wood tracheid traits before specific-locus amplified fragment sequencing (SLAF-seq) was employed to perform a genome wide association study (GWAS). Subsequently, the results were used to screen single nucleotide polymorphism (SNP) loci and candidate genes that exhibited a significant correlation with the studied traits. We developed 4,058,883 SLAF-tags and 12,275,765 SNP loci, and our analyses identified a total of 96 SNP loci that showed significant correlations with three earlywood tracheid traits using a mixed linear model (MLM). Next, candidate genes were screened in the 100 kb zone (50 kb upstream, 50 kb downstream) of each of the SNP loci, whereby 67 candidate genes were obtained in earlywood tracheid traits, including 34 genes of known function and 33 genes of unknown function. We provide the most significant SNP for each trait-locus combination and candidate genes occurring within the GWAS hits. These resources provide a foundation for the development of markers that could be used in wood traits improvement and candidate genes for the development of earlywood tracheid in P. crassifolia.

High-density genetic linkage mapping reveals low stability of QTLs across environments for economic traits in Eucalyptus

Introduction

Eucalyptus urophylla, E. tereticornis and their hybrids are the most important commercial forest tree species in South China where they are grown for pulpwood and solid wood production. Construction of a fine-scale genetic linkage map and detecting quantitative trait loci (QTL) for economically important traits linked to these end-uses will facilitate identification of the main candidate genes and elucidate the regulatory mechanisms.

Method

A high-density consensus map (a total of 2754 SNPs with 1359.18 cM) was constructed using genotyping by sequencing (GBS) on clonal progenies of E. urophylla × tereticornis hybrids. QTL mapping of growth and wood property traits were conducted in three common garden experiments, resulting in a total of 108 QTLs. A total of 1052 candidate genes were screened by the efficient combination of QTL mapping and transcriptome analysis.

Results

Only ten QTLs were found to be stable across two environments, and only one (qSG10Stable mapped on chromosome 10, and associated with lignin syringyl-to-guaiacyl ratio) was stable across all three environments. Compared to other QTLs, qSG10Stable explained a very high level of phenotypic variation (18.4-23.6%), perhaps suggesting that QTLs with strong effects may be more stably inherited across multiple environments. Screened candidate genes were associated with some transcription factor families, such as TALE, which play an important role in the secondary growth of plant cell walls and the regulation of wood formation.

Discussion

While QTLs such as qSG10Stable, found to be stable across three sites, appear to be comparatively uncommon, their identification is likely to be a key to practical QTL-based breeding. Further research involving clonally-replicated populations, deployed across multiple target planting sites, will be required to further elucidate QTL-by-environment interactions.

Genetic dissection of growth trajectories in forest trees: From FunMap to FunGraph

Growth is the developmental process involving important genetic components. Functional mapping (FunMap) has been used as an approach to map quantitative trait loci (QTLs) governing growth trajectories by incorporating growth equations. FunMap is based on reductionism thinking, with a power to identify a small set of significant QTLs from the whole pool of genome-wide markers. Yet, increasing evidence shows that a complex trait is controlled by all genes the organism may possibly carry. Here, we describe and demonstrate a different mapping approach that encapsulates all markers into genetic interaction networks. This approach, symbolized as FunGraph, combines functional mapping, evolutionary game theory, and prey-predator theory into mathematical graphs, allowing the observed genetic effect of a locus to be decomposed into its independent component (resulting from this locus' intrinsic capacity) and dependent component (due to extrinsic regulation by other loci). Using FunGraph, we can visualize and trace the roadmap of how each locus interact with every other locus to impact growth. In a population-based association study of Euphrates poplar, we use FunGraph to identify the previously neglected genetic interaction effects that contribute to the genetic architecture of juvenile stem growth. FunGraph could open up a novel gateway to comprehend the global genetic control mechanisms of complex traits.

Genetic Dissection of Growth and Eco-Physiological Traits Associated with Altitudinal Adaptation in Sakhalin Fir (Abies sachalinensis) Based on QTL Mapping

(1) Background: The genetic basis of local adaptation in conifers remains poorly understood because of limited research evidence and the lack of suitable genetic materials. Sakhalin fir (Abies sachalinensis) is an ideal organism for elucidating the genetic basis of local adaptation because its altitudinal adaptation has been demonstrated, and suitable materials for its linkage mapping are available. (2) Method: We constructed P336 and P236 linkage maps based on 486 and 516 single nucleotide polymorphisms, respectively, that were derived from double digest restriction site-associated DNA sequences. We measured the growth and eco-physiological traits associated with morphology, phenology, and photosynthesis, which are considered important drivers of altitudinal adaptation. (3) Results: The quantitative trait loci (QTLs) for growth traits, phenology, needle morphology, and photosynthetic traits were subsequently detected. Similar to previous studies on conifers, most traits were controlled by multiple QTLs with small or moderate effects. Notably, we detected that one QTL for the crown area might be a type-A response regulator, a nuclear protein responsible for the cytokinin-induced shoot elongation. (4) Conclusion: The QTLs detected in this study include potentially important genomic regions linked to altitudinal adaptation in Sakhalin fir.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

Construction of a high-density genetic map and QTL mapping of leaf traits and plant growth in an interspecific F1 population of Catalpa bungei × Catalpa duclouxii Dode

BACKGROUND

Catalpa bungei is an important tree species used for timber in China and widely cultivated for economic and ornamental purposes. A high-density linkage map of C. bungei would be an efficient tool not only for identifying key quantitative trait loci (QTLs) that affect important traits, such as plant growth and leaf traits, but also for other genetic studies.

RESULTS

Restriction site-associated DNA sequencing (RAD-seq) was used to identify molecular markers and construct a genetic map. Approximately 280.77 Gb of clean data were obtained after sequencing, and in total, 25,614,295 single nucleotide polymorphisms (SNPs) and 2,871,647 insertions-deletions (InDels) were initially identified in the genomes of 200 individuals of a C. bungei (7080) × Catalpa duclouxii (16-PJ-3) F₁ population and their parents. Finally, 9072 SNP and 521 InDel markers that satisfied the requirements for constructing a genetic map were obtained. The integrated genetic map contained 9593 pleomorphic markers in 20 linkage groups and spanned 3151.63 cM, with an average distance between adjacent markers of 0.32 cM. Twenty QTLs for seven leaf traits and 13 QTLs for plant height at five successive time points were identified using our genetic map by inclusive composite interval mapping (ICIM). Q16-60 was identified as a QTL for five leaf traits, and three significant QTLs (Q9-1, Q18-66 and Q18-73) associated with plant growth were detected at least twice. Genome annotation suggested that a cyclin gene participates in leaf trait development, while the growth of C. bungei may be influenced by CDC48C and genes associated with phytohormone synthesis.

CONCLUSIONS

This is the first genetic map constructed in C. bungei and will be a useful tool for further genetic study, molecular marker-assisted breeding and genome assembly.

Identification and functional analysis of LecRLK genes in Taxodium 'Zhongshanshan'

Background

Lectin receptor-like protein kinases (LecRLKs) can transform external stimuli into intracellular signals and play important regulatory roles in plant development and response to environmental stressors. However, research on the LecRLK gene family of conifers has seldom been reported.

Methods

Putative LecRLK genes were identified in the transcriptome of Taxodium 'Zhongshanshan'. The classification, domain structures, subcellular localization prediction, and expression patterns of LecRLK genes, as well as co-expressed genes, were analyzed using bioinformatics methods. Fifteen representative genes were further selected for qRT-PCR analysis in six tissues and under five different environmental stressor conditions.

Results

In total, 297 LecRLK genes were identified, including 155 G-type, 140 L-type, and 2 C-type. According to the classification, G-type and L-type LecRLK genes both can be organized into seven groups. The domain architecture of G-type proteins were more complex compared with that of L- and C-type proteins. Conservative motifs were found in G-type and L-type diverse lectin domains. Prediction and transient expression experiments to determine subcellular localization showed that LecRLKs were mainly concentrated in the cell membrane system, and some members were located at multiple sites at the same time. RNA-seq-based transcriptomics analysis suggested functional redundancy and divergence within each group. Unigenes co-expressed with LecRLKs in the transcriptome were found to be enriched in pathways related to signal transduction and environmental adaptation. Furthermore, qRT-PCR analysis of representative genes showed evidence of functional divergence between different groups.

Conclusions

This is the first study to conduct an identification and expression analysis of the LecRLK gene family in Taxodium. These results provide a basis for future studies on the evolution and function of this important gene family in Taxodium.