Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Biased Retention of Environment-Responsive Genes Following Genome Fractionation

The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.

An ancient whole-genome duplication in barnacles contributes to their diversification and intertidal sessile life adaptation

INTRODUCTION

Whole-genome duplication (WGD) is one of the most sudden and dramatic events rarely reported in invertebrates, but its occurrence can lead to physiological, morphological, and behavioral diversification. WGD has also never been reported in barnacles, which is one of the most unique groups of crustaceans with extremely speciallized morphology (calcareous shells) and habits (intertidal sessile lifestyle).

OBJECTIVES

To investigate whether WGD has occurred in barnacles and examine its potential role in driving the adaptive evolution and diversification of barnacles.

METHODS

Based on a newly sequenced and assembled chromosome-level barnacle genome, a novel WGD event has been identified in barnacles through a comprehensive analysis of interchromosomal synteny, the Hox gene cluster, and synonymous substitution distribution.

RESULTS

We provide ample evidences for WGD in the barnacle genomes. Comparative genomic analysis indicates that this WGD event predates the divergence of Thoracicalcarea, occurring more than 247 million years ago. The retained ohnologs from the WGD are primarily enriched in various pathways related to environmental information processing, shedding light on the adaptive evolution and diversification of intertidal sessile lifestyle. In addition, transcriptomic analyses show that most of these ohnologs were differentially expressed following the ebb of tide. And the cytochrome P450 ohnologs with differential expression patterns are subject to subfunctionalization and/or neofunctionalization for intertidal adaptation. Besides WGD, parallel evolution underlying intertidal adaptation has also occurred in barnacles.

CONCLUSION

This study revealed an ancient WGD event in the barnacle genomes, which is potentially associated with the origin and diversification of thoracican barnacles, and may have contributed to the adaptive evolution of their intertidal sessile lifestyle.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima

PREMISE

The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy.

METHODS

Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium.

RESULTS

Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance.

CONCLUSIONS

These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Characterization, evolution, and abiotic stress responses of leucine-rich repeat receptor-like protein kinases (LRR-RLK) in Liriodendron chinense

BACKGROUND

Liriodendron chinense is susceptible to extinction due to the increasing severity of abiotic stresses resulting from global climate change, consequently impacting its growth, development, and geographic distribution. However, the L. chinense remains pivotal in both socio-economic and ecological realms. The LRR-RLK (leucine-rich repeat receptor-like protein kinase) genes, constituting a substantial cluster of receptor-like kinases in plants, are crucial for plant growth and stress regulation and are unexplored in the L. chinense.

RESULT

233 LchiLRR-RLK genes were discovered, unevenly distributed across 17 chromosomes and 24 contigs. Among these, 67 pairs of paralogous genes demonstrated gene linkages, facilitating the expansion of the LchiLRR-RLK gene family through tandem (35.82%) and segmental (64.18%) duplications. The synonymous and nonsynonymous ratios showed that the LchiLRR-RLK genes underwent a purifying or stabilizing selection during evolution. Investigations in the conserved domain and protein structures revealed that the LchiLRR-RLKs are highly conserved, carrying conserved protein kinase and leucine-rich repeat-like domians that promote clustering in different groups implicating gene evolutionary conservation. A deeper analysis of LchiLRR-RLK full protein sequences phylogeny showed 13 groups with a common ancestor protein. Interspecies gene collinearity showed more orthologous gene pairs between L. chinense and P. trichocarpa, suggesting various similar biological functions between the two plant species. Analysis of the functional roles of the LchiLRR-RLK genes using the qPCR demonstrated that they are involved in cold, heat, and salt stress regulation, especially, members of subgroups VIII, III, and Xa.

CONCLUSION

Conclusively, the LRR-RLK genes are conserved in L. chinense and function to regulate the temperature and salt stresses, and this research provides new insights into understanding LchiLRR-RLK genes and their regulatory effects in abiotic stresses.

Genome-Wide Identification and Characterization of RdHSP Genes Related to High Temperature in Rhododendron delavayi

Heat shock proteins (HSPs) are molecular chaperones that play essential roles in plant development and in response to various environmental stresses. Understanding R. delavayi HSP genes is of great importance since R. delavayi is severely affected by heat stress. In the present study, a total of 76 RdHSP genes were identified in the R. delavayi genome, which were divided into five subfamilies based on molecular weight and domain composition. Analyses of the chromosome distribution, gene structure, and conserved motif of the RdHSP family genes were conducted using bioinformatics analysis methods. Gene duplication analysis showed that 15 and 8 RdHSP genes were obtained and retained from the WGD/segmental duplication and tandem duplication, respectively. Cis-element analysis revealed the importance of RdHSP genes in plant adaptations to the environment. Moreover, the expression patterns of RdHSP family genes were investigated in R. delavayi treated with high temperature based on our RNA-seq data, which were further verified by qRT-PCR. Further analysis revealed that nine candidate genes, including six RdHSP20 subfamily genes (RdHSP20.4, RdHSP20.8, RdHSP20.6, RdHSP20.3, RdHSP20.10, and RdHSP20.15) and three RdHSP70 subfamily genes (RdHSP70.15, RdHSP70.21, and RdHSP70.16), might be involved in enhancing the heat stress tolerance. The subcellular localization of two candidate RdHSP genes (RdHSP20.8 and RdHSP20.6) showed that two candidate RdHSPs were expressed and function in the chloroplast and nucleus, respectively. These results provide a basis for the functional characterization of HSP genes and investigations on the molecular mechanisms of heat stress response in R. delavayi.

Integrative gene duplication and genome-wide analysis as an approach to facilitate wheat reverse genetics: An example in the TaCIPK family

INTRODUCTION

Reverse genetic studies conducted in the plant with a complex or polyploidy genome enriched with large gene families (like wheat) often meet challenges in identifying the key candidate genes related to important traits and prioritizing the genes for functional experiments.

OBJECTIVE

To overcome the above-mentioned challenges of reverse genetics, this work aims to establish an efficient multi-species strategy for genome-wide gene identification and prioritization of the key candidate genes.

METHODS

We established the integrative gene duplication and genome-wide analysis (iGG analysis) as a strategy for pinpointing key candidate genes deserving functional research. The iGG captures the evolution, and the expansion/contraction of large gene families across phylogeny-related species and integrates spatial-temporal expression information for gene function inference. Transgenic approaches were also employed to functional validation.

RESULTS

As a proof-of-concept for the iGG analysis, we took the wheat calcineurin B-like protein-interacting protein kinases (CIPKs) family as an example. We identified CIPKs from seven monocot species, established the orthologous relationship of CIPKs between rice and wheat, and characterized Triticeae-specific CIPK duplicates (e.g., CIPK4 and CIPK17). Integrated with our analysis of CBLs and CBL-CIPK interaction, we revealed that divergent expressions of TaCBLs and TaCIPKs could play an important role in keeping the stoichiometric balance of CBL-CIPK. Furthermore, we validated the function of TaCIPK17-A2 in the regulation of drought tolerance by using transgenic approaches. Overexpression of TaCIPK17 enhanced antioxidant capacity and improved drought tolerance in wheat.

CONCLUSION

The iGG analysis leverages evolutionary and comparative genomics of crops with large genomes to rapidly highlight the duplicated genes potentially associated with speciation, domestication and/or particular traits that deserve reverse-genetic functional studies. Through the identification of Triticeae-specific TaCIPK17 duplicates and functional validation, we demonstrated the effectiveness of the iGG analysis and provided a new target gene for improving drought tolerance in wheat.

Detection of colinear blocks and synteny and evolutionary analyses based on utilization of MCScanX

As different taxa evolve, gene order often changes slowly enough that chromosomal 'blocks' with conserved gene orders (synteny) are discernible. The MCScanX toolkit ( https://github.com/wyp1125/MCScanX ) was published in 2012 as freely available software for the detection of such 'colinear blocks' and subsequent synteny and evolutionary analyses based on genome-wide gene location and protein sequence information. Owing to its simplicity and high efficiency for colinear block detection, MCScanX provides a powerful tool for conducting diverse synteny and evolutionary analyses. Moreover, the detection of colinear blocks has been embraced as an integral step for pangenome graph construction. Here, new application trends of MCScanX are explored, striving to better connect this increasingly used tool to other tools and accelerate insight generation from exponentially growing sequence data. We provide a detailed protocol that covers how to install MCScanX on diverse platforms, tune parameters, prepare input files from data from the National Center for Biotechnology Information, run MCScanX and its visualization and evolutionary analysis tools, and connect MCScanX with external tools, including MCScanX-transposed, Circos and SynVisio. This protocol is easily implemented by users with minimal computational background and is adaptable to new data of interest to them. The data and utility programs for this protocol can be obtained from http://bdx-consulting.com/mcscanx-protocol .

Genome-wide identification and expression analysis of GRAS gene family in Eucalyptus grandis

BACKGROUND

The GRAS gene family is a class of plant-specific transcription factors with important roles in many biological processes, such as signal transduction, disease resistance and stress tolerance, plant growth and development. So far, no information available describes the functions of the GRAS genes in Eucalyptus grandis.

RESULTS

A total of 82 GRAS genes were identified with amino acid lengths ranging from 267 to 817 aa, and most EgrGRAS genes had one exon. Members of the GRAS gene family of Eucalyptus grandis are divided into 9 subfamilies with different protein structures, while members of the same subfamily have similar gene structures and conserved motifs. Moreover, these EgrGRAS genes expanded primarily due to segmental duplication. In addition, cis-acting element analysis showed that this family of genes was involved involved in the signal transduction of various plant hormones, growth and development, and stress response. The qRT-PCR data indicated that 18 EgrGRAS genes significantly responded to hormonal and abiotic stresses. Among them, the expression of EgrGRAS13, EgrGRAS68 and EgrGRAS55 genes was significantly up-regulated during the treatment period, and it was hypothesised that members of the EgrGRAS family play an important role in stress tolerance.

CONCLUSIONS

In this study, the phylogenetic relationship, conserved domains, cis-elements and expression patterns of GRAS gene family of Eucalyptus grandis were analyzed, which filled the gap in the identification of GRAS gene family of Eucalyptus grandis and laid the foundation for analyzing the function of EgrGRAS gene in hormone and stress response.

Long Noncoding RNAs in Response to Hyperosmolarity Stress, but Not Salt Stress, Were Mainly Enriched in the Rice Roots

Due to their immobility and possession of underground parts, plants have evolved various mechanisms to endure and adapt to abiotic stresses such as extreme temperatures, drought, and salinity. However, the contribution of long noncoding RNAs (lncRNAs) to different abiotic stresses and distinct rice seedling parts remains largely uncharacterized beyond the protein-coding gene (PCG) layer. Using transcriptomics and bioinformatics methods, we systematically identified lncRNAs and characterized their expression patterns in the roots and shoots of wild type (WT) and ososca1.1 (reduced hyperosmolality-induced [Ca2+]i increase in rice) seedlings under hyperosmolarity and salt stresses. Here, 2937 candidate lncRNAs were identified in rice seedlings, with intergenic lncRNAs representing the largest category. Although the detectable sequence conservation of lncRNAs was low, we observed that lncRNAs had more orthologs within the Oryza. By comparing WT and ososca1.1, the transcription level of OsOSCA1.1-related lncRNAs in roots was greatly enhanced in the face of hyperosmolality stress. Regarding regulation mode, the co-expression network revealed connections between trans-regulated lncRNAs and their target PCGs related to OsOSCA1.1 and its mediation of hyperosmolality stress sensing. Interestingly, compared to PCGs, the expression of lncRNAs in roots was more sensitive to hyperosmolarity stress than to salt stress. Furthermore, OsOSCA1.1-related hyperosmolarity stress-responsive lncRNAs were enriched in roots, and their potential cis-regulated genes were associated with transcriptional regulation and signaling transduction. Not to be ignored, we identified a motif-conserved and hyperosmolarity stress-activated lncRNA gene (OSlncRNA), speculating on its origin and evolutionary history in Oryza. In summary, we provide a global perspective and a lncRNA resource to understand hyperosmolality stress sensing in rice roots, which helps to decode the complex molecular networks involved in plant sensing and adaptation to stressful environments.

Arogenate dehydratase isoforms strategically deregulate phenylalanine biosynthesis in Akebia trifoliata

Arogenate dehydratase (ADT) is key for phenylalanine (Phe) biosynthesis in plants. To examine ADT components and function in Akebia trifoliata, a representative of Ranunculaceae, we first identified eight ADTs (AktADT1-8, encoding sequences varying from 1032 to 1962 bp) in the A. trifoliata reference genome and five proteins (AktADT1, AktADT4, AktADT7, AktADT8 and AktADT8s) with moonlighting prephenate dehydratase (PDT) activity and Km values varying from 0.43 to 2.17 mM. Structurally, two basic residue combinations (Val314/Ala317 and Ala314/Val317) in the PAC domain are essential for the moonlighting PDT activity of ADTs. Functionally, AktADT4 and AktADT8 successfully restored the wild-type phenotype of pha2, a knockout mutant of Saccharomyces cerevisiae. In addition, AktADTs are ubiquitously expressed, but their expression levels are tissue specific, and the half maximal inhibitory concentration (IC50) of Phe for AktADTs ranged from 49.81 to 331.17 μM. Both AktADT4 and AktADT8 and AktADT8s localized to chloroplast stromules and the cytosol, respectively, while the remaining AktADTs localized to the chloroplast stroma. These findings suggest that various strategies exist for regulating Phe biosynthesis in A. trifoliata. This provides a reasonable explanation for the high Phe content and insights for further genetic improvement of the edible fruits of A. trifoliata.

Inverse and Proportional Trans Modulation of Gene Expression in Human Aneuploidies

Genomic imbalance in aneuploidy is often detrimental to organisms. To gain insight into the molecular basis of aneuploidies in humans, we analyzed transcriptome data from several autosomal and sex chromosome aneuploidies. The results showed that in human aneuploid cells, genes located on unvaried chromosomes are inversely or proportionally trans-modulated, while a subset of genes on the varied chromosomes are compensated. Less genome-wide modulation is found for sex chromosome aneuploidy compared with autosomal aneuploidy due to X inactivation and the retention of dosage sensitive regulators on both sex chromosomes to limit the effective dosage change. We also found that lncRNA and mRNA can have different responses to aneuploidy. Furthermore, we analyzed the relationship between dosage-sensitive transcription factors and their targets, which illustrated the modulations and indicates genomic imbalance is related to stoichiometric changes in components of gene regulatory complexes.In summary, this study demonstrates the existence of trans-acting effects and compensation mechanisms in human aneuploidies and contributes to our understanding of gene expression regulation in unbalanced genomes and disease states.

Identification and Characterization of the Glutathione S-Transferase Gene Family in Blueberry (Vaccinium corymbosum) and Their Potential Roles in Anthocyanin Intracellular Transportation

The glutathione S-transferases (GSTs, EC 2.5.1.18) constitute a versatile enzyme family with pivotal roles in plant stress responses and detoxification processes. Recent discoveries attributed the additional function of facilitating anthocyanin intracellular transportation in plants to GSTs. Our study identified 178 VcGST genes from 12 distinct subfamilies in the blueberry genome. An uneven distribution was observed among these genes across blueberry's chromosomes. Members within the same subfamily displayed homogeneity in gene structure and conserved protein motifs, whereas marked divergence was noted among subfamilies. Functional annotations revealed that VcGSTs were significantly enriched in several gene ontology and KEGG pathway categories. Promoter regions of VcGST genes predominantly contain light-responsive, MYB-binding, and stress-responsive elements. The majority of VcGST genes are subject to purifying selection, with whole-genome duplication or segmental duplication serving as key processes that drive the expansion of the VcGST gene family. Notably, during the ripening of the blueberry fruit, 100 VcGST genes were highly expressed, and the expression patterns of 24 of these genes demonstrated a strong correlation with the dynamic content of fruit anthocyanins. Further analysis identified VcGSTF8, VcGSTF20, and VcGSTF22 as prime candidates of VcGST genes involved in the anthocyanin intracellular transport. This study provides a reference for the exploration of anthocyanin intracellular transport mechanisms and paves the way for investigating the spectrum of GST functions in blueberries.

Identification and characterization of the CRK gene family in the wheat genome and analysis of their expression profile in response to high temperature-induced male sterility

Cysteine-rich receptor-like kinases (CRKs) play many important roles during plant development, including defense responses under both biotic and abiotic stress, reactive oxygen species (ROS) homeostasis, callose deposition and programmed cell death (PCD). However, there are few studies on the involvement of the CRK family in male sterility due to heat stress in wheat (Triticum aestivum L.). In this study, a genome-wide characterization of the CRK family was performed to investigate the structural and functional attributes of the wheat CRKs in anther sterility caused by heat stress. A total of 95 CRK genes were unevenly distributed on 18 chromosomes, with the most genes distributed on chromosome 2B. Paralogous homologous genes with Ka/Ks ratios less than 1 may have undergone strong purifying selection during evolution and are more functionally conserved. The collinearity analysis results of CRK genes showed that wheat and Arabidopsis (A. thaliana), foxtail millet, Brachypodium distachyon (B. distachyon), and rice have three, 12, 15, and 11 pairs of orthologous genes, respectively. In addition, the results of the network interactions of genes and miRNAs showed that five miRNAs were in the hub of the interactions map, namely tae-miR9657b-5p, tae-miR9780, tae-miR9676-5p, tae-miR164, and tae-miR531. Furthermore, qRT-PCR validation of the six TaCRK genes showed that they play key roles in the development of the mononuclear stage anthers, as all six genes were expressed at highly significant levels in heat-stressed male sterile mononuclear stage anthers compared to normal anthers. We hypothesized that the TaCRK gene is significant in the process of high-temperature-induced sterility in wheat based on the combination of anther phenotypes, paraffin sections, and qRT-PCR data. These results improve our understanding of their relationship.

Genome-Wide Identification and Expression Profile Analysis of the Phenylalanine Ammonia-Lyase Gene Family in Hevea brasiliensis

The majority of the world's natural rubber comes from the rubber tree (Hevea brasiliensis). As a key enzyme for synthesizing phenylpropanoid compounds, phenylalanine ammonia-lyase (PAL) has a critical role in plant satisfactory growth and environmental adaptation. To clarify the characteristics of rubber tree PAL family genes, a genome-wide characterization of rubber tree PALs was conducted in this study. Eight PAL genes (HbPAL1-HbPAL8), which spread over chromosomes 3, 7, 8, 10, 12, 13, 14, 16, and 18, were found to be present in the genome of H. brasiliensis. Phylogenetic analysis classified HbPALs into groups I and II, and the group I HbPALs (HbPAL1-HbPAL6) displayed similar conserved motif compositions and gene architectures. Tissue expression patterns of HbPALs quantified by quantitative real-time PCR (qPCR) proved that distinct HbPALs exhibited varying tissue expression patterns. The HbPAL promoters contained a plethora of cis-acting elements that responded to hormones and stress, and the qPCR analysis demonstrated that abiotic stressors like cold, drought, salt, and H2O2-induced oxidative stress, as well as hormones like salicylic acid, abscisic acid, ethylene, and methyl jasmonate, controlled the expression of HbPALs. The majority of HbPALs were also regulated by powdery mildew, anthracnose, and Corynespora leaf fall disease infection. In addition, HbPAL1, HbPAL4, and HbPAL7 were significantly up-regulated in the bark of tapping panel dryness rubber trees relative to that of healthy trees. Our results provide a thorough comprehension of the characteristics of HbPAL genes and set the groundwork for further investigation of the biological functions of HbPALs in rubber trees.

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.)

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Genome organization and botanical diversity

The rich diversity of angiosperms, both the planet's dominant flora and the cornerstone of agriculture, is integrally intertwined with a distinctive evolutionary history. Here, we explore the interplay between angiosperm genome organization and botanical diversity, empowered by genomic approaches ranging from genetic linkage mapping to analysis of gene regulation. Commonality in the genetic hardware of plants has enabled robust comparative genomics that has provided a broad picture of angiosperm evolution and implicated both general processes and specific elements in contributing to botanical diversity. We argue that the hardware of plant genomes-both in content and in dynamics-has been shaped by selection for rather substantial differences in gene regulation between plants and animals such as maize and human, organisms of comparable genome size and gene number. Their distinctive genome content and dynamics may reflect in part the indeterminate development of plants that puts strikingly different demands on gene regulation than in animals. Repeated polyploidization of plant genomes and multiplication of individual genes together with extensive rearrangement and differential retention provide rich raw material for selection of morphological and/or physiological variations conferring fitness in specific niches, whether natural or artificial. These findings exemplify the burgeoning information available to employ in increasing knowledge of plant biology and in modifying selected plants to better meet human needs.

Anthocyanin gene enrichment in the distal region of cotton chromosome A07: mechanisms of reproductive organ coloration

Introduction

The biosynthesis of secondary metabolites like anthocyanins is often governed by metabolic gene clusters (MGCs) in the plant ancestral genome. However, the existence of gene clusters specifically regulating anthocyanin accumulation in certain organs is not well understood.

Methods and results

In this study, we identify MGCs linked to the coloration of cotton reproductive organs, such as petals, spots, and fibers. Through genetic analysis and map-based cloning, we pinpointed key genes on chromosome A07, such as PCC/GhTT19, which is involved in anthocyanin transport, and GbBM and GhTT2-3A, which are associated with the regulation of anthocyanin and proanthocyanidin biosynthesis. Our results demonstrate the coordinated control of anthocyanin and proanthocyanidin pathways, highlighting the evolutionary significance of MGCs in plant adaptation. The conservation of these clusters in cotton chromosome A07 across species underscores their importance in reproductive development and color variation. Our study sheds light on the complex biosynthesis and transport mechanisms for plant pigments, emphasizing the role of transcription factors and transport proteins in pigment accumulation.

Discussion

This research offers insights into the genetic basis of color variation in cotton reproductive organs and the potential of MGCs to enhance our comprehension of plant secondary metabolism.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

Tandemly duplicated MYB genes are functionally diverged in the regulation of anthocyanin biosynthesis in soybean

Gene duplications have long been recognized as a driving force in the evolution of genes, giving rise to novel functions. The soybean (Glycine max) genome is characterized by a large number of duplicated genes. However, the extent and mechanisms of functional divergence among these duplicated genes in soybean remain poorly understood. In this study, we revealed that 4 MYB genes (GmMYBA5, GmMYBA2, GmMYBA1, and Glyma.09g235000)-presumably generated by tandem duplication specifically in the Phaseoleae lineage-exhibited a stronger purifying selection in soybean compared to common bean (Phaseolus vulgaris). To gain insights into the diverse functions of these tandemly duplicated MYB genes in anthocyanin biosynthesis, we examined the expression, transcriptional activity, induced metabolites, and evolutionary history of these 4 MYB genes. Our data revealed that Glyma.09g235000 is a pseudogene, while the remaining 3 MYB genes exhibit strong transcriptional activation activity, promoting anthocyanin biosynthesis in different soybean tissues. GmMYBA5, GmMYBA2, and GmMYBA1 induced anthocyanin accumulation by upregulating the expression of anthocyanin pathway-related genes. Notably, GmMYBA5 showed a lower capacity for gene induction compared to GmMYBA2 and GmMYBA1. Metabolomics analysis further demonstrated that GmMYBA5 induced distinct anthocyanin accumulation in Nicotiana benthamiana leaves and soybean hairy roots compared to GmMYBA2 and GmMYBA1, suggesting their functional divergence leading to the accumulation of different metabolites accumulation following gene duplication. Together, our data provide evidence of functional divergence within the MYB gene cluster following tandem duplication, which sheds light on the potential evolutionary directions of gene duplications during legume evolution.

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid-flavonoid pathway in Thymus mongolicus

BACKGROUND

Thymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.

RESULTS

Chromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid-flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.

CONCLUSION

This study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

Genome-Wide Analysis of the TCP Transcription Factor Gene Family in Pepper (Capsicum annuum L.)

TCP transcription factors play a key role in regulating various developmental processes, particularly in shoot branching, flower development, and leaf development, and these factors are exclusively found in plants. However, comprehensive studies investigating TCP transcription factors in pepper (Capsicum annuum L.) are lacking. In this study, we identified 27 CaTCP members in the pepper genome, which were classified into Class I and Class II through phylogenetic analysis. The motif analysis revealed that CaTCPs in the same class exhibit similar numbers and distributions of motifs. We predicted that 37 previously reported miRNAs target 19 CaTCPs. The expression levels of CaTCPs varied in various tissues and growth stages. Specifically, CaTCP16, a member of Class II (CIN), exhibited significantly high expression in flowers. Class I CaTCPs exhibited high expression levels in leaves, while Class II CaTCPs showed high expression in lateral branches, especially in the CYC/TB1 subclass. The expression profile suggests that CaTCPs play specific roles in the developmental processes of pepper. We provide a theoretical basis that will assist in further functional validation of the CaTCPs.

Genome-Wide Identification and Transcriptome Analysis of the Hsp70 Gene Family in Monodonta labio Reveals Its Role in Response to Nanoplastics Stress

For marine invertebrates, the disruption of organismal physiology and behavior by nanoplastics (NPs) has been extensively reported. Heat shock proteins (Hsps) are important for redundant protein breakdown, environmental changes, and intracellular protein transport. An exhaustive identification of Hsp70 genes and an experiment where different concentrations of NPs were stressed were performed to study how Hsp70 genes respond to NPs stress in Monodonta labio. Our results identified 15 members of Hsp70 within the genome of M. labio and provided insights into their responses to different concentrations of acute NP stress. Phylogenetic analyses revealed extensive amplification of the Hsp70 genes from the Hsc70 subfamily, with gene duplication events. As a result of NP stress, five of fifteen genes showed significant upregulation or downregulation. Three Hsp70 genes were highly expressed at an NP concentration of 0.1 mg/L, and no genes were downregulated. At 10 mg/L, they showed significant upregulation of two genes and significant downregulation of two genes. At 1 mg/L treatment, three genes were significantly downregulated, and no genes were significantly upregulated. Moreover, a purifying selection was revealed using a selection test conducted on duplicate gene pairs, indicating functional redundancy. This work is the first thorough examination of the Hsp70s in Archaeogastropoda. The findings improve knowledge of Hsp70s in molluscan adaptation to NP stress and intertidal living and offer essential data for the biological study of M. labio.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

Polyploid speciation in Zea (Poaceae): cytogenetic insights

MAIN CONCLUSION

The analysis of meiotic pairing affinities and genomic formulae in species and hybrids of Zea allowed us to speculate an evolutionary model to recreate the ancient polyploidization of maize and allied species. The meiotic pairing affinities and the genomic formulae analysis in Zea species and hybrids obtained in new and previous crosses, together with the molecular data known in the genus, allowed us to speculate an evolutionary model to attempt to recreate the ancient polyploidization process of Zea species. We propose that x = 5 semispecies are the ancestors of all modern species of the genus. The complex evolutionary process that originated the different taxa could be included hybridization between sympatric diploid ancestral semispecies (2n = 10) and recurrent duplication of the hybrid chromosome number, resulting in distinct auto- and allopolyploids. After the merger and doubling of independent genomes would have undergone cytological and genetical diploidization, implying revolutionary changes in genome organization and genic balance processes. Based on the meiotic behaviour of the 2n = 30 hybrids, that showed homoeology between the A subgenomes of all parental species, we propose that this subgenome A would be pivotal in all the species and would have conserved the rDNA sequences and the pairing regulator locus (PrZ). In the hypothetical model postulated here, the ancestral semispecies with the pivotal subgenome A would have had a wide geographic distribution, co-occurring and hybridizing with the semispecies harbouring B subgenomes, thus enabling sympatric speciation.

Genome-Wide Identification and Characterization of the TIFY Gene Family and Their Expression Patterns in Response to MeJA and Aluminum Stress in Centipedegrass (Eremochloa ophiuroides)

The TIFY family is a group of novel plant-specific transcription factors involved in plant development, signal transduction, and responses to stress and hormones. TIFY genes have been found and functionally characterized in a number of plant species. However, there is no information about this family in warm-season grass plants. The current study identified 24 TIFY genes in Eremochloa ophiuroides, a well-known perennial warm-season grass species with a high tolerance to aluminum toxicity and good adaptability to the barren acidic soils. All of the 24 EoTIFYs were unevenly located on six out of nine chromosomes and could be classified into two subfamilies (ZIM/ZML and JAZ), consisting of 3 and 21 genes, respectively, with the JAZ subfamily being further divided into five subgroups (JAZ I to JAZ V). The amino acids of 24 EoTIFYs showed apparent differences between the two subfamilies based on the analysis of gene structures and conserved motifs. MCScanX analysis revealed the tandem duplication and segmental duplication of several EoTIFY genes occurred during E. ophiuroides genome evolution. Syntenic analyses of TIFY genes between E. ophiuroides and other five plant species (including A. thaliana, O. sativa, B. distachyon, S. biocolor, and S. italica) provided valuable clues for understanding the potential evolution of the EoTIFY family. qRT-PCR analysis revealed that EoTIFY genes exhibited different spatial expression patterns in different tissues. In addition, the expressions of EoTIFY genes were highly induced by MeJA and all of the EoTIFY family members except for EoJAZ2 displayed upregulated expression by MeJA. Ten EoTIFY genes (EoZML1, EoZML1, EoJAZ1, EoJAZ3, EoJAZ5, EoJAZ6, EoJAZ8, EoJAZ9, EoJAZ10, and EoJAZ21) were observed to be highly expressed under both exogenous MeJA treatment and aluminum stress, respectively. These results suggest that EoTIFY genes play a role in the JA-regulated pathway of plant growth and aluminum resistance as well. The results of this study laid a foundation for further understanding the function of TIFY genes in E. ophiuroides, and provided useful information for future aluminum tolerance related breeding and gene function research in warm-season grass plants.

Genome-wide identification and transcriptome profiling expression analysis of the U-box E3 ubiquitin ligase gene family related to abiotic stress in maize (Zea mays L.)

BACKGROUND

The U-box gene family encodes E3 ubiquitin ligases involved in plant hormone signaling pathways and abiotic stress responses. However, there has yet to be a comprehensive analysis of the U-box gene family in maize (Zea mays L.) and its responses to abiotic stress.

RESULTS

In this study, 85 U-box family proteins were identified in maize and were classified into four subfamilies based on phylogenetic analysis. In addition to the conserved U-box domain, we identified additional functional domains, including Pkinase, ARM, KAP and Tyr domains, by analyzing the conserved motifs and gene structures. Chromosomal localization and collinearity analysis revealed that gene duplications may have contributed to the expansion and evolution of the U-box gene family. GO annotation and KEGG pathway enrichment analysis identified a total of 105 GO terms and 21 KEGG pathways that were notably enriched, including ubiquitin-protein transferase activity, ubiquitin conjugating enzyme activity and ubiquitin-mediated proteolysis pathway. Tissue expression analysis showed that some ZmPUB genes were specifically expressed in certain tissues and that this could be due to their functions. In addition, RNA-seq data for maize seedlings under salt stress revealed 16 stress-inducible plant U-box genes, of which 10 genes were upregulated and 6 genes were downregulated. The qRT-PCR results for genes responding to abiotic stress were consistent with the transcriptome analysis. Among them, ZmPUB13, ZmPUB18, ZmPUB19 and ZmPUB68 were upregulated under all three abiotic stress conditions. Subcellular localization analysis showed that ZmPUB19 and ZmPUB59 were located in the nucleus.

CONCLUSIONS

Overall, our study provides a comprehensive analysis of the U-box gene family in maize and its responses to abiotic stress, suggesting that U-box genes play an important role in the stress response and providing insights into the regulatory mechanisms underlying the response to abiotic stress in maize.

Investigation of regulatory divergence between homoeologs in the recently formed allopolyploids, Tragopogon mirus and T. miscellus (Asteraceae)

Polyploidy is an important evolutionary process throughout eukaryotes, particularly in flowering plants. Duplicated gene pairs (homoeologs) in allopolyploids provide additional genetic resources for changes in molecular, biochemical, and physiological mechanisms that result in evolutionary novelty. Therefore, understanding how divergent genomes and their regulatory networks reconcile is vital for unraveling the role of polyploidy in plant evolution. Here, we compared the leaf transcriptomes of recently formed natural allotetraploids (Tragopogon mirus and T. miscellus) and their diploid parents (T. porrifolius X T. dubius and T. pratensis X T. dubius, respectively). Analysis of 35 400 expressed loci showed a significantly higher level of transcriptomic additivity compared to old polyploids; only 22% were non-additively expressed in the polyploids, with 5.9% exhibiting transgressive expression (lower or higher expression in the polyploids than in the diploid parents). Among approximately 7400 common orthologous regions (COREs), most loci in both allopolyploids exhibited expression patterns that were vertically inherited from their diploid parents. However, 18% and 20.3% of the loci showed novel expression bias patterns in T. mirus and T. miscellus, respectively. The expression changes of 1500 COREs were explained by cis-regulatory divergence (the condition in which the two parental subgenomes do not interact) between the diploid parents, whereas only about 423 and 461 of the gene expression changes represent trans-effects (the two parental subgenomes interact) in T. mirus and T. miscellus, respectively. The low degree of both non-additivity and trans-effects on gene expression may present the ongoing evolutionary processes of the newly formed Tragopogon polyploids (~80-90 years).

Genome-wide identification of B-box zinc finger (BBX) gene family in Medicago sativa and their roles in abiotic stress responses

BACKGROUND

B-box (BBX) family is a class of zinc finger transcription factors (TFs) that play essential roles in regulating plant growth, development, as well as abiotic stress. However, no systematic analysis of BBX genes has yet been conducted in alfalfa (Medica go sativa L.), and their functions have not been elucidated up to now.

RESULTS

In this study, 28 MsBBX genes were identified from the alfalfa genome, which were clustered into 4 subfamilies according to an evolutionary tree of BBX proteins. Exon-intron structure and conserved motif analysis reflected the evolutionary conservation of MsBBXs in alfalfa. Collinearity analysis showed that segmental duplication promoted the expansion of the MsBBX family. Analysis of cis-regulatory elements suggested that the MsBBX genes possessed many growth/development-, light-, phytohormone-, and abiotic stress-related elements. MsBBX genes were differentially expressed in leaves, flowers, pre-elongated stems, elongated stems, roots and nodules, and most MsBBXs were remarkably induced by drought, salt and various plant growth regulators (ABA, JA, and SA). Further functional verification demonstrated that overexpressing of the MsBBX11 gene clearly promoted salt tolerance in transgenic Arabidopsis by regulating growth and physiological processes of seedlings.

CONCLUSIONS

This research provides insights into further functional research and regulatory mechanisms of MsBBX family genes under abiotic stress of alfalfa.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Gene expression bias between the subgenomes of allopolyploid hybrids is an emergent property of the kinetics of expression

Hybridization coupled to polyploidy, or allopolyploidy, has dramatically shaped the evolution of flowering plants, teleost fishes, and other lineages. Studies of recently formed allopolyploid plants have shown that the two subgenomes that merged to form that new allopolyploid do not generally express their genes equally. Instead, one of the two subgenomes expresses its paralogs more highly on average. Meanwhile, older allopolyploidy events tend to show biases in duplicate losses, with one of the two subgenomes retaining more genes than the other. Since reduced expression is a pathway to duplicate loss, understanding the origins of expression biases may help explain the origins of biased losses. Because we expect gene expression levels to experience stabilizing selection, our conceptual frameworks for how allopolyploid organisms form tend to assume that the new allopolyploid will show balanced expression between its subgenomes. It is then necessary to invoke phenomena such as differences in the suppression of repetitive elements to explain the observed expression imbalances. Here we show that, even for phenotypically identical diploid progenitors, the inherent kinetics of gene expression give rise to biases between the expression levels of the progenitor genes in the hybrid. Some of these biases are expected to be gene-specific and not give rise to global differences in progenitor gene expression. However, particularly in the case of allopolyploids formed from progenitors with different genome sizes, global expression biases favoring one subgenome are expected immediately on formation. Hence, expression biases are arguably the expectation upon allopolyploid formation rather than a phenomenon needing explanation. In the future, a deeper understanding of the kinetics of allopolyploidy may allow us to better understand both biases in duplicate losses and hybrid vigor.

Genome-wide identification and expression analysis of the glutamate receptor gene family in sweet potato and its two diploid relatives

Plant glutamate receptor (GLR) homologs are crucial calcium channels that play an important role in plant development, signal transduction, and response to biotic and abiotic stresses. However, the GLR gene family has not yet been thoroughly and systematically studied in sweet potato. In this study, a total of 37 GLR genes were identified in the cultivated hexaploid sweet potato (Ipomoea batatas), and 32 GLR genes were discovered in each of the two diploid relatives (Ipomoea trifida and Ipomoea triloba) for the first time. Based on their evolutionary relationships to those of Arabidopsis, these GLRs were split into five subgroups. We then conducted comprehensive analysis to explore their physiological properties, protein interaction networks, promoter cis-elements, chromosomal placement, gene structure, and expression patterns. The results indicate that the homologous GLRs of the cultivated hexaploid sweet potato and its two relatives are different. These variations are reflected in their functions related to plant growth, hormonal crosstalk, development of tuberous roots, resistance to root rot, and responses to abiotic stress factors, all of which are governed by specific individual GLR genes. This study offers a comprehensive analysis of GLR genes in sweet potato and its two diploid relatives. It also provides a theoretical basis for future research into their regulatory mechanisms, significantly influencing the field of molecular breeding in sweet potatoes.

Comprehensive Analysis of the Catalase (CAT) Gene Family and Expression Patterns in Rubber Tree (Hevea brasiliensis) under Various Abiotic Stresses and Multiple Hormone Treatments

Catalase (CAT) is one of the key enzymes involved in antioxidant defense systems and mainly scavenges H2O2 and plays a vital role in plant growth, development, and various adverse stresses. To date, a systematic study of the CAT gene family in rubber tree has not been reported. In this study, five HbCAT gene family members were identified from the rubber tree genome, and these were mainly clustered into two subfamilies. Gene structure and motif analysis showed that exon-intron and motif patterns were conserved across different plant species. Sequence analysis revealed that HbCAT proteins contain one active catalytic site, one heme-ligand signature sequence, three conserved amino acid residues (His, Tyr, and Asn), and one peroxisome-targeting signal 1 (PTS1) sequence. Fragment duplication is a selection pressure for the evolution of the HbCAT family based on Ka/Ks values. Analysis of cis-acting elements in the promoters indicated that HbCAT gene expression might be regulated by abscisic acid (ABA), salicylic acid (SA), and MYB transcription factors; furthermore, these genes might be involved in plant growth, development, and abiotic stress responses. A tissue-specific expression analysis showed that HbCATs gradually increased with leaf development and were highly expressed in mature leaves. Gene expression profiling exhibited the differential expression of the HbCATs under cold, heat, drought, and NaCl stresses. Our results provide comprehensive information about the HbCAT gene family, laying the foundation for further research on its function in rubber tree.

Tandem gene duplications contributed to high-level azole resistance in a rapidly expanding Candida tropicalis population

Invasive diseases caused by the globally distributed commensal yeast Candida tropicalis are associated with mortality rates of greater than 50%. Notable increases of azole resistance have been observed in this species, particularly within Asia-Pacific regions. Here, we carried out a genetic population study on 1571 global C. tropicalis isolates using multilocus sequence typing (MLST). In addition, whole-genome sequencing (WGS) analysis was conducted on 629 of these strains, comprising 448 clinical invasive strains obtained in this study and 181 genomes sourced from public databases. We found that MLST clade 4 is the predominant azole-resistant clone. WGS analyses demonstrated that dramatically increasing rates of azole resistance are associated with a rapid expansion of cluster AZR, a sublineage of clade 4. Cluster AZR isolates exhibited a distinct high-level azole resistance, which was induced by tandem duplications of the ERG11A395T gene allele. Ty3/gypsy-like retrotransposons were found to be highly enriched in this population. The alarming expansion of C. tropicalis cluster AZR population underscores the urgent need for strategies against growing threats of antifungal resistance.

Characteristics and expression of lncRNA and transposable elements in Drosophila aneuploidy

Aneuploidy can globally affect the expression of the whole genome, which is detrimental to organisms. Dosage-sensitive regulators usually have multiple intermolecular interactions, and changes in their stoichiometry are responsible for the dysregulation of the regulatory network. Currently, studies on noncoding genes in aneuploidy are relatively rare. We studied the characteristics and expression profiles of long noncoding RNAs (lncRNAs) and transposable elements (TEs) in aneuploid Drosophila. It is found that lncRNAs and TEs are affected by genomic imbalance and appear to be more sensitive to an inverse dosage effect than mRNAs. Several dosage-sensitive lncRNAs and TEs were detected for their expression patterns during embryogenesis, and their biological functions in the ovary and testes were investigated using tissue-specific RNAi. This study advances our understanding of the noncoding sequences in imbalanced genomes and provides a novel perspective for the study of aneuploidy-related human diseases such as cancer.

Identification and analysis of BAHD superfamily related to malonyl ginsenoside biosynthesis in Panax ginseng

Introduction

The BAHD (benzylalcohol O-acetyl transferase, anthocyanin O-hydroxycinnamoyl transferase, N-hydroxycinnamoyl anthranilate benzoyl transferase and deacetylvindoline 4-O-acetyltransferase), has various biological functions in plants, including catalyzing the biosynthesis of terpenes, phenolics and esters, participating in plant stress response, affecting cell stability, and regulating fruit quality.

Methods

Bioinformatics methods, real-time fluorescence quantitative PCR technology, and ultra-high-performance liquid chromatography combined with an Orbitrap mass spectrometer were used to explore the relationship between the BAHD gene family and malonyl ginsenosides in Panax ginseng.

Results

In this study, 103 BAHD genes were identified in P. ginseng, mainly distributed in three major clades. Most PgBAHDs contain cis-acting elements associated with abiotic stress response and plant hormone response. Among the 103 genes, 68 PgBAHDs are WGD (whole-genome duplication) genes. The significance of malonylation in biosynthesis has garnered considerable attention in the study of malonyltransferases. The phylogenetic tree results showed 34 PgBAHDs were clustered with genes that have malonyl characterization. Among them, seven PgBAHDs (PgBAHD4, 45, 65, 74, 90, 97, and 99) showed correlations > 0.9 with crucial enzyme genes involved in ginsenoside biosynthesis and > 0.8 with malonyl ginsenosides. These seven genes were considered potential candidates involved in the biosynthesis of malonyl ginsenosides.

Discussion

These results help elucidate the structure, evolution, and functions of the P. ginseng BAHD gene family, and establish the foundation for further research on the mechanism of BAHD genes in ginsenoside biosynthesis.

Transposition, duplication, and divergence of the telomerase RNA underlies the evolution of Mimulus telomeres

Telomeres are nucleoprotein complexes with a crucial role of protecting chromosome ends. It consists of simple repeat sequences and dedicated telomere-binding proteins. Because of its vital functions, components of the telomere, for example its sequence, should be under strong evolutionary constraint. But across all plants, telomere sequences display a range of variation and the evolutionary mechanism driving this diversification is largely unknown. Here, we discovered in Monkeyflower (Mimulus) the telomere sequence is even variable between species. We investigated the basis of Mimulus telomere sequence evolution by studying the long noncoding telomerase RNA (TR), which is a core component of the telomere maintenance complex and determines the telomere sequence. We conducted total RNA-based de novo transcriptomics from 16 Mimulus species and analyzed reference genomes from 6 species, and discovered Mimulus species have evolved at least three different telomere sequences: (AAACCCT)n, (AAACCCG)n, and (AAACCG)n. Unexpectedly, we discovered several species with TR duplications and the paralogs had functional consequences that could influence telomere evolution. For instance, M. lewisii had two sequence-divergent TR paralogs and synthesized a telomere with sequence heterogeneity, consisting of AAACCG and AAACCCG repeats. Evolutionary analysis of the M. lewisii TR paralogs indicated it had arisen from a transposition-mediate duplication process. Further analysis of the TR from multiple Mimulus species showed the gene had frequently transposed and inserted into new chromosomal positions during Mimulus evolution. From our results, we propose the TR transposition, duplication, and divergence model to explain the evolutionary sequence turnovers in Mimulus and potentially all plant telomeres.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Genomes of Endotrypanum monterogeii from Panama and Zelonia costaricensis from Brazil: Expansion of Multigene Families in Leishmaniinae Parasites That Are Close Relatives of Leishmania spp

The Leishmaniinae subfamily of the Trypanosomatidae contains both genus Zelonia (monoxenous) and Endotrypanum (dixenous). They are amongst the nearest known relatives of Leishmania, which comprises many human pathogens widespread in the developing world. These closely related lineages are models for the genomic biology of monoxenous and dixenous parasites. Herein, we used comparative genomics to identify the orthologous groups (OGs) shared among 26 Leishmaniinae species to investigate gene family expansion/contraction and applied two phylogenomic approaches to confirm relationships within the subfamily. The Endotrypanum monterogeii and Zelonia costaricensis genomes were assembled, with sizes of 29.9 Mb and 38.0 Mb and 9.711 and 12.201 predicted protein-coding genes, respectively. The genome of E. monterogeii displayed a higher number of multicopy cell surface protein families, including glycoprotein 63 and glycoprotein 46, compared to Leishmania spp. The genome of Z. costaricensis presents expansions of BT1 and amino acid transporters and proteins containing leucine-rich repeat domains, as well as a loss of ABC-type transporters. In total, 415 and 85 lineage-specific OGs were identified in Z. costaricensis and E. monterogeii. The evolutionary relationships within the subfamily were confirmed using the supermatrix (3384 protein-coding genes) and supertree methods. Overall, this study showed new expansions of multigene families in monoxenous and dixenous parasites of the subfamily Leishmaniinae.

Systematic analysis of the Rboh gene family in seven gramineous plants and its roles in response to arbuscular mycorrhizal fungi in maize

BACKGROUND

Plant respiratory burst oxidase homolog (Rboh) gene family produces reactive oxygen species (ROS), and it plays key roles in plant-microbe interaction. Most Rboh gene family-related studies mainly focused on dicotyledonous plants; however, little is known about the roles of Rboh genes in gramineae.

RESULTS

A total of 106 Rboh genes were identified in seven gramineae species, including Zea mays, Sorghum bicolor, Brachypodium distachyon, Oryza sativa, Setaria italica, Hordeum vulgare, and Triticum aestivum. The Rboh protein sequences showed high similarities, suggesting that they may have conserved functions across different species. Duplication mode analysis detected whole-genome/segmental duplication (WGD)/(SD) and dispersed in the seven species. Interestingly, two local duplication (LD, including tandem and proximal duplication) modes were found in Z. mays, S. italica and H. vulgare, while four LD were detected in T. aestivum, indicating that these genes may have similar functions. Collinearity analysis indicated that Rboh genes are at a stable evolution state in all the seven species. Besides, Rboh genes from Z. mays were closely related to those from S. bicolor, consistent with the current understanding of plant evolutionary history. Phylogenetic analysis showed that the genes in the subgroups I and II may participate in plant-AM fungus symbiosis. Cis-element analysis showed that different numbers of elements are related to fungal induction in the promoter region. Expression profiles of Rboh genes in Z. mays suggested that Rboh genes had distinct spatial expression patterns. By inoculation with AM fungi, our transcriptome analysis showed that the expression of Rboh genes varies upon AM fungal inoculation. In particularly, ZmRbohF was significantly upregulated after inoculation with AM fungi. pZmRbohF::GUS expression analyses indicated that ZmRbohF was induced by arbuscular mycorrhizal fungi in maize. By comparing WT and ZmRbohF mutant, we found ZmRbohF had limited impact on the establishment of maize-AM fungi symbiosis, but play critical roles in regulating the proper development of arbuscules.

CONCLUSIONS

This study provides a comprehensive analysis of the evolution relationship of Rboh genes in seven gramineae species. Results showed that several Rboh genes regulate maize-AM fungal symbiosis process. This study provides valuable information for further studies of Rboh genes in gramineae.

Comprehensive characterization and expression analysis of enzymatic antioxidant gene families in passion fruit (Passiflora edulis)

Passion fruit, a valuable tropical fruit, faces climate-related growth challenges. Antioxidant enzymes are vital for both stress protection and growth regulation in plants. We first provided systemic analysis of enzymatic antioxidant gene families in passion fruit, identifying 90 members including 11 PeSODs, 45 PeAPXs, 8 PeCATs, 7 PeGPXs, 6 PeMDHARs, 8 PeDHARs, and 5 PeGRs. Gene members in each gene family with same subcellular localization showed closer phylogenetic relationship. Many antioxidant genes exhibited tissue- or developmental stage-specific expression patterns during floral and fruit development, with some widely expressed. Their co-expressed genes were linked to photosynthesis and energy metabolism, suggesting roles in protecting highly proliferating tissues from oxidative damage. Potential genes for enhancing temperature stress resistance were identified. The involvement of diverse regulatory factors including miRNAs, transcription factors, and CREs might contribute to the complex roles of antioxidant genes. This study informs future research on antioxidant genes and passion fruit breeding.

VDR promotes testosterone synthesis in mouse Leydig cells via regulation of cholesterol side chain cleavage cytochrome P450 (Cyp11a1) expression

BACKGROUND

The vitamin D receptor (VDR) mediates the pleiotropic biological actions that include osteoporosis, immune responses and androgen synthesis.VDR is widely expressed in testis cells such as Leydig cells, Sertoli cells, and sperm. The levels of steroids are critical for sexual development. In the early stage of steroidogenesis, cholesterol is converted to pregnenolone (precursor of most steroid hormones) by cholesterol side-chain lyase (CYP11A1), which eventually synthesizes the male hormone testosterone.

OBJECTIVE

This study aims to reveal how VDR regulates CYP11A1 expression and affects testosterone synthesis in murine Leydig cells.

METHODS

The levels of VDR, CYP11A1 were determined by quantitative real-time polymerase chain reaction (RT-qPCR) or western blot. Targeted relationship between VDR and Cyp11a1 was evaluated by dual-luciferase reporter assay. The levels of testosterone concentrations in cell culture media serum by enzyme-linked immunosorbent assay (ELISA).

RESULTS

Phylogenetic and motif analysis showed that the Cyp11a1 family had sequence loss, which may have special biological functions during evolution. The results of promoter prediction showed that vitamin D response element (VDRE) existed in the upstream promoter region of murine Cyp11a1. Dual-luciferase assay confirmed that VDR could bind candidate VDREs in upstream region of Cyp11a1, and enhance gene expression. Tissue distribution and localizatio analysis showed that Cyp11a1 was mainly expressed in testis, and dominantly existed in murine Leydig cells. Furthermore, over-expression VDR and CYP11A1 significantly increased testosterone synthesis in mice Leydig cells.

CONCLUSIONS

Active vitamin D3 (VD3) and Vdr interference treatment showed that VD3/VDR had a positive regulatory effect on Cyp11a1 expression and testosterone secretion. VDR promotes testosterone synthesis in male mice by up-regulating Cyp11a1 expression, which played an important role for male reproduction.

Genome-wide identification and expression analysis of the BURP domain-containing genes in Malus domestica

The conserved BURP-containing proteins are specific to plants and play a crucial role in plant growth, development, and response to abiotic stresses. However, less is known about the systematic characterization of BURP-containing proteins in apple. This study aimed to identify and analyze all BURP-containing genes in the apple genome, as well as to examine their expression patterns through various bioinformatics methods. Eighteen members of BURP-containing genes were identified in apple, six members lacked signal peptides, and the secondary structure was mainly a Random coil of BURP-containing genes. Gene structure and Motif analysis showed that proteins have similar structures and are conserved at the C-terminal. Cis-acting element analysis revealed that the proteins contain phytohormone and stress response elements, and chromosomal localization revealed that the family is unevenly distributed across eight chromosomes, with duplication of fragments leading to the expansion of family proteins. Tissue expression showed that MdPG3 and MdPG4 were expressed in different tissues and different varieties, MdRD2 and MdRD7 were highly expressed in 'M74' fruits and MdRD7 in 'M49' leaves, while MdUSP1 was highly expressed in 'GD' roots. The quantitative real-time PCR analysis showed that the expressions of six and seven genes were significantly up-regulated under NaCl and PEG treatments, respectively, whereas MdRD7 was significantly up-regulated under NaCl and PEG treatment over time. This study offers a comprehensive identification and expression analysis of BURP-containing proteins in apple. The findings provide a theoretical foundation for further exploration of the functions of this protein family.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01393-7.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

Complex origin, evolution, and diversification of non-canonically organized OVATE-OFP and OVATE-Like OFP gene pair across Embryophyta

Ovate Family Proteins (OFP) is a plant-specific gene family of negative transcriptional regulators. Till-date, a handful of in-silico studies have provided glimpses into family size, expansion patterns, and genic features across all major plant lineages. A major lacuna exists in understanding origin of organisation complexity of members such as those arranged in a head-to-head manner which may lead to transcriptional co-regulation via a common bi-directional promoter. To address this gap, we investigated the origin, organization and evolution of two head-to-head arranged gene pairs of homologs of AtOFP2-AtOFP17, and, AtOFP4-AtOFP20 across Archaeplastida. The ancestral forms of AtOFP2, AtOFP4, AtOFP17, and AtOFP20 are likely to have evolved in last common ancestors of Embryophyta (land plants) given their complete absence in Rhodophyta and Chlorophyta. The OFP gene family originated and expanded in Bryophyta, including protein variants with complete (OVATE-OFP) or partial (OVATE-Like OFP) OVATE domain; with head-to-head organization present only in Spermatophyta (gymnosperms and angiosperms). Ancestral State Reconstruction revealed the origin of head-to-head organized gene pair in gymnosperms, with both genes being OVATE-OFP (homologs of AtOFP2/4). Phylogenetic reconstruction and copy number analysis suggests the presence of a single copy of the head-to-head arranged pair of OFP2/4 (OVATE)-OFP17/20 (OVATE-Like) in all angiosperms except Brassicaceae, and a duplication event in last common ancestor of core Brassicaceae approximately 32-54 MYA leading to origin of AtOFP2-AtOFP17 and AtOFP4-AtOFP20 as paralogs. Synteny analysis of genomic regions harbouring homologs of AtOFP2-AtOFP17, AtOFP4-AtOFP20 and AtOFP2/4-AtOFP17/20 across angiosperms suggested ancestral nature of AtOFP2-AtOFP17 gene pair. The present study thus establishes the orthology and evolutionary history of two non-canonically organised gene pairs with variation in their OVATE domain. The non-canonical organisation, atleast in Brassicaceae, has the potential of generating complex transcriptional regulation mediated via a common bi-directional promoter. The study thus lays down a framework to understand evolution of gene and protein structure, transcriptional regulation and function across a phylogenetic lineage through comparative analyses.

Superoxide dismutase gene family in cassava revealed their involvement in environmental stress via genome-wide analysis

Superoxide dismutase (SOD) is a crucial metal-containing enzyme that plays a vital role in catalyzing the dismutation of superoxide anions, converting them into molecular oxygen and hydrogen peroxide, essential for enhancing plant stress tolerance. We identified 8 SOD genes (4 CSODs, 2 FSODs, and 2 MSODs) in cassava. Bioinformatics analyses provided insights into chromosomal location, phylogenetic relationships, gene structure, conserved motifs, and gene ontology annotations. MeSOD genes were classified into two groups through phylogenetic analysis, revealing evolutionary connections. Promoters of these genes harbored stress-related cis-elements. Duplication analysis indicated the functional significance of MeCSOD2/MeCSOD4 and MeMSOD1/MeMSOD2. Through qRT-PCR, MeCSOD2 responded to salt stress, MeMSOD2 to drought, and cassava bacterial blight. Silencing MeMSOD2 increased XpmCHN11 virulence, indicating MeMSOD2 is essential for cassava's defense against XpmCHN11 infection. These findings enhance our understanding of the SOD gene family's role in cassava and contribute to strategies for stress tolerance improvement.

Chromosome-level genome assembly of Niphotrichum japonicum provides new insights into heat stress responses in mosses

With a diversity of approximately 22,000 species, bryophytes (hornworts, liverworts, and mosses) represent a major and diverse lineage of land plants. Bryophytes can thrive in many extreme environments as they can endure the stresses of drought, heat, and cold. The moss Niphotrichum japonicum (Grimmiaceae, Grimmiales) can subsist for extended periods under heat and drought conditions, providing a good candidate for studying the genetic basis underlying such high resilience. Here, we de novo assembled the genome of N. japonicum using Nanopore long reads combined with Hi-C scaffolding technology to anchor the 191.61 Mb assembly into 14 pseudochromosomes. The genome structure of N. japonicum's autosomes is mostly conserved and highly syntenic, in contrast to the sparse and disordered genes present in its sex chromosome. Comparative genomic analysis revealed the presence of 10,019 genes exclusively in N. japonicum. These genes may contribute to the species-specific resilience, as demonstrated by the gene ontology (GO) enrichment. Transcriptome analysis showed that 37.44% (including 3,107 unique genes) of the total annotated genes (26,898) exhibited differential expression as a result of heat-induced stress, and the mechanisms that respond to heat stress are generally conserved across plants. These include the upregulation of HSPs, LEAs, and reactive oxygen species (ROS) scavenging genes, and the downregulation of PPR genes. N. japonicum also appears to have distinctive thermal mechanisms, including species-specific expansion and upregulation of the Self-incomp_S1 gene family, functional divergence of duplicated genes, structural clusters of upregulated genes, and expression piggybacking of hub genes. Overall, our study highlights both shared and species-specific heat tolerance strategies in N. japonicum, providing valuable insights into the heat tolerance mechanism and the evolution of resilient plants.

Genome-wide identification of GH9 gene family and the assessment of its role during fruit abscission zone formation in Vaccinium ashei

KEY MESSAGE

The genomic location and stage-specific expression pattern of GH9 genes reveal their critical roles during fruit abscission zone formation in Vaccinium ashei. Glycosyl hydrolase family 9 (GH9) cellulases play a crucial role in both cellulose synthesis and hydrolysis during plant growth and development. Despite this importance, there is currently no study on the involvement of GH9-encoding genes, specifically VaGH9s, in abscission zone formation of rabbiteye blueberries (Vaccinium ashei). In this study, we identified a total of 61 VaGH9s in the genome, which can be classified into 3 subclasses based on conserved motifs and domains, gene structures, and phylogenetic analyses. Our synteny analysis revealed that VaGH9s are more closely related to the GH9s of Populus L. than to those of Arabidopsis, Vitis vinifera, and Citrus sinensis. In silico structural analysis predicted that most of VaGH9s are hydrophilic, and localized in cell membrane and/or cell wall, and the variable sets of cis-acting regulatory elements and functional diversity with four categories of stress response, hormone regulation, growth and development, and transcription factor-related elements are present in the promoter sequence of VaGH9s genes. Transcriptomic analysis showed that there were 22 differentially expressed VaGH9s in fruit abscission zone tissue at the veraison stage, and the expression of VaGH9B2 and VaGH9C10 was continuously increased during fruit maturation, which were in parallel with the increasing levels of cellulase activity and oxidative stress indicators, suggesting that they are involved in the separation stage of fruit abscission in Vaccinium ashei. Our work identified 22 VaGH9s potentially involved in different stages of fruit abscission and would aid further investigation into the molecular regulation of abscission in rabbiteye blueberries fruit.