Genetic Variation of Growth Traits and Genotype-by-Environment Interactions in Clones of Catalpa bungei and Catalpa fargesii f. duclouxii

Effects of nitrate- and ammonium- nitrogen on anatomical and physiological responses of Catalpa bungei under full and partial root-zone drought

Catalpa bungei is a precious timber species distributed in North China where drought often occurs. To clarify adaptive responses of C. bungei to partial- and full- root-zone drought under the influence of nitrogen forms, a two-factor experiment was conducted in which well-watered (WW), partial root-zone drought in horizontal direction (H-PRD) and in vertical direction (V-PRD), and full root-zone drought (FRD) were combined with nitrate-nitrogen (NN) and ammonium-nitrogen (AN) treatments. C. bungei responded to FRD by sharply closing stomata, decreasing gas exchange rate and increasing leaf instantaneous water use efficiency (WUEi). Under FRD condition, the growth of seedlings was severely inhibited and the effect of N forms was covered up by the drastic drought effect. In comparison, stomata conductance and gas exchanges were moderately inhibited by PRDs. WUEi in V-PRD treatment was superior to H-PRD due to the active stomata regulation resulting from a higher ABA level and active transcription of genes in abscisic acid (ABA) signaling pathway under V-PRD. Under both PRDs and FRD, nitrate benefited antioxidant defense, stomata regulation and leaf WUEi. Under V-PRD, WUEi in nitrate treatment was superior to that in ammonium treatment due to active stomata regulation by signaling network of nitric oxide (NO), Ca2+ and ABA. Under FRD, WUEi was higher in nitrate treatment due to the favoring photosynthetic efficiency resulting from active NO signal and antioxidant defense. The interactive effect of water and N forms was significant on wood xylem development. Superoxide dismutase (SOD) and catalase (CAT) largely contributes to stress tolerance and xylem development.

QTL mapping of drought-related traits in the hybrids of Populus deltoides 'Danhong'×Populus simonii 'Tongliao1'

BACKGROUND

Poplar trees provide a large amount of wood material, but many parts of the world are arid or semi-arid areas because of insufficient annual precipitation, which seriously affects the growth of poplar trees. Populus simonii 'Tongliao1' shows strong tolerance to stress environments, and Populus deltoides 'Danhong' shows a stronger growth rate in a suitable environment. To identify drought tolerance-related QTLs and genes, an F₁ population derived from the cross between the 'Danhong' and 'Tongliao 1' Populus was assessed under drought stress.

RESULTS

We measured drought-related traits such as the relative height growth, relative diameter growth, leaf senescence number, specific leaf area, and leaf relative water content in the population under control and drought environments. The results showed that drought stress reduced the plant height relative growth, ground diameter relative growth, specific leaf area and leaf relative water content and increased the number of leaf drops. A total of 208 QTLs were identified by QTL mapping analysis, and they consisted of 92, 63 and 53 QTLs under control, drought stress treatment and drought index conditions, respectively. A molecular identification marker for drought tolerance, np2841, which was associated with a QTL (qDLRWC-LG10-1) for relative leaf water content, was initially developed. We mined 187 candidate genes for QTL regions of five traits under a drought environment. The reference genome annotation for Populus trichocarpa and a homologous gene analysis of Arabidopsis thaliana identified two candidate genes, Potri.003G171300 and Potri.012G123900, with significant functions in response to drought stress. We identified five key regulatory genes (Potri.006G273500, Potri.007G111500, Potri.007G111600, Potri.007G111700, and Potri.007G111800) related to drought tolerance through the poplar coexpression network.

CONCLUSION

In this study, our results indicate that the QTLs can effectively enhance the drought tolerance of poplar. It is a step closer towards unravelling the genetic basis of poplar drought tolerance-related traits, and to providing validated candidate genes and molecular markers for future genetic improvement.

The Current and Future Potential Geographical Distribution and Evolution Process of Catalpa bungei in China

Catalpa bungei C. A. Mey. (C. bungei) is one of the recommended native species for ecological management in China. It is a fast-growing tree of high economic and ecological importance, but its rare resources, caused by anthropogenic destruction and local climatic degradation, have not satisfied the requirements. It has been widely recommended for large-scale afforestation of ecological management and gradually increasing in recent years, but the impact mechanism of climate change on its growth has not been studied yet. Studying the response of species to climate change is an important part of national afforestation planning. Based on combinations of climate, topography, soil variables, and the multiple model ensemble (MME) of CMIP6, this study explored the relationship between C. bungei and climate change, then constructed Maxent to predict its potential distribution under SSP126 and SSP585 and analyzed its dominant environmental factors. The results showed that C. bungei is widely distributed in Henan, Hebei, Hubei, Anhui, Jiangsu, and Shaanxi provinces and others where it covers an area of 2.96 × 106 km2. Under SSP126 and SSP585, its overall habitat area will increase by more than 14.2% in 2080–2100, which mainly indicates the transformation of unsuitable areas into low suitable areas. The center of its distribution will migrate to the north with a longer distance under SSP585 than that under SSP126, and it will transfer from the junction of Shaanxi and Hubei province to the north of Shaanxi province under SSP585 by 2100. In that case, C. bungei shows a large-area degradation trend in the south of the Yangtze River Basin but better suitability in the north of the Yellow River Basin, such as the Northeast Plain, the Tianshan Mountains, the Loess Plateau, and others. Temperature factors have the greatest impact on the distribution of C. bungei. It is mainly affected by the mean temperature of the coldest quarter, followed by precipitation of the wettest month, mean diurnal range, and precipitation of the coldest quarter. Our results hence demonstrate that the increase of the mean temperature of the coldest quarter becomes the main reason for its degradation, which simultaneously means a larger habitat boundary in Northeast China. The findings provide scientific evidence for the ecological restoration and sustainable development of C. bungei in China.

A Bivariate Mapping Model Identifies Major Covariation QTLs for Biomass Allocation Between Leaf and Stem Growth of Catalpa bungei

Biomass allocation plays a critical role in plant morphological formation and phenotypic plasticity, which greatly impact plant adaptability and competitiveness. While empirical studies on plant biomass allocation have focused on molecular biology and ecology approaches, detailed insight into the genetic basis of biomass allocation between leaf and stem growth is still lacking. Herein, we constructed a bivariate mapping model to identify covariation QTLs governing carbon (C) allocation between the leaves and stem as well as the covariation of traits within and between organs in a full-sib mapping population of C. bungei. A total of 123 covQTLs were detected for 23 trait pairs, including six leaf traits (leaf length, width, area, perimeter, length/width ratio and petiole length) and five stem traits (height, diameter at breast height, wood density, stemwood volume and stemwood biomass). The candidate genes were further identified in tissue-specific gene expression data, which provided insights into the genetic architecture underlying C allocation for traits or organs. The key QTLs related to growth and biomass allocation, which included UVH1, CLPT2, GAD/SPL, COG1 and MTERF4, were characterised and verified via gene function annotation and expression profiling. The integration of a bivariate Quantitative trait locus mapping model and gene expression profiling will enable the elucidation of genetic architecture underlying biomass allocation and covariation growth, in turn providing a theoretical basis for forest molecular marker-assisted breeding with specific C allocation strategies for adaptation to heterogeneous environments.

Transcriptomic, Proteomic, and Metabolic Profiles of Catalpa bungei Tension Wood Reveal New Insight Into Lignin Biosynthesis Involving Transcription Factor Regulation

Lignin is a complex polymer in plant cell walls whose proportion is second only to that of cellulose and plays an important role in the mechanical properties of wood and stress resistance of plants. Here, we induced tension wood (TW) formation in Catalpa bungei by artificial bending and analyzed the lignin metabolism of the TW. LC-MS analysis showed that a significantly higher content of coniferyl aldehyde was observed in the TW cell wall than in the opposite wood (OW) and normal wood (NW) cell walls. TW had significantly lower contents of coniferyl alcohol than OW and NW. Raman spectroscopy results indicated that TW had lower total lignin than OW and NW. The transcription and translation levels of most of the differentially expressed genes (DEGs) involved in lignin monomer biosynthesis indicated upregulation in TW/OW and TW/NW. We found no significant difference in the transcription levels of three collision gases (CADs) between TW and OW or between NW, but their translation levels were significantly downregulated in TW, suggesting post-transcriptional control for CAD. We predicted and analyzed transcription factors that could target DEGs involved in lignin monomer biosynthesis in TW. Based on the analysis of the relationships of targeting and coexpression, we found that NAC (evm.model.group1.695) could potentially target 4CLs and CCoAOMT, that HD-Zip (evm.model.group7.1157) had potential targeting relationships with CCoAOMT, F5H, and CCR, and that their expression levels were significantly positive. It is speculated that the upregulation of NAC and HD-ZIP transcription factors activates the expression of downstream target genes, which leads to a significant increase in coniferyl aldehyde in TW. However, the decrease in total lignin in TW may be caused by the significant downregulation of CAD translation and the significant decrease in precursors (coniferyl alcohol). Whether the expression of CAD genes is regulated by post-transcriptional control and affects TW lignin metabolism needs further study.

Genetic Evaluation and Combined Selection for the Simultaneous Improvement of Growth and Wood Properties in Catalpa bungei Clones

Catalpa bungei is an important timber tree. Improvements in growth and wood quality are important goals of C. bungei breeding, and it is necessary to understand the genetic parameters of specific target traits and genetic correlation between growth traits and wood properties for tree breeding. In this study, the genetic parameters of height, diameter at breast height (DBH) and wood properties were estimated and genetic and phenotypic correlations between growth traits and wood properties were evaluated in C. bungei. Finally, different selection scenarios were used to evaluate and select optimal clones. The results showed that there were significant differences in growth and wood properties among clones. The wood hardness (0.66–0.79), basic density (0.89), air-dried density (0.89) and compression strength parallel to the grain of wood (CSP) (0.84) had high repeatability. The variance component proportions indicated that the variation in wood properties came mainly from different genotypes (clones) rather than from different individuals of the same clone. The DBH showed a significant negative genetic correlation with the hardness of radial section (HRS) (−643), basic density (−0.531) and air-dry density (−0.495). This unfavorable relationship makes it difficult to improve growth and wood quality simultaneously in C. bungei. We selected the optimal clones under different scenarios, and we obtained 7.75–9.06% genetic gains for growth in the scenario in which height and DBH were the target traits. Genetic gains of 7.43–14.94% were obtained for wood properties by selecting optimal clones in the scenario in which wood properties were the target traits. Approximately 5% and 4% genetic gains were obtained for growth and wood properties, respectively, for the combined selection. This study provides new insights into the genetic improvement of wood quality in C. bungei.

The complete chloroplast genome of Catalpa fargesii Bur. f. duclouxii (Dode) Gilmour

In this study, we reported and characterized the complete chloroplast genome sequence of Catalpa fargesii Bur. f. duclouxii (Dode) Gilmour. The chloroplast genome was determined to be 158,250 bp in length. It contained large single-copy (LSC) and small single-copy (SSC) regions of 84,929 bp and 12,663 bp, respectively, which were separated by a pair of 30,329 bp inverted repeat (IR) regions. The genome is predicted to contain 121 genes, including 78 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. The overall GC content of the genome is 38.1%. A phylogenetic tree reconstructed by 12 chloroplast genomes reveals that C. fargesii is mostly related to Catalpa. ovata and Catalpa. speciosa. This study identified the unique characteristics of the C. fargesii cp genome, which will provide a theoretical basis for species identification and biological research.

Genome-wide analysis of lncRNA and mRNA expression and endogenous hormone regulation during tension wood formation in Catalpa bungei

BACKGROUND

Phytohormones are the key factors regulating vascular development in plants, and they are also involved in tension wood (TW) formation. Although the theory of hormone distribution in TW formation is widely supported, the effects of endogenous hormones on TW formation have not yet been assessed. In this study, TW formation was induced in Catalpa bungei by artificial bending. The phytohormone content of TW, opposite wood (OW) and normal wood (NW) was determined using liquid chromatography-mass spectrometry (LC-MS), and transcriptome sequencing was performed. The hormone content and related gene expression data were comprehensively analyzed.

RESULTS

The results of analyses of the plant hormone contents indicated significantly higher levels of cis-zeatin (cZ), indoleacetic acid (IAA) and abscisic acid (ABA) in TW than in OW. Genes involved in the IAA and ABA synthesis pathways, such as ALDH (evm.

MODEL

group5.1511) and UGT (evm.

MODEL

scaffold36.20), were significantly upregulated in TW. and the expression levels of ARF (evm.

MODEL

group5.1332), A-ARR (evm.

MODEL

group0.1600), and TCH4 (evm.

MODEL

group2.745), which participate in IAA, cZ and Brassinolide (BR) signal transduction, were significantly increased in TW. In particular, ARF expression may be regulated by long noncoding RNAs (lncRNAs) and the HD-ZIP transcription factor ATHB-15.

CONCLUSIONS

We constructed a multiple hormone-mediated network of C. bungei TW formation based on hormone levels and transcriptional expression profiles were identified during TW formation.

Phenotypic Variability and Genetic Diversity in a Pinus koraiensis Clonal Trial in Northeastern China

Combining phenotypic and genetic characteristics in a genetic variation study is of paramount importance to effectively orient the selection of producers’ elite trees in a seed orchard. In total, 28 phenotypic characteristics and 16 microsatellite loci were used to analyze the clonal genetic variation, to characterize the genetic diversity, and to refine the genetic classifications of 110 Pinus koraiensis clones grown in the Naozhi orchard in northeastern China. All clones were significantly different in most traits. Most of the phenotypic characteristics showed great genetic variation among clones, while the genotypic differentiation was weak between the selection sites of clones. The SSR markers showed a relatively high level of genetic diversity (Na = 4.67 ± 0.43, Ne = 2.916 ± 0.18, I = 1.15 ± 0.07, Ho = 0.69 ± 0.04, He = 0.62 ± 0.02, and mean polymorphic information content (PIC) of 0.574), with higher heterozygosity as an indication of a lower probability of inbreeding events in the orchard. Despite weak correlation coefficients between dissimilarity matrices (r(A/B), range equal to 0.022, p-value < 0.001), the genetic and phenotypic classifications congruently subdivided all the clones into three major groups. The patterns of phenotypic trait variations and genetic diversity are valuable to effectively select materials in breeding programs of P. koraiensis.

Transcriptomics and Proteomics Reveal the Cellulose and Pectin Metabolic Processes in the Tension Wood (Non-G-Layer) of Catalpa bungei

: Catalpa bungei is an economically important tree with high-quality wood and highly valuable to the study of wood formation. In this work, the xylem microstructure of C. bungei tension wood (TW) was observed, and we performed transcriptomics, proteomics and Raman spectroscopy of TW, opposite wood (OW) and normal wood (NW). The results showed that there was no obvious gelatinous layer (G-layer) in the TW of C. bungei and that the secondary wall deposition in the TW was reduced compared with that in the OW and NW. We found that most of the differentially expressed mRNAs and proteins were involved in carbohydrate polysaccharide synthesis. Raman spectroscopy results indicated that the cellulose and pectin content and pectin methylation in the TW were lower than those in the OW and NW, and many genes and proteins involved in the metabolic pathways of cellulose and pectin, such as galacturonosyltransferase (GAUT), polygalacturonase (PG), endoglucanase (CLE) and β-glucosidase (BGLU) genes, were significantly upregulated in TW. In addition, we found that the MYB2 transcription factor may regulate the pectin degradation genes PG1 and PG3, and ARF, ERF, SBP and MYB1 may be the key transcription factors regulating the synthesis and decomposition of cellulose. In contrast to previous studies on TW with a G-layer, our results revealed a change in metabolism in TW without a G-layer, and we inferred that the change in the pectin type, esterification and cellulose characteristics in the TW of C. bungei may contribute to high tensile stress. These results will enrich the understanding of the mechanism of TW formation.

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.

Deciphering Genetic Architecture of Adventitious Root and Related Shoot Traits in Populus Using QTL Mapping and RNA-Seq Data

Understanding the genetic architecture of adventitious root and related shoot traits will facilitate the cultivation of superior genotypes. In this study, we measured 12 adventitious root and related shoot traits of 434 F₁ genotypes originating from Populus deltoides ‘Danhong’ × Populus simonii ‘Tongliao1’ and conducted an integrative analysis of quantitative trait locus (QTL) mapping and RNA-Seq data to dissect their genetic architecture and regulatory genes. Extensive segregation, high repeatability, and significant correlation relationship were detected for the investigated traits. A total of 150 QTLs were associated with adventitious root traits, explaining 3.1–6.1% of phenotypic variation (PVE); while 83 QTLs were associated with shoot traits, explaining 3.1–19.8% of PVE. Twenty-five QTL clusters and 40 QTL hotspots were identified for the investigated traits. Ten QTL clusters were overlapped in both adventitious root traits and related shoot traits. Transcriptome analysis identified 10,172 differentially expressed genes (DEGs) among two parents, three fine rooting and three poor-rooting genotypes, 143 of which were physically located within the QTL intervals. K-means cluster and weighted gene co-expression network analysis showed that PtAAAP19 (Potri.004G111400) encoding amino acid transport protein was tightly associated with adventitious roots and highly expressed in fine-rooting genotypes. Compare with ‘Danhong’, 153 bp deletion in the coding sequence of PtAAAP19 in ‘Tongliao1’ gave rise to lack one transmembrane domain, which might cause the variation of adventitious roots. Taken together, this study deciphered the genetic basis of adventitious root and related shoot traits and provided potential function genes for genetic improvement of poplar breeding.