Gene Duplication Research Articles

Pangenome – its Aspect and Prospect

Pangenome is a very new discipline showing a collection of unique and variable genomes of any species in one model. This discipline is a combination of three subjects of Biology like Computer Science and Applied Mathematics. All genes from all species or strains are called Pangenomes. Pangenome was first discovered during the preparation of protein-based vaccine of the G group of Streptococcus sp. which is the main cause of neonatal life-threatening disease by a team of scientists led by Tettelin. Conventional method of Vaccine preparation was not successful due to the existence of five different serotypes in this organism. Then Tettelin with groups of scientists tried to use the genome sequence information for preparing protein- based vaccine in G group of Streptococcus. During their studies, two gap free genomes were found in this G group of bacteria. After sequencing studies of many strains, they found several regions of diversity (genomic diversity) among different isolates of same strain. After sequencing a large number of genomes, they realized that genes belong to isolates present in the species but not present in genomes. Then the Pangenome concept came to their mind and genes were classified into Core genome, Accessory or Dispensable genome and Species-specific genomes. Bacteria genomes have a distinct number of genes but there are differences in the genomes of different strains of a single species that is mainly due to horizontal gene transfer of genes from other bacteria. The Pangenome concept and its discovery was discussed. Plant genomes are dynamic with whole genome duplication, tandem duplication, transposable element duplication, chromosomal rearrangements or structural variations like deletion, inversion, translocation and recombination between species. Several methods have already been tried for the improvement of crop plants like Marker Assisted Selection (MAS), QTL mapping, Single Nucleotide Polymorphisms (SNP), Copy number variations and Genome Wide Association Studies (GWAS) by comparing with Reference Genome Assemblies representing the genome sequence of the genome of an individual or organism. But the Reference genome is not able to locate all genetic variations of a species. To overcome this limitation the concept of Pangenomes may help in the crop improvement which has been discussed. Pangenomic studies have already been done in many crop plants like Rice, Wheat, Maize, Brassica, Soybean etc. Some of all these aspects have been discussed.

Targeting rapid TKI-induced AXL upregulation overcomes adaptive ERK reactivation and exerts antileukemic effects in FLT3/ITD acute myeloid leukemia.

Acute myeloid leukemia (AML) patients with the FMS-related receptor tyrosine kinase 3 internal tandem duplication (FLT3/ITD) mutation have a poorer prognosis, and treatment with FLT3 tyrosine kinase inhibitors (TKIs) has been hindered by resistance mechanisms. One such mechanism is known as adaptive resistance, in which downstream signaling pathways are reactivated after initial inhibition. Past work has shown that FLT3/ITD cells undergo adaptive resistance through the reactivation of extracellular signal-regulated kinase (ERK) signaling within 24 h of sustained FLT3 inhibition. We investigated the mechanism(s) responsible for this ERK reactivation and hypothesized that targeting tyrosine-protein kinase receptor UFO (AXL), another receptor tyrosine kinase that has been implicated in cancer resistance, may overcome the adaptive ERK reactivation. Experiments revealed that AXL is upregulated and activated in FLT3/ITD cell lines mere hours after commencing TKI treatment. AXL inhibition combined with FLT3 inhibition to decrease the ERK signal rebound and to exert greater anti-leukemia effects than with either treatment alone. Finally, we observed that TKI-induced AXL upregulation occurs in patient samples, and combined inhibition of both AXL and FLT3 increased efficacy in our in vivo models. Taken together, these data suggest that AXL plays a role in adaptive resistance in FLT3/ITD AML and that combined AXL and FLT3 inhibition might improve FLT3/ITD AML patient outcomes.

Genome-wide analysis of glutathione S-transferase genes in four Prunus species and the function of PmGSTF2, activated by PmMYBa1, in regulating anthocyanin accumulation in Prunus mume

Glutathione S-transferases (GSTs) are proteases with multiple physiological functions and play an important role in plant responses to abiotic stresses. Nevertheless, there is a paucity of systematic research on GST genes in Prunus genus. Here, 330 GST genes in four Prunus species were identified for the first time and classified into eight subgroups based on protein sequence and conserved structure, among which Tau subfamily genes had the largest number. The amino acid lengths of GST-encoded proteins in the four species ranged from 66 to 1152 aa, most of which were soluble proteins and located in the cytoplasm and chloroplasts. The GST family was propelled by tandem duplications, yet robust purifying selection constrained its divergence. Conserved motif and domain analysis revealed that the majority of PmGSTs exhibited a highly conserved GST-N structure. The expression pattern of PmGSTs exhibited tissue specificity and spatiotemporal specificity. qRT-PCR validated the transcriptome results and 11 genes were differentially expressed in varieties with different flower and stem colors. In addition, we discovered an anthocyanin-related gene PmGSTF2, which can effectively restore the anthocyanin and proanthocyanidin deficiency-related phenotypes of the Arabidopsis tt19 mutant. Recombinant PmGSTF2 enhanced the water solubility of cyanidin and cyanidin-3-O-glucoside in vitro. Moreover, PmMYBa1 could directly bind to the promoter of PmGSTF2 and activate its expression. The findings revealed that GSTs were preserved in Prunus species and that PmGSTF2 was critical in regulating anthocyanin accumulation.

Graph-based self-supervised learning for repeat detection in metagenomic assembly.

Repetitive DNA (repeats) poses significant challenges for accurate and efficient genome assembly and sequence alignment. This is particularly true for metagenomic data, in which genome dynamics such as horizontal gene transfer, gene duplication, and gene loss/gain complicate accurate genome assembly from metagenomic communities. Detecting repeats is a crucial first step in overcoming these challenges. To address this issue, we propose GraSSRep, a novel approach that leverages the assembly graph's structure through graph neural networks (GNNs) within a self-supervised learning framework to classify DNA sequences into repetitive and nonrepetitive categories. Specifically, we frame this problem as a node classification task within a metagenomic assembly graph. In a self-supervised fashion, we rely on a high-precision (but low-recall) heuristic to generate pseudolabels for a small proportion of the nodes. We then use those pseudolabels to train a GNN embedding and a random forest classifier to propagate the labels to the remaining nodes. In this way, GraSSRep combines sequencing features with predefined and learned graph features to achieve state-of-the-art performance in repeat detection. We evaluate our method using simulated and synthetic metagenomic data sets. The results on the simulated data highlight GraSSRep's robustness to repeat attributes, demonstrating its effectiveness in handling the complexity of repeated sequences. Additionally, experiments with synthetic metagenomic data sets reveal that incorporating the graph structure and the GNN enhances the detection performance. Finally, in comparative analyses, GraSSRep outperforms existing repeat detection tools with respect to precision and recall.

Genome-Wide Identification and Functional Validation of Actin Depolymerizing Factor (ADF) Gene Family in Gossypium hirsutum L.

The Actin Depolymerizing Factor (ADF) protein, highly conserved among eukaryotes, is essential for plant growth, development, and stress responses. Cotton, a vital economic crop with applications spanning oilseed, textiles, and military sectors, has seen a limited exploration of its ADF gene family. This research has identified 118 unique ADF sequences across four principal cotton species: Gossypium hirsutum L., Gossypium barbadense Linn, Gossypium raimondii, and Asiatic cotton. The study found that the structural domains and physicochemical properties of these proteins are largely uniform across species. The ADF genes were classified into four subfamilies with a notable expansion in groups III and IV due to tandem and chromosomal duplication events. A thorough analysis revealed a high degree of conservation in gene structure, including exon counts and the lengths of introns and exons, with the majority of genes containing three exons, aligning with the characteristics of the ADF family. RNA-seq analysis uncovered a spectrum of responses by GhADFs to various abiotic stresses with GhADF19 showing the most significant reaction. Virus-induced gene silencing (VIGS) experiments were conducted to assess the role of GhADF19 in plant growth under abiotic stress. The results demonstrated that plants with silenced GhADF19 exhibited significantly slower growth rates and lower dry weights when subjected to cold, salt, and drought stress compared to the control group. This marked reduction in growth and dry weight under stress conditions highlights the potential importance of GhADF19 in stress tolerance mechanisms.

A chromosomal-level genome assembly of Corcyra cephalonica Stainton (Lepidoptera: Pyralidae)

The family Pyralidae is a large and important group within the superfamily Pyraloidea, consisting of more than 5,900 described species. The rice moth Corcyra cephalonica is a noteworthy pest affecting stored grains and products. However, the absence of high-quality genomes hinders our comprehension of the ecology of C. cephalonica. Herein, we propose a chromosome-level genome of C. cephalonica achieved through integrating Nanopore, Illumina, and Hi-C reads. The genome of 453.37 Mb in size consists of 98 scaffolds and 154 contigs, with a scaffold N50 size of 16.67 Mb. The majority of contigs (99.36%, 450.45 Mb) were effectively anchored to 29 chromosomes. The BUSCO analysis of the assembly revealed a completeness of 99.0% (n = 1,367), identifying 98.6% single-copy BUSCO genes and 0.4% duplicated BUSCO genes. Furthermore, the genome contains 48.74% (220.96 Mb) repeat elements and encompasses 14,135 predicted protein-coding genes. This study offers a valuable genomic foundation for advancing our comprehension of the biology of C. cephalonica.

Odorant receptors tuned to isothiocyanates in Drosophila melanogaster and their evolutionary expansion in herbivorous relatives.

Plants release complex volatile compounds to attract mutualists, deter herbivores, and deceive pollinators. Here, we used herbivorous specialist flies that feed on mustard plants ( Scaptomyza spp.) and microbe-feeding species ( Drosophila melanogaster and Scaptomyza spp.) to investigate how plant-derived electrophilic toxins such as isothiocyanates (ITCs) affect insects, and how flies detect these compounds through olfaction. In survival assays, D. melanogaster exposed to volatile allyl isothiocyanate (AITC), a toxin derived from many Brassicales plants, were acutely intoxicated, demonstrating the high toxicity of this volatile compound to non-specialized insects. Through single sensillum recordings (SSR) from olfactory organs and behavioral assays, we found that the Odorant receptor 42a (Or42a) is necessary for AITC detection and behavioral aversion. Comparative transcriptome and RNA FISH studies across the drosophilid genus Scaptomyza revealed lineage-specific triplication of Or42a in the Brassicales specialists and a doubling of Or42a -positive-olfactory sensory neurons. Heterologous expression experiments showed that Or42a paralogs in Brassicales-specialists exhibited broadened sensitivity to ITCs in a paralog-specific manner. Finally, AlphaFold2 modeling followed by site-directed mutagenesis and SSR identified two critical amino acid substitutions that conferred Or42a heighten sensitivity to Brassicales-derived ITCs. Our findings suggest that ITCs, which are toxic to most insects, can be detected and avoided by non-specialists like D. melanogaster through olfaction. In Brassicales specialists, these same Ors experienced gene duplication events that resulted in an expanded sensitivity to ITC compounds. Thus, the insect's olfactory system can rapidly adapt to toxic ecological niches provided by chemically-defended host plants through co-option of chemosensory capabilities already present in their ancestors.

Genome-wide identification and expression analysis of the ZRT, IRT-like protein (ZIP) family in Nicotiana tabacum.

Iron (Fe) and Zinc (Zn) are essential micronutrients for plant growth and development. ZIP (ZRT, IRT-like protein) transporters, known for their role in the regulation of Zinc and Iron uptake, are pivotal in facilitating the absorption, transport, and maintenance of Fe/Zn homeostasis in plants. Nicotiana tabacum has been widely used as a model plant for gene function analysis; however, the tobacco ZIP genes have not been identified systematically. In this study, we have identified a comprehensive set of 32 NtZIP genes, which were phylogenetically categorized into three distinct clades. The gene structures, characterized by their exon/intron organization, and the protein motifs are relatively conserved, particularly among genes within the same clade. These NtZIP genes exhibit an uneven distribution across 12 chromosomes. The gene localization analysis revealed the presence of 11 pairs of homeologous locus genes and 7 pairs of tandem duplication genes within the genome. To further explore the functionality of these genes, real-time quantitative reverse transcription PCR was employed to assess their expression levels in roots subjected to metal deficiency. The results indicated that certain NtZIP genes are specifically upregulated in response to either Fe or Zn deficiency. Additionally, the presence of specific cis-elements within their promoter regions, such as the E-box associated with Fe deficiency response and the ZDRE box linked to Zn deficiency response, was identified. This study lays a foundational groundwork for future research into the biological functions of NtZIP genes in tobacco in micronutrient regulation and homeostasis.

Hidden evolutionary constraints dictate the retention of coronavirus accessory genes.

Coronaviruses exhibit many mechanisms of genetic innovation, including the acquisition of accessory genes that originate by capture of cellular genes or through duplication of existing viral genes. Accessory genes influence viral host range and cellular tropism, but little is known about how selection acts on these variable regions of virus genomes. We used experimental evolution of mouse hepatitis virus (MHV) encoding a cellular AKAP7 phosphodiesterase and an inactive native phosphodiesterase, NS2 to model the evolutionary fate of accessory genes. After courses of serial infection, the gene encoding inactive NS2, ORF2, unexpectedly remained intact, suggesting it is under cryptic constraint uncoupled from the function of NS2. In contrast, AKAP7 was retained under strong selection but rapidly lost under relaxed selection. Experimental evolution also led to altered viral replication in a cell type-specific manner and changed the relative proportions of subgenomic viral RNA in plaque-purified viral isolates, revealing additional mechanisms of adaptation. Guided by the retention of ORF2 and similar patterns in related betacoronaviruses, we analyzed ORF8 of SARS-CoV-2, which arose via gene duplication and contains premature stop codons in several globally successful lineages. As with MHV ORF2, the coding-defective SARS-CoV-2 ORF8 gene remains largely intact, mirroring patterns observed during MHV experimental evolution, challenging assumptions on the dynamics of gene loss in virus genomes and extending these findings to viruses currently adapting to humans.

Challenges in Assembling the Dated Tree of Life.

The assembly of a comprehensive and dated Tree of Life (ToL) remains one of the most formidable challenges in evolutionary biology. The complexity of life's history, involving both vertical and horizontal transmission of genetic information, defies its representation by a simple bifurcating phylogeny. With the advent of genome and metagenome sequencing, vast amounts of data have become available. However, employing this information for phylogeny and divergence time inference has introduced significant theoretical and computational hurdles. This perspective addresses some key methodological challenges in assembling the dated ToL, namely, the identification and classification of homologous genes, accounting for gene tree-species tree mismatch due to population-level processes along with duplication, loss, and horizontal gene transfer, and the accurate dating of evolutionary events. Ultimately, the success of this endeavor requires new approaches that integrate knowledge databases with optimized phylogenetic algorithms capable of managing complex evolutionary models.

Genetic inheritance studies in rice (Oryza sativa L.) for grain protein, quality and yield

The present study was undertaken to investigate the genetic parameters of yield component and quality traits of the F3 population from BPT 5204 and JAK-686 cross in rice. High heritability and high genetic advance as a percentage of the mean were recorded for productive tillers per plant, grains per panicle, panicle length and grain yield per plant indicating the effectiveness of direct phenotypic selection for improvement of these traits. For traits such as grains per panicle, 1000 grain weight, amylose content and grain protein content, the magnitude of additive (d) effects surpassed dominance (h) gene effects. This implies that phenotypic selection using methods like pureline selection, mass selection would be beneficial. Gel consistency exhibited complimentary gene interaction. Furthermore, duplicate gene interaction was prevalent in traits like days to 50% flowering, productive tillers per plant, plant height, panicle length, grain yield per plant, kernel length, kernel breadth, and kernel L/B ratio, making selection in later generations more effective. The inheritance studies highlighted the importance of selecting desired recombinants from the segregating population to enhance rice quality as well as grain yield. Keywords: Rice, Grain protein content, Quality, Variability, Heritability, Gene action

Genome-wide analysis of the WRKY gene family and their response to low-temperature stress in elephant grass.

Elephant grass (Cenchrus purpureus) is a perennial forage grass characterized by tall plants, high biomass and wide adaptability. Low-temperature stress severely limits elephant grass biomass and geographic distribution. WRKY is one of the largest families of plant-specific transcription factors and plays important roles in plant resistance to low-temperature. However, the understanding of the WRKY family in grasses is limited. In this study, we conducted a genome-wide characterization of WRKY proteins in elephant grass, including gene structure, phylogeny, expression, conserved motif organization, and functional annotation, to identify key CpWRKY candidates involved in cold tolerance. In this study, a total of 176 WRKY genes were identified in elephant grass. It was found that 172 were unevenly distributed across its 14 chromosomes, while the remaining 4 genes were not anchored to any chromosome. The genes were classified into three groups based on their WRKY conserved domains and zinc finger motifs. There were 12, 8, 19, 27, 12, 18 and 80 CpWRKYs belonging to group I, group IIa, group IIb, group IIc, group IId, group IIe and group III, respectively. We hypothesized that the ancient subgroup IIc WRKY gene is the ancestor of all WRKY genes in elephant grass. Most CpWRKYs in the same group have similar structure and motif composition. A total of 169 duplicate gene pairs were identified, suggesting that segmental duplication might have contributed to the expansion of the CpWRKY gene family. Ka/Ks analysis revealed that most of the CpWRKYs were subjected to purifying selection during the evolution. It was also found that six genes (CpWRKY51, CpWRKY81, CpWRKY100, CpWRKY101, CpWRKY140 and CpWRKY143) exhibited higher expression in roots compare to leaves, and were significantly induced by low temperature stress. Among them, CpWRKY81 had the highest expression under low-temperature stress, and its over-expression significantly enhanced the cold tolerance in yeast. In this study, we characterized WRKY genes in elephant grass and further investigated their physicochemical properties, evolution, and expression patterns under low-temperature stress. This research provides valuable resources for identifying key CpWRKY genes that contribute to cold tolerance in elephant grass.

Genome-Wide Identification and Expression Analysis of the BTB Gene Superfamily Provides Insight into Sex Determination and Early Gonadal Development of Alligator sinensis.

The BTB gene superfamily is widely distributed among higher eukaryotes and plays a significant role in numerous biological processes. However, there is limited knowledge about the structure and function of BTB genes in the critically endangered species Alligator sinensis, which is endemic to China. A total of 170 BTB genes were identified from the A. sinensis genome, classified into 13 families, and unevenly distributed across 16 chromosomes. Analysis of gene duplication events yielded eight pairs of tandem duplication genes and six pairs of segmental duplication genes. Phylogenetics shows that the AsBTB genes are evolutionarily conserved. The cis-regulatory elements in the AsBTB family promoter region reveal their involvement in multiple biological processes. Protein interaction network analysis indicates that the protein interactions of the AsBTB genes are centered around CLU-3, mainly participating in the regulation of biological processes through the ubiquitination pathway. The expression profile and protein interaction network analysis of AsBTB genes during sex differentiation and early gonadal development indicate that AsBTB genes are widely expressed in this process and involves numerous genes and pathways for regulation. This study provides a basis for further investigation of the role of the BTB gene in sex differentiation and gonadal development in A. sinensis.

The distribution of beneficial mutational effects between two sister yeast species poorly explains natural outcomes of vineyard adaptation.

Domesticated strains of Saccharomyces cerevisiae have adapted to resist copper and sulfite, two chemical stressors commonly used in winemaking. S. paradoxus has not adapted to these chemicals despite being consistently present in sympatry with S. cerevisiae in vineyards. This contrast could be driven by a number of factors including niche differences or differential access to resistance mutations between species. In this study, we used a comparative mutagenesis approach to test whether S. paradoxus is mutationally constrained with respect to acquiring greater copper and sulfite resistance. For both species, we assayed the rate, effect size, and pleiotropic costs of resistance mutations and sequenced a subset of 150 mutants. We found that the distributions of mutational effects displayed by the two species were similar and poorly explained the natural pattern. We also found that chromosome VIII aneuploidy and loss of function mutations in PMA1 confer copper resistance in both species, whereas loss of function mutations in REG1 were only a viable route to copper resistance in S. cerevisiae. We also observed a de novo duplication of the CUP1 gene in S. paradoxus but not in S. cerevisiae. For sulfite, loss of function mutations in RTS1 and KSP1 confer resistance in both species, but mutations in RTS1 have larger effects in S. paradoxus. Our results show that even when available mutations are largely similar, species can differ in the adaptive paths available to them. They also demonstrate that assays of the distribution of mutational effects may lack predictive insight concerning adaptive outcomes.

Are there conserved biosynthetic genes in lichens? Genome-wide assessment of terpene biosynthetic genes suggests ubiquitous distribution of the squalene synthase cluster

Lichen-forming fungi (LFF) are prolific producers of functionally and structurally diverse secondary metabolites, most of which are taxonomically exclusive and play lineage-specific roles. To date, widely distributed, evolutionarily conserved biosynthetic pathways in LFF are not known. However, this idea stems from polyketide derivatives, since most biochemical research on lichens has concentrated on polyketide synthases (PKSs). Here, we present the first systematic identification and comparison of terpene biosynthetic genes of LFF using all the available Lecanoromycete reference genomes and 22 de novo sequenced ones (111 in total, representing 60 genera and 23 families). We implemented genome mining and gene networking approaches to identify and group the biosynthetic gene clusters (BGCs) into networks of similar BGCs. Our large-scale analysis led to the identification of 724 terpene BGCs with varying degrees of pairwise similarity. Most BGCs in the dataset were unique with no similarity to a previously known fungal or bacterial BGC or among each other. Remarkably, we found two BGCs that were widely distributed in LFF. Interestingly, both conserved BGCs contain the same core gene, i.e., putatively a squalene/phytoene synthase (SQS), involved in sterol biosynthesis. This indicates that early gene duplications, followed by gene losses/gains and gene rearrangement are the major evolutionary factors shaping the composition of these widely distributed SQS BGCs across LFF. We provide an in-depth overview of these BGCs, including the transmembrane, conserved, variable and LFF-specific regions. Our study revealed that lichenized fungi do have a highly conserved BGC, providing the first evidence that a biosynthetic gene may constitute essential genes in lichens.

Dosage constraint of the ribosome-associated molecular chaperone drives the evolution and fates of its duplicates in bacteria.

Gene duplication events presumably occur in tig, which encodes the ribosome-associated molecular chaperone trigger factor (TF). However, TF is singly copied in virtually every bacterium, and these exceptions with multiple TF homologs always retain only one complete copy while other homologs lack the N-terminal domain. Here, we reveal the manner and mechanism underlying the evolution and fates of TF duplicates in bacteria. We discovered that the mutation-to-loss or retention-to-sub/neofunctionalization of TF duplicates is associated with the dosage constraint of N-terminal complete TF. The dosage constraint of TF is attributed to its characteristic ribosome binding and substrate-holding activities, causing a decrease in protein productivity and profile changes in cellular proteome.

Comparative analysis of mitochondrial genomes of invasive weed Mikania micrantha and its indigenous congener Mikania cordata

Mikania micrantha and Mikania cordata are two distinct species in China. The former is notorious as one of the top 100 worst invasive species, whereas the latter is an indigenous species harmless to native plants or the environment. They form an ideal congener pair for comparative studies aimed at deeply understanding the invasion mechanisms of the exotic weed. In this study, we have assembled and annotated the mitogenomes of both species using Illumina and PacBio sequencing data and compared their characteristic differences. The complete mitogenome of M. micrantha is a double-stranded DNA with a length of 336,564 bp, while the mitogenome of M. cordata exhibits a branching structure, consisting of two small circular molecules and six linear molecules, with a combined length totaling 335,444 bp. Compared to M. cordata, M. micrantha has less SSRs, tandem repeats, dispersed repeats, mitochondrial protein coding genes (PCGs). The two plants show similar codon usage patterns. This comparative study has revealed the structure and function of the mitogenomes of the two species and laid a solid foundation for investigating the effects of gene loss and duplication on the development of invasive traits in M. micrantha.

Genome-Wide Analysis of Potato CCT Family Genes and Its Response to Auxin Substances

The control of flowering time plays an important role in the growth and development of potato tubers. The CCT (CO, COL and TOC1) gene family is involved in the flowering process of plants. In this study, a total of 32 StCCT family genes were identified and further classified into five subfamilies, including COL (17 members), PRR (4 members), ZIM (3 members), ASML2 (6 members) and TCR1 (2 members), based on their phylogenetic relationship. An analysis of the gene structure, motif compositions and conserved domain provided support for this classification. The StCCT genes were unevenly distributed on 12 chromosomes of the potato plant. In total, six gene duplication events were observed, which played a crucial role in the expansion of the StCCT family genes in the potato. The expression profiles exhibited diverse expression patterns of the StCCT genes in six tissues (leaf, shoot, root, tuber, stolon, and flower), StCCT32 is only expressed in flowers, while StCCT19 and StCCT8 are highly expressed in flowers and tubers, respectively. The StCCT genes exhibit different expression patterns in response to IAA and TIBA treatments at different concentrations across three tissues (leaf, stem, and tuber). After IAA and TIBA treatments, it was found that the expression level of StCCT7 was low in leaves and stems but significantly increased in tubers. Collectively, this study provided valuable information for the further study of potato formation and development and provided candidate genes for molecular breeding in the potato.

Systematic analysis of the ARF gene family in Fagopyrum dibotrys and its potential roles in stress tolerance.

The auxin response factor (ARF) is a plant-specific transcription factor that regulates the expression of auxin response genes by binding directly to their promoters. They play an important role in the regulation of plant growth and development, as well as in the response to biotic and abiotic stresses. However, the identification and functional analysis of ARFs in Fagopyrum dibotrys are still unclear. In this study, a total of 26 FdARF genes were identified using bioinformatic methods. Their chromosomal location, gene structure, physical and chemical properties of their encoded protein, subcellular location, phylogenetic tree, conserved motifs and cis-acting elements in FdARF promoters were analyzed. The results showed that 26 FdARF genes were unevenly distributed on 8 chromosomes, with the largest distribution on chromosome 4 and the least distribution on chromosome 3. Most FdARF proteins are located in the nucleus, except for the proteins FdARF7 and FdARF21 located to the cytoplasm and nucleus, while FdARF14, FdARF16, and FdARF25 proteins are located outside the chloroplast and nucleus. According to phylogenetic analysis, 26 FdARF genes were divided into 6 subgroups. Duplication analysis indicates that the expansion of the FdARF gene family was derived from segmental duplication rather than tandem duplication. The prediction based on cis-elements of the promoter showed that 26 FdARF genes were rich in multiple stress response elements, suggesting that FdARFs may be involved in the response to abiotic stress. Expression profiling analysis showed that most of the FdARF genes were expressed in the roots, stems, leaves, and tubers of F. dibotrys, but their expression exhibits a certain degree of tissue specificity. qRT-PCR analysis revealed that most members of the FdARF gene were up- or down-regulated in response to abiotic stress. The results of this study expand our understanding of the functional role of FdARFs in response to abiotic stress and lay a theoretical foundation for further exploration of other functions of FdARF genes.

Evolution and functional divergence of glycosyltransferase genes shaped the quality and cold tolerance of tea plants.

Plant glycosyltransferases (UGTs) play a key role in plant growth and metabolism. Here, we examined the evolutionary landscape among UGTs in 28 fully sequenced species from early algae to angiosperms. Our findings revealed a distinctive expansion and contraction of UGTs in the G and H groups in tea (Camellia sinensis), respectively. Whole-genome duplication and tandem duplication events jointly drove the massive expansion of UGTs, and the interplay of natural and artificial selection has resulted in marked functional divergence within the G group of the sinensis-type tea population. In Cluster II of group G, differences in substrate selection (e.g., Abscisic Acid) of the enzymes encoded by UGT genes led to their functional diversification, and these genes influence tolerance to abiotic stresses such as low temperature and drought via different modes of positive and negative regulation, respectively. UGTs in Cluster III of the G group have diverse aroma substrate preferences, which contributes a diverse aroma spectrum of the sinensis-type tea population. All Cluster III genes respond to low-temperature stress, whereas UGTs within Cluster III-1, shaped by artificial selection, are unresponsive to drought. This suggests that artificial selection of tea plants focused on improving quality and cold tolerance as primary targets.