Investigation of salt tolerance in cotton germplasm by analyzing agro-physiological traits and ERF genes expression

The development of genotypes that can tolerate high levels of salt is crucial for the efficient use of salt-affected land and for enhancing crop productivity worldwide. Therefore, incorporating salinity tolerance is a critical trait that crops must possess. Salt resistance is a complex character, controlled by multiple genes both physiologically and genetically. To examine the genetic foundation of salt tolerance, we assessed 16 F1 hybrids and their eight parental lines under normal and salt stress (15 dS/m) conditions. Under salt stress conditions significant reduction was observed for plant height (PH), bolls/plant (NBP), boll weight (BW), seed cotton yield (SCY), lint% (LP), fiber length (FL), fiber strength (FS), potassium to sodium ratio (K+/Na+), potassium contents (K+), total soluble proteins (TSP), carotenoids (Car) and chlorophyll contents. Furthermore, the mean values for hydrogen peroxide (H2O2), sodium contents (Na+), catalase (CAT), superoxide dismutase (SOD), peroxidase (POD), and fiber fineness (FF) were increased under salt stress. Moderate to high heritability and genetic advancement was observed for NBP, BW, LP, SCY, K+/Na+, SOD, CAT, POD, Car, TSP, FL, and FS. Mean performance and multivariate analysis of 24 cotton genotypes based on various agro-physiological and biochemical parameters suggested that the genotypes FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, along with crosses FBS-FALCON × JSQ-White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 achieved the maximum values for K+/Na+, K+, TSP, POD, Chlb, CAT, Car, LP, FS, FL, PH, NBP, BW, and SCY under salt stress and declared as salt resistant genotypes. The above-mentioned genotypes also showed relatively higher expression levels of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 at 15 dS/m and proved the role of these ERF genes in salt tolerance in cotton. These findings suggest that these genotypes have the potential for the development of salt-tolerant cotton varieties with desirable fiber quality traits.

Molecular characterization and expression pattern of Rubisco activase gene GhRCAβ2 in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Rubisco activase (RCA) is a pivotal enzyme that can catalyse the activation of Rubisco in carbon assimilation pathway. Many studies have shown that RCA may be a potential target for genetic manipulation aimed at enhancing photosynthetic efficiency and crop yield.

OBJECTIVE

To understand the biological function of the GhRCAβ2 gene in upland cotton, we cloned the coding sequence (CDS) of the GhRCAβ2 gene and investigated its sequence features, evolutionary relationship, subcellular localization, promoter sequence and expression pattern.

METHODS

The bioinformatics tools were used to analyze the sequence features of GhRCAβ2 protein. Transient transformation of Arabidopsis mesophyll protoplasts was performed to determine the subcellular localization of the GhRCAβ2 protein. The expression pattern of the GhRCAβ2 gene was examined by analyzing transcriptome data and using the quantitative real-time PCR (qRT-PCR).

RESULTS

The full-length CDS of GhRCAβ2 was 1317 bp, and it encoded a protein with a chloroplast transit peptide. The GhRCAβ2 had two conserved ATP-binding domains, and did not have the C-terminal extension (CTE) domain that was unique to the RCA α-isoform in plants. Evolutionarily, GhRCAβ2 was clustered in Group A, and had a close evolutionary relationship with the soybean RCA. Western blot analysis demonstrated that GhRCAβ2 was immunoreactive to the RCA antibody displaying a molecular weight similar to that of the RCA β-isoform. The GhRCAβ2 protein was found in chloroplast, aligning with its role as a vital enzyme in the process of photosynthesis. The GhRCAβ2 gene had a leaf tissue-specific expression pattern, and the yellow-green leaf mutant exhibited a decreased expression of GhRCAβ2 in comparison to the wild-type cotton plants. The GhRCAβ2 promoter contained several cis-acting elements that respond to light, phytohormones and stress, suggesting that the expression of GhRCAβ2 may be regulated by these factors. An additional examination of stress response indicated that GhRCAβ2 expression was influenced by cold, heat, salt, and drought stress. Notably, diverse expression pattern was observed across different stress conditions. Additionally, low phosphorus and low potassium stress may result in a notable reduction in the expression of GhRCAβ2 gene.

CONCLUSION

Our findings will establish a basis for further understanding the function of the GhRCAβ2 gene, as well as providing valuable genetic knowledge to improve cotton photosynthetic efficiency and yield under challenging environmental circumstances.

GhCOL2 Positively Regulates Flowering by Activating the Transcription of GhHD3A in Upland Cotton (Gossypium hirsutum L.)

Plants have evolved sophisticated signaling networks to adjust flowering time, ensuring successful reproduction. Two crucial flowering regulators, FLOWERING LOCUS T (FT) and CONSTANS (CO), play pivotal roles in regulating flowering across various species. Previous studies have indicated that suppressing Gossypium hirsutum CONSTANS-LIKE 2 (GhCOL2), a homolog of Arabidopsis CO, leads to delayed flowering in cultivated cotton. However, the underlying regulatory mechanisms remain unknown. In this study, a yeast one-hybrid and dual-LUC expression assays were used to elucidate the molecular mechanism through which GhCOL2 regulates the transcription of GhHD3A. RT-qPCR was used to examine the expression of GhCOL2 and GhHD3A. Our findings reveal that GhCOL2 directly binds to CCACA cis-elements and atypical CORE (TGTGTATG) cis-elements in the promoter regions of HEADING DATE 3 A (HD3A), thereby activating GhHD3A transcription. Notably, GhCOL2 and GhHD3A exhibited high expression levels in the adult stage and low levels in the juvenile stage. Interestingly, the expression of GhCOL2 and GhHD3A varied significant between the two cotton varieties (Tx2094 and Maxxa). In summary, our study enhances the understanding of the molecular mechanism by which cotton GhCOL2-GhHD3A regulates flowering at the molecular level. Furthermore, it contributes to a broader comprehension of the GhCOL2-GhHD3A model in G. hirsutum.

Cotton host resistance as a tool for managing Rotylenchulus reniformis in Louisiana

The reniform nematode, Rotylenchulus reniformis, is a major yield-limiting pest of upland cotton (Gossypium hirsutum) in the United States that has been steadily increasing in incidence in many states. Reniform nematode-resistant cotton cultivars have recently become commercially available for cotton producers; however, few field trials have evaluated their efficacy as a nematode management tool. The aim of this study was to evaluate reniform nematode population development, plant growth, and seed cotton yield of reniform nematode-resistant cotton cultivars in two nematode-infested fields in Louisiana. Replicated small-plot field trials were conducted in St. Joseph, LA (NERS field) and Winnsboro, LA (MRRS field) during the 2022 and 2023 growing seasons. In 2022, cultivars evaluated included: (1) DP 1646 B2XF (susceptible/tolerant), (2) DP 2141NR B3XF (resistant), (3) PHY 332 W3FE (resistant), (4) PHY 411 W3FE (resistant), and (5) PHY 443 W3FE (resistant). In 2023, an additional susceptible cotton cultivar, PHY 340 W3FE, was also included. All nematode-resistant cotton cultivars evaluated provided suppression of reniform nematode population development relative to that of the susceptible cotton cultivars, with suppression of nematode soil population densities at harvest ranging from 49 - 81% relative to DP 1646 B2XF. The resistant cultivar PHY 411 W3FE provided the most consistent suppression of reniform nematode population development, reducing reniform nematode soil population densities at harvest in both field locations and both trial years. In contrast, DP 2141NR B3XF only reduced soil population densities at harvest in the NERS field in 2023. Despite relatively consistent nematode suppression and improvements in plant vigor ratings and canopy coverage associated with the resistant cotton cultivars, a yield increase was only observed with PHY 332 W3FE and PHY 411 W3FE planted at the NERS field in 2023. Despite strong resistance to reniform nematode in the evaluated cotton cultivars, nematode soil population densities still increased during the growing season in plots planted with resistant cotton cultivars, emphasizing the need for additional management tactics to use alongside host resistance. This study indicates that new reniform nematode-resistant cotton cultivars show promising potential to reduce nematode population development during the growing season in Louisiana.

Gossypium purpurascens genome provides insight into the origin and domestication of upland cotton

INTRODUCTION

Allotetraploid upland cotton (Gossypium hirsutum L.) is native to the Mesoamerican and Caribbean regions, had been improved in the southern United States by the mid-eighteenth century, was then dispersed worldwide. However, a Hainan Island Native Cotton (HIC) has long been grown extensively on Hainan Island, China.

OBJECTIVES

Explore HIC's evolutionary relationship and genomic diversity with other tetraploid cottons, its origin and whether it was used for YAZHOUBU (Yazhou cloth, World Intangible Cultural Heritage) weaving, and the role of structural variations (SVs) in upland cotton domestication.

METHODS

We assembled a high-quality genome of one HIC plant. We performed phylogenetic analysis, divergence time estimation, principal component analysis and population differentiation estimation using cotton assemblies and/or resequencing data. SVs were detected by whole-genome comparison. A F2 population was used for linkage analysis and to study effects of SVs. Buoyancy and salt water tolerance tests for seeds were conducted.

RESULTS

We found that the HIC belongs to G. purpurascens. G. purpurascens is best classified as a primitive race of G. hirsutum. The potential for long range transoceanic dispersal of G. purpurascens seeds was proved. A set of SVs, selective sweep regions between G. hirsutum races and cultivars, and quantitative trait loci (QTLs) of eleven agronomic traits were obtained. SVs, especially large-scale SVs, were found to have important effects on cotton domestication and improvement. Of them, eight large-scale inversions strongly associated with yield and fiber quality have probably undergone artificial selection in domestication.

CONCLUSION

G. purpurascens including HIC is a primitive race of G. hirsutum, probably disperse to Hainan from Central America by floating on ocean currents, may have been partly domesticated, planted and was likely used for YAZHOUBU weaving in Hainan much earlier than the Pre-Columbian period. SV plays an important role in cotton domestication and improvement.

GhVOZ1-AVP1 module positively regulates salt tolerance in upland cotton (Gossypium hirsutum L.)

Vascular Plant One‑zinc Finger (VOZ) transcription factor can respond to a variety of abiotic stresses, however its function in cotton and the molecular mechanisms of response to salt tolerance remained unclear. In this study, we found that GhVOZ1 is highly expressed in stamen and stem of cotton under normal conditions. The expression of GhVOZ1 increased significantly after 3 h of salt treatment in three-leaf staged upland cotton. Overexpressed transgenic lines of GhVOZ1 in Arabidopsis and upland cotton were treated with salt stress and we found that GhVOZ1 could respond positively to salt stress. GhVOZ1 can regulate Arabidopsis Vacuolar Proton Pump Pyrophosphatase (H+-PPase) gene (AVP1) expression through specific binding to GCGTCTAAAGTACGC site on GhAVP1 promoter, which was examined through Dual-luciferase assay and Electrophoretic mobility shift assay (EMSA). AVP1 expression was significantly increased in Arabidopsis with GhVOZ1 overexpression, while GhAVP1 expression was decreased in virus induced gene silenced (VIGS) cotton plants of GhVOZ1. Knockdown of GhAVP1 expression in cotton plants by VIGS showed decreased superoxide dismutase (SOD) and peroxidase (POD) activities, whereas an increased malondialdehyde (MDA) content and ultimately decreased salt tolerance. The GhVOZ1-AVP1 module could maintain sodium ion homeostasis through cell ion transport and positively regulate the salt tolerance in cotton, providing new ideas and insights for the study of salt tolerance.

Genomic and epigenomic insights into the mechanism of cold response in upland cotton (Gossypium hirsutum)

Functional genome research, including gene transcriptional and posttranslational modifications of histones, can benefit greatly from a high-quality genome assembly. Histone modification plays a significant role in modulating the responses to abiotic stress in plants. However, there are limited reports on the involvement of dynamic changes in histone modification in cold stress response in upland cotton. In this study, the genome of an elite accession, YM11, with considerable cold stress tolerance was de novo assembled, which yielded a genome of 2343.06 Mb with a contig N50 of 88.96 Mb, and a total of 73,821 protein-coding gene models were annotated. Comparisons among YM11 and five Gossypium allopolyploid cotton assemblies highlighted a large amount of structural variations and presence/absence variations. We analyzed transcriptome and metabolome changes in YM11 seedlings subjected to cold stress. Using the CUT&Tag method, genome-wide H3K4me3 and H3K9ac modification patterns and effect of histone changes on gene expression were profiled during cold stress. Significant and consistently changing histone modifications and the gene expressions were screened, of which transcription factors (TFs) were highlighted. Our results suggest a positive correlation between the changes in H3K4me3, H3K9ac modifications and cold stress-responsive gene activation. This genome assembly and comprehensive analysis of genome-wide histone modifications and gene expression provide insights into the genomic variation and epigenetic responses to cold stress in upland cotton.

Development and validation of allele-specific PCR-based SNP typing in a gene on chromosome D03 conferring resistance to Fusarium wilt race 4 in Upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum) is the most important fiber crop for the global textile industry. Fusarium oxysporum f. sp. vasinfectum (FOV) is one of the most destructive soil-borne fungal pathogens in cotton. Among eight pathogenic races and other strains, FOV race 4 (FOV4) is the most virulent race in US cotton production. A single nucleotide polymorphism (SNP) in a glutamate receptor-like gene (GhGLR4.8) on chromosome D03 was previously identified and validated to confer resistance to FOV race 7, and targeted genome sequencing demonstrated that it was also associated with resistance to FOV4. The objective of this study was to develop an easy and convenient PCR-based marker assay. To target the resistance SNP, a forward primer for the SNP with a mismatch in the 3rd position was designed for both the resistance (R) and susceptibility (S) alleles, respectively, with addition of 20-mer T7 promoter primer to the 5' end of the forward primer for the R allele. The two forward primers, in combination with each of five common reverse primers, were targeted to amplify amplicons of 50-260 bp in size with R and S alleles differing in 20 bp. Results showed that each of three common reverse primers in combination with the two forward primers produced polymorphic markers between R and S plants that were consistent with the targeted genome sequencing results. The polymorphism was distinctly resolved using both polyacrylamide and agarose gel electrophoreses. In addition, a sequence comparative analysis between the resistance gene and homologous sequences in sequenced tetraploid and diploid A and D genome species showed that none of the species possessed the resistance gene allele, suggesting its recent origin from a natural point mutation. The allele-specific PCR-based SNP typing method based on a three-primer combination provides a fast and convenient marker-assisted selection method to search and select for FOV4-resistant Upland cotton.

Quantitative trait loci and candidate genes for yield-related traits of upland cotton revealed by genome-wide association analysis under drought conditions

BACKGROUND

Due to the influence of extreme weather, the environment in China's main cotton-producing areas is prone to drought stress conditions, which affect the growth and development of cotton and lead to a decrease in cotton yield.

RESULTS

In this study, 188 upland cotton germplasm resources were phenotyped for data of 8 traits (including 3 major yield traits) under drought conditions in three environments for two consecutive years. Correlation analysis revealed significant positive correlations between the three yield traits. Genetic analysis showed that the estimated heritability of the seed cotton index (SC) under drought conditions was the highest (80.81%), followed by that of boll weight (BW) (80.64%) and the lint cotton index (LC) (70.49%) With genome-wide association study (GWAS) analysis, a total of 75 quantitative trait loci (QTLs) were identified, including two highly credible new QTL hotspots. Three candidate genes (Gh_D09G064400, Gh_D10G261000 and Gh_D10G254000) located in the two new QTL hotspots, QTL51 and QTL55, were highly expressed in the early stage of fiber development and showed significant correlations with SC, LC and BW. The expression of three candidate genes in two extreme materials after drought stress was analyzed by qRT-PCR, and the expression of these two materials in fibers at 15, 20 and 25 DPA. The expression of these three candidate genes was significantly upregulated after drought stress and was significantly higher in drought-tolerant materials than in drought-sensitive materials. In addition, the expression levels of the three candidate genes were higher in the early stage of fiber development (15 DPA), and the expression levels in drought-tolerant germplasm were higher than those in drought-sensitive germplasm. These three candidate genes may play an important role in determining cotton yield under drought conditions.

CONCLUSIONS

This study is helpful for understanding the regulatory genes affecting cotton yield under drought conditions and provides germplasm and candidate gene resources for breeding high-yield cotton varieties under these conditions.

Targeted development of diagnostic SNP markers for resistance to Fusarium wilt race 4 in Upland cotton (Gossypium hirsutum)

Fusarium wilt caused by the soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) has become one of the most important emerging diseases in US cotton production. Numerous QTLs have been reported for resistance to FOV; however, no major FOV4-resistance QTL or gene has been identified and used in breeding Upland cotton (Gossypium hirsutum) for FOV4 resistance. In this study, a panel of 223 Chinese Upland cotton accessions was evaluated for FOV4 resistance based on seedling mortality rate (MR) and stem and root vascular discoloration (SVD and RVD). SNP markers were developed based on targeted genome sequencing using AgriPlex Genomics. The chromosome region at 2.130-2.292 Mb on D03 was significantly correlated with both SVD and RVD but not with MR. Based on the two most significant SNP markers, accessions homozygous for AA or TT SNP genotype averaged significantly lower SVD (0.88 vs. 2.54) and RVD (1.46 vs. 3.02) than those homozygous for CC or GG SNP genotype. The results suggested that a gene or genes within the region conferred resistance to vascular discoloration caused by FOV4. The Chinese Upland accessions had 37.22% homozygous AA or TT SNP genotype and 11.66% heterozygous AC or TG SNP genotype, while 32 US elite public breeding lines all had the CC or GG SNP genotype. Among 463 obsolete US Upland accessions, only 0.86% possessed the AA or TT SNP genotype. This study, for the first time, has developed diagnostic SNPs for marker-assisted selection and identified FOV4-resistant Upland germplasms with the SNPs.

Analysis of PAT1 subfamily members in the GRAS family of upland cotton and functional characterization of GhSCL13-2A in Verticillium dahliae resistance

GhSCL13-2A, a member of the PAT1 subfamily in the GRAS family, positively regulates cotton resistance to Verticillium dahliae by mediating the jasmonic acid and salicylic acid signaling pathways and accumulation of reactive oxygen species. Verticillium wilt (VW) is a devastating disease of upland cotton (Gossypium hirsutum) that is primarily caused by the soil-borne fungus Verticillium dahliae. Scarecrow-like (SCL) proteins are known to be involved in plant abiotic and biotic stress responses, but their roles in cotton defense responses are still unclear. In this study, a total of 25 GhPAT1 subfamily members in the GRAS family were identified in upland cotton. Gene organization and protein domain analysis showed that GhPAT1 members were highly conserved. GhPAT1 genes were widely expressed in various tissues and at multiple developmental stages, and they were responsive to jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) signals. Furthermore, GhSCL13-2A was induced by V. dahliae infection. V. dahliae resistance was enhanced in Arabidopsis thaliana by ectopic overexpression of GhSCL13-2A, whereas cotton GhSCL13-2A knockdowns showed increased susceptibility. Levels of reactive oxygen species (ROS) and JA were also increased and SA content was decreased in GhSCL13-2A knockdowns. At the gene expression level, PR genes and SA signaling marker genes were down-regulated and JA signaling marker genes were upregulated in GhSCL13-2A knockdowns. GhSCL13-2A was shown to be localized to the cell membrane and the nucleus. Yeast two-hybrid and luciferase complementation assays indicated that GhSCL13-2A interacted with GhERF5. In Arabidopsis, V. dahliae resistance was enhanced by GhERF5 overexpression; in cotton, resistance was reduced in GhERF5 knockdowns. This study revealed a positive role of GhSCL13-2A in V. dahliae resistance, establishing it as a strong candidate gene for future breeding of V. dahliae-resistant cotton cultivars.

Genome-wide identification and characterization of polycomb repressive complex 2 core components in upland cotton (Gossypium hirsutum L.)

BACKGROUND

The evolutionarily conserved Polycomb Repressive Complex 2 (PRC2) plays a vital role in epigenetic gene repression by depositing tri-methylation on lysine residue K27 of histone H3 (H3K27me3) at the target loci, thus participating in diverse biological processes. However, few reports about PRC2 are available in plant species with large and complicated genomes, like cotton.

RESULTS

Here, we performed a genome-wide identification and comprehensive analysis of cotton PRC2 core components, especially in upland cotton (Gossypium hirsutum). Firstly, a total of 8 and 16 PRC2 core components were identified in diploid and tetraploid cotton species, respectively. These components were classified into four groups, E(z), Su(z)12, ESC and p55, and the members in the same group displayed good collinearity, similar gene structure and domain organization. Next, we cloned G. hirsutum PRC2 (GhPRC2) core components, and found that most of GhPRC2 proteins were localized in the nucleus, and interacted with each other to form multi-subunit complexes. Moreover, we analyzed the expression profile of GhPRC2 genes. The transcriptome data and quantitative real-time PCR (qRT-PCR) assays indicated that GhPRC2 genes were ubiquitously but differentially expressed in various tissues, with high expression levels in reproductive organs like petals, stamens and pistils. And the expressions of several GhPRC2 genes, especially E(z) group genes, were responsive to various abiotic and biotic stresses, including drought, salinity, extreme temperature, and Verticillium dahliae (Vd) infection.

CONCLUSION

We identified PRC2 core components in upland cotton, and systematically investigated their classifications, phylogenetic and synteny relationships, gene structures, domain organizations, subcellular localizations, protein interactions, tissue-specific and stresses-responsive expression patterns. Our results will provide insights into the evolution and composition of cotton PRC2, and lay the foundation for further investigation of their biological functions and regulatory mechanisms.

Combined transcriptome and metabolome analysis revealed pathways involved in improved salt tolerance of Gossypium hirsutum L. seedlings in response to exogenous melatonin application

BACKGROUND

Salinization is major abiotic stress limiting cotton production. Melatonin (MT) has been implicated in salt stress tolerance in multiple crops including upland cotton. Here, we explored the transcriptomic and metabolomic response of a salt-tolerant self-bred high-yielding cotton line SDS-01, which was exogenously sprayed with four MT concentrations (50, 100, 200, and 500 μM).

RESULTS

Here we found that MT improves plant biomass and growth under salt stress. The combined transcriptome sequencing and metabolome profiling approach revealed that photosynthetic efficiency is improved by increasing the expressions of chlorophyll metabolism and antenna proteins in MT-treated seedlings. Additionally, linoleic acid and flavonoid biosynthesis were improved after MT treatment. The Na+/K+ homeostasis-related genes were increasingly expressed in salt-stressed seedlings treated with MT as compared to the ones experiencing only salt stress. Melatonin treatment activated a cascade of plant-hormone signal transduction and reactive oxygen scavenging genes to alleviate the detrimental effects of salt stress. The global metabolome profile revealed an increased accumulation of flavonoids, organic acids, amino acids and derivatives, saccharides, and phenolic acids in MT-treated seedlings. Interestingly, N, N'-Diferuloylputrescine a known antioxidative compound was highly accumulated after MT treatment.

CONCLUSION

Collectively, our study concludes that MT is a salt stress regulator in upland cotton and alleviates salt-stress effects by modulating the expressions of photosynthesis (and related pathways), flavonoid, ROS scavenging, hormone signaling, linoleic acid metabolism, and ion homeostasis-related genes.

Identification of superior parents with high fiber quality using molecular markers and phenotypes based on a core collection of upland cotton (Gossypium hirsutum L.)

The combination of molecular markers and phenotypes to select superior parents has become the goal of modern breeders. In this study, 491 upland cotton (Gossypium hirsutum L.) accessions were genotyped using the CottonSNP80K array and then a core collection (CC) was constructed. Superior parents with high fiber quality were identified using molecular markers and phenotypes based on the CC. The Nei diversity index, Shannon's diversity index, and polymorphism information content among chromosomes for 491 accessions ranged from 0.307 to 0.402, 0.467 to 0.587, and 0.246 to 0.316, with mean values of 0.365, 0.542, and 0.291, respectively. A CC containing 122 accessions was established and was categorized into eight clusters based on the K2P genetic distances. From the CC, 36 superior parents (including duplicates) were selected, which contained the elite alleles of markers and ranked in the top 10% of phenotypic values for each fiber quality trait. Among the 36 materials, eight were for fiber length, four were for fiber strength, nine were for fiber micronaire, five were for fiber uniformity, and ten were for fiber elongation. In particular, the nine materials, 348 (Xinluzhong34), 319 (Xinluzhong3), 325 (Xinluzhong9), 397 (L1-14), 205 (XianIII9704), 258 (9D208), 464 (DP201), 467 (DP150), and 465 (DP208), possessed the elite alleles of markers for at least two traits and could be given priority in breeding applications for a more synchronous improvement of fiber quality. The work provides an efficient method for superior parent selection and will facilitate the application of molecular design breeding to cotton fiber quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01300-0.

Identification and expression analysis of arabinogalactan protein genes in cotton reveal the function of GhAGP15 in Verticillium dahliae resistance

Verticillium wilt, primarily caused by the fungal pathogen Verticillium dahliae, is a serious disease in cotton. Arabinogalactan proteins (AGPs), a class of hydroxyproline-rich glycoproteins, have been widely implicated in plant growth and environmental adaptation. The purpose of this study is to identify and characterize AGP members in cotton plants and explore their roles in responding to environmental stressors. In total, 65 GhAGP members were identified in upland cotton (Gossypium hirsutum), along with 43, 35, and 37 AGP members that were also identified in G. barbadense, G. arboreum, and G. raimondii, respectively. According to gene structure and protein domains analysis, GhAGP genes in upland cotton are highly conserved. Meanwhile, tandem duplication events have occurred frequently throughout cotton's evolutionary history. Expression analysis showed that GhAGP genes were widely expressed during growth and development and in response to abiotic stressors. Many cis-elements related to hormonal responses and environmental stressors were detected in GhAGP promoter regions. GhAGP genes participate in responding to cold, drought, and salt stress, and were sensitive to ET signaling. Furthermore, the expression level of GhAGP15 was elevated during V. dahliae infection and resistance against V. dahliae in upland cotton was significantly weakened by silencing GhAGP15 using a virus-induced gene silencing (VIGS) approach. Our results further suggest that the function of GhAGP15 in V. dahliae resistance might be involved in regulation of the JA, SA, and reactive oxygen species (ROS) pathways. The comprehensive analysis of AGP genes in cotton performed in this study provides a basic framework for further functional research of these genes.

The MADS transcription factor GhFYF is involved in abiotic stress responses in upland cotton (Gossypium hirsutum L.)

Cotton is an important textile industry raw material crops, which plays a critical role in the development of society. MADS transcription factors (TFs) play a key role about the flowering time, flower development, and abiotic stress responses in plants, but little is known about their functions on abiotic stress in cotton. In this study, a MIKC^C subfamily gene from cotton, GhFYF (FOREVER YOUNG FLOWER), was isolated and characterized. Our data showed that GhFYF localized to the nucleus. A β-glucuronidase (GUS) activity assay revealed that the promoter of GhFYF was mainly expressed in the flower and seed of ProGhFYF::GUS transgenic A. thaliana plants. The GUS staining of flowers and seeds was deepened after drought, salt treatment, and the expression level of the GUS gene and corresponding stress genes AtERD10, AtAnnexin1 are up-regulated in the inflorescence. Overexpression GhFYF in A. thaliana could promote the seed germination and growth under different salt concentrations, and determin the proline content. Yeast two-hybrid (Y2H) assays showed that GhFYF interacted with the HAD-like protein GhGPP2, which has responds to abiotic stress. Our findings indicate that GhFYF is involved in abiotic stress responses, especially for salt stress. This work establishes a solid foundation for further functional analysis of the GhFYF gene in cotton.

vasinfectum) race 4 in a backcross inbred line population of Upland cotton

KEY MESSAGE

A backcross inbred line population of cotton was evaluated for Fusarium wilt race 4 resistance at different days after inoculation (DAI). Both constitutively expressed and developmentally regulated QTLs were detected. The soil-borne fungus Fusarium oxysporum f. sp. vasinfectum (FOV) race 4 (FOV4) causes Fusarium wilt including seedling mortality in cotton. A backcross inbred line (BIL) population of 181 lines, derived from a bi-parental cross of moderately resistant non-recurrent Hai 7124 (Gossypium barbadense) and recurrent parent CCRI 36 (G. hirsutum), was evaluated under temperature-controlled conditions for FOV4 resistance with artificial inoculations. Based on three replicated tests evaluated at 7, 14, 21, and 28 days after inoculation (DAI), only 2-5 BILs showed lower disease severity ratings (DSR) than the parents while 22-50 BILs were more susceptible, indicating transgressive segregation toward susceptibility. Although DSR were overall congruent between DAI, there were many BILs displaying different responses to FOV4 across DAI. Genetic mapping using 7709 SNP markers identified 42 unique QTLs for four evaluation parameters- disease incidence (DI), DSR, mortality rate (MR), and area under disease progress curve (AUDPC), including 26 for two or more parameters. All five QTLs for AUDPC were co-localized with QTLs for DI, DSR, and/or MR at one or two DAI, indicating the unnecessary use of AUDPC in QTL mapping for FOV4 resistance. Those common QTLs explained the significant positive associations between parameters observed. Ten common QTLs with negative or positive additive effects were detected between DAI. DAI-specific and consistent QTLs were detected between DAI in cotton for the first time, suggesting the existence of both constitutively expressed and developmentally regulated QTLs for FOV4 resistance and the importance of evaluating genetic populations for FOV4 resistance at different growth stages.

Single-molecule real-time transcript sequencing of developing cotton anthers facilitates genome annotation and fertility restoration candidate gene discovery

Heterosis refers to the superior phenotypes observed in hybrids. Cytoplasmic male sterility (CMS) system plays an important role in cotton heterosis utilization. However, the global gene expression patterns of CMS-D2 and its interaction with the restorer gene Rf1 remain unclear. Here, the full-length transcript sequencing was performed in anthers of the CMS-D2 restorer line using PacBio single-molecule real-time sequencing technology. Combining PacBio SMRT long-read isoforms and Illumina RNA-seq data, 107,066 isoforms from 44,338 loci were obtained, including 10,086 novel isoforms of novel genes and 66,419 new isoforms of known genes. Totally 56,572 alternative splicing (AS) events, 1146 lncRNAs, 61 fusion transcripts and 10,466 genes exhibited alternative polyadenylation (APA), and 60,995 novel isoforms with predicted open reading frames (ORFs) were further identified. Furthermore, the specifically expressed genes in restorer line were selected and confirmed by qRT-PCR. These findings provide a basis for upland cotton genome annotation and transcriptome research, and will help to reveal the molecular mechanism of interaction between Rf1 and CMS-D2 cytoplasm.

Mitigation of salt stress response in upland cotton (Gossypium hirsutum) by exogenous melatonin

As a pleiotropic signal molecule, melatonin is ubiquitous throughout the animal and plant kingdoms and plays important roles in the regulation of plant growth, development, and responses to environmental stresses. In this study, we quantified the endogenous melatonin levels in upland cotton (Gossypium hirsutum L.), using high-performance liquid chromatography-tandem mass spectrometry. The melatonin concentrations in root, stem, and leaf were 150.60, 37.92, and 40.58 ng g fresh weight^- 1, respectively. The effects of exogenous melatonin (1 µM) on plant growth, photosynthesis, antioxidant enzyme activity, and ion homeostasis in upland cotton seedlings exposed to 100 mM NaCl treatment were determined. Pretreatment (prior to exposure to salt stress) of seedlings with exogenous melatonin significantly alleviated plant growth inhibition by salt stress and maintained an improved photosynthetic capacity. The application of melatonin also significantly reduced the salt-induced oxidative damage, possibly through the accumulation of osmotic regulatory substances and the activation of antioxidant enzymes. We also showed that exogenous melatonin regulated the expression of stress-responsive and ion-channel genes under salinity, which could contribute to improved salt tolerance in cotton. Taken together, our study provides evidence that cotton contains endogenous melatonin, and it may have unraveled crucial evidence of the role of melatonin in cotton against salt stress.

Correlation analysis of the transcriptome and metabolome reveals the role of the flavonoid biosynthesis pathway in regulating axillary buds in upland cotton (Gossypium hirsutum L.)

Flavonoids are involved in axillary bud development in upland cotton. The phenylpropanoid and flavonoid biosynthesis pathways regulate axillary bud growth by promoting the transport of auxin in upland cotton. In cotton production, simplified cultivation and mechanical harvesting are emerging trends that depend on whether the cotton plant type meets production requirements. The axillary bud is an important index of cotton plant-type traits, and the molecular mechanism of axillary bud development in upland cotton has not yet been completely studied. Here, a combined investigation of transcriptome and metabolome analyses in G. hirsutum CCRI 117 at the fourth week (stage 1), fifth week (stage 2) and sixth week (stage 3) after seedling emergence was performed. The metabolome results showed that the total lipid, amino acid and organic acid contents in the first stalk node decreased during axillary bud development. The abundance of 71 metabolites was altered between stage 2 and stage 1, and 32 metabolites exhibited significantly altered abundance between stage 3 and stage 2. According to the correlation analysis of metabolome and transcriptome profiles, we found that phenylpropanoid and flavonoid biosynthesis pathways exhibit high enrichment degrees of both differential metabolites and differential genes in three stages. Based on the verification of hormone, soluble sugar and flavonoid detection, we propose a model for flavonoid-mediated regulation of axillary bud development in upland cotton, revealing that the decrease in secondary metabolites of phenylpropanoid and flavonoid biosynthesis is an essential factor to promote the transport of auxin and subsequently promote the growth of axillary buds. Our findings provide novel insights into the regulation of phenylpropanoid and flavonoid biosynthesis in axillary bud development and could prove useful for cultivating machine-harvested cotton varieties with low axillary buds.

Production of tocotrienols in seeds of cotton (Gossypium hirsutum L.) enhances oxidative stability and offers nutraceutical potential

Upland cotton (Gossypium hirsutum L.) is an economically important multi-purpose crop cultivated globally for fibre, seed oil and protein. Cottonseed oil also is naturally rich in vitamin E components (collectively known as tocochromanols), with α- and γ-tocopherols comprising nearly all of the vitamin E components. By contrast, cottonseeds have little or no tocotrienols, tocochromanols with a wide range of health benefits. Here, we generated transgenic cotton lines expressing the barley (Hordeum vulgare) homogentisate geranylgeranyl transferase coding sequence under the control of the Brassica napus seed-specific promoter, napin. Transgenic cottonseeds had ~twofold to threefold increases in the accumulation of total vitamin E (tocopherols + tocotrienols), with more than 60% γ-tocotrienol. Matrix assisted laser desorption ionization-mass spectrometry imaging showed that γ-tocotrienol was localized throughout the transgenic embryos. In contrast, the native tocopherols were distributed unequally in both transgenic and non-transgenic embryos. α- Tocopherol was restricted mostly to cotyledon tissues and γ-tocopherol was more enriched in the embryonic axis tissues. Production of tocotrienols in cotton embryos had no negative impact on plant performance or yield of other important seed constituents including fibre, oil and protein. Advanced generations of two transgenic events were field grown, and extracts of transgenic seeds showed increased antioxidant activity relative to extracts from non-transgenic seeds. Furthermore, refined cottonseed oil from the two transgenic events showed 30% improvement in oxidative stability relative to the non-transgenic cottonseed oil. Taken together, these materials may provide new opportunities for cottonseed co-products with enhanced vitamin E profile for improved shelf life and nutrition.

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

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

Upland cotton germplasm

Bacterial blight (BB), caused by Xanthomonas citri pv. malvacearum (Xcm), is a destructive disease to cotton production in many countries. In the U.S., Xcm race 18 is the most virulent and widespread race and can cause serious yield losses. Planting BB-resistant cotton cultivars is the most effective method of controlling this disease. In this study, 335 U.S. Upland cotton accessions were evaluated for resistance to race 18 using artificial inoculations by scratching cotyledons on an individual plant basis in a greenhouse. The analysis of variance detected significant genotypic variation in disease incidence, and 50 accessions were resistant including 38 lines with no symptoms on either cotyledons or true leaves. Many of the resistant lines were developed in the MAR (multi-adversity resistance) breeding program at Texas A&M University, whereas others were developed before race 18 was first reported in the U.S. in 1973, suggesting a broad base of resistance to race 18. A genome-wide association study (GWAS) based on 26,301 single nucleotide polymorphic (SNP) markers detected 11 quantitative trait loci (QTL) anchored by 79 SNPs, including three QTL on each of the three chromosomes A01, A05 and D02, and one QTL on each of D08 and D10. This study has identified a set of obsolete Upland germplasm with resistance to race 18 and specific chromosomal regions delineated by SNPs for resistance. The results will assist in breeding cotton for BB resistance and facilitate further genomic studies in fine mapping resistance genes to enhance the understanding of the genetic basis of BB resistance in cotton.

Assessment of heterosis based on parental genetic distance estimated with SSR and SNP markers in upland cotton (Gossypium hirsutum L.)

Background

Heterosis has been extensively utilized in different crops and made a significant contribution to global food security. Genetic distance (GD) is one of the valuable criteria for selecting parents in hybrid breeding. The objectives of this study were to estimate the GD between parents using both simple sequence repeat (SSR) markers and single nucleotide polymorphism (SNP) markers and to investigate the efficiency of the prediction of hybrid performance based on GD. The experiment comprised of four male parents, 282 female parents and 1128 F₁, derived from NCII mating scheme. The hybrids, their parents and two check cultivars were evaluated for two years. Performance of F₁, mid-parent heterosis (MPH), and best parent heterosis (BPH) were evaluated for ten agronomic and fiber quality traits, including plant height, boll weight, boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index.

Results

Heterosis was observed in all hybrids and, the traits like plant height, boll number, boll weight and lint percentage exhibited higher heterosis than the fiber quality traits. Correlations were significant between parental and F₁ performances. The F₁ performances between three hybrid sets (Elite×Elite, Exotic×Elite, and Historic×Elite) showed significant differences in eight traits, including boll number, lint percentage, fiber length, fiber strength, fiber uniformity, fiber elongation ratio, micronaire, and spinning consistent index. The correlation of the GD assessed by both SSR and SNP markers was significantly positive. The cluster analysis based on GD results estimated using SNP showed that all the female parents divided into five groups and the F₁ performance between these five groups showed significant differences in four traits, including lint percentage, micronaire, fiber strength, and fiber elongation ratio. The correlation between GD and F₁ performance, MPH and BPH were significant for lint percentage and micronaire.

Conclusions

Our results suggested that GD between parents could be helpful in heterosis prediction for certain traits. This study reveals that molecular marker analysis can serve as a basis for assigning germplasm into heterotic groups and to provide guidelines for parental selection in hybrid cotton breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07431-6.

Extensive haplotypes are associated with population differentiation and environmental adaptability in Upland cotton (Gossypium hirsutum)

Three extensive eco-haplotypes associated with population differentiation and environmental adaptability in Upland cotton were identified, with A06_85658585, A08_43734499 and A06_113104285 considered the eco-loci for environmental adaptability. Population divergence is suggested to be the primary force driving the evolution of environmental adaptability in various species. Chromosome inversion increases reproductive isolation between subspecies and accelerates population divergence to adapt to new environments. Although modern cultivated Upland cotton (Gossypium hirsutum L.) has spread worldwide, the noticeable phenotypic differences still existed among cultivars grown in different areas. In recent years, the long-distance migration of cotton cultivation areas throughout China has demanded that breeders better understand the genetic basis of environmental adaptability in Upland cotton. Here, we integrated the genotypes of 419 diverse accessions, long-term environment-associated variables (EAVs) and environment-associated traits (EATs) to evaluate subgroup differentiation and identify adaptive loci in Upland cotton. Two highly divergent genomic regions were found on chromosomes A06 and A08, which likely caused by extensive chromosome inversions. The subgroups could be geographically classified based on distinct haplotypes in the divergent regions. A genome-wide association study (GWAS) also confirmed that loci located in these regions were significantly associated with environmental adaptability in Upland cotton. Our study first revealed the cause of population divergence in Upland cotton, as well as the consequences of variation in its environmental adaptability. These findings provide new insights into the genetic basis of environmental adaptability in Upland cotton, which could accelerate the development of molecular markers for adaptation to climate change in future cotton breeding.

Genome-wide analysis of chloride channel-encoding gene family members and identification of CLC genes that respond to Cl-/salt stress in upland cotton

Chloride channels (CLCs) are kinds of anion transport protein family members that are mainly distributed in cell endomembrane systems of prokaryotic and eukaryotic organisms and mediate anion (Cl^-, as a representative) transport and homeostasis. Some CLC genes have been reported to be involved in Cl^-/salt tolerance of plants exposed to NaCl stress. Through BLAST in cotton database, a total of 22 CLCs were identified in genomes A and D in upland cotton (Gossypium hirsutum L.), and except for GhCLC6 and GhCLC17, they formed highly similar homologous genes pairs. According to the prediction in PlantCARE database, many cis-acting elements related to abiotic stress responses, including ABREs, AREs, GT-1s, G-boxes, MYBs, MYCs, etc., were found in the promoters of GhCLCs. qRT-PCR revealed that most GhCLC gene expression was upregulated in the roots and leaves of cotton seedlings under salt stress, and those of homologous GhCLC4/15, GhCLC5/16, and GhCLC7/18 displayed more obvious expression. Furthermore, according to leaf virus-induced gene silencing (VIGS) assay and compared with the salt-stressed GhCLC4/15- and GhCLC7/18-silenced cotton plants, the salt-stressed GhCLC5/16-silenced plants displayed relatively better growth with significant increases in both Cl^- content and Cl^-/NO₃^- ratio in the roots and drop of the same parameters in the leaves. These results indicate that homologous GhCLC5/16, with the highest NaCl^-induced upregulation of expression and the maximum number of MYC cis-acting elements, might be the key members contributing to cotton Cl^-/salt tolerance by regulating the transport, interaction and homeostasis of Cl^- and NO₃^-.

Temporal salt stress-induced transcriptome alterations and regulatory mechanisms revealed by PacBio long-reads RNA sequencing in Gossypium hirsutum

Background

Cotton (Gossypium hirsutum) is considered a fairly salt tolerant crop however, salinity can still cause significant economic losses by affecting the yield and deteriorating the fiber quality. We studied a salt-tolerant upland cotton cultivar under temporal salt stress to unfold the salt tolerance molecular mechanisms. Biochemical response to salt stress (400 mM) was measured at 0 h, 3 h, 12 h, 24 h and 48 h post stress intervals and single-molecule long-read sequencing technology from Pacific Biosciences (PacBio) combined with the unique molecular identifiers approach was used to identify differentially expressed genes (DEG).

Results

Antioxidant enzymes including, catalase (CAT), peroxidase (POD), superoxide dismutase (SOD) were found significantly induced under temporal salt stress, suggesting that reactive oxygen species scavenging antioxidant machinery is an essential component of salt tolerance mechanism in cotton. We identified a wealth of novel transcripts based on the PacBio long reads sequencing approach. Prolonged salt stress duration induces high number of DEGs. Significant numbers of DEGs were found under key terms related to stress pathways such as “response to oxidative stress”, “response to salt stress”, “response to water deprivation”, “cation transport”, “metal ion transport”, “superoxide dismutase”, and “reductase”. Key DEGs related to hormone (abscisic acid, ethylene and jasmonic acid) biosynthesis, ion homeostasis (CBL-interacting serine/threonine-protein kinase genes, calcium-binding proteins, potassium transporter genes, potassium channel genes, sodium/hydrogen exchanger or antiporter genes), antioxidant activity (POD, SOD, CAT, glutathione reductase), transcription factors (myeloblastosis, WRKY, Apetala 2) and cell wall modification were found highly active in response to salt stress in cotton. Expression fold change of these DEGs showed both positive and negative responses, highlighting the complex nature of salt stress tolerance mechanisms in cotton.

Conclusion

Collectively, this study provides a good insight into the regulatory mechanism under salt stress in cotton and lays the foundation for further improvement of salt stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07260-z.

GWAS reveals consistent QTL for drought and salt tolerance in a MAGIC population of 550 lines derived from intermating of 11 Upland cotton (Gossypium hirsutum) parents

Cotton is grown in arid and semi-arid regions where abiotic stresses such as drought and salt are prevalent. There is a lack of studies that simultaneously address the genetic and genomic basis of tolerance to drought and salt stress. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs) together with their 11 Upland cotton parents with a total of 473,516 polymorphic SNP markers was used to identify quantitative trait loci (QTL) for drought tolerance (DT) and salt tolerance (ST) at the seedling stage based on two replicated greenhouse tests. Transgressive segregation occurred in the MAGIC-RILs, indicating that tolerant and sensitive alleles recombined for tolerance to the abiotic stress during the intermating process for the population development. A total of 20 QTL were detected for DT including 13 and 7 QTL based on plant height (PH) and dry shoot weight (DSW), respectively; and 23 QTL were detected for ST including 12 and 11 QTL for PH and DSW, respectively. There were several chromosomes with QTL clusters for abiotic stress tolerance including four QTL on chromosome A13 and three QTL on A01 for DT, and four QTL on D08 and three QTL on A11 for ST. Nine QTL (21% of the 43 QTL) detected were in common between DT and ST, indicating a common genetic basis for DT and ST. The narrow chromosomal regions for most of the QTL detected in this study allowed identification of 53 candidate genes associated with responses to salt and drought stress and abiotic stimulus. The QTL identified for both DT and ST have significantly augmented the repertoire of QTL for abiotic stress tolerance that can be used for marker-assisted selection to develop cultivars with resilience to drought and/or salt and further genomic studies towards the identification of drought and salt tolerance genes in cotton.

GhCIPK6a increases salt tolerance in transgenic upland cotton by involving in ROS scavenging and MAPK signaling pathways

BACKGROUND

Salt stress is one of the most damaging abiotic stresses in production of Upland cotton (Gossypium hirsutum). Upland cotton is defined as a medium salt-tolerant crop. Salinity hinders root development, shoots growth, and reduces the fiber quality.

RESULTS

Our previous study verified a GhCIPK6a gene response to salt stress in G. hirsutum. The homologs of GhCIPK6a were analyzed in A₂ (G. arboreum), D₅ (G. raimondii), and AD₁ (G. hirsutum) genomes. GhCIPK6a localized to the vacuole and cell membrane. The GhCBL1-GhCIPK6a and GhCBL8-GhCIPK6a complexes localized to the nucleus and cytomembrane. Overexpression of GhCIPK6a enhanced expression levels of co-expressed genes induced by salt stress, which scavenged ROS and involved in MAPK signaling pathways verified by RNA-seq analysis. Water absorption capacity and cell membrane stability of seeds from GhCIPK6a overexpressed lines was higher than that of wild-type seeds during imbibed germination stage. The seed germination rates and seedling field emergence percentages of GhCIPK6a overexpressed lines were higher than that of control line under salt stress. Moreover, overexpressing of GhCIPK6a in cotton increased lint percentage, and fiber length uniformity under salt stress.

CONCLUSIONS

We verified the function of GhCIPK6a by transformation and RNA-seq analysis. GhCIPK6a overexpressed lines exhibited higher tolerance to abiotic stresses, which functioned by involving in ROS scavenging and MAPK pathways. Therefore, GhCIPK6a has the potential for cotton breeding to improve stress-tolerance.

An RTM-GWAS procedure reveals the QTL alleles and candidate genes for three yield-related traits in upland cotton

BACKGROUND

Cotton (Gossypium spp.) fiber yield is one of the key target traits, and improved fiber yield has always been thought of as an important objective in the breeding programs and production. Although some studies had been reported for the understanding of genetic bases for cotton yield-related traits, the detected quantitative trait loci (QTL) for the traits is still very limited. To uncover the whole-genome QTL controlling three yield-related traits in upland cotton (Gossypium hirsutum L.), phenotypic traits were investigated under four planting environments and 9244 single-nucleotide polymorphism linkage disequilibrium block (SNPLDB) markers were developed in an association panel consisting of 315 accessions.

RESULTS

A total of 53, 70 and 68 significant SNPLDB loci associated with boll number (BN), boll weight (BW) and lint percentage (LP), were respectively detected through a restricted two-stage multi-locus multi-allele genome-wide association study (RTM-GWAS) procedure in multiple environments. The haplotype/allele effects of the significant SNPLDB loci were estimated and the QTL-allele matrices were organized for offering the abbreviated genetic composition of the population. Among the significant SNPLDB loci, six of them were simultaneously identified in two or more single planting environments and were thought of as the stable SNPLDB loci. Additionally, a total of 115 genes were annotated in the nearby regions of the six stable SNPLDB loci, and 16 common potential candidate genes controlling target traits of them were predicted by two RNA-seq data. One of 16 genes (GH_D06G2161) was mainly expressed in the early ovule-development stages, and the stable SNPLDB locus (LDB_19_62926589) was mapped in its promoter region.

CONCLUSION

This study identified the QTL alleles and candidate genes that could provide important insights into the genetic basis of yield-related traits in upland cotton and might facilitate breeding cotton varieties with high yield.

Genome-wide association study reveals the genetic basis of fiber quality traits in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Fiber quality is an important economic trait of cotton, and its improvement is a major goal of cotton breeding. To better understand the genetic mechanisms responsible for fiber quality traits, we conducted a genome-wide association study to identify and mine fiber-quality-related quantitative trait loci (QTLs) and genes.

RESULTS

In total, 42 single nucleotide polymorphisms (SNPs) and 31 QTLs were identified as being significantly associated with five fiber quality traits. Twenty-five QTLs were identified in previous studies, and six novel QTLs were firstly identified in this study. In the QTL regions, 822 genes were identified and divided into four clusters based on their expression profiles. We also identified two pleiotropic SNPs. The SNP locus i52359Gb was associated with fiber elongation, strength, length and uniformity, while i11316Gh was associated with fiber strength and length. Moreover, these two SNPs were nonsynonymous and located in genes Gh_D09G2376 and Gh_D06G1908, respectively. RT-qPCR analysis revealed that these two genes were preferentially expressed at one or more stages of cotton fiber development, which was consistent with the RNA-seq data. Thus, Gh_D09G2376 and Gh_D06G1908 may be involved in fiber developmental processes.

CONCLUSIONS

The findings of this study provide insights into the genetic bases of fiber quality traits, and the identified QTLs or genes may be applicable in cotton breeding to improve fiber quality.

Mitochondrial proteomic analysis reveals that proteins relate to oxidoreductase activity play a central role in pollen fertility in cotton

Cotton (Gossypium hirsutum L.) is an important economic crop. Cytoplasm male sterility (CMS) has been used to develop hybrid system and to produce hybrid seeds in cotton, but the molecular mechanism of CMS remains unclear. Mitochondria are semi-autonomous organelles, which play an important role in the reproduction of flowering plants. Male sterility has been proved associated with mitochondrial dysfunction in plants. In present study, a new strategy of proteomic sequencing data-independent acquisition (DIA) was used to analysis protein abundance across CMS lines 2074A (cytoplasm of Gossypium harknessii, D_2-2) and 2074S (cytoplasm of G. hirsutum, AD₁), and their maintainer 2074B. Comparing with transcriptome results showed that there is little consistence between proteome and transcriptome. A total of 2095 protein species were identified in three materials, and 186 and 161 differentially proteins were detected in the comparisons of 2074A vs 2074B, and 2074S vs 2074B, respectively. Among them, 49 and 50 proteins were specific existed in anther, and mainly participated in oxidoreductase activity, carbohydrate metabolism, fatty acid metabolism, cell aging, wax or cutin deposition and signal transduction. Gh_A07G0770 and Gh_D05G1908 were specific up-regulated in sterility lines, and the other genes Gh_D08G1196, Gh_D12G1971, Gh_A11G1250, Gh_D08G0388 were down-regulated, which presented similar expression tendency verified by qRT-PCR, transcriptome and proteome results. These six genes related to lipid synthesis, response to oxidative stress and cell aging, suggested them being involved in CMS occurrence. Using virus-induced gene silencing (VIGS) system, sterility obtained demonstrated the silencing Gh_A11G1250 in maintainer 2074B led to partial anthers abortion. Gh_A11G1250 encoded a mitochondrial localization of peroxisomal-like protein, participated in response to reactive oxygen species (ROS). Twenty-two proteins interacting with Gh_A11G1250 mainly related to chlorophyll biosynthetic process, photoperiodism and flowering, which showed different expression pattern between the male sterile line 2074A and maintainer 2074B. This novel research based on mitochondrial proteomics comparison confirmed that DAPs related to oxidative stress are critical to pollen abortion. BIOLOGICAL SIGNIFICANCE: Cytoplasm male sterility (CMS) system is utilized widely for hybrid production in cotton. However, the genetic and molecular mechanisms of CMS still need to be further elucidated. Up till now, fewer comprehensive comparisons of the mitochondrial proteomes from cotton CMS line and maintainer line have been reported. In this study, we performed a novel comparison of mitochondrial protein profiles in two CMS lines and their common maintainer line. Based on our results, we found a potential protein related to oxidative stress led to the anthers abortion. These results accumulate data to interpret the molecular mechanisms of CMS in cotton.

hirsutum L.)

BACKGROUND

The axillary bud is an important index of cotton plant-type traits, and the molecular mechanism of axillary bud development in upland cotton has not yet been reported. We obtained a mutant (designated mZ571) with a high-budding phenotype in axillary bud development from the low-budding phenotype variety G. hirsutum Z571 (CCRI 9A02), which provided ideal materials for the study of complex regulatory networks of axillary bud development. In this study, RNA sequencing was carried out to detect gene expression levels during three stages of axillary buds in Z571 (LB, low budding) and mZ571 mutant (HB, high budding).

RESULTS

A total of 7162 DEGs were identified in the three groups (HB-E vs. LB-E, HB-G1 vs. LB-G1, HB-G2 vs. LB-G2), including 4014 downregulated and 3184 upregulated DEGs. Additionally, 221 DEGs were commonly identified in all three groups, accounting for approximately 3.09% of the total DEGs. These DEGs were identified, annotated and classified. A significant number of DEGs were related to hormone metabolism, hormone signal transduction, and starch and sucrose metabolism. In addition, 45, 22 and 9 DEGs involved in hormone metabolic pathways and 67, 22 and 19 DEGs involved in hormone signal transduction pathwayspathway were identified in HB-E vs. LB-E, HB-G1 vs. LB-G1, and HB-G2 vs. LB-G2, respectively, suggesting that endogenous hormones are the primary factors influencing cotton axillary bud growth. Hormone and soluble sugar content measurements revealed that mZ571 exhibited higher concentrations of zeatin, gibberellins and soluble sugar in all three stages, which confirmed that these hormone metabolism-, hormone signal transduction- and starch metabolism-related genes showed interaction effects contributing to the divergence of axillary bud growth between mZ571 and Z571.

CONCLUSIONS

Our results confirmed the importance of endogenous hormones and sugars in the development of axillary buds, and we found that mZ571 plants, with a high-budding phenotype of axillary buds, exhibited higher endogenous hormone and sugar concentrations. Overall, we present a model for the emergence and development of cotton axillary buds that provides insights into the complexity and dynamic nature of the regulatory network during axillary bud emergence and development.

Comprehensive analysis of the Gossypium hirsutum L. respiratory burst oxidase homolog (Ghrboh) gene family

BACKGROUND

Plant NADPH oxidase (NOX), also known as respiratory burst oxidase homolog (rboh), encoded by the rboh gene, is a key enzyme in the reactive oxygen species (ROS) metabolic network. It catalyzes the formation of the superoxide anion (O2•-), a type of ROS. In recent years, various studies had shown that members of the plant rboh gene family were involved in plant growth and developmental processes as well as in biotic and abiotic stress responses, but little is known about its functional role in upland cotton.

RESULTS

In the present study, 26 putative Ghrboh genes were identified and characterized. They were phylogenetically classified into six subfamilies and distributed at different densities across 18 of the 26 chromosomes or scaffolds. Their exon-intron structures, conserved domains, synteny and collinearity, gene family evolution, regulation mediated by cis-acting elements and microRNAs (miRNAs) were predicted and analyzed. Additionally, expression profiles of Ghrboh gene family were analyzed in different tissues/organs and at different developmental stages and under different abiotic stresses, using RNA-Seq data and real-time PCR. These profiling studies indicated that the Ghrboh genes exhibited temporal and spatial specificity with respect to expression, and might play important roles in cotton development and in stress tolerance through modulating NOX-dependent ROS induction and other signaling pathways.

CONCLUSIONS

This comprehensive analysis of the characteristics of the Ghrboh gene family determined features such as sequence, synteny and collinearity, phylogenetic and evolutionary relationship, expression patterns, and cis-element- and miRNA-mediated regulation of gene expression. Our results will provide valuable information to help with further gene cloning, evolutionary analysis, and biological function analysis of cotton rbohs.

Introgression Leads to Genomic Divergence and Responsible for Important Traits in Upland Cotton

Understanding the genetic diversity and population structure of germplasms is essential when selecting parents for crop breeding. The genomic changes that occurred during the domestication and improvement of Upland cotton (Gossypium hirsutum) remains poorly understood. Besides, the available genetic resources from cotton cultivars are limited. By applying restriction site-associated DNA marker sequencing (RAD-seq) technology to 582 tetraploid cotton accessions, we confirmed distinct genomic regions on chromosomes A06 and A08 in Upland cotton cultivar subgroups. Based on the pedigree, reported QTLs, introgression analyses, and genome-wide association study (GWAS), we suggest that these divergent regions might have resulted from the introgression of exotic lineages of G. hirsutum landraces and their wild relatives. These regions were the typical genomic signatures that might be responsible for maturity and fiber quality on chromosome A06 and chromosome A08, respectively. Moreover, these genomic regions are located in the putative pericentromeric regions, implying that their application will be challenging. In the study, based on high-density SNP markers, we reported two genomic signatures on chromosomes A06 and A08, which might originate from the introgression events in the Upland cotton population. Our study provides new insights for understanding the impact of historic introgressions on population divergence and important agronomic traits of modern Upland cotton cultivars.

Dissecting the genetic architecture of seed-cotton and lint yields in Upland cotton using genome-wide association mapping

Seed-cotton yield (SY) and lint yield (LY) are the most important yield traits of cotton. Thus, it is critical to dissect their genetic architecture. Upland cotton (Gossypium hirsutum) is widely grown worldwide. In this study, a genome-wide association mapping was performed based on the CottonSNP80K array to dissect the genetic architecture of SY and LY in Upland cotton. Twenty-three significant associations were detected within four environments, including 11 associated with SY and 12 associated with LY. Seven single nucleotide polymorphisms (SNPs), TM234, TM237, TM247, TM255, TM256, TM263, and TM264, were co-associated with the two traits, which may indicate pleiotropy or intergenic tight linkages. Five SNPs, TM13332, TM39771, TM57119, TM81653, and TM81660, were coincided with those of previous reports and could be used in marker-assisted selection. Combining functional annotations with expression analyses of the genes identified within 400 kb of the significantly associated SNPs, we hypothesize that the three genes, Gh_D05G1077 and Gh_D13G1571 for SY, and Gh_A11G0775 for LY, may have the potential to increase cotton yield. The results would provide useful information for understanding the genetic basis of yield traits in Upland cotton and for facilitating its high-yield breeding through molecular design.

High-density linkage map construction and QTL analyses for fiber quality, yield and morphological traits using CottonSNP63K array in upland cotton (Gossypium hirsutum L.)

Background

Improving fiber quality and yield are the primary research objectives in cotton breeding for enhancing the economic viability and sustainability of Upland cotton production. Identifying the quantitative trait loci (QTL) for fiber quality and yield traits using the high-density SNP-based genetic maps allows for bridging genomics with cotton breeding through marker assisted and genomic selection. In this study, a recombinant inbred line (RIL) population, derived from cross between two parental accessions, which represent broad allele diversity in Upland cotton, was used to construct high-density SNP-based linkage maps and to map the QTLs controlling important cotton traits.

Results

Molecular genetic mapping using RIL population produced a genetic map of 3129 SNPs, mapped at a density of 1.41 cM. Genetic maps of the individual chromosomes showed good collinearity with the sequence based physical map. A total of 106 QTLs were identified which included 59 QTLs for six fiber quality traits, 38 QTLs for four yield traits and 9 QTLs for two morphological traits. Sub-genome wide, 57 QTLs were mapped in A sub-genome and 49 were mapped in D sub-genome. More than 75% of the QTLs with favorable alleles were contributed by the parental accession NC05AZ06. Forty-six mapped QTLs each explained more than 10% of the phenotypic variation. Further, we identified 21 QTL clusters where 12 QTL clusters were mapped in the A sub-genome and 9 were mapped in the D sub-genome. Candidate gene analyses of the 11 stable QTL harboring genomic regions identified 19 putative genes which had functional role in cotton fiber development.

Conclusion

We constructed a high-density genetic map of SNPs in Upland cotton. Collinearity between genetic and physical maps indicated no major structural changes in the genetic mapping populations. Most traits showed high broad-sense heritability. One hundred and six QTLs were identified for the fiber quality, yield and morphological traits. Majority of the QTLs with favorable alleles were contributed by improved parental accession. More than 70% of the mapped QTLs shared the similar map position with previously reported QTLs which suggest the genetic relatedness of Upland cotton germplasm. Identification of QTL clusters could explain the correlation among some fiber quality traits in cotton. Stable and major QTLs and QTL clusters of traits identified in the current study could be the targets for map-based cloning and marker assisted selection (MAS) in cotton breeding. The genomic region on D12 containing the major stable QTLs for micronaire, fiber strength and lint percentage could be potential targets for MAS and gene cloning of fiber quality traits in cotton.

A high density SLAF-seq SNP genetic map and QTL for seed size, oil and protein content in upland cotton

BACKGROUND

Cotton is a leading natural fiber crop. Beyond its fiber, cottonseed is a valuable source of plant protein and oil. Due to the much higher value of cotton fiber, there is less consideration of cottonseed quality despite its potential value. Though some QTL controlling cottonseed quality have been identified, few of them that warrant further study are known. Identifying stable QTL controlling seed size, oil and protein content is necessary for improvement of cottonseed quality.

RESULTS

In this study, a recombinant inbred line (RIL) population was developed from a cross between upland cotton cultivars/lines Yumian 1 and M11. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to construct a genetic map that covered 3353.15 cM with an average distance between consecutive markers of 0.48 cM. The seed index, together with kernel size, oil and protein content were further used to identify QTL. In total, 58 QTL associated with six traits were detected, including 13 stable QTL detected in all three environments and 11 in two environments.

CONCLUSION

A high resolution genetic map including 7033 SNP loci was constructed through specific locus amplified fragment sequencing technology. A total of 13 stable QTL associated with six cottonseed quality traits were detected. These stable QTL have the potential for fine mapping, identifying candidate genes, elaborating molecular mechanisms of cottonseed development, and application in cotton breeding programs.

Data describing Upland cotton cultivars and advanced breeding lines used in Colombia

In the last century, more than a hundred of cultivars were used in the cotton production system in Colombia. Breeding for cultivars adapted to tropical environments had been the main purpose of the Colombian agricultural research institutions dedicated to cotton. Data describing yield and fiber quality traits of these cultivars (and the introduced ones mainly from USA) is scattered across grey literature which reduces chances of discovering, accessing and assessing this information.This data article contains databases describing i) Colombian and introduced Upland cotton cultivars used in Colombia and ii) ramulosis-resistance scores of lines developed by the Colombian breeding program. The first database was constructed from data extracted from grey literature mainly produced by ICA and CORPOICA (rebranded today as AGROSAVIA), the Colombian agricultural research agencies. The second one describes the Cereté lines (LCER) database. These advanced breeding lines were developed for improved yield performance in tropical environments, specifically monsoon and savanna climates. The LCER dataset also describes the ramulosis field resistance of these cultivars. Ramulosis is an endemic disease in South America caused by Colletotrichum gossypii var. cephalosporioides. The data in this article supports and augments information presented in the research articles [1]: and [3].

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.

Stress response of NAD+-dependent formate dehydrogenase in Gossypium hirsutum L. grown under copper toxicity

Cotton (Gossypium hirsutum L.), which is not directly involved in the food chain, appears to be a suitable candidate to remove heavy metals from the food chain and to be a commercial plant which could be planted in contaminated soils. The key point of this approach is selection of the right genotype, which has heavy metal resistance or hyperaccumulation properties. Therefore, in the present study, two G. hirsutum genotypes, Erşan-92 and N-84S, were grown under copper stress and investigated to obtain further insights about the heavy metal tolerance mechanisms of plants by focusing on the expression of NAD⁺-dependent formate dehydrogenase (FDH). In accordance with the results, which were obtained from RT-PCR analysis and activity measurements, in the Erşan-92 root tissue, FDH activity increased significantly with increasing metal concentrations and a 6.35-fold higher FDH activity was observed in the presence of 100-μM Cu. As opposed to Erşan-92, the maximum FDH activity in the roots of N-84S, which were untreated with copper as the control plants, was measured as 0.0141-U mg^-1 g^-1 FW, and the activity decreased significantly with the increasing metal concentrations. The metallothionein (GhMT3a) transcript level of the plants grown in a medium containing different Cu concentrations showed nearly the same pattern as that of the FDH gene transcription. It was observed that while the tolerance of N-84S in the lower Cu concentration reduces remarkably, Erşan-92 continues to struggle up to 100-μM Cu. The results of the SOD analysis also confirm this activity of Erşan-92 against the Cu stress.

A combined small RNA and transcriptome sequencing analysis reveal regulatory roles of miRNAs during anther development of Upland cotton carrying cytoplasmic male sterile Gossypium harknessii (D2) cytoplasm

BACKGROUND

Cytoplasmic male sterility (CMS) in flowering plants is usually caused by incompatibility between mitochondrial and nuclear genomes, and can be restored by nuclear genes known as restorer-of-fertility (Rf). Although the CMS/Rf system is useful and convenient for economic production of commercial hybrid seed, the molecular mechanisms of CMS occurrence and fertility restoration in cotton are unclear.

RESULTS

Here, a combined small RNA and transcriptome sequencing analysis was performed on floral buds at the meiosis stage in three-line hybrid cotton system, and differentially expressed microRNAs (DEMs) and their target genes were identified and further analyzed for a possible involvement in CMS and fertility restoration. Totally 10 and 30 differentially expressed miRNA-target gene pairs were identified in A-B and A-R comparison group, respectively. A putative regulatory network of CMS occurrence and fertility restoration-related miRNA-target pairs during anther development were then constructed. The RLM-RACE analysis showed that gra-miR7505b regulates a PPR gene (Gh_D05G3392) by cleaving precisely at the 643 nt and 748 nt sites. The further analysis indicated that the sequence variation in the binding regions of Gh_D05G3392 and Gh_D05G3356 may cause a lower cleavage efficiency of the PPR genes by miR7505b and miR7505 in R line, respectively, leading to the up-regulation of the PPR genes and fertility restoration. These results have established their genetic involvement in fertility restoration in the CMS-D2 system.

CONCLUSION

Our combined miRNA and transcriptome analysis in three-line hybrid cotton system provides new insights into the molecular mechanisms of CMS occurrence and fertility restoration, which will contribute to further hybrid breeding in cotton.

Simple sequence repeat markers associated/linked with agronomic traits, as core primers, are eminently suitable for DNA fingerprinting in Upland cotton

Analyzing the genetic differences among crop germplasm resources scientifically and accurately is very important for the selection of core accessions, the identification of new cultivars, and the determination of seed purity. However, phenotypic selection per se is not sufficient to identify genetically distinct accessions. In this study, 26 out of 83 simple sequence repeat markers associated/linked with cotton important agronomic traits derived from our previous and other published research, corresponding to the 26 chromosomes of Upland cotton (Gossypium hirsutum L.), were selected as core primers for DNA fingerprinting construction. The 26 markers showed clear band patterns, good repeatability and high polymorphism. The average alleles, gene diversity index and polymorphism information content were 3.12, 0.4312 and 0.3830, respectively. Using TM-1, a genetic standard line for Upland cotton, as the control, DNA fingerprinting pattern and DNA barcodes were obtained based on the core primers. There was a significant positive correlation between genetic distance matrix determined using 26 core primers and that determined using more primers (335) derived from previous research, further suggesting that the core primers were eminently suitable for DNA fingerprinting in Upland cotton. This study provides a molecular basis for assessing identification, authenticity and seed purity of cotton cultivars.

RNA interference of GhPEPC2 enhanced seed oil accumulation and salt tolerance in Upland cotton

Phosphoenolpyruvate carboxylase (PEPCase) mainly produces oxaloacetic acid for tricarboxylic acid (TCA) cycle. Here we reported that GhPEPC2 silencing with PEPC2-RNAi vector could regulate oil and protein accumulation in cottonseeds. In GhPEPC2 transgenic plants, PEPCase activities in immature embryos were significantly reduced, and the oil content in seed kernel was increased 7.3 percentages, whereas total proteins decreased 5.65 percentages. Compared to wild type, agronomical traits of transgenic plant were obviously unaffected. Furthermore, gene expression profile of GhPEPC2 transgenic seeds were investigated using RNA-seq, most lipid synthesis related genes were up-regulated, but amino acid metabolic related genes were down-regulated. In addition, the GhPEPC2 transgenic cotton seedlings were stressed using sodium salts at seedling stage, and the salt tolerance was significantly enhanced. Our observations of GhPEPC2 in cotton would shade light on understanding the regulation of oil content, protein accumulation and salt tolerance enhancement in other plants.

Exploration of miRNAs and target genes of cytoplasmic male sterility line in cotton during flower bud development

Cytoplasmic male sterility (CMS) lines provide crucial material to harness heterosis for crop plants, which serves as an important strategy for hybrid seed production. However, the molecular mechanism remains obscure. Although microRNAs (miRNAs) play important roles in vegetative growth and reproductive growth, there are few reports on miRNAs regulating the development of male sterility in Upland cotton. In present study, 12 small RNA libraries were constructed and sequenced for two development stages of flower buds from a CMS line and its maintainer line. Based on the results, 256 novel miRNAs were allocated to 141 new miRNA families, and 77 known miRNAs belonging to 54 conserved miRNA families were identified as well. Comparative analysis revealed that 61 novel and 10 conserved miRNAs were differentially expressed. Further transcriptome analysis identified 232 target genes for these miRNAs, which participated in cellular developmental process, cell death, pollen germination, and sexual reproduction. In addition, expression patterns of typical miRNA and the negatively regulated target genes, such as PPR, ARF, AP2, and AFB, were verified by qRT-PCR in cotton flower buds. These targets were previously reported to be related to reproduction development and male sterility, suggesting that miRNAs might act as regulators of CMS occurrence. Some miRNAs displayed specific expression profiles in special developmental stages of CMS line and its fertile hybrid (F₁). Present study offers new information on miRNAs and their related target genes in exploiting CMS mechanism, and revealing the miRNA regulatory networks in Upland cotton.

Detection of favorable alleles for yield and yield components by association mapping in upland cotton

Association mapping based on linkage disequilibrium provides a promising tool for dissecting the genetic basis underlying complex traits. To reveal the genetic variations of yield and yield components traits in upland cotton, 403 upland cotton accessions were collected and analyzed by 560 genome-wide simple sequence repeats (SSRs). A diverse panel consisting of 403 upland cotton accessions was grown in six different environments, and the yield and yield component traits were measured, and 560 SSR markers covering the whole genome were mapped. Association studies were performed to uncover the genotypic and phenotypic variations using a mixed linear model. Favorable alleles and typical accessions for yield traits were identified. A total of 201 markers were polymorphic, revealing 394 alleles. The average gene diversity and polymorphism information content were 0.556 and 0.483, respectively. Based on a population structure analysis, 403 accessions were divided into two subgroups. A mixed linear model analysis of the association mapping detected 43 marker loci according to the best linear unbiased prediction and in at least three of the six environments(- lgP > 1.30, P < 0.05). Among the 43 associated markers, five were associated with more than two traits simultaneously and nine were coincident with those identified previously. Based on phenotypic effects, favorable alleles and typical accessions that contained the elite allele loci related to yield traits were identified and are widely used in practical breeding. This study detected favorable quantitative trait loci's alleles and typical accessions for yield traits, these are excellent genetic resources for future high-yield breeding by marker-assisted selection in upland cotton in China.

QTL Mapping for Fiber Quality and Yield Traits Based on Introgression Lines Derived from Gossypium hirsutum × G. tomentosum

The tetraploid species Gossypium hirsutum is cultivated widely throughout the world with high yield and moderate fiber quality, but its genetic basis is narrow. A set of 107 introgression lines (ILs) was developed with an interspecific cross using G. hirsutumacc. 4105 as the recurrent parent and G. tomentosum as the donor parent. A specific locus amplified fragment sequencing (SLAF-seq) strategy was used to obtain high-throughput single nucleotide polymorphism (SNP) markers. In total, 3157 high-quality SNP markers were obtained and further used for identification of quantitative trait loci (QTLs) for fiber quality and yield traits evaluated in multiple environments. In total, 74 QTLs were detected that were associated with five fiber quality traits (30 QTLs) and eight yield traits (44 QTLs), with 2.02-30.15% of the phenotypic variance explained (PVE), and 69 markers were found to be associated with these thirteen traits. Eleven chromosomes in the A sub-genome (At) harbored 47 QTLs, and nine chromosomes in the D sub-genome (Dt) harbored 27 QTLs. More than half (44 QTLs = 59.45%) showed positive additive effects for fiber and yield traits. Five QTL clusters were identified, with three in the At, comprised of thirteen QTLs, and two in the Dt comprised of seven QTLs. The ILs developed in this study and the identified QTLs will facilitate further molecular breeding for improvement of Upland cotton in terms of higher yield with enhanced fiber quality.

Overexpression of GhDof1 improved salt and cold tolerance and seed oil content in Gossypium hirsutum

A homologous GhDof1, which belongs to a large family of plant-specific transcription factor DOF, was isolated from Upland cotton (Gossypium hirsutum L.). GhDof1 protein was located in the nucleus of onion epidermal cells, the core domain of transcriptional activity existed in the C-terminal, and the activity elements of GhDof1 promoter existed in the regions of -645∼ -469bp and -286∼ -132bp of transcriptional start codon. GhDof1 constitutively expressed in leaves, roots and stems, accumulated highest in leaves. The salinity and cold treatments induced GhDof1 transcript accumulation. The GhDof1-overexpressed cotton showed significantly higher salt and cold tolerance over the wild-type plants. Under salt stress, the root growth of overexpressed GhDof1 lines was promoted. The expression levels of stress-responsive genes, GhP5CS, GhSOD and GhMYB, were differently up-regulated in transgenic lines. Oil contents increased in some transgenic plants, and protein contents reduced compared with transformed receptor. These results suggested that GhDof1 was a functional transcription factor for improving the abiotic tolerance and seed oil content in Upland cotton.

Comparative genomic analysis of the PKS genes in five species and expression analysis in upland cotton

Plant type III polyketide synthase (PKS) can catalyse the formation of a series of secondary metabolites with different structures and different biological functions; the enzyme plays an important role in plant growth, development and resistance to stress. At present, the PKS gene has been identified and studied in a variety of plants. Here, we identified 11 PKS genes from upland cotton (Gossypium hirsutum) and compared them with 41 PKS genes in Populus tremula, Vitis vinifera, Malus domestica and Arabidopsis thaliana. According to the phylogenetic tree, a total of 52 PKS genes can be divided into four subfamilies (I-IV). The analysis of gene structures and conserved motifs revealed that most of the PKS genes were composed of two exons and one intron and there are two characteristic conserved domains (Chal_sti_synt_N and Chal_sti_synt_C) of the PKS gene family. In our study of the five species, gene duplication was found in addition to Arabidopsis thaliana and we determined that purifying selection has been of great significance in maintaining the function of PKS gene family. From qRT-PCR analysis and a combination of the role of the accumulation of proanthocyanidins (PAs) in brown cotton fibers, we concluded that five PKS genes are candidate genes involved in brown cotton fiber pigment synthesis. These results are important for the further study of brown cotton PKS genes. It not only reveals the relationship between PKS gene family and pigment in brown cotton, but also creates conditions for improving the quality of brown cotton fiber.

Fine mapping and candidate gene analysis of the virescent gene v 1 in Upland cotton (Gossypium hirsutum)

The young leaves of virescent mutants are yellowish and gradually turn green as the plants reach maturity. Understanding the genetic basis of virescent mutants can aid research of the regulatory mechanisms underlying chloroplast development and chlorophyll biosynthesis, as well as contribute to the application of virescent traits in crop breeding. In this study, fine mapping was employed, and a recessive gene (v ₁) from a virescent mutant of Upland cotton was narrowed to an 84.1-Kb region containing ten candidate genes. The GhChlI gene encodes the cotton Mg-chelatase I subunit (CHLI) and was identified as the candidate gene for the virescent mutation using gene annotation. BLAST analysis showed that the GhChlI gene has two copies, Gh_A10G0282 and Gh_D10G0283. Sequence analysis indicated that the coding region (CDS) of GhChlI is 1269 bp in length, with three predicted exons and one non-synonymous nucleotide mutation (G1082A) in the third exon of Gh_D10G0283, with an amino acid (AA) substitution of arginine (R) to lysine (K). GhChlI-silenced TM-1 plants exhibited a lower GhChlI expression level, a lower chlorophyll content, and the virescent phenotype. Analysis of upstream regulatory elements and expression levels of GhChlI showed that the expression quantity of GhChlI may be normal, and with the development of the true leaf, the increase in the Gh_A10G0282 dosage may partially make up for the deficiency of Gh_D10G0283 in the v ₁ mutant. Phylogenetic analysis and sequence alignment revealed that the protein sequence encoded by the third exon of GhChlI is highly conserved across diverse plant species, in which AA substitutions among the completely conserved residues frequently result in changes in leaf color in various species. These results suggest that the mutation (G1082A) within the GhChlI gene may cause a functional defect of the GhCHLI subunit and thus the virescent phenotype in the v₁ mutant. The GhChlI mutation not only provides a tool for understanding the associations of CHLI protein function and the chlorophyll biosynthesis pathway but also has implications for cotton breeding.