The relationship between atmospheric particulate matter, leaf surface microstructure, and the phyllosphere microbial diversity of Ulmus L.

BACKGROUND: Plants can retain atmospheric particulate matter (PM) through their unique foliar microstructures, which has a profound impact on the phyllosphere microbial communities. Yet, the underlying mechanisms linking atmospheric particulate matter (PM) retention by foliar microstructures to variations in the phyllosphere microbial communities remain a mystery. In this study, we conducted a field experiment with ten Ulmus lines. A series of analytical techniques, including scanning electron microscopy, atomic force microscopy, and high-throughput amplicon sequencing, were applied to examine the relationship between foliar surface microstructures, PM retention, and phyllosphere microbial diversity of Ulmus L. RESULTS: We characterized the leaf microstructures across the ten Ulmus lines. Chun exhibited a highly undulated abaxial surface and dense stomatal distribution. Langya and Xingshan possessed dense abaxial trichomes, while Lieye, Zuiweng, and Daguo had sparsely distributed, short abaxial trichomes. Duomai, Qingyun, and Lang were characterized by sparse stomata and flat abaxial surfaces, whereas Jinye had sparsely distributed but extensive stomata. The mean leaf retention values for total suspended particulate (TSP), PM2.5, PM2.5-10, PM10-100, and PM> 100 were 135.76, 6.60, 20.10, 90.98, and 13.08 µg·cm- 2, respectively. Trichomes substantially contributed to PM2.5 retention, while larger undulations enhanced PM2.5-10 retention, as evidenced by positive correlations between PM2.5 and abaxial trichome density and between PM2.5-10 and the adaxial raw microroughness values. Phyllosphere microbial diversity patterns varied among lines, with bacteria dominated by Sediminibacterium and fungi by Mycosphaerella, Alternaria, and Cladosporium. Redundancy analysis confirmed that dense leaf trichomes facilitated the capture of PM2.5-associated fungi, while bacteria were less impacted by PM and struggled to adhere to leaf microstructures. Long and dense trichomes provided ideal microhabitats for retaining PM-borne microbes, as evidenced by positive feedback loops between PM2.5, trichome characteristics, and the relative abundances of microorganisms like Trichoderma and Aspergillus. CONCLUSIONS: Based on our findings, a three-factor network profile was constructed, which provides a foundation for further exploration into how different plants retain PM through foliar microstructures, thereby impacting phyllosphere microbial communities.

Nanoplastic toxicity induces metabolic shifts in Populus × euramericana cv. '74/76' revealed by multi-omics analysis.

There is increasing global concern regarding the pervasive issue of plastic pollution. We investigated the response of Populus × euramericana cv. '74/76' to nanoplastic toxicity via phenotypic, microanatomical, physiological, transcriptomic, and metabolomic approaches. Polystyrene nanoplastics (PS-NPs) were distributed throughout the test plants after the application of PS-NPs. Nanoplastics principally accumulated in the roots; minimal fractions were translocated to the leaves. In leaves, however, PS-NPs easily penetrated membranes and became concentrated in chloroplasts, causing thylakoid disintegration and chlorophyll degradation. Finally, oxidant damage from the influx of PS-NPs led to diminished photosynthesis, stunted growth, and etiolation and/or wilting. By integrating dual-omics data, we found that plants could counteract mild PS-NP-induced oxidative stress through the antioxidant enzyme system without initiating secondary metabolic defense mechanisms. In contrast, severe PS-NP treatments promoted a shift in metabolic pattern from primary metabolism to secondary metabolic defense mechanisms, an effect that was particularly pronounced during the upregulation of flavonoid biosynthesis. Our findings provide a useful framework from which to further clarify the roles of key biochemical pathways in plant responses to nanoplastic toxicity. Our work also supports the development of effective strategies to mitigate the environmental risks of nanoplastics by biologically immobilizing them in contaminated lands.

Genome-wide molecular evolution analysis of the GRF and GIF gene families in Plantae (Archaeplastida).

BACKGROUND: Plant growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) interact with each other and collectively have important regulatory roles in plant growth, development, and stress responses. Therefore, it is of great significance to explore the systematic evolution of GRF and GIF gene families. However, our knowledge and understanding of the role of GRF and GIF genes during plant evolution has been fragmentary. RESULTS: In this study, a large number of genomic and transcriptomic datasets of algae, mosses, ferns, gymnosperms and angiosperms were used to systematically analyze the evolution of GRF and GIF genes during the evolution of plants. The results showed that GRF gene first appeared in the charophyte Klebsormidium nitens, whereas the GIF genes originated relatively early, and these two gene families were mainly expanded by segmental duplication events after plant terrestrialization. During the process of evolution, the protein sequences and functions of GRF and GIF family genes are relatively conservative. As cooperative partner, GRF and GIF genes contain the similar types of cis-acting elements in their promoter regions, which enables them to have similar transcriptional response patterns, and both show higher levels of expression in reproductive organs and tissues and organs with strong capacity for cell division. Based on protein-protein interaction analysis and verification, we found that the GRF-GIF protein partnership began to be established in pteridophytes and is highly conserved across different terrestrial plants. CONCLUSIONS: These results provide a foundation for further exploration of the molecular evolution and biological functions of GRF and GIF genes.

GRP78 promotes bone metastasis of prostate cancer by regulating bone microenvironment through Sonic hedgehog signaling.

Bone metastasis is the leading cause of tumor-related deaths in patients with prostate cancer (PCa). The interactions between PCa and the bone microenvironment form a vicious cycle. However, the complex molecular mechanism by which PCa regulates the bone microenvironment remains unclear. To determine the role of glucose-regulated protein (GRP78) in bone metastasis and growth, we established intracardiac injection and tibial injection models, and performed their histological staining. To assess the effect of GRP78 on the differentiation of osteoblasts and osteoclasts, we performed cell co-culture, enzyme-linked immunosorbent assay, alizarin red staining, and tartrate-resistant acid phosphatase staining. We found that GRP78 is upregulated in PCa tissues and that its upregulation is associated with PCa progression in patients. Functional experiments showed that GRP78 overexpression in PCa cells considerably promotes bone metastasis and induces bone microstructure changes. Silencing GRP78 substantially inhibits the migration and invasion of PCa cells in vitro and bone metastasis and tumor growth in vivo. Mechanistically, GRP78 promotes the migration and invasion of PCa cells via the Sonic hedgehog (Shh) signaling pathway. Cell co-culture showed that GRP78 promotes the differentiation of osteoblasts and osteoclasts through Shh signaling. Our findings suggest that tumor-bone matrix interactions owing to GRP78-activated paracrine Shh signaling by PCa cells regulate the differentiation of osteoblasts and osteoclasts. This process promotes bone metastasis and the proliferation of PCa cells in the bone microenvironment. Targeting the GRP78/Shh axis can serve as a therapeutic strategy to prevent bone metastasis and improve the quality of life of patients with PCa.

Genome-wide identification of the Pyrus R2R3-MYB gene family and PhMYB62 regulation analysis in Pyrus hopeiensis flowers at low temperature.

The R2R3-MYB gene family play an important role in plant growth, development and stress responses. In this study, a total of 122 PcoR2R3-MYB genes were identified and grouped into 26 clades in pear. And these PcoMYBs were unevenly distributed among 17 chromosomes. The sequence characteristics, conversed motifs, exon/intron structures, classification, duplication events and cis-acting elements were also investigated. The gene duplication events showed that segmental duplication may play key roles in expansion of the PcoMYB gene family. Pyrus hopeiensis, which is a valuable wild resource, has strong cold resistance. An integrative analyses of miRNA and mRNA showed that PhMYB62 was involved in regulating low-temperature stress in P. hopeiensis flower organs. Subcellular localization analysis showed that PhMYB62 protein was specifically localized to the nucleus. The result of DAP-seq showed that PhMYB62 responded to low-temperature stress in P. hopeiensis by regulating TFs, which were associated with plant stress resistance, and POD, GAUT12, AUX28 and CHS genes. Subsequently, yeast one-hybrid verified that PhMYB62 could bind and activate the promoter of POD gene. The current study would provide a comprehensive information for further functional research on the stress-responsive R2R3-MYB gene candidates in pear, and may help to identify the genes associated with cold resistance for the cultivation of cold-resistant pear varieties.

Genome-wide identification and expression analysis of PtJAZ gene family in poplar (Populus trichocarpa).

BACKGROUND: The jasmonate ZIM domain (JAZ) protein is a key repressor of the jasmonate signal transduction pathway, which plays an important role in plant growth and development and defense responses. In this study, based on the published whole-genome data, we identified members of the JAZ gene family in Populus trichocarpa. Through a series of bioinformatic approaches, their expression patterns under various stress conditions have been analyzed to explore and excavate the endogenous resistance genes of poplar and provide a theoretical basis for breeding new varieties of poplar resistance. RESULTS: A total of 13 PtJAZ genes have been identified in P. trichocarpa and designated as PtJAZ1-PtJAZ13. Those 13 PtJAZ genes were unevenly distributed on nine chromosomes, and they could be divided into four subfamilies. The gene structures and motif composition of the members derived from the same subfamily were similar. Collinearity analysis demonstrated that, compared with Arabidopsis thaliana and Oryza sativa, the most collinear pairs (13) were found in P. trichocarpa and Eucalyptus robusta. Cis-acting element analysis suggested that the promoter regions of PtJAZs contained a large number of hormones and stress response elements, of which abscisic acid (ABA) and methyl jasmonate (MeJA) hormone response elements were the most abundant. The PtJAZ genes not only had diverse expression patterns in different tissues, but they also responded to various abiotic and biotic stress conditions. The co-expression network and GO and KEGG analyses showed that JAZ genes were closely related to insect resistance. CONCLUSIONS: In this study, applying bioinformatic methods, 13 PtJAZ gene family members from P. trichocarpa were identified and comprehensively analyzed. By further studying the function of the poplar JAZ gene family, the aim is to select genes with better insect resistance and stress resistance so as to lay a solid foundation for the subsequent breeding of new poplar varieties.

Identification of Differentially Expressed lncRNAs in Response to Blue Light and Expression Pattern Analysis of Populus tomentosa Hybrid Poplar 741.

Poplar is an important shelterbelt, timber stand, and city tree species that has been the focus of forestry research. The regulatory role of the long non-coding RNA molecule (lncRNA; length > 200 nt) has been a research hotspot in plants. In this study, seedlings of 741 poplar were irradiated with LED blue and white light, and the Illumina HiSeq 2000 sequencing platform was used to identify lncRNAs. |logFC| > 1 and p < 0.05 were considered to indicate differentially expressed lncRNAs, and nine differentially expressed lncRNAs were screened, the target genes of which were predicted, and three functionally annotated target genes were obtained. The differentially expressed lncRNAs were identified as miRNA targets. Six lncRNAs were determined to be target sites for twelve mRNAs in six miRNA families. LncRNAs and their target genes, including lncRNA MSTRG.20413.1-ptc-miR396e-5p-GRF9, were verified using quantitative real-time polymerase chain reaction analysis, and the expression patterns were analyzed. The analysis showed that the ptc-miR396e-5p expression was downregulated, while lncRNA MSTRG.20413.1 and GRF9 expression was upregulated, after blue light exposure. These results indicate that lncRNAs interact with miRNAs to regulate gene expression and affect plant growth and development.

Effect of T-DNA Integration on Growth of Transgenic Populus × euramericana cv. Neva Underlying Field Stands.

Multigene cotransformation has been widely used in the study of genetic improvement in crops and trees. However, little is known about the unintended effects and causes of multigene cotransformation in poplars. To gain insight into the unintended effects of T-DNA integration during multigene cotransformation in field stands, here, three lines (A1-A3) of Populus × euramericana cv. Neva (PEN) carrying Cry1Ac-Cry3A-BADH genes and three lines (B1-B3) of PEN carrying Cry1Ac-Cry3A-NTHK1 genes were used as research objects, with non-transgenic PEN as the control. Experimental stands were established at three common gardens in three locations and next generation sequencing (NGS) was used to identify the insertion sites of exogenous genes in six transgenic lines. We compared the growth data of the transgenic and control lines for four consecutive years. The results demonstrated that the tree height and diameter at breast height (DBH) of transgenic lines were significantly lower than those of the control, and the adaptability of transgenic lines in different locations varied significantly. The genotype and the experimental environment showed an interaction effect. A total of seven insertion sites were detected in the six transgenic lines, with B3 having a double-site insertion and the other lines having single copies. There are four insertion sites in the gene region and three insertion sites in the intergenic region. Analysis of the bases near the insertion sites showed that AT content was higher than the average chromosome content in four of the seven insertion sites within 1000 bp. Transcriptome analysis suggested that the differential expression of genes related to plant hormone transduction and lignin synthesis might be responsible for the slow development of plant height and DBH in transgenic lines. This study provides an integrated analysis of the unintended effects of transgenic poplar, which will benefit the safety assessment and reasonable application of genetically modified trees.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species.

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

Integrated physiological and metabolomic analyses reveal changes during the natural senescence of Quercus mongolica leaves.

Quercus mongolica is a common landscape, afforestation, and construction timber species in northern China with high ecological, economic, and ornamental value. Leaf senescence is a complex process that has important implications for plant growth and development. To explore changes of metabolites during the ageing of Quercus mongolica leaves, we investigated physiological responses and metabolite composition in ageing leaves harvested from 15-20-year-old Quercus mongolica. Leaf samples of Q. mongolica were collected when they were still green (at maturity) (stage 1), during early senescence (stage 2), and during late senescence (stage 3). These leaves were then subjected to physiological index and metabolome sequencing analyses. The physiological analysis showed that the leaves of Q. mongolica changed from green to yellow during senescence, which induced significant accumulation of soluble sugar and significant reductions in the concentration of soluble protein and chlorophyll. Peroxidase and catalase were the main antioxidant enzymes mitigating leaf senescence. Metabolomic analysis identified 797 metabolites during leaf senescence. Compared to stage 1, 70 differential metabolites were screened in stage 2 and 72 were screened in stage 3. Differential metabolites in the two senescent stages were principally enriched in amino acid metabolism, lipid metabolism and secondary metabolite biosynthesis. The contents of N-oleoylethanolamine and N, N-dimethylglycine were significantly increased only in stage 2, while the contents of trifolin, astragalin, valine, isoleucine, leucine, and citric acid were significantly increased only in stage 3. Histidine, homoserine, tryptophan, tyrosine, phenylalanine, proline, norleucine, N-glycyl-L-leucine, linoleic acid, linolenic acid, gallic acid, 3-indoleacrylic acid, 3-amino-2-naphthoic acid, 3-hydroxy-3-methylpentane-1,5-dioic acid, 2,3,4-trihydroxybenzoic acid, trifolin, astragalin, DL-2-aminoadipic acid, pinoresinol dimethyl ether, dimethylmatairesinol, and lysophosphatidylcholine increased during both stage 2 and stage 3. Increasing contents of these metabolites may constitute the main mechanism by which Q. mongolica leaves adapt to senescence.

Metabolomics and Transcriptomics Analyses Reveals the Molecular Regulatory Mechanisms of Walnut (Juglans regia L.) Embryos in Response to Shade Treatment.

The walnut is an important nut that has numerous uses worldwide. However, due to dwarf and close plantation methods as well as continuous cloudy or rainy days that occur during periods of walnut oil accumulation, the walnut fruit exhibits varying degrees of stress under low-light conditions. However, the effects of shade on metabolites and genes in walnut embryos remain unclear in the literature. The purpose of this study is to investigate the lipid biosynthesis process that occurs in walnut embryos under shade treatment via the use of metabolomics and transcriptomics analyses. The results indicate that the oil content decreases significantly under shaded conditions, while the protein content increases significantly. The expression levels of fatty acid desaturase 2 (FAD2) and stearoyl-ACP-desaturase (SAD) involved in the lipid biosynthesis mechanism were significantly reduced in the shaded group, which resulted in reductions in oleic (C18:1), linoleic (C18:2), and α-linolenic (C18:3) acids. The reduced oil content was consistent with the downregulation of genes associated with the lipid biosynthesis mechanism. In the amino acid biosynthesis process, the upregulated cysteine synthase (cscK) and anthranilate synthase beta subunit 2 (trpG) genes promoted the accumulation of L-aspartic acid and L-citrulline. The increase in protein content was consistent with the upregulation of genes related to amino acid biosynthesis. Thus, our study provides new insights into the regulatory mechanisms of shade underlying overall walnut fruit quality.

Progress in the understanding of WRKY transcription factors in woody plants.

The WRKY transcription factor (TF) family, named for its iconic WRKY domain, is among the largest and most functionally diverse TF families in higher plants. WRKY TFs typically interact with the W-box of the target gene promoter to activate or inhibit the expression of downstream genes; these TFs are involved in the regulation of various physiological responses. Analyses of WRKY TFs in numerous woody plant species have revealed that WRKY family members are broadly involved in plant growth and development, as well as responses to biotic and abiotic stresses. Here, we review the origin, distribution, structure, and classification of WRKY TFs, along with their mechanisms of action, the regulatory networks in which they are involved, and their biological functions in woody plants. We consider methods currently used to investigate WRKY TFs in woody plants, discuss outstanding problems, and propose several new research directions. Our objective is to understand the current progress in this field and provide new perspectives to accelerate the pace of research that enable greater exploration of the biological functions of WRKY TFs.

Conjoint Analysis of Genome-Wide lncRNA and mRNA Expression during the Salicylic Acid Response in Populus × euramericana.

Long noncoding RNAs (lncRNAs) participate in a wide range of biological processes, but lncRNAs in plants remain largely unknown; in particular, we lack a systematic identification of plant lncRNAs involved in hormone responses. To explore the molecular mechanism of the response of poplar to salicylic acid (SA), the changes in protective enzymes, which are closely related to plant resistance induced by exogenous SA, were studied, and the expression of mRNA and lncRNA were determined by high-throughput RNA sequencing. The results showed that the activities of phenylalanine ammonia lyase (PAL) and polyphenol oxidase (PPO), in the leaves of Populus × euramericana, were significantly increased by exogenous SA application. High-throughput RNA sequencing showed that 26,366 genes and 5690 lncRNAs were detected under the different treatment conditions: SA and H2O application. Among these, 606 genes and 49 lncRNAs were differentially expressed. According to target prediction, lncRNAs and target genes involved in light response, stress response, plant disease resistance, and growth and development, were differentially expressed in SA-treated leaves. Interaction analysis showed that lncRNA-mRNA interactions, following exogenous SA, were involved in the response of poplar leaves to the external environment. Our study provides a comprehensive view of Populus × euramericana lncRNAs and offers insights into the potential functions and regulatory interactions of SA-responsive lncRNAs, thus forming the foundation for future functional analysis of SA-responsive lncRNAs in Populus × euramericana.

Metabolomics and Transcriptomics Provide Insights into Lipid Biosynthesis in the Embryos of Walnut (Juglans regia L.).

Walnut (Juglans regia L.) is an important woody oilseed tree species due to its commercial value. However, the regulation mechanism of walnut oil accumulation is still poorly understood, which restricted the breeding and genetic improvement of high-quality oil-bearing walnuts. In order to explore the metabolic mechanism that regulates the synthesis of walnut oil, we used transcriptome sequencing technology and metabolome technology to comprehensively analyze the key genes and metabolites involved in oil synthesis of the walnut embryo at 60, 90, and 120 days after pollination (DAP). The results showed that the oil and protein contents increased gradually during fruit development, comprising 69.61% and 18.32% of the fruit, respectively, during ripening. Conversely, the contents of soluble sugar and starch decreased gradually during fruit development, comprising 2.14% and 0.84%, respectively, during ripening. Transcriptome sequencing generated 40,631 unigenes across 9 cDNA libraries. We identified 51 and 25 candidate unigenes related to the biosynthesis of fatty acid and the biosynthesis of triacylglycerol (TAG), respectively. The expression levels of the genes encoding Acetyl-CoA carboxylase (ACCase), long-chain acyl-CoA synthetases (LACS), 3-oxoacyl-ACP synthase II (KASII), and glycerol-3-phosphate acyl transfer (GPAT) were upregulated at 60 DAP relative to the levels at 90 and 120 DAP, while the stearoyl-ACP-desaturase (SAD) and fatty acid desaturase 2 (FAD2) genes were highly abundantly expressed during all walnut developmental periods. We found that ABSCISIC ACID INSENSEITIVE3 (ABI3), WRINKLEDl (WRI1), LEAFY COTYLEDON1 (LEC1), and FUSCA3 (FUS3) may be key transcription factors involved in lipid synthesis. Additionally, the metabolomics analysis detected 706 metabolites derived from 18 samples, among which, 4 are implicated in the TAG synthesis, 2 in the glycolysis pathway, and 5 in the tricarboxylic acid cycle (TCA cycle) pathway. The combined analysis of the related genes and metabolites in TAG synthesis showed that phospholipid:diacylglycerol acyltransferase (PDAT) genes were highly abundantly expressed across walnut fruit developmental periods, and their downstream metabolite TAG gradually accumulated with the progression of fruit development. The FAD2 gene showed consistently higher expression during fruit development, and its downstream metabolites 18:2-PC and 18:3-PC gradually accumulated. The ACCase, LACS, SAD, FAD2, and PDAT genes may be crucial genes required for walnut oil synthesis. Our data will enrich public databases and provide new insights into functional genes related to lipid metabolism in walnut.

Genome-Wide Characterization of the Mitogen-Activated Protein Kinase Gene Family and Their Expression Patterns in Response to Drought and Colletotrichum Gloeosporioides in Walnut (Juglans regia).

Mitogen-activated protein kinases (MAPKs) are a family of Ser/Thr (serine/threonine) protein kinases that play very important roles in plant responses to biotic and abiotic stressors. However, the MAPK gene family in the important crop walnut (Juglans regia L.) has been less well studied compared with other species. We discovered 25 JrMAPK members in the Juglans genome in this study. The JrMAPK gene family was separated into four subfamilies based on phylogenetic analysis, and members of the same subgroup had similar motifs and exons/introns. A variety of cis-acting elements, mainly related to the light response, growth and development, stress response, and hormone responses, were detected in the JrMAPK gene promoters. Collinearity analysis showed that purification selection was the main driving force in JrMAPK gene evolution, and segmental and tandem duplications played key roles in the expansion of the JrMAPK gene family. The RNA-Seq (RNA Sequencing) results indicated that many of the JrMAPK genes were expressed in response to different levels of Colletotrichum gloeosporioides infection. JrMAPK1, JrMAPK3, JrMAPK4, JrMAPK5, JrMAPK6, JrMAPK7, JrMAPK9, JrMAPK11, JrMAPK12, JrMAPK13, JrMAPK17, JrMAPK19, JrMAPK20, and JrMAPK21 were upregulated at the transcriptional level in response to the drought stress treatment. The results of this study will help in further investigations of the evolutionary history and biological functions of the MAPK gene family in walnut.

Complete chloroplast genomes and comparative analysis of Ligustrum species.

In this study, we assembled and annotated the chloroplast (cp) genomes of four Ligustrum species, L. sinense, L. obtusifolium, L. vicaryi, and L. ovalifolium 'Aureum'. Including six other published Ligustrum species, we compared various characteristics such as gene structure, sequence alignment, codon preference, and nucleic acid diversity, and performed positive-selection genes screening and phylogenetic analysis. The results showed that the cp genome of Ligustrum was 162,185-166,800 bp in length, with a circular tetrad structure, including a large single-copy region (86,885-90,106 bp), a small single-copy region (11,446-11,499 bp), and a pair of IRa and IRb sequences with the same coding but in opposite directions (31,608-32,624 bp). This structure is similar to the cp genomes of most angiosperms. We found 132-137 genes in the cp genome of Ligustrum, including 89-90 protein-coding genes, 35-39 tRNAs, and 8 rRNAs. The GC content was 37.93-38.06% and varied among regions, with the IR region having the highest content. The single-nucleotide (A/T)n was dominant in simple-sequence repeats of the Ligustrum cp genome, with an obvious A/T preference. Six hotspot regions were identified from multiple sequence alignment of Ligustrum; the ycf1 gene region and the clpP1 exon region can be used as potential DNA barcodes for the identification and phylogeny of the genus Ligustrum. Branch-site model and Bayes empirical Bayes (BEB) analysis showed that four protein-coding genes (accD, clpP, ycf1, and ycf2) were positively selected, and BEB analysis showed that accD and rpl20 had positively selected sites. A phylogenetic tree of Oleaceae species was constructed based on the whole cp genomes, and the results were consistent with the traditional taxonomic results. The phylogenetic results showed that genus Ligustrum is most closely related to genus Syringa. Our study provides important genetic information to support further investigations of the phylogenetic development and adaptive evolution of Ligustrum species.

Cold Resistance of Euonymus japonicus Beihaidao Leaves and Its Chloroplast Genome Structure and Comparison with Celastraceae Species.

Euonymus japonicus Beihaidao is one of the most economically important ornamental species of the Euonymus genus. There are approximately 97 genera and 1194 species of plants worldwide in this family (Celastraceae). Using E. japonicus Beihaidao, we conducted a preliminary study of the cold resistance of this species, evaluated its performance during winter, assembled and annotated its chloroplast genome, and performed a series of analyses to investigate its gene structure GC content, sequence alignment, and nucleic acid diversity. Our objectives were to understand the evolutionary relationships of the genus and to identify positive selection genes that may be related to adaptations to environmental change. The results indicated that E. japonicus Beihaidao leaves have certain cold resistance and can maintain their viability during wintering. Moreover, the chloroplast genome of E. japonicus Beihaidao is a typical double-linked ring tetrad structure, which is similar to that of the other four Euonymus species, E. hamiltonianus, E. phellomanus, E. schensianus, and E. szechuanensis, in terms of gene structure, gene species, gene number, and GC content. Compared to other Celastraceae species, the variation in the chloroplast genome sequence was lower, and the gene structure was more stable. The phylogenetic relationships of 37 species inferred that members of the Euonymus genus do not form a clade and that E. japonicus Beihaidao is closely related to E. japonicus and E. fortunei. A total of 11 functional positive selected genes were identified, which may have played an important role in the process of Celastraceae species adapting to environmental changes. Our study provides important genetic information to support further investigations into the phylogenetic development and adaptive evolution of Celastraceae species.

Complete chloroplast genome structure of four Ulmus species and Hemiptelea davidii and comparative analysis within Ulmaceae species.

In this study, the chloroplast (cp) genomes of Hemiptelea davidii, Ulmus parvifolia, Ulmus lamellosa, Ulmus castaneifolia, and Ulmus pumila 'zhonghuajinye' were spliced, assembled and annotated using the Illumina HiSeq PE150 sequencing platform, and then compared to the cp genomes of other Ulmus and Ulmaceae species. The results indicated that the cp genomes of the five sequenced species showed a typical tetrad structure with full lengths ranging from 159,113 to 160,388 bp. The large single copy (LSC), inverted repeat (IR), and small single copy (SSC) lengths were in the range of 87,736-88,466 bp, 26,317-26,622 bp and 18,485-19,024 bp, respectively. A total of 130-131 genes were annotated, including 85-86 protein-coding genes, 37 tRNA genes and eight rRNA genes. The GC contents of the five species were similar, ranging from 35.30 to 35.62%. Besides, the GC content was different in different region and the GC content in IR region was the highest. A total of 64-133 single sequence repeat (SSR) loci were identified among all 21 Ulmaceae species. The (A)n and (T)n types of mononucleotide were highest in number, and the lengths were primarily distributed in 10-12 bp, with a clear AT preference. A branch-site model and a Bayes Empirical Bayes analysis indicated that the rps15 and rbcL had the positive selection sites. Besides, the analysis of mVISTA and sliding windows got a lot of hotspots such as trnH/psbA, rps16/trnQ, trnS/trnG, trnG/trnR and rpl32/trnL, which could be utilized as potential markers for the species identification and phylogeny reconstruction within Ulmus in the further studies. Moreover, the evolutionary tree of Ulmaceae species based on common protein genes, whole cp genome sequences and common genes in IR region of the 23 Ulmaceae species were constructed using the ML method. The results showed that these Ulmaceae species were divided into two branches, one that included Ulmus, Zelkova and Hemiptelea, among which Hemiptelea was the first to differentiate and one that included Celtis, Trema, Pteroceltis, Gironniera and Aphananthe. Besides, these variations found in this study could be used for the classification, identification and phylogenetic study of Ulmus species. Our study provided important genetic information to support further investigations into the phylogenetic development and adaptive evolution of Ulmus and Ulmaceae species.

Effects of three regeneration methods on the growth and bacterial community diversity of Populus × euramericana.

To study the effects of different regeneration methods on the growth and bacterial community diversity of Populus × euramericana cv. '74/76' (poplar 107), we investigated the growth of poplar 107 trees under three regeneration methods in 2017 and 2020, and sequenced the 16S rDNA V5-V7 regions in stem endophytic, root endophytic, and rhizosphere soil bacteria present in samples from the three regeneration methods using the Illumina high-throughput sequencing platform. The growth analysis showed that stump grafting regeneration (ST) and stump sprouting regeneration (SP) presented similar tree height and diameter at breast height (DBH), which were significantly lower by planted seedling regeneration (CK). The high-throughput sequencing results showed that the rhizosphere soil bacteria appeared to be significantly more diverse and rich than the root and stem endophytic bacteria. Cluster analysis showed that the similarity of bacterial community structure among the rhizosphere soil, root, and stem was small. Thus, the three sample types showed significant differences in bacteria. While comparing the two years, 2020 was significantly more diverse and rich than 2017. With the increase in stand age, the abundance of Proteobacteria increased and the abundance of Acidobacteria decreased. Among the three regeneration methods, ST significantly increased the diversity of stem endophytic bacteria. Chthoniobacter was enriched in SP, which promoted the decomposition of organic matter, and more plant growth promoting rhizobacteria (PGPR) were accumulated in the rhizosphere of SP and ST. The composition of the bacterial community was similar in the three regeneration methods, but the community composition was different. Regeneration and transformation of poplar plantations can be better carried out by stump grafting and stump sprouting.

Integrated physiological and transcriptomic analyses reveal the molecular mechanism behind the response to cultivation in Quercus mongolica.

Quercus mongolica, a common tree species for building and landscaping in northern China, has great commercial and ecological value. The seedlings of Q. mongolica grow poorly and develop chlorosis when introduced from high-altitude mountains to low-altitude plains. Effective cultivation measures are key to improving the quality of seedlings. To investigate the complex responses of Q. mongolica to different cultivation measures, we compared the adaptability of 3-year-old Q. mongolica seedlings to pruning (P), irrigation (W), and fertilization [F (nitro compound fertilizer with 16N-16P-16K)]. Physiological measurements and transcriptome sequencing were performed on leaves collected under the P treatments (control, cutting, removal of all lateral branches, and removal of base branches to one-third of seedling height), the W treatments (0, 1, 2, 3, 4, or 5 times in sequence), and the F treatments (0, 2, 4, and 6 g/plant). Analyses of the physiological data showed that P was more effective than W or F for activating intracellular antioxidant systems. By contrast, W and F were more beneficial than P for inducing the accumulation of soluble sugar. OPLS-DA identified superoxide dismutase, malondialdehyde, and peroxidase as critical physiological indices for the three cultivation measures. Transcriptome analyses revealed 1,012 differentially expressed genes (DEGs) in the P treatment, 1,035 DEGs in the W treatment, and 1,175 DEGs in the F treatment; these DEGs were mainly enriched in Gene Ontology terms related to the stress response and signal transduction. Weighted gene coexpression network analyses indicated that specific gene modules were significantly correlated with MDA (one module) and soluble sugar (four modules). Functional annotation of the hub genes differentially expressed in MDA and soluble sugar-related modules revealed that Q. mongolica responded and adapted to different cultivation measures by altering signal transduction, hormone levels, reactive oxygen species, metabolism, and transcription factors. The hub genes HOP3, CIPK11, WRKY22, and BHLH35 in the coexpression networks may played a central role in responses to the cultivation practices. These results reveal the mechanism behind the response of Q. mongolica to different cultivation measures at the physiological and molecular levels and provide insight into the response of plants to cultivation measures.