Search

Determining how a new hybrid lineage can achieve reproductive isolation is a key to understanding the process and mechanisms of homoploid hybrid speciation. Here, we evaluated the degree and nature of reproductive isolation between the ecologically successful hybrid species Pinus densata and its parental species P. tabuliformis and P. yunnanensis. We performed interspecific crosses among the three species to assess their crossability. We then conducted reciprocal transplantation experiments to evaluate their fitness differentiation, and to examine how natural populations representing different directions of introgression differ in adaptation. The crossing experiments revealed weak genetic barriers among the species. The transplantation trials showed manifest evidence of local adaptation as the three species all performed best in their native habitats. Pinus densata populations from the western edge of its distribution have evolved a strong local adaptation to the specific habitat in that range; populations representing different directions of introgressants with the two parental species all showed fitness disadvantages in this P. densata habitat. These observations illustrate that premating isolation through selection against immigrants from other habitat types or postzygotic isolation through selection against backcrosses between the three species is strong. Thus, ecological selection in combination with endogenous components and geographic isolation has likely played a significant role in the speciation of P. densata. [ABSTRACT FROM AUTHOR]

In this note we prove a stability estimate for an inverse heat source problem. Based on the obtained stability estimate, we present a generalized Tikhonov regularization method and obtain the error estimate. Numerical experiment shows that the generalized Tikhonov regularization works well. [ABSTRACT FROM AUTHOR]

Southwest China is a biodiversity hotspot characterized by complex topography, heterogeneous regional climates and rich flora. The processes and driving factors underlying this hotspot remain to be explicitly tested across taxa to gain a general understanding of the evolution of biodiversity and speciation in the region. In this study, we examined the role played by historically neutral processes, geography and environment in producing the current genetic diversity of the subtropical pine Pinus yunnanensis. We used genetic and ecological methods to investigate the patterns of genetic differentiation and ecological niche divergence across the distribution range of this species. We found both continuous genetic differentiation over the majority of its range, and discrete isolated local clusters. The discrete differentiation between two genetic groups in the west and east peripheries is consistent with niche divergence and geographical isolation of these groups. In the central area of the species’ range, population structure was shaped mainly by neutral processes and geography rather than by ecological selection. These results show that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, and illustrate the importance of ecological factors in forming or maintaining genetic divergence across a complex landscape. Our findings differ from other phylogenetic studies that identified the historical drainage system in the region as the primary factor shaping population structure, and highlight the heterogeneous contributions that geography and environment have made to genetic diversity among taxa in southwest China. [ABSTRACT FROM AUTHOR]

Southwest China is a biodiversity hotspot characterized by complex topography, heterogeneous regional climates and rich flora. The processes and driving factors underlying this hotspot remain to be explicitly tested across taxa to gain a general understanding of the evolution of biodiversity and speciation in the region. In this study, we examined the role played by historically neutral processes, geography and environment in producing the current genetic diversity of the subtropical pine Pinus yunnanensis. We used genetic and ecological methods to investigate the patterns of genetic differentiation and ecological niche divergence across the distribution range of this species. We found both continuous genetic differentiation over the majority of its range, and discrete isolated local clusters. The discrete differentiation between two genetic groups in the west and east peripheries is consistent with niche divergence and geographical isolation of these groups. In the central area of the species’ range, population structure was shaped mainly by neutral processes and geography rather than by ecological selection. These results show that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, and illustrate the importance of ecological factors in forming or maintaining genetic divergence across a complex landscape. Our findings differ from other phylogenetic studies that identified the historical drainage system in the region as the primary factor shaping population structure, and highlight the heterogeneous contributions that geography and environment have made to genetic diversity among taxa in southwest China. [ABSTRACT FROM AUTHOR]

Pinus densata is an ecologically successful homoploid hybrid that inhabits vast areas of heterogeneous terrain on the south-eastern Tibetan Plateau as a result of multiple waves of colonization. Its region of origin, route of colonization onto the plateau and the directions of introgression with its parental species have previously been defined, but little is known about the isolation and divergence history of its populations. In this study, we surveyed nucleotide polymorphism over eight nuclear loci in 19 representative populations of P. densata and its parental species. Using this information and coalescence simulations, we assessed the historical changes in its population size, gene flow and divergence in time and space. The results indicate a late Miocene origin for P. densata associated with the recent uplift of south-eastern Tibet. The subsequent differentiation between geographical regions of this species began in the late Pliocene and was induced by regional topographical changes and Pleistocene glaciations. The ancestral P. densata population had a large effective population size but the central and western populations were established by limited founders, suggesting that there were severe bottlenecks during the westward migration out of the ancestral hybrid zone. After separating from their ancestral populations, population expansion occurred in all geographical regions especially in the western range. Gene flow in P. densata was restricted to geographically neighbouring populations, resulting in significant differentiation between regional groups. The new information on the divergence and demographic history of P. densata reported herein enhances our understanding of its speciation process on the Tibetan Plateau. [ABSTRACT FROM AUTHOR]

In this paper, a new method is put forward for simultaneous reliability and reliability-sensitivity analyses based on the information-reuse of sparse grid numerical integration (SGNI). First, the reliability analysis is conducted on the basis of fractional exponential moments-based maximum entropy method (FEM-MEM), where the SGNI is employed for FEM assessments. The reliability index can be evaluated by integrating over the distribution derived by FEM-MEM. Then, the reliability-sensitivity analysis is carried out by reusing the output samples in previous reliability analysis and updating the corresponding weights, where no additional model evaluations are required. Then, the FEM-MEM is applied again to derive the conditional distribution and reliability index. By comparing the conditional and original reliability indexes, one can define the reliability-sensitivity index to identify the importance of each random variable to reliability. Since only one-round of model evaluations are necessary in the proposed method, the computational efficiency is highly desirable. Four numerical examples are investigated to check the effectiveness of the proposed method, where pertinent results obtained from Monte Carlo simulations (MCS) and Sobol’s index are compared. The results demonstrate the proposed method is accurate and efficient for simultaneous reliability and reliability-sensitivity analyses. [ABSTRACT FROM AUTHOR]

Populus tomentosa, an indigenous tree species, is widely distributed and cultivated over 1,000,000 km2 in China, contributing significantly to forest production, ecological conservation and urban–rural greening. Although a reference genome is available for P. tomentosa, the intricate interspecific hybrid origins, chromosome structural variations (SVs) and sex determination mechanisms remain confusion and unclear due to its broad and even overlapping geographical distribution, extensive morphological variations and cross infiltration among white poplar species. We conducted a haplotype‐resolved de novo assembly of P. tomentosa elite individual GM107, which comprises subgenomes a and b with a total genome size of 714.9 Mb. We then analysed the formation of hybrid species and the phylogenetic evolution and sex differentiation across the entire genus. Phylogenomic analyses suggested that GM107 likely originated from a hybridisation event between P. alba (♀) and P. davidiana (♂) approximately 3.8 Mya. A total of 1551 chromosome SVs were identified between the two subgenomes. More noteworthily, a distinctive inversion structure spanning 2.15–2.95 Mb was unveiled among Populus, Tacamahaca, Turaga, Aigeiros poplar species and Salix, highlighting a unique evolutionary feature. Intriguingly, a novel sex genotype of the ZY type, which represents a crossover between XY and ZW systems, was identified and confirmed through both natural and artificial hybrids populations. These novel insights offer significant theoretical value for the study of the species' evolutionary origins and serve as a valuable resource for ecological genetics and forest biotechnology. [ABSTRACT FROM AUTHOR]

Background: Cupressus gigantea, a rare and endangered tree species with remarkable medicinal value, is endemic to the Tibetan Plateau. Yet, little is known about the underlying genetics of the unique ecological adaptability of this extremely long-lived conifer with a large genome size. Here, we present its first de novo and multi-tissue transcriptome in-depth characterization. Results: We performed Illumina paired-end sequencing and RNA libraries assembly derived from terminal buds, male and female strobili, biennial leaves, and cambium tissues taken from adult C. gigantea. In total, large-scale high-quality reads were assembled into 101,092 unigenes, with an average sequence length of 1029 bp, and 6848 unigenes (6.77%) were mapped against the KEGG databases to identify 292 pathways. A core set of 41,373 genes belonging to 2412 orthologous gene families shared between C. gigantea and nine other plants was revealed. In addition, we identified 2515 small to larger-size gene families containing in total 9223 genes specific to C. gigantea, and enriched for gene ontologies relating to biotic interactions. We identified an important terpene synthases gene family expansion with its 121 putative members. Conclusions: This study presents the first comprehensive transcriptome characterization of C. gigantea. Our results will facilitate functional genomic studies to support genetic improvement and conservation programs for this endangered conifer. [ABSTRACT FROM AUTHOR]

Paeonia ostii is a traditional ornamental and medicinal species that has attracted considerable interest for its high oil value. To facilitate the effective and rational cultivation and application of P. ostii in China, it is necessary to determine its potential spatial habitat distribution and environmental requirements. Using high-resolution environmental data for current and future climate scenarios, the potential suitable area and climatic requirements of P. ostii were modelled. Among the 11 environmental variables investigated, growing degree days, precipitation of the wettest month, mean temperature of the coldest quarter, global UV-B radiation, annual precipitation, and soil pH played major roles in determining the suitability of a habitat for the cultivation of P. ostii. Under the current environmental conditions in China, a total area of 20.31 × 105 km2 is suitable for growing P. ostii, accounting for 21.16% of the country's total land area. Under the two future climate scenario/year combinations (i.e., representative concentration pathways [RCPs], RCP2.6 and RCP8.5 in 2050), this species would increase its suitable area at high latitudes while decrease at low latitudes. These results present valuable information and a theoretical reference point for identifying the suitable cultivation areas of P. ostii. [ABSTRACT FROM AUTHOR]

Summary: Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated.First, we presented a high‐quality chromosome‐scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein‐coding genes. Second, our metabolome analyses uncovered pelargonidin‐3‐O‐glucoside, cyanidin‐3‐O‐arabinoside, and cyanidin‐3‐O‐glucoside as the main pigments involved in leaf color transition. Third, gene co‐expression further identified the MYB‐bHLH‐WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation.Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co‐expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein–protein and DNA–protein interaction assays.Our high‐quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus. [ABSTRACT FROM AUTHOR]

In this paper, stochastic dynamic response and reliability assessment of non-linear structures driven by non-stationary non-Gaussian ground motions are investigated. First, the non-stationary non-Gaussian ground motions are simulated by using the translation process theory. To circumvent the problem of a large number of iterations, an improved upgrading model is developed to obtain a power spectral density which is compatible with the target non-Gaussian probability distribution and is as close as possible to the target non-Gaussian power spectral density. To conduct stochastic dynamic response and reliability analyses of non-linear structures driven by non-stationary non-Gaussian ground motions, the probability density evolution method is employed, whereby a new high-dimensional points selection strategy based on the clustering algorithm is put forward. The new strategy of points selection is quite flexible and simple while taking into account both efficiency and accuracy. The efficacy and applicability of the proposed method are verified through a numerical example and a practical engineering example, in which the influence of non-Gaussianity on the probabilistic seismic response and reliability of the structure is discussed. • Nonstationary non-Gaussian ground motions are simulated with accuracy. • High-dimensional points selection is performed based on the clustering method. • Probabilistic dynamic analysis of nonlinear structures are carried out based on PDEM. • Numerical examples verify the efficacy and applicability of the proposed method. [ABSTRACT FROM AUTHOR]

Color change during flower opening is common; however, little is understood on the biochemical and molecular basis related. Lilac (Syringa oblata), a well-known woody ornamental plant with obvious petal color changes, is an ideal model. Here, we presented chromosome-scale genome assembly for lilac, resolved the flavonoids metabolism, and identified key genes and potential regulatory networks related to petal color change. The genome assembly is 1.05 Gb anchored onto 23 chromosomes, with a BUSCO score of 96.6%. Whole-genome duplication (WGD) event shared within Oleaceae was revealed. Metabolome quantification identified delphinidin-3-O-rutinoside (Dp3Ru) and cyanidin-3-O-rutinoside (Cy3Ru) as the major pigments; gene co-expression networks indicated WRKY an essential regulation factor at the early flowering stage, ERF more important in the color transition period (from violet to light nearly white), while the MBW complex participated in the entire process. Our results provide a foundation for functional study and molecular breeding in lilac. A high-quality genome assembly for the lilac (Syringa oblata) provides insight into the evolution of this shrub and potential mechanisms underlying petal color change. [ABSTRACT FROM AUTHOR]

Summary: Understanding the origin and distribution of genetic diversity across landscapes is critical for predicting the future of organisms in changing climates. This study investigated how adaptive and demographic forces have shaped diversity and population structure in Pinus densata, a keystone species on Qinghai‐Tibetan Plateau (QTP).We examined the distribution of genomic diversity across the range of P. densata using exome capture sequencing. We applied spatially explicit tests to dissect the impacts of allele surfing, geographic isolation and environmental gradients on population differentiation and forecasted how this genetic legacy may limit the persistence of P. densata in future climates.We found that allele surfing from range expansion could explain the distribution of 39% of the c. 48 000 genotyped single nucleotide polymorphisms (SNPs). Uncorrected, these allele frequency clines severely confounded inferences of selection. After controlling for demographic processes, isolation‐by‐environment explained 9.2–19.5% of the genetic structure, with c. 4.0% of loci being affected by selection. Allele surfing and genotype–environment associations resulted in genomic mismatch under projected climate scenarios.We illustrate that significant local adaptation, when coupled with reduced diversity as a result of demographic history, constrains potential evolutionary response to climate change. The strong signal of genomic vulnerability in P. densata may be representative for other QTP endemics. [ABSTRACT FROM AUTHOR]

miRNAs are important regulatory components involving in many biological processes, including plant development, vegetative and reproductive growth, and stress response. However, identification and characterization of miRNAs still remain limited for conifer species. In this study, with deep sequencing, we obtained 1,314,450 unique reads with 18–30 nt length from a stress-tolerant conifer, Sabina chinensis. We identified 37 conserved and 103 novel miRNAs, their unique characteristics were further analyzed, and 10 randomly selected were validated by qRT-PCR. Through miRNA target predictions and annotations, we found miRNA may have several targets as well a target could be regulated by several miRNAs, and a total of 2,397 mRNAs were predicted to be targets of the 140 miRNAs. These targets included not only important transcription factors such as auxin response factors, but also indispensable non-transcriptional factor proteins. Pathway-based analysis showed that S. chinensis miRNAs are involved in 172 metabolic pathways, of which 3 were discovered in adaptation-related pathways, indicating their possible relevance to the species’ stress-tolerance characteristics. This study is expected to lay the foundation for exploring the regulative roles of miRNAs in development, growth, and response to environmental stresses of S. chinensis. [ABSTRACT FROM AUTHOR]

Picea asperata and P. crassifolia have sympatric ranges and are closely related, but the differences between these species at the plastome level are unknown. To better understand the patterns of variation among Picea plastomes, the complete plastomes of P. asperata and P. crassifolia were sequenced. Then, the plastomes were compared with the complete plastomes of P. abies and P. morrisonicola, which are closely and distantly related to the focal species, respectively. We also used these sequences to construct phylogenetic trees to determine the relationships among and between the four species as well as additional taxa from Pinaceae and other gymnosperms. Analysis of our sequencing data allowed us to identify 438 single nucleotide polymorphism (SNPs) point mutation events, 95 indel events, four inversion events, and seven highly variable regions, including six gene spacer regions (psbJ-petA, trnT-psaM, trnS-trnD, trnL-rps4, psaC-ccsA, and rps7-trnL) and one gene (ycf1). The highly variable regions are appropriate targets for future use in the phylogenetic reconstructions of closely related, sympatric species of Picea as well as Pinaceae in general. [ABSTRACT FROM AUTHOR]

Pollination dynamics was studied in a first-generation Platycladusorientalis (L.) Franco seed orchard with pedigree reconstruction using eight nuclear and four chloroplast SSRs. The pedigree reconstruction assigned 371 of 448 studied seeds to one of the orchard’s 192 candidate male parents and showed a high level of outcrossing and pollen contamination in the orchard’s seed crop. While the orchard’s seed population showed greater allelic richness compared with the parental population, a few alleles present in the parental population were missing in the seed crop. Additionally, we detected no significant correlation between male reproductive energy (pollen yield) and male reproductive success; however, uneven parental contribution was also observed. Pollen management practices were recommended to ensure the maintenance of genetic diversity in the seed crops and increase in genetic gain. [ABSTRACT FROM AUTHOR]

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids., Competing Interests: Conflict of interest statement. None declared., (© The Author(s) 2024. Published by Oxford University Press on behalf of American Society of Plant Biologists.)

Long non-coding RNAs (lncRNAs) play essential roles in various biological processes, such as chromatin remodeling, post-transcriptional regulation, and epigenetic modifications. Despite their critical functions in regulating plant growth, root development, and seed dormancy, the identification of plant lncRNAs remains a challenge due to the scarcity of specific and extensively tested identification methods. Most mainstream machine learning-based methods used for plant lncRNA identification were initially developed using human or other animal datasets, and their accuracy and effectiveness in predicting plant lncRNAs have not been fully evaluated or exploited. To overcome this limitation, we retrained several models, including CPAT, PLEK, and LncFinder, using plant datasets and compared their performance with mainstream lncRNA prediction tools such as CPC2, CNCI, RNAplonc, and LncADeep. Retraining these models significantly improved their performance, and two of the retrained models, LncFinder-plant and CPAT-plant, alongside their ensemble, emerged as the most suitable tools for plant lncRNA identification. This underscores the importance of model retraining in tackling the challenges associated with plant lncRNA identification. Finally, we developed a pipeline (Plant-LncPipe) that incorporates an ensemble of the two best-performing models and covers the entire data analysis process, including reads mapping, transcript assembly, lncRNA identification, classification, and origin, for the efficient identification of lncRNAs in plants. The pipeline, Plant-LncPipe, is available at: https://github.com/xuechantian/Plant-LncRNA-pipline., Competing Interests: The authors declare no competing interests., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nanjing Agricultural University.)

Xanthoceras sorbifolium (yellowhorn) is a woody oil plant with super stress resistance and excellent oil characteristics. The yellowhorn oil can be used as biofuel and edible oil with high nutritional and medicinal value. However, genetic studies on yellowhorn are just in the beginning, and fundamental biological questions regarding its very long-chain fatty acid (VLCFA) biosynthesis pathway remain largely unknown. In this study, we reconstructed the VLCFA biosynthesis pathway and annotated 137 genes encoding relevant enzymes. We identified four oleosin genes that package triacylglycerols (TAGs) and are specifically expressed in fruits, likely playing key roles in yellowhorn oil production. Especially, by examining time-ordered gene co-expression network (TO-GCN) constructed from fruit and leaf developments, we identified key enzymatic genes and potential regulatory transcription factors involved in VLCFA synthesis. In fruits, we further inferred a hierarchical regulatory network with MYB-related ( XS03G0296800 ) and B3 ( XS02G0057600 ) transcription factors as top-tier regulators, providing clues into factors controlling carbon flux into fatty acids. Our results offer new insights into key genes and transcriptional regulators governing fatty acid production in yellowhorn, laying the foundation for efforts to optimize oil content and fatty acid composition. Moreover, the gene expression patterns and putative regulatory relationships identified here will inform metabolic engineering and molecular breeding approaches tailored to meet biofuel and bioproduct demands., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Shi, Ma, Wang, Liu, Yan, Tian, Li, Bao, Chen, Zhao, Xiang, Jia, Nie, Guan and Mao.)

Terpenes and terpenoids are key natural compounds for plant defense, development, and composition of plant oil. The synthesis and accumulation of a myriad of volatile terpenoid compounds in these plants may dramatically alter the quality and flavor of the oils, which provide great commercial utilization value for oil-producing plants. Terpene synthases ( TPSs ) are important enzymes responsible for terpenic diversity. Investigating the differentiation of the TPS gene family could provide valuable theoretical support for the genetic improvement of oil-producing plants. While the origin and function of TPS genes have been extensively studied, the exact origin of the initial gene fusion event - it occurred in plants or microbes - remains uncertain. Furthermore, a comprehensive exploration of the TPS gene differentiation is still pending. Here, phylogenetic analysis revealed that the fusion of the TPS gene likely occurred in the ancestor of land plants, following the acquisition of individual C- and N- terminal domains. Potential mutual transfer of TPS genes was observed among microbes and plants. Gene synteny analysis disclosed a differential divergence pattern between TPS-c and TPS-e/f subfamilies involved in primary metabolism and those (TPS-a/b/d/g/h subfamilies) crucial for secondary metabolites. Biosynthetic gene clusters (BGCs) analysis suggested a correlation between lineage divergence and potential natural selection in structuring terpene diversities. This study provides fresh perspectives on the origin and evolution of the TPS gene family., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Yan, Zhou, Liu, Zhao, Tian, Shi, Bao, Li, Jia, Nie, Guo, Kong, Porth and Mao.)

Amur honeysuckle (Lonicera maackii) is a widely used medicinal plant of the Caprifoliaceae family that produces chlorogenic acid. Research on this plant mainly focuses on its ornamental value and medicinal compounds, but a reference genome sequence and molecular resources for accelerated breeding are currently lacking. Herein, nanopore sequencing and high-throughput chromosome conformation capture (Hi-C) allowed a chromosome-level genome assembly of L. maackii (2n = 18). A global view of the gene regulatory network involved in the biosynthesis of chlorogenic acid and the dynamics of fruit coloration in L. maackii was established through metabolite profiling and transcriptome analyses. Moreover, we identified the genes encoding hydroxycinnamoyl-CoA quinate transferase (LmHQT) and hydroxycinnamoyl-CoA shikimic/quinate transferase (LmHCT), which localized to the cytosol and nucleus. Heterologous overexpression of these genes in Nicotiana benthamiana leaves resulted in elevated chlorogenic acid contents. Importantly, HPLC analyses revealed that LmHCT and LmHQTs recombinant proteins modulate the accumulation of chlorogenic acid (CGA) using quinic acid and caffeoyl CoA as substrates, highlighting the importance of LmHQT and LmHCT in CGA biosynthesis. These results confirmed that LmHQTs and LmHCT catalyze the biosynthesis of CGA in vitro. The genomic data presented in this study will offer a valuable resource for the elucidation of CGA biosynthesis and facilitating selective molecular breeding., Competing Interests: Conflict of interest statement. The authors declare that there is no conflict of interests., (© American Society of Plant Biologists 2023. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Coriaria nepalensis Wall. (Coriariaceae) is a nitrogen-fixing shrub which forms root nodules with the actinomycete Frankia. Oils and extracts of C. nepalensis have been reported to be bacteriostatic and insecticidal, and C. nepalensis bark provides a valuable tannin resource. Here, by combining PacBio HiFi sequencing and Hi-C scaffolding techniques, we generated a haplotype-resolved chromosome-scale genome assembly for C. nepalensis. This genome assembly is approximately 620 Mb in size with a contig N50 of 11 Mb, with 99.9% of the total assembled sequences anchored to 40 pseudochromosomes. We predicted 60,862 protein-coding genes of which 99.5% were annotated from databases. We further identified 939 tRNAs, 7,297 rRNAs, and 982 ncRNAs. The chromosome-scale genome of C. nepalensis is expected to be a significant resource for understanding the genetic basis of root nodulation with Frankia, toxicity, and tannin biosynthesis., (© 2023. The Author(s).)

Wood decay resistance (WDR) is marking the value of wood utilization. Many trees of the Lauraceae have exceptional WDR, as evidenced by their use in ancient royal palace buildings in China. However, the genetics of WDR remain elusive. Here, through comparative genomics, we revealed the unique characteristics related to the high WDR in Lauraceae trees. We present a 1.27-Gb chromosome-level assembly for Lindera megaphylla (Lauraceae). Comparative genomics integrating major groups of angiosperm revealed Lauraceae species have extensively shared gene microsynteny associated with the biosynthesis of specialized metabolites such as isoquinoline alkaloids, flavonoid, lignins and terpenoid, which play significant roles in WDR. In Lauraceae genomes, tandem and proximal duplications (TD/PD) significantly expanded the coding space of key enzymes of biosynthesis pathways related to WDR, which may enhance the decay resistance of wood by increasing the accumulation of these compounds. Among Lauraceae species, genes of WDR-related biosynthesis pathways showed remarkable expansion by TD/PD and conveyed unique and conserved motifs in their promoter and protein sequences, suggesting conserved gene collinearity, gene expansion and gene regulation supporting the high WDR. Our study thus reveals genomic profiles related to biochemical transitions among major plant groups and the genomic basis of WDR in the Lauraceae., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Tian, Guo, Yan, Shi, Nie, Zhao, Bao, Li, Kong, Su, Mao and Lin.)

Few incidents of ancient allopolyploidization (polyploidization by hybridization or merging diverged genomes) were previously revealed, although there is significant evidence for the accumulation of whole genome duplications (WGD) in plants. Here, we focused on Ericales, one of the largest and most diverse angiosperm orders with significant ornamental and economic value. Through integrating 24 high-quality whole genome data selected from ~ 200 Superasterids genomes/species and an algorithm of topology-based gene-tree reconciliation, we explored the evolutionary history of in Ericales with ancient complex. We unraveled the allopolyploid origin of Ericales and detected extensive lineage-specific gene loss following the polyploidization. Our study provided a new hypothesis regarding the origin of Ericales and revealed an instructive perspective of gene loss as a pervasive source of genetic variation and adaptive phenotypic diversity in Ericales., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Nie, Tian, Kong, Zhao, Chen, Jiao, El-Kassaby, Porth, Yang, Zhao and Mao.)

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle . Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy , has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle , while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii . This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nanjing Agricultural University.)

Among the three genomes in plant cells, the mitochondrial genome (mitogenome) is the least studied due to complex recombination and intergenomic transfer. In gymnosperms only ∼20 mitogenomes have been released thus far, which hinders a systematic investigation into the tempo and mode of mitochondrial DNA evolution in seed plants. Here, we report the complete mitogenome sequence of Platycladus orientalis (Cupressaceae). This mitogenome is assembled as two circular-mapping chromosomes with a size of ∼2.6 Mb and which contains 32 protein-coding genes, three rRNA and seven tRNA genes, and 1,068 RNA editing sites. Repetitive sequences, including dispersed repeats, transposable elements (TEs), and tandem repeats, made up 23% of the genome. Comparative analyses with 17 other mitogenomes representing the five gymnosperm lineages revealed a 30-fold difference in genome size, 80-fold in repetitive content, and 230-fold in substitution rate. We found dispersed repeats are highly associated with mitogenome expansion ( r = 0.99), and most of them were accumulated during recent duplication events. Syntenic blocks and shared sequences between mitogenomes decay rapidly with divergence time ( r = 0.53), with the exceptions of Ginkgo and Cycads which retained conserved genome structure over long evolutionary time. Our phylogenetic analysis supports a sister group relationship of Cupressophytes and Gnetophytes; both groups are unique in that they lost 8-12 protein-coding genes, of which 4-7 intact genes are likely transferred to nucleus. These two clades also show accelerated and highly variable substitution rates relative to other gymnosperms. Our study highlights the dynamic and enigmatic evolution of gymnosperm mitogenomes., Competing Interests: Author R-GZ was employed by the company Ori (Shandong) Gene Science and Technology Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Liu, Zhao, Zhang, Mao and Wang.)

Populus has a wide ecogeographical range spanning the Northern Hemisphere, and interspecific hybrids are common. Populus tomentosa Carr. is widely distributed and cultivated in the eastern region of Asia, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. Reference genomes are available for several Populus species, however, our goals were to produce a very high quality de novo chromosome-level genome assembly in P. tomentosa genome that could serve as a reference for evolutionary and ecological studies of hybrid speciation throughout the genus. Here, combining long-read sequencing and Hi-C scaffolding, we present a high-quality, haplotype-resolved genome assembly. The genome size was 740.2 Mb, with a contig N50 size of 5.47 Mb and a scaffold N50 size of 46.68 Mb, consisting of 38 chromosomes, as expected with the known diploid chromosome number (2n = 2x = 38). A total of 59,124 protein-coding genes were identified. Phylogenomic analyses revealed that P. tomentosa is comprised of two distinct subgenomes, which we deomonstrate is likely to have resulted from hybridization between Populus adenopoda as the female parent and Populus alba var. pyramidalis as the male parent, with an origin of approximately 3.93 Ma. Although highly colinear, significant structural variation was found between the two subgenomes. Our study provides a valuable resource for ecological genetics and forest biotechnology., (© 2021 John Wiley & Sons Ltd.)

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species., (© 2021. The Author(s).)

LTR retrotransposons (LTR-RTs) are ubiquitous and represent the dominant repeat element in plant genomes, playing important roles in functional variation, genome plasticity and evolution. With the advent of new sequencing technologies, a growing number of whole-genome sequences have been made publicly available, making it possible to carry out systematic analyses of LTR-RTs. However, a comprehensive and unified annotation of LTR-RTs in plant groups is still lacking. Here, we constructed a plant intact LTR-RTs dataset, which is designed to classify and annotate intact LTR-RTs with a standardized procedure. The dataset currently comprises a total of 2,593,685 intact LTR-RTs from genomes of 300 plant species representing 93 families of 46 orders. The dataset is accompanied by sequence, diverse structural and functional annotation, age determination and classification information associated with the LTR-RTs. This dataset will contribute valuable resources for investigating the evolutionary dynamics and functional implications of LTR-RTs in plant genomes., (© 2021. The Author(s).)

Azaleas (Ericaceae) comprise one of the most diverse ornamental plants, renowned for their cultural and economic importance. We present a chromosome-scale genome assembly for Rhododendron simsii, the primary ancestor of azalea cultivars. Genome analyses unveil the remnants of an ancient whole-genome duplication preceding the radiation of most Ericaceae, likely contributing to the genomic architecture of flowering time. Small-scale gene duplications contribute to the expansion of gene families involved in azalea pigment biosynthesis. We reconstruct entire metabolic pathways for anthocyanins and carotenoids and their potential regulatory networks by detailed analysis of time-ordered gene co-expression networks. MYB, bHLH, and WD40 transcription factors may collectively regulate anthocyanin accumulation in R. simsii, particularly at the initial stages of flower coloration, and with WRKY transcription factors controlling progressive flower coloring at later stages. This work provides a cornerstone for understanding the underlying genetics governing flower timing and coloration and could accelerate selective breeding in azalea.

Platycladus orientalis is an ecologically important native conifer in Northern China and exotic species in many parts of the world; however, knowledge about the species' genetics and genome are very limited. The availability of well-developed battery of genetic markers, with large genome coverage, is a prerequisite for the species genetic dissection of adaptive attributes and efficient selective breeding. Here, we present a genome-wide genotyping method with double-digestion restriction site associated DNA sequencing (ddRAD-seq) that is effective in generating large number of Mendelian markers for genome mapping and other genetic applications. Using 139 megagametophytes collected from a single mother tree, we assembled 397,226 loci, of which 108,683 (27.4%) were polymorphic. After stringent filtering for 1:1 segregation ratio and missing rate of <20%, the remaining 23,926 loci (22% of the polymorphic loci) were ordered into 11 linkage groups (LGs) and distributed across 7,559 unique positions, with a total map length of 1,443 cM and an average spacing of 0.2 cM between adjacent unique positions. The 11 LGs correspond to the species' 11 haploid genome chromosome number. This genetic map is among few high-density maps available for conifers to date, and represents the first genetic map for P. orientalis The information generated serves as a solid foundation not only for marker-assisted breeding efforts, but also for comparative conifer genomic studies., (Copyright © 2019 Jin et al.)

Background: Malania oleifera, a member of the Olacaceae family, is an IUCN red listed tree, endemic and restricted to the Karst region of southwest China. This tree's seed is valued for its high content of precious fatty acids (especially nervonic acid). However, studies on its genetic makeup and fatty acid biogenesis are severely hampered by a lack of molecular and genetic tools., Findings: We generated 51 Gb and 135 Gb of raw DNA sequences, using Pacific Biosciences (PacBio) single-molecule real-time and 10× Genomics sequencing, respectively. A final genome assembly, with a scaffold N50 size of 4.65 Mb and a total length of 1.51 Gb, was obtained by primary assembly based on PacBio long reads plus scaffolding with 10× Genomics reads. Identified repeats constituted ∼82% of the genome, and 24,064 protein-coding genes were predicted with high support. The genome has low heterozygosity and shows no evidence for recent whole genome duplication. Metabolic pathway genes relating to the accumulation of long-chain fatty acid were identified and studied in detail., Conclusions: Here, we provide the first genome assembly and gene annotation for M. oleifera. The availability of these resources will be of great importance for conservation biology and for the functional genomics of nervonic acid biosynthesis., (© The Author(s) 2019. Published by Oxford University Press.)

Background: Long non-coding RNAs (lncRNAs) play widespread roles in gene regulation and cellular processes. However, the functional roles of lncRNAs in hepatocellular carcinoma (HCC) are not yet well elucidated. The aim of the present study was to measure the levels of lncRNA PCAT-1 expression in HCC and evaluate its clinical significance in the development and progression of HCC., Methods: We examined the expression of PCAT-1 in 117 HCC tissues and adjacent non-tumor tissues using quantitative real-time-PCR and analyzed its correlation with the clinical parameters., Results: Our data showed that PCAT-1 expression in HCC tissues was significantly increased compared with adjacent non-tumor tissues (P<0.05). Up-regulated expression of PCAT-1 was significantly associated with TNM stage and metastasis (P<0.05), but not other clinical parameters. Moreover, Kaplan-Meier survival analysis showed that a high expression level of PCAT-1 resulted in a significantly poor overall survival of HCC patients. The multivariate Cox regression analysis demonstrated that PCAT-1 expression level was an independent prognostic factor for the overall survival rate of HCC patients., Conclusions: Our data suggested that the increased expression of PCAT-1 was associated with advanced clinical parameters and poor overall survival of HCC patients, indicating that PCAT-1 up-regulation may serve as a novel biomarker of poor prognosis in HCC patients.

Objective: To investigate the prevalence and antibiotic resistance of Shigella isolated from children with diarrhea for the guidance of clinical treatment and prevention and control of bacillary dysentery., Method: A total of 156 strains of Shigella were isolated from feces of children with diarrhea in Zhejiang Xiaoshan Hospital from January 2008 to December 2010. The antimicrobial resistance of the strains was detected by disk diffusion method and the extended-spectrum beta-lactamases (ESBLs) in these isolates were determined using phenotypic confirmatory test; the isolates of ESBLs producing Shigella sonnei were analyzed by REP-PCR., Result: Among 156 strains of Shigella isolated, the most common groups were Shigella sonnei (130 strains, accounting for 83.3%) and Shigella fleaneri (26 strains, accounting for 16.7%), and 81 (51.9%) strains were identified as ESBLs producers, and the positive rates in 2008, 2009 and 2010 were 32.0%, 41.4% and 59.8%, respectively. The results of antibiotic susceptibility test displayed that the resistance rates of ESBLs producing Shigella to ampicillin, cotrimoxazole, cefotaxime, piperacillin were higher than 90%. However, the resistance rates to cefepime, ceftazidime, levofloxacin and ciprofloxacin were low; The resistance of ESBLs producing strains to piperacillin (100% vs. 77.3%), cefotaxime (100% vs. 0), ceftazidime (14.8% vs. 0), cefepime (28.4% vs. 0), cotrimoxazole (95.1% vs. 86.7%) was significantly higher than that of non-ESBLs producing strains (χ(2) = 20.605, 156.000, 12.037, 24.979, 45.040, respectively; P < 0.05). No isolate was resistant to piperacillin/tazobactam and imipenem. There were 7 genotypes among 74 ESBLs producing Shigella sonnei, respectively type A (50), type B (12), type C (8), type D (1), type E (1), type F (1), and type G (1)., Conclusion: The isolation rate of ESBLs-producing isolate was high in Shigella from pediatric patients with diarrhea, and the number is going up year by year, and these ESBLs producing Shigella sonnei strains in genotype A are dominant in recent years, Piperacillin/tazobactam is the drug of choice for children with ESBLs producing Shigella infection.