Genomic loci associated with leaf abscission contribute to machine picking and environmental adaptability in upland cotton (Gossypium hirsutum L.).

INTRODUCTION: Defoliation by applying defoliants before machine picking is an important agricultural practice that enhances harvesting efficiency and leads to increased raw cotton purity. However, the fundamental characteristics of leaf abscission and the underlying genetic basis in cotton are not clearly understood. OBJECTIVES: In this study, we aimed to (1) reveal the phenotypic variations in cotton leaf abscission, (2) discover the whole-genome differentiation sweeps and genetic loci related to defoliation, (3) identify and verify the functions of key candidate genes associated with defoliation, and (4) explore the relationship between haplotype frequency of loci and environmental adaptability. METHODS: Four defoliation-related traits of 383 re-sequenced Gossypium hirsutum accessions were investigated in four environments. The genome-wide association study (GWAS), linkage disequilibrium (LD) interval genotyping and functional identification were conducted. Finally, the haplotype variation related to environmental adaptability and defoliation traits was revealed. RESULTS: Our findings revealed the fundamental phenotypic variations of defoliation traits in cotton. We showed that defoliant significantly increased the defoliation rate without incurring yield and fiber quality penalties. The strong correlations between defoliation traits and growth period traits were observed. A genome-wide association study of defoliation traits identified 174 significant SNPs. Two loci (RDR7 on A02 and RDR13 on A13) that significantly associated with the relative defoliation rate were described, and key candidate genes GhLRR and GhCYCD3;1, encoding a leucine-rich repeat (LRR) family protein and D3-type cell cyclin 1 protein respectively, were functional verified by expression pattern analysis and gene silencing. We found that combining of two favorable haplotypes (HapRDR7 and HapRDR13) improved sensitivity to defoliant. The favorable haplotype frequency generally increased in high latitudes in China, enabling adaptation to the local environment. CONCLUSION: Our findings lay an important foundation for the potentially broad application of leveraging key genetic loci in breeding machine-pickable cotton.

Physiological and transcriptional analyses reveal formation of memory under recurring drought stresses in seedlings of cotton (Gossypium hirsutum).

Plants are frequently subjected to a range of environmental stresses, including drought, salinity, cold, pathogens, and herbivore attacks. To survive in such conditions, plants have evolved a novel adaptive mechanism known as 'stress memory'. The formation of stress memories necessitates coordinated responses at the cellular, genetic/genomic, and epigenetic levels, involving altered physiological responses, gene activation, hyper-induction and chromatin modification. Cotton (Gossypium spp.) is an important economic crop with numerous applications and high economic value. In this study, we establish G. hirsutum drought memory following cycles of mild drought and re-watering treatments and analyzed memory gene expression patterns. Our findings reveal the physiological, biochemical, and molecular mechanisms underlying drought stress memory formation in G. hirsutum. Specifically, H3K4me3, a histone modification, plays a crucial role in regulating [+ /+ ] transcriptional memory. Moreover, we investigated the intergenerational inheritance of drought stress memory in G. hirsutum. Collectively, our data provides theoretical guidance for cotton breeding.

Identification and Characterization of the Growth-Regulating Factors-Interacting Factors in Cotton.

Growth-regulating factors-interacting factors (GIFs) are a type of transcription co-activators in plants, playing crucial roles in plants' growth, development, and stress adaptation. Here, a total of 35 GIF genes were identified and clustered into two groups by phylogenetic analysis in four cotton genus. The gene structure and conserved domain analysis proved the conservative characteristics of GIF genes in cotton. The function of GIF genes was evaluated in two cotton accessions, Ji A-1-7 (33xi) and King, which have larger and smaller lateral root numbers, respectively. The results showed that the expression of GhGIF4 in Ji A-1-7 (33xi) was higher than that in King. The enzyme activity and microstructure assay showed a higher POD activity, lower MDA content, and more giant cells of the lateral root emergence part phenotype in Ji A-1-7 (33xi) than in King. A mild waterlogging assay showed the GIF genes were down-regulated in the waterlogged seedling. Further confirmation of the suppression of GhGIF4 in cotton plants further confirmed that GhGIF4 could reduce the lateral root numbers in cotton. This study could provide a basis for future studies of the role of GIF genes in upland cotton.

Integrating Genome-wide association and whole transcriptome analysis to reveal genetic control of leaf traits in Gossypium arboreum L.

Leaves are important organs for crop photosynthesis and transpiration, and their morphological characteristics can directly reflect the growth state of plants. Accurate measurement of leaf traits and mining molecular markers are of great significance to the study of cotton growth. Here, we performed a Genome-wide association study on 7 leaf traits in 213 Asian cotton accessions. 32 significant SNPs and 44 genes were identified. A field experiment showed significant difference in leaf hair and leaf area between DPL971 and its natural mutant DPL972. We also compared the leaf transcriptome difference between DPL971 and DPL972, and found a batch of differentially expressed genes and non-coding RNAs (including lncRNAs, microRNAs, and circRNAs). After integrating the GWAS and transcriptome results, we finally selected two coding genes (Ga03G2383 and Ga05G3412) and two microRNAs (hbr-miR156, unconservative_Chr03_contig343_2364) as the candidate for leaf traits. Those findings will provide important genomic resources for cotton leaf improvement breeding.

Genome-wide association study for seedling biomass-related traits in Gossypium arboreum L.

BACKGROUND: Seedling stage plant biomass is usually used as an auxiliary trait to study plant growth and development or stress adversities. However, few molecular markers and candidate genes of seedling biomass-related traits were found in cotton. RESULT: Here, we collected 215 Gossypium arboreum accessions, and investigated 11 seedling biomass-related traits including the fresh weight, dry weight, water content, and root shoot ratio. A genome-wide association study (GWAS) utilizing 142,5003 high-quality SNPs identified 83 significant associations and 69 putative candidate genes. Furthermore, the transcriptome profile of the candidate genes emphasized higher expression of Ga03G1298, Ga09G2054, Ga10G1342, Ga11G0096, and Ga11G2490 in four representative cotton accessions. The relative expression levels of those five genes were further verified by qRT-PCR. CONCLUSIONS: The significant SNPs, candidate genes identified in this study are expected to lay a foundation for studying the molecular mechanism for early biomass development and related traits in Asian cotton.

Genetic Factors Underlying Single Fiber Quality in A-Genome Donor Asian Cotton (Gossypium arboreum).

The study of A-genome Asian cotton as a potential fiber donor in Gossypium species may offer an enhanced understanding of complex genetics and novel players related to fiber quality traits. Assessment of individual fibers providing classified fiber quality information to the textile industry is Advanced Fiber Information System (AFIS) in the recent technological era. Keeping the scenario, a diverse collection of 215 Asiatic cotton accessions were evaluated across three agro-ecological zones of China. Genome-Wide Association Studies (GWAS) was performed to detect association signals related to 17 AFIS fiber quality traits grouped into four categories viz: NEPs, fiber length, maturity, and fineness. Significant correlations were found within as well as among different categories of various traits related to fiber quality. Fiber fineness has shown a strong correlation to all other categories, whereas these categories are shown interrelationships via fiber-fineness. A total of 7,429 SNPs were found in association with 17 investigated traits, of which 177 were selected as lead SNPs. In the vicinity of these lead SNPs, 56 differentially expressed genes in various tissues/development stages were identified as candidate genes. This compendium connecting trait-SNP-genes may allow further prioritization of genes in GWAS loci to enable mechanistic studies. These identified quantitative trait nucleotides (QTNs) may prove helpful in fiber quality improvement in Asian cotton through marker-assisted breeding as well as in reviving eroded genetic factors of G. hirsutum via introgression breeding.

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum).

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (- log10(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

GWAS Mediated Elucidation of Heterosis for Metric Traits in Cotton (Gossypium hirsutum L.) Across Multiple Environments.

For about a century, plant breeding has widely exploited the heterosis phenomenon-often considered as hybrid vigor-to increase agricultural productivity. The ensuing F1 hybrids can substantially outperform their progenitors due to heterozygous combinations that mitigate deleterious mutations occurring in each genome. However, only fragmented knowledge is available concerning the underlying genes and processes that foster heterosis. Although cotton is among the highly valued crops, its improvement programs that involve the exploitation of heterosis are still limited in terms of significant accomplishments to make it broadly applicable in different agro-ecological zones. Here, F1 hybrids were derived from mating a diverse Upland Cotton germplasm with commercially valuable cultivars in the Line × Tester fashion and evaluated across multiple environments for 10 measurable traits. These traits were dissected into five different heterosis types and specific combining ability (SCA). Subsequent genome-wide predictions along-with association analyses uncovered a set of 298 highly significant key single nucleotide polymorphisms (SNPs)/Quantitative Trait Nucleotides (QTNs) and 271 heterotic Quantitative Trait Nucleotides (hQTNs) related to agronomic and fiber quality traits. The integration of a genome wide association study with RNA-sequence analysis yielded 275 candidate genes in the vicinity of key SNPs/QTNs. Fiber micronaire (MIC) and lint percentage (LP) had the maximum number of associated genes, i.e., each with 45 related to QTNs/hQTNs. A total of 54 putative candidate genes were identified in association with HETEROSIS of quoted traits. The novel players in the heterosis mechanism highlighted in this study may prove to be scientifically and biologically important for cotton biologists, and for those breeders engaged in cotton fiber and yield improvement programs.

The genomic basis of geographic differentiation and fiber improvement in cultivated cotton.

Large-scale genomic surveys of crop germplasm are important for understanding the genetic architecture of favorable traits. The genomic basis of geographic differentiation and fiber improvement in cultivated cotton is poorly understood. Here, we analyzed 3,248 tetraploid cotton genomes and confirmed that the extensive chromosome inversions on chromosomes A06 and A08 underlies the geographic differentiation in cultivated Gossypium hirsutum. We further revealed that the haplotypic diversity originated from landraces, which might be essential for understanding adaptative evolution in cultivated cotton. Introgression and association analyses identified new fiber quality-related loci and demonstrated that the introgressed alleles from two diploid cottons had a large effect on fiber quality improvement. These loci provided the potential power to overcome the bottleneck in fiber quality improvement. Our study uncovered several critical genomic signatures generated by historical breeding effects in cotton and a wealth of data that enrich genomic resources for the research community.

Large-fragment insertion activates gene GaFZ (Ga08G0121) and is associated with the fuzz and trichome reduction in cotton (Gossypium arboreum).

Cotton seeds are typically covered by lint and fuzz fibres. Natural 'fuzzless' mutants are an ideal model system for identifying genes that regulate cell initiation and elongation. Here, using a genome-wide association study (GWAS), we identified a ~ 6.2 kb insertion, larINDELFZ , located at the end of chromosome 8, composed of a ~ 5.0 kb repetitive sequence and a ~ 1.2 kb fragment translocated from chromosome 12 in fuzzless Gossypium arboreum. The presence of larINDELFZ was associated with a fuzzless seed and reduced trichome phenotypes in G. arboreum. This distant insertion was predicted to be an enhancer, located ~ 18 kb upstream of the dominant-repressor GaFZ (Ga08G0121). Ectopic overexpression of GaFZ in Arabidopsis thaliana and G. hirsutum suggested that GaFZ negatively modulates fuzz and trichome development. Co-expression and interaction analyses demonstrated that GaFZ might impact fuzz fibre/trichome development by repressing the expression of genes in the very-long-chain fatty acid elongation pathway. Thus, we identified a novel regulator of fibre/trichome development while providing insights into the importance of noncoding sequences in cotton.

Genome-wide dissection of hybridization for fiber quality- and yield-related traits in upland cotton.

An evaluation of combining ability can facilitate the selection of suitable parents and superior F1 hybrids for hybrid cotton breeding, although the molecular genetic basis of combining ability has not been fully characterized. In the present study, 282 female parents were crossed with four male parents in accordance with the North Carolina II mating scheme to generate 1128 hybrids. The parental lines were genotyped based on restriction site-associated DNA sequencing and 306 814 filtered single nucleotide polymorphisms were used for genome-wide association analysis involving the phenotypes, general combining ability (GCA) values, and specific combining ability values of eight fiber quality- and yield-related traits. The main results were: (i) all parents could be clustered into five subgroups based on population structure analyses and the GCA performance of the female parents had significant differences between subgroups; (ii) 20 accessions with a top 5% GCA value for more than one trait were identified as elite parents for hybrid cotton breeding; (iii) 120 significant single nucleotide polymorphisms, clustered into 66 quantitative trait loci, such as the previously reported Gh_A07G1769 and GhHOX3 genes, were found to be significantly associated with GCA; and (iv) identified quantitative trait loci for GCA had a cumulative effect on GCA of the accessions. Overall, our results suggest that pyramiding the favorable loci for GCA may improve the efficiency of hybrid cotton breeding.

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.

Integration of conventional and advanced molecular tools to track footprints of heterosis in cotton.

BACKGROUND: Heterosis, a multigenic complex trait extrapolated as sum total of many phenotypic features, is widely utilized phenomenon in agricultural crops for about a century. It is mainly focused on establishing vigorous cultivars with the fact that its deployment in crops necessitates the perspective of genomic impressions on prior selection for metric traits. In spite of extensive investigations, the actual mysterious genetic basis of heterosis is yet to unravel. Contemporary crop breeding is aimed at enhanced crop production overcoming former achievements. Leading cotton improvement programs remained handicapped to attain significant accomplishments. RESULTS: In mentioned context, a comprehensive project was designed involving a large collection of cotton accessions including 284 lines, 5 testers along with their respective F1 hybrids derived from Line × Tester mating design were evaluated under 10 diverse environments. Heterosis, GCA and SCA were estimated from morphological and fiber quality traits by L × T analysis. For the exploration of elite marker alleles related to heterosis and to provide the material carrying such multiple alleles the mentioned three dependent variables along with trait phenotype values were executed for association study aided by microsatellites in mixed linear model based on population structure and linkage disequilibrium analysis. Highly significant 46 microsatellites were discovered in association with the fiber and yield related traits under study. It was observed that two-thirds of the highly significant associated microsatellites related to fiber quality were distributed on D sub-genome, including some with pleiotropic effect. Newly discovered 32 hQTLs related to fiber quality traits are one of prominent findings from current study. A set of 96 exclusively favorable alleles were discovered and C tester (A971Bt) posited a major contributor of these alleles primarily associated with fiber quality. CONCLUSIONS: Hence, to uncover hidden facts lying within heterosis phenomenon, discovery of additional hQTLs is required to improve fibre quality. To grab prominent improvement in influenced fiber quality and yield traits, we suggest the A971 Bt cotton cultivar as fundamental element in advance breeding programs as a parent of choice.

Dissection of complicate genetic architecture and breeding perspective of cottonseed traits by genome-wide association study.

BACKGROUND: Cottonseed is one of the most important raw materials for plant protein, oil and alternative biofuel for diesel engines. Understanding the complex genetic basis of cottonseed traits is requisite for achieving efficient genetic improvement of the traits. However, it is not yet clear about their genetic architecture in genomic level. GWAS has been an effective way to explore genetic basis of quantitative traits in human and many crops. This study aims to dissect genetic mechanism seven cottonseed traits by a GWAS for genetic improvement. RESULTS: A genome-wide association study (GWAS) based on a full gene model with gene effects as fixed and gene-environment interaction as random, was conducted for protein, oil and 5 fatty acids using 316 accessions and ~ 390 K SNPs. Totally, 124 significant quantitative trait SNPs (QTSs), consisting of 16, 21, 87 for protein, oil and fatty acids (palmitic, linoleic, oleic, myristic, stearic), respectively, were identified and the broad-sense heritability was estimated from 71.62 to 93.43%; no QTS-environment interaction was detected for the protein, the palmitic and the oleic contents; the protein content was predominantly controlled by epistatic effects accounting for 65.18% of the total variation, but the oil content and the fatty acids except the palmitic were mainly determined by gene main effects and no epistasis was detected for the myristic and the stearic. Prediction of superior pure line and hybrid revealed the potential of the QTSs in the improvement of cottonseed traits, and the hybrid could achieve higher or lower genetic values compared with pure lines. CONCLUSIONS: This study revealed complex genetic architecture of seven cottonseed traits at whole genome-wide by mixed linear model approach; the identified genetic variants and estimated genetic component effects of gene, gene-gene and gene-environment interaction provide cotton geneticist or breeders new knowledge on the genetic mechanism of the traits and the potential molecular breeding design strategy.

Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield.

Upland cotton is the most important natural-fiber crop. The genomic variation of diverse germplasms and alleles underpinning fiber quality and yield should be extensively explored. Here, we resequenced a core collection comprising 419 accessions with 6.55-fold coverage depth and identified approximately 3.66 million SNPs for evaluating the genomic variation. We performed phenotyping across 12 environments and conducted genome-wide association study of 13 fiber-related traits. 7,383 unique SNPs were significantly associated with these traits and were located within or near 4,820 genes; more associated loci were detected for fiber quality than fiber yield, and more fiber genes were detected in the D than the A subgenome. Several previously undescribed causal genes for days to flowering, fiber length, and fiber strength were identified. Phenotypic selection for these traits increased the frequency of elite alleles during domestication and breeding. These results provide targets for molecular selection and genetic manipulation in cotton improvement.

Association mapping analysis of fiber yield and quality traits in Upland cotton (Gossypium hirsutum L.).

Fiber yield and quality are the most important traits for Upland cotton (Gossypium hirsutum L.). Identifying high yield and good fiber quality genes are the prime concern of researchers in cotton breeding. Association mapping offers an alternative and powerful method for detecting those complex agronomic traits. In this study, 198 simple sequence repeats (SSRs) were used to screen markers associated with fiber yield and quality traits with 302 elite Upland cotton accessions that were evaluated in 12 locations representing the Yellow River and Yangtze River cotton growing regions of China. Three subpopulations were found after the estimation of population structure. The pair-wise kinship values varied from 0 to 0.867. Only 1.59% of the total marker locus pairs showed significant linkage disequilibrium (LD, p < 0.001). The genome-wide LD decayed within the genetic distance of ~30 to 32 cM at r 2 = 0.1, and decreased to ~1 to 2 cM at r 2 = 0.2, indicating the potential for association mapping. Analysis based on a mixed linear model detected 57 significant (p < 0.01) marker-trait associations, including seven associations for fiber length, ten for fiber micronaire, nine for fiber strength, eight for fiber elongation, five for fiber uniformity index, five for fiber uniformity ratio, six for boll weight and seven for lint percent, for a total of 35 SSR markers, of which 11 markers were associated with more than one trait. Among marker-trait associations, 24 associations coincided with the previously reported quantitative trait loci (QTLs), the remainder were newly identified QTLs/genes. The QTLs identified in this study will potentially facilitate improvement of fiber yield and quality in the future cotton molecular breeding programs.

Association mapping for epistasis and environmental interaction of yield traits in 323 cotton cultivars under 9 different environments.

Improving yield is a major objective for cotton breeding schemes, and lint yield and its three component traits (boll number, boll weight and lint percentage) are complex traits controlled by multiple genes and various environments. Association mapping was performed to detect markers associated with these four traits using 651 simple sequence repeats (SSRs). A mixed linear model including epistasis and environmental interaction was used to screen the loci associated with these four yield traits by 323 accessions of Gossypium hirsutum L. evaluated in nine different environments. 251 significant loci were detected to be associated with lint yield and its three components, including 69 loci with individual effects and all involved in epistasis interactions. These significant loci explain ∼ 62.05% of the phenotypic variance (ranging from 49.06% ∼ 72.29% for these four traits). It was indicated by high contribution of environmental interaction to the phenotypic variance for lint yield and boll numbers, that genetic effects of SSR loci were susceptible to environment factors. Shared loci were also observed among these four traits, which may be used for simultaneous improvement in cotton breeding for yield traits. Furthermore, consistent and elite loci were screened with -Log10 (P-value) >8.0 based on predicted effects of loci detected in different environments. There was one locus and 6 pairs of epistasis for lint yield, 4 loci and 10 epistasis for boll number, 15 loci and 2 epistasis for boll weight, and 2 loci and 5 epistasis for lint percentage, respectively. These results provided insights into the genetic basis of lint yield and its components and may be useful for marker-assisted breeding to improve cotton production.

[Genetic analysis of fuzzless in cotton germplasm].

The present study was conducted to evaluate genetic analysis of fuzzless seed trait in cotton. One hundred and two upland cotton (G. hirsutum) and eighty-five island cotton (G. barbadense) were used to cross with the same lines, TM-1 (G. hirsutum) and Xinhai 13 (G. barbadense), respectively. Two different F1 populations obtained were assessed to specify the dominant and recessive inheritance of fiber fuzziness in these lines. Three F1 populations (Kuguangzi × TM-1, Luwuxu × TM-1, and SA65 × TM-1) displaying recessive fiber fuzziness inheritance were selected to construct the F2 population for a further genetic study of fuzzless seed trait. The results of this study indicated that (1) the same materials showed different quantities of fuzzy fiber in different environments. Less fuzzy fiber was found in Xinjian and Hainan compared to Anyang. Thus, the quantity of fuzzy cotton seed depends on ecological environment. (2) In upland cotton, the inheritance of fiber fuzziness was dominant for 26 accessions (25.49%), incompletely dominant for 8 accessions (7.84%), and recessive for 22 accessions (21.57%). The inheritance of fiber fuzziness in island cotton was dominant for 5 accessions (5.88%), incompletely dominant for 16 accessions (18.82%), and recessive for 9 accessions (10.59%). Analysis of F2 population indicated that the fiber fuzziness of Kuguangzi was controlled by two recessive complementary effect alleles. The fiber fuzziness of Luwuxu was controlled by two recessive additive effect alleles, and a single recessive gene controlled the same trait for SA65. Fiber fuzziness evaluation in cotton germplasm provides the genetic and basic information for cotton fiber development study and breeding.