Caffeic acid O-methyltransferase (COMT) is an important protein that participates in lignin synthesis and is associated with the ratio of G-/S-type lignin in plants. COMTs are associated with the wood properties of forest trees; however, little known about the COMT family in Catalpa bungei, a valuable timber tree species in China . We performed a comprehensive analysis of COMT genes in the C. bungei genome by describing the gene structure and phylogenetic relationships of each family member using bioinformatics-based methods. A total of 23 putative COMT genes were identified using the conserved domain sequences and amino acid sequences of COMTs from Arabidopsis thaliana and Populus trichocarpa as probes. Phylogenetic analysis showed that 23 CbuCOMTs can be divided into three groups based on their structural characteristics; five conserved domains were found in the COMT family. Promoter analysis indicated that the CbuCOMT promoters included various cis-acting elements related to growth and development. Real-time quantitative polymerase chain reaction (PCR) analysis showed differential expression among CbuCOMTs. CbuCOMT2, 7, 8, 9, 10, 12, 13, 14, 21, and 23 were mainly expressed in xylem. Only CbuCOMT23 was significantly downregulated in tension wood and upregulated in opposite wood compared to normal wood. Our study provides new information about the CbuCOMT gene family and will facilitate functional characterisation in further research.

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.

Clones of Catalpa bungei and Catalpa fargesii f. duclouxii were studied over several years in central China to explore genetic variation in growth traits and to identify clones of high wood yield and high stability. The genetic parameters for height, diameter at breast height (DBH), and stem volume of clones, were estimated. The effect of clone × year on the increment of stem volume in the two species was analyzed by genotype and genotype × environment (GGE) biplot methods. Significant differences in growth traits among clones and between species were found. The growth of C. bungei exceeded that of C. fargesii f. duclouxii after 4 years. Furthermore, from the 5th year, the repeatability and genetic variation coefficient (GCV) of the C. bungei clones were higher than those of the C. fargesii f. duclouxii clones in most cases. The phenotypic variation coefficient (PCV) of the C. fargesii f. duclouxii clones was significantly lower than that of the C. bungei clones. The repeatability of stem volume was intermediate or high in the two species. ANOVA revealed significant effects of the clone by year interaction in these two species. GGE biplot analysis revealed that wood yield and stability were largely independent in C. bungei; clones 22-03, 19-27, and 20-01 were the optimal clones in this species. In contrast, the optimal clones 63 and 128 of C. fargesii f. duclouxii combined the desired characteristics of high yield and high stability. In conclusion, our results indicated that the height and stem volume of C. bungei was under strong genetic control, whereas that of C. fargesii f. duclouxii was influenced by the environment more than by genetic effects. Genetic improvement by clone selection can be expected to be effective, as the repeatability of stem volume was high. Francis and Kannenberg’s method and GGE biplot analysis were used in combination to evaluate the clones. C. bungei clone 22-03 and C. fargesii f. duclouxii clones 63 and 128 were identified as the optimal clones, which exhibited both a high increment of stem volume and high stability.

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.

Sophora moorcroftiana is a perennial leguminous low shrub endemic to the Yarlung Zangbo River basin in Tibet with irreplaceable economic and ecological value. To determine the drivers of evolution in this species, 225 individuals belonging to 15 populations from different geographic locations were sampled, and population genetics was studied using high-throughput genotyping-by-sequencing (GBS). Based on genetic diversity analysis, phylogenetic analysis, principal component analysis, and structure analysis, 15 natural populations were clustered into the following five subgroups: subgroup I (Shigatse subgroup) was located in the upper reaches of the Yarlung Zangbo River with a relatively high level of population genetic variation (means for PIC, Shannon and PI were 0.173, 0.326 and 0.0000305, respectively), and gene flow within the subgroup was also high (mean value for Nm was 4.67). Subgroup II (including Pop 7 and Pop 8; means for PIC, Shannon and PI were 0.182, 0.345 and 0.0000321, respectively), located in the middle reaches of the Yarlung Zangbo River had relatively high levels of gene flow with the populations distributed in the upper and lower reaches. The Nm between subgroup II with subgroups I and III was 3.271 and 2.894, respectively. Considering all the genetic diversity indices Pop 8 had relatively high genetic diversity. Subgroup III (the remaining mixed subgroup of Lhasa and Shannan) was located in the middle reaches of the Yarlung Zangbo River and the means for PIC, Shannon and PI were 0.172, 0.324 and 0.0000303, respectively. Subgroup IV (Nyingchi subgroup), located in the lower reaches of the Yarlung Zangbo River basin, showed a further genetic distance from the other subgroups and the means for PIC, Shannon and PI were 0.147, 0.277 and 0.0000263, respectively. Subgroup V (Nyingchi Gongbu Jiangda subgroup), located in the upper reaches of the Niyang River, had the lowest level of genetic variation (means for PIC, Shannon and PI were 0.106, 0.198 and 0.0000187, respectively) and gene flow with other populations (mean value for Nm was 0.42). According to the comprehensive analysis, the S. moorcroftiana populations generally expanded from upstream to downstream and displayed a high level of genetic differentiation in the populations in the upper and lower reaches. There were high levels of gene exchange between the central populations with upstream and downstream populations, and wind-induced seed dispersal was an important factor in the formation of this gene exchange mode.

Long non-coding RNAs (lncRNAs) have crucial roles in various biological regulatory processes. However, the study of lncRNAs is limited in woody plants. Catalpa bungei is a valuable ornamental tree with a long cultivation history in China, and a deeper understanding of the floral transition mechanism in C. bungei would be interesting from both economic and scientific perspectives. In this study, we categorized C. bungei buds from early flowering (EF) and normal flowering (NF) varieties into three consecutive developmental stages. These buds were used to systematically study lncRNAs during floral transition using high-throughput sequencing to identify molecular regulatory networks. Quantitative real-time PCR was performed to study RNA expression changes in different stages. In total, 12,532 lncRNAs and 26,936 messenger RNAs (mRNAs) were detected. Moreover, 680 differentially expressed genes and 817 differentially expressed lncRNAs were detected during the initiation of floral transition. The results highlight the mRNAs and lncRNAs that may be involved in floral transition, as well as the many lncRNAs serving as microRNA precursors. We predicted the functions of lncRNAs by analysing the relationships between lncRNAs and mRNAs. Seven lncRNA-mRNA interaction pairs may participate in floral transition. This study is the first to identify lncRNAs and their potential functions in floral transition, providing a starting point for detailed determination of the functions of lncRNAs in C. bungei.