Search

SEP genes are famous for their function in the morphological novelty of bisexual flowers. Although the diverse functions of SEP genes were reported, only the regulatory mechanisms underlying floral organ development have been addressed. In this study, we identified SEP-like genes in Gossypium and found that SEP3 genes were duplicated in diploid cotton varieties. GhSEP4.1 and GhSEP4.2 were abundantly transcribed in the shoot apical meristem (SAM), but only GhSEP4.2 was expressed in the leaf vasculature. The expression pattern of GhSEPs in floral organs was conserved with that of homologs in Arabidopsis, except for GhSEP2 that was preponderantly expressed in ovules and fibers. The overexpression and silencing of each single GhSEP gene suggested their distinct role in promoting flowering via direct binding to GhAP1 and GhLFY genomic regions. The curly leaf and floral defects in overexpression lines with a higher expression of GhSEP genes revealed the concentration-dependent target gene regulation of GhSEP proteins. Moreover, GhSEP proteins were able to dimerize and interact with flowering time regulators. Together, our results suggest the dominant role of GhSEP4.2 in leaves to promote flowering via GhAP1-A04, and differently accumulated GhSEP proteins in the SAM alternately participate in forming the dynamic tetramer complexes to target at the different loci of GhAP1 and GhLFY to maintain reproductive growth. The regulatory roles of cotton SEP genes reveal their conserved and diversified functions.

Summary Verticillium wilt (VW) is a destructive disease that results in great losses in cotton yield and quality. Identifying genetic variation that enhances crop disease resistance is a primary objective in plant breeding. Here we reported a GWAS of cotton VW resistance in a natural‐variation population, challenged by different pathogenicity stains and different environments, and found 382 SNPs significantly associated with VW resistance. The associated signal repeatedly peaked in chromosome Dt11 (68 798 494‐69 212 808) containing 13 core elite alleles undescribed previously. The core SNPs can make the disease reaction type from susceptible to tolerant or resistant in accessions with alternate genotype compared to reference genotype. Of the genes associated with the Dt11 signal, 25 genes differentially expressed upon Verticillium dahliae stress, with 21 genes verified in VW resistance via gene knockdown and/or overexpression experiments. We firstly discovered that a gene cluster of L‐type lectin‐domain containing receptor kinase (GhLecRKs‐V.9) played an important role in VW resistance. These results proved that the associated Dt11 region was a major genetic locus responsible for VW resistance. The frequency of the core elite alleles (FEA) in modern varieties was significantly higher than the early/middle varieties (12.55% vs 4.29%), indicating that the FEA increased during artificial selection breeding. The current developmental resistant cultivars, JND23 and JND24, had fixed these core elite alleles during breeding without yield penalty. These findings unprecedentedly provided genomic variations and promising alleles for promoting cotton VW resistance improvement.

Cotton produces natural fiber for the textile industry. The genetic effects of genomic structural variations underlying agronomic traits remain unclear. Here, we generate two high-quality genomes of Gossypium hirsutum cv. NDM8 and Gossypium barbadense acc. Pima90, and identify large-scale structural variations in the two species and 1,081 G. hirsutum accessions. The density of structural variations is higher in the D-subgenome than in the A-subgenome, indicating that the D-subgenome undergoes stronger selection during species formation and variety development. Many structural variations in genes and/or regulatory regions potentially influencing agronomic traits were discovered. Of 446 significantly associated structural variations, those for fiber quality and Verticillium wilt resistance are located mainly in the D-subgenome and those for yield mainly in the A-subgenome. Our research provides insight into the role of structural variations in genotype-to-phenotype relationships and their potential utility in crop improvement., High-quality genomes of two cultivated tetraploid cottons Gossypium hirsutum cv. NDM8 and Gossypium barbadense acc. Pima90 and resequencing of 1,081 G. hirsutum accessions provide insights into the role of structural variations.

Cotton (Gossypium hirsutum) is constantly attacked by pathogens and insects. The most efficient control strategy is to develop resistant varieties using broad-spectrum gene resources. Several resistance loci harboured by superior varieties have been identified through genome-wide association studies. However, the key genes and/or loci have not been functionally identified. In this study, we identified a locus significantly associated with Verticillium wilt (VW) resistance, and within a 145.5-kb linkage disequilibrium, two non-specific lipid transfer protein genes (named GhnsLTPsA10) were highly expressed under Verticillium pathogen stress. The expression of GhnsLTPsA10 significantly increased in roots upon Verticillium dahliae stress but significantly decreased in leaves under insect attack. Furthermore, GhnsLTPsA10 played antagonistic roles in positively regulating VW and Fusarium wilt resistance and negatively mediating aphid and bollworm resistance in transgenic Arabidopsis and silenced cotton. By combining transcriptomic, histological and physiological analyses, we determined that GhnsLTPsA10-mediated phenylpropanoid metabolism further affected the balance of the downstream metabolic flux of flavonoid and lignin biosynthesis. The divergent expression of GhnsLTPsA10 in roots and leaves coordinated resistance of cotton against fungal pathogens and insects via the redirection of metabolic flux. In addition, GhnsLTPsA10 contributed to reactive oxygen species accumulation. Therefore, in this study, we elucidated the novel function of GhnsLTP and the molecular association between disease resistance and insect resistance, balanced by GhnsLTPsA10. This broadens our knowledge of the biological function of GhnsLTPsA10 in crops and provides a useful locus for genetic improvement of cotton.

Lipids are major and essential constituents of plant cells and provide energy for various metabolic processes. However, the function of the lipid signal in defence against Verticillium dahliae, a hemibiotrophic pathogen, remains unknown. Here, we characterized 19 conserved stearoyl‐ACP desaturase family proteins from upland cotton (Gossypium hirsutum). We further confirmed that GhSSI2 isoforms, including GhSSI2‐A, GhSSI2‐B, and GhSSI2‐C located on chromosomes A10, D10, and A12, respectively, played a dominant role to the cotton 18:1 (oleic acid) pool. Suppressing the expression of GhSSI2s reduced the 18:1 level, which autoactivated the hypersensitive response (HR) and enhanced cotton Verticillium wilt and Fusarium wilt resistance. We found that low 18:1 levels induced phenylalanine ammonia‐lyase‐mediated salicylic acid (SA) accumulation and activated a SA‐independent defence response in GhSSI2s‐silenced cotton, whereas suppressing expression of GhSSI2s affected PDF1.2‐dependent jasmonic acid (JA) perception but not the biosynthesis and signalling cascade of JA. Further investigation showed that structurally divergent resistance‐related genes and nitric oxide (NO) signal were activated in GhSSI2s‐silenced cotton. Taken together, these results indicate that SA‐independent defence response, multiple resistance‐related proteins, and elevated NO level play an important role in GhSSI2s‐regulated Verticillium wilt resistance. These findings broaden our knowledge regarding the lipid signal in disease resistance and provide novel insights into the molecular mechanism of cotton fungal disease resistance., GhSACPD gene family isoforms GhSSI2‐A, GhSSI2‐B, and GhSSI2‐C displayed an important role in regulating 18:1 fatty acid level and negatively regulated cotton disease resistance to Verticillium wilt.

Additional file 3: Table S3. The correlation (R2) of gene expression (log10 (FPKM+ 1)) between two biological replicates.

Additional file 2: Table S2. G. hirsutum fiber transcriptomes assembly statistics.

Additional file 1: Table S1. Primers used for qPCR.

Background Dirigent (DIR) proteins mediate regioselectivity and stereoselectivity during lignan biosynthesis and are also involved in lignin, gossypol and pterocarpan biosynthesis. This gene family plays a vital role in enhancing stress resistance and in secondary cell-wall development, but systematical understanding is lacking in cotton. Results In this study, 107 GbDIRs and 107 GhDIRs were identified in Gossypium barbadense and Gossypium hirsutum, respectively. Most of these genes have a classical gene structure without intron and encode proteins containing a signal peptide. Phylogenetic analysis showed that cotton DIR genes were classified into four distinct subfamilies (a, b/d, e, and f). Of these groups, DIR-a and DIR-e were evolutionarily conserved, and segmental and tandem duplications contributed equally to their formation. In contrast, DIR-b/d mainly expanded by recent tandem duplications, accompanying with a number of gene clusters. With the rapid evolution, DIR-b/d-III was a Gossypium-specific clade involved in atropselective synthesis of gossypol. RNA-seq data highlighted GhDIRs in response to Verticillium dahliae infection and suggested that DIR gene family could confer Verticillium wilt resistance. We also identified candidate DIR genes related to fiber development in G. barbadense and G. hirsutum and revealed their differential expression. To further determine the involvement of DIR genes in fiber development, we overexpressed a fiber length-related gene GbDIR78 in Arabidopsis and validated its function in trichomes and hypocotyls. Conclusions These findings contribute novel insights towards the evolution of DIR gene family and provide valuable information for further understanding the roles of DIR genes in cotton fiber development as well as in stress responses.

As the largest cultivated fiber crop in the world, cotton (Gossypium hirsutum) is often exposed to various biotic stresses during its growth periods. Verticillium wilt caused by Verticillium dahliae is a severe disease in cotton, and the molecular mechanism of cotton resistance for Verticillium wilt needs to be further investigated. Here, we revealed that the cotton genome contains nine types of GST genes. An evolutionary analysis showed that a newly identified cluster (including Gh_A09G1508, Gh_A09G1509 and Gh_A09G1510) located on chromosome 09 of the A-subgenome was under positive selection pressure during the formation of an allotetraploid. Transcriptome analysis showed that this cluster participates in Verticillium wilt resistance. Because the Gh_A09G1509 gene showed the greatest differential expression in the resistant cultivar under V. dahliae stress, we overexpressed this gene in tobacco and found that its overexpression resulted in enhanced Verticillium wilt resistance. Suppression of the gene cluster via virus-induced gene silencing made cotton plants of the resistant cultivar Nongda601 significantly susceptible. These results demonstrated that the GST cluster played an important role in Verticillium wilt resistance. Further investigation showed that the encoded enzymes of the cluster were essential for the delicate equilibrium between the production and scavenging of H2 O2 during V. dahliae stress.

Background Evaluating phenotypic traits is very important for the selection of elite lines in Gossypium hirsutum L. Cotton breeders are interested in using diverse genotypes in hybridization that can segregate for traits of interested with the possibility of selection and genetic gain. Information on phenotypic and molecular diversity helps the breeders for parental selection. Methods In this study, 719 global collections of G. hirsutum L. were evaluated for five fibre-related traits during two consecutive years in eight different environments. A series of phenotypic data for fibre quality traits were obtained and the elite accessions were further screened using principal component analysis and phylogenetic tree construction based on single nucleotide polymorphism markers. Results We found that fibre quality traits showed a wide range of variation among the G. hirsutum accessions over 2 years. In general, accessions from outside China tended to have higher fibre length (FL) and fibre strength (FS) than did Chinese accessions. Among different regional accessions in China, North/Northwest accessions tended to have the highest FL, FS and best fibre macronaire. By assessing five fibre quality traits over 2 years with genotypic data, 31 elite germplasms reaching double-thirty quality values (FL ≥ 30 mm and FS ≥ 30 cN·tex− 1) were selected. Conclusions This study provided a detailed phenotypic diversity description of a population representing a wide range of upland cotton germplasm. Our findings provide useful information about possible elite fibre quality parents for cotton breeding programs.

Background Cyclophilins (CYPs), belonging to the peptidyl prolyl cis/trans isomerase (PPIase) superfamily, play important roles during plant responses to biotic and abiotic stresses. Results Here, a total of 79 CYPs were identified in the genome of Gossypium hirsutum. Of which, 65 GhCYPs only contained one cyclophilin type PPIase domain, others 14 GhCYPs contain additional domains. A number of cis-acting elements related to phytohormone signaling were predicated in the upstream of GhCYPs ORF. The expression analysis revealed that GhCYPs were induced in response to cold, hot, salt, PEG and Verticillium dahliae infection. In addition, the functional importance of GhCYP-3 in Verticillium wilt resistance was also presented in this study. GhCYP-3 showed both cytoplasmic and nuclear localization. Overexpression of GhCYP-3 in Arabidopsis significantly improved Verticillium wilt resistance of the plants. Recombinant GhCYP-3 displayed PPIase activity and evident inhibitory effects on V. dahliae in vitro. Moreover, the extracts from GhCYP-3 transgenic Arabidopsis displayed significantly inhibit activity to conidia germinating and hyphal growth of V. dahliae. Conclusions Our study identified the family members of cotton CYP genes using bioinformatics tools. Differential expression patterns of GhCYPs under various abiotic stress and V. dahliae infection conditions provide a comprehensive understanding of the biological functions of candidate genes. Moreover, GhCYP-3 involved in the resistance of cotton to V. dahliae infection presumably through antifungal activity. Electronic supplementary material The online version of this article (10.1186/s12870-019-1848-1) contains supplementary material, which is available to authorized users.

Salt tolerance in cotton is highly imperative for improvement in the response to decreasing farmland and soil salinization. However, little is known about the genetic basis underlying salt tolerance in cotton, especially the seedling stage. In this study, we evaluated two salt-tolerance-related traits of a natural population comprising 713 upland cotton (Gossypium hirsutum L.) accessions worldwide at the seedling stage and performed a genome-wide association study (GWAS) to identify marker-trait associations under salt stress using the Illumina Infinium CottonSNP63K array. A total of 23 single nucleotide polymorphisms (SNPs) that represented seven genomic regions on chromosomes A01, A10, D02, D08, D09, D10, and D11 were significantly associated with the two salt-tolerance-related traits, relative survival rate (RSR) and salt tolerance level (STL). Of these, the two SNPs i46598Gh and i47388Gh on D09 were simultaneously associated with the two traits. Based on all loci, we screened 280 possible candidate genes showing different expression levels under salt stress. Most of these genes were involved in transcription factors, transporters and enzymes and were previously reported as being involved in plant salt tolerance, such as NAC, MYB, NXH, WD40, CDPK, LEA, and CIPK. We further validated six putative candidate genes by qRT-PCR and found a differential expression level between salt-tolerant and salt-sensitive varieties. Our findings provide valuable information for enhancing the understanding of complicated mechanisms of salt tolerance in G. hirsutum seedlings and cotton salt tolerance breeding by molecular marker-assisted selection.

A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield. Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

Summary Genetic improvement of fibre quality is one of the main breeding goals for the upland cotton, Gossypium hirsutum, but there are difficulties with precise selection of traits. Therefore, it is important to improve the understanding of the genetic basis of phenotypic variation. In this study, we conducted phenotyping and genetic variation analyses of 719 diverse accessions of upland cotton based on multiple environment tests and a recently developed Cotton 63K Illumina Infinium SNP array and performed a genome‐wide association study (GWAS) of fibre quality traits. A total of 10 511 polymorphic SNPs distributed in 26 chromosomes were screened across the cotton germplasms, and forty‐six significant SNPs associated with five fibre quality traits were detected. These significant SNPs were scattered over 15 chromosomes and were involved in 612 unique candidate genes, many related to polysaccharide biosynthesis, signal transduction and protein translocation. Two major haplotypes for fibre length and strength were identified on chromosomes Dt11 and At07. Furthermore, by combining GWAS and transcriptome analysis, we identified 163 and 120 fibre developmental genes related to length and strength, respectively, of which a number of novel genes and 19 promising genes were screened. These results provide new insight into the genetic basis of fibre quality in G. hirsutum and provide candidate SNPs and genes to accelerate the improvement of upland cotton.

Fiber strength is an essential trait of fiber property in cotton, and it is quantitatively inherited. Identification of stable quantitative trait loci (QTL) contributing to fiber strength would provide the key basis for marker-assisted selection (MAS) in cotton breeding. In this study, four interspecific hybridization populations were established with a common G. barbadense parent Pima 90-53 and two G. hirsutum parents (CCRI 8 and Handan 208), each of which had fiber strength characteristic. Based on the phenotypic data of fiber strength from seven environments, a stable QTL, qFS-c9-1, was detected and validated on c9 in a marker interval between SSR markers NAU2395 and NAU1092. The QTL explaining 14.4-17.9 % of the phenotypic variation was firstly detected in two populations (CCRI 8 × Pima 90-53, BC1F1 and BC1F2) and its derived lines in four environments. And it accounting for 12.1-14.8 % of the phenotypic variation was further confirmed in two populations (Handan 208 × Pima 90-53, BC1F1, and F2) under one environment. In silico mapping using three sequenced cotton genomes indicated that homologous genes, anchored by NAU2395 and NAU1092, were aligned to the G. arboreum genome within a physical distance between 81.10 Mbps and 87.07 Mbps. In that interval, several genes were confirmed in literatures to associate with fiber development. Among these genes, seven genes were further selected for an expression analysis through fiber development transcriptome database, revealing unique expression patterns across different stages of fiber development between CCRI 8 and Pima 90-53. The genes underlying qFS-c9-1 were favorable to fine mapping and cloning. The current study results provided valuable evidence for mapping stable QTL of fiber strength utilizing multiple populations and environments, as well as map-based cloning the candidate gene underlying the QTL for future prospective research directions.