Transfer learning improves predictions in lignin content of Chinese fir based on Raman spectra.

Lignin in biomass plays significant role in substitution of synthetic polymer and reduction of energy expenditure, and the lignin content was usually determined by wet chemical methods. However, the methods' heavy workload, low efficiency, huge consumption of chemicals and use of toxic reagents render them unsuitable for sustainable development and environmental protection. Chinese fir, a prevalent angiosperm tree, holds immense importance for various industries. Since our previous work found that Raman spectroscopy could accurately predict the lignin content in poplar, we propose that the lignin content of Chinese fir can be estimated by similar strategy. The results suggested that the peak at 2895 cm-1 is the optimal choice of internal standard peak and algorithm of XGBoost demonstrates the highest accuracy among all algorithms. Furthermore, transfer learning was successfully introduced to enhance the accuracy and robustness of the model. Ultimately, we report that a machine learning algorithm, combining transfer learning with XGBoost or LightGBM, offers an accurate, high-efficiency and environmental friendly method for predicting the lignin content of Chinese fir using Raman spectra.

Identification and characterization of the WOX Gene Family revealed two WUS Clade Members associated with embryo development in Cunninghamia lanceolata.

The WUSCHEL-related homeobox (WOX) gene family is vital for plant development and stress response. In this study, we conducted a comprehensive analysis of WOX genes in Cunninghamia lanceolata (C. lanceolata) and subsequently explored the potential roles of two ClWOX genes within the WUS clade. In total, six ClWOX genes were identified through a full-length transcriptome analysis. These genes, exhibiting conserved structural and functional motifs, were assigned to the ancient clade and Modern/WUS clade, respectively, through a phylogenetic analysis. Our expression analysis indicated that these ClWOX genes were highly expressed in the middle and late developmental stages of zygotic embryos in C. lanceolata. Moreover, only ClWOX5 and ClWOX6 within the Modern/WUS clade exhibited transcriptional activity, and their expressions were also induced in response to auxin and wounding. Overexpression of ClWOX5 and ClWOX6 in Arabidopsis caused a partially sterile phenotype, resulting in a very low seed setting rate. Transcriptomic analysis revealed that expressions of many embryo-defective (EMB) genes, phytohormone-related genes, and transcription factors (TFs) were dramatically altered in ClWOX5 and ClWOX6 transgenic plants, which suggested that ClWOX5 and ClWOX6 may play specific important roles in embryo development via complex gene networks. In addition, overexpression of ClWOX5 and ClWOX6 in leaf segments promoted shoot regeneration in tobacco, indicating that ClWOX5 and ClWOX6 can promote plant regeneration and could be used to improve genetic transformation. In conclusion, these results help to elucidate the function of the WOX gene and provide a valuable basis for future studies of the developmental regulation and applications of WOX genes in C. lanceolata.

The transcriptional landscape of Populus pattern/effector-triggered immunity and how PagWRKY18 involved in it.

Plants trigger a robust immune response by activating massive transcriptome reprogramming through crosstalk between PTI and ETI. However, how PTI and ETI contribute to the quantitative or/and qualitative output of immunity and how they work together when both are being activated were unclear. In this study, we performed a comprehensive overview of pathogen-triggered transcriptomic reprogramming by analyzing temporal changes in the transcriptome up to 144 h after Colletotrichum gloeosporioides inoculated in Populus. Moreover, we constructed a hierarchical gene regulatory network of PagWRKY18 and its potential target genes to explore the underlying regulatory mechanisms of PagWRKY18 that are not yet clear. Interestingly, we confirmed that PagWRKY18 protein can directly bind the W-box elements in the promoter of a transmembrane leucine-rich repeat receptor-like kinase, PagSOBIR1 gene, to trigger PTI. At the same time, PagWRKY18 functions in disease tolerance by modulation of ROS homeostasis and induction of cell death via directly targeting PagGSTU7 and PagPR4 respectively. Furthermore, PagPR4 can interact with PagWRKY18 to inhibit the expression of PagPR4 genes, forming a negative feedback loop. Taken together, these results suggest that PagWRKY18 may be involved in regulating crosstalk between PTI and ETI to activate a robust immune response and maintain intracellular homeostasis.

PPGR: a comprehensive perennial plant genomes and regulation database.

Perennial woody plants hold vital ecological significance, distinguished by their unique traits. While significant progress has been made in their genomic and functional studies, a major challenge persists: the absence of a comprehensive reference platform for collection, integration and in-depth analysis of the vast amount of data. Here, we present PPGR (Resource for Perennial Plant Genomes and Regulation; https://ngdc.cncb.ac.cn/ppgr/) to address this critical gap, by collecting, integrating, analyzing and visualizing genomic, gene regulation and functional data of perennial plants. PPGR currently includes 60 species, 847 million protein-protein/TF (transcription factor)-target interactions, 9016 transcriptome samples under various environmental conditions and genetic backgrounds. Noteworthy is the focus on genes that regulate wood production, seasonal dormancy, terpene biosynthesis and leaf senescence representing a wealth of information derived from experimental data, literature mining, public databases and genomic predictions. Furthermore, PPGR incorporates a range of multi-omics search and analysis tools to facilitate browsing and application of these extensive datasets. PPGR represents a comprehensive and high-quality resource for perennial plants, substantiated by an illustrative case study that demonstrates its capacity in unraveling gene functions and shedding light on potential regulatory processes.

AC-UNet: an improved UNet-based method for stem and leaf segmentation in Betula luminifera.

Plant phenotypic traits play an important role in understanding plant growth dynamics and complex genetic traits. In phenotyping, the segmentation of plant organs, such as leaves and stems, helps in automatically monitoring growth and improving screening efficiency for large-scale genetic breeding. In this paper, we propose an AC-UNet stem and leaf segmentation algorithm based on an improved UNet. This algorithm aims to address the issues of feature edge information loss and sample breakage in the segmentation of plant organs, specifically in Betula luminifera. The method replaces the backbone feature extraction network of UNet with VGG16 to reduce the redundancy of network information. It adds a multi-scale mechanism in the splicing part, an optimized hollow space pyramid pooling module, and a cross-attention mechanism in the expanding network part at the output end to obtain deeper feature information. Additionally, Dice_Boundary is introduced as a loss function in the back-end of the algorithm to circumvent the sample distribution imbalance problem. The PSPNet model achieves mIoU of 58.76%, mPA of 73.24%, and Precision of 66.90%, the DeepLabV3 model achieves mIoU of 82.13%, mPA of 91.47%, and Precision of 87.73%, on the data set. The traditional UNet model achieves mIoU of 84.45%, mPA of 91.11%, and Precision of 90.63%, and the Swin-UNet model achieves . The mIoU is 79.02%, mPA is 85.99%, and Precision is 88.73%. The AC-UNet proposed in this article achieved excellent performance on the Swin-UNet dataset, with mIoU, mPA, and Precision of 87.50%, 92.71%, and 93.69% respectively, which are better than the selected PSPNet, DeepLabV3, traditional UNet, and Swin-UNet. Commonly used semantic segmentation algorithms. Experiments show that the algorithm in this paper can not only achieve efficient segmentation of the stem and leaves of Betula luminifera but also outperforms the existing state-of-the-art algorithms in terms of both speed. This can provide more accurate auxiliary support for the subsequent acquisition of plant phenotypic traits.

Orthogroup and phylotranscriptomic analyses identify transcription factors involved in the plant cold response: A case study of Arabidopsis BBX29.

C-repeat binding factors (CBFs) are well-known transcription factors (TFs) that regulate plant cold acclimation. RNA sequencing (RNA-seq) data from diverse plant species provide opportunities to identify other TFs involved in the cold response. However, this task is challenging because gene gain and loss has led to an intertwined community of co-orthologs and in-paralogs between and within species. Using orthogroup (closely related homologs) analysis, we identified 10,549 orthogroups in five representative eudicots. A phylotranscriptomic analysis of cold-treated seedlings from eudicots identified 35 high-confidence conserved cold-responsive transcription factor orthogroups (CoCoFos). These 35 CoCoFos included the well-known cold-responsive regulators CBFs, HSFC1, ZAT6/10, and CZF1 among others. We used Arabidopsis BBX29 for experimental validation. Expression and genetic analyses showed that cold-induction of BBX29 is CBF- and abscisic acid-independent, and BBX29 is a negative regulator of cold tolerance. Integrative RNA-seq and Cleavage Under Targets and Tagmentation followed by sequencing analyses revealed that BBX29 represses a set of cold-induced TFs (ZAT12, PRR9, RVE1, MYB96, etc.). Altogether, our analysis yielded a library of eudicot CoCoFos and demonstrated that BBX29 is a negative regulator of cold tolerance in Arabidopsis.

A CNN-LSTM-att hybrid model for classification and evaluation of growth status under drought and heat stress in chinese fir (Cunninghamia lanceolata).

BACKGROUND: Cunninghamia lanceolata (Chinese fir), is one of the most important timber trees in China. With the global warming, to develop new resistant varieties to drought or heat stress has become an essential task for breeders of Chinese fir. However, classification and evaluation of growth status of Chinese fir under drought or heat stress are still labor-intensive and time-consuming. RESULTS: In this study, we proposed a CNN-LSTM-att hybrid model for classification of growth status of Chinese fir seedlings under drought and heat stress, respectively. Two RGB image datasets of Chinese fir seedling under drought and heat stress were generated for the first time, and utilized in this study. By comparing four base CNN models with LSTM, the Resnet50-LSTM was identified as the best model in classification of growth status, and LSTM would dramatically improve the classification performance. Moreover, attention mechanism further enhanced performance of Resnet50-LSTM, which was verified by Grad-CAM. By applying the established Resnet50-LSTM-att model, the accuracy rate and recall rate of classification was up to 96.91% and 96.79% for dataset of heat stress, and 96.05% and 95.88% for dataset of drought, respectively. Accordingly, the R2 value and RMSE value for evaluation on growth status under heat stress were 0.957 and 0.067, respectively. And, the R2 value and RMSE value for evaluation on growth status under drought were 0.944 and 0.076, respectively. CONCLUSION: In summary, our proposed model provides an important tool for stress phenotyping in Chinese fir, which will be a great help for selection and breeding new resistant varieties in future.

Melatonin-mediated development and abiotic stress tolerance in plants.

Melatonin is a multifunctional molecule that has been widely discovered in most plants. An increasing number of studies have shown that melatonin plays essential roles in plant growth and stress tolerance. It has been extensively applied to alleviate the harmful effects of abiotic stresses. In view of its role in regulating aspects of plant growth and development, we ponder and summarize the scientific discoveries about seed germination, root development, flowering, fruit maturation, and senescence. Under abiotic and biotic stresses, melatonin brings together many pathways to increase access to treatments for the symptoms of plants and to counteract the negative effects. It has the capacity to tackle regulation of the redox, plant hormone networks, and endogenous melatonin. Furthermore, the expression levels of several genes and the contents of diverse secondary metabolites, such as polyphenols, terpenoids, and alkaloids, were significantly altered. In this review, we intend to examine the actions of melatonin in plants from a broader perspective, explore the range of its physiological functions, and analyze the relationship between melatonin and other metabolites and metabolic pathways.

The complete chloroplast genome of Rhododendron huadingense (Ericaceae).

Rhododendron huadingense is an important horticultural plant that belongs to the Ericaceae family. In this study, the chloroplast genome sequence of R. huadingense is reported. The chloroplast genome of R. huadingense was 198,952 bp in length and had an angiosperm-typical quadripartite structure with a large single-copy (LSC) region of 108,557 bp, a small single-copy (SSC) region of 53 bp, and two inverted repeat regions (IRs) of 45,171 bp. One hundred and thirteen unique genes including 79 protein-coding genes, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes were identified in the chloroplast genome. Further phylogenetic analysis revealed a close relationship between R. huadingense and R. molle. The complete chloroplast genome of R. huadingense provides valuable genetic information for the phylogeny, varieties breeding and sustainable utilization of this species.

Genome-wide identification and characterization of OVATE family proteins in Betula luminifera reveals involvement of BlOFP3 and BlOFP5 genes in leaf development.

Ovate family proteins (OFP) are plant-specific transcription factors involved in regulating morphologies of the lateral organs, plant growth and development. However, the functional roles of OFP genes in Betula luminifera, an important timber tree species, are not well studied. In this study, we identified 20 BlOFP genes and analyzed their phylogenetic relationship, gene structure, conserved motifs, and cis-elements. Further, expression analysis indicates that BlOFP genes were up-regulated in leaves on the one-year-old branch compared to leaves on the current-year branch and bract, except BlOFP7, BlOFP11, BlOFP14 and BlOFP12. The overexpression of BlOFP3 and BlOFP5 in Arabidopsis thaliana not only resulted in a slower growth rate but also produced sawtooth shape, flatter and darker green rosette leaves. Further investigation showed that the leaf thickness of the transgenic plants was more than double that of the wild type, which was caused by the increasement in the number and size of palisade tissue cells. Furthermore, the expression analysis also indicated that the expressions of several genes related to leaf development were significantly changed in the transgene plants. These results suggested the significant roles of BlOFP3 and BlOFP5 in leaf development. Moreover, protein-protein interaction studies showed that BlOFP3 interacts with BlKNAT5, and BlOFP5 interacts with BlKNAT5, BlBLH6 and BlBLH7. In conclusion, our study demonstrates that BlOFP3 and BlOFP5 were involved in leaf shape and thickness regulation by forming a complex with BlKNAT5, BlBLH6 and BlBLH7. In addition, our study serves as a guide for future functional genomic studies of OFP genes of the B. luminifera.

Comparative physiological analyses and the genetic basis reveal heat stress responses mechanism among different Betula luminifera populations.

Betula luminifera is a subtropical fast-growing timber species with high economic value. However, along with global warming, heat stress become one of the main environmental variables that limit the productivity of B. luminifera, and the response of diverse geographic populations to high temperatures is still unclear. In order to offer a comprehensive understanding of the behavior of B. luminifera under heat stress, the physiological responses of six B. luminifera populations (across the core distribution area) were described in this work in an integrated viewpoint. The results showed that a multi-level physiological regulatory network may exist in B. luminifera, the first response was the activity of resistant enzymes [e.g., peroxidase (POD)] at a preliminary stage of 2 h heat stress, and then the proline (osmoregulation substance) content began to increase after 24 h of continuous high-temperature treatment. In addition, photosynthesis was stronlgly affected by heat stress, and the net photosynthetic rate (Pn ) showed a downward trend under heat treatment in all six B. luminifera populations. Interestingly, although the physiological change patterns of the six B. luminifera populations were relatively consistent for the same parameter, there were obvious differences among different populations. Comprehensive analysis revealed that the physiological response of Rongshui (RS) was the most stable, and this was the representative B. luminifera population. Illumina RNA-seq analysis was applied to reveal the specific biological process of B. luminifera under heat stress using the RS population, and a total of 116,484 unigenes were obtained. The differentially expressed genes (DEGs) between different time periods under heat stress were enriched in 34 KEGG pathways, and the limonene and pinene degradation pathway was commonly enriched in all pairwise comparisons. Moreover, transcription factors including bHLH (basic helix-loop-helix), MYB, WRKY, and NAC (NAM, ATAF1/2, and CUC2) were identified. In this study, the physiological response and tolerance mechanisms of B. luminifera under high temperature stress were revealed, which can conducive to the basis of B. luminifera selection and resistance assessment for cultivation and breeding.

Identification of Nitrogen Starvation-Responsive miRNAs to Reveal the miRNA-Mediated Regulatory Network in Betula luminifera.

Because of the immobility, plants encounter a series of stresses, such as varied nutrient concentrations in soil, which regulate plant growth, development, and phase transitions. Nitrogen (N) is one of the most limiting factors for plants, which was exemplified by the fact that low nitrogen (LN) has a great adverse effect on plant growth and development. In the present study, we explored the potential role of microRNAs (miRNAs) in response to LN stress in Betula luminifera. We identified 198 miRNAs using sRNA sequencing, including 155 known and 43 novel miRNAs. Among them, 98 known miRNAs and 31 novel miRNAs were differentially expressed after 0.5 h or 24 h of LN stress. Based on degradome data, 122 differential expressed miRNAs (DEmiRNAs) including 102 known miRNAs and 20 novel miRNAs targeted 203 genes, comprising 321 miRNA-target pairs. A big proportion of target genes were transcription factors and functional proteins, and most of the Gene Ontology terms were enriched in biological processes; moreover, one Kyoto Encyclopedia of Genes and Genomes term "ascorbate and aldarate metabolism" was significantly enriched. The expression patterns of six miRNAs and their corresponding target genes under LN stress were monitored. According to the potential function for targets of DEmiRNAs, a proposed regulatory network mediated by miRNA-target pairs under LN stress in B. luminifera was constructed. Taken together, these findings provide useful information to elucidate miRNA functions and establish a framework for exploring N signaling networks mediated by miRNAs in B. luminifera. It may provide new insights into the genetic engineering of the high use efficiency of N in forestry trees.

The chromosome-scale genome of Phoebe bournei reveals contrasting fates of terpene synthase (TPS)-a and TPS-b subfamilies.

Terpenoids, including aromatic volatile monoterpenoids and sesquiterpenoids, function in defense against pathogens and herbivores. Phoebe trees are remarkable for their scented wood and decay resistance. Unlike other Lauraceae species investigated to date, Phoebe species predominantly accumulate sesquiterpenoids instead of monoterpenoids. Limited genomic data restrict the elucidation of terpenoid variation and functions. Here, we present a chromosome-scale genome assembly of a Lauraceae tree, Phoebe bournei, and identify 72 full-length terpene synthase (TPS) genes. Genome-level comparison shows pervasive lineage-specific duplication and contraction of TPS subfamilies, which have contributed to the extreme terpenoid variation within Lauraceae species. Although the TPS-a and TPS-b subfamilies were both expanded via tandem duplication in P. bournei, more TPS-a copies were retained and constitutively expressed, whereas more TPS-b copies were lost. The TPS-a genes on chromosome 8 functionally diverged to synthesize eight highly accumulated sesquiterpenes in P. bournei. The essential oil of P. bournei and its main component, ß-caryophyllene, exhibited antifungal activities against the three most widespread canker pathogens of trees. The TPS-a and TPS-b subfamilies have experienced contrasting fates over the evolution of P. bournei. The abundant sesquiterpenoids produced by TPS-a proteins contribute to the excellent pathogen resistance of P. bournei trees. Overall, this study sheds light on the evolution and adaptation of terpenoids in Lauraceae and provides valuable resources for boosting plant immunity against pathogens in various trees and crops.

Full-Length Transcriptome Sequencing and Comparative Transcriptomic Analyses Provide Comprehensive Insight Into Molecular Mechanisms of Cellulose and Lignin Biosynthesis in Cunninghamia lanceolata.

Cunninghamia lanceolata is an essential timber species that provide 20%-30% raw materials for China's timber industry. Although a few transcriptomes have been published in C. lanceolata, full-length mRNA transcripts and regulatory mechanisms behind the cellulose and lignin biosynthesis have not been thoroughly investigated. Here, PacBio Iso-seq and RNA-seq analyses were adapted to identify the full-length and differentially expressed transcripts along a developmental gradient from apex to base of C. lanceolata shoots. A total of 48,846 high-quality full-length transcripts were obtained, of which 88.0% are completed transcriptome based on benchmarking universal single-copy orthologs (BUSCO) assessment. Along stem developmental gradient, 18,714 differentially expressed genes (DEGs) were detected. Further, 28 and 125 DEGs were identified as enzyme-coding genes of cellulose and lignin biosynthesis, respectively. Moreover, 57 transcription factors (TFs), including MYB and NAC, were identified to be involved in the regulatory network of cellulose and lignin biosynthesis through weighted gene co-expression network analysis (WGCNA). These TFs are composed of a comparable regulatory network of secondary cell wall formation in angiosperms, revealing a similar mechanism may exist in gymnosperms. Further, through qRT-PCR, we also investigated eight specific TFs involved in compression wood formation. Our findings provide a comprehensive and valuable source for molecular genetics breeding of C. lanceolata and will be beneficial for molecular-assisted selection.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir).

BACKGROUND: R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line. RESULTS: The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms. CONCLUSIONS: This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

Nitrogen addition affects eco-physiological interactions between two tree species dominating in subtropical forests.

Nitrogen (N) deposition affects plant growth and interspecific interaction. This study aimed to explore the effect of N deposition on the growth and eco-physiological interactions between two tree species dominating in subtropical forests. A greenhouse experiment was conducted for 6 months in which the conifer Cunninghamia lanceolata and the broadleaved Phoebe chekiangensis were grown in monocultures and in a mixture under two levels of N addition: 0 and 45 kg ha-1 yr-1. The plant growth, root architecture, biomass distribution, element contents in plants and soil, and photosynthetic physiology were determined. The height and crown width of both seedlings tended to be higher in the mixture than in the monoculture when grown without N addition. P. chekiangensis was superior to C. lanceolata in resource acquisition and showed a greater net photosynthetic rate, plant height, crown width, total biomass, and belowground biomass distribution. In the mixture, N addition increased the net photosynthetic rate and decreased the height, ground diameter, and crown width of both species. Belowground biomass distribution was decreased in C. lanceolata but increased in P. chekiangensis under N addition. The P contents in both seedlings were higher in the mixture than in monocultures. Results showed N addition aggravated the competition and weakened the growth of both species in the mixture, largely determined by the competition for resources through the changing root architecture and biomass allocation. Our results provide new insights into the mechanisms of interspecific interaction in response to increasing N deposition in silvicultural practice.

Colchicine-Induced Polyploidy in Rhododendron fortunei Lindl.

Polyploidy in Rhododendron fortunei has great potential to improve its horticultural and commercial value, and to also meet market demands. In this study, a feasible method for polyploid induction in R. fortunei via colchicine treatment was established, and the obtained polyploid plants were identified and characterized. As a result, the stem bases of tissue-cultured plantlets treated with 0.1% colchicine for 24 h showed the highest polyploid induction with a rate of 36.67%. By flow cytometric analysis, 69 tetraploids and 29 octoploids were identified in the regenerated plants that were examined. Phenotypic analysis indicated that the leaves of tetraploid and octoploid plants were smaller, rounder and thicker with more abundant and longer epidermal hairs than those of diploids. Furthermore, the stomata of polyploids were larger and sparser than those of diploids. An increase in chlorophyll content was also detected in polyploids, which resulted in darker green leaves. In conclusion, our study established an effective method to induce polyploidy in R. fortunei, which could be used to develop new genetic resources for breeding R. fortunei and other Rhododendron species in the future.

Genome survey of Chinese fir (Cunninghamia lanceolata): Identification of genomic SSRs and demonstration of their utility in genetic diversity analysis.

Chinese fir (Cunninghamia lanceolata) is an important coniferous species that accounts for 20-30% of the total commercial timber production in China. Though traditional breeding of Chinese fir has achieved remarkable success, molecular-assisted breeding has made little progress due to limited availability of genomic information. In this study, a survey of Chinese fir genome was performed using the Illumina HiSeq Xten sequencing platform. K-mer analysis indicated that Chinese fir has a large genome of approximately 11.6 Gb with 74.89% repetitive elements and is highly heterozygous. Meanwhile, its genome size was estimated to be 13.2 Gb using flow cytometry. A total of 778.02 Gb clean reads were assembled into 10,982,272 scaffolds with an N50 of 1.57 kb. In total, 362,193 SSR loci were detected with a frequency of 13.18 kb. Dinucleotide repeats were the most abundant (up to 73.6% of the total SSRs), followed by trinucleotide and tetranucleotide repeats. Forty-six polymorphic pairs were developed, and 298 alleles were successfully amplified from 199 Chinese fir clones. The average PIC value was 0.53, indicating that the identified genomic SSR (gSSR) markers have a high degree of polymorphism. In addition, these breeding resources were divided into three groups, and a limited gene flow existed among these inferred groups.

RNA-Seq analysis of differential gene expression in Betula luminifera xylem during the early stages of tension wood formation.

BACKGROUND: Betula luminifera H. Winkler, which is widely distributed in southern China, is an economically important broadleaf tree species. However, little genomic information of B. luminifera is available, and little is known about the molecular mechanisms of wood formation in this species. Meanwhile, few efforts have focused on investigating the early transcriptional changes during tension wood formation in woody plants. RESULTS: A reference transcriptome dataset was first generated containing 45,700 Unigenes, and 35,135 (76.9%) Unigenes were annotated by a BLAST similarity search against four public databases. Then, based on an anatomical investigation, the global gene expression changes during the early stages of tension wood formation were analyzed. Gene expression profiling showed that a total of 13,273 Unigenes were differentially regulated during the early stages of tension wood formation. Most genes involved in cellulose and lignin biosynthesis were highlighted to reveal their biological importance in tension wood formation. In addition, the transcription levels of many genes involved in the auxin response pathway were significantly changed during the early stages of tension wood formation. Furthermore, 18 TFs co-expressed with key enzymes of cellulose synthesis were identified. CONCLUSIONS: Our results revealed the transcriptional changes associated with TW formation and identified potential key genes in the regulation of this process. These results will help to dissect the molecular mechanism of wood formation and provide key candidate genes for marker-assisted selection in B. luminifera.

Molecular Characterization of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family in Betula luminifera.

As a major family of plant-specific transcription factors, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes play vital regulatory roles in plant growth, development and stress responses. In this study, 18 SPL genes were identified and cloned from Betula luminifera. Two zinc finger-like structures and a nuclear location signal (NLS) segments were existed in the SBP domains of all BlSPLs. Phylogenetic analysis showed that these genes were clustered into nine groups (group I-IX). The intron/exon structure and motif composition were highly conserved within the same group. 12 of the 18 BlSPLs were experimentally verified as the targets of miR156, and two cleavage sites were detected in these miR156-targeted BlSPL genes. Many putative cis-elements, associated with light, stresses and phytohormones response, were identified in the promoter regions of BlSPLs, suggesting that BlSPL genes are probably involved in important physiological processes and developmental events. Tissue-specific expression analysis showed that miR156-targeted BlSPLs exhibited a more differential expression pattern, while most miR156-nontargeted BlSPLs tended to be constitutively expressed, suggesting the distinct roles of miR156-targeted and nontargeted BlSPLs in development and growth of B. luminifera. Further expression analysis revealed that miR156-targeted BlSPLs were dramatically up-regulated with age, whereas mature BlmiR156 level was apparently declined with age, indicating that miR156/SPL module plays important roles in vegetative phase change of B. luminifera. Moreover, yeast two-hybrid assay indicated that several miR156-targeted and nontargeted BlSPLs could interact with two DELLA proteins (BlRGA and BlRGL), which suggests that certain BlSPLs take part in the GA regulated processes through protein interaction with DELLA proteins. All these results provide an important basis for further exploring the biological functions of BlSPLs in B. luminifera.