Cotton functional genomics reveals global insight into genome evolution and fiber development

Unveiling the power of MYB transcription factors: Master regulators of multi-stress responses and development in cotton

Plants, being immobile, are subject to environmental stresses more than other creatures, necessitating highly effective stress tolerance systems. Transcription factors (TFs) play a crucial role in the adaptation mechanism as they can be activated by diverse signals and ultimately control the expression of stress-responsive genes. One of the most prominent plant TFs family is MYB (myeloblastosis), which is involved in secondary metabolites, developmental mechanisms, biological processes, cellular architecture, metabolic pathways, and stress responses. Extensive research has been conducted on the involvement of MYB TFs in crops, while their role in cotton remains largely unexplored. We also utilized genome-wide data to discover potential 440 MYB genes and investigated their plausible roles in abiotic and biotic stress conditions, as well as in different tissues across diverse transcriptome databases. This review primarily summarized the structure and classification of MYB TFs biotic and abiotic stress tolerance and their role in secondary metabolism in different crops, especially in cotton. However, it intends to identify gaps in current knowledge and emphasize the need for further research to enhance our understanding of MYB roles in plants.

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis

BACKGROUND

GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism.

RESULTS

A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner.

CONCLUSIONS

Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Origin, evolution, and diversification of the wall-associated kinase gene family in plants

KEY MESSAGE

The study of the origin, evolution, and diversification of the wall-associated kinase gene family in plants facilitates their functional investigations in the future. Wall-associated kinases (WAKs) make up one subfamily of receptor-like kinases (RLKs), and function directly in plant cell elongation and responses to biotic and abiotic stresses. The biological functions of WAKs have been extensively characterized in angiosperms; however, the origin and evolutionary history of the WAK family in green plants remain unclear. Here, we performed a comprehensive analysis of the WAK family to reveal its origin, evolution, and diversification in green plants. In total, 1061 WAK genes were identified in 37 species from unicellular algae to multicellular plants, and the results showed that WAK genes probably originated before bryophyte differentiation and were widely distributed in land plants, especially angiosperms. The phylogeny indicated that the land plant WAKs gave rise to five clades and underwent lineage-specific expansion after species differentiation. Cis-acting elements and expression patterns analyses of WAK genes in Arabidopsis and rice demonstrated the functional diversity of WAK genes in these two species. Many gene gains and losses have occurred in angiosperms, leading to an increase in the number of gene copies. The evolutionary trajectory of the WAK family during polyploidization was uncovered using Gossypium species. Our results provide insights into the evolution of WAK genes in green plants, facilitating their functional investigations in the future.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

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.

The hexokinase Gene Family in Cotton: Genome-Wide Characterization and Bioinformatics Analysis

Hexokinase (HXK) is involved in hexose phosphorylation, sugar sensing, and signal transduction, all of which regulate plant growth and adaptation to stresses. Gossypium hirsutum L. is one of the most important fiber crops in the world, however, little is known about the HXKs gene family in G. hirsutum L. We identified 17 GhHXKs from the allotetraploid G. hirsutum L. genome (AADD). G. raimondii (DD) and G. arboreum (AA) are the diploid progenitors of G. hirsutum L. and contributed equally to the At_genome and Dt_genome GhHXKs genes. The chromosomal locations and exon-intron structures of GhHXK genes among cotton species are conservative. Phylogenetic analysis grouped the HXK proteins into four and three groups based on whether they were monocotyledons and dicotyledons, respectively. Duplication event analysis demonstrated that HXKs in G. hirsutum L. primarily originated from segmental duplication, which prior to diploid hybridization. Experiments of qRT-PCR, transcriptome and promoter cis-elements demonstrated that GhHXKs' promoters have auxin and GA responsive elements that are highly expressed in the fiber initiation and elongation stages, while the promoters contain ABA-, MeJA-, and SA-responsive elements that are highly expressed during the synthesis of the secondary cell wall. We performed a comprehensive analysis of the GhHXK gene family is a vital fiber crop, which lays the foundation for future studies assessing its role in fiber development.

Biochemical evidence of epicuticular wax compounds involved in cotton-whitefly interaction

Sucking insects require a surface of plants on which the legs and the eggs of insects will adhere and to which insect mouthparts will access. The primary plant protection against insects is their surface property, which hinders the attachment of the insect’s legs and eggs. The epicuticular waxes chemistry influences the fine structure of the cuticular surface. In current study, an attempt was made to investigate the variation of chemical compounds in epicuticular waxes of four cotton species that classify them resistant or susceptible i.e., Gossypium abroreum, G. hirsutum, G. arboreum wax deficient mutant (GaWM3) and G. harknessi which were evaluated for their interaction with whitefly and CLCuV transmission. Gossypium hirsutum an insect and CLCuV susceptible cotton variety, was found to have four compounds namely Trichloroacetic acid, hexadecylester, P-xylenolpthalein, 2-cyclopentene-1-ol, 1-phenyl-and Phenol, 2,5-bis [1,1- dimethyl] which could interact with chitin of whitefly while only two compounds in Gossypium arboreum an insect and CLCuV resistant cotton variety could interact with chitin of whitefly. Similarly, GaWM3 and Gossypium harkasnessi were found to have only a single compound. Number of whiteflies found on leaves of G. hirsutum was much higher as compared to other cotton species. Keeping this fact in mind a wax biosynthetic gene CER3, from Arabidopsis thaliana was transformed into G. hirsutum and the plants were evaluated for their resistance against whitefly and CLCuV transmission. In microscopic analysis transgenic plants clearly showed higher amounts of leaf waxes as compared to non-transgenics. The least whitefly population and CLCuV titer of <10,000 units was found in transgenic plants compared to non-transgenic cotton where it was ≈4.5X106 units that confirmed the role of wax in insect interaction and ultimately to CLCuV transmission. This study provides novel insight on wax related compounds involved in cotton-whitefly interaction, which potentially can help in developing more efficient control strategies for this destructive pest.

Transformation and functional verification of Cry5Aa in cotton

Most of the cotton bollworm-resistant genes applied in cotton are more than 20 years and they all belong to Cry1Ab/c family, but the insect-resistant effects of Cry5Aa on cotton were rarely reported. The possible risk of resistance is increasing. The study synthesized a novel bollworm-resistant gene Cry5Aa artificially based on preferences of cotton codon. The new gene was transferred to cotton through the method of pollen tube pathway. The transgenic strains were identified by kanamycin test in field and laboratory PCR analysis. Meanwhile, an insect resistance test was conducted by artificial bollworm feeding with transgenic leaves and GK19 was used as a control in this study. Results showed that rate of positive transgenic strains with kanamycin resistance in the first generation (T1), the second generation (T2) and the third generation (T3) respectively were 7.76%, 73.1% and 95.5%. However, PCR analysis showed that the positive strain rate in T1, T2 and T3 were 2.35%, 55.8% and 94.5%, respectively. The resistant assay of cotton bollworm showed that the mortality rate of the second, third and fourth instar larva feed by the transgenic cotton leaves, were 85.42%, 73.35% and 62.79%, respectively. There was a significant difference between transgenic plant of Cry5Aa and GK19 in insect resistance. Finally, we also conducted the further analysis of gene expression patterns, gene flow and the effect on non-target pest in the study. The results showed that Cry5Aa gene had less environmental impact, and Cry5Aa has been transferred successfully and expressed stably in cotton. Therefore, the novel bollworm resistance gene can partially replace the current insect-resistance gene of Lepidoptera insects.

Sequencing Multiple Cotton Genomes Reveals Complex Structures and Lays Foundation for Breeding

Cotton is a major fiber plant, which provides raw materials for clothing, protecting humans from the harsh environment of cold or hot weathers, enriching the culture and custom of human societies. Due to its importance, the diploid and tetraploid genomes of different cotton plants have been repeatedly sequenced to obtain their complete and fine genome sequences. These valuable genome data sets revealed the evolutionary past of the cotton plants, which were recursively affected by polyploidization, with a decaploidization contributing to the formation of the genus Gossypium, and a neo-tetraploidization contributing to the formation of nowadays widely cultivated cotton plants. Post-polyploidization genome instability resulted in numerous structural changes of the genomes, such as gene loss, DNA inversion and translocation, illegitimate recombination, and accumulation of repetitive sequences, and functional innovation accompanied by elevated evolutionary rates of genes. Many these changes have been asymmetric between subgnomes of the tetraploid cottons, rendering their divergent profiles of biological regulation and function. The availability of whole-genome sequences has now paved the way to identify and clone functional genes, e.g., those relating to fiber development, and to enhance breeding efforts to cultivate cottons to produce high-yield and high-quality fibers, and to resist environmental and biological stress.

Genome sequence of Gossypium herbaceum and genome updates of Gossypium arboreum and Gossypium hirsutum provide insights into cotton A-genome evolution

Upon assembling the first Gossypium herbaceum (A1) genome and substantially improving the existing Gossypium arboreum (A2) and Gossypium hirsutum ((AD)1) genomes, we showed that all existing A-genomes may have originated from a common ancestor, referred to here as A0, which was more phylogenetically related to A1 than A2. Further, allotetraploid formation was shown to have preceded the speciation of A1 and A2. Both A-genomes evolved independently, with no ancestor-progeny relationship. Gaussian probability density function analysis indicates that several long-terminal-repeat bursts that occurred from 5.7 million years ago to less than 0.61 million years ago contributed compellingly to A-genome size expansion, speciation and evolution. Abundant species-specific structural variations in genic regions changed the expression of many important genes, which may have led to fiber cell improvement in (AD)1. Our findings resolve existing controversial concepts surrounding A-genome origins and provide valuable genomic resources for cotton genetic improvement.

Heterologous Expression of GbTCP4, a Class II TCP Transcription Factor, Regulates Trichome Formation and Root Hair Development in Arabidopsis

Two class I family teosinte branched1/cycloidea/proliferating cell factor1 (TCP) proteins from allotetraploid cotton are involved in cotton fiber cell differentiation and elongation and root hair development. However, the biological function of most class II TCP proteins is unclear. This study sought to reveal the characteristics and functions of the sea-island cotton class II TCP gene GbTCP4 by biochemical, genetic, and molecular biology methods. GbTCP4 protein localizes to nuclei, binding two types of TCP-binding cis-acting elements, including the one in its promoter. Expression pattern analysis revealed that GbTCP4 is widely expressed in tissues, with the highest level in flowers. GbTCP4 is expressed at different fiber development stages and has high transcription in fibers beginning at 5 days post anthesis (DPA). GbTCP4 overexpression increases primary root hair length and density and leaf and stem trichomes in transgenic Arabidopsis relative to wild-type plants (WT). GbTCP4 binds directly to the CAPRICE (CPC) promoter, increasing CPC transcript levels in roots and reducing them in leaves. Compared with WT plants, lignin content in the stems of transgenic Arabidopsis overexpressing GbTCP4 increased, and AtCAD5 gene transcript levels increased. These results suggest that GbTCP4 regulates trichome formation and root hair development in Arabidopsis and may be a candidate gene for regulating cotton fiber elongation.

Analysis of the chromatin binding affinity of retrotransposases reveals novel roles in diploid and tetraploid cotton

LTR-retrotransposable elements are major components of diploid (Gossypium arboreum) and tetraploid (Gossypium hirsutum) cotton genomes that have undergone dramatic increases in copy number during the course of evolution. However, little is known about the biological functions of LTR-retrotransposable elements in cotton. Here, we show that a copia-like LTR-retrotransposable element has maintained considerable activity in both G. arboreum and G. hirsutum. We identified two functional domains of the retrotransposon and analyzed their expression levels in various cotton tissues, including leaves, ovules, and germinating seeds. ChIP-qPCR (chromatin immunoprecipitation followed by quantitative PCR), using a copia-specific antibody, established that copia-like proteins primarily bind to the first exons of several protein-coding genes in cotton cells. This finding suggests that retrotransposons play a novel, important role in regulating the transcriptional activities of protein-coding genes with various biological activities.