Genome-wide identification of oxidosqualene cyclase genes regulating natural rubber in Taraxacum kok-saghyz.

MAIN CONCLUSION: Nine TkOSC genes have been identified by genome-wide screening. Among them, TkOSC4-6 might be more crucial for natural rubber biosynthesis in Taraxacum kok-saghyz roots. Taraxacum kok-saghyz Rodin (TKS) roots contain large amounts of natural rubber, inulin, and valuable metabolites. Oxidosqualene cyclase (OSC) is a key member for regulating natural rubber biosynthesis (NRB) via the triterpenoid biosynthesis pathway. To explore the functions of OSC on natural rubber producing in TKS, its gene family members were identified in TKS genome via genome-wide screening. Nine TkOSCs were identified, which were mainly distributed in the cytoplasm. Their family genes experienced a neutral selection during the evolution process. Overall sequence homology analysis OSC proteins revealed 80.23% similarity, indicating a highly degree of conservation. Pairwise comparisons showed a multiple sequence similarity ranging from 57% to 100%. Protein interaction prediction revealed that TkOSCs may interact with baruol synthase, sterol 1,4-demethylase, lupeol synthase and squalene epoxidase. Phylogenetic analysis showed that OSC family proteins belong to two branches. TkOSC promoter regions contain cis-acting elements related to plant growth, stress response, hormones response and light response. Protein accumulation analysis demonstrated that TkOSC4, TkOSC5 and TkOSC6 proteins had strong expression levels in the root, latex and plumular axis. Comparison of gene expression patterns showed TkOSC1, TkOSC4, TkOSC5, TkOSC6, TkOSC7, TkOSC8 and TkOSC9 might be important in regulating NRB. Combination of gene and protein results revealed TkOSC4-6 might be more crucial, and the data might contribute to a more profound understanding of the roles of OSCs for NRB in TKS roots.

Genome-wide epigenetic dynamics of tea leaves under mechanical wounding stress during oolong tea postharvest processing.

Understanding the epigenetic responses to mechanical wounding stress during the postharvest processing of oolong tea provides insight into the reprogramming of the tea genome and its impact on tea quality. Here, we characterized the 5mC DNA methylation and chromatin accessibility landscapes of tea leaves subjected to mechanical wounding stress during the postharvest processing of oolong tea. Analysis of the differentially methylated regions and preferentially accessible promoters revealed many overrepresented TF-binding motifs, highlighting sets of TFs that are likely important for the quality of oolong tea. Within these sets, we constructed a chromatin accessibility-mediated gene regulatory network specific to mechanical wounding stress. In combination with the results of the TF-centred yeast one-hybrid assay, we identified potential binding sites of CsMYC2 and constructed a gene regulatory network centred on CsMYC2, clarifying the potential regulatory role of CsMYC2 in the postharvest processing of oolong tea. Interestingly, highly accessible chromatin and hypomethylated cytosine were found to coexist in the promoter region of the indole biosynthesis gene (tryptophan synthase ß-subunit, CsTSB) under wounding stress, which indicates that these two important epigenetic regulatory mechanisms are jointly involved in regulating the synthesis of indole during the postharvest processing of oolong tea. These findings improve our understanding of the epigenetic regulatory mechanisms involved in quality formation during the postharvest processing of oolong tea.

The first Chromosomal-level genome assembly of Sageretia thea using Nanopore long reads and Pore-C technology.

Sageretia thea, a notable species within the mock buckthorn genus, is recognized for its intriguing biogeographical distribution and diverse medicinal properties. Despite this significance, genomic studies on S. thea are still in the nascent stages. We present the first chromosome-level genome assembly of S. thea that was generated using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies complemented by Pore-C chromatin conformation capture. The genome assembly had a size of 197.8 Mb with 12 chromosomal scaffolds and a scaffold N50 length of 15.9 Mb. A total of 25,434 protein-coding genes were identified and functionally annotated, and the gene model indicated 96.5% complete eukaryotic BUSCOs. Additionally, orthologous gene profiling and synteny analysis were performed to elucidate the evolutionary relationships within the Rhamnaceae family and Rosales. This high-quality chromosomal genome is the first genomic view of S. thea, which will serve as the basis for future studies on its biological and medicinal properties, and evolutionary history.

Genome editing for improvement of biotic and abiotic stress tolerance in cereals.

Global agricultural production must quadruple by 2050 to fulfil the needs of a growing global population, but climate change exacerbates the difficulty. Cereals are a very important source of food for the world population. Improved cultivars are needed, with better resistance to abiotic stresses like drought, salt, and increasing temperatures, and resilience to biotic stressors like bacterial and fungal infections, and pest infestation. A popular, versatile, and helpful method for functional genomics and crop improvement is genome editing. Rapidly developing genome editing techniques including clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein (Cas) are very important. This review focuses on how CRISPR/Cas9 genome editing might enhance cereals' agronomic qualities in the face of climate change, providing important insights for future applications. Genome editing efforts should focus on improving characteristics that confer tolerance to conditions exacerbated by climate change (e.g. drought, salt, rising temperatures). Improved water usage efficiency, salt tolerance, and heat stress resilience are all desirable characteristics. Cultivars that are more resilient to insect infestations and a wide range of biotic stressors, such as bacterial and fungal diseases, should be created. Genome editing can precisely target genes linked to disease resistance pathways to strengthen cereals' natural defensive systems.

Revolutionizing soybean genomics: How CRISPR and advanced sequencing are unlocking new potential.

Soybean Glycine max L., paleopolyploid genome, poses challenges to its genetic improvement. However, the development of reference genome assemblies and genome sequencing has completely changed the field of soybean genomics, allowing for more accurate and successful breeding techniques as well as research. During the single-cell revolution, one of the most advanced sequencing tools for examining the transcriptome landscape is single-cell RNA sequencing (scRNA-seq). Comprehensive resources for genetic improvement of soybeans may be found in the SoyBase and other genomics databases. CRISPR-Cas9 genome editing technology provides promising prospects for precise genetic modifications in soybean. This method has enhanced several soybean traits, including as yield, nutritional value, and resistance to both biotic and abiotic stresses. With base editing techniques that allow for precise DNA modifications, the use of CRISPR-Cas9 is further increased. With the availability of the reference genome for soybeans and the following assembly of wild and cultivated soybeans, significant chromosomal rearrangements and gene duplication events have been identified, offering new perspectives on the complex genomic structure of soybeans. Furthermore, major single nucleotide polymorphisms (SNPs) linked to stachyose and sucrose content have been found through genome-wide association studies (GWAS), providing important tools for enhancing soybean carbohydrate profiles. In order to open up new avenues for soybean genetic improvement, future research approaches include investigating transcriptional divergence processes, enhancing genetic resources, and incorporating CRISPR-Cas9 technologies.

Wheat improvement through advances in single nucleotide polymorphism (SNP) detection and genotyping with a special emphasis on rust resistance.

KEY MESSAGE: Single nucleotide polymorphism (SNP) markers in wheat and their prospects in breeding with special reference to rust resistance. Single nucleotide polymorphism (SNP)-based markers are increasingly gaining momentum for screening and utilizing vital agronomic traits in wheat. To date, more than 260 million SNPs have been detected in modern cultivars and landraces of wheat. This rapid SNP discovery was made possible through the release of near-complete reference and pan-genome assemblies of wheat and its wild relatives, coupled with whole genome sequencing (WGS) of thousands of wheat accessions. Further, genotyping customized SNP sites were facilitated by a series of arrays (9 to 820Ks), a cost effective substitute WGS. Lately, germplasm-specific SNP arrays have been introduced to characterize novel traits and detect closely linked SNPs for marker-assisted breeding. Subsequently, the kompetitive allele-specific PCR (KASP) assay was introduced for rapid and large-scale screening of specific SNP markers. Moreover, with the advances and reduction in sequencing costs, ample opportunities arise for generating SNPs artificially through mutations and in combination with next-generation sequencing and comparative genomic analyses. In this review, we provide historical developments and prospects of SNP markers in wheat breeding with special reference to rust resistance where over 50 genetic loci have been characterized through SNP markers. Rust resistance is one of the most essential traits for wheat breeding as new strains of the Puccinia fungus, responsible for rust diseases, evolve frequently and globally.

Haplotype-based pangenomes reveal genetic variations and climate adaptations in moso bamboo populations.

Moso bamboo (Phyllostachys edulis), an ecologically and economically important forest species in East Asia, plays vital roles in carbon sequestration and climate change mitigation. However, intensifying climate change threatens moso bamboo survival. Here we generate high-quality haplotype-based pangenome assemblies for 16 representative moso bamboo accessions and integrated these assemblies with 427 previously resequenced accessions. Characterization of the haplotype-based pangenome reveals extensive genetic variation, predominantly between haplotypes rather than within accessions. Many genes with allele-specific expression patterns are implicated in climate responses. Integrating spatiotemporal climate data reveals more than 1050 variations associated with pivotal climate factors, including temperature and precipitation. Climate-associated variations enable the prediction of increased genetic risk across the northern and western regions of China under future emissions scenarios, underscoring the threats posed by rising temperatures. Our integrated haplotype-based pangenome elucidates moso bamboo's local climate adaptation mechanisms and provides critical genomic resources for addressing intensifying climate pressures on this essential bamboo. More broadly, this study demonstrates the power of long-read sequencing in dissecting adaptive traits in climate-sensitive species, advancing evolutionary knowledge to support conservation.

An Improved Chromosome-scale Genome Assembly and Population Genetics resource for Populus tremula.

Aspen (Populus tremula L.) is a keystone species and a model system for forest tree genomics. We present an updated resource comprising a chromosome-scale assembly, population genetics and genomics data. Using the resource, we explore the genetic basis of natural variation in leaf size and shape, traits with complex genetic architecture. We generated the genome assembly using long-read sequencing, optical and high-density genetic maps. We conducted whole-genome resequencing of the Umeå Aspen (UmAsp) collection. Using the assembly and re-sequencing data from the UmAsp, Swedish Aspen (SwAsp) and Scottish Aspen (ScotAsp) collections we performed genome-wide association analyses (GWAS) using Single Nucleotide Polymorphisms (SNPs) for 26 leaf physiognomy phenotypes. We conducted Assay of Transposase Accessible Chromatin sequencing (ATAC-Seq), identified genomic regions of accessible chromatin, and subset SNPs to these regions, improving the GWAS detection rate. We identified candidate long non-coding RNAs in leaf samples, quantified their expression in an updated co-expression network, and used this to explore the functions of candidate genes identified from the GWAS. A GWAS found SNP associations for seven traits. The associated SNPs were in or near genes annotated with developmental functions, which represent candidates for further study. Of particular interest was a ~177-kbp region harbouring associations with several leaf phenotypes in ScotAsp. We have incorporated the assembly, population genetics, genomics, and GWAS data into the PlantGenIE.org web resource, including updating existing genomics data to the new genome version, to enable easy exploration and visualisation. We provide all raw and processed data to facilitate reuse in future studies.

Atlas of telomeric repeat diversity in Arabidopsis thaliana.

BACKGROUND: Telomeric repeat arrays at the ends of chromosomes are highly dynamic in composition, but their repetitive nature and technological limitations have made it difficult to assess their true variation in genome diversity surveys. RESULTS: We have comprehensively characterized the sequence variation immediately adjacent to the canonical telomeric repeat arrays at the very ends of chromosomes in 74 genetically diverse Arabidopsis thaliana accessions. We first describe several types of distinct telomeric repeat units and then identify evolutionary processes such as local homogenization and higher-order repeat formation that shape diversity of chromosome ends. By comparing largely isogenic samples, we also determine repeat number variation of the degenerate and variant telomeric repeat array at both the germline and somatic levels. Finally, our analysis of haplotype structure uncovers chromosome end-specific patterns in the distribution of variant telomeric repeats, and their linkage to the more proximal non-coding region. CONCLUSIONS: Our findings illustrate the spectrum of telomeric repeat variation at multiple levels in A. thaliana-in germline and soma, across all chromosome ends, and across genetic groups-thereby expanding our knowledge of the evolution of chromosome ends.

Unravelling DNA methylation dynamics during developmental stages in Quercus ilex subsp. ballota [Desf.] Samp.

BACKGROUND: DNA methylation is a critical factor influencing plant growth, adaptability, and phenotypic plasticity. While extensively studied in model and crop species, it remains relatively unexplored in holm oak and other non-domesticated forest trees. This study conducts a comprehensive in-silico mining of DNA methyltransferase and demethylase genes within the holm oak genome to enhance our understanding of this essential process in these understudied species. The expression levels of these genes in adult and seedling leaves, as well as embryos, were analysed using quantitative real-time PCR (qRT-PCR). Global DNA methylation patterns were assessed through methylation-sensitive amplified polymorphism (MSAP) techniques. Furthermore, specific methylated genomic sequences were identified via MSAP sequencing (MSAP-Seq). RESULT: A total of 13 DNA methyltransferase and three demethylase genes were revealed in the holm oak genome. Expression levels of these genes varied significantly between organs and developmental stages. MSAP analyses revealed a predominance of epigenetic over genetic variation among organs and developmental stages, with significantly higher global DNA methylation levels observed in adult leaves. Embryos exhibited frequent demethylation events, while de novo methylation was prevalent in seedling leaves. Approximately 35% of the genomic sequences identified by MSAP-Seq were methylated, predominantly affecting nuclear genes and intergenic regions, as opposed to repetitive sequences and chloroplast genes. Methylation was found to be more pronounced in the exonic regions of nuclear genes compared to their promoter and intronic regions. The methylated genes were predominantly associated with crucial biological processes such as photosynthesis, ATP synthesis-coupled electron transport, and defence response. CONCLUSION: This study opens a new research direction in analysing variability in holm oak by evaluating the epigenetic events and mechanisms based on DNA methylation. It sheds light on the enzymatic machinery governing DNA (de)methylation, and the changes in the expression levels of methylases and demethylases in different organs along the developmental stages. The expression level was correlated with the DNA methylation pattern observed, showing the prevalence of de novo methylation and demethylation events in seedlings and embryos, respectively. Several methylated genes involved in the regulation of transposable element silencing, lipid biosynthesis, growth and development, and response to biotic and abiotic stresses are highlighted. MSAP-seq integrated with whole genome bisulphite sequencing and advanced sequencing technologies, such as PacBio or Nanopore, will bring light on epigenetic mechanisms regulating the expression of specific genes and its correlation with the phenotypic variability and the differences in the response to environmental cues, especially those related to climate change.

Genome-wide investigation of the PIF gene family in alfalfa (Medicago sativa L.) expression profiles during development and stress.

BACKGROUND: Phytochrome-interacting factors (PIFs) plays an important role in plants as hubs for intracellular signaling regulation. The PIF gene family has been identified and characterized in many plants, but alfalfa (Medicago sativa L.), an important perennial high-quality legume forage, has not been reported on the PIF gene family. RESULTS: In this study, we presented the identification and characterization of five MsPIF genes in alfalfa (Medicago sativa L.). Phylogenetic analysis indicated that PIFs from alfalfa and other four plant species could be divided into three groups supported by similar motif analysis. The collinearity analysis of the MsPIF gene family showed the presence of two gene pairs, and the collinearity analysis with AtPIFs showed three gene pairs, indicating that the evolutionary process of this family is relatively conservative. Analysis of cis-acting elements in promoter regions of MsPIF genes indicated that various elements were related to light, abiotic stress, and plant hormone responsiveness. Gene expression analyses demonstrated that MsPIFs were primarily expressed in the leaves and were induced by various abiotic stresses. CONCLUSION: This study conducted genome-wide identification, evolution, synteny analysis, and expression analysis of the PIFs in alfalfa. Our study lays a foundation for the study of the biological functions of the PIF gene family and provides a useful reference for improving abiotic stress resistance in alfalfa.

A near-complete chromosome-level genome assembly of looseleaf lettuce (Lactuca sativa var. crispa).

Lettuce (Lactuca sativa L., Asteraceae) is one of the most important vegetable crops, known for its various horticultural types and significant morphological variation. The first reference genome of lettuce, a crisphead type (L. sativa var. capitata cv. Salinas), was previously released. Here, we reported a near-complete chromosome-level reference genome for looseleaf lettuce (L. sativa var. crispa). PacBio high-fidelity sequencing, Oxford Nanopore, and Hi-C technologies were employed to produce genome assembly. The final assembly is 2.59 Gb in length with a contig N50 of 205.47 Mb, anchored onto nine chromosomes, containing 14 recognizable telomeres and only 11 gaps. Repetitive sequences account for 77.11% of the genome, and 41,375 protein-coding genes were predicted, with 99.10% of these assigned functional annotations. This chromosome-level genome enriched genomic resources for various horticultural types of lettuce and will facilitate the characterization of morphological variation and genetic improvement in lettuce.

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion.

BACKGROUND: Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions. RESULTS: Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP". CONCLUSION: The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

Genome-Wide Analysis and Expression Profiling of Soybean RbcS Family in Response to Plant Hormones and Functional Identification of GmRbcS8 in Soybean Mosaic Virus.

Rubisco small subunit (RbcS), a core component with crucial effects on the structure and kinetic properties of the Rubisco enzyme, plays an important role in response to plant growth, development, and various stresses. Although Rbcs genes have been characterized in many plants, their muti-functions in soybeans remain elusive. In this study, a total of 11 GmRbcS genes were identified and subsequently divided into three subgroups based on a phylogenetic relationship. The evolutionary analysis revealed that whole-genome duplication has a profound effect on GmRbcSs. The cis-acting elements responsive to plant hormones, development, and stress-related were widely found in the promoter region. Expression patterns based on the RT-qPCR assay exhibited that GmRbcS genes are expressed in multiple tissues, and notably Glyma.19G046600 (GmRbcS8) exhibited the highest expression level compared to other members, especially in leaves. Moreover, differential expressions of GmRbcS genes were found to be significantly regulated by exogenous plant hormones, demonstrating their potential functions in diverse biology processes. Finally, the function of GmRbcS8 in enhancing soybean resistance to soybean mosaic virus (SMV) was further determined through the virus-induced gene silencing (VIGS) assay. All these findings establish a strong basis for further elucidating the biological functions of RbcS genes in soybeans.

Genome-Wide Identification of microRNAs Associated with Starch Biosynthesis and Endosperm Development in Foxtail Millet.

Foxtail millet is one of the oldest crops, and its endosperm contains up to 70% of starch. Grain filling is an important starch accumulation process associated with foxtail millet yield and quality. However, the molecular mechanisms of grain filling in foxtail millet are relatively unclear. Here, we investigate the genes and regulated miRNAs associated with starch synthesis and metabolism in foxtail millet using high-throughput small RNA, mRNA and degradome sequencing. The regulation of starch synthesis and quality is carried out mainly at the 15 DAA to 35 DAA stage during grain filling. The DEGs between waxy and non-waxy foxtail millet were significant, especially for GBSS. Additionally, ptc-miR169i_R+2_1ss21GA, fve-miR396e_L-1R+1, mtr-miR162 and PC-5p-221_23413 regulate the expression of genes associated with the starch synthesis pathway in foxtail millet. This study provides new insights into the molecular mechanisms of starch synthesis and quality formation in foxtail millet.

Genome-Wide Identification of NAC Gene Family Members of Tree Peony (Paeonia suffruticosa Andrews) and Their Expression under Heat and Waterlogging Stress.

An important family of transcription factors (TFs) in plants known as NAC (NAM, ATAF1/2, and CUC2) is crucial for the responses of plants to environmental stressors. In this study, we mined the NAC TF family members of tree peony (Paeonia suffruticosa Andrews) from genome-wide data and analyzed their response to heat and waterlogging stresses in conjunction with transcriptome data. Based on tree peony's genomic information, a total of 48 PsNAC genes were discovered. Based on how similar their protein sequences were, these PsNAC genes were divided into 14 branches. While the gene structures and conserved protein motifs of the PsNAC genes within each branch were largely the same, the cis-acting elements in the promoter region varied significantly. Transcriptome data revealed the presence of five PsNAC genes (PsNAC06, PsNAC23, PsNAC38, PsNAC41, PsNAC47) and one PsNAC gene (PsNAC37) in response to heat and waterlogging stresses, respectively. qRT-PCR analysis reconfirmed the response of these five PsNAC genes to heat stress and one PsNAC gene to waterlogging stress. This study lays a foundation for the study of the functions and regulatory mechanisms of NAC TFs in tree peony. Meanwhile, the NAC TFs of tree peony in response to heat and waterlogging stress were excavated, which is of great significance for the selection and breeding of new tree peony varieties with strong heat and waterlogging tolerance.

Exploring Evolutionary Pathways and Abiotic Stress Responses through Genome-Wide Identification and Analysis of the Alternative Oxidase (AOX) Gene Family in Common Oat (Avena sativa).

The alternative oxidase (AOX), a common terminal oxidase in the electron transfer chain (ETC) of plants, plays a crucial role in stress resilience and plant growth and development. Oat (Avena sativa), an important crop with high nutritional value, has not been comprehensively studied regarding the AsAOX gene family. Therefore, this study explored the responses and potential functions of the AsAOX gene family to various abiotic stresses and their potential evolutionary pathways. Additionally, we conducted a genome-wide analysis to explore the evolutionary conservation and divergence of AOX gene families among three Avena species (Avena sativa, Avena insularis, Avena longiglumis) and four Poaceae species (Avena sativa, Oryza sativa, Triticum aestivum, and Brachypodium distachyon). We identified 12 AsAOX, 9 AiAOX, and 4 AlAOX gene family members. Phylogenetic, motif, domain, gene structure, and selective pressure analyses revealed that most AsAOXs, AiAOXs, and AlAOXs are evolutionarily conserved. We also identified 16 AsAOX segmental duplication pairs, suggesting that segmental duplication may have contributed to the expansion of the AsAOX gene family, potentially preserving these genes through subfunctionalization. Chromosome polyploidization, gene structural variations, and gene fragment recombination likely contributed to the evolution and expansion of the AsAOX gene family as well. Additionally, we hypothesize that AsAOX2 may have potential function in resisting wounding and heat stresses, while AsAOX4 could be specifically involved in mitigating wounding stress. AsAOX11 might contribute to resistance against chromium and waterlogging stresses. AsAOX8 may have potential fuction in mitigating ABA-mediated stress. AsAOX12 and AsAOX5 are most likely to have potential function in mitigating salt and drought stresses, respectively. This study elucidates the potential evolutionary pathways of the AsAOXs gene family, explores their responses and potential functions to various abiotic stresses, identifies potential candidate genes for future functional studies, and facilitates molecular breeding applications in A. sativa.

Genome-Wide Identification of NAC Family Genes and Their Expression Analyses in Response to Osmotic Stress in Cannabis sativa L.

NAC (NAM, ATAF1/2, and CUC2) transcription factors are unique and essential for plant growth and development. Although the NAC gene family has been identified in a wide variety of plants, its chromosomal location and function in Cannabis sativa are still unknown. In this study, a total of 69 putative CsNACs were obtained, and chromosomal location analysis indicated that the CsNAC genes mapped unevenly to 10 chromosomes. Phylogenetic analyses showed that the 69 CsNACs could be divided into six subfamilies. Additionally, the CsNAC genes in group IV-a are specific to Cannabis sativa and contain a relatively large number of exons. Promoter analysis revealed that most CsNAC promoters contained cis-elements related to plant hormones, the light response, and abiotic stress. Furthermore, transcriptome expression profiling revealed that 24 CsNAC genes in two Cannabis sativa cultivars (YM1 and YM7) were significantly differentially expressed under osmotic stress, and these 12 genes presented differential expression patterns across different cultivars according to quantitative real-time PCR (RT-qPCR) analysis. Among these, the genes homologous to the CsNAC18, CsNAC24, and CsNAC61 genes have been proven to be involved in the response to abiotic stress and might be candidate genes for further exploration to determine their functions. The present study provides a comprehensive insight into the sequence characteristics, structural properties, evolutionary relationships, and expression patterns of NAC family genes under osmotic stress in Cannabis sativa and provides a basis for further functional characterization of CsNAC genes under osmotic stress to improve agricultural traits in Cannabis sativa.

Pangenome Identification and Analysis of Terpene Synthase Gene Family Members in Gossypium.

Terpene synthases (TPSs), key gatekeepers in the biosynthesis of herbivore-induced terpenes, are pivotal in the diversity of terpene chemotypes across and within plant species. Here, we constructed a gene-based pangenome of the Gossypium genus by integrating the genomes of 17 diploid and 10 tetraploid species. Within this pangenome, 208 TPS syntelog groups (SGs) were identified, comprising 2 core SGs (TPS5 and TPS42) present in all 27 analyzed genomes, 6 softcore SGs (TPS11, TPS12, TPS13, TPS35, TPS37, and TPS47) found in 25 to 26 genomes, 131 dispensable SGs identified in 2 to 24 genomes, and 69 private SGs exclusive to a single genome. The mutational load analysis of these identified TPS genes across 216 cotton accessions revealed a great number of splicing variants and complex splicing patterns. The nonsynonymous/synonymous Ka/Ks value for all 52 analyzed TPS SGs was less than one, indicating that these genes were subject to purifying selection. Of 208 TPS SGs encompassing 1795 genes, 362 genes derived from 102 SGs were identified as atypical and truncated. The structural analysis of TPS genes revealed that gene truncation is a major mechanism contributing to the formation of atypical genes. An integrated analysis of three RNA-seq datasets from cotton plants subjected to herbivore infestation highlighted nine upregulated TPSs, which included six previously characterized TPSs in G. hirsutum (AD1_TPS10, AD1_TPS12, AD1_TPS40, AD1_TPS42, AD1_TPS89, and AD1_TPS104), two private TPSs (AD1_TPS100 and AD2_TPS125), and one atypical TPS (AD2_TPS41). Also, a TPS-associated coexpression module of eight genes involved in the terpenoid biosynthesis pathway was identified in the transcriptomic data of herbivore-infested G. hirsutum. These findings will help us understand the contributions of TPS family members to interspecific terpene chemotypes within Gossypium and offer valuable resources for breeding insect-resistant cotton cultivars.

Genome-wide identification of CaWD40 proteins reveals the involvement of a novel complex (CaAN1-CaDYT1-CaWD40-91) in anthocyanin biosynthesis and genic male sterility in Capsicum annuum.

BACKGROUND: The WD40 domain, one of the most abundant in eukaryotic genomes, is widely involved in plant growth and development, secondary metabolic biosynthesis, and mediating responses to biotic and abiotic stresses. WD40 repeat (WD40) protein has been systematically studied in several model plants but has not been reported in the Capsicum annuum (pepper) genome. RESULTS: Herein, 269, 237, and 257 CaWD40 genes were identified in the Zunla, CM334, and Zhangshugang genomes, respectively. CaWD40 sequences from the Zunla genome were selected for subsequent analysis, including chromosomal localization, phylogenetic relationships, sequence characteristics, motif compositions, and expression profiling. CaWD40 proteins were unevenly distributed on 12 chromosomes, encompassing 19 tandem duplicate gene pairs. The 269 CaWD40s were divided into six main branches (A to F) with 17 different types of domain distribution. The CaWD40 gene family exhibited diverse expression patterns, and several genes were specifically expressed in flowers and seeds. Yeast two-hybrid (Y2H) and dual-luciferase assay indicated that CaWD40-91 could interact with CaAN1 and CaDYT1, suggesting its involvement in anthocyanin biosynthesis and male sterility in pepper. CONCLUSIONS: In summary, we systematically characterized the phylogeny, classification, structure, and expression of the CaWD40 gene family in pepper. Our findings provide a valuable foundation for further functional investigations on WD40 genes in pepper.