Genome-wide analysis of cotton C2H2-zinc finger transcription factor family and their expression analysis during fiber development

Overexpression of wheat C2H2 zinc finger protein transcription factor TaZAT8-5B enhances drought tolerance and root growth in Arabidopsis thaliana

MAIN CONCLUSION

TaZAT8-5B, a C2H2 zinc finger protein transcription factor, positively regulates drought tolerance in transgenic Arabidopsis. It promotes root growth under drought stress via the Aux/IAA-ARF module in the auxin signaling pathway. C2H2 zinc finger proteins (C2H2-ZFPs) represent the largest but relatively unexplored family of transcription factors in plants. This is particularly evident in wheat, where the functions of only a few C2H2-ZFP genes have been confirmed. In this study, we identified a novel C2H2-ZFP gene, TaZAT8-5B. This gene shows high expression in roots and flowers and is significantly induced by heat, drought, and salt stress. Under drought stress, overexpressing TaZAT8-5B in Arabidopsis resulted in increased proline content and superoxide dismutase (SOD) activity in leaves. It also led to reduced stomatal aperture and water loss, while inducing the expression of P5CS1, RD29A, and DREB1A. Consequently, it alleviated drought stress-induced malondialdehyde (MDA) accumulation and improved drought tolerance. Additionally, TaZAT8-5B promoted lateral root initiation under mannitol stress and enhanced both lateral and primary root growth under long-term drought stress. Moreover, TaZAT8-5B was induced by indole-3-acetic acid (IAA). Overexpressing TaZAT8-5B under drought stress significantly inhibited the expression of auxin signaling negative regulatory genes IAA12 and IAA14. Conversely, downstream genes (ARF7, LBD16, LBD18, and CDKA1) of IAA14 and IAA12 were upregulated in TaZAT8-5B overexpressing plants compared to wild-type (WT) plants. These findings suggest that TaZAT8-5B regulates root growth and development under drought stress via the Aux/IAA-ARF module in the auxin signaling pathway. In summary, this study elucidates the role of TaZAT8-5B in enhancing drought tolerance and its involvement in root growth and development through the auxin signaling pathway. These findings offer new insights into the functional analysis of homologous genes of TaZAT8-5B, particularly in Gramineae species.

The 4Fs of cotton: genome editing of cotton for fiber, food, feed, and fuel to achieve zero hunger

Cotton is globally known for its high-priority cellulose-rich natural fiber. In addition to providing fiber for the textile industry, it is an important source material for edible oil, livestock feed, and fuel products. Global warming and the growing population are the major challenges to the world's agriculture and the potential risks to food security. In this context, improving output traits in cotton is necessary to achieve sustainable cotton production. During the last few years, high throughput omics techniques have aided in identifying crucial genes associated with traits of cotton fiber, seed, and plant architecture which could be targeted with more precision and efficiency through the CIRPSR/Cas-mediated genome editing technique. The various CRISPR/Cas systems such as CRISPR/Cas9, CRISPR/nCas9, and CRISPR/Cas12a have been employed to edit cotton genes associated with a wide range of traits including fiber length, flowering, leaf colour, rooting, seed oil, plant architecture, gossypol content, somatic embryogenesis, and biotic and abiotic stresses tolerance, highlighting its effectiveness in editing the cotton genome. Thus, CRISPR/Cas-mediated genome editing has emerged as a technique of choice to tailor crop phenotypes for better yield potential and environmental resilience. The review covers a comprehensive analysis of cotton phenotypic traits and their improvement with the help of the latest genome editing tools to improve fiber, food, feed, and fuel-associated genes of cotton to ensure food security.

Genome-Wide Identification and Expression Analysis of the Cys2His2 Zinc Finger Protein Gene Family in Flammulina filiformis

Zinc finger proteins (ZFPs) are essential transcription factors in eukaryotes, particularly the extensively studied C2H2 family, which is known for its involvement in various biological processes. This research provides a thorough examination and analysis of the C2H2-ZFP gene family in Flammulina filiformis. Using bioinformatics tools, 58 FfC2H2-ZFP genes spread across 11 chromosomes were identified and scrutinized in detail for their gene structures, protein characteristics, and phylogenetic relationships. The study of phylogenetics and synteny sheds light on the evolutionary relationships among C2H2-ZFPs in F. filiformis and other fungi, revealing a complex evolutionary past. The identification of conserved cis-regulatory elements in the gene promoter regions suggests intricate functionalities, particularly in the developmental and stress response pathways. By utilizing RNA-seq and qRT-PCR techniques, the expression patterns of these genes were explored across different developmental stages and tissues of F. filiformis, unveiling distinct expression profiles. Notably, significant expression variations were observed in the stipe elongation region and pilei of various sizes, indicating potential roles in fruiting body morphogenesis. This study enhances our knowledge of the C2H2-ZFP gene family in F. filiformis and lays the groundwork for future investigations into their regulatory mechanisms and applications in fungal biology and biotechnology.

Genome-wide identification of the Q-type C2H2 zinc finger protein gene family and expression analysis under abiotic stress in lotus (Nelumbo nucifera G.)

BACKGROUND

Lotus (Nelumbo nucifera G.) is an important aquatic plant with high ornamental, economic, cultural and ecological values, but abiotic stresses seriously affect its growth and distribution. Q-type C2H2 zinc finger proteins (ZFPs) play an important role in plant growth development and environmental stress responses. Although the Q-type C2H2 gene family has been identified in some plants, limited reports has been carried out it in lotus.

RESULTS

In this study, we identified 45 Q-type NnZFP members in lotus. Based on the phylogenetic tree, these Q-type NnZFP gene family members were divided into 4 groups, including C1-1i, C1-2i, C1-3i and C1-4i. Promoter cis-acting elements analysis indicated that most Q-type NnZFP gene family members in lotus were associated with response to abiotic stresses. Through collinearity analyses, no tandem duplication gene pairs and 14 segmental duplication gene pairs were identified, which showed that duplication events might play a key role in the expansion of the Q-type NnZFP gene family. The synteny results suggested that 54 and 28 Q-type NnZFP genes were orthologous to Arabidopsis and rice, respectively. The expression patterns of these Q-type NnZFP genes revealed that 30 Q-type NnZFP genes were expressed in at least one lotus tissue. Nn5g30550 showed relatively higher expression levels in all tested tissues. 12 genes were randomly selected with at least one gene from each phylogenetic clade, and the expression of these selected genes were confirmed by qRT-PCR (quantitative real-time polymerase chain reaction). The results indicated that Q-type NnZFP genes were extensively involved in cadmium, drought, salt and cold stresses responses. Among them, 11 genes responded to at least three different stress treatments, especially Nn2g12894, which induced by all four treatments.

CONCLUSIONS

These results could increase our understanding of the characterization of the Q-type NnZFP gene family and provide relevant information for further functional analysis of Q-type NnZFP genes in plant development, and abiotic stress tolerance in lotus.

Identification and expression of the Di19 gene family in response to abiotic stress in common bean (Phaseolus vulgaris L.)

Drought-induced 19 (Di19) protein plays critical biological functions in response to adversity as well as in plant growth and development. Exploring the role and mechanism of Di19 in abiotic stress responses is of great significance for improving plant tolerance. In this study, six Di19 genes were identified in the common bean (Phaseolus vulgaris L.), which were mainly derived from segmental duplications. These genes share conserved exon/intron structures and were classified into three subfamilies based on their phylogenetic relationships. The composition and arrangement of conserved motifs were consistent with their phylogenetic relationships. Many hormone- and stress-responsive elements were distributed in the promoters region of PvDi19 genes. Variations in histidine residues in the Cys2/His2 (C2H2) zinc-finger domains resulted in an atypical tertiary structure of PvDi19-5. Gene expression analysis showed rapid induction of PvDi19-1 in roots by 10% PEG treatment, and PvDi19-2 in leaves by 20% PEG treatment, respectively. Most PvDi19s exhibited insensitivity to saline-alkali stress, except for PvDi19-6, which was notably induced during later stages of treatment. The most common bean Di19 genes were inhibited or not regulated by cadmium stress, but the expression of PvDi19-6 in roots was significantly upregulated when subjected to lower concentrations of cadmium (5 mmol). Moreover, Di19s exhibited greater sensitivity to severe cold stress (6°C). These findings enhance our understanding of the role of PvDi19s in common bean abiotic stress responses and provide a basis for future genetic enhancements in common bean stress tolerance.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Genome-wide analysis of the C2H2-ZFP gene family in Stevia rebaudiana reveals involvement in abiotic stress response

Stevia (Stevia rebaudiana Bertoni) is a natural sweetener plant that accumulates highly sweet steviol glycosides (SGs) especially in leaves. Stevia is native to humid areas and does not have a high tolerance to drought which is the most serious abiotic stress restricting its production worldwide. C2H2 zinc finger proteins (C2H2-ZFPs) are a group of well-known transcription factors that involves in various developmental, physiological and biochemical activities as well as in response to abiotic stresses. Here we analyzed C2H2-ZFP gene family in stevia and identified a total of 185 putative SrC2H2-ZF proteins from the genome sequence of S. rebaudiana. We further characterized the identified C2H2-ZF domains and their organization, additional domains and motifs and analyzed their physicochemical properties, localization and gene expression patterns. The cis-element analysis suggested multiple roles of SrC2H2-ZFPs in response to light, phytohormone, and abiotic stresses. In silico analysis revealed that the stevia C2H2-ZFP genes are interactively expressed in different tissues and developmental stages and some C2H2-ZFP genes are involved in response to drought stress. This study provides a background for future exploration of the functional, and regulatory aspects of the C2H2-ZFP gene family in S. rebaudiana.

Genome-wide identification and expression analysis of the C2H2-zinc finger transcription factor gene family and screening of candidate genes involved in floral development in Coptis teeta Wall. (Ranunculaceae)

Background: C2H2-zinc finger transcription factors comprise one of the largest and most diverse gene superfamilies and are involved in the transcriptional regulation of flowering. Although a large number of C2H2 zinc-finger proteins (C2H2-ZFPs) have been well characterized in a number of model plant species, little is known about their expression and function in Coptis teeta. C. teeta displays two floral phenotypes (herkogamy phenotypes). It has been proposed that the C2H2-zinc finger transcription factor family may play a crucial role in the formation of floral development and herkogamy observed in C. teeta. As such, we performed a genome-wide analysis of the C2H2-ZFP gene family in C. teeta. Results: The complexity and diversity of C. teeta C2H2 zinc finger proteins were established by evaluation of their physicochemical properties, phylogenetic relationships, exon-intron structure, and conserved motifs. Chromosome localization showed that 95 members of the C2H2 zinc-finger genes were unevenly distributed across the nine chromosomes of C. teeta, and that these genes were replicated in tandem and segmentally and had undergone purifying selection. Analysis of cis-acting regulatory elements revealed a possible involvement of C2H2 zinc-finger proteins in the regulation of phytohormones. Transcriptome data was then used to compare the expression levels of these genes during the growth and development of the two floral phenotypes (F-type and M-type). These data demonstrate that in groups A and B, the expression levels of 23 genes were higher in F-type flowers, while 15 genes showed higher expressions in M-type flowers. qRT-PCR analysis further revealed that the relative expression was highly consistent with the transcriptome data. Conclusion: These data provide a solid basis for further in-depth studies of the C2H2 zinc finger transcription factor gene family in this species and provide preliminary information on which to base further research into the role of the C2H2 ZFPs gene family in floral development in C. teeta.

Genome-wide identification of the C2H2 zinc finger gene family and expression analysis under salt stress in sweetpotato

Introduction

The higher plant transcription factor C2H2 zinc finger protein (C2H2-ZFP) is essential for plant growth, development, and stress response. There are limited studies on C2H2-ZFP genes in sweetpotato, despite a substantial number of C2H2-ZFP genes having been systematically found in plants.

Methods

In this work, 178 C2H2-ZFP genes were found in sweetpotato, distributed randomly on 15 chromosomes, and given new names according to where they were located. These members of the zinc finger gene family are separated into six branches, as shown by the phylogenetic tree. 24 tandem repeats of IbZFP genes and 46 fragment repeats were identified, and a homology study revealed that IbZFP genes linked more regions with wild relative species of sweetpotato as well as rhizome plants like potato and cassava. And we analyzed the expression patterns of IbZFP genes during the early development of sweetpotato storage roots (SRs) and salt stress using transcriptome data, and identified 44 IbZFP genes that exhibited differences in expression levels during the early expansion of sweetpotato SRs in different varieties, and 92 IbZFP genes that exhibited differences in expression levels under salt stress in salt tolerant and salt sensitive sweetpotato varieties. Additionally, we cloned six IbZFP genes in sweetpotato and analyzed their expression patterns in different tissues, their expression patterns under abiotic stress and hormone treatment, and subcellular localization.

Results and discussion

The results showed that the IbZFP genes had tissue specificity in sweetpotato and were induced to varying degrees by drought and salt stress. ABA and GA3 treatments also affected the expression of the IbZFP genes. We selected IbZFP105, which showed significant differences in expression levels under salt stress and ABA treatment, to be heterologously expressed in Arabidopsis thaliana. We found that IbZFP105 OE lines exhibited higher tolerance to salt stress and ABA stress. This indicates that IbZFP105 can enhance the salt tolerance of plants. These results systematically identified the evolution and expression patterns of members of the C2H2-ZFP gene family in sweetpotato, providing a theoretical basis for studying the role of IbZFP genes in the development of sweetpotato SRs and in resistance to stress.

Genome-wide identification of the C2H2-Zinc finger gene family and functional validation of CsZFP7 in citrus nucellar embryogenesis

KEY MESSAGE

Genome-wide identification of C2H2-ZF gene family in the poly- and mono-embryonic citrus species and validation of the positive role of CsZFP7 in sporophytic apomixis. The C2H2 zinc finger (C2H2-ZF) gene family is involved in plant vegetative and reproductive development. Although a large number of C2H2 zinc-finger proteins (C2H2-ZFPs) have been well characterized in some horticultural plants, little is known about the C2H2-ZFPs and their function in citrus. In this work, we performed a genome-wide sequence analysis and identified 97 and 101 putative C2H2-ZF gene family members in the genomes of sweet orange (C. sinensis, poly-embryonic) and pummelo (C. grandis, mono-embryonic), respectively. Phylogenetic analysis categorized citrus C2H2-ZF gene family into four clades, and their possible functions were inferred. According to the numerous regulatory elements on promoter, citrus C2H2-ZFPs can be divided into five different regulatory function types that indicate functional differentiation. RNA-seq data revealed 20 differentially expressed C2H2-ZF genes between poly-embryonic and mono-embryonic ovules at two stages of citrus nucellar embryogenesis, among them CsZFP52 specifically expressed in mono-embryonic pummelo ovules, while CsZFP7, 37, 44, 45, 67 and 68 specifically expressed in poly-embryonic sweet orange ovules. RT-qPCR further validated that CsZFP7 specifically expressed at higher levels in poly-embryonic ovules, and down-regulation of CsZFP7 in the poly-embryonic mini citrus (Fortunella hindsii) increased rate of mono-embryonic seeds compared with the wild type, indicating the regulatory potential of CsZFP7 in nucellar embryogenesis of citrus. This work provided a comprehensive analysis of C2H2-ZF gene family in citrus, including genome organization and gene structure, phylogenetic relationships, gene duplications, possible cis-elements on promoter regions and expression profiles, especially in the poly- and mono-embryogenic ovules, and suggested that CsZFP7 is involved in nucellar embryogenesis.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Genome-wide identification, comprehensive characterization of transcription factors, cis-regulatory elements, protein homology, and protein interaction network of DREB gene family in Solanum lycopersicum

Tomato is a drought-sensitive crop which has high susceptibility to adverse climatic changes. Dehydration-responsive element-binding (DREB) are significant plant transcription factors that have a vital role in regulating plant abiotic stress tolerance by networking with DRE/CRT cis-regulatory elements in response to stresses. In this study, bioinformatics analysis was performed to conduct the genome-wide identification and characterization of DREB genes and promoter elements in Solanum lycopersicum. In genome-wide coverage, 58 SlDREB genes were discovered on 12 chromosomes that justified the criteria of the presence of AP2 domain as conserved motifs. Intron-exon organization and motif analysis showed consistency with phylogenetic analysis and confirmed the absence of the A3 class, thus dividing the SlDREB genes into five categories. Gene expansion was observed through tandem duplication and segmental duplication gene events in SlDREB genes. Ka/Ks values were calculated in ortholog pairs that indicated divergence time and occurrence of purification selection during the evolutionary period. Synteny analysis demonstrated that 32 out of 58 and 47 out of 58 SlDREB genes were orthologs to Arabidopsis and Solanum tuberosum, respectively. Subcellular localization predicted that SlDREB genes were present in the nucleus and performed primary functions in DNA binding to regulate the transcriptional processes according to gene ontology. Cis-acting regulatory element analysis revealed the presence of 103 motifs in 2.5-kbp upstream promoter sequences of 58 SlDREB genes. Five representative SlDREB proteins were selected from the resultant DREB subgroups for 3D protein modeling through the Phyre2 server. All models confirmed about 90% residues in the favorable region through Ramachandran plot analysis. Moreover, active catalytic sites and occurrence in disorder regions indicated the structural and functional flexibility of SlDREB proteins. Protein association networks through STRING software suggested the potential interactors that belong to different gene families and are involved in regulating similar functional and biological processes. Transcriptome data analysis has revealed that the SlDREB gene family is engaged in defense response against drought and heat stress conditions in tomato. Overall, this comprehensive research reveals the identification and characterization of SlDREB genes that provide potential knowledge for improving abiotic stress tolerance in tomato.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton

Interspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056-1255 ENAs showed transgressive expression in fibres at 5-15 dpa (days post-anthesis) of some BILs, 520 of which were located in FL-quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL-QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR-cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Genome-Wide Analysis of C2H2 Zinc Finger Gene Family and Its Response to Cold and Drought Stress in Sorghum [Sorghum bicolor (L.) Moench]

C2H2 zinc finger protein (C2H2-ZFP) is one of the most important transcription factor families in higher plants. In this study, a total of 145 C2H2-ZFPs was identified in Sorghum bicolor and randomly distributed on 10 chromosomes. Based on the phylogenetic tree, these zinc finger gene family members were divided into 11 clades, and the gene structure and motif composition of SbC2H2-ZFPs in the same clade were similar. SbC2H2-ZFP members located in the same clade contained similar intron/exon and motif patterns. Thirty-three tandem duplicated SbC2H2-ZFPs and 24 pairs of segmental duplicated genes were identified. Moreover, synteny analysis showed that sorghum had more collinear regions with monocotyledonous plants such as maize and rice than did dicotyledons such as soybean and Arabidopsis. Furthermore, we used quantitative RT-PCR (qRT-PCR) to analyze the expression of C2H2-ZFPs in different organs and demonstrated that the genes responded to cold and drought. For example, Sobic.008G088842 might be activated by cold but is inhibited in drought in the stems and leaves. This work not only revealed an important expanded C2H2-ZFP gene family in Sorghum bicolor but also provides a research basis for determining the role of C2H2-ZFPs in sorghum development and abiotic stress resistance.

Genome-wide identification of C2H2-type zinc finger gene family members and their expression during abiotic stress responses in orchardgrass (Dactylis glomerata)

The C2H2-type zinc finger protein (ZFP) family is one of the largest transcription factor families in the plant kingdom and its members are involved in plant growth, development, and stress responses. As an economically valuable perennial graminaceous forage crop, orchardgrass (Dactylis glomerata) is an important feedstuff resource owing to its high yield and quality. In this study, 125 C2H2-type ZFPs in orchardgrass (Dg-ZFPs) were identified and further classified by phylogenetic analysis. The members with similar gene structures were generally clustered into the same groups, with proteins containing the conserved QALGGH motif being concentrated in groups VIII and IX. Gene ontology and miRNA target analyses indicated that Dg-ZFPs likely perform diverse biological functions through their gene interactions. The RNA-seq data revealed differentially expressed genes across tissues and development phases, suggesting that some Dg-ZFPs might participate in growth and development regulation. Abiotic stress responses of Dg-ZFP genes were verified by qPCR and Saccharomyces cerevisiae transformation, revealing that Dg-ZFP125 could enhance the tolerance of yeasts to osmotic and salt stresses. Our study performed a novel systematic analysis of Dg-ZFPs in orchardgrass, providing a reference for this gene family in other grasses and revealing new insights for enhancing gene utilization.

Molecular structures and functional exploration of NDA family genes respond tolerant to alkaline stress in Gossypium hirsutum L

Background

The internal NAD(P)H dehydrogenase (NDA) gene family was a member of the NAD(P)H dehydrogenase (ND) gene family, mainly involved in the non-phosphorylated respiratory pathways in mitochondria and played crucial roles in response to abiotic stress.

Methods

The whole genome identification, structure analysis and expression pattern of NDA gene family were conducted to analyze the NDA gene family.

Results

There were 51, 52, 26, and 24 NDA genes identified in G. hirsutum, G. barbadense, G. arboreum and G. raimondii, respectively. According to the structural characteristics of genes and traits of phylogenetic tree, we divided the NDA gene family into 8 clades. Gene structure analysis showed that the NDA gene family was relatively conservative. The four Gossypium species had good collinearity, and segmental duplication played an important role in the evolution of the NDA gene family. Analysis of cis-elements showed that most GhNDA genes contained cis-elements related to light response and plant hormones (ABA, MeJA and GA). The analysis of the expression patterns of GhNDA genes under different alkaline stress showed that GhNDA genes were actively involved in the response to alkaline stress, possibly through different molecular mechanisms. By analyzing the existing RNA-Seq data after alkaline stress, it was found that an NDA family gene GhNDA32 was expressed, and then theGhNDA32 was silenced by virus-induced gene silencing (VIGS). By observing the phenotype, we found that the wilting degree of silenced plants was much higher than that of the control plant after alkaline treatment, suggesting that GhNDA32 gene was involved in the response to alkaline stress.

Conclusions

In this study, GhNDAs participated in response to alkaline stress, especially NaHCO₃ stress. It was of great significance for the future research on the molecular mechanism of NDA gene family in responding to abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00372-8.

Genome-wide analysis of ZAT gene family revealed GhZAT6 regulates salt stress tolerance in G. hirsutum

High salt environments can induce stress in different plants. The genes containing the ZAT domain constitute a family that belongs to a branch of the C₂H₂ family, which plays a vital role in responding to abiotic stresses. In this study, we identified 169 ZAT genes from seven plant species, including 44 ZAT genes from G. hirsutum. Phylogenetic tree analysis divided ZAT genes in six groups with conserved gene structure, protein motifs. Two C₂H₂ domains and an EAR domain and even chromosomal distribution on At and Dt sub-genome chromosomes of G. hirsutum was observed. GhZAT6 was primarily expressed in the root tissue and responded to NaCl and ABA treatments. Subcellular localization found that GhZAT6 was located in the nucleus and demonstrated transactivation activity during a transactivation activity assay. Arabidopsis transgenic lines overexpressing the GhZAT6 gene showed salt tolerance and grew more vigorously than WT on MS medium supplemented with 100 mmol NaCl. Additionally, the silencing of the GhZAT6 gene in cotton plants showed more obvious leaf wilting than the control plants, which were subjected to 400 mmol NaCl treatment. Next, the expressions of GhAPX1, GhFSD1, GhFSD2, and GhSOS3 were significantly lower in the GhZAT6-silenced plants treated with NaCl than the control. Based on these findings, GhZAT6 may be involved in the ABA pathway and mediate salt stress tolerance by regulating ROS-related gene expression.

Advances in the Regulation of Epidermal Cell Development by C2H2 Zinc Finger Proteins in Plants

Plant epidermal cells, such as trichomes, root hairs, salt glands, and stomata, play pivotal roles in the growth, development, and environmental adaptation of terrestrial plants. Cell fate determination, differentiation, and the formation of epidermal structures represent basic developmental processes in multicellular organisms. Increasing evidence indicates that C2H2 zinc finger proteins play important roles in regulating the development of epidermal structures in plants and plant adaptation to unfavorable environments. Here, we systematically summarize the molecular mechanism underlying the roles of C2H2 zinc finger proteins in controlling epidermal cell formation in plants, with an emphasis on trichomes, root hairs, and salt glands and their roles in plant adaptation to environmental stress. In addition, we discuss the possible roles of homologous C2H2 zinc finger proteins in trichome development in non-halophytes and salt gland development in halophytes based on bioinformatic analysis. This review provides a foundation for further study of epidermal cell development and abiotic stress responses in plants.

Genome-wide analysis, identification, evolution and genomic organization of dehydration responsive element-binding (DREB) gene family in Solanum tuberosum

Background

The dehydration responsive element-binding (DREB) gene family plays a crucial role as transcription regulators and enhances plant tolerance to abiotic stresses. Although the DREB gene family has been identified and characterized in many plants, knowledge about it in Solanum tuberosum (Potato) is limited.

Results

In the present study, StDREB gene family was comprehensively analyzed using bioinformatics approaches. We identified 66 StDREB genes through genome wide screening of the Potato genome based on the AP2 domain architecture and amino acid conservation analysis (Valine at position 14th). Phylogenetic analysis divided them into six distinct subgroups (A1–A6). The categorization of StDREB genes into six subgroups was further supported by gene structure and conserved motif analysis. Potato DREB genes were found to be distributed unevenly across 12 chromosomes. Gene duplication proved that StDREB genes experienced tandem and segmental duplication events which led to the expansion of the gene family. The Ka/Ks ratios of the orthologous pairs also demonstrated the StDREB genes were under strong purification selection in the course of evolution. Interspecies synteny analysis revealed 45 and 36 StDREB genes were orthologous to Arabidopsis and Solanum lycopersicum, respectively. Moreover, subcellular localization indicated that StDREB genes were predominantly located within the nucleus and the StDREB family’s major function was DNA binding according to gene ontology (GO) annotation.

Conclusions

This study provides a comprehensive and systematic understanding of precise molecular mechanism and functional characterization of StDREB genes in abiotic stress responses and will lead to improvement in Solanum tuberosum.

Genomic insight into the divergence and adaptive potential of a forgotten landrace G. hirsutum L. purpurascens

Wild progenitors are an excellent source for strengthening the genetic basis and accumulation of desirable variation lost because of directional selection and adaptation in modern cultivars. Here, we re-evaluate a landrace of Gossypium hirsutum, formerly known as Gossypium purpurascens. Our study seeks to understand the genomic structure, variation, and breeding potential of this landrace, providing potential insights into the biogeographic history and genomic changes likely associated with domestication. A core set of accessions, including current varieties, obsolete accessions, G. purpurascens, and other geographical landraces, are subjected to genotyping along with multilocation phenotyping. Population fixation statistics suggests a marked differentiation between G. purpurascens and three other groups, emphasizing the divergent genomic behavior of G. purpurascens. Phylogenetic analysis establishes the primitive nature of G. purpurascens, identifying it as a vital source of functional variation, the inclusion of which in the upland cotton (cultivated G. hirsutum) gene pool may broaden the genetic basis of modern cultivars. Genome-wide association results indicate multiple loci associated with domestication regions corresponding to flowering and fiber quality. Moreover, the conserved nature of G. purpurascens can also provide insights into the evolutionary process of G. hirsutum.

Genome-wide identification of the DUF668 gene family in cotton and expression profiling analysis of GhDUF668 in Gossypium hirsutum under adverse stress

Background

Domain of unknown function 668 (DUF668) may play a crucial role in the plant growth and developmental response to adverse stress. However, our knowledge of the function of the DUF668 gene family is limited.

Results

Our study was conducted based on the DUF668 gene family identified from cotton genome sequencing. Phylogenetic analysis showed that the DUF668 family genes can be classified into four subgroups in cotton. We identified 32 DUF668 genes, which are distributed on 17 chromosomes and most of them located in the nucleus of Gossypium hirsutum. Gene structure and motif analyses revealed that the members of the DUF668 gene family can be clustered in G. hirsutum into two broad groups, which are relatively evolutionarily conserved. Transcriptome data analysis showed that the GhDUF668 genes are differentially expressed in different tissues under various stresses (cold, heat, drought, salt, and Verticillium dahliae), and expression is generally increased in roots and stems. Promoter and expression analyses indicated that Gh_DUF668–05, Gh_DUF668–08, Gh_DUF668–11, Gh_DUF668–23 and Gh_DUF668–28 in G. hirsutum might have evolved resistance to adverse stress. Additionally, qRT-PCR revealed that these 5 genes in four cotton lines, KK1543 (drought resistant), Xinluzao 26 (drought sensitive), Zhongzhimian 2 (disease resistant) and Simian 3 (susceptible), under drought and Verticillium wilt stress were all significantly induced. Roots had the highest expression of these 5 genes before and after the treatment. Among them, the expression levels of Gh_DUF668–08 and Gh_DUF668–23 increased sharply at 6 h and reached a maximum at 12 h under biotic and abiotic stress, which showed that they might be involved in the process of adverse stress resistance in cotton.

Conclusion

The significant changes in GhDUF668 expression in the roots after adverse stress indicate that GhDUF668 is likely to increase plant resistance to stress. This study provides an important theoretical basis for further research on the function of the DUF668 gene family and the molecular mechanism of adverse stress resistance in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07716-w.

Genome-Wide Identification and Expression Patterns of the C2H2-Zinc Finger Gene Family Related to Stress Responses and Catechins Accumulation in Camellia sinensis [L.] O. Kuntze

The C2H2-zinc finger protein (C2H2-ZFP) is essential for the regulation of plant development and widely responsive to diverse stresses including drought, cold and salt stress, further affecting the late flavonoid accumulation in higher plants. Tea is known as a popular beverage worldwide and its quality is greatly dependent on the physiological status and growing environment of the tea plant. To date, the understanding of C2H2-ZFP gene family in Camellia sinensis [L.] O. Kuntze is not yet available. In the present study, 134 CsC2H2-ZFP genes were identified and randomly distributed on 15 chromosomes. The CsC2H2-ZFP gene family was classified into four clades and gene structures and motif compositions of CsC2H2-ZFPs were similar within the same clade. Segmental duplication and negative selection were the main forces driving the expansion of the CsC2H2-ZFP gene family. Expression patterns suggested that CsC2H2-ZFPs were responsive to different stresses including drought, salt, cold and methyl jasmonate (MeJA) treatment. Specially, several C2H2-ZFPs showed a significant correlation with the catechins content and responded to the MeJA treatment, which might contribute to the tea quality and specialized astringent taste. This study will lay the foundations for further research of C2H2-type zinc finger proteins on the stress responses and quality-related metabolites accumulation in C. sinensis.

Identification of Chlorophyll Metabolism- and Photosynthesis-Related Genes Regulating Green Flower Color in Chrysanthemum by Integrative Transcriptome and Weighted Correlation Network Analyses

Green chrysanthemums are difficult to breed but have high commercial value. The molecular basis for the green petal color in chrysanthemum is not fully understood. This was investigated in the present study by RNA sequencing analysis of white and green ray florets collected at three stages of flower development from the F₁ progeny of the cross between Chrysanthemum × morifolium “Lüdingdang” with green-petaled flowers and Chrysanthemum vistitum with white-petaled flowers. The chlorophyll content was higher and chloroplast degradation was slower in green pools than in white pools at each developmental stage. Transcriptome analysis revealed that genes that were differentially expressed between the two pools were enriched in pathways related to chlorophyll metabolism and photosynthesis. We identified the transcription factor genes CmCOLa, CmCOLb, CmERF, and CmbHLH as regulators of the green flower color in chrysanthemum by differential expression analysis and weighted gene co-expression network analysis. These findings can guide future efforts to improve the color palette of chrysanthemum flowers through genetic engineering.

Identification of hub genes through co-expression network of major QTLs of fiber length and strength traits in multiple RIL populations of cotton

The present study demonstrated a de novo correlation among fiber quality genes in multiple RIL populations including sGK9708 × 0-153, LMY22 × LY343 and Lumianyan28 × Xinluzao24. The current study was conducted to identify the major common QTLs including fiber length and strength, and to identify the co-expression networks of fiber length and strength QTLs harbored genes to target the hub genes. The RNA-seq data of sGK9708 × 0-153 population highlighted 50 and 48 candidate genes of fiber length and fiber strength QTLs. A total of 29 and 21 hub genes were identified in fiber length and strength co-expression network modules. The absolute values of correlation coefficient close to 1 resulted highly positive correlation among hub genes. Results also suggested that the gene correlation significantly influence the gene expression at different fiber development stages. These results might provide useful reference for further experiments in multiple RIL populations and suggest potential candidate genes for functional studies in cotton.

Genome-Wide Association Analysis Reveals Loci and Candidate Genes Involved in Fiber Quality Traits Under Multiple Field Environments in Cotton (Gossypium hirsutum)

Fiber length, fiber strength, and fiber micronaire are the main fiber quality parameters in cotton. Thus, mining the elite and stable loci/alleles related to fiber quality traits and elucidating the relationship between the two may accelerate genetic improvement of fiber quality in cotton. Here, genome-wide association analysis (GWAS) was performed for fiber quality parameters based on phenotypic data, and 56,010 high-quality single nucleotide polymorphisms (SNPs) using 242 upland cotton accessions under 12 field environments were obtained. Phenotypic analysis exhibited that fiber length (FL) had a positive correlation with fiber strength (FS) and had a negative correlation with fiber micronaire (Mic). Genetic analysis also indicated that FL, FS, and Mic had high heritability of more than 80%. A total of 67 stable quantitative trait loci (QTLs) were identified through GWAS analysis, including 31 for FL, 21 for FS, and 22 for Mic. Of them, three pairs homologous QTLs were detected between A and D subgenomes, and seven co-located QTLs with two fiber quality parameters were found. Compared with the reported QTLs, 34 co-located with previous studies, and 33 were newly revealed. Integrated with transcriptome analysis, we selected 256, 244, and 149 candidate genes for FL, FS, and Mic, respectively. Gene Ontology (GO) analysis showed that most of the genes located in QTLs interval of the three fiber quality traits were involved in sugar biosynthesis, sugar metabolism, microtubule, and cytoskeleton organization, which played crucial roles in fiber development. Through correlation analysis between haplotypes and phenotypes, three genes (GH_A05G1494, GH_D11G3097, and GH_A05G1082) predominately expressed in fiber development stages were indicated to be potentially responsible for FL, FS, and Mic, respectively. The GH_A05G1494 encoded a protein containing SGS-domain, which is related to tubulin-binding and ubiquitin-protein ligase binding. The GH_D11G3097 encoded 20S proteasome beta subunit G1, and was involved in the ubiquitin-dependent protein catabolic process. The GH_A05G1082 encoded RAN binding protein 1 with a molecular function of GTPase activator activity. These results provide new insights and candidate loci/genes for the improvement of fiber quality in cotton.

Bioinformatics of Metalloproteins and Metalloproteomes

Trace metals are inorganic elements that are required for all organisms in very low quantities. They serve as cofactors and activators of metalloproteins involved in a variety of key cellular processes. While substantial effort has been made in experimental characterization of metalloproteins and their functions, the application of bioinformatics in the research of metalloproteins and metalloproteomes is still limited. In the last few years, computational prediction and comparative genomics of metalloprotein genes have arisen, which provide significant insights into their distribution, function, and evolution in nature. This review aims to offer an overview of recent advances in bioinformatic analysis of metalloproteins, mainly focusing on metalloprotein prediction and the use of different metals across the tree of life. We describe current computational approaches for the identification of metalloprotein genes and metal-binding sites/patterns in proteins, and then introduce a set of related databases. Furthermore, we discuss the latest research progress in comparative genomics of several important metals in both prokaryotes and eukaryotes, which demonstrates divergent and dynamic evolutionary patterns of different metalloprotein families and metalloproteomes. Overall, bioinformatic studies of metalloproteins provide a foundation for systematic understanding of trace metal utilization in all three domains of life.