Positive roles of the Ca2+ sensors GbCML45 and GbCML50 in improving cotton Verticillium wilt resistance.

As a universal second messenger, cytosolic calcium (Ca2+) functions in multifaceted intracellular processes, including growth, development and responses to biotic/abiotic stresses in plant. The plant-specific Ca2+ sensors, calmodulin and calmodulin-like (CML) proteins, function as members of the second-messenger system to transfer Ca2+ signal into downstream responses. However, the functions of CMLs in the responses of cotton (Gossypium spp.) after Verticillium dahliae infection, which causes the serious vascular disease Verticillium wilt, remain elusive. Here, we discovered that the expression level of GbCML45 was promoted after V. dahliae infection in roots of cotton, suggesting its potential role in Verticillium wilt resistance. We found that knockdown of GbCML45 in cotton plants decreased resistance while overexpression of GbCML45 in Arabidopsis thaliana plants enhanced resistance to V. dahliae infection. Furthermore, there was physiological interaction between GbCML45 and its close homologue GbCML50 by using yeast two-hybrid and bimolecular fluorescence assays, and both proteins enhanced cotton resistance to V. dahliae infection in a Ca2+-dependent way in a knockdown study. Detailed investigations indicated that several defence-related pathways, including salicylic acid, ethylene, reactive oxygen species and nitric oxide signalling pathways, as well as accumulations of lignin and callose, are responsible for GbCML45- and GbCML50-modulated V. dahliae resistance in cotton. These results collectively indicated that GbCML45 and GbCML50 act as positive regulators to improve cotton Verticillium wilt resistance, providing potential targets for exploitation of improved Verticillium wilt-tolerant cotton cultivars by genetic engineering and molecular breeding.

GhHAM regulates GoPGF-dependent gland development and contributes to broad-spectrum pest resistance in cotton.

Cotton is a globally cultivated crop, producing 87% of the natural fiber used in the global textile industry. The pigment glands, unique to cotton and its relatives, serve as a defense structure against pests and pathogens. However, the molecular mechanism underlying gland formation and the specific role of pigment glands in cotton's pest defense are still not well understood. In this study, we cloned a gland-related transcription factor GhHAM and generated the GhHAM knockout mutant using CRISPR/Cas9. Phenotypic observations, transcriptome analysis, and promoter-binding experiments revealed that GhHAM binds to the promoter of GoPGF, regulating pigment gland formation in cotton's multiple organs via the GoPGF-GhJUB1 module. The knockout of GhHAM significantly reduced gossypol production and increased cotton's susceptibility to pests in the field. Feeding assays demonstrated that more than 80% of the cotton bollworm larvae preferred ghham over the wild type. Furthermore, the ghham mutants displayed shorter cell length and decreased gibberellins (GA) production in the stem. Exogenous application of GA3 restored stem cell elongation but not gland formation, thereby indicating that GhHAM controls gland morphogenesis independently of GA. Our study sheds light on the functional differentiation of HAM proteins among plant species, highlights the significant role of pigment glands in influencing pest feeding preference, and provides a theoretical basis for breeding pest-resistant cotton varieties to address the challenges posed by frequent outbreaks of pests.

Strigolactones positively regulate Verticillium wilt resistance in cotton via crosstalk with other hormones.

Verticillium wilt caused by Verticillium dahliae is a serious vascular disease in cotton (Gossypium spp.). V. dahliae induces the expression of the CAROTENOID CLEAVAGE DIOXYGENASE 7 (GauCCD7) gene involved in strigolactone (SL) biosynthesis in Gossypium australe, suggesting a role for SLs in Verticillium wilt resistance. We found that the SL analog rac-GR24 enhanced while the SL biosynthesis inhibitor TIS108 decreased cotton resistance to Verticillium wilt. Knock-down of GbCCD7 and GbCCD8b genes in island cotton (Gossypium barbadense) decreased resistance, whereas overexpression of GbCCD8b in upland cotton (Gossypium hirsutum) increased resistance to Verticillium wilt. Additionally, Arabidopsis (Arabidopsis thaliana) SL mutants defective in CCD7 and CCD8 putative orthologs were susceptible, whereas both Arabidopsis GbCCD7- and GbCCD8b-overexpressing plants were more resistant to Verticillium wilt than wild-type (WT) plants. Transcriptome analyses showed that several genes related to the jasmonic acid (JA)- and abscisic acid (ABA)-signaling pathways, such as MYELOCYTOMATOSIS 2 (GbMYC2) and ABA-INSENSITIVE 5, respectively, were upregulated in the roots of WT cotton plants in responses to rac-GR24 and V. dahliae infection but downregulated in the roots of both GbCCD7- and GbCCD8b-silenced cotton plants. Furthermore, GbMYC2 suppressed the expression of GbCCD7 and GbCCD8b by binding to their promoters, which might regulate the homeostasis of SLs in cotton through a negative feedback loop. We also found that GbCCD7- and GbCCD8b-silenced cotton plants were impaired in V. dahliae-induced reactive oxygen species (ROS) accumulation. Taken together, our results suggest that SLs positively regulate cotton resistance to Verticillium wilt through crosstalk with the JA- and ABA-signaling pathways and by inducing ROS accumulation.

CONSTITUTIVE EXPRESSER OF PATHOGENESIS-RELATED GENES 5 is an RNA-binding protein controlling plant immunity via an RNA processing complex.

Plant innate immunity is capable of combating diverse and ever evolving pathogens. The plasticity of innate immunity could be boosted by RNA processing. Arabidopsis thaliana CONSTITUTIVE EXPRESSER OF PATHOGENESIS-RELATED GENES 5 (CPR5), a key negative immune regulator, is a component of the nuclear pore complex. Here we further identified CPR5 as a component of RNA processing complexes. Through genetic screening, we found that RNA splicing activator NineTeen Complex and RNA polyadenylation factor CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR, coordinately function downstream of CPR5 to activate plant immunity. CPR5 and these two regulators form a complex that is localized in nuclear speckles, an RNA processing organelle. Intriguingly, we found that CPR5 is an RNA-binding protein belonging to the Transformer 2 (Tra2) subfamily of the serine/arginine-rich family. The RNA recognition motif of CPR5 protein binds the Tra2-targeted RNA sequence in vitro and is functionally replaceable by those of Tra2 subfamily proteins. In planta, it binds RNAs of CPR5-regulated alternatively spliced genes (ASGs) identified by RNA-seq. ARGONAUTE 1 (AGO1) is one of the ASGs and, consistent with this, the ago1 mutant suppresses the cpr5 phenotype. These findings reveal that CPR5 is an RNA-binding protein linking RNA processing with plant immunity.

lncRNA7 and lncRNA2 modulate cell wall defense genes to regulate cotton resistance to Verticillium wilt.

In plants, long noncoding RNAs (lncRNAs) regulate disease resistance against fungi and other pathogens. However, the specific mechanism behind this regulation remains unclear. In this study, we identified disease resistance-related lncRNAs as well as their regulating genes and assessed their functions by infection of cotton (Gossypium) chromosome segment substitution lines with Verticillium dahliae. Our results demonstrated that lncRNA7 and its regulating gene Pectin methylesterase inhibitor 13 (GbPMEI13) positively regulated disease resistance via the silencing approach, while ectopic overexpression of GbPMEI13 in Arabidopsis (Arabidopsis thaliana) promoted growth and enhanced resistance to V. dahliae. In contrast, lncRNA2 and its regulating gene Polygalacturonase 12 (GbPG12) negatively regulated resistance to V. dahliae. We further found that fungal disease-related agents, including the pectin-derived oligogalacturonide (OG), could downregulate the expression of lncRNA2 and GbPG12, leading to pectin accumulation. Conversely, OG upregulated the expression of lncRNA7, which encodes a plant peptide phytosulfokine (PSK-α), which was confirmed by lncRNA7 overexpression and Ultra Performance Liquid Chromatography Tandem Mass Spectrometry (UPLC-MS) experiments. We showed that PSK-α promoted 3-Indoleacetic acid (IAA) accumulation and activated GbPMEI13 expression through Auxin Response Factor 5. Since it is an inhibitor of pectin methylesterase (PME), GbPMEI13 promotes pectin methylation and therefore increases the resistance to V. dahliae. Consistently, we also demonstrated that GbPMEI13 inhibits the mycelial growth and spore germination of V. dahliae in vitro. In this study, we demonstrated that lncRNA7, lncRNA2, and their regulating genes modulate cell wall defense against V. dahliae via auxin-mediated signaling, providing a strategy for cotton breeding.

The use of ribosome-nascent chain complex-seq to reveal the translated mRNA profile and the role of ASN1 in resistance to Verticillium wilt in cotton.

We combined traditional mRNA-seq and RNC-seq together to reveal post-transcriptional regulation events impacting gene expression and interactions between the serious fungal pathogen Verticillium dahliae and a susceptible host, Gossypium hirsutum TM-1. After screening the differentially expressed and translated genes, V. dahliae infection was observed to influence gene transcription and translation in its host. Interestingly, the asparagine synthase (ASN1) gene transcripts increased significantly with the increase of infection time, while the rate of ASN1 protein accumulation in host TM-1 was distinctly lower than that in resistant hosts. We knocked down the ASN1 gene in resistant plants (ZZM2), and found that Verticillium-resistance was significantly reduced upon knockdown of ASN1. Our study revealed both transcriptional and post-transcriptional regulation of gene expression in TM-1 cotton plants infected by V. dahliae, and showed that ASN1 functions in the V. dahliae resistance process. These insights support breeding of disease resistance in cotton.

Large-fragment insertion activates gene GaFZ (Ga08G0121) and is associated with the fuzz and trichome reduction in cotton (Gossypium arboreum).

Cotton seeds are typically covered by lint and fuzz fibres. Natural 'fuzzless' mutants are an ideal model system for identifying genes that regulate cell initiation and elongation. Here, using a genome-wide association study (GWAS), we identified a ~ 6.2 kb insertion, larINDELFZ , located at the end of chromosome 8, composed of a ~ 5.0 kb repetitive sequence and a ~ 1.2 kb fragment translocated from chromosome 12 in fuzzless Gossypium arboreum. The presence of larINDELFZ was associated with a fuzzless seed and reduced trichome phenotypes in G. arboreum. This distant insertion was predicted to be an enhancer, located ~ 18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). Ectopic overexpression of GaFZ in Arabidopsis thaliana and G. hirsutum suggested that GaFZ negatively modulates fuzz and trichome development. Co-expression and interaction analyses demonstrated that GaFZ might impact fuzz fibre/trichome development by repressing the expression of genes in the very-long-chain fatty acid elongation pathway. Thus, we identified a novel regulator of fibre/trichome development while providing insights into the importance of noncoding sequences in cotton.

A copy number variant at the HPDA-D12 locus confers compact plant architecture in cotton.

Improving yield is a primary mission for cotton (Gossypium hirsutum) breeders; development of cultivars with suitable architecture for high planting density (HPDA) can increase yield per unit area. We characterized a natural cotton mutant, AiSheng98 (AS98), which exhibits shorter height, shorter branch length, and more acute branch angle than wild-type. A copy number variant at the HPDA locus on Chromosome D12 (HPDA-D12), encoding a dehydration-responsive element-binding (DREB) transcription factor, GhDREB1B, strongly affects plant architecture in the AS98 mutant. We found an association between a tandem duplication of a c. 13.5 kb segment in HPDA-D12 and elevated GhDREB1B expression resulting in the AS98 mutant phenotype. GhDREB1B overexpression confers a significant decrease in plant height and branch length, and reduced branch angle. Our results suggest that fine-tuning GhDREB1B expression may be a viable engineering strategy for modification of plant architecture favorable to high planting density in cotton.

Comprehensive Analysis of Respiratory Burst Oxidase Homologs (Rboh) Gene Family and Function of GbRboh5/18 on Verticillium Wilt Resistance in Gossypium barbadense.

Respiratory burst oxidase homologs (Rbohs) play a predominant role in reactive oxygen species (ROS) production, which is crucial in plant growth, differentiation, as well as their responses to biotic and abiotic stresses. To date, however, there is little knowledge about the function of cotton Rboh genes. Here, we identified a total of 87 Rbohs from five sequenced Gossypium species (the diploids Gossypium arboreum, Gossypium raimondii, and Gossypium australe, and the allotetraploids Gossypium hirsutum and Gossypium barbadense) via BLAST searching their genomes. Phylogenetic analysis of the putative 87 cotton Rbohs revealed that they were divided into seven clades. All members within the same clade are generally similar to each other in terms of gene structure and conserved domain arrangement. In G. barbadense, the expression levels of GbRbohs in the CladeD were induced in response to a fungal pathogen and to hormones (i.e., jasmonic acid and abscisic acid), based upon which the main functional member in CladeD was discerned to be GbRboh5/18. Further functional and physiological analyses showed that the knock-down of GbRboh5/18 expression attenuates plant resistance to Verticillium dahliae infection. Combined with the molecular and biochemical analyses, we found less ROS accumulation in GbRboh5/18-VIGS plants than in control plants after V. dahliae infection. Overexpression of GbRboh5/18 in G. barbadense resulted in more ROS accumulation than in control plants. These results suggest that GbRboh5/18 enhances the cotton plants' resistance against V. dahliae by elevating the levels of ROS accumulation. By integrating phylogenetic, molecular, and biochemical approaches, this comprehensive study provides a detailed overview of the number, phylogeny, and evolution of the Rboh gene family from five sequenced Gossypium species, as well as elucidating the function of GbRboh5/18 for plant resistance against V. dahliae. This study sheds fresh light on the molecular evolutionary properties and function of Rboh genes in cotton, and provides a reference for improving cotton's responses to the pathogen V. dahliae.

Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis.

The diploid wild cotton species Gossypium australe possesses excellent traits including resistance to disease and delayed gland morphogenesis, and has been successfully used for distant breeding programmes to incorporate disease resistance traits into domesticated cotton. Here, we sequenced the G. australe genome by integrating PacBio, Illumina short read, BioNano (DLS) and Hi-C technologies, and acquired a high-quality reference genome with a contig N50 of 1.83 Mb and a scaffold N50 of 143.60 Mb. We found that 73.5% of the G. australe genome is composed of various repeat sequences, differing from those of G. arboreum (85.39%), G. hirsutum (69.86%) and G. barbadense (69.83%). The G. australe genome showed closer collinear relationships with the genome of G. arboreum than G. raimondii and has undergone less extensive genome reorganization than the G. arboreum genome. Selection signature and transcriptomics analyses implicated multiple genes in disease resistance responses, including GauCCD7 and GauCBP1, and experiments revealed induction of both genes by Verticillium dahliae and by the plant hormones strigolactone (GR24), salicylic acid (SA) and methyl jasmonate (MeJA). Experiments using a Verticillium-resistant domesticated G. barbadense cultivar confirmed that knockdown of the homologues of these genes caused a significant reduction in resistance against Verticillium dahliae. Moreover, knockdown of a newly identified gland-associated gene GauGRAS1 caused a glandless phenotype in partial tissues using G. australe. The G. australe genome represents a valuable resource for cotton research and distant relative breeding as well as for understanding the evolutionary history of crop genomes.

Knockdown of Gh_A05G1554 (GhDHN_03) and Gh_D05G1729 (GhDHN_04) Dehydrin genes, Reveals their potential role in enhancing osmotic and salt tolerance in cotton.

In this investigation, whole-genome identification and functional characterization of the cotton dehydrin genes was carried out. A total of 16, 7, and 7 dehydrin proteins were identified in G. hirsutum, G. arboreum and G. raimondii, respectively. Through RNA sequence data and RT-qPCR validation, Gh_A05G1554 (GhDHN_03) and Gh_D05G1729 (GhDHN_04) were highly upregulated, and knockdown of the two genes, significantly reduced the ability of the plants to tolerate the effects of osmotic and salt stress. The VIGS-plants recorded significantly higher concentration levels of oxidants, hydrogen peroxide (H2O2) and malondialdehyde (MDA), furthermore, the four stress responsive genes GhLEA2, Gh_D12G2017 (CDKF4), Gh_A07G0747 (GPCR) and a transcription factor, trihelix, Gh_A05G2067, were significantly downregulated in VIGS-plants, but upregulated in wild types under osmotic and salt stress condition. The result indicated that dehydrin proteins are vital for plants and can be exploited in developing a more osmotic and salt stress-resilient germplasm to boost and improve cotton production.

Functional analysis of the GbDWARF14 gene associated with branching development in cotton.

Plant architecture, including branching pattern, is an important agronomic trait of cotton crops. In recent years, strigolactones (SLs) have been considered important plant hormones that regulate branch development. In some species such as Arabidopsis, DWARF14 is an unconventional receptor that plays an important role in the SL signaling pathway. However, studies on SL receptors in cotton are still lacking. Here, we cloned and analysed the structure of the GbD14 gene in Gossypium barbadense and found that it contains the domains necessary for a SL receptor. The GbD14 gene was expressed primarily in the roots, leaves and vascular bundles, and the GbD14 protein was determined via GFP to localize to the cytoplasm and nucleus. Gene expression analysis revealed that the GbD14 gene not only responded to SL signals but also was differentially expressed between cotton plants whose types of branching differed. In particular, GbD14 was expressed mainly in the axillary buds of normal-branching cotton, while it was expressed the most in the leaves of nulliplex-branch cotton. In cotton, the GbD14 gene can be induced by SL and other plant hormones, such as indoleacetic acid, abscisic acid, and jasmonic acid. Compared with wild-type Arabidopsis, GbD14-overexpressing Arabidopsis responded more rapidly to SL signals. Moreover, we also found that GbD14 can rescue the multi-branched phenotype of Arabidopsis Atd14 mutants. Our results indicate that the function of GbD14 is similar to that of AtD14, and GbD14 may be a receptor for SL in cotton and involved in regulating branch development. This research provides a theoretical basis for a profound understanding of the molecular mechanism of branch development and ideal plant architecture for cotton breeding improvements.

Overexpression of the Key Virulence Factor 1,3-1,4-ß-d-Glucanase in the Endophytic Bacterium Bacillus halotolerans Y6 To Improve Verticillium Resistance in Cotton.

Verticillium wilt, caused by Verticillium dahliae, results in a dramatic loss of cotton yields in China. There is great potential for biocontrol to manage this destructive crop disease. In this study, we obtained the endophytic bacterium Bacillus halotolerans Y6 from Verticillium wilt-resistant cotton Gossypium barbadense Xinhai15; this bacterium possesses strong antagonistic abilities that inhibit V. dahliae spore germination and mycelial growth. The results of the enzyme activity assay, heterologous expression, and gene knockdown showed that the key virulence factor of Y6 for antagonizing V. dahliae was ß -glucanase Bgy6. To facilitate field tests of biological control, we constructed the homologous Bgy6-overexpression strain OY6. Compared with the wild-type Y6 strain, the ß-glucanase activity of OY6 was increased by 91.79%, and the inhibition rate of OY6 against V. dahliae V991 exceeded 96.7%. Moreover, the spores of V. dahliae V991 treated with OY6 showed more mucus and larger holes on the surface, as observed by scanning electron microscopy. Potting test results illustrated that both OY6 and Y6 could improve the resistance of upland cotton to Verticillium wilt. With the inoculation of V. dahliae V991 for 45 days, the disease index of G. hirsutum TM-1 treated with OY6 was only 8.33, which was significantly lower than that in plants treated with the wild-type strain Y6 (17.86) or the controls without bacteria (35.94). Our research provides a new idea for the control of Verticillium wilt in upland cotton via transforming endophytic bacteria of Verticillium wilt-resistant cotton and proposes a new solution to prevent and control Verticillium wilt.

Identifying favorable alleles for improving key agronomic traits in upland cotton.

BACKGROUND: Gossypium hirsutum L. is grown worldwide and is the largest source of natural fiber crop. We focus on exploring the favorable alleles (FAs) for upland cotton varieties improvement, and further understanding the history of accessions selection and acumination of favorable allele during breeding. RESULTS: The genetic basis of phenotypic variation has been studied. But the accumulation of favorable alleles in cotton breeding history in unknown, and potential favorable alleles to enhance key agronomic traits in the future cotton varieties have not yet been identified. Therefore, 419 upland cotton accessions were screened, representing a diversity of phenotypic variations of 7362 G. hirsutum, and 15 major traits were investigated in 6 environments. These accessions were categorized into 3 periods (early, medium, and modern) according to breeding history. All accessions were divided into two major groups using 299 polymorphic microsatellite markers: G1 (high fiber yield and quality, late maturity) and G2 (low fiber yield and quality, early maturity). The proportion of G1 genotype gradually increased from early to modern breeding periods. Furthermore, 21 markers (71 alleles) were significantly associated (-log P > 4) with 15 agronomic traits in multiple environments. Seventeen alleles were identified as FAs; these alleles accumulated more in the modern period than in other periods, consistent with their phenotypic variation trends in breeding history. Our results demonstrate that the favorable alleles accumulated through breeding effects, especially for common favorable alleles. However, the potential elite accessions could be rapidly screened by rare favorable alleles. CONCLUSION: In our study, genetic variation and genome-wide associations for 419 upland cotton accessions were analyzed. Two favorable allele types were identified during three breeding periods, providing important information for yield/quality improvement of upland cotton germplasm.

A tuber mustard AP2/ERF transcription factor gene, BjABR1, functioning in abscisic acid and abiotic stress responses, and evolutionary trajectory of the ABR1 homologous genes in Brassica species.

The AP2/ERF superfamily of transcription factors is one of the largest transcription factor families in plants and plays an important role in plant development processes and stress responses. In this study, BjABR1, an AP2/ERF superfamily gene, from tuber mustard (Brassica juncea var. tumida Tsen et Lee), sharing high amino acid sequence similarity with the AtABR1 (Arabidopsis thaliana AP2-like abscisic acid repressor 1) gene, were performed functional research, and the ABR1 homologous genes in Brassica species were identified and performed phylogenetic analysis. The promoter sequence of BjABR1 contained many phytohormone- and stress-related cis-elements; ABA (abscisic acid) and abiotic stresses can induce BjABR1 expression in tuber mustard; overexpression of BjABR1 in Arabidopsis can alleviate plant sensitivity to ABA and salt and osmotic stresses, and the alleviation may be due to changes in stress/ABA-induced gene expression. These results indicated that BjABR1 functions in ABA and abiotic stress responses. By BLAST searches against the genome database of five Brassica species (three diploids, B. rapa, B. nigra, and B. oleracea, and two allotetraploid, B. juncea and B. napus) using the protein sequence of AtABR1, 3, 3, 3, 6, and 5 ABR1 homologous genes in B. nigra, B. rapa, B. oleracea, B. juncea, and B. napus were identified, respectively, and they shared high sequence similarity. By sequence analysis, annotation mistakes of the protein-coding regions of two ABR1 homologous genes, GSBRNA2T00134741001 and BjuB007684, were found and corrected. Then, the evolution analysis of these ABR1 homologous genes showed that the ancestor of the three diploid species had three ABR1 homologous genes and each diploid inherited all the three genes from their ancestor; then, allotetraploid B. juncea inherited all the six genes from B. rapa and B. nigra with no gene lost, while allotetraploid B. napus inherited all the three genes from B. oleracea and two genes from B. rapa with one gene lost, indicating that ABR1 homologous genes possessed greater hereditary conservation in Brassica species. The ABR1 homologous genes between B. rapa and B. oleracea shared much higher sequence similarity compared to that of B. nigra in diploid species, indicating that ABR1 homologous genes in B. nigra had experienced more rapid evolution, and B. rapa and B. oleracea may share closer relationship compared to B. nigra. Moreover, the spatial and temporal expression analysis of six ABR1 homologous genes of tuber mustard showed that they possessed different expression models. These results imply that ABR1 homologous genes are important to Brassica plants, and they may possess similar function in ABA and abiotic stress responses but play a role in different tissues and growing stages of plant. This study will provide the foundation to the functional research of ABR1 homologous genes in the Brassica species and help to reveal and understand the evolution mechanisms of Brassica species.

SSR-Linkage map of interspecific populations derived from Gossypium trilobum and Gossypium thurberi and determination of genes harbored within the segregating distortion regions.

Wild cotton species have significant agronomic traits that can be introgressed into elite cultivated varieties. The use of a genetic map is important in exploring, identification and mining genes which carry significant traits. In this study, 188 F2mapping individuals were developed from Gossypium thurberi (female) and Gossypium trilobum (male), and were genotyped by using simple sequence repeat (SSR) markers. A total of 12,560 simple sequence repeat (SSR) markers, developed by Southwest University, thus coded SWU were screened out of which only 994 were found to be polymorphic, and 849 markers were linked in all the 13 chromosomes. The map had a length of 1,012.458 cM with an average marker distance of 1.193 cM. Segregation distortion regions (SDRs) were observed on Chr01, Chr02, Chr06, Chr07 Chr09, Chr10 and Chr11 with a large proportion of the SDR regions segregating towards the heterozygous allele. There was good syntenic block formation that revealed good collinearity between the genetic and physical map of G. raimondii, compared to the Dt_sub genome of the G. hirsutum and G. barbadense. A total of 2,496 genes were mined within the SSR related regions. The proteins encoding the mined genes within the SDR had varied physiochemical properties; their molecular weights ranged from 6.586 to 252.737 kDa, charge range of -39.5 to 52, grand hydropathy value (GRAVY) of -1.177 to 0.936 and isoelectric (pI) value of 4.087 to 12.206. The low GRAVY values detected showed that the proteins encoding these genes were hydrophilic in nature, a property common among the stress responsive genes. The RNA sequence analysis revealed more of the genes were highly upregulated in various stages of fiber development for instance; Gorai.002G241300 was highly up regulated at 5, 10, 20 and 25 day post anthesis (DPA). Validation through RT-qPCR further revealed that these genes mined within the SDR regions might be playing a significant role under fiber development stages, therefore we infer that Gorai.007G347600 (TFCA), Gorai.012G141600 (FOLB1), Gorai.006G024500 (NMD3), Gorai.002G229900 (LST8) and Gorai.002G235200 (NSA2) are significantly important in fiber development and in turn the quality, and further researches needed to be done to elucidate their exact roles in the fiber development process. The construction of the genetic map between the two wild species paves away for the mapping of quantitative trait loci (QTLs) since the average distance between the markers is small, and mining of genes on the SSR regions will provide an insight in identifying key genes that can be introgressed into the cultivated cotton cultivars.

Changes in DNA methylation assessed by genomic bisulfite sequencing suggest a role for DNA methylation in cotton fruiting branch development.

Cotton plant architecture, including fruit branch formation and flowering pattern, influences plant light exploitation, cotton yield and planting cost. DNA methylation has been widely observed at different developmental stages in both plants and animals and is associated with regulation of gene expression, chromatin remodelling, genome protection and other functions. Here, we investigated the global epigenetic reprogramming during the development of fruiting branches and floral buds at three developmental stages: the seedling stage, the pre-squaring stage and the squaring stage. We first identified 22 cotton genes which potentially encode DNA methyltransferases and demethylases. Among them, the homologous genes of CMT, DRM2 and MET1 were upregulated at pre-squaring and squaring stages, suggesting that DNA methylation is involved in the development of floral buds and fruit branches. Although the global methylation at all of three developmental stages was not changed, the CHG-type methylation of non-expressed genes was higher than those of expressed genes. In addition, we found that the expression of the homologous genes of the key circadian rhythm regulators, including CRY, LHY and CO, was associated with changes of DNA methylation at three developmental stages.

Identification of circularRNAs and their targets in Gossypium under Verticillium wilt stress based on RNA-seq.

Circular RNAs (circRNAs), a class of recently discovered non-coding RNAs, play a role in biological and developmental processes. A recent study showed that circRNAs exist in plants and play a role in their environmental stress responses. However, cotton circRNAs and their role in Verticillium wilt response have not been identified up to now. In this study, two CSSLs (chromosome segment substitution lines) of G.barbadense introgressed into G. hirsutum, CSSL-1 and CSSL-4 (a resistant line and a susceptible line to Verticillium wilt, respectively), were inoculated with V. dahliae for RNA-seq library construction and circRNA analysis. A total of 686 novel circRNAs were identified. CSSL-1 and CSSL-4 had similar numbers of circRNAs and shared many circRNAs in common. However, CSSL-4 differentially expressed approximately twice as many circRNAs as CSSL-1, and the differential expression levels of the common circRNAs were generally higher in CSSL-1 than in CSSL-4. Moreover, two C-RRI comparisons, C-RRI-vs-C-RRM and C-RRI-vs-C-RSI, possessed a large proportion (approximately 50%) of the commonly and differentially expressed circRNAs. These results indicate that the differentially expressed circRNAs may play roles in the Verticillium wilt response in cotton. A total of 280 differentially expressed circRNAs were identified. A Gene Ontology analysis showed that most of the 'stimulus response' term source genes were NBS family genes, of which most were the source genes from the differentially expressed circRNAs, indicating that NBS genes may play a role in Verticillium wilt resistance and might be regulated by circRNAs in the disease-resistance process in cotton.

Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species.

BACKGROUND: Nucleotide binding site (NBS) genes encode a large family of disease resistance (R) proteins in plants. The availability of genomic data of the two diploid cotton species, Gossypium arboreum and Gossypium raimondii, and the two allotetraploid cotton species, Gossypium hirsutum (TM-1) and Gossypium barbadense allow for a more comprehensive and systematic comparative study of NBS-encoding genes to elucidate the mechanisms of cotton disease resistance. RESULTS: Based on the genome assembly data, 246, 365, 588 and 682 NBS-encoding genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. The distribution of NBS-encoding genes among the chromosomes was nonrandom and uneven, and was tended to form clusters. Gene structure analysis showed that G. arboreum and G. hirsutum possessed a greater proportion of CN, CNL, and N genes and a lower proportion of NL, TN and TNL genes compared to that of G. raimondii and G. barbadense, while the percentages of RN and RNL genes remained relatively unchanged. The percentage changes among them were largest for TNL genes, about 7 times. Exon statistics showed that the average exon numbers per NBS gene in G. raimondii and G. barbadense were all greater than that in G. arboretum and G. hirsutum. Phylogenetic analysis revealed that the TIR-NBS genes of G. barbadense were closely related with that of G. raimondii. Sequence similarity analysis showed that diploid cotton G. arboreum possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. hirsutum, while diploid G. raimondii possessed a larger proportion of NBS-encoding genes similar to that of allotetraploid cotton G. barbadense. The synteny analysis showed that more NBS genes in G. raimondii and G. arboreum were syntenic with that in G. barbadense and G. hirsutum, respectively. CONCLUSIONS: The structural architectures, amino acid sequence similarities and synteny of NBS-encoding genes between G. arboreum and G. hirsutum, and between G. raimondii and G. barbadense were the highest among comparisons between the diploid and allotetraploid genomes, indicating that G. hirsutum inherited more NBS-encoding genes from G. arboreum, while G. barbadense inherited more NBS-encoding genes from G. raimondii. This asymmetric evolution of NBS-encoding genes may help to explain why G. raimondii and G. barbadense are more resistant to Verticillium wilt, whereas G. arboreum and G. hirsutum are more susceptible to Verticillium wilt. The disease resistances of the allotetraploid cotton were related to their NBS-encoding genes especially in regard from which diploid progenitor they were derived, and the TNL genes may have a significant role in disease resistance to Verticillium wilt in G. raimondii and G. barbadense.

Genome-wide identification of the TIFY gene family in three cultivated Gossypium species and the expression of JAZ genes.

TIFY proteins are plant-specific proteins containing TIFY, JAZ, PPD and ZML subfamilies. A total of 50, 54 and 28 members of the TIFY gene family in three cultivated cotton species-Gossypium hirsutum, Gossypium barbadense and Gossypium arboretum-were identified, respectively. The results of phylogenetic analysis showed that these TIFY genes were divided into eight clusters. The different clusters of gene family members often have similar gene structures, including the number of exons. The results of quantitative reverse transcription polymerase chain reaction (qRT-PCR) showed that different JAZ genes displayed distinct expression patterns in the leaves of upland cotton under treatment with Gibberellin (GA), methyl jasmonate (MeJA), Jasmonic acid (JA) and abscisic acid (ABA). Different groups of JAZ genes exhibited different expression patterns in cotton leaves infected with Verticillium dahliae. The results of the comparative analysis of TIFY genes in the three cultivated species will be useful for understanding the involvement of these genes in development and stress resistance in cotton.