Genome-wide expression analysis suggests glutaredoxin genes response to various stresses in cotton

Genome-Wide Identification of the CYP716 Gene Family in Platycodon grandiflorus (Jacq.) A. DC. and Its Role in the Regulation of Triterpenoid Saponin Biosynthesis

The main type of saponins occurring in the root of Platycodon grandiflorus (Jacq.) A. DC. are oleanolic acid glycosides. The CYP716 gene family plays a major role in catalyzing the conversion of β-amyrin into oleanolic acid. However, studies on the CYP716 genes in P. grandiflorus are limited, and its evolutionary history remains poorly understood. In this study, 22 PgCYP716 genes were identified, distributed among seven subfamilies. Cis-acting elements of the PgCYP716 promoters were mainly involved in plant hormone regulation and responses to abiotic stresses. PgCYP716A264, PgCYP716A391, PgCYP716A291, and PgCYP716BWv3 genes were upregulated in the root and during saponin accumulation, as shown by RNA-seq analysis, suggesting that these four genes play an important role in saponin synthesis. The results of subcellular localization indicated that these four genes encoded membrane proteins. Furthermore, the catalytic activity of these four genes was proved in the yeast, which catalyzed the conversion of β-amyrin into oleanolic acid. We found that the content of β-amyrin, platycodin D, platycoside E, platycodin D3, and total saponins increased significantly when either of the four genes was over expressed in the transgenic hair root. In addition, the expression of PgSS, PgGPPS2, PgHMGS, and PgSE was also upregulated while these four genes were overexpressed. These data support that these four PgCYP716 enzymes oxidize β-amyrin to produce oleanolic acid, ultimately promoting saponin accumulation by activating the expression of upstream pathway genes. Our results enhanced the understanding of the functional variation among the PgCYP716 gene family involved in triterpenoid biosynthesis and provided a theoretical foundation for improving saponin content and enriching the saponin biosynthetic pathway in P. grandiflorus.

Genome-wide identification of walnut (Juglans regia) PME gene family members and expression analysis during infection with Cryptosphaeria pullmanensis pathogens

Walnuts exhibit a higher resistance to diseases, though they are not completely immune. This study focuses on the Pectin methylesterase (PME) gene family to investigate whether it is involved in disease resistance in walnuts. These 21 genes are distributed across 12 chromosomes, with four pairs demonstrating homology. Variations in conserved motifs and gene structures suggest diverse functions within the gene family. Phylogenetic and collinear gene pairs of the PME family indicate that the gene family has evolved in a relatively stable way. The cis-acting elements and gene ontology enrichment of these genes, underscores their potential role in bolstering walnuts' defense mechanisms. Transcriptomic analyses were conducted under conditions of Cryptosphaeria pullmanensis infestation and verified by RT-qPCR. The results showed that certain JrPME family genes were activated in response, leading to the hypothesis that some members may confer resistance to the disease.

GhVIM28, a negative regulator identified from VIM family genes, positively responds to salt stress in cotton

The VIM (belonged to E3 ubiquitin ligase) gene family is crucial for plant growth, development, and stress responses, yet their role in salt stress remains unclear. We analyzed phylogenetic relationships, chromosomal localization, conserved motifs, gene structure, cis-acting elements, and gene expression patterns of the VIM gene family in four cotton varieties. Our findings reveal 29, 29, 17, and 14 members in Gossypium hirsutum (G.hirsutum), Gossypium barbadense (G.barbadense), Gossypium arboreum (G.arboreum), and Gossypium raimondii (G. raimondii), respectively, indicating the maturity and evolution of this gene family. motifs among GhVIMs genes were observed, along with the presence of stress-responsive, hormone-responsive, and growth-related elements in their promoter regions. Gene expression analysis showed varying patterns and tissue specificity of GhVIMs genes under abiotic stress. Silencing GhVIM28 via virus-induced gene silencing revealed its role as a salt-tolerant negative regulator. This work reveals a mechanism by which the VIM gene family in response to salt stress in cotton, identifying a potential negative regulator, GhVIM28, which could be targeted for enhancing salt tolerance in cotton. The objective of this study was to explore the evolutionary relationship of the VIM gene family and its potential function in salt stress tolerance, and provide important genetic resources for salt tolerance breeding of cotton.

Identification and expression profiling of GAPDH family genes involved in response to Sclerotinia sclerotiorum infection and phytohormones in Brassica napus

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is a crucial enzyme in glycolysis, an essential metabolic pathway for carbohydrate metabolism across all living organisms. Recent research indicates that phosphorylating GAPDH exhibits various moonlighting functions, contributing to plant growth and development, autophagy, drought tolerance, salt tolerance, and bacterial/viral diseases resistance. However, in rapeseed (Brassica napus), the role of GAPDHs in plant immune responses to fungal pathogens remains unexplored. In this study, 28 genes encoding GAPDH proteins were revealed in B. napus and classified into three distinct subclasses based on their protein structural and phylogenetic relationships. Whole-genome duplication plays a major role in the evolution of BnaGAPDHs. Synteny analyses revealed orthologous relationships, identifying 23, 26, and 26 BnaGAPDH genes with counterparts in Arabidopsis, Brassica rapa, and Brassica oleracea, respectively. The promoter regions of 12 BnaGAPDHs uncovered a spectrum of responsive elements to biotic and abiotic stresses, indicating their crucial role in plant stress resistance. Transcriptome analysis characterized the expression profiles of different BnaGAPDH genes during Sclerotinia sclerotiorum infection and hormonal treatment. Notably, BnaGAPDH17, BnaGAPDH20, BnaGAPDH21, and BnaGAPDH22 exhibited sensitivity to S. sclerotiorum infection, oxalic acid, hormone signals. Intriguingly, under standard physiological conditions, BnaGAPDH17, BnaGAPDH20, and BnaGAPDH22 are primarily localized in the cytoplasm and plasma membrane, with BnaGAPDH21 also detectable in the nucleus. Furthermore, the nuclear translocation of BnaGAPDH20 was observed under H2O2 treatment and S. sclerotiorum infection. These findings might provide a theoretical foundation for elucidating the functions of phosphorylating GAPDH.

Capsicum chinense Jacq.-derived glutaredoxin (CcGRXS12) alters redox status of the cells to confer resistance against pepper mild mottle virus (PMMoV-I)

KEY MESSAGE

The CcGRXS12 gene protects plants from cellular oxidative damage that are caused by both biotic and abiotic stresses. The protein possesses GSH-disulphide oxidoreductase property but lacks Fe-S cluster assembly mechanism. Glutaredoxins (Grxs) are small, ubiquitous and multi-functional proteins. They are present in different compartments of plant cells. A chloroplast targeted Class I GRX (CcGRXS12) gene was isolated from Capsicum chinense during the pepper mild mottle virus (PMMoV) infection. Functional characterization of the gene was performed in Nicotiana benthamiana transgenic plants transformed with native C. chinense GRX (Nb:GRX), GRX-fused with GFP (Nb:GRX-GFP) and GRX-truncated for chloroplast sequences fused with GFP (Nb:Δ2MGRX-GFP). Overexpression of CcGRXS12 inhibited the PMMoV-I accumulation at the later stage of infection, accompanied with the activation of salicylic acid (SA) pathway pathogenesis-related (PR) transcripts and suppression of JA/ET pathway transcripts. Further, the reduced accumulation of auxin-induced Glutathione-S-Transferase (pCNT103) in CcGRXS12 overexpressing lines indicated that the protein could protect the plants from the oxidative stress caused by the virus. PMMoV-I infection increased the accumulation of pyridine nucleotides (PNs) mainly due to the reduced form of PNs (NAD(P)H), and it was high in Nb:GRX-GFP lines compared to other transgenic lines. Apart from biotic stress, CcGRXS12 protects the plants from abiotic stress conditions caused by H2O2 and herbicide paraquat. CcGRXS12 exhibited GSH-disulphide oxidoreductase activity in vitro; however, it was devoid of complementary Fe-S cluster assembly mechanism found in yeast. Overall, this study proves that CcGRXS12 plays a crucial role during biotic and abiotic stress in plants.

The integration of TRX/GRX systems and phytohormonal signalling pathways in plant stress and development

Plant acclimation to changing environmental conditions involves the interaction of different signalling molecules, including reactive oxygen species and hormones. Redox regulation exerted by thioredoxin (TRX) and glutaredoxin (GRX), two oxidoreductases, is emerging as a specific point of control mediating signal transduction pathways associated with plant growth and stress response. Phytohormones are messengers that coordinate plant cell activities to regulate growth, defence, and productivity, although their cross-talk with components of the redox system is less known. The present review focuses on our current knowledge of the interplay that occurs between TRX and GRX systems and phytohormonal signalling pathways in connection with the control of plant development and stress responses. Here, we consider the regulation that phytohormones exert on TRX and GRX systems, as well as the involvement of these redox proteins in the control of phytohormone-mediated signalling pathways.

Raffinose degradation-related gene GhAGAL3 was screened out responding to salinity stress through expression patterns of GhAGALs family genes

A-galactosidases (AGALs), the oligosaccharide (RFO) catabolic genes of the raffinose family, play crucial roles in plant growth and development and in adversity stress. They can break down the non-reducing terminal galactose residues of glycolipids and sugar chains. In this study, the whole genome of AGALs was analyzed. Bioinformatics analysis was conducted to analyze members of the AGAL family in Gossypium hirsutum, Gossypium arboreum, Gossypium barbadense, and Gossypium raimondii. Meanwhile, RT-qPCR was carried out to analyze the expression patterns of AGAL family members in different tissues of terrestrial cotton. It was found that a series of environmental factors stimulated the expression of the GhAGAL3 gene. The function of GhAGAL3 was verified through virus-induced gene silencing (VIGS). As a result, GhAGAL3 gene silencing resulted in milder wilting of seedlings than the controls, and a significant increase in the raffinose content in cotton, indicating that GhAGAL3 responded to NaCl stress. The increase in raffinose content improved the tolerance of cotton. Findings in this study lay an important foundation for further research on the role of the GhAGAL3 gene family in the molecular mechanism of abiotic stress resistance in cotton.

Identification and adaptive evolution analysis of glutaredoxin genes in Populus spp

Glutaredoxin (GRX) is a class of small redox proteins widely involved in cellular redox homeostasis and the regulation of various cellular processes. The role of GRX gene in the differentiation of Populus spp. is rarely reported. We compared the similarities and differences of GRX genes among four sections of poplar using bioinformatics, corrected the annotations of some GRX genes, and focused on analysing their transcript profiling and adaptive evolution in Populus spp. A total of 219 GRX genes were identified in four sections of poplar, among which annotations for 13 genes were corrected. Differences in GRX genes were found between sect. Turanga, represented by P. euphratica, and other poplar sections. Most notably, P. euphratica had the smallest number of duplication events for GRX genes (n = 9) and no tandem duplications, whereas there were >25 duplication events for all other poplars. Furthermore, we detected 18 pairs of GRX genes under positive selection pressure in various sections of poplar, and identified two groups of GRX genes in the Salicaceae that potentially underwent positive selection. Expression profiling results showed that the PtrGRX34 and its orthologous genes were upregulated under stress treatments. In summary, the GRX gene family underwent expansion during poplar differentiation, and some genes underwent rapid evolution during this process, which may be beneficial for Populus spp. to adapt to environmental changes. This study may provide more insights into the molecular mechanisms of Populus spp. adaptation to environmental changes and the adaptive evolution of GRX genes.

Glutaredoxin in Rice Growth, Development, and Stress Resistance: Mechanisms and Research Advances

Rice (Oryza sativa L.) is a staple food for more than half of the global population. Various abiotic and biotic stresses lead to accumulation of reactive oxygen species in rice, which damage macromolecules and signaling pathways. Rice has evolved a variety of antioxidant systems, including glutaredoxin (GRX), that protect against various stressors. A total of 48 GRX gene loci have been identified on 11 of the 12 chromosomes of the rice genome; none were found on chromosome 9. GRX proteins were classified into four categories according to their active sites: CPYC, CGFS, CC, and GRL. In this paper, we summarized the recent research advances regarding the roles of GRX in rice development regulation and response to stresses, and discussed future research perspectives related to rice production. This review could provide information for rice researchers on the current status of the GRX and serve as guidance for breeding superior varieties.

Isolation, Purification, and Application of Protoplasts and Transient Expression Systems in Plants

Protoplasts, derived from plant cells, exhibit remarkable totipotency and hold significant value across a wide spectrum of biological and biotechnological applications. These versatile applications encompass protein subcellular localization and interaction analysis, gene expression regulation, functional characterization, gene editing techniques, and single-cell sequencing. Protoplasts' usability stems from their inherent accessibility and their ability to efficiently incorporate exogenous genes. In this review, we provide a comprehensive overview, including details on isolation procedures and influencing factors, purification and viability assessment methodologies, and the utilization of the protoplast transient expression system. The aim is to provide a comprehensive overview of current applications and offer valuable insights into protoplast isolation and the establishment of transient expression systems in a diverse range of plant species, thereby serving as a valuable resource for the plant science community.

The Mining of Genetic Loci and the Analysis of Candidate Genes to Identify the Physical and Chemical Markers of Anti-Senescence in Rice

Premature senescence is a common occurrence in rice production, and seriously affects rice plants' nutrient utilization and growth. A total of 120 recombinant inbred lines (RILs) were obtained from successive self-crossing of F12 generations derived from Huazhan and Nekken2. The superoxide dismutase (SOD) activity, malondialdehyde (MDA), content and catalase (CAT) activity related to the anti-senescence traits and enzyme activity index of rice were measured for QTL mapping using 4858 SNPs. Thirteen QTLs related to anti-senescence were found, among which the highest LOD score was 5.70. Eighteen anti-senescence-related genes were found in these regions, and ten of them differed significantly between the parents. It was inferred that LOC_Os01g61500, LOC_Os01g61810, and LOC_Os04g40130 became involved in the regulation of the anti-senescence molecular network upon upregulation of their expression levels. The identified anti-senescence-related QTLs and candidate genes provide a genetic basis for further research on the mechanism of the molecular network that regulates premature senescence.

GhCYS2 governs the tolerance against cadmium stress by regulating cell viability and photosynthesis in cotton

Cysteine, an early sulfur-containing compound in plants, is of significant importance in sulfur metabolism. CYS encodes cysteine synthetase that further catalyzes cysteine synthesis. In this investigation, CYS genes, identified from genome-wide analysis of Gossypium hirsutum bioinformatically, led to the discovery of GhCYS2 as the pivotal gene responsible for Cd2+ response. The silencing of GhCYS2 through virus-induced gene silencing (VIGS) rendered plants highly susceptible to Cd2+ stress. Silencing GhCYS2 in plants resulted in diminished levels of cysteine and glutathione while leading to the accumulation of MDA and ROS within cells, thereby impeding the regular process of photosynthesis. Consequently, the stomatal aperture of leaves decreased, epidermal cells underwent distortion and deformation, intercellular connections are dramatically disrupted, and fissures manifested between cells. Ultimately, these detrimental effected culminating in plant wilting and a substantial reduction in biomass. The association established between Cd2+ and cysteine in this investigation offered a valuable reference point for further inquiry into the functional and regulatory mechanisms of cysteine synthesis genes.

Genome-wide analysis of the P450 gene family in tea plant (Camellia sinensis) reveals functional diversity in abiotic stress

BACKGROUND

Cytochrome P450 (Cytochrome P450s) genes are involved in the catalysis of various reactions, including growth, development, and secondary metabolite biosynthetic pathways. However, little is known about the characteristics and functions of the P450 gene family in Camellia sinensis (C. sinensis).

RESULTS

To reveal the mechanisms of tea plant P450s coping with abiotic stresses, analyses of the tea plant P450 gene family were conducted using bioinformatics-based methods. In total, 273 putative P450 genes were identified from the genome database of C. sinensis. The results showed that P450s were well-balanced across the chromosomes I to XV of entire genome, with amino acid lengths of 268-612 aa, molecular weights of 30.95-68.5 kDa, and isoelectric points of 4.93-10.17. Phylogenetic analysis divided CsP450s into 34 subfamilies, of which CYP71 was the most abundant. The predicted subcellular localization results showed that P450 was distributed in a variety of organelles, with chloroplasts, plasma membrane,,and cytoplasm localized more frequently. The promoter region of CsP450s contained various cis-acting elements related to phytohormones and stress responses. In addition, ten conserved motifs (Motif1-Motif10) were identified in the CsP450 family proteins, with 27 genes lacking introns and only one exon. The results of genome large segment duplication showed that there were 37 pairs of genes with tandem duplication. Interaction network analysis showed that CsP450 could interact with multiple types of target genes, and there are protein interactions within the family. Tissue expression analysis showed that P450 was highly expressed in roots and stems. Moreover, qPCR analysis of the relative expression level of the gene under drought and cold stress correlated with the sequencing results.

CONCLUSIONS

This study lays the foundation for resolving the classification and functional study of P450 family genes and provides a reference for the molecular breeding of C. sinensis.

Analysis of the histidine kinase gene family and the role of GhHK8 in response to drought tolerance in cotton

As an important member of the two-component system (TCS), histidine kinases (HKs) play important roles in various plant developmental processes and signal transduction in response to a wide range of biotic and abiotic stresses. So far, the HK gene family has not been investigated in Gossypium. In this study, a total of 177 HK gene family members were identified in cotton. They were further divided into seven groups, and the protein characteristics, genetic relationship, gene structure, chromosome location, collinearity, and cis-elements identification were comprehensively analyzed. Whole genome duplication (WGD) / segmental duplication may be the reason why the number of HK genes doubled in tetraploid Gossypium species. Expression analysis revealed that most cotton HK genes were mainly expressed in the reproductive organs and the fiber at initial stage. Gene expression analysis revealed that HK family genes are involved in cotton abiotic stress, especially drought stress and salt stress. In addition, gene interaction networks showed that HKs were involved in the regulation of cotton abiotic stress, especially drought stress. VIGS experiments have shown that GhHK8 is a negative regulatory factor in response to drought stress. Our systematic analysis provided insights into the characteristics of the HK genes in cotton and laid a foundation for further exploring their potential in drought stress resistance in cotton.

Genome-wide identification and expression pattern profiling of the ATP-binding cassette gene family in tea plant (Camelliasinensis)

The ATP-binding cassette (ABC) gene family is one of the largest and oldest protein families, consisting of ATP-driven transporters facilitating substrate transportation across cell membranes. However, little is known about the evolution and biological function of the ABC gene family in tea plants. In this study, we performed a genome-wide identification and expression analysis of genes encoding ABC transporter proteins in Camellia sinensis. Our analysis of 170 ABC genes revealed that CsABCs were unevenly distributed across 15 chromosomes, with an amino acid length ranging from 188 to 2489 aa, molecular weight ranging from 20.29 to 277.34 kDa, and an isoelectric point ranging from 4.89 to 10.63. Phylogenetic analysis showed that CsABCs were divided into eight subfamilies, among which the ABCG subfamily was the most abundant. Furthermore, the subcellular localization of CsABCs indicated that they were present in various organelles. Collinearity analysis between the tea plant and Arabidopsis thaliana genomes revealed that the CsABC genes were homologous to the AtABC genes. Large gene fragment duplication analysis identified ten gene pairs as tandem repeats, and interaction network analysis demonstrated that CsABCs interacted with various types of target genes, with protein interactions also occurring within the family. Tissue expression analysis indicated that CsABCs were highly expressed in roots, stems, and leaves and were easily induced by drought and cold stress. Moreover, qRT-PCR analysis of the relative expression level of the gene under drought and cold stress correlated with the sequencing results. Identifying ABC genes in tea plants lays a foundation for the classification and functional analysis of ABC family genes, which can facilitate molecular breeding and the development of new tea varieties.

Genome wide identification and characterization of fertility associated novel CircRNAs as ceRNA reveal their regulatory roles in sheep fecundity

Reproductive traits play a vital role in determining the production efficiency of sheep. Maximizing the production is of paramount importance for breeders worldwide due to the growing population. Circular RNAs (circRNAs) act as miRNA sponges by absorbing miRNA activity through miRNA response elements (MREs) and participate in ceRNA regulatory networks (ceRNETs) to regulate mRNA expression. Despite of extensive research on role of circRNAs as miRNA sponges in various species, their specific regulatory roles and mechanism in sheep ovarian tissue are still not well understood. In this study, we performed whole genome sequencing of circRNAs, miRNA and mRNA employing bioinformatic techniques on ovine tissues of two contrasting sheep breeds "Small tail Han (X_LC) and Dolang sheep (D_LC)", which results into identification of 9,878 circRNAs with a total length of 23,522,667 nt and an average length of 2,381.32 nt. Among them, 44 differentially expressed circRNAs (DECs) were identified. Moreover, correlation between miRNA-mRNA and lncRNA-miRNA provided us with to prediction of miRNA binding sites on nine differentially expressed circRNAs and 165 differentially expressed mRNAs using miRanda. miRNA-mRNA and lncRNA-miRNA pairs with negative correlation were selected to determine the ceRNA score along with positively correlated pairs from lncRNA and mRNA network. Integration of ceRNA score and positively correlated pairs exhibit a significant ternary relationship among circRNAs-miRNA-mRNA demonestrated by ceRNA, comprising of 50 regulatory pairs sharring common nodes and predicted potential differentially expressed circRNAs-miRNAs-mRNAs regulatory axis. Based on functional enrichment analysis shortlisted key ceRNA regulatory pairs associated with reproduction including circRNA_3257-novel579_mature-EPHA3, circRNA_8396-novel130_mature-LOC101102473, circRNA_4140- novel34_mature > novel661_mature-KCNK9, and circRNA_8312-novel339_mature-LOC101110545. Furthermore, expression profiling, functional enrichments and qRT-PCR analysis of key target genes infer their implication in reproduction and metabolism. ceRNA target mRNAs evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations following genomic organizations will provide new insights underlying molecular mechanisms of reproduction, and establish a solid foundation for future research. Graphical abstract summarizing the scheme of study.

Four ciliate-specific expansion events occurred during actin gene family evolution of eukaryotes

Actin gene family is a divergent and ancient eukaryotic cellular cytoskeletal gene family, and participates in many essential cellular processes. Ciliated protists offer us an excellent opportunity to investigate gene family evolution, since their gene families evolved faster in ciliates than in other eukaryotes. Nonetheless, actin gene family is well studied in few model ciliate species but little is known about its evolutionary patterns in ciliates. Here, we analyzed the evolutionary pattern of eukaryotic actin gene family based on genomes/transcriptomes of 36 species covering ten ciliate classes, as well as those of nine non-ciliate eukaryotic species. Results showed: (1) Except for conventional actins and actin-related proteins (Arps) shared by various eukaryotes, at least four ciliate-specific subfamilies occurred during evolution of actin gene family. Expansions of Act2 and ArpC were supposed to have happen in the ciliate common ancestor, while expansions of ActI and ActII may have occurred in the ancestor of Armophorea, Muranotrichea, and Spirotrichea. (2) The number of actin isoforms varied greatly among ciliate species. Environmental adaptability, whole genome duplication (WGD) or segmental duplication events, distinct spatial and temporal patterns of expression might play driving forces for the increasement of isoform numbers. (3) The 'birth and death' model of evolution could explain the evolution of actin gene family in ciliates. And actin genes have been generally under strong negative selection to maintain protein structures and physiological functions. Collectively, we provided meaningful information for understanding the evolution of eukaryotic actin gene family.

Comprehensive evolutionary analysis of growth-regulating factor gene family revealing the potential molecular basis under multiple hormonal stress in Gramineae crops

Growth-regulating factors (GRFs) are plant-specific transcription factors that contain two highly conserved QLQ and WRC domains, which control a range of biological functions, including leaf growth, floral organ development, and phytohormone signaling. However, knowledge of the evolutionary patterns and driving forces of GRFs in Gramineae crops is limited and poorly characterized. In this study, a total of 96 GRFs were identified from eight crops of Brachypodium distachyon, Hordeum vulgare, Oryza sativa L. ssp. indica, Oryza rufipogon, Oryza sativa L. ssp. japonica, Setaria italic, Sorghum bicolor and Zea mays. Based on their protein sequences, the GRFs were classified into three groups. Evolutionary analysis indicated that the whole-genome or segmental duplication plays an essential role in the GRFs expansion, and the GRFs were negatively selected during the evolution of Gramineae crops. The GRFs protein function as transcriptional activators with distinctive structural motifs in different groups. In addition, the expression of GRFs was induced under multiple hormonal stress, including IAA, BR, GA3, 6BA, ABA, and MeJ treatments. Specifically, OjGRF11 was significantly induced by IAA at 6 h after phytohormone treatment. Transgenic experiments showed that roots overexpressing OjGRF11 were more sensitive to IAA and affect root elongation. This study will broaden our insights into the origin and evolution of the GRF family in Gramineae crops and will facilitate further research on GRF function.

Systematic analysis and expression of Gossypium 2ODD superfamily highlight the roles of GhLDOXs responding to alkali and other abiotic stress in cotton

BACKGROUND

2-oxoglutarate-dependent dioxygenase (2ODD) is the second largest family of oxidases involved in various oxygenation/hydroxylation reactions in plants. Many members in the family regulate gene transcription, nucleic acid modification/repair and secondary metabolic synthesis. The 2ODD family genes also function in the formation of abundant flavonoids during anthocyanin synthesis, thereby modulating plant development and response to diverse stresses.

RESULTS

Totally, 379, 336, 205, and 204 2ODD genes were identified in G. barbadense (Gb), G. hirsutum (Gh), G. arboreum (Ga), and G. raimondii (Gb), respectively. The 336 2ODDs in G. hirsutum were divided into 15 subfamilies according to their putative functions. The structural features and functions of the 2ODD members in the same subfamily were similar and evolutionarily conserved. Tandem duplications and segmental duplications served essential roles in the large-scale expansion of the cotton 2ODD family. Ka/Ks values for most of the gene pairs were less than 1, indicating that 2ODD genes undergo strong purifying selection during evolution. Gh2ODDs might act in cotton responses to different abiotic stresses. GhLDOX3 and GhLDOX7, two members of the GhLDOX subfamily from Gh2ODDs, were significantly down-regulated in transcription under alkaline stress. Moreover, the expression of GhLDOX3 in leaves was significantly higher than that in other tissues. These results will provide valuable information for further understanding the evolution mechanisms and functions of the cotton 2ODD genes in the future.

CONCLUSIONS

Genome-wide identification, structure, and evolution and expression analysis of 2ODD genes in Gossypium were carried out. The 2ODDs were highly conserved during evolutionary. Most Gh2ODDs were involved in the regulation of cotton responses to multiple abiotic stresses including salt, drought, hot, cold and alkali.

Genetic and evolutionary dissection of melatonin response signaling facilitates the regulation of plant growth and stress responses

The expansion of gene families during evolution could generate functional diversity among their members to regulate plant growth and development. Melatonin, a phylogenetically ancient molecule, is vital for many aspects of a plant's life. Understanding the functional diversity of the molecular players involved in melatonin biosynthesis, signaling, and metabolism will facilitate the regulation of plant phenotypes. However, the molecular mechanism of melatonin response signaling elements in regulating this network still has many challenges. Here, we provide an in-depth analysis of the functional diversity and evolution of molecular components in melatonin signaling pathway. Genetic analysis of multiple mutants in plant species will shed light on the role of gene families in melatonin regulatory pathways. Phylogenetic analysis of these genes was performed, which will facilitate the identification of melatonin-related genes for future study. Based on the abovementioned signal networks, the mechanism of these genes was summarized to provide reference for studying the regulatory mechanism of melatonin in plant phenotypes. We hope that this work will facilitate melatonin research in higher plants and finely tuned spatio-temporal regulation of melatonin signaling.

Systematic identification of TPS genes in Gossypium and their characteristics in response to flooding stress

Terpene synthases (TPS) is a key enzyme in the synthesis of plant terpenoids. Studies on TPSs have not been reported in Gossypium barbadense and Gossypium arboreum. 260 TPSs were identified in Gossypium, including 71 in Gossypium hirsutum, 75 in Gossypium. barbadense, 60 in Gossypium. arboreum, and 54 in Gossypium raimondii. We systematically analyzed the TPS gene family of Gossypium from three aspects: gene structure, evolutionary process and gene function. (1) Gene structure: Based on the protein structure of two conserved domains (PF01397 and PF03936), the TPS gene family is divided into five clades: TPS -a, -b, -c, -e/f and -g. (2) Evolution: Whole genome duplication and segmental duplication are the main modes of TPS gene amplification. (3) Function: The abundance of cis-acting elements may reveal the functional diversity of TPSs in cotton. TPS gene has tissue specific expression in cotton. The hypomethylation of the exon of TPSs may help to enhance the adaptability of cotton to flooding stress. In conclusion, this study can broaden the understanding of structure-evolution-function of the TPS gene family, and provide reference for the mining and verification of new genes.

Genome-wide analysis of the SBT gene family involved in drought tolerance in cotton

The subtilisin-like proteases (SBTs) are a large family of serine peptidases that are unique to plants. Previous studies have shown that SBTs are associated with developmental processes and environmental responses. However, comprehensive identification and systematic analysis of the SBT family have not been conducted in cotton. We used bioinformatics methods to analyze the structural characteristics, phylogenetic relationships, gene structures, expression modes, evolutionary relationships, selection pressures and stress responses of SBT gene family members in upland cotton. In this study, we identified 120 and 112 SBTs in the tetraploid cotton species G. hirsutum and G. barbadense, while 67 and 69 SBTs were identified in the diploid species G. arboreum and G. raimondii, respectively; these SBTs were divided into five distinct subfamilies. We identified the SBT gene GhSBT27A, and explore its function through virus-induced gene silencing and transmission electron microscopy. These results suggested that the GhSBT27A gene was involved in the response to drought stress. These results lay a foundation for further study on the drought stress mechanism of cotton.

GhAAO2 was observed responding to NaHCO3 stress in cotton compared to AAO family genes

BACKGROUND

Abscisic acid (ABA) is an important stress hormone, the changes of abscisic acid content can alter plant tolerance to stress, abscisic acid is crucial for studying plant responses to abiotic stress. The abscisic acid aldehyde oxidase (AAO) plays a vital role in the final step in the synthesis of abscisic acid, therefore, understanding the function of AAO gene family is of great significance for plants to response to abiotic stresses.

RESULT

In this study, 6, 8, 4 and 4 AAO genes were identified in four cotton species. According to the structural characteristics of genes and the traits of phylogenetic tree, we divided the AAO gene family into 4 clades. Gene structure analysis showed that the AAO gene family was relatively conservative. The analysis of cis-elements showed that most AAO genes contained cis-elements related to light response and plant hormones. Tissue specificity analysis under NaHCO₃ stress showed that GhAAO2 gene was differentially expressed in both roots and leaves. After GhAAO2 gene silencing, the degree of wilting of seedlings was lighter than that of the control group, indicating that GhAAO2 could respond to NaHCO₃ stress.

CONCLUSIONS

In this study, the AAO gene family was analyzed by bioinformatics, the response of GhAAO gene to various abiotic stresses was preliminarily verified, and the function of the specifically expressed gene GhAAO2 was further verified. These findings provide valuable information for the study of potential candidate genes related to plant growth and stress.

Identification, Classification and Characterization Analysis of FBXL Gene in Cotton

F-box/LR (FBXL), Leucine-rich repeats in F-box proteins, belongs to the Skp1-Cullin1-F-box protein (SCF) E3 ligase family. FBXL genes play important roles in plant growth, such as plant hormones, responses to environmental stress, and floral organ development. Here, a total of 518 FBXL genes were identified and analyzed in six plant species. Phylogenetic analysis showed that AtFBXLs, VvFBXLs, and GrFBXLs were clustered into three subfamilies (Ⅰ-Ⅲ). Based on the composition of the F-box domain and carboxyl-terminal amino acid sequence, FBXL proteins were classified into three types (Type-A/-B/-C). Whole-genome duplication (WGD) along with tandem duplications and segmental contributed to the expansion of this gene family. The result indicates that four cotton species are also divided into three subfamilies. FBXLs in cotton were classified into three clades by phylogenetic and structural analyses. Furthermore, expression analyses indicated that the expression patterns of GhFBXLs in different cotton tissues were different. The highly expressed of GH_A07G2363 in 5-8 mm anthers, indicates that this gene might play a role in the reproductive process, providing candidate genes for future studies on cotton fertility materials. This study provides an original functional opinion and a useful interpretation of the FBXL protein family in cotton.

A Comprehensive Analysis of the DUF4228 Gene Family in Gossypium Reveals the Role of GhDUF4228-67 in Salt Tolerance

Soil salinization conditions seriously restrict cotton yield and quality. Related studies have shown that the DUF4228 proteins are pivotal in plant resistance to abiotic stress. However, there has been no systematic identification and analysis of the DUF4228 gene family in cotton and their role in abiotic stress. In this study, a total of 308 DUF4228 genes were identified in four Gossypium species, which were divided into five subfamilies. Gene structure and protein motifs analysis showed that the GhDUF4228 proteins were conserved in each subfamily. In addition, whole genome duplication (WGD) events and allopolyploidization might play an essential role in the expansion of the DUF4228 genes. Besides, many stress-responsive (MYB, MYC) and hormone-responsive (ABA, MeJA) related cis-elements were detected in the promoters of the DUF4228 genes. The qRT-PCR results showed that GhDUF4228 genes might be involved in the response to abiotic stress. VIGS assays and the measurement of relative water content (RWC), Proline content, POD activity, and malondialdehyde (MDA) content indicated that GhDUF4228-67 might be a positive regulator of cotton response to salt stress. The results in this study systematically characterized the DUF4228s in Gossypium species and will provide helpful information to further research the role of DUF4228s in salt tolerance.

Genome-Wide Identification and Expression Analysis of the Zinc Finger Protein Gene Subfamilies under Drought Stress in Triticum aestivum

The zinc finger protein (ZFP) family is one of plants’ most diverse family of transcription factors. These proteins with finger-like structural domains have been shown to play a critical role in plant responses to abiotic stresses such as drought. This study aimed to systematically characterize Triticum aestivum ZFPs (TaZFPs) and understand their roles under drought stress. A total of 9 TaC2H2, 38 TaC3HC4, 79 TaCCCH, and 143 TaPHD were identified, which were divided into 4, 7, 12, and 14 distinct subgroups based on their phylogenetic relationships, respectively. Segmental duplication dominated the evolution of four subfamilies and made important contributions to the large-scale amplification of gene families. Syntenic relationships, gene duplications, and Ka/Ks result consistently indicate a potential strong purifying selection on TaZFPs. Additionally, TaZFPs have various abiotic stress-associated cis-acting regulatory elements and have tissue-specific expression patterns showing different responses to drought and heat stress. Therefore, these genes may play multiple functions in plant growth and stress resistance responses. This is the first comprehensive genome-wide analysis of ZFP gene families in T. aestivum to elucidate the basis of their function and resistance mechanisms, providing a reference for precise manipulation of genetic engineering for drought resistance in T. aestivum.

Expression Analysis of a Novel Oxidoreductase Glutaredoxin 2 in Black Tiger Shrimp, Penaeus monodon

Glutaredoxin (Grx) is a glutathione-dependent oxidoreductase that is an important component of the redox system in organisms. However, there is a serious lack of sequence information and functional validation related to Grx in crustaceans. In this study, a novel Grx was identified in Penaeus monodon (PmGrx2). The full-length cDNA of PmGrx2 is 998 bp, with an open reading frame (ORF) of 441 bp, encoding 119 amino acids. Sequence alignment showed that PmGrx2 had the highest identity with Grx2 of Penaeus vannamei at 96.64% and clustered with Grx2 of other crustaceans. Quantitative real-time PCR (qRT-PCR) analysis showed that PmGrx2 was expressed in all examined tissues, with higher expression levels in the stomach and testis. PmGrx2 was continuously expressed during development and had the highest expression level in the zygote stage. Both ammonia-N stress and bacterial infection could differentially induce the expression of PmGrx2 in hepatopancreas and gills. When PmGrx2 was inhibited, the expression of antioxidant enzymes was suppressed, the degree of apoptosis increased, and the GSH content decreased with the prolongation of ammonia-N stress. Inhibition of PmGrx2 resulted in shrimp being exposed to a greater risk of oxidative damage. In addition, an SNP locus was screened on the exons of PmGrx2 that was significantly associated with an ammonia-N-stress-tolerance trait. This study suggests that PmGrx2 is involved in redox regulation and plays an important role in shrimps’ resistance to marine environmental stresses.

Genome-wide identification, evolution and function analysis of UGTs superfamily in cotton

Glycosyltransferases mainly catalyse the glycosylation reaction in living organisms and widely exists in plants. UGTs have been identified from G. raimondii, G. arboreum and G. hirsutum. However, Genome-wide systematic analysis of UGTs superfamily have not been studied in G. barbadense. 752 UGTs were identified from four cotton species and grouped into 18 clades, of which R was newly discovered clades. Most UGTs were clustered at both ends of the chromosome and showed a heterogeneous distribution. UGT proteins were widely distributed in cells, with the highest distribution in chloroplasts. UGTs of the same clade shared similar intron/exon structural features. During evolution, the gene family has undergone strong selection for purification. UGTs were significantly enriched in “transcriptional activity (GO:0016758)” and “metabolic processes (GO:0008152)”. Genes from the same clade differed in function under various abiotic stresses. The analysis of cis-acting element and qRT–PCR may indicate that GHUGTs play important roles in plant growth, development and abiotic stress. We further found that GHUGT74-2 plays an important role under submergence. The study broadens the understanding of UGTs in terms of gene characteristics, evolutionary processes, and gene function in cotton and provides a new way to systematically and globally understand the structure–function relationship of multigene families in the evolutionary process.

Overexpression of cotton GhNAC072 gene enhances drought and salt stress tolerance in transgenic Arabidopsis

BACKGROUND

Crops face several environmental stresses (biotic and abiotic), thus resulting in severe yield losses. Around the globe abiotic stresses are the main contributors of plant damages, primarily drought and salinity. Many genes and transcription factors are involved in abiotic and biotic stress responses. NAC TF (Transcription Factors) improves tolerance to stresses by controlling the physiological and enzyme activities of crops.

RESULTS

In current research, GhNAC072 a highly upregulated TF in RNA-Seq was identified as a hub gene in the co-expression network analysis (WGCNA). This gene was transformed to Arabidopsis thaliana to confirm its potential role in drought and salt stress tolerance. Significant variations were observed in the morpho-physiological traits with high relative leaf water contents, chlorophyll contents, higher germination and longer root lengths of the overexpressed lines and low excised leaf loss and ion leakage as compared to the wildtype plants. Besides, overexpressed lines have higher amounts of antioxidants and low oxidant enzyme activities than the wildtype during the period of stress exposure.

CONCLUSIONS

In summary, the above analysis showed that GhNAC072 might be the true candidate involved in boosting tolerance mechanisms under drought and salinity stress.

Reactive oxygen species in plants: an invincible fulcrum for biotic stress mitigation

Climate change-associated environmental vagaries have amplified the incidence of pests and pathogens on plants, thus imparting the increased quest for management strategies. Plants respond to stresses through intricate signaling networks that regulate diverse cellular mechanisms. Reactive oxygen species (ROS) are cardinal towards the maintenance of normal plant activities as well as improving stress management. Plants that exhibit a fine balance between ROS levels and its management apparently mitigate stresses better. There have been very many compendiums on signaling and management of ROS during several abiotic stresses. However, expansion of knowledge related to ROS induction and homeostasis during biotic stresses is pertinent. Hence, considering its importance, we provide insights in this review on how plants signal and manage ROS upon an oxidative burst during their interaction with pathogens and herbivores. Substantial degree of molecular changes and pivotal roles of ROS have been detected during phyto-pathogen/herbivore interactions, opening novel platforms to understand signaling/management of events under varied biotic stresses. It is interesting to know that, though plants react to biotic stresses through oxidative burst, receptors and elicitors involved in the signal transduction differ across stresses. The review provides explicit details about the specific signaling of ROS production in plants under pathogen and herbivore attack. Furthermore, we also provide an update about tackling the accumulated ROS under biotic stresses as another pivotal step. ROS signaling and homeostasis can be exploited as critical players and a fulcrum to tackle biotic stresses, thus paving the way for futuristic combinatorial stress management strategies. KEY POINTS: • The review is a comprehension of redox signaling and management in plants during herbivory and pathogen infection • Reactive oxygen species (ROS) is an important factor during normal plant activities as well as in their response to stresses. Diverse modes of ROS signaling and management have been observed during both biotic stresses independently • Exploration of plant biology in multi-stress resistant plants like the crop wild relatives could pave the way for combinatorial management of stress for a better tomorrow.

LBD, de Sá MEL, Grossi-de-Sa MF, Mehta A. Overexpression of cotton genes GhDIR4 and GhPRXIIB in Arabidopsis thaliana improves plant resistance to root-knot nematode (Meloidogyne incognita) infection

Gossypium hirsutum L. represents the best cotton species for fiber production, thus computing the largest cultivated area worldwide. Meloidogyne incognita is a root-knot nematode (RKN) and one of the most important species of Meloidogyne genus, which has a wide host range, including cotton plants. Phytonematode infestations can only be partially controlled by conventional agricultural methods, therefore, more effective strategies to improve cotton resistance to M. incognita disease are highly desirable. The present study employed functional genomics to validate the involvement of two previously identified candidate genes, encoding dirigent protein 4-GhDIR4 and peroxiredoxin-2-GhPRXIIB, in cotton defense against M. incognita. Transgenic A. thaliana plant lines overexpressing GhDIR4 and GhPRXIIB genes were generated and displayed significantly improved resistance against M. incognita infection in terms of female nematode abundance in the roots when compared to wild-type control plants. For our best target-gene GhDIR4, an in-silico functional analysis, including multiple sequence alignment, phylogenetic relationship, and search for specific protein motifs unveiled potential orthologs in other relevant crop plants, including monocots and dicots. Our findings provide valuable information for further understanding the roles of GhDIR and GhPRXIIB genes in cotton defense response against RKN nematode.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03282-4.

Genome-Wide Expression Profiling and Networking Reveals an Imperative Role of IMF-Associated Novel CircRNAs as ceRNA in Pigs

Intramuscular fat (IMF) deposition is a biological process that has a strong impact on the nutritional and sensorial properties of meat, with relevant consequences on human health. Pork loins determine the effects of marbling on the sensory attributes and meat quality properties, which differ among various pig breeds. This study explores the crosstalk of non-coding RNAs with mRNAs and analyzes the potential pathogenic role of IMF-associated competing endogenous RNA (ceRNA) in IMF tissues, which offer a framework for the functional validation of key/potential genes. A high-throughput whole-genome transcriptome analysis of IMF tissues from longissimus dorsi muscles of Large White (D_JN) and Laiwu (L_JN) pigs resulted in the identification of 283 differentially expressed circRNAs (DECs), including two key circRNAs (circRNA-23437, circRNA-08840) with potential binding sites for multiple miRNAs regulating the whole network. The potential ceRNA mechanism identified the DEC target miRNAs-mRNAs involved in lipid metabolism, fat deposition, meat quality, and metabolic syndrome via the circRNA-miRNA-mRNA network, concluding that ssc-mir-370 is the most important target miRNA shared by both key circRNAs. TGM2, SLC5A6, ECI1, FASN, PER1, SLC25A34, SOD1, and COL5A3 were identified as hub genes through an intensive protein-protein interaction (PPI) network analysis of target genes acquired from the ceRNA regulatory network. Functional enrichments, pathway examinations, and qRT-PCR analyses infer their implications in fat/cholesterol metabolism, insulin secretion, and fatty acid biosynthesis. Here, circRNAs and miRNA sequencing accompanied by computational techniques were performed to analyze their expressions in IMF tissues from the longissimus dorsi muscles of two pig breeds. Their target gene evolutionary trajectories, expression profiling, functional enrichments, subcellular localizations, and structural advances with high-throughput protein modeling, following genomic organizations, will provide new insights into the underlying molecular mechanisms of adipocyte differentiation and IMF deposition and a much-needed qualitative framework for future research to improve meat quality and its role as a biomarker to treat lipid metabolic syndromes.

Structure and character analysis of cotton response regulator genes family reveals that GhRR7 responses to draught stress

Background

Cytokinin signal transduction is mediated by a two-component system (TCS). Two-component systems are utilized in plant responses to hormones as well as to biotic and abiotic environmental stimuli. In plants, response regulatory genes (RRs) are one of the main members of the two-component system (TCS).

Method

From the aspects of gene structure, evolution mode, expression type, regulatory network and gene function, the evolution process and role of RR genes in the evolution of the cotton genome were analyzed.

Result

A total of 284 RR genes in four cotton species were identified. Including 1049 orthologous/paralogous gene pairs were identified, most of which were whole genome duplication (WGD). The RR genes promoter elements contain phytohormone responses and abiotic or biotic stress-related cis-elements. Expression analysis showed that RR genes family may be negatively regulate and involved in salt stress and drought stress in plants. Protein regulatory network analysis showed that RR family proteins are involved in regulating the DNA-binding transcription factor activity (COG5641) pathway and HP kinase pathways. VIGS analysis showed that the GhRR7 gene may be in the same regulatory pathway as GhAHP5 and GhPHYB, ultimately negatively regulating cotton drought stress by regulating POD, SOD, CAT, H₂O₂ and other reactive oxygen removal systems.

Conclusion

This study is the first to gain insight into RR gene members in cotton. Our research lays the foundation for discovering the genes related to drought and salt tolerance and creating new cotton germplasm materials for drought and salt tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-022-00394-2.

Genome-wide identification, characterization and expression analysis of glutaredoxin gene family (Grxs) in Phaseolus vulgaris

Glutaredoxins (Grxs) are ubiquitous oxidoreductase proteins implicated in development and abiotic stress response mainly through maintaining redox homoeostasis. Here, we conducted the first systematic analysis of the Grx gene family (PvGrx) in the most popular legume Phaseolus vulgaris (common bean). A total of 50 PvGrx genes were identified, and divided into four classes (CC-type, CGFS-type, CPYC-type and Grl-type) based on the phylogenetic analysis. The different classes have different introns-exons structures and conserved motifs, indicating functional divergence in the PvGrx family. Both tandem and segmental duplications were found to be involved in the expansion of PvGrx family that underwent a purifying selection by excluding the deleterious loss-of-function mutations. Cis-acting regulatory elements and gene ontology analyses predicted their role of distinctive members in abiotic stress response and hormonal signalling. RNA-seq based expression analysis revealed their differential expression pattern during plant development. On the other hand, RT q-PCR analysis revealed that target PvGrx isoforms were associated with nodule organogenesis and symbiosis based on their expression profiles. In addition, a battery of PvGrx candidates were markedly upregulated by different abiotic stressors suggesting their broad spectrum of functions. These findings serve as a reference for functional analysis and genetic improvement in P. vulgaris and related legume species.

Genome-wide expression analysis of carboxylesterase (CXE) gene family implies GBCXE49 functional responding to alkaline stress in cotton

BACKGROUND

Carboxylesterase (CXE) is a type of hydrolase with α/β sheet hydrolase activity widely found in animals, plants and microorganisms, which plays an important role in plant growth, development and resistance to stress.

RESULTS

A total of 72, 74, 39, 38 CXE genes were identified in Gossypium barbadense, Gossypium hirsutum, Gossypium raimondii and Gossypium arboreum, respectively. The gene structure and expression pattern were analyzed. The GBCXE genes were divided into 6 subgroups, and the chromosome distribution of members of the family were mapped. Analysis of promoter cis-acting elements showed that most GBCXE genes contain cis-elements related to plant hormones (GA, IAA) or abiotic stress. These 6 genes we screened out were expressed in the root, stem and leaf tissues. Combined with the heat map, GBCXE49 gene was selected for subcellular locate and confirmed that the protein was expressed in the cytoplasm.

CONCLUSIONS

The collinearity analysis of the CXE genes of the four cotton species in this family indicated that tandem replication played an indispensable role in the evolution of the CXE gene family. The expression patterns of GBCXE gene under different stress treatments indicated that GBCXE gene may significantly participate in the response to salt and alkaline stress through different mechanisms. Through the virus-induced gene silencing technology (VIGS), it was speculated that GBCXE49 gene was involved in the response to alkaline stress in G. barbadense.

Exploration of glutathione reductase for abiotic stress response in bread wheat (Triticum aestivum L.)

A total of seven glutathione reductase (GR) genes were identified in Triticum aestivum, which were used for comparative structural characterization, phylogenetic analysis and expression profiling with the GR genes of other cereal plants. The modulated gene expression and enzyme activity revealed the role of GRs in abiotic stress response in T. aestivum. Glutathione reductase (GR) is an enzymatic antioxidant that converts oxidized glutathione (GSSG) into reduced glutathione (GSH) through the ascorbate-glutathione cycle. In this study, a total of seven GR genes forming two homeologous groups were identified in the allohexaploid genome of bread wheat (Triticum aestivum). Besides, we identified three GR genes in each Aegilops tauschii, Brachypodium distachyon, Triticum urartu and Sorghum bicolor, which were used for comparative characterization. Phylogenetic analysis revealed the clustering of GR proteins into two groups; class I and class II, which were predicted to be localized in cytoplasm and chloroplast, respectively. The exon-intron and conserved motif patterns were almost conserved in each group, in which a maximum of 10 and 17 exons were present in chloroplastic and cytoplasmic GRs, respectively. The protein structure analysis confirmed the occurrence of conserved pyridine nucleotide disulfide oxidoreductase (Pyr_redox) and pyridine nucleotide disulfide oxidoreductase dimerization (Pyr_redox_dim) domains in each GR. The active site of GR proteins consisted of two conserved cysteine residues separated by four amino acid residues. Promoter analysis revealed the occurrence of growth and stress-related cis-active elements. Tissue-specific expression profiling suggested the involvement of GRs in both vegetative and reproductive tissue development in various plants. The differential expression of TaGR genes and enhanced GR enzyme activity suggested their roles under drought, heat, salt and arsenic stress. Interaction of GRs with other proteins and chemical compounds of the ascorbate-glutathione cycle revealed their coordinated functioning. The current study will provide a foundation for the validation of the precise role of each GR gene in future studies.

Identification and Structure Analysis of KCS Family Genes Suggest Their Reponding to Regulate Fiber Development in Long-Staple Cotton Under Salt-Alkaline Stress

Plant 3-ketoacyl-CoA synthase (KCS) gene family catalyzed a β ketoacyl-CoA synthase, which was the rate-limiting enzyme for the synthesis of very long chain fatty acids (VLCFAs). Gossypium barbadense was well-known not only for high-quality fiber, which was perceived as a cultivated species of Gossypium. In this study, a total of 131 KCS genes were identified in four cotton species, there were 38, 44, 26, 23 KCS genes in the G. barbadense, the G. hirsutum, the G. arboreum and G. raimondii, respectively. The gene structure and expression pattern were analyzed. GBKCS genes were divided into six subgroups, the chromosome distribution of members of the family were mapped. The prediction of cis-acting elements of the GBKCS gene promoters suggested that the GBKCS genes may be involved in hormone signaling, defense and the stress response. Collinearity analysis on the KCS genes of the four cotton species were formulated. Tandem duplication played an indispensable role in the evolution of the KCS gene family. Specific expression analysis of 20 GBKCS genes indicated that GBKCS gene were widely expressed in the first 25 days of fiber development. Among them, GBKCS3, GBKCS8, GBKCS20, GBKCS34 were expressed at a high level in the initial long-term level of the G. barbadense fiber. This study established a foundation to further understanding of the evolution of KCS genes and analyze the function of GBKCS genes.

Identification of SNAT Family Genes Suggests GhSNAT3D Functional Reponse to Melatonin Synthesis Under Salinity Stress in Cotton

Serotonin N-acetyltransferase (SNAT) is a key enzyme in the biosynthesis of melatonin, and plays an important role in the regulation of melatonin synthesis. The study of SNAT is of great significance to understand the function of melatonin. In this study, we analyzed the structural characteristics, phylogenetic relationship, gene structure, expression pattern, evolutionary relationship and stress response of the members of the SNAT gene family in upland cotton through bioinformatics. A putative Serotonin n-acetyltransferase gene GhSNAT3D was identified, and preliminarily function of GhSNAT3D was verified by virus-induced gene silencing. We identified a total of 52 SNAT genes in the whole genome of G. hirsutum, and part of the GhSNATs were regulated by exogenous melatonin. The content of melatonin, antioxidant enzyme activity and Ca²⁺ content of GhSNAT3D gene silenced plants decreased, and the salt tolerance of GhSNAT3D gene silenced plants was reduced. Exogenous melatonin supplementation restored the salt tolerance of GhSNAT3D gene silenced plants. GhSNAT3D may interact with GhSNAT25D and ASMT to regulate melatonin synthesis. This study provided an important basis for further study on the regulation of melatonin in cotton against abiotic stress.

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 expression analysis of phospholipase A1 (PLA1) gene family suggests phospholipase A1-32 gene responding to abiotic stresses in cotton

Cotton is the most important crop for the production of natural fibres used in the textile industry. High salinity, drought, cold and high temperature represent serious abiotic stresses, which seriously threaten cotton production. Phospholipase A_S has an irreplaceable role in lipid signal transmission, growth and development and stress events. Phospholipase A can be divided into three families: PLA1, PLA2 and pPLA. Among them, the PLA1 family is rarely studied in plants. In order to study the potential functions of the PLA1 family in cotton, the bioinformatics analysis of the PLA1 family was correlated with cotton adversity, and tissue-specific analysis was performed. Explore the structure-function relationship of PLA1 members. It is found that the expression of GbPLA1-32 gene is affected by a variety of environmental stimuli, indicating that it plays a very important role in stress and hormone response, and closely associates the cotton adversity with this family. Through further functional verification, we found that virus-induced GbPLA1-32 gene silencing (VIGS) caused Gossypium barbadense to be sensitive to salt stress. This research provides an important basis for further research on the molecular mechanism of cotton resistance to abiotic stress.

Genome-wide identification and characteristic analysis of the downstream melatonin metabolism gene GhM2H in Gossypium hirsutum L

Background

Melatonin 2-hydroxylase (M2H) is the first enzyme in the catabolism pathway of melatonin, which catalyzes the production of 2-hydroxymelatonin (2-OHM) from melatonin. The content of 2-hydroxymelatonin in plants is much higher than that of melatonin. So M2H may be a key enzyme in the metabolic pathway of melatonin.

Method

We conducted a systematic analysis of the M2H gene family in Gossypium hirsutum based on the whole genome sequence by integrating the structural characteristics, phylogenetic relationships, expression profile, and biological stress of the members of the Gossypium hirsutum M2H gene family.

Result

We identified 265 M2H genes in the whole genome of Gossypium hirsutum, which were divided into 7 clades (clades I-VII) according to phylogenetic analysis. Most M2H members in each group had similar motif composition and gene structure characteristics. More than half of GhM2H members contain ABA-responsive elements and MeJA-responsive elements. Under different stress conditions, the expression levels of the gene changed, indicating that GhM2H members were involved in the regulation of abiotic stress. Some genes in the GhM2H family were involved in regulating melatonin levels in cotton under salt stress, and some genes were regulated by exogenous melatonin.

Conclusion

This study is helpful to explore the function of GhM2H, the downstream metabolism gene of melatonin in cotton, and lay the foundation for better exploring the molecular mechanism of melatonin improving cotton's response to abiotic stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40659-021-00358-y.

Identification and Analysis of GhEXO Gene Family Indicated That GhEXO7_At Promotes Plant Growth and Development Through Brassinosteroid Signaling in Cotton (Gossypium hirsutum L.)

Brassinosteroids (BRs), an efficient plant endogenous hormone, significantly promotes plant nutrient growth adapting to biological and abiotic adversities. BRs mainly promote plant cell elongation by regulating gene expression patterns. EXORDIUM (EXO) genes have been characterized as the indicators of BR response genes. Cotton, an ancient crop, is of great economic value and its fibers can be made into all kinds of fabrics. However, EXO gene family genes have not been full identified in cotton. 175 EXO genes were identified in nine plant species, of which 39 GhEXO genes in Gossypium hirsutum in our study. A phylogenetic analysis grouped all of the proteins encoded by the EXO genes into five major clades. Sequence identification of conserved amino acid residues among monocotyledonous and dicotyledonous species showed a high level of conservation across the N and C terminal regions. Only 25% the GhEXO genes contain introns besides conserved gene structure and protein motifs distribution. The 39 GhEXO genes were unevenly distributed on the 18 At and Dt sub-genome chromosomes. Most of the GhEXO genes were derived from gene duplication events, while only three genes showed evidence of tandem duplication. Homologous locus relationships showed that 15 GhEXO genes are located on collinear blocks and that all orthologous/paralogous gene pairs had Ka > Ks values, indicating purifying selection pressure. The GhEXO genes showed ubiquitous expression in all eight tested cotton tissues and following exposure to three phytohormones, IAA, GA, and BL. Furthermore, GhEXO7_At was mainly expressed in response to BL treatment, and was predominantly expressed in the fibers. GhEXO7_At was found to be a plasma membrane protein, and its ectopic expression in Arabidopsis mediated BR-regulated plant growth and development with altered expression of DWF4, CPD, KCS1, and EXP5. Additionally, the functions of GhEXO7_At were confirmed by virus-induced gene silencing (VIGS) in cotton. This study will provide important genetic resources for future cotton breeding programs.

Genome-wide identification and function analysis of HMAD gene family in cotton (Gossypium spp.)

BACKGROUND

The abiotic stress such as soil salinization and heavy metal toxicity has posed a major threat to sustainable crop production worldwide. Previous studies revealed that halophytes were supposed to tolerate other stress including heavy metal toxicity. Though HMAD (heavy-metal-associated domain) was reported to play various important functions in Arabidopsis, little is known in Gossypium.

RESULTS

A total of 169 G. hirsutum genes were identified belonging to the HMAD gene family with the number of amino acids ranged from 56 to 1011. Additionally, 84, 76 and 159 HMAD genes were identified in each G. arboreum, G. raimondii and G. barbadense, respectively. The phylogenetic tree analysis showed that the HMAD gene family were divided into five classes, and 87 orthologs of HMAD genes were identified in four Gossypium species, such as genes Gh_D08G1950 and Gh_A08G2387 of G. hirsutum are orthologs of the Gorai.004G210800.1 and Cotton_A_25987 gene in G. raimondii and G. arboreum, respectively. In addition, 15 genes were lost during evolution. Furthermore, conserved sequence analysis found the conserved catalytic center containing an anion binding (CXXC) box. The HMAD gene family showed a differential expression levels among different tissues and developmental stages in G. hirsutum with the different cis-elements for abiotic stress.

CONCLUSIONS

Current study provided important information about HMAD family genes under salt-stress in Gossypium genome, which would be useful to understand its putative functions in different species of cotton.

Genome-wide identification of CK gene family suggests functional expression pattern against Cd2+ stress in Gossypium hirsutum L

Choline kinase (CK) gene plays an important role in plants growth, development and resistance to stress. It mainly regulates the biosynthesis of phosphatidylcholine. This study aims to explore the structure-function relationship, and to provide a framework for functional validation and biochemical characterization of various CK genes. Our analysis showed that 87 CK genes were identified in cotton and 7 diploid plants, of which 43 genes encode CK proteins in 4 cotton species, and 13 genes were identified in Gossypium hirsutum. Most of GhCK genes are affected by the abiotic stress conditions, indicating the importance of CK proteins for plant development and response to abiotic stress. RT-qPCR analysis showed the tissue specificity of GhCK genes in response to Cd²⁺ and other abiotic stresses. Under Cd²⁺ stress, the expression level of GhCK gene family members has undergone different changes. The expression level of GhCK5 was enhanced, indicating that Cd²⁺ stress caused the increase of phosphatidylcholine content, which in turn reacted on the plant cell membrane, finally reached the absorption of Cd²⁺ into plant cells to repair Cd²⁺ the purpose of contaminated soil. This study will further broaden our understanding of the association between evolution and function of the GhCK gene family.

OPT gene family analysis of potato (Solanum tuberosum) responding to heavy metal stress: Comparative omics and co-expression networks revealed the underlying core templates and specific response patterns

Oligopeptides transporter (OPT) can maintain intracellular metal homeostat, however, their evolutionary characteristics, as well as their expression patterns in heavy metal exposure, remain unclear. Compared with previous OPT family identification, we identified 94 OPT genes (including 21 in potato) in potato and 4 other plants by HMMER program based on OPT domain (PF03169) for the first time. Secondly, conserved and special OPTs were found through comprehensive analysis. Thirdly, spatio-temporal tissue specific expression patterns and co-expression frameworks of potato OPT genes under different heavy metal stress were constructed. These data can provide excellent gene resources for food security and soil remediation.

Cotton DMP gene family: characterization, evolution, and expression profiles during development and stress

Members of DOMAIN OF UNKNOWN FUNCTION 679 membrane protein (DMP) gene family, a type of plant-specific membrane proteins, have been proposed to function in various physiological processes such as reproductive development and senescence in plants. Here, a total of 174 DMP genes were identified and analyzed in 16 plant species (including 58 DMPs in four cotton species). Phylogenetic analysis showed that these DMPs could be clustered into five subfamilies (I-V). 137 duplicated cotton gene pairs were identified and most duplicate events were formed by whole-genome duplication (WGD)/segmental duplications. Expression analysis revealed that most of cotton DMPs were mainly expressed in the reproductive organs (the sepal, petal, pistil and anther) and the fiber of secondary cell wall stage. GhDMPs promoter regions containing the different cis-elements also showed different responses to abiotic stress. In addition, gene interaction networks showed that DMPs, as an endomembrane system, were involved in plant senescence process and flower reproductive development. We speculated GhDMP8-A/-D, GbDMP8-A/-D could be used as some candidate gene for inducing cotton haploid. This genome-wide study provides a systematic analysis of the cotton DMP gene family, and further insights towards understanding the potential functions of candidate genes.

Genome-wide identification, characterization and expression profile of glutaredoxin gene family in relation to fruit ripening and response to abiotic and biotic stresses in banana (Musa acuminata)

Glutaredoxins (GRXs) are disulfide oxidoreductases that are involved in various biological processes. However, little information on the role of GRXs in the regulation of fruit ripening and the response to stress is available. In this study, we isolated 64 GRX genes from banana genome. Their encoded GRX proteins could be classified into four classes: CC, CGFS, CPYC and GRL types. The distribution and synteny of these GRXs on chromosomes, the gene structures, the promoter sequences, and the possible protein subcellular localizations were characterized. Molecular interaction network analysis suggested that MaGRX might interact with glutathione reductase (GR), sulfiredoxin, peroxiredoxin (Prx), and NADPH-dependent thioredoxin reductase C (NTRC), contributing to the antioxidative defense of banana fruit. MicroRNA prediction showed that MaGRX genes might be targeted by different miRNAs. Transcriptome analysis characterized the expression profiles of different MaGRX genes during banana fruit ripening, and in response to different storage stresses. The results suggested that CC-type, CPYC-type and GRL-type MaGRXs might be more active than CGFS-type MaGRXs during banana fruit ripening and the response to stress. Moreover, MaGRX6/7/9/11/17/23/28 and MaGRL3/16/19 might play important roles in regulating fruit ripening or in response to low and high temperature, or Fusarium proliferatum infection.