Unraveling genomic regions and candidate genes for multiple disease resistance in upland cotton using meta-QTL analysis.

Verticillium wilt (VW), Fusarium wilt (FW) and Root-knot nematode (RKN) are the main diseases affecting cotton production. However, many reported quantitative trait loci (QTLs) for cotton resistance have not been used for agricultural practices because of inconsistencies in the cotton genetic background. The integration of existing cotton genetic resources can facilitate the discovery of important genomic regions and candidate genes involved in disease resistance. Here, an improved and comprehensive meta-QTL analysis was conducted on 487 disease resistant QTLs from 31 studies in the last two decades. A consensus linkage map with genetic overall length of 3006.59 cM containing 8650 markers was constructed. A total of 28 Meta-QTLs (MQTLs) were discovered, among which nine MQTLs were identified as related to resistance to multiple diseases. Candidate genes were predicted based on public transcriptome data and enriched in pathways related to disease resistance. This study used a method based on the integration of Meta-QTL, known genes and transcriptomics to reveal major genomic regions and putative candidate genes for resistance to multiple diseases, providing a new basis for marker-assisted selection of high disease resistance in cotton breeding.

GhSINA1, a SEVEN in ABSENTIA ubiquitin ligase, negatively regulates fiber development in Upland cotton.

Protein ubiquitination is essential for plant growth and responses to the environment. The SEVEN IN ABSENTIA (SINA) ubiquitin ligases have been extensively studied in plants, but information on their roles in fiber development is limited. Here, we identified GhSINA1 in Upland cotton (Gossypium hirsutum), which has a conserved RING finger domain and SINA domain. Quantitative real-time PCR (qRT-PCR) analysis showed that GhSINA1 was preferentially expressed during fiber initiation and elongation, especially during initiation in the fuzzless-lintless cotton mutant. Subcellular localization experiments indicated that GhSINA1 localized to the nucleus. In vitro ubiquitination analysis revealed that GhSINA1 has E3 ubiquitin ligase activity. Ectopic overexpression of GhSINA1 in Arabidopsis thaliana reduced the number and length of root hairs and trichomes. Yeast two-hybrid (Y2H), firefly luciferase complementation imaging (LCI), and bimolecular fluorescence complementation (BiFC) assays demonstrated that the GhSINA1 proteins could interact with each other to form homodimers and heterodimers. Overall, these results suggest that GhSINA1 may act as a negative regulator in cotton fiber development through homodimerization and heterodimerization.

Genetic structure, gene flow pattern, and association analysis of superior germplasm resources in domesticated upland cotton (Gossypium hirsutum L.).

Gene flow patterns and the genetic structure of domesticated crops like cotton are not well understood. Furthermore, marker-assisted breeding of cotton has lagged far behind that of other major crops because the loci associated with cotton traits such as fiber yield and quality have scarcely been identified. In this study, we used 19 microsatellites to first determine the population genetic structure and patterns of gene flow of superior germplasm resources in upland cotton. We then used association analysis to identify which markers were associated with 15 agronomic traits (including ten yield and five fiber quality traits). The results showed that the upland cotton accessions have low levels of genetic diversity (polymorphism information content = 0.427), although extensive gene flow occurred among different ecological and geographic regions. Bayesian clustering analysis indicated that the cotton resources used in this study did not belong to obvious geographic populations, which may be the consequence of a single source of domestication followed by frequent genetic introgression mediated by human transference. A total of 82 maker-trait associations were examined in association analysis and the related ratios for phenotypic variations ranged from 3.04% to 47.14%. Interestingly, nine SSR markers were detected in more than one environmental condition. In addition, 14 SSR markers were co-associated with two or more different traits. It was noteworthy that NAU4860 and NAU5077 markers detected at least in two environments were simultaneously associated with three fiber quality traits (uniformity index, specific breaking strength and micronaire value). In conclusion, these findings provide new insights into the population structure and genetic exchange pattern of cultivated cotton accessions. The quantitative trait loci of domesticated cotton identified will also be very useful for improvement of yield and fiber quality of cotton in molecular breeding programs.

Comparative Chloroplast Genomics of Gossypium Species: Insights Into Repeat Sequence Variations and Phylogeny.

Cotton is one of the most economically important fiber crop plants worldwide. The genus Gossypium contains a single allotetraploid group (AD) and eight diploid genome groups (A-G and K). However, the evolution of repeat sequences in the chloroplast genomes and the phylogenetic relationships of Gossypium species are unclear. Thus, we determined the variations in the repeat sequences and the evolutionary relationships of 40 cotton chloroplast genomes, which represented the most diverse in the genus, including five newly sequenced diploid species, i.e., G. nandewarense (C1-n), G. armourianum (D2-1), G. lobatum (D7), G. trilobum (D8), and G. schwendimanii (D11), and an important semi-wild race of upland cotton, G. hirsutum race latifolium (AD1). The genome structure, gene order, and GC content of cotton species were similar to those of other higher plant plastid genomes. In total, 2860 long sequence repeats (>10 bp in length) were identified, where the F-genome species had the largest number of repeats (G. longicalyx F1: 108) and E-genome species had the lowest (G. stocksii E1: 53). Large-scale repeat sequences possibly enrich the genetic information and maintain genome stability in cotton species. We also identified 10 divergence hotspot regions, i.e., rpl33-rps18, psbZ-trnG (GCC), rps4-trnT (UGU), trnL (UAG)-rpl32, trnE (UUC)-trnT (GGU), atpE, ndhI, rps2, ycf1, and ndhF, which could be useful molecular genetic markers for future population genetics and phylogenetic studies. Site-specific selection analysis showed that some of the coding sites of 10 chloroplast genes (atpB, atpE, rps2, rps3, petB, petD, ccsA, cemA, ycf1, and rbcL) were under protein sequence evolution. Phylogenetic analysis based on the whole plastomes suggested that the Gossypium species grouped into six previously identified genetic clades. Interestingly, all 13 D-genome species clustered into a strong monophyletic clade. Unexpectedly, the cotton species with C, G, and K-genomes were admixed and nested in a large clade, which could have been due to their recent radiation, incomplete lineage sorting, and introgression hybridization among different cotton lineages. In conclusion, the results of this study provide new insights into the evolution of repeat sequences in chloroplast genomes and interspecific relationships in the genus Gossypium.