Fine-mapping and candidate gene analysis of qFS-Chr. D02, a QTL for fibre strength introgressed from a semi-wild cotton into Gossypium hirsutum

An aldehyde dehydrogenase gene, GhALDH7B4_A06, positively regulates fiber strength in upland cotton (Gossypium hirsutum L.)

High fiber strength (FS) premium cotton has significant market demand. Consequently, enhancing FS is a major objective in breeding quality cotton. However, there is a notable lack of known functionally applicable genes that can be targeted for breeding. To address this issue, our study used specific length-amplified fragment sequencing combined with bulk segregant analysis to study FS trait in an F2 population. Subsequently, we integrated these results with previous quantitative trait locus mapping results regarding fiber quality, which used simple sequence repeat markers in F2, F2:3, and recombinant inbred line populations. We identified a stable quantitative trait locus qFSA06 associated with FS located on chromosome A06 (90.74-90.83 Mb). Within this interval, we cloned a gene, GhALDH7B4_A06, which harbored a critical mutation site in coding sequences that is distinct in the two parents of the tested cotton line. In the paternal parent Ji228, the gene is normal and referred to as GhALDH7B4_A06O; however, there is a nonsense mutation in the maternal parent Ji567 that results in premature termination of protein translation, and this gene is designated as truncated GhALDH7B4_A06S. Validation using recombinant inbred lines and gene expression analysis revealed that this mutation site is correlated with cotton FS. Virus-induced gene silencing of GhALDH7B4 in cotton caused significant decreases in FS and fiber micronaire. Conversely, GhALDH7B4_A06O overexpression in Arabidopsis boosted cell wall component contents in the stem. The findings of our study provide a candidate gene for improving cotton fiber quality through molecular breeding.

Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L.)

KEY MESSAGE

GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3's expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.

Genome-wide association study for boll weight in Gossypium hirsutum races

High yield has always been an essential target in almost all of the cotton breeding programs. Boll weight (BW) is a key component of cotton yield. Numerous linkage mapping and genome-wide association studies (GWAS) have been performed to understand the genetic mechanism of BW, but information on the markers/genes controlling BW remains limited. In this study, we conducted a GWAS for BW using 51,268 high-quality single-nucleotide polymorphisms (SNPs) and 189 Gossypium hirsutum accessions across five different environments. A total of 55 SNPs significantly associated with BW were detected, of which 29 and 26 were distributed in the A and D subgenomes, respectively. Five SNPs were simultaneously detected in two environments. For TM5655, TM8662, TM36371, and TM50258, the BW grouped by alleles of each SNP was significantly different. The ± 550 kb regions around these four key SNPs contained 262 genes. Of them, Gh_A02G1473, Gh_A10G1765, and Gh_A02G1442 were expressed highly at 0 to 1 days post-anthesis (dpa), - 3 to 0 dpa, and - 3 to 0 dpa in ovule of TM-1, respectively. They were presumed as the candidate genes for fiber cell differentiation, initiation, or elongation based on gene annotation of their homologs. Overall, these results supplemented valuable information for dissecting the genetic architecture of BW and might help to improve cotton yield through molecular marker-assisted selection breeding and molecular design breeding.

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture

BACKGROUND

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production.

GENOMICS APPROACHES IN COTTON

In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research.

CONCLUSIONS

This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

Genome-Wide Association Study of Lint Percentage in Gossypium hirsutum L. Races

Lint percentage is one of the most essential yield components and an important economic index for cotton planting. Improving lint percentage is an effective way to achieve high-yield in cotton breeding worldwide, especially upland cotton (Gossypium hirsutum L.). However, the genetic basis controlling lint percentage has not yet been systematically understood. Here, we performed a genome-wide association mapping for lint percentage using a natural population consisting of 189 G. hirsutum accessions (188 accessions of G. hirsutum races and one cultivar TM-1). The results showed that 274 single-nucleotide polymorphisms (SNPs) significantly associated with lint percentage were detected, and they were distributed on 24 chromosomes. Forty-five SNPs were detected at least by two models or at least in two environments, and their 5 Mb up- and downstream regions included 584 makers related to lint percentage identified in previous studies. In total, 11 out of 45 SNPs were detected at least in two environments, and their 550 Kb up- and downstream region contained 335 genes. Through RNA sequencing, gene annotation, qRT-PCR, protein-protein interaction analysis, the cis-elements of the promotor region, and related miRNA prediction, Gh_D12G0934 and Gh_A08G0526 were selected as key candidate genes for fiber initiation and elongation, respectively. These excavated SNPs and candidate genes could supplement marker and gene information for deciphering the genetic basis of lint percentage and facilitate high-yield breeding programs of G. hirsutum ultimately.

Transcriptome Revealed GhPP2C43-A Negatively Regulates Salinity Tolerance in an Introgression Line from a Semi-wild Upland Cotton

Salt damage is one of the major threats to sustainable cotton production owing to the limited arable land in China mainly occupied by the production of staple food crops. Salt-stress tolerant cotton varieties are lacking in production and, the mechanisms underpinning salt-stress tolerance in cotton remain enigmatic. Here, DM37, an intraspecific introgression line from G. hirsutum race yucatanense acc TX-1046 into the G. hirsutum acc TM-1 background, was found to be highly tolerant to salt stress. Its seed germination rate and germination potential were significantly higher than the recipient TM-1 under salt stress. Physiological analysis showed DM37 had higher proline content and Peroxidase activity, as well as lower Na+/K+ ratios at the seedling stage, consistent with higher seedling survival rate after durable salt stress. Furthermore, comparative transcriptome analysis revealed that responsive patterns to salt stress in DM37 were different from TM-1. Weighted Correlation Network Analysis (WGCNA) demonstrated that co-expression modules associated with salt stress in DM37 also differed from TM-1. Out of them, GhPP2C43-A, a phosphatase gene, exhibited negative regulation of salt-stress tolerance verified by VIGS and transgenic Arabidopsis. Gene expression showed GhPP2C43-A in TM-1 was induced by durable salt stress but not in DM37 probably attributing to the variation of cis-element in its promoter, thereby being conferred different salt-stress tolerance. Our result would provide new genes/germplasms from semi-wild cotton in salt-stress tolerant cotton breeding. This study would give us new insights into the mechanisms underpinning the salt-stress tolerance in cotton.

Fine mapping and candidate gene analysis of qFL-A12-5: a fiber length-related QTL introgressed from Gossypium barbadense into Gossypium hirsutum

The fiber length-related qFL-A12-5 identified in CSSLs introgressed from Gossypium barbadense into Gossypium hirsutum was fine-mapped to an 18.8 kb region on chromosome A12, leading to the identification of the GhTPR gene as a potential regulator of cotton fiber length. Fiber length is a key determinant of fiber quality in cotton, and it is a key target of artificial selection for breeding and domestication. Although many fiber length-related quantitative trait loci have been identified, there are few reports on their fine mapping or candidate gene validation, thus hampering efforts to understand the mechanistic basis of cotton fiber development. Our previous study identified the qFL-A12-5 associated with superior fiber quality on chromosome A12 in the chromosome segment substitution line (CSSL) MBI7747 (BC₄F_3:5). A single segment substitution line (CSSL-106) screened from BC₆F₂ was backcrossed to construct a larger segregation population with its recurrent parent CCRI45, thus enabling the fine mapping of 2852 BC₇F₂ individuals using denser simple sequence repeat markers to narrow the qFL-A12-5 to an 18.8 kb region of the genome, in which six annotated genes were identified in Gossypium hirsutum. Quantitative real-time PCR and comparative analyses led to the identification of GH_A12G2192 (GhTPR) encoding a tetratricopeptide repeat-like superfamily protein as a promising candidate gene for qFL-A12-5. A comparative analysis of the protein-coding regions of GhTPR among Hai1, MBI7747, and CCRI45 revealed two non-synonymous mutations. The overexpression of GhTPR resulted in longer roots in Arabidopsis, suggesting that GhTPR may regulate cotton fiber development. These results provide a foundation for future efforts to improve cotton fiber length.

Fine-mapping and candidate gene analysis of qFL-c10-1 controlling fiber length in upland cotton (Gossypium hirsutum L.)

A fiber length QTL, qFL-c10-1, was fine-mapped to a 96.5-kb region containing one gene that has not been characterized in plants. Fiber length is an important component of cotton fiber quality, which is associated with other quality properties such as fiber strength, fiber maturity, and fineness. In our previous studies, a stable QTL qFL-c10-1 controlling fiber length had been identified on chromosome A10 in an upland cotton recombinant inbred line (RIL) population from a cross between Jimian5 and DH962. To fine-map qFL-c10-1, an F₂ population with 1081 individual plants from a cross between a recombinant line DJ61 and Jimian5 was established. Using linkage analysis and progeny recombination experiment, qFL-c10-1 was mapped into a 96.5-kb genomic region that just contained one proper transcript Ghir_A10G022020 (described as GhFL10), an undescribed gene in plants. One 214-bp deletion was identified in the promoter region of DJ61 compared with Jimian5. Quantitative real-time PCR (qRT-PCR) and comparative analysis of parental sequences suggested that GhFL10 was the most promising candidate gene for qFL-c10-1. According to RNA-seq, yeast two-hybrid assay and bimolecular fluorescence complementation (BiFC), we speculate that GhFL10 interacts with NF-YA transcription factors to negatively regulate fiber elongation.

Analysis of transcriptome data and quantitative trait loci enables the identification of candidate genes responsible for fiber strength in Gossypium barbadense

Gossypium barbadense possesses a superior fiber quality because of its fiber length and strength. An in-depth analysis of the underlying genetic mechanism could aid in filling the gap in research regarding fiber strength and could provide helpful information for Gossypium barbadense breeding. Three quantitative trait loci related to fiber strength were identified from a Gossypium barbadense recombinant inbred line (PimaS-7 × 5917) for further analysis. RNA sequencing was performed in the fiber tissues of PimaS-7 × 5917 0–35 days postanthesis. Four specific modules closely related to the secondary wall-thickening stage were obtained using the weighted gene coexpression network analysis. In total, 55 genes were identified as differentially expressed from 4 specific modules. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for enrichment analysis, and Gbar_D11G032910, Gbar_D08G020540, Gbar_D08G013370, Gbar_D11G033670, and Gbar_D11G029020 were found to regulate fiber strength by playing a role in the composition of structural constituents of cytoskeleton and microtubules during fiber development. Quantitative real-time PCR results confirmed the accuracy of the transcriptome data. This study provides a quick strategy for exploring candidate genes and provides new insights for improving fiber strength in cotton.

Genome-Wide Dissection of the Genetic Basis for Drought Tolerance in Gossypium hirsutum L. Races

Drought seriously threats the growth and development of Gossypium hirsutum L. To dissect the genetic basis for drought tolerance in the G. hirsutum L. germplasm, a population, consisting of 188 accessions of G. hirsutum races and a cultivar (TM-1), was genotyped using the Cotton80KSNP biochip, and 51,268 high-quality single-nucleotide polymorphisms (SNPs) were obtained. Based on the phenotypic data of eight drought relative traits from four environments, we carried out association mapping with five models using GAPIT software. In total, thirty-six SNPs were detected significantly associated at least in two environments or two models. Among these SNPs, 8 and 28 (including 24 SNPs in 5 peak regions) were distributed in the A and D subgenome, respectively; eight SNPs were found to be distributed within separate genes. An SNP, TM73079, located on chromosome D10, was simultaneously associated with leaf fresh weight, leaf wilted weight, and leaf dry weight. Another nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported quantitative trait loci for drought tolerance. Moreover, 520 genes located 200 kb up- and down-stream of 36 SNPs were obtained and analyzed based on gene annotation and transcriptome sequencing. The results showed that three candidate genes, Gh_D08G2462, Gh_A03G0043, and Gh_A12G0369, may play important roles in drought tolerance. The current GWAS represents the first investigation into mapping QTL for drought tolerance in G. hirsutum races and provides important information for improving cotton cultivars.