Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants

Terpene synthases GhTPS6 and GhTPS47 participate in resistance to Verticillium dahliae in upland cotton

Terpene synthases (TPSs) are enzymes responsible for catalyzing the production of diverse terpenes, the largest class of secondary metabolites in plants. Here, we identified 107 TPS gene loci encompassing 92 full-length TPS genes in upland cotton (Gossypium hirsutum L.). Phylogenetic analysis showed they were divided into six subfamilies. Segmental duplication and tandem duplication events contributed greatly to the expansion of TPS gene family, particularly the TPS-a and TPS-b subfamilies. Expression profile analysis screened out that GhTPSs may mediate the interaction between cotton and Verticillium dahliae. Three-dimensional structures and subcellular localizations of the two selected GhTPSs, GhTPS6 and GhTPS47, which belong to the TPS-a subfamily, demonstrated similarity in protein structures and nucleus and cytoplasm localization. Virus-induced gene silencing (VIGS) of the two GhTPSs yielded plants characterized by increased wilting and chlorosis, more severe vascular browning, and higher disease index than control plants. Additionally, knockdown of GhTPS6 and GhTPS47 led to the down-regulation of cotton terpene synthesis following V. dahliae infection, indicating that these two genes may positively regulate resistance to V. dahliae through the modulation of disease-resistant terpene biosynthesis. Overall, our study represents a comprehensive analysis of the G. hirsutum TPS gene family, revealing their potential roles in defense responses against Verticillium wilt.

Phyllostachys edulis argonaute genes function in the shoot architecture

Argonaute (AGO) proteins are the core components of the RNA-induced silencing complexes (RISC) in the cytoplasm and nucleus, and are necessary for the development of plant shoot meristem, which gives rise to the above-ground plant body. In this study, we identified 23 Phyllostachys edulis AGO genes (PhAGOs) that were distributed unequally on the 14 unmapped scaffolds. Gene collinearity and phylogeny analysis showed that the innovation of PhAGO genes was mainly due to dispersed duplication and whole-genome duplication, which resulted in the enlarged PhAGO family. PhAGO genes were expressed in a temporal-spatial expression pattern, and they encoded proteins differently localized in the cytoplasm and/or nucleus. Overexpression of the PhAGO2 and PhAGO4 genes increased the number of tillers or leaves in Oryza sativa and affected the shoot architecture of Arabidopsis thaliana. These results provided insight into the fact that PhAGO genes play important roles in plant development.

The chromosome-level genome assembly of Cananga odorata provides insights into its evolution and terpenoid biosynthesis

Cananga odorata is known as a natural perfume tree of the Annonaceae family in Magnoliales. However, its phylogenetic position and the molecular mechanisms involved in the biosynthesis of the floral volatile organic compounds (VOCs) remain unclear. Here, by combining a variety of sequencing platforms, we present a telomere-to-telomere (T2T) genome of C. odorata with 735.83 Mb, which represents the highest integrity and assembly quality of genome in magnoliid plants reported to date. Phylogenetic analysis based on multiple datasets and approaches showed that C. odorata, as a member of magnoliids, is sister to eudicots, after their divergence from monocots. Metabolomic of VOCs in the essential oil and flowers scent showed that sesquiterpenes, especially β-caryophyllene, were the major compounds. Two CoTPS21 homologues derived from tandem duplication events were highly expressed during flower development and were identified as the key sesquiterpene synthases for the production of β-caryophyllene. In addition, CoSPL3 and CoSPL9 were considered as potential transcription factors for activating the expression of CoTPS21 homologues. Our results shed light on the molecular mechanisms underlying the biosynthesis of the unique floral fragrance in C. odorata and provide new insights into the phylogenetic position of magnoliids.

Genome-wide identification of R2R3-MYB family in Eriobotrya japonica and functional analysis of EjMYB5 involved in proanthocyanidin biosynthesis

Loquat (Eriobotrya japonica Lindl.) is a popular fruit and medicinal plant. Proanthocyanidins (PAs), as one of the main types of flavonoids, are the key components of loquat fruit quality and medicinal properties. However, the identification of transcription factors (TFs) involved in PA accumulation in loquat remains limited. R2R3-MYB TFs play key regulatory role in PA accumulation in plants. In this study, 190 R2R3-MYB TFs were identified in loquat genome. Combined with transcriptome data, R2R3-MYB TF EjMYB5 involved in PA accumulation in loquat was isolated. EjMYB5 was transcriptional activator localized to nucleus. Expression of EjMYB5 was closely related to PA accumulation in loquat fruits. Heterogenous overexpression of EjMYB5 in tomato (Solanum lycopersicum) inhibited anthocyanin accumulation and promoted PA accumulation. Additionally, transient overexpression of EjMYB5 in tobacco (Nicotiana benthamiana) leaves promoted PA accumulation by upregulating flavonoid biosynthesis genes (NtDFR, NtANS, and NtLAR). Transcriptome analysis of EjMYB5-overexpressing tomato fruits suggested that EjMYB5 was involved in several biological pathways, including lipid metabolism, MAPK signaling, phenylpropanoid biosynthesis, and flavonoid biosynthesis. Collectively, our findings provided basic data for further analysis the function of R2R3-MYB TFs in loquat, and revealed that EjMYB5 functioned as PA accumulation in loquat.

Comprehensive genome-wide identification of Snf2 gene family and their expression profile under salt stress in six Brassica species of U's triangle model

MAIN CONCLUSION

We comprehensively identified and analyzed the Snf2 gene family. Some Snf2 genes were involved in responding to salt stress based on the RNA-seq and qRT-PCR analysis. Sucrose nonfermenting 2 (Snf2) proteins are core components of chromatin remodeling complexes that not only alter DNA accessibility using the energy of ATP hydrolysis, but also play a critical regulatory role in growth, development, and stress response in eukaryotes. However, the comparative study of Snf2 gene family in the six Brassica species in U's triangle model remains unclear. Here, a total of 405 Snf2 genes were identified, comprising 53, 50, and 46 in the diploid progenitors: Brassica rapa (AA, 2n = 20), Brassica nigra (BB, 2n = 16), and Brassica oleracea (CC, 2n = 18), and 93, 91, and 72 in the allotetraploid: Brassica juncea (AABB, 2n = 36), Brassica napus (AACC, 2n = 38), and Brassica carinata (BBCC, 2n = 34), respectively. These genes were classified into six clades and further divided into 18 subfamilies based on their conserved motifs and domains. Intriguingly, these genes showed highly conserved chromosomal distributions and gene structures, indicating that few dynamic changes occurred during the polyploidization. The duplication modes of the six Brassica species were diverse, and the expansion of most Snf2 in Brassica occurred primarily through dispersed duplication (DSD) events. Additionally, the majority of Snf2 genes were under purifying selection during polyploidization, and some Snf2 genes were associated with various abiotic stresses. Both RNA-seq and qRT-PCR analysis showed that the expression of BnaSnf2 genes was significantly induced under salt stress, implying their involvement in salt tolerance response in Brassica species. The results provide a comprehensive understanding of the Snf2 genes in U's triangle model species, which will facilitate further functional analysis of the Snf2 genes in Brassica plants.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

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 .

Chromosome-level genome assembly of Pedicularis kansuensis illuminates genome evolution of facultative parasitic plant

Parasitic plants have a heterotrophic lifestyle, in which they withdraw all or part of their nutrients from their host through the haustorium. Despite the release of many draft genomes of parasitic plants, the genome evolution related to the parasitism feature of facultative parasites remains largely unknown. In this study, we present a high-quality chromosomal-level genome assembly for the facultative parasite Pedicularis kansuensis (Orobanchaceae), which invades both legume and grass host species in degraded grasslands on the Qinghai-Tibet Plateau. This species has the largest genome size compared with other parasitic species, and expansions of long terminal repeat retrotransposons accounting for 62.37% of the assembly greatly contributed to the genome size expansion of this species. A total of 42,782 genes were annotated, and the patterns of gene loss in P. kansuensis differed from other parasitic species. We also found many mobile mRNAs between P. kansuensis and one of its host species, but these mobile mRNAs could not compensate for the functional losses of missing genes in P. kansuensis. In addition, we identified nine horizontal gene transfer (HGT) events from rosids and monocots, as well as one single-gene duplication events from HGT genes, which differ distinctly from that of other parasitic species. Furthermore, we found evidence for HGT through transferring genomic fragments from phylogenetically remote host species. Taken together, these findings provide genomic insights into the evolution of facultative parasites and broaden our understanding of the diversified genome evolution in parasitic plants and the molecular mechanisms of plant parasitism.

Genome-scale identification, expression and evolution analysis of B-box members in Dendrobium huoshanense

B-box (BBX) proteins have been recognized as vital determinants in plant development, morphogenesis, and adaptive responses to a myriad of environmental stresses. These zinc-finger proteins play a pivotal role in various biological processes. Their influence spans photomorphogenesis, the regulation of flowering, and imparting resilience to a wide array of challenges, encompassing both biotic and abiotic factors. Chromosome localization, gene structure and conserved motifs, phylogenetic analysis, collinearity analysis, expression profiling, fluorescence quantitative analysis, and tobacco transient transformation methods were used for functional localization and expression pattern analysis of the DhBBX gene. A total of 23 DhBBX members were identified from Dendrobium huoshanense. Subsequent phylogenetic evaluations effectively segregated these genes into five discrete evolutionary subsets. The predictions of subcellular localizations revealed that all these proteins were localized in the nucleus. The genetic composition and patterns showed that the majority of these genes consisted of several exons, with a few variations that could be attributed to transposon insertion. A comprehensive analysis using qRT-PCR was conducted to unravel the expression patterns of these genes in D. huoshanense, with a specific concentration on their responses to various hormone treatments and cold stress. Subcellular localization reveals that DhBBX21 and DhBBX9 are located in the nucleus. Our results provide a deep comprehension of the complex regulatory mechanisms of BBXs in response to various environmental and hormonal stimuli. These discoveries encourage further detailed and focused investigations into the operational dynamics of the BBX gene family in a wider range of plant species.

Malate dehydrogenase in plants: evolution, structure, and a myriad of functions

Malate dehydrogenase (MDH) catalyzes the interconversion of oxaloacetate and malate coupled to the oxidation/reduction of coenzymes NAD(P)H/NAD(P)+. While most animals have two isoforms of MDH located in the cytosol and mitochondria, all major groups of land plants have at least six MDHs localized to the cytosol, mitochondria, plastids, and peroxisomes. This family of enzymes participates in important reactions in plant cells including photosynthesis, photorespiration, lipid metabolism, and NH4+ metabolism. MDH also helps to regulate the energy balance in the cell and may help the plant cope with various environmental stresses. Despite their functional diversity, all of the plant MDH enzymes share a similar structural fold and act as dimers. In this review, we will introduce readers to our current understanding of the plant MDHs, including their evolution, structure, and function. The focus will be on the MDH enzymes of the model plant Arabidopsis thaliana.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

Phylogenetic insight into ABCE gene subfamily in plants

ATP-BINDING CASSETTE SUBFAMILY E MEMBER (ABCE) proteins are one of the most conserved proteins across eukaryotes and archaea. Yeast and most animals possess a single ABCE gene encoding the critical translational factor ABCE1. In several plant species, including Arabidopsis thaliana and Oryza sativa, two or more ABCE gene copies have been identified, however information related to plant ABCE gene family is still missing. In this study we retrieved ABCE gene sequences of 76 plant species from public genome databases and comprehensively analyzed them with the reference to A. thaliana ABCE2 gene (AtABCE2). Using bioinformatic approach we assessed the conservation and phylogeny of plant ABCEs. In addition, we performed haplotype analysis of AtABCE2 and its paralogue AtABCE1 using genomic sequences of 1,135 A. thaliana ecotypes. Plant ABCE proteins showed overall high sequence conservation, sharing at least 78% of amino acid sequence identity with AtABCE2. We found that over half of the selected species have two to eight ABCE genes, suggesting that in plants ABCE genes can be classified as a low-copy gene family, rather than a single-copy gene family. The phylogenetic trees of ABCE protein sequences and the corresponding coding sequences demonstrated that Brassicaceae and Poaceae families have independently undergone lineage-specific split of the ancestral ABCE gene. Other plant species have gained ABCE gene copies through more recent duplication events. We also noticed that ploidy level but not ancient whole genome duplications experienced by a species impacts ABCE gene family size. Deeper analysis of AtABCE2 and AtABCE1 from 1,135 A. thaliana ecotypes revealed four and 35 non-synonymous SNPs, respectively. The lower natural variation in AtABCE2 compared to AtABCE1 is in consistence with its crucial role for plant viability. Overall, while the sequence of the ABCE protein family is highly conserved in the plant kingdom, many plants have evolved to have more than one copy of this essential translational factor.

Genome-wide identification and expression analysis of the Dof gene family reveals their involvement in hormone response and abiotic stresses in sunflower (Helianthus annuus L.)

DNA binding with one finger (Dof), plant-specific zinc finger transcription factors, can participate in various physiological and biochemical processes during the life of plants. As one of the most important oil crops in the world, sunflower (Helianthus annuus L.) has significant economic and ornamental value. However, a systematic analysis of H. annuus Dof (HaDof) members and their functions has not been extensively conducted. In this study, we identified 50 HaDof genes that are unevenly distributed on 17 chromosomes of sunflower. We present a comprehensive overview of the HaDof genes, including their chromosome locations, phylogenetic analysis, and expression profile characterization. Phylogenetic analysis classified the 366 Dof members identified from 11 species into four groups (further subdivided into nine subfamilies). Segmental duplications are predominantly contributed to the expansion of sunflower Dof genes, and all segmental duplicate gene pairs are under purifying selection due to strong evolutionary constraints. Furthermore, we observed differential expression patterns for HaDof genes in normal tissues as well as under hormone treatment or abiotic stress conditions by analyzing RNA-seq data from previous studies and RT-qPCR data in our current study. The expression of HaDof04 and HaDof43 were not detected in any samples, which implied that they may be gradually undergoing pseudogenization process. Some HaDof genes, such as HaDof25 and HaDof30, showed responsiveness to exogenous plant hormones, such as kinetin, brassinosteroid, auxin or strigolactone, while others like HaDof15 and HaDof35 may participate in abiotic stress resistance of sunflower seedling. Our study represents the initial step towards understanding the phylogeny and expression characterization of sunflower Dof family genes, which may provide valuable reference information for functional studies on hormone response, abiotic stress resistance, and molecular breeding in sunflower and other species.

Evolutionary conservation of receptor compensation for stem cell homeostasis in Solanaceae plants

Stem cell homeostasis is pivotal for continuous and programmed formation of organs in plants. The precise control of meristem proliferation is mediated by the evolutionarily conserved signaling that encompasses complex interactions among multiple peptide ligands and their receptor-like kinases. Here, we identified compensation mechanisms involving the CLAVATA1 (CLV1) receptor and its paralogs, BARELY ANY MERISTEMs (BAMs), for stem cell proliferation in two Solanaceae species, tomato and groundcherry. Genetic analyses of higher-order mutants deficient in multiple receptor genes, generated via CRISPR-Cas9 genome editing, reveal that tomato SlBAM1 and SlBAM2 compensate for slclv1 mutations. Unlike the compensatory responses between orthologous receptors observed in Arabidopsis, tomato slclv1 mutations do not trigger transcriptional upregulation of four SlBAM genes. The compensation mechanisms within receptors are also conserved in groundcherry, and critical amino acid residues of the receptors associated with the physical interaction with peptide ligands are highly conserved in Solanaceae plants. Our findings demonstrate that the evolutionary conservation of both compensation mechanisms and critical coding sequences between receptor-like kinases provides a strong buffering capacity during stem cell homeostasis in tomato and groundcherry.

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77 266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

The auxin response factor (ARF) gene family in Cyclocarya paliurus: genome-wide identification and their expression profiling under heat and drought stresses

Auxin response factors (ARFs), as the main components of auxin signaling, play a crucial role in various processes of plant growth and development, as well as in stress response. So far, there have been no reports on the genome-wide identification of the ARF transcription factor family in Cyclocarya paliurus, a deciduous tree plant in the family Juglaceae. In this study, a total of 34 CpARF genes were identified based on whole genome sequence, and they were unevenly distributed on 16 chromosomes, with the highest distribution on chromosome 6. Domain analysis of CpARF proteins displayed that 31 out of 34 CpARF proteins contain a typical B3 domain (DBD domain), except CpARF12/ CpARF14/CpARF31, which all belong to Class VI. And 20 CpARFs (58.8%) contain an auxin_IAA binding domain, and are mainly distributed in classes I, and VI. Phylogenetic analysis showed that CpARF was divided into six classes (I-VI), each containing 4, 4, 1, 8, 4, and 13 members, respectively. Gene duplication analysis showed that there are 14 segmental duplications and zero tandem repeats were identified in the CpARF gene family of the C. paliurus genome. The Ka/Ks ratio of duplicate gene pairs indicates that CpARF genes are subjected to strong purification selection pressure. Synteny analysis showed that C. paliurus shared the highest homology in 74 ARF gene pairs with Juglans regia, followed by 73, 51, 25, and 11 homologous gene pairs with Populus trichocarpa, Juglans cathayensis, Arabidopsis, and rice, respectively. Promoter analysis revealed that 34 CpARF genes had cis-elements related to hormones, stress, light, and growth and development except for CpARF12. The expression profile analysis showed that almost all CpARF genes were differentially expressed in at least one tissue, and several CpARF genes displayed tissue-specific expression. Furthermore, 24 out of the 34 CpARF genes have significantly response to drought stress (P < 0.05), and most of them (16) being significantly down-regulated under moderate drought treatment. Meanwhile, the majority of CpARF genes (28) have significantly response to drought stress (P < 0.05), and most of them (26) are significantly down-regulated under severe drought treatment. Furthermore, 32 out of the 34 CpARF genes have significantly response to high, middle, and low salt stress under salt treatment (P < 0.05). Additionally, subcellular localization analysis confirmed that CpARF16 and CpARF32 were all localized to nucleus. Thus, our findings expand the understanding of the function of CpARF genes and provide a basis for further functional studies on CpARF genes in C. paliurus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01474-1.

Upgraded durian genome reveals the role of chromosome reshuffling during ancestral karyotype evolution, lignin biosynthesis regulation, and stress tolerance

Durian (Durio zibethinus) is a tropical fruit that has a unique flavor and aroma. It occupies a significant phylogenetic position within the Malvaceae family. Extant core-eudicot plants are reported to share seven ancestral karyotypes that have undergone reshuffling, resulting in an abundant genomic diversity. However, the ancestral karyotypes of the Malvaceae family, as well as the evolution trajectory leading to the 28 chromosomes in durian, remain poorly understood. Here, we report the high-quality assembly of the durian genome with comprehensive comparative genomic analyses. By analyzing the collinear blocks between cacao and durian, we inferred 11 Malvaceae ancestral karyotypes. These blocks were present in a single-copy form in cacao and mainly in triplicates in durian, possibly resulting from a recent whole genome triplication (WGT) event that led to hexaploidization of the durian genome around 20 (17-24) million years ago. A large proportion of the duplicated genes in durian, such as those involved in the lignin biosynthesis module for phenylpropane biosynthesis, are derived directly from whole genome duplication, which makes it an important force in reshaping its genomic architecture. Transcriptome studies have revealed that genes involved in feruloyl-CoA formations were highly preferentially expressed in fruit peels, indicating that the thorns produced on durian fruit may comprise guaiacyl and syringyl lignins. Among all the analyzed transcription factors (TFs), members of the heat shock factor family (HSF) were the most significantly upregulated under heat stress. All subfamilies of genes encoding heat shock proteins (HSPs) in the durian genome appear to have undergone expansion. The potential interactions between HSF Dzi05.397 and HSPs were examined and experimentally verified. Our study provides a high-quality durian genome and reveals the reshuffling mechanism of ancestral Malvaceae chromosomes to produce the durian genome. We also provide insights into the mechanism underlying lignin biosynthesis and heat stress tolerance.

Chromosome-level genome assembly and characterization of the Calophaca sinica genome

Calophaca sinica is a rare plant endemic to northern China which belongs to the Fabaceae family and possesses rich nutritional value. To support the preservation of the genetic resources of this plant, we have successfully generated a high-quality genome of C. sinica (1.06 Gb). Notably, transposable elements (TEs) constituted ~73% of the genome, with long terminal repeat retrotransposons (LTR-RTs) dominating this group of elements (~54% of the genome). The average intron length of the C. sinica genome was noticeably longer than what has been observed for closely related species. The expansion of LTR-RTs and elongated introns emerged had the largest influence on the enlarged genome size of C. sinica in comparison to other Fabaceae species. The proliferation of TEs could be explained by certain modes of gene duplication, namely, whole genome duplication (WGD) and dispersed duplication (DSD). Gene family expansion, which was found to enhance genes associated with metabolism, genetic maintenance, and environmental stress resistance, was a result of transposed duplicated genes (TRD) and WGD. The presented genomic analysis sheds light on the genetic architecture of C. sinica, as well as provides a starting point for future evolutionary biology, ecology, and functional genomics studies centred around C. sinica and closely related species.

Evolutionary origin and gradual accumulation with plant evolution of the LACS family

BACKGROUND

LACS (long-chain acyl-CoA synthetase) genes are widespread in organisms and have multiple functions in plants, especially in lipid metabolism. However, the origin and evolutionary dynamics of the LACS gene family remain largely unknown.

RESULTS

Here, we identified 1785 LACS genes in the genomes of 166 diverse plant species and identified the clades (I, II, III, IV, V, VI) of six clades for the LACS gene family of green plants through phylogenetic analysis. Based on the evolutionary history of plant lineages, we found differences in the origins of different clades, with Clade IV originating from chlorophytes and representing the origin of LACS genes in green plants. The structural characteristics of different clades indicate that clade IV is relatively independent, while the relationships between clades (I, II, III) and clades (V, VI) are closer. Dispersed duplication (DSD) and transposed duplication (TRD) are the main forces driving the evolution of plant LACS genes. Network clustering analysis further grouped all LACS genes into six main clusters, with genes within each cluster showing significant co-linearity. Ka/Ks results suggest that LACS family genes underwent purifying selection during evolution. We analyzed the phylogenetic relationships and characteristics of six clades of the LACS gene family to explain the origin, evolutionary history, and phylogenetic relationships of different clades and proposed a hypothetical evolutionary model for the LACS family of genes in plants.

CONCLUSIONS

Our research provides genome-wide insights into the evolutionary history of the LACS gene family in green plants. These insights lay an important foundation for comprehensive functional characterization in future research.

Complete telomere-to-telomere genomes uncover virulence evolution conferred by chromosome fusion in oomycete plant pathogens

Variations in chromosome number are occasionally observed among oomycetes, a group that includes many plant pathogens, but the emergence of such variations and their effects on genome and virulence evolution remain ambiguous. We generated complete telomere-to-telomere genome assemblies for Phytophthora sojae, Globisporangium ultimum, Pythium oligandrum, and G. spinosum. Reconstructing the karyotype of the most recent common ancestor in Peronosporales revealed that frequent chromosome fusion and fission drove changes in chromosome number. Centromeres enriched with Copia-like transposons may contribute to chromosome fusion and fission events. Chromosome fusion facilitated the emergence of pathogenicity genes and their adaptive evolution. Effectors tended to duplicate in the sub-telomere regions of fused chromosomes, which exhibited evolutionary features distinct to the non-fused chromosomes. By integrating ancestral genomic dynamics and structural predictions, we have identified secreted Ankyrin repeat-containing proteins (ANKs) as a novel class of effectors in P. sojae. Phylogenetic analysis and experiments further revealed that ANK is a specifically expanded effector family in oomycetes. These results revealed chromosome dynamics in oomycete plant pathogens, and provided novel insights into karyotype and effector evolution.

Jasmonate ZIM Domain Protein (JAZ) Gene SLJAZ15 Increases Resistance to Orobanche aegyptiaca in Tomato

Orobanche aegyptiaca Pers. is a holoparasitic plant that severely reduces tomato (Solanum lycopersicum L.) production in China. However, there is a lack of effective control methods and few known sources of genetic resistance. In this study, we focused on key genes in the JAZ family, comparing the JAZ family in Arabidopsis thaliana (L. Heynh.) to the tomato genome. After identifying the JAZ family members in S. lycopersicum, we performed chromosomal localization and linear analysis with phylogenetic relationship analysis of the JAZ family. We also analyzed the gene structure of the JAZ gene family members in tomato and the homology of the JAZ genes among the different species to study their relatedness. The key genes for O. aegyptiaca resistance were identified using VIGS (virus-induced gene silencing), and the parasitization rate of silenced tomato plants against O. aegyptiaca increased by 47.23-91.13%. The genes were localized in the nucleus by subcellular localization. Heterologous overexpression in A. thaliana showed that the key gene had a strong effect on the parasitization process of O. aegyptiaca, and the overexpression of the key gene reduced the parasitization rate of O. aegyptiaca 1.69-fold. Finally, it was found that the SLJAZ15 gene can positively regulate the hormone content in tomato plants and affect plant growth and development, further elucidating the function of this gene.

Genome-wide identification and expression analysis of protein arginine methyltransferase and JmjC domain-containing family in apple

Histone methylation is an important type of histone modification that regulates gene expression in plants. In this study, we identified 14 arginine methylation-related genes (Protein Arginine Methyltransferase, MdPRMT) and 32 demethylation-related genes (JmjC Domain-Containing Family, MdJMJ) in apple. Furthermore, we investigated the phylogenetic relationship, chromosome distribution, gene structure, motif analysis, promoter sequence analysis, and expression patterns of MdPRMT and MdJMJ genes. Homology analysis showed a high degree of conservation and homology between PRMT and JMJ genes in Arabidopsis and apple. We identified the types of duplicated genes in the MdJMJ and MdPRMT gene families, found a large number of whole-genome duplicates (WGD) gene pairs and a small number of tandem duplicates (TD) pairs, transposed duplication (TRD) gene pairs as well as proximal duplicates (PD) pairs, and discussed the possible evolutionary pathways of the gene families from the perspective of duplicated genes. Homology analysis showed a high degree of conservation and homology between PRMT and JMJ genes in Arabidopsis and apple. In addition, the promoter regions of MdPRMT and MdJMJ contain numerous cis-acting elements involved in plant growth and development, hormone response, and stress responses. Based on the transcriptional profiles of MdPRMT and MdJMJ in different tissues and developmental stages, it was found that MdPRMT and MdJMJ may play multiple roles in apple growth and development, for example, MdJMJ21 may be involved in the regulation of apple endosperm formation. MdPRMT and MdJMJ exhibit different expression patterns in response to hormone signaling in apple, MdJMJ3, MdJMJ18, MdJMJ30, MdPRMT2, MdPRMT13, and MdPRMT14 may play roles in apple response to drought stress, while the expression of MdJMJ13, MdPRMT3, MdPRMT4, and MdPRMT6 is affected by cold stress. Our study provides a foundation for determining the functional roles of MdPRMT and MdJMJ genes in apple.

Genome-wide analysis of NAC transcription factors and exploration of candidate genes regulating selenium metabolism in Broussonetia papyrifera

MAIN CONCLUSION

Genome-wide identification revealed 79 BpNAC genes belonging to 16 subfamilies, and their gene structures and evolutionary relationships were characterized. Expression analysis highlighted their importance in plant selenium stress responses. Paper mulberry (Broussonetia papyrifera), a deciduous arboreal plant of the Moraceae family, is distinguished by its leaves, which are abundant in proteins, polysaccharides, and flavonoids, positioning it as a novel feedstock. NAC transcription factors, exclusive to plant species, are crucial in regulating growth, development, and response to biotic and abiotic stress. However, extensive characterization of the NAC family within paper mulberry is lacking. In this study, 79 BpNAC genes were identified from the paper mulberry genome, with an uneven distribution across 13 chromosomes. A comprehensive, genome-wide analysis of BpNACs was performed, including investigating gene structures, promoter regions, and chromosomal locations. Phylogenetic tree analysis, alongside comparisons with Arabidopsis thaliana NACs, allowed for categorizing these genes into 16 subfamilies in alignment with gene structure and motif conservation. Collinearity analysis suggested a significant homologous relationship between the NAC genes of paper mulberry and those in Morus notabilis, Ficus hispida, Antiaris toxicaria, and Cannabis sativa. Integrating transcriptome data and Se content revealed that 12 BpNAC genes were associated with selenium biosynthesis. Subsequent RT-qPCR analysis corroborated the correlation between BpNAC59, BpNAC62 with sodium selenate, and BpNAC55 with sodium selenite. Subcellular localization experiments revealed the nuclear functions of BpNAC59 and BpNAC62. This study highlights the potential BpNAC transcription factors involved in selenium metabolism, providing a foundation for strategically breeding selenium-fortified paper mulberry.

An integrated multi-omics approach reveals polymethoxylated flavonoid biosynthesis in Citrus reticulata cv. Chachiensis

Citrus reticulata cv. Chachiensis (CRC) is an important medicinal plant, its dried mature peels named "Guangchenpi", has been used as a traditional Chinese medicine to treat cough, indigestion, and lung diseases for several hundred years. However, the biosynthesis of the crucial natural products polymethoxylated flavonoids (PMFs) in CRC remains unclear. Here, we report a chromosome-scale genome assembly of CRC with the size of 314.96 Mb and a contig N50 of 16.22 Mb. Using multi-omics resources, we discover a putative caffeic acid O-methyltransferase (CcOMT1) that can transfer a methyl group to the 3-hydroxyl of natsudaidain to form 3,5,6,7,8,3',4'-heptamethoxyflavone (HPMF). Based on transient overexpression and virus-induced gene silencing experiments, we propose that CcOMT1 is a candidate enzyme in HPMF biosynthesis. In addition, a potential gene regulatory network associated with PMF biosynthesis is identified. This study provides insights into PMF biosynthesis and may assist future research on mining genes for the biosynthesis of plant-based medicines.

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.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

Genome-Wide Identification and Evolutionary Analysis of Receptor-like Kinase Family Genes Provides Insights into Anthracnose Resistance of Dioscorea alata

Dioscorea alata, commonly known as "greater yam", is a vital crop in tropical and subtropical regions of the world, yet it faces significant threats from anthracnose disease, mainly caused by Colletotrichum gloeosporioides. However, exploring disease resistance genes in this species has been challenging due to the difficulty of genetic mapping resulting from the loss of the flowering trait in many varieties. The receptor-like kinase (RLK) gene family represents essential immune receptors in plants. In this study, genomic analysis revealed 467 RLK genes in D. alata. The identified RLKs were distributed unevenly across chromosomes, likely due to tandem duplication events. However, a considerable number of ancient whole-genome or segmental duplications dating back over 100 million years contributed to the diversity of RLK genes. Phylogenetic analysis unveiled at least 356 ancient RLK lineages in the common ancestor of Dioscoreaceae, which differentially inherited and expanded to form the current RLK profiles of D. alata and its relatives. The analysis of cis-regulatory elements indicated the involvement of RLK genes in diverse stress responses. Transcriptome analysis identified RLKs that were up-regulated in response to C. gloeosporioides infection, suggesting their potential role in resisting anthracnose disease. These findings provide novel insights into the evolution of RLK genes in D. alata and their potential contribution to disease resistance.

Genome-wide characterization of the NBLRR gene family provides evolutionary and functional insights into blast resistance in pearl millet (Cenchrus americanus (L.) Morrone)

MAIN CONCLUSION

The investigation is the first report on genome-wide identification and characterization of NBLRR genes in pearl millet. We have shown the role of gene loss and purifying selection in the divergence of NBLRRs in Poaceae lineage and candidate CaNBLRR genes for resistance to Magnaporthe grisea infection. Plants have evolved multiple integral mechanisms to counteract the pathogens' infection, among which plant immunity through NBLRR (nucleotide-binding site, leucine-rich repeat) genes is at the forefront. The genome-wide mining in pearl millet (Cenchrus americanus (L.) Morrone) revealed 146 CaNBLRRs. The variation in the branch length of NBLRRs showed the dynamic nature of NBLRRs in response to evolving pathogen races. The orthology of NBLRRs showed a predominance of many-to-one orthologs, indicating the divergence of NBLRRs in the pearl millet lineage mainly through gene loss events followed by gene gain through single-copy duplications. Further, the purifying selection (Ka/Ks < 1) shaped the expansion of NBLRRs within the lineage of pear millet and other members of Poaceae. Presence of cis-acting elements, viz. TCA element, G-box, MYB, SARE, ABRE and conserved motifs annotated with P-loop, kinase 2, RNBS-A, RNBS-D, GLPL, MHD, Rx-CC and LRR suggests their putative role in disease resistance and stress regulation. The qRT-PCR analysis in pearl millet lines showing contrasting responses to Magnaporthe grisea infection identified CaNBLRR20, CaNBLRR33, CaNBLRR46 CaNBLRR51, CaNBLRR78 and CaNBLRR146 as putative candidates. Molecular docking showed the involvement of three and two amino acid residues of LRR domains forming hydrogen bonds (histidine, arginine and threonine) and salt bridges (arginine and lysine) with effectors. Whereas 14 and 20 amino acid residues of CaNBLRR78 and CaNBLRR20 showed hydrophobic interactions with 11 and 9 amino acid residues of effectors, Mg.00g064570.m01 and Mg.00g006570.m01, respectively. The present investigation gives a comprehensive overview of CaNBLRRs and paves the foundation for their utility in pearl millet resistance breeding through understanding of host-pathogen interactions.

Functional Divergence in Orthologous Transcription Factors: Insights from AtCBF2/3/1 and OsDREB1C

Despite traditional beliefs of orthologous genes maintaining similar functions across species, growing evidence points to their potential for functional divergence. C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) are critical in cold acclimation, with their overexpression enhancing stress tolerance but often constraining plant growth. In contrast, a recent study unveiled a distinctive role of rice OsDREB1C in elevating nitrogen use efficiency (NUE), photosynthesis, and grain yield, implying functional divergence within the CBF/DREB1 orthologs across species. Here, we delve into divergent molecular mechanisms of OsDREB1C and AtCBF2/3/1 by exploring their evolutionary trajectories across rice and Arabidopsis genomes, regulatomes, and transcriptomes. Evolutionary scrutiny shows discrete clades for OsDREB1C and AtCBF2/3/1, with the Poaceae-specific DREB1C clade mediated by a transposon event. Genome-wide binding profiles highlight OsDREB1C's preference for GCCGAC compared to AtCBF2/3/1's preference for A/GCCGAC, a distinction determined by R12 in the OsDREB1C AP2/ERF domain. Cross-species multiomic analyses reveal shared gene orthogroups (OGs) and underscore numerous specific OGs uniquely bound and regulated by OsDREB1C, implicated in NUE, photosynthesis, and early flowering, or by AtCBF2/3/1, engaged in hormone and stress responses. This divergence arises from gene gains/losses (∼16.7% to 25.6%) and expression reprogramming (∼62.3% to 66.2%) of OsDREB1C- and AtCBF2/3/1-regulated OGs during the extensive evolution following the rice-Arabidopsis split. Our findings illustrate the regulatory evolution of OsDREB1C and AtCBF2/3/1 at a genomic scale, providing insights on the functional divergence of orthologous transcription factors following gene duplications across species.

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.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

MYB transcription factors encoded by diversified tandem gene clusters cause varied Morella rubra fruit color

Chinese bayberry (Morella rubra) is a fruit tree with a remarkable variation in fruit color, ranging from white to dark red as determined by anthocyanin content. In dark red "Biqi" (BQ), red "Dongkui" (DK), pink "Fenhong" (FH), and white "Shuijing" (SJ), we identified an anthocyanin-related MYB transcription factor-encoding gene cluster of four members, i.e. MrMYB1.1, MrMYB1.2, MrMYB1.3, and MrMYB2. Collinear analysis revealed that the MYB tandem cluster may have occurred in a highly conserved region of many eudicot genomes. Two alleles of MrMYB1.1 were observed; MrMYB1.1-1 (MrMYB1.1n) was a full-length allele and homozygous in "BQ", MrMYB1.1-2 (MrMYB1.1d) was a nonfunctional allele with a single base deletion and homozygous in "SJ", and MrMYB1.1n/MrMYB1.1d were heterozygous in "DK" and "FH". In these four cultivars, expression of MrMYB1.1, MrMYB1.2, and MrMYB2 was enhanced during ripening. Both alleles were equally expressed in MrMYB1.1n/MrMYB1.1d heterozygous cultivars as revealed by a cleaved amplified polymorphic sequence marker. Expression of MrMYB1.3 was restricted to some dark red cultivars only. Functional characterization revealed that MrMYB1.1n and MrMYB1.3 can induce anthocyanin accumulation while MrMYB1.1d, MrMYB1.2, and MrMYB2 cannot. DNA-protein interaction assays indicated that MrMYB1.1n and MrMYB1.3 can directly bind to and activate the promoters of anthocyanin-related genes via interaction with a MYC-like basic helix-loop-helix protein MrbHLH1. We concluded that the specific genotype of MrMYB1.1 alleles, as well as the exclusive expression of MrMYB1.3 in some dark red cultivars, contributes to fruit color variation. The study provides insights into the mechanisms for regulation of plant anthocyanin accumulation by MYB tandem clusters.

Genome-Wide Analysis of MYB Gene Family in Chrysanthemum ×morifolium Provides Insights into Flower Color Regulation

MYBs constitute the second largest transcription factor (TF) superfamily in flowering plants with substantial structural and functional diversity, which have been brought into focus because they affect flower colors by regulating anthocyanin biosynthesis. Up to now, the genomic data of several Chrysanthemum species have been released, which provides us with abundant genomic resources for revealing the evolution of the MYB gene family in Chrysanthemum species. In the present study, comparative analyses of the MYB gene family in six representative species, including C. lavandulifolium, C. seticuspe, C. ×morifolium, Helianthus annuus, Lactuca sativa, and Arabidopsis thaliana, were performed. A total of 1104 MYBs, which were classified into four subfamilies and 35 lineages, were identified in the three Chrysanthemum species (C. lavandulifolium, C. seticuspe, and C. ×morifolium). We found that whole-genome duplication and tandem duplication are the main duplication mechanisms that drove the occurrence of duplicates in CmMYBs (particularly in the R2R3-MYB subfamily) during the evolution of the cultivated chrysanthemums. Sequence structure and selective pressure analyses of the MYB gene family revealed that some of R2R3-MYBs were subjected to positive selection, which are mostly located on the distal telomere segments of the chromosomes and contain motifs 7 and 8. In addition, the gene expression analysis of CmMYBs in different organs and at various capitulum developmental stages of C. ×morifolium indicated that CmMYBS2, CmMYB96, and CmMYB109 might be the negative regulators for anthocyanin biosynthesis. Our results provide the phylogenetic context for research on the genetic and functional evolution of the MYB gene family in Chrysanthemum species and deepen our understanding of the regulatory mechanism of MYB TFs on the flower color of C. ×morifolium.

Genome-wide identification and analysis of ascorbate peroxidase (APX) gene family in hemp (Cannabis sativa L.) under various abiotic stresses

Ascorbate peroxidase (APX) plays a critical role in molecular mechanisms such as plant development and defense against abiotic stresses. As an important economic crop, hemp (Cannabis sativa L.) is vulnerable to adverse environmental conditions, such as drought, cold, salt, and oxidative stress, which lead to a decline in yield and quality. Although APX genes have been characterized in a variety of plants, members of the APX gene family in hemp have not been completely identified. In this study, we (1) identified eight members of the CsAPX gene family in hemp and mapped their locations on the chromosomes using bioinformatics analysis; (2) examined the physicochemical characteristics of the proteins encoded by these CsAPX gene family members; (3) investigated their intraspecific collinearity, gene structure, conserved domains, conserved motifs, and cis-acting elements; (4) constructed a phylogenetic tree and analyzed interspecific collinearity; and (5) ascertained expression differences in leaf tissue subjected to cold, drought, salt, and oxidative stresses using quantitative real-time-PCR (qRT-PCR). Under all four stresses, CsAPX6, CsAPX7, and CsAPX8 consistently exhibited significant upregulation, whereas CsAPX2 displayed notably higher expression levels under drought stress than under the other stresses. Taken together, the results of this study provide basic genomic information on the expression of the APX gene family and pave the way for studying the role of APX genes in abiotic stress.

Genome composition and GC content influence loci distribution in reduced representation genomic studies

BACKGROUND

Genomic architecture is a key evolutionary trait for living organisms. Due to multiple complex adaptive and neutral forces which impose evolutionary pressures on genomes, there is a huge variability of genomic features. However, their variability and the extent to which genomic content determines the distribution of recovered loci in reduced representation sequencing studies is largely unexplored.

RESULTS

Here, by using 80 genome assemblies, we observed that whereas plants primarily increase their genome size by expanding their intergenic regions, animals expand both intergenic and intronic regions, although the expansion patterns differ between deuterostomes and protostomes. Loci mapping in introns, exons, and intergenic categories obtained by in silico digestion using 2b-enzymes are positively correlated with the percentage of these regions in the corresponding genomes, suggesting that loci distribution mostly mirrors genomic architecture of the selected taxon. However, exonic regions showed a significant enrichment of loci in all groups regardless of the used enzyme. Moreover, when using selective adaptors to obtain a secondarily reduced loci dataset, the percentage and distribution of retained loci also varied. Adaptors with G/C terminals recovered a lower percentage of selected loci, with a further enrichment of exonic regions, while adaptors with A/T terminals retained a higher percentage of loci and slightly selected more intronic regions than expected.

CONCLUSIONS

Our results highlight how genome composition, genome GC content, RAD enzyme choice and use of base-selective adaptors influence reduced genome representation techniques. This is important to acknowledge in population and conservation genomic studies, as it determines the abundance and distribution of loci.

TIP aquaporins in Cyperus esculentus: genome-wide identification, expression profiles, subcellular localizations, and interaction patterns

BACKGROUND

Tonoplast intrinsic proteins (TIPs), which typically mediate water transport across vacuolar membranes, play an essential role in plant growth, development, and stress responses. However, their characterization in tigernut (Cyperus esculentus L.), an oil-bearing tuber plant of the Cyperaceae family, is still in the infancy.

RESULTS

In this study, a first genome-wide characterization of the TIP subfamily was conducted in tigernut, resulting in ten members representing five previously defined phylogenetic groups, i.e., TIP1-5. Although the gene amounts are equal to that present in two model plants Arabidopsis and rice, the group composition and/or evolution pattern were shown to be different. Except for CeTIP1;3 that has no counterpart in both Arabidopsis and rice, complex orthologous relationships of 1:1, 1:2, 1:3, 2:1, and 2:2 were observed. Expansion of the CeTIP subfamily was contributed by whole-genome duplication (WGD), transposed, and dispersed duplications. In contrast to the recent WGD-derivation of CeTIP3;1/-3;2, synteny analyses indicated that TIP4 and - 5 are old WGD repeats of TIP2, appearing sometime before monocot-eudicot divergence. Expression analysis revealed that CeTIP genes exhibit diverse expression profiles and are subjected to developmental and diurnal fluctuation regulation. Moreover, when transiently overexpressed in tobacco leaves, CeTIP1;1 was shown to locate in the vacuolar membrane and function in homo/heteromultimer, whereas CeTIP2;1 is located in the cell membrane and only function in heteromultimer. Interestingly, CeTIP1;1 could mediate the tonoplast-localization of CeTIP2;1 via protein interaction, implying complex regulatory patterns.

CONCLUSIONS

Our findings provide a global view of CeTIP genes, which provide valuable information for further functional analysis and genetic improvement through manipulating key members in tigernut.

Genome-wide identification and stress response analysis of BcaCPK gene family in amphidiploid Brassica carinata

BACKGROUND

Calcium-dependent protein kinases (CPKs) are crucial for recognizing and transmitting Ca2+ signals in plant cells, playing a vital role in growth, development, and stress response. This study aimed to identify and detect the potential roles of the CPK gene family in the amphidiploid Brassica carinata (BBCC, 2n = 34) using bioinformatics methods.

RESULTS

Based on the published genomic information of B. carinata, a total of 123 CPK genes were identified, comprising 70 CPK genes on the B subgenome and 53 on the C subgenome. To further investigate the homologous evolutionary relationship between B. carinata and other plants, the phylogenetic tree was constructed using CPKs in B. carinata and Arabidopsis thaliana. The phylogenetic analysis classified 123 family members into four subfamilies, where gene members within the same subfamily exhibited similar conserved motifs. Each BcaCPK member possesses a core protein kinase domain and four EF-hand domains. Most of the BcaCPK genes contain 5 to 8 introns, and these 123 BcaCPK genes are unevenly distributed across 17 chromosomes. Among these BcaCPK genes, 120 replicated gene pairs were found, whereas only 8 genes were tandem duplication, suggesting that dispersed duplication mainly drove the family amplification. The results of the Ka/Ks analysis indicated that the CPK gene family of B. carinata was primarily underwent purification selection in evolutionary selection. The promoter region of most BcaCPK genes contained various stress-related cis-acting elements. qRT-PCR analysis of 12 selected CPK genes conducted under cadmium and salt stress at various points revealed distinct expression patterns among different family members in response to different stresses. Specifically, the expression levels of BcaCPK2.B01a, BcaCPK16.B02b, and BcaCPK26.B02 were down-regulated under both stresses, whereas the expression levels of other members were significantly up-regulated under at least one stress.

CONCLUSION

This study systematically identified the BcaCPK gene family in B. carinata, which contributes to a better understanding the CPK genes in this species. The findings also serve as a reference for analyzing stress responses, particularly in relation to cadmium and salt stress in B. carinata.

GATA transcription factor in common bean: A comprehensive genome-wide functional characterization, identification, and abiotic stress response evaluation

The GATA transcription factors (TFs) have been extensively studied for its regulatory role in various biological processes in many plant species. The functional and molecular mechanism of GATA TFs in regulating tolerance to abiotic stress has not yet been studied in the common bean. This study analyzed the functional identity of the GATA gene family in the P. vulgaris genome under different abiotic and phytohormonal stress. The GATA gene family was systematically investigated in the P. vulgaris genome, and 31 PvGATA TFs were identified. The study found that 18 out of 31 PvGATA genes had undergone duplication events, emphasizing the role of gene duplication in GATA gene expansion. All the PvGATA genes were classified into four significant subfamilies, with 8, 3, 6, and 13 members in each subfamily (subfamilies I, II, III, and IV), respectively. All PvGATA protein sequences contained a single GATA domain, but subfamily II members had additional domains such as CCT and tify. A total of 799 promoter cis-regulatory elements (CREs) were predicted in the PvGATAs. Additionally, we used qRT-PCR to investigate the expression profiles of five PvGATA genes in the common bean roots under abiotic conditions. The results suggest that PvGATA01/10/25/28 may play crucial roles in regulating plant resistance against salt and drought stress and may be involved in phytohormone-mediated stress signaling pathways. PvGATA28 was selected for overexpression and cloned into N. benthamiana using Agrobacterium-mediated transformation. Transgenic lines were subjected to abiotic stress, and results showed a significant tolerance of transgenic lines to stress conditions compared to wild-type counterparts. The seed germination assay suggested an extended dormancy of transgenic lines compared to wild-type lines. This study provides a comprehensive analysis of the PvGATA gene family, which can serve as a foundation for future research on the function of GATA TFs in abiotic stress tolerance in common bean plants.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Allopolyploid origin and diversification of the Hawaiian endemic mints

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

Insights into cucumber (Cucumis sativus) genetics: Genome-wide discovery and computational analysis of the Calreticulin Domain-Encoding gene (CDEG) family

Cucumber is an essential vegetable crop throughout the world. Cucumber development is vital for accomplishing both quality and productivity requirements. Meanwhile, numerous factors have resulted in substantial cucumber losses. However, the calreticulin domain-encoding genes (CDEGs) in cucumber were not well-characterized and had little function. In the genome-wide association study (GWAS), we recognized and characterized the CDEGs in Cucumis sativus (cucumber). Through a comprehensive study of C. sativus, our research has unveiled the presence of three unique genes, denoted as CsCRTb, CsCRT3, and CsCNX1, unevenly distributed on three chromosomes in the genome of C. sativus. In accordance to the phylogenetic investigation, these genes may be categorized into three subfamilies. Based on the resemblance with AtCDE genes, we reorganized the all CsCDE genes in accordance with international nomenclature. The expression analysis and cis-acting components revealed that each of CsCDE gene promoter region enclosed number of cis-elements connected with hormone and stress response. According to subcellular localization studies demonstrated that, they were found in deferent locations of the cell such as endoplasmic reticulum, plasma membrane, golgi apparatus, and vacuole, according to subcellular localization studies. Chromosomal distribution analysis and synteny analysis demonstrated the probability of segmental or tandem duplications within the cucumber CDEG gene family. Additionally, miRNAs displayed diverse modes of action, including mRNA cleavage and translational inhibition. We used the RNA seq data to analyze the expression of CDEG genes in response to cold stress and also improved cold tolerance, which was brought on by treating cucumber plants to an exogenous chitosan oligosaccharide spray. Our investigation revealed that these genes responded to this stress in a variety of ways, demonstrating that they may adapt quickly to environmental changes in cucumber plants. This study provides a base for further understanding in reference to CDE gene family and reveals that genes play significant functions in cucumber stress responses.

Haplotype-resolved genome assembly of Phanera championii reveals molecular mechanisms of flavonoid synthesis and adaptive evolution

Phanera championii is a medicinal liana plant that has successfully adapted to hostile karst habitats. Despite extensive research on its medicinal components and pharmacological effects, the molecular mechanisms underlying the biosynthesis of critical flavonoids and its adaptation to karst habitats remain elusive. In this study, we performed high-coverage PacBio and Hi-C sequencing of P. championii, which revealed its high heterozygosity and phased the genome into two haplotypes: Hap1 (384.60 Mb) and Hap2 (383.70 Mb), encompassing a total of 58 612 annotated genes. Comparative genomes analysis revealed that P. championii experienced two whole-genome duplications (WGDs), with approximately 59.59% of genes originating from WGD events, thereby providing a valuable genetic resource for P. championii. Moreover, we identified a total of 112 genes that were strongly positively selected. Additionally, about 81.60 Mb of structural variations between the two haplotypes. The allele-specific expression patterns suggested that the dominant effect of P. championii was the elimination of deleterious mutations and the promotion of beneficial mutations to enhance fitness. Moreover, our transcriptome and metabolome analysis revealed alleles in different tissues or different haplotypes collectively regulate the synthesis of flavonoid metabolites. In summary, our comprehensive study highlights the significance of genomic and morphological adaptation in the successful adaptation of P. championii to karst habitats. The high-quality phased genomes obtained in this study serve as invaluable genomic resources for various applications, including germplasm conservation, breeding, evolutionary studies, and elucidation of pathways governing key biological traits of P. championii.

Drying without dying: A genome database for desiccation-tolerant plants and evolution of desiccation tolerance

Desiccation is typically fatal, but a small number of land plants have evolved vegetative desiccation tolerance (VDT), allowing them to dry without dying through a process called anhydrobiosis. Advances in sequencing technologies have enabled the investigation of genomes for desiccation-tolerant plants over the past decade. However, a dedicated and integrated database for these valuable genomic resources has been lacking. Our prolonged interest in VDT plant genomes motivated us to create the "Drying without Dying" database, which contains a total of 16 VDT-related plant genomes (including 10 mosses) and incorporates 10 genomes that are closely related to VDT plants. The database features bioinformatic tools, such as blast and homologous cluster search, sequence retrieval, Gene Ontology term and metabolic pathway enrichment statistics, expression profiling, co-expression network extraction, and JBrowser exploration for each genome. To demonstrate its utility, we conducted tailored PFAM family statistical analyses, and we discovered that the drought-responsive ABA transporter AWPM-19 family is significantly tandemly duplicated in all bryophytes but rarely so in tracheophytes. Transcriptomic investigations also revealed that response patterns following desiccation diverged between bryophytes and angiosperms. Combined, the analyses provided genomic and transcriptomic evidence supporting a possible divergence and lineage-specific evolution of VDT in plants. The database can be accessed at http://desiccation.novogene.com. We expect this initial release of the "Drying without Dying" plant genome database will facilitate future discovery of VDT genetic resources.

Adaptive evolution of chloroplast division mechanisms during plant terrestrialization

Despite extensive research, the origin and evolution of the chloroplast division machinery remain unclear. Here, we employ recently sequenced genomes and transcriptomes of Archaeplastida clades to identify the core components of chloroplast division and reconstruct their evolutionary histories, respectively. Our findings show that complete division ring structures emerged in Charophytes. We find that Glaucophytes experienced strong selection pressure, generating diverse variants adapted to the changing terrestrial environments. By integrating the functions of chloroplast division genes (CDGs) annotated in a workflow developed using large-scale multi-omics data, we further show that dispersed duplications acquire more species-specific functions under stronger selection pressures. Notably, PARC6, a dispersed duplicate CDG, regulates leaf color and plant growth in Solanum lycopersicum, demonstrating neofunctionalization. Our findings provide an integrated perspective on the functional evolution of chloroplast division machinery and highlight the potential of dispersed duplicate genes as the primary source of adaptive evolution of chloroplast division.

Identification, Evolutionary Dynamics, and Gene Expression Patterns of the ACP Gene Family in Responding to Salt Stress in Brassica Genus

Acyl carrier proteins (ACPs) have been reported to play a crucial role in responding to biotic and abiotic stresses, regulating growth and development. However, the biological function of the ACP gene family in the Brassica genus has been limited until now. In this study, we conducted a comprehensive analysis and identified a total of 120 ACP genes across six species in the Brassica genus. Among these, there were 27, 26, and 30 ACP genes in the allotetraploid B. napus, B. juncea, and B. carinata, respectively, and 14, 13, and 10 ACP genes in the diploid B. rapa, B. oleracea, and B. nigra, respectively. These ACP genes were further classified into six subclades, each containing conserved motifs and domains. Interestingly, the majority of ACP genes exhibited high conservation among the six species, suggesting that the genome evolution and polyploidization processes had relatively minor effects on the ACP gene family. The duplication modes of the six Brassica species were diverse, and the expansion of most ACPs in Brassica occurred primarily through dispersed duplication (DSD) events. Furthermore, most of the ACP genes were under purifying selection during the process of evolution. Subcellular localization experiments demonstrated that ACP genes in Brassica species are localized in chloroplasts and mitochondria. Cis-acting element analysis revealed that most of the ACP genes were associated with various abiotic stresses. Additionally, RNA-seq data revealed differential expression levels of BnaACP genes across various tissues in B. napus, with particularly high expression in seeds and buds. qRT-PCR analysis further indicated that BnaACP genes play a significant role in salt stress tolerance. These findings provide a comprehensive understanding of ACP genes in Brassica plants and will facilitate further functional analysis of these genes.

Genome-wide identification analysis of the 4-Coumarate: CoA ligase (4CL) gene family expression profiles in Juglans regia and its wild relatives J. Mandshurica resistance and salt stress

Persian walnut (Juglans regia) and Manchurian walnut (Juglans mandshurica) belong to Juglandaceae, which are vulnerable, temperate deciduous perennial trees with high economical, ecological, and industrial values. 4-Coumarate: CoA ligase (4CL) plays an essential function in plant development, growth, and stress. Walnut production is challenged by diverse stresses, such as salinity, drought, and diseases. However, the characteristics and expression levels of 4CL gene family in Juglans species resistance and under salt stress are unknown. Here, we identified 36 Jr4CL genes and 31 Jm4CL genes, respectively. Based on phylogenetic relationship analysis, all 4CL genes were divided into three branches. WGD was the major duplication mode for 4CLs in two Juglans species. The phylogenic and collinearity analyses showed that the 4CLs were relatively conserved during evolution, but the gene structures varied widely. 4CLs promoter region contained multiply cis-acting elements related to phytohormones and stress responses. We found that Jr4CLs may be participated in the regulation of resistance to anthracnose. The expression level and some physiological of 4CLs were changed significantly after salt treatment. According to qRT-PCR results, positive regulation was found to be the main mode of regulation of 4CL genes after salt stress. Overall, J. mandshurica outperformed J. regia. Therefore, J. mandshurica can be used as a walnut rootstock to improve salt tolerance. Our results provide new understanding the potential functions of 4CL genes in stress tolerance, offer the theoretical genetic basis of walnut varieties adapted to salt stress, and provide an important reference for breeding cultivated walnuts for stress tolerance.

Genome-wide identification of GATA transcription factors in tetraploid potato and expression analysis in differently colored potato flesh

The GATA gene family belongs to a kind of transcriptional regulatory protein featuring a zinc finger motif, which is essential for plant growth and development. However, the identification of the GATA gene family in tetraploid potato is still not performed. In the present research, a total of 88 GATA genes in the tetraploid potato C88.v1 genome were identified by bioinformatics methods. These StGATA genes had an uneven distribution on 44 chromosomes, and the corresponding StGATA proteins were divided into four subfamilies (I-IV) based on phylogenetic analysis. The cis-elements of StGATA genes were identified, including multiple cis-elements related to light-responsive and hormone-responsive. The collinearity analysis indicates that segmental duplication is a key driving force for the expansion of GATA gene family in tetraploid potato, and that the GATA gene families of tetraploid potato and Arabidopsis share a closer evolutionary relationship than rice. The transcript profiling analysis showed that all 88 StGATA genes had tissue-specific expression, indicating that the StGATA gene family members participate in the development of multiple potato tissues. The RNA-seq analysis was also performed on the tuber flesh of two potato varieties with different color, and 18 differentially expressed GATA transcription factor genes were screened, of which eight genes were validated through qRT-PCR. In this study, we identified and characterized StGATA transcription factors in tetraploid potato for the first time, and screened differentially expressed genes in potato flesh with different color. It provides a theoretical basis for further understanding the StGATA gene family and its function in anthocyanin biosynthesis.