Duplicated Gene Pairs Research Articles

Genome-Wide Identification of the bHLH Gene Family in Callerya speciosa Reveals Its Potential Role in the Regulation of Isoflavonoid Biosynthesis

Callerya speciosa (Champ. ex Benth.) Schot is a significant leguminous plant valued for its edible tuberous roots, which are a plentiful source of isoflavonoids. Basic helix–loop–helix (bHLH) transcription factors (TFs) have been reported to regulate secondary metabolism in plants, especially flavonoid biosynthesis. However, the bHLH genes in C. speciosa have not yet been reported, and their regulatory role in isoflavonoid biosynthesis remains unexplored. Here, 146 CsbHLH genes were identified in the C. speciosa genome, classifying them into 23 subfamilies based on the gene structures and phylogenetic relationships. All the CsbHLH proteins contained both motifs 1 and 2, whereas motif 8 was only distributed in subgroup III (d + e). Collinearity analysis demonstrated that fragmental replications are the primary driver of CsbHLH evolution, with the majority of duplicated CsbHLH gene pairs experiencing selective pressure. Nine candidate CsbHLH genes were found to play a potential role in regulating isoflavonoid biosynthesis through a combination of gene-to-metabolite correlation analysis and weighted gene co-expression network analysis (WGCNA). Additionally, the cis-regulatory elements and response to MeJA of these nine genes were characterized and confirmed through quantitative real-time PCR (qRT-PCR) analysis. Among them, three CsbHLHs (CsbHLH9, CsbHLH89, and CsbHLH95) were selected for further investigation. Yeast two-hybrid (Y2H), dual-luciferase (LUC) assays, bimolecular fluorescence complementation (BiFC) assays, and transient transformation demonstrated that CsbHLH9 acted as a transcriptional activator through its interaction with CsMYB36 and binding to the promoters of isoflavonoid biosynthesis genes in a MeJA-induced manner, such as CsIFR2, CsI3′H2, and CsCHS4, to promote isoflavonoid (calycosin, calycosin-7-o-glucoside, and formononetin) accumulation. Our results establish a basis for the functional analysis of bHLH genes and investigations into the molecular mechanisms underlying isoflavonoid biosynthesis in C. speciosa.

Genome-Wide Identification, Characterization and Expression Patterns of the DBB Transcription Factor Family Genes in Wheat

Double B-box (DBB) proteins are plant-specific transcription factors (TFs) that play crucial roles in plant growth and stress responses. This study investigated the classification, structure, conserved motifs, chromosomal locations, cis-elements, duplication events, expression levels, and protein interaction network of the DBB TF family genes in common wheat (Triticum aestivum L.). In all, twenty-seven wheat DBB genes (TaDBBs) with two conserved B-box domains were identified and classified into six subgroups based on sequence features. A collinearity analysis of the DBB family genes among wheat, Arabidopsis, and rice revealed some duplicated gene pairs and highly conserved genes in wheat. An expression pattern analysis indicated that wheat TaDBBs were involved in plant growth, responses to drought stress, light/dark, and abscisic acid treatment. A large number of cis-acting regulatory elements related to light response are enriched in the predicted promoter regions of 27 TaDBBs. Furthermore, some of TaDBBs can interact with COP1 or HY5 based on the STRING database prediction and yeast two-hybrid (Y2H) assay, indicating the potential key roles of TaDBBs in the light signaling pathway. Conclusively, our study revealed the potential functions and regulatory mechanisms of TaDBBs in plant growth and development under drought stress, light, and abscisic acid.

In silico analysis of R2R3-MYB transcription factors in the basal eudicot model, Aquilegia coerulea

R2R3-MYBs are an important group of transcription factors that regulate crucial developmental processes across the plant kingdom; yet no comprehensive analysis of the R2R3-MYBs in the early-diverging eudicot clade of Ranunculaceae has been conducted so far. In the present study, Aquilegia coerulea is chosen to understand the extent of conservation and divergence of R2R3-MYBs as a representative of the family by analysing the genomic distribution, organization, gene structure, physiochemical properties, protein architecture, evolution and possible mode of expansion. Genome-wide analysis showed the presence of 82 putative homologues classified into 21 subgroups, based on phylogenetic analysis of full-length protein sequences. The domain has remained largely conserved across all homologues with few differences from the characterized Arabidopsis thaliana R2R3-MYBs. The topology of the phylogenetic tree remains the same when full-length protein sequences are used, indicating that the evolution of R2R3-MYBs is driven by the domain region only. This is supported by the presence of similar structures of exon–intron and conserved motifs within the same subgroup. Furthermore, comparisons of the AqcoeR2R3-MYB members with monocots and core-eudicots revealed the evolutionary expansion of a few functional clades, such as A. thaliana R2R3-MYB subgroup 6 (SG6), the upstream regulatory factors of floral pigment biosynthesis and floral color. The reconstructed evolutionary history of SG6-like genes across angiosperms highlights the occurrence of independent duplication events in the genus Aquilegia. AqcoeR2R3-MYB genes are present in all seven chromosomes of A. coerulea, most of which result from local and segmental duplications. Selection analysis of these duplicated gene pairs indicates purifying selection except one, and the physiochemical analyses of R2R3-MYBs reveal differences among the MYBs signifying their functional diversification. This study paves the way for further investigation of paralogous copies and their probable role in the evolution of different floral traits in A. coerulea. It lays the foundation for functional genomic studies of R2R3-MYBs in the basal eudicots and facilitates comparative studies among angiosperms. The work also provides a framework for deciphering novel genetic regulatory pathways that govern the diversity of floral morphology.

Genome-Wide Analysis of the MAPKKK Gene Family Under Abiotic Stresses in Moso Bamboo (Phyllostachys edulis)

Mitogen-activated protein kinase kinase kinases (MAPKKKs) are the upstream components of MAPK cascades and are involved in mediating stress responses and developmental processes. Although MAPKKK genes have been investigated in many plants, the identification and characterization of MAPKKKs in moso bamboo were still limited. Here, 134 MAPKKKs were identified as unevenly distributed on 23 chromosomes (except for chromosome 1) of moso bamboo and divided into three subfamilies by phylogenetic analysis. The gene structure and conserved motif of PeMAPKKKs were investigated. The expansions of PeMAPKKKs were driven by whole-genome duplication (WGD) or segmental duplication events. The duplicated gene pairs were under purifying selection based on the Ka/Ks ratios, suggesting they underwent functional conservation. Most PeMAPKKKs contained cis-elements related to development, hormones, and stress responses. Tissue expression patterns showed that PeMAPKKKs had multiple expression patterns. The qPCR analysis showed distinct expression patterns of PeMAPKKKs under drought, salt, and cold stress conditions. Taken together, this study provides a solid foundation for future functional characterizations of MAPKKKs and identifies candidate stress-responsive genes for further study in moso bamboo.

Genome-Wide Identification of the Cyclic Nucleotide-Gated Ion Channel Gene Family and Expression Profiles Under Low-Temperature Stress in Luffa cylindrica L.

Cyclic nucleotide-gated ion channels (CNGCs) are cell membrane channel proteins for calcium ions. They have been reported to play important roles in survival and in the responses to environmental factors in various plants. However, little is known about the CNGC family and its functions in luffa (Luffa cylindrica L.). In this study, a bioinformatics-based method was used to identify members of the CNGC gene family in L. cylindrica. In total, 20 LcCNGCs were detected, and they were grouped into five subfamilies (I, II, Ⅲ, IV-a, and IV-b) in a phylogenetic analysis with CNGCs from Arabidopsis thaliana (20 AtCNGCs) and Momordica charantia (17 McCNGCs). The 20 LcCNGC genes were unevenly distributed on 11 of the 13 chromosomes in luffa, with none on Chromosomes 1 and 5. The members of each subfamily encoded proteins with highly conserved functional domains. An evolutionary analysis of CNGCs in luffa revealed three gene losses and a motif deletion. An examination of gene replication events during evolution indicated that two tandemly duplicated gene pairs were the primary driving force behind the evolution of the LcCNGC gene family. PlantCARE analyses of the LcCNGC promoter regions revealed various cis-regulatory elements, including those responsive to plant hormones (abscisic acid, methyl jasmonate, and salicylic acid) and abiotic stresses (light, drought, and low temperature). The presence of these cis-acting elements suggested that the encoded CNGC proteins may be involved in stress responses, as well as growth and development. Transcriptome sequencing (RNA-seq) analyses revealed tissue-specific expression patterns of LcCNGCs in various plant parts (roots, stems, leaves, flowers, and fruit) and the upregulation of some LcCNGCs under low-temperature stress. To confirm the accuracy of the RNA-seq data, 10 cold-responsive LcCNGC genes were selected for verification by quantitative real-time polymerase chain reaction (RT-qPCR) analysis. Under cold conditions, LcCNGC4 was highly upregulated (>50-fold increase in its transcript levels), and LcCNGC3, LcCNGC6, and LcCNGC13 were upregulated approximately 10-fold. Our findings provide new information about the evolution of the CNGC family in L. cylindrica and provide insights into the functions of the encoded CNGC proteins.

Genome-wide identification and analysis of anthocyanin synthesis-related R2R3-MYB genes in Fragaria pentaphylla

BackgroundMYB transcription factors regulate anthocyanin biosynthesis across numerous plant species. However, comprehensive genome-wide investigations regarding the R2R3-MYB gene family and its involvement in regulating anthocyanin biosynthesis in the red and white fruit color morphs of Fragaria pentaphylla remain scarce.ResultsA total of 101 FpR2R3-MYB genes were identified from the F. pentaphylla genome and were divided into 34 subgroups based on phylogenetic analysis. Gene structure (exon/intron) and protein motifs were particularly conserved among the FpR2R3-MYB genes, especially members within the same subgroup. The FpR2R3-MYB genes were distributed over seven F. pentaphylla chromosomes. Analysis of gene duplication events revealed five pairs of tandem duplication genes and 16 pairs of segmental duplication genes, suggesting that segmental duplications are the major pattern for expansion of the FpR2R3-MYB gene family expansion in F. pentaphylla. Cis-regulatory elements of the FpR2R3-MYB promoters were involved in cellular development, phytohormones, environmental stress and photoresponse. Based on the analysis of the FpR2R3-MYB gene family and transcriptome sequencing (RNA-seq) data, FpMYB9 was identified as a key transcription factor involved in the regulation of anthocyanin synthesis in F. pentaphylla fruits. The expression of FpMYB9 increases significantly during the ripening stage of red fruits, as confirmed by reverse transcription quantitative real-time PCR. In addition, subcellular localization experiments further confirmed the nuclear presence of FpMYB9, supporting its role as a transcription factor involved in anthocyanin biosynthesis.ConclusionOur results showed that the FpR2R3-MYB genes are highly conserved and play important roles in the anthocyanin biosynthesis in F. pentaphylla fruits. Our results also provide a compelling basis for further understanding of the regulatory mechanism underlying the role of FpMYB9 in anthocyanin formation in F. pentaphylla fruits.

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.

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.

Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms.

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. In this study, we analyzed the RFPKM between homoeologs in 3 angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these 3 species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage-balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in 1 copy of homoeologs is inclined toward nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage-balance constraints, in the evolutionary retention of homoeologs in plants.

Genome-wide investigation of caffeic acid O-methyltransferases and expression analysis under different abiotic stresses and anthracnose infection in pepper (Capsicum annuum L.)

ABSTRACT The caffeic acid O-methyl transferases (COMTs) facilitate the biosynthesis of lignin and play important roles in many processes, including growth, development and response to stresses. Pepper (Capsicum annuum L.) is a major horticultural crop with significant nutritional and medicinal values. Nevertheless, a comprehensive assessment of the COMT genes in pepper and their involvement in drought, salt tolerance and anthracnose resistance remains unestablished. In the present study, 80 pepper COMT genes were identified and classified into four major groups based on domain structure and phylogenetics. Additionally, the sequence analysis, gene structures, and duplication events were analysed. The genes were unevenly distributed across 10 chromosomes, and 10 pairs of tandem and 13 pairs of segmental duplicated genes were detected in the COMT gene family. Examination of cis-acting regulatory elements detected the presence of phytohormone-specific, stress-receptive and light responsive elements in the pepper COMT promoters. Moreover, qRT-PCR analysis revealed significant upregulated expression of CanCOMT48 to drought stress, CanCOMT14 to salinity stress and CanCOMT18 to anthracnose infection by Colletotrichum truncatum, suggesting their important functions in different stresses. Our results provide significant insights into the expression and probable function of selective COMTs which could be used as potential candidate genes for pepper tolerance breeding.

Genome-wide analysis of apple CNGC family allows the identification of MdCNGC15A negatively regulating apple salt tolerance

Calcium (Ca2+) is essential for signal conduction and plant growth. Cyclic nucleotide-gated channels (CNGCs) are Ca2+ transporters that regulate Ca2+ signalling and homeostasis by modulating its transmembrane transport, thereby influencing plant development as well as the biotic and abiotic stress responses. Although identified in numerous plant species, the CNGC family has not been characterized in apple until now. Here, 20 MdCNGCs were identified from the apple genome and were randomly distributed on 13 chromosomes. Phylogenetic analysis classified these MdCNGCs into five groups (I, Ⅱ, Ⅲ, Ⅳ-a, and Ⅳ-b), with five pairs of segmental duplicated genes being detected via collinearity analysis. Sequence alignment and analyses of gene structures, conserved motifs, and 3D structures indicated high structural conservation, particularly within groups. Yeast two-hybrid (Y2H) assays demonstrated interactions between most MdCNGCs and the Ca2+ receptor MdCaM7.1, except for MdCNGC1B and MdCNGC15A. Promoter analysis and expression profiling revealed significant responses to abiotic stress, particularly salt stress, in some MdCNGCs. Silencing MdCNGC15A significantly enhanced apple plants salt tolerance, while its overexpression in apple calli significantly decreased tolerance, as shown by transgenic analysis. Collectively, our results demonstrate the crucial role of MdCNGCs in abiotic stress responses and provide valuable insights for future functional and regulatory studies in apples.

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

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

Genome-wide identification and expression analysis of the SPL gene family and its response to abiotic stress in barley (Hordeum vulgare L.)

BackgroundSquamosa promoter-binding protein-like (SPL) is a plant-specific transcription factor that is widely involved in the regulation of plant growth and development, including flower and grain development, stress responses, and secondary metabolite synthesis. However, this gene family has not been comprehensively evaluated in barley, the most adaptable cereal crop with a high nutritional value.ResultsIn this study, a total of 15 HvSPL genes were identified based on the Hordeum vulgare genome. These genes were named HvSPL1 to HvSPL15 based on the chromosomal distribution of the HvSPL genes and were divided into seven groups (I, II, III, V, VI, VII, and VIII) based on the phylogenetic tree analysis. Chromosomal localization revealed one pair of tandem duplicated genes and one pair of segmental duplicated genes. The HvSPL genes exhibited the highest collinearity with the monocotyledonous plant, Zea mays (27 pairs), followed by Oryza sativa (18 pairs), Sorghum bicolor (16 pairs), and Arabidopsis thaliana (3 pairs), and the fewest homologous genes with Solanum lycopersicum (1 pair). The distribution of the HvSPL genes in the evolutionary tree was relatively scattered, and HvSPL proteins tended to cluster with SPL proteins from Z. mays and O. sativa, indicating a close relationship between HvSPL and SPL proteins from monocotyledonous plants. Finally, the spatial and temporal expression patterns of the 14 HvSPL genes from different subfamilies were determined using quantitative real-time polymerase chain reaction (qRT-PCR). Based on the results, the HvSPL gene family exhibited tissue-specific expression and played a regulatory role in grain development and abiotic stress. HvSPL genes are highly expressed in various tissues during seed development. The expression levels of HvSPL genes under the six abiotic stress conditions indicated that many genes responded to stress, especially HvSPL8, which exhibited high expression under multiple stress conditions, thereby warranting further attention.ConclusionIn this study, 15 SPL gene family members were identified in the genome of Hordeum vulgare, and the phylogenetic relationships, gene structure, replication events, gene expression, and potential roles of these genes in millet development were studied. Our findings lay the foundation for exploring the HvSPL genes and performing molecular breeding of barley.

Genome-Wide Characterization of Alfin-like Genes in Brassica napus and Functional Analyses of BnaAL02 and BnaAL28 in Response to Nitrogen and Phosphorus Deficiency.

Alfin-like proteins (ALs) form a plant-specific transcription factor (TF) gene family involved in the regulation of plant growth and development, and abiotic stress response. In this study, 30 ALs were identified in Brassica napus ecotype 'Zhongshuang 11' genome (BnaALs), and unevenly distributed on 15 chromosomes. Structural characteristic analysis showed that all of the BnaALs contained two highly conserved domains: the N terminal DUF3594 domain and the C-terminal PHD-finger domain. The BnaALs were classified into four groups (Group I-IV), supported by conserved intron-exon and protein motif structures in each group. The allopolyploid event between B. oleracea and B. rapa ancestors and the small-scale duplication events in B. napus both contributed to the large BnaALs expansion. The promoter regions of BnaALs contained multiple abiotic stress cis-elements. The BnaALs in I-IV groups were mainly expressed in cotyledon, petal, root, silique, and seed tissues, and the duplicated gene pairs shared highly similar expression patterns. RNA-seq and RT-qPCR analysis showed that BnaALs were obviously induced by low nitrogen (LN) and low phosphorus (LP) treatments in roots. Overexpressing BnaAL02 and BnaAL28 in Arabidopsis demonstrated their functions in response to LN and LP stresses. BnaAL28 enhanced primary roots' (PRs) length and lateral roots' (LRs) number under LP and LN conditions, where BnaAL02 can inhibit LR numbers under the two conditions. They can promote root hair (RH) elongation under LP conditions; however, they suppressed RH elongation under LN conditions. Our result provides new insight into the functional dissection of this family in response to nutrient stresses in plants.

In Silico Analysis: Molecular Characterization and Evolutionary Study of CLCN Gene Family in Buffalo.

Chloride channels (ClCs) have received global interest due to their significant role in the regulation of ion homeostasis, fluid transport, and electrical excitability of tissues and organs in different mammals and contributing to various functions, such as neuronal signaling, muscle contraction, and regulating the electrolytes' balance in kidneys and other organs. In order to define the chloride voltage-gated channel (CLCN) gene family in buffalo, this study used in silico analyses to examine physicochemical properties, evolutionary patterns, and genome-wide identification. We identified eight CLCN genes in buffalo. The ProtParam tool analysis identified a number of important physicochemical properties of these proteins, including hydrophilicity, thermostability, in vitro instability, and basic nature. Based on their evolutionary relationships, a phylogenetic analysis divided the eight discovered genes into three subfamilies. Furthermore, a gene structure analysis, motif patterns, and conserved domains using TBtool demonstrated the significant conservation of this gene family among selected species over the course of evolution. A comparative amino acid analysis using ClustalW revealed similarities and differences between buffalo and cattle CLCN proteins. Three duplicated gene pairs were identified, all of which were segmental duplications except for CLCN4-CLCN5, which was a tandem duplication in buffalo. For each gene pair, the Ka/Ks test ratio findings showed that none of the ratios was more than one, indicating that these proteins were likely subject to positive selection. A synteny analysis confirmed a conserved pattern of genomic blocks between buffalo and cattle. Transcriptional control in cells relies on the binding of transcription factors to specific sites in the genome. The number of transcription factor binding sites (TFBSs) was higher in cattle compared to buffalo. Five main recombination breakpoints were identified at various places in the recombination analysis. The outcomes of our study provide new knowledge about the CLCN gene family in buffalo and open the door for further research on candidate genes in vertebrates through genome-wide studies.

Molecular characterization of cassava zinc finger-homeodomain (ZF-HD) transcription factors reveals their role in disease resistance

The production of cassava (Manihot esculenta Crantz) is constantly threatened by cassava bacterial blight (CBB), caused by Xanthomonas phaseoli pv. manihotis (Xpm). Zinc finger-homeodomain (ZF-HD) belongs to a family of homozygous heterotypic cassette genes widely implicated in various developmental and physiological processes in plants. Despite their importance, a comprehensive analysis of ZF-HD genes, particularly those involved in disease resistance, has not been performed for cassava. In the present study, we utilized bioinformatics methods to identify 21 ZF-HD genes distributed across 11 chromosomes of cassava genome, with the majority exhibiting gene structure without introns. Phylogenetic analysis categorized these genes into two major groups (MIF and ZHD) with five subgroups. We observed fourteen pairs of duplicated genes, suggesting that segmental duplication has likely facilitated the expansion of the cassava ZF-HD gene family. Comparative orthologous analyses between cassava and other plant species shed light on the evolutionary trajectory of this gene family. Promoter analyses revealed multiple hormone- and stress-related elements, indicative of a functional role in stress responses. Expression profiling through RNA-seq and quantitative real-time PCR (qRT-PCR) demonstrated that certain cassava ZF-HD genes are up-regulated in response to Xpm infection, suggesting their involvement in defense mechanisms. Notably, MeZHD7 gene was identified via virus induced gene silencing (VIGS) as potentially crucial in conferring resistance against CBB. Results from subcellular localization experiments indicated that MeZHD7 was localized in the nucleus. The Luciferase reporter assay demonstrated an interaction between MeZHD7 and MeMIF5. These findings may lay the foundation for further cloning and functional analyses of cassava ZF-HD genes, particularly those associated with pathogen resistance.

Altered interactions between cis-regulatory elements partially resolve BLADE-ON-PETIOLE genetic redundancy in Capsella rubella.

Duplicated genes are thought to follow one of three evolutionary trajectories that resolve their redundancy: neofunctionalization, subfunctionalization, or pseudogenization. Differences in expression patterns have been documented for many duplicated gene pairs and interpreted as evidence of subfunctionalization and a loss of redundancy. However, little is known about the functional impact of such differences and about their molecular basis. Here, we investigate the genetic and molecular basis for the partial loss of redundancy between the two BLADE-ON-PETIOLE genes BOP1 and BOP2 in red shepherd's purse (Capsella rubella) compared to Arabidopsis (Arabidopsis thaliana). While both genes remain almost fully redundant in A. thaliana, BOP1 in C. rubella can no longer ensure wild-type floral organ numbers and suppress bract formation, due to an altered expression pattern in the region of the cryptic bract primordium. We use two complementary approaches, transgenic rescue of A. thaliana atbop1 atbop2 double mutants and deletions in the endogenous AtBOP1 promoter, to demonstrate that several BOP1 promoter regions containing conserved noncoding sequences interact in a nonadditive manner to control BOP1 expression in the bract primordium and that changes in these interactions underlie the evolutionary divergence between C. rubella and A. thaliana BOP1 expression and activity. Similarly, altered interactions between cis-regulatory regions underlie the divergence in functional promoter architecture related to the control of floral organ abscission by BOP1. These findings highlight the complexity of promoter architecture in plants and suggest that changes in the interactions between cis-regulatory elements are key drivers for evolutionary divergence in gene expression and the loss of redundancy.

Genome-wide analysis of the HD-Zip IV gene family in Camellia sinensis and functional characterization of CsHDZIV3 under trichome development

The HD-Zip IV transcription factors, which belong to the homeodomain leucine zipper (HD-ZIP) family, are essential for the development of epidermal cells in plants. However, there is limited information about HD-Zip IV in tea plant (Camellia sinensis L.). Here we identified 10 HD-Zip IV members (CsHDZIV1∼10) in tea plant, which are distributed in 8 chromosomes with two segmentally duplicated gene pairs. The exon number of CsHDZIVs ranges from 9 to 11, with most genes having two conserved motifs in the 3’UTR. CsHDZIV1∼9 have almost identical motifs and tertiary structures. The promoter region of CsHDZIVs was predicted to harbor light, phytohormone, stress, growth and developmental elements. 9 CsHDZIVs are divided into 5 subfamilies and located in the nucleus. They exhibited expression trends across 12 tea plant tissues and responded to PEG, cold, NaCl and shading treatments. CsHDZIV3 overexpression partially and completely restored trichome development defects of gl2 leaves and stems, respectively. CsHDZIV3 showed histochemical staining in Arabidopsis trichomes and protein interaction with itself. It might regulate trichome development through forming homodimers. These findings contribute to functional study of tea plant HD-Zip IV and provide guiding significance for the investigation of trichome development and the improvement of trichome traits.

Identifying potential anthocyanin biosynthesis regulator in Chinese cherry by comprehensive genome-wide characterization of the R2R3-MYB transcription factor gene family

BackgroundChinese cherry [Cerasus pseudocerasus (Lindl.) G.Don] (syn. Prunus pseudocerasus Lindl.) is an economically important fruiting cherry species with a diverse range of attractive colors, spanning from the lightest yellow to the darkest black purple. However, the MYB transcription factors involved in anthocyanin biosynthesis underlying fruit color variation in Chinese cherry remain unknown.ResultsIn this study, we characterized the R2R3-MYB gene family of Chinese cherry by genome-wide identification and compared it with those of 10 Rosaceae relatives and Arabidopsis thaliana. A total of 1490 R2R3-MYBs were classified into 43 subfamilies, which included 29 subfamilies containing both Rosaceae MYBs and AtMYBs. One subfamily (S45) contained only Rosaceae MYBs, while three subfamilies (S12, S75, and S77) contained only AtMYBs. The variation in gene numbers within identical subfamilies among different species and the absence of certain subfamilies in some species indicated the species-specific expansion within MYB gene family in Chinese cherry and its relatives. Segmental and tandem duplication events primarily contributed to the expansion of Chinese cherry R2R3-CpMYBs. The duplicated gene pairs underwent purifying selection during evolution after duplication events. Phylogenetic relationships and transcript profiling revealed that CpMYB10 and CpMYB4 are involved in the regulation of anthocyanin biosynthesis in Chinese cherry fruits. Expression patterns, transient overexpression and VIGS results confirmed that CpMYB10 promotes anthocyanin accumulation in the fruit skin, while CpMYB4 acts as a repressor, inhibiting anthocyanin biosynthesis of Chinese cherry.ConclusionsThis study provides a comprehensive and systematic analysis of R2R3-MYB gene family in Chinese cherry and Rosaceae relatives, and identifies two regulators, CpMYB10 and CpMYB4, involved in anthocyanin biosynthesis in Chinese cherry. These results help to develop and utilize the potential functions of anthocyanins in Chinese cherry.