High-density genetic variation maps reveal the correlation between asymmetric interspecific introgressions and improvement of agronomic traits in Upland and Pima cotton varieties developed in Xinjiang, China

Suppressing a β-1,3-glucanase gene expression increases the seed and fibre yield in cotton

Seeds are initiated from the carpel margin meristem (CMM) and high seed yield is top one of breeding objectives for many crops. β-1,3-glucanases play various roles in plant growth and developmental processes; however, whether it participates in CMM development and seed formation remains largely unknown. Here, we identified a β-1,3-glucanase gene (GLU19) as a determinant of CMM callose deposition and seed yield in cotton. GLU19 was differentially expressed in carpel tissues between Gossypium barbadense (Gb) and Gossypium hirsutum (Gh). Based on resequencing data, one interspecies-specific InDel in the promoter of GLU19 was further detected. The InDel was involved in the binding site of the CRABS CLAW (CRC) transcription factor, a regulator of carpel development. We found that the CRC binding affinity to the GLU19 promoter of G. barbadense was higher than that of G. hirsutum. Since G. barbadense yields fewer seeds than G. hirsutum, we speculated that stronger CRC binding to the GLU19 promoter activated higher expression of GLU19 which in turn suppressed seed production. Consistent with this hypothesis was that the overexpression of GhGLU19 caused reduced seed number, boll weight and less callose formation in CMM. Conversely, GhGLU19-knockdown (GhGLU19-KD) cotton led to the opposite phenotypes. By crossing GhGLU19-KD lines with several G. hirsutum and G. barbadense cotton accessions, all F1 and F2 plants carrying GhGLU19-KD transgenic loci exhibited higher seed yield than control plants without the locus. The increased seed effect was also found in the down-regulation of Arabidopsis orthologs lines, indicating that this engineering strategy may improve the seed yield in other crops.

Genetic Analysis of Cotton Fiber Traits in Gossypium Hybrid Lines

Cotton plays a crucial role in the progress of the textile industry and the betterment of human life by providing natural fibers. In our study, we explored the genetic determinants of cotton architecture and fiber yield and quality by crossbreeding Gossypium hirsutum and Gossypium barbadense, creating a recombinant inbred line (RIL) population. Utilizing SNP markers, we constructed an extensive genetic map encompassing 7,730 markers over 2,784.2 cM. We appraised two architectural and seven fiber traits within six environments, identifying 58 QTLs, of which 49 demonstrated stability across these environments. These encompassed QTLs for traits such as lint percentage (LP), boll weight (BW), fiber strength (STRENGTH), seed index (SI), and micronaire (MIC), primarily located on chromosomes chr-A07, chr-D06, and chr-D07. Notably, chr-D07 houses a QTL region affecting SI, corroborated by multiple studies. Within this region, the genes BZIP043 and SEP2 were identified as pivotal, with SEP2 particularly showing augmented expression in developing ovules. These discoveries contribute significantly to marker-assisted selection, potentially elevating both the yield and quality of cotton fiber production. These findings provide valuable insights into marker-assisted breeding strategies, offering crucial information to enhance fiber yield and quality in cotton production.

Combined genome and transcriptome analysis of elite fiber quality in Gossypium barbadense

Gossypium barbadense, which is one of several species of cotton, is well known for its superior fiber quality. However, the genetic basis of its high-quality fiber remains largely unexplored. Here, we resequenced 269 G. barbadense accessions. Phylogenetic structure analysis showed that the set of accessions was clustered into 3 groups: G1 and G2 mainly included modern cultivars from Xinjiang, China, and G3 was related to widely introduced accessions in different regions worldwide. A genome-wide association study of 5 fiber quality traits across multiple field environments identified a total of 512 qtls (main-effect QTLs) and 94 qtlEs (QTL-by-environment interactions) related to fiber quality, of which 292 qtls and 57 qtlEs colocated with previous studies. We extracted the genes located in these loci and performed expression comparison, local association analysis, and introgression segment identification. The results showed that high expression of hormone-related genes during fiber development, introgressions from Gossypium hirsutum, and the recombination of domesticated elite allelic variation were 3 major contributors to improve the fiber quality of G. barbadense. In total, 839 candidate genes with encoding region variations associated with elite fiber quality were mined. We confirmed that haplotype GB_D03G0092H traced to G. hirsutum introgression, with a 1-bp deletion leading to a frameshift mutation compared with GB_D03G0092B, significantly improved fiber quality. GB_D03G0092H is localized in the plasma membrane, while GB_D03G0092B is in both the nucleus and plasma membrane. Overexpression of GB_D03G0092H in Arabidopsis (Arabidopsis thaliana) significantly improved the elongation of longitudinal cells. Our study systematically reveals the genetic basis of the superior fiber quality of G. barbadense and provides elite segments and gene resources for breeding high-quality cotton cultivars.

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 and haplotype provide insights into the population differentiation and breeding improvement of Gossypium barbadense

INTRODUCTION

Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare.

OBJECTIVES

We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton.

METHODS

Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3-79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton.

RESULTS

One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian.

CONCLUSION

Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1

BACKGROUND

Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear.

RESULTS

In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR).

CONCLUSIONS

Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Genome-wide artificial introgressions of Gossypium barbadense into G. hirsutum reveal superior loci for simultaneous improvement of cotton fiber quality and yield traits

INTRODUCTION

The simultaneous improvement of fiber quality and yield for cotton is strongly limited by the narrow genetic backgrounds of Gossypium hirsutum (Gh) and the negative genetic correlations among traits. An effective way to overcome the bottlenecks is to introgress the favorable alleles of Gossypium barbadense (Gb) for fiber quality into Gh with high yield.

OBJECTIVES

This study was to identify superior loci for the improvement of fiber quality and yield.

METHODS

Two sets of chromosome segment substitution lines (CSSLs) were generated by crossing Hai1 (Gb, donor-parent) with cultivar CCRI36 (Gh) and CCRI45 (Gh) as genetic backgrounds, and cultivated in 6 and 8 environments, respectively. The kmer genotyping strategy was improved and applied to the population genetic analysis of 743 genomic sequencing data. A progeny segregating population was constructed to validate genetic effects of the candidate loci.

RESULTS

A total of 68,912 and 83,352 genome-wide introgressed kmers were identified in the CCRI36 and CCRI45 populations, respectively. Over 90 % introgressions were homologous exchanges and about 21 % were reverse insertions. In total, 291 major introgressed segments were identified with stable genetic effects, of which 66(22.98 %), 64(21.99 %), 35(12.03 %), 31(10.65 %) and 18(6.19 %) were beneficial for the improvement of fiber length (FL), strength (FS), micronaire, lint-percentage (LP) and boll-weight, respectively. Thirty-nine introgression segments were detected with stable favorable additive effects for simultaneous improvement of 2 or more traits in Gh genetic background, including 6 could increase FL/FS and LP. The pyramiding effects of 3 pleiotropic segments (A07:C45Clu-081, D06:C45Clu-218, D02:C45Clu-193) were further validated in the segregating population.

CONCLUSION

The combining of genome-wide introgressions and kmer genotyping strategy showed significant advantages in exploring genetic resources. Through the genome-wide comprehensive mining, a total of 11 clusters (segments) were discovered for the stable simultaneous improvement of FL/FS and LP, which should be paid more attention in the future.

Genetic diversity and population structure analyses and genome-wide association studies of photoperiod sensitivity in cotton (Gossypium hirsutum L.)

KEY MESSAGE

Genetic diversity and population structure analyses showed progressively narrowed diversity in US Upland cotton compared to land races. GWAS identified genomic regions and candidate genes for photoperiod sensitivity in cotton. Six hundred fifty-seven accessions that included elite cotton germplasm (DIV panel), lines of a public cotton breeding program (FB panel), and tropical landrace accessions (TLA panel) of Gossypium hirsutum L. were genotyped with cottonSNP63K array and phenotyped for photoperiod sensitivity under long day-length conditions. The genetic diversity analysis using 26,952 polymorphic SNPs indicated a progressively narrowed diversity from the landraces (0.230) to the DIV panel accessions (0.195) and FB panel (0.116). Structure analysis in the US germplasm identified seven subpopulations representing all four major regions of the US cotton belt. Three subpopulations were identified within the landrace accessions. The highest fixation index (FST) of 0.65 was found between landrace accessions of Guatemala and the Plains-type cultivars from Southwest cotton region while the lowest FST values were between the germplasms of Mid-South and Southeastern regions. Genome wide association studies (GWAS) of photoperiod response using 600 phenotyped accessions identified 14 marker trait associations spread across eight Upland cotton chromosomes. Six of these marker trait associations, on four chromosomes (A10, D04, D05, and D06), showed significant epistatic interactions. Targeted genomic analysis identified regions with 19 candidate genes including Transcription factor Vascular Plant One-Zinc Finger 1 (VOZ1) and Protein Photoperiod-Independent Early Flowering 1 (PIE1) genes. Genetic diversity and genome wide analyses of photoperiod sensitivity in diverse cotton germplasms will enable the use of genomic tools to systematically utilize the tropical germplasm and its beneficial alleles for broadening the genetic base in Upland cotton.

Genome-wide characterization, evolutionary analysis, and expression pattern analysis of the trihelix transcription factor family and gene expression analysis under MeJA treatment in Panax ginseng

Panax ginseng is a well-known medicinal plant with several pharmacological uses in China. The trihelix family transcription factors, also known as GT factors, can be involved in the regulation of growth and developmental processes in plants. There have been no in-depth reports or systematic studies about the trihelix transcription factor in ginseng. In this study, the structure, chromosomal localization, gene duplication, phylogeny, functional differentiation, expression patterns and coexpression interactions of trihelix transcripts were analysed using bioinformatics methods based on the ginseng transcriptome database. Thirty-two trihelix transcription factor genes were identified in ginseng, and these genes were alternatively spliced to obtain 218 transcripts. These transcripts were unevenly distributed on different chromosomes of ginseng, and phylogenetic analysis classified the PgGT transcripts into five subgroups. Gene Ontology (GO) analysis classified PgGT transcripts into eight functional subclasses, indicating that they are functionally diverse. The expression pattern analysis of 218 PgGT transcripts revealed that their expression was tissue-specific and spatiotemporally-specific in 14 different tissues of 4-year-old ginseng, 4 different ages of ginseng roots, and 42 farmers' cultivars of 4-year-old ginseng roots. Despite the differences in the expression patterns of these transcripts, coexpression network analysis revealed that these transcripts could be expressed synergistically in ginseng. In addition, two randomly selected PgGT transcripts in each of the five different subfamilies were subjected to methyl jasmonate treatment at different times, and PgGT was able to respond to the regulation of methy1 jasmonate. These results provide a theoretical basis and gene resources for an in-depth study of the function of trihelix genes in other plants.

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.

CottonMD: a multi-omics database for cotton biological study

Cotton is an important economic crop, and many loci for important traits have been identified, but it remains challenging and time-consuming to identify candidate or causal genes/variants and clarify their roles in phenotype formation and regulation. Here, we first collected and integrated the multi-omics datasets including 25 genomes, transcriptomes in 76 tissue samples, epigenome data of five species and metabolome data of 768 metabolites from four tissues, and genetic variation, trait and transcriptome datasets from 4180 cotton accessions. Then, a cotton multi-omics database (CottonMD, http://yanglab.hzau.edu.cn/CottonMD/) was constructed. In CottonMD, multiple statistical methods were applied to identify the associations between variations and phenotypes, and many easy-to-use analysis tools were provided to help researchers quickly acquire the related omics information and perform multi-omics data analysis. Two case studies demonstrated the power of CottonMD for identifying and analyzing the candidate genes, as well as the great potential of integrating multi-omics data for cotton genetic breeding and functional genomics research.

Genome-wide association study and transcriptome analysis reveal key genes controlling fruit branch angle in cotton

Fruit branch angle (FBA), a pivotal component of cotton plant architecture, is vital for field and mechanical harvesting. However, the molecular mechanism of FBA formation is poorly understood in cotton. To uncover the genetic basis for FBA formation in cotton, we performed a genome-wide association study (GWAS) of 163 cotton accessions with re-sequencing data. A total of 55 SNPs and 18 candidate genes were significantly associated with FBA trait. By combining GWAS and transcriptome analysis, four genes underlying FBA were identified. An FBA-associated candidate gene Ghi_A09G08736, which is homologous to SAUR46 in Arabidopsis thaliana, was detected in our study. In addition, transcriptomic evidence was provided to show that gravity and light were implicated in the FBA formation. This study provides new insights into the genetic architecture of FBA that informs architecture breeding in cotton.

Genetic Mapping and Analysis of a Compact Plant Architecture and Precocious Mutant in Upland Cotton

With the promotion and popularization of machine cotton-picking, more and more attention has been paid to the selection of early-maturity varieties with compact plant architecture. The type of fruit branch is one of the most important factors affecting plant architecture and early maturity of cotton. Heredity analysis of the cotton fruit branch is beneficial to the breeding of machine-picked cotton. Phenotype analysis showed that the types of fruit branches in cotton are controlled by a single recessive gene. Using an F2 population crossed with Huaxin102 (normal branch) and 04N-11 (nulliplex branch), BSA (Bulked Segregant Analysis) resequencing analysis and GhNB gene cloning in 04N-11, and allelic testing, showed that fruit branch type was controlled by the GhNB gene, located on chromosome D07. Ghnb5, a new recessive genotype of GhNB, was found in 04N-11. Through candidate gene association analysis, SNP 20_15811516_SNV was found to be associated with plant architecture and early maturity in the Xinjiang natural population. The GhNB gene, which is related to early maturity and the plant architecture of cotton, is a branch-type gene of cotton. The 20_15811516_SNV marker, obtained from the Xinjiang natural population, was used for the assisted breeding of machine-picked cotton varieties.

Construction of a high-density genetic map and identification of QTLs related to agronomic and physiological traits in an interspecific (Gossypium hirsutum × Gossypium barbadense) F2 population

Background

Advances in genome sequencing technology, particularly restriction-site associated DNA sequence (RAD-seq) and whole-genome resequencing, have greatly aided the construction of cotton interspecific genetic maps based on single nucleotide polymorphism (SNPs), Indels, and other types of markers. High-density genetic maps can improve accuracy of quantitative trait locus (QTL) mapping, narrow down location intervals, and facilitate identification of the candidate genes.

Result

In this study, 249 individuals from an interspecific F2 population (TM-1 and Hai7124) were re-sequenced, yielding 6303 high-confidence bin markers spanning 5057.13 cM across 26 cotton chromosomes. A total of 3380 recombination hot regions RHRs were identified which unevenly distributed on the 26 chromosomes. Based on this map, 112 QTLs relating to agronomic and physiological traits from seedling to boll opening stage were identified, including 15 loci associated with 14 traits that contained genes harboring nonsynonymous SNPs. We analyzed the sequence and expression of these ten candidate genes and discovered that GhRHD3 (GH_D10G0500) may affect fiber yield while GhGPAT6 (GH_D04G1426) may affect photosynthesis efficiency.

Conclusion

Our research illustrates the efficiency of constructing a genetic map using binmap and QTL mapping on the basis of a certain size of the early-generation population. High-density genetic map features high recombination exchanges in number and distribution. The QTLs and the candidate genes identified based on this high-density genetic map may provide important gene resources for the genetic improvement of cotton.

A Calmodulin-Like Gene (GbCML7) for Fiber Strength and Yield Improvement Identified by Resequencing Core Accessions of a Pedigree in Gossypium barbadense

Sea Island cotton (Gossypium barbadense) is world-renowned for its superior natural fiber. Although fiber strength is one of the most important fiber quality traits, genes contributing to fiber strength are poorly understood. Production of sea island cotton also is inextricably linked to improving its relatively low yield, thus enhancing the importance of joint improvement of both fiber quality and yield. We used genomic variation to uncover the genetic evidence of trait improvement resulting from pedigree breeding of Sea Island cotton. This pedigree was aimed at improving fiber strength and yielded an elite cultivar, XH35. Using a combination of genome-wide association study (GWAS) and selection screens, we detected 82 putative fiber-strength-related genes. Expression analysis confirmed a calmodulin-like gene, GbCML7, which enhanced fiber strength in a specific haplotype. This gene is a major-effect gene, which interacts with a minor-effect gene, GbTUA3, facilitating the enhancement of fiber strength in a synergistic fashion. Moreover, GbCML7 participates in the cooperative improvement of fiber strength, fiber length, and fiber uniformity, though a slight compromise exists between the first two of these traits and the latter. Importantly, GbCML7 is shown to boost yield in some backgrounds by increasing multiple yield components to varying degrees, especially boll number. Our work provides valuable genomic evidence and a key genetic factor for the joint improvement of fiber quality and yield in Sea Island cotton.

CottonGVD: A Comprehensive Genomic Variation Database for Cultivated Cottons

Cultivated cottons are the most important economic crop, which produce natural fiber for the textile industry. In recent years, the genetic basis of several essential traits for cultivated cottons has been gradually elucidated by decoding their genomic variations. Although an abundance of resequencing data is available in public, there is still a lack of a comprehensive tool to exhibit the results of genomic variations and genome-wide association study (GWAS). To assist cotton researchers in utilizing these data efficiently and conveniently, we constructed the cotton genomic variation database (CottonGVD; http://120.78.174.209/ or http://db.cngb.org/cottonGVD). This database contains the published genomic information of three cultivated cotton species, the corresponding population variations (SNP and InDel markers), and the visualized results of GWAS for major traits. Various built-in genomic tools help users retrieve, browse, and query the variations conveniently. The database also provides interactive maps (e.g., Manhattan map, scatter plot, heatmap, and linkage disequilibrium block) to exhibit GWAS and expression GWAS results. Cotton researchers could easily focus on phenotype-associated loci visualization, and they are interested in and screen for candidate genes. Moreover, CottonGVD will continue to update by adding more data and functions.

Whole-genome resequencing of 240 Gossypium barbadense accessions reveals genetic variation and genes associated with fiber strength and lint percentage

KEY MESSAGE

Genetic variation in a G. barbadense population was revealed using resquencing. GWAS on G.barbadense population identified several candidate genes associated with fiber strength and lint percentage. Gossypium barbadense is the second-largest cultivated cotton species planted in the world, which is characterized by high fiber quality. Here, we described the global pattern of genetic polymorphisms for 240 G. barbadense accessions based on the whole-genome resequencing. A total of 3,632,231 qualified single-nucleotide polymorphisms (SNPs) and 221,354 insertion-deletions (indels) were obtained. We conducted a genome-wide association study (GWAS) on 12 traits under four environments. Two traits with more stable associated variants, fiber strength and lint percentage, were chosen for further analysis. Three putative candidate genes, HD16 orthology (GB_D11G3437), WDL2 orthology (GB_D11G3460) and TUBA1 orthology (GB_D11G3471), on chromosome D11 were found to be associated with fiber strength, and one gene orthologous to Arabidopsis Receptor-like protein kinase HERK 1 (GB_A07G1034) was predicated to be the candidate gene for the lint percentage improvement. The identified genes may serve as promising targets for genetic engineering to accelerate the breeding process for G. barbadense and the high-density genome variation map constructed in this work may facilitate our understanding of the genetic architecture of cotton traits.

Divergent improvement of two cultivated allotetraploid cotton species

Interspecific genomic variation can provide a genetic basis for local adaptation and domestication. A series of studies have presented its role of interspecific haplotypes and introgressions in adaptive traits, but few studies have addressed their role in improving agronomic character. Two allotetraploid Gossypium species, Gossypium barbadense (Gb) and G. hirsutum (Gh) originating from the Americas, are cultivated independently. Here, through sequencing and the comparison of one GWAS panel in 229 Gb accessions and two GWAS panels in 491 Gh accessions, we found that most associated loci or functional haplotypes for agronomic traits were highly divergent, representing the strong divergent improvement between Gb and Gh. Using a comprehensive interspecific haplotype map, we revealed that six interspecific introgressions from Gh to Gb were significantly associated with the phenotypic performance of Gb, which could explain 5%-40% of phenotypic variation in yield and fibre qualities. In addition, three introgressions overlapped with six associated loci in Gb, indicating that these introgression regions were under further selection and stabilized during improvement. A single interspecific introgression often possessed yield-increasing potential but decreased fibre qualities, or the opposite, making it difficult to simultaneously improve yield and fibre qualities. Our study not only has proved the importance of interspecific functional haplotypes or introgressions in the divergent improvement of Gb and Gh, but also supports their potential value in further human-mediated hybridization or precision breeding.

Genomic mapping and identification of candidate genes encoding nulliplex-branch trait in sea-island cotton (Gossypium barbadense L.) by multi-omics analysis

Nulliplex branch is a key architectural trait in sea-island cotton (Gossypium barbadense L.), but its genetic basis is not well understood. Here we investigated the genetic basis of the nulliplex-branch trait in cotton by combining newly created bulked segregant analysis (BSA)-seq data, published RNA-seq data, and published whole-genome resequencing (WGR) data. We delimited the nulliplex-branch locus (qD07-NB) to D07, region 14.8-17.1 Mb, using various BSA methods and markers. We integrated our BSA data with WGR data of sea-island cotton varieties and detected a missense single-nucleotide polymorphism in the candidate gene (Gbar_D07G011870) of qD07-NB. This gene was under positive selection during sea-island cotton breeding in the Xinjiang Uygur Autonomous Region, China. Notably, the nulliplex-branch varieties possessed a better fiber quality than the long-branch varieties, and a set of high-quality molecular markers was identified for molecular breeding of the nulliplex-branch trait in cotton. We combined BSA-seq and RNA-seq data to compare gene expression profiles between two elite sea-island cotton varieties during three developmental stages. We identified eleven relevant candidate genes, five downregulated and six upregulated, in the qD07-NB locus. This research will expand our understanding of the genetic basis of the nulliplex-branch trait and provide guidance for architecture-focused breeding in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01229-w.

Evaluation and genome-wide association study of resistance to bacterial blight race 18 in U.S. 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.

Genome-Wide Association Analysis Reveals Loci and Candidate Genes Involved in Fiber Quality Traits Under Multiple Field Environments in Cotton (Gossypium hirsutum)

Fiber length, fiber strength, and fiber micronaire are the main fiber quality parameters in cotton. Thus, mining the elite and stable loci/alleles related to fiber quality traits and elucidating the relationship between the two may accelerate genetic improvement of fiber quality in cotton. Here, genome-wide association analysis (GWAS) was performed for fiber quality parameters based on phenotypic data, and 56,010 high-quality single nucleotide polymorphisms (SNPs) using 242 upland cotton accessions under 12 field environments were obtained. Phenotypic analysis exhibited that fiber length (FL) had a positive correlation with fiber strength (FS) and had a negative correlation with fiber micronaire (Mic). Genetic analysis also indicated that FL, FS, and Mic had high heritability of more than 80%. A total of 67 stable quantitative trait loci (QTLs) were identified through GWAS analysis, including 31 for FL, 21 for FS, and 22 for Mic. Of them, three pairs homologous QTLs were detected between A and D subgenomes, and seven co-located QTLs with two fiber quality parameters were found. Compared with the reported QTLs, 34 co-located with previous studies, and 33 were newly revealed. Integrated with transcriptome analysis, we selected 256, 244, and 149 candidate genes for FL, FS, and Mic, respectively. Gene Ontology (GO) analysis showed that most of the genes located in QTLs interval of the three fiber quality traits were involved in sugar biosynthesis, sugar metabolism, microtubule, and cytoskeleton organization, which played crucial roles in fiber development. Through correlation analysis between haplotypes and phenotypes, three genes (GH_A05G1494, GH_D11G3097, and GH_A05G1082) predominately expressed in fiber development stages were indicated to be potentially responsible for FL, FS, and Mic, respectively. The GH_A05G1494 encoded a protein containing SGS-domain, which is related to tubulin-binding and ubiquitin-protein ligase binding. The GH_D11G3097 encoded 20S proteasome beta subunit G1, and was involved in the ubiquitin-dependent protein catabolic process. The GH_A05G1082 encoded RAN binding protein 1 with a molecular function of GTPase activator activity. These results provide new insights and candidate loci/genes for the improvement of fiber quality in cotton.

Genome-wide association analysis reveals loci and candidate genes involved in fiber quality traits in sea island cotton (Gossypium barbadense)

BACKGROUND

Sea island cotton (Gossypium barbadense) has markedly superior high quality fibers, which plays an important role in the textile industry and acts as a donor for upland cotton (G. hirsutum) fiber quality improvement. The genetic characteristics analysis and the identification of key genes will be helpful to understand the mechanism of fiber development and breeding utilization in sea island cotton.

RESULTS

In this study, 279 sea island cotton accessions were collected from different origins for genotyping and phenotyping fiber quality traits. A set of 6303 high quality single nucleotide polymorphisms (SNPs) were obtained by high-density CottonSNP80K array. The population characteristics showed that the sea island cotton accessions had wide genetic diversity and were clustered into three groups, with Group1 closely related to Menoufi, an original sea island cotton landrace, and Group2 and Group3 related to widely introduced accessions from Egypt, USA and Former Soviet Union. Further, we used 249 accessions and evaluated five fiber quality traits under normal and salt environments over 2 years. Except for fiber uniformity (FU), fiber length (FL) and fiber elongation (FE) were significantly decreased in salt conditions, while fiber strength (FS) and fiber micronaire (MIC) were increased. Based on 6303 SNPs and genome-wide association study (GWAS) analysis, a total of 34 stable quantitative trait loci (QTLs) were identified for the five fiber quality traits with 25 detected simultaneously under normal and salt environments. Gene Ontology (GO) analysis indicated that candidate genes in the 25 overlapped QTLs were enriched mostly in "cellular and biological process". In addition, "xylem development" and "response to hormone" pathways were also found. Haplotype analyses found that GB_A03G0335 encoding an E3 ubiquitin-protein ligase in QTL TM6004 had SNP variation (A/C) in gene region, was significantly correlated with FL, FS, FU, and FE, implying a crucial role in fiber quality.

CONCLUSIONS

The present study provides a foundation for genetic diversity of sea island cotton accessions and will contribute to fiber quality improvement in breeding practice.