Genome-wide association study identified genetic variations and candidate genes for plant architecture component traits in Chinese upland cotton

Key message

Thirty significant associations between 22 SNPs and five plant architecture component traits in Chinese upland cotton were identified via GWAS. Four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits. A candidate gene, Gh_D03G0922, might be responsible for plant height in upland cotton.

Abstract

A compact plant architecture is increasingly required for mechanized harvesting processes in China. Therefore, cotton plant architecture is an important trait, and its components, such as plant height, fruit branch length and fruit branch angle, affect the suitability of a cultivar for mechanized harvesting. To determine the genetic basis of cotton plant architecture, a genome-wide association study (GWAS) was performed using a panel composed of 355 accessions and 93,250 single nucleotide polymorphisms (SNPs) identified using the specific-locus amplified fragment sequencing method. Thirty significant associations between 22 SNPs and five plant architecture component traits were identified via GWAS. Most importantly, four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits, and these SNPs were harbored in one linkage disequilibrium block. Furthermore, 21 candidate genes for plant architecture were predicted in a 0.95-Mb region including the four peak SNPs. One of these genes (Gh_D03G0922) was near the significant SNP D03_31584163 (8.40 kb), and its Arabidopsis homologs contain MADS-box domains that might be involved in plant growth and development. qRT-PCR showed that the expression of Gh_D03G0922 was upregulated in the apical buds and young leaves of the short and compact cotton varieties, and virus-induced gene silencing (VIGS) proved that the silenced plants exhibited increased PH. These results indicate that Gh_D03G0922 is likely the candidate gene for PH in cotton. The genetic variations and candidate genes identified in this study lay a foundation for cultivating moderately short and compact varieties in future Chinese cotton-breeding programs.

     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

A genome‐wide association study of early‐maturation traits in upland cotton based on the CottonSNP80K array

Genome‐wide association studies (GWASs) efficiently identify genetic loci controlling traits at a relatively high resolution. In this study, variations in major early‐maturation traits, including seedling period (SP), bud period (BP), flower and boll period (FBP), and growth period (GP), of 169 upland cotton accessions were investigated, and a GWAS of early maturation was performed based on a CottonSNP80K array. A total of 49,650 high‐quality single‐nucleotide polymorphisms (SNPs) were screened, and 29 significant SNPs located on chromosomes A6, A7, A8, D1, D2, and D9, were repeatedly identified as associated with early‐maturation traits, in at least two environments or two algorithms. Of these 29 significant SNPs, 1, 12, 11, and 5 were related to SP, BP, FBP, and GP, respectively. Six peak SNPs, TM47967, TM13732, TM20937, TM28428, TM50283, and TM72552, exhibited phenotypic contributions of approximately 10%, which could allow them to be used for marker‐assisted selection. One of these, TM72552, as well as four other SNPs, TM72554, TM72555, TM72558, and TM72559, corresponded to the quantitative trait loci previously reported. In total, 274 candidate genes were identified from the genome sequences of upland cotton and were categorized based on their functional annotations. Finally, our studies identified Gh_D01G0340 and Gh_D01G0341 as potential candidate genes for improving cotton early maturity.     

A genome-wide association study of early-maturation traits in upland cotton based on the CottonSNP8oK array

Genome-wide association studies(GWASs)efficiently identify genetic loci controlling traits at a relatively high resolution.In this study,variations in major early-maturation traits,including seedling period(SP),bud period(BP),flower and boll period(FBP),and growth period(GP),of 169 upland cotton accessions were investigated,and a GWAS of early maturation was performed based on a CottonSNP80K array.A total of49,650 high-quality single-nucleotide polymorphisms(SNPs)were screened,and 29 significant SNPs located on chromosomes A6,A7,A8,D1,D2,and D9,were repeatedly identified as associated with early-maturation traits,in at least two environments or two algorithms.Of these 29 significant SNPs,1,12,11,and 5 were related to SP,BP,FBP,and GP,respectively.Six peak SNPs,TM47967,TM13732,TM20937,TM28428,TM50283,and TM72552,exhibited phenotypic contributions of approximately 10%,which could allow them to be used for marker-assisted selection.One of these,TM72552,as well as four other SNPs,TM72554,TM72555,TM72558,and TM72559,corresponded to the quantitative trait loci previously reported.In total,274 candidate genes were identified from the genome sequences of upland cotton and were categorized based on their functional annotations.Finally,our studies identified Gh_D01G0340 and Gh_D01G0341 as potential candidate genes for improving cotton early maturity.     

Genetic dissection of lint yield and fiber quality traits of <Emphasis Type="Italic">G. hirsutum</Emphasis> in <Emphasis Type="Italic">G. barbadense</Emphasis> background

Gossypium hirsutum L. is a widely cultivated species characterized by its high yield and wide environmental adaptability, while Gossypium barbadense is well known for its superior fiber quality. In the present report, we, for the first time, developed G. hirsutum chromosome segment introgression lines (ILs) in a G. barbadense background (GhILs_Gb) and genetically dissected the inheritance of lint yield and fiber quality of G. hirsutum in G. barbadense background. The GhILs_Gb contains introgressed segments spanning 4121.20 cM, which represents 82.20% of the tetraploid cotton genome, with an average length of 18.65 cM. A total of 39 quantitative trait loci (QTLs) for six traits are identified in this IL population planted in Xinjiang. Four QTL clusters are detected. Of them, however, three clusters have deleterious effects on fiber length and strength and boll weight, and only one cluster on Chr. D9 can be used in marker-assisted selection (MAS) to increase lint percentage and decrease micronaire value in G. barbadense. QTL mapping showed that most of yield-related QTLs detected have positive effects and increase lint yield in G. barbadense, while most of fiber quality-related QTLs have deleterious effects except for micronaire. It suggested that G. hirsutum evolved to have a high lint yield. Several lines improved in lint percentage and boll size in G. barbadense by introgressed one fragment of G. hirsutum have been developed from the GhILs_Gb. The ILs developed, and the analyses presented here will enhance the understanding of the genetics of lint yield and fiber quality in G. hirsutum and facilitate further molecular breeding to improve lint yield in G. barbadense.     

Fine-mapping qFS07.1 controlling fiber strength in upland cotton (<Emphasis Type="Italic">Gossypium hirsutum</Emphasis> L.)

Key message

qFS07.1 controlling fiber strength was fine-mapped to a 62.6-kb region containing four annotated genes. RT-qPCR and sequence of candidate genes identified an LRR RLK gene as the most likely candidate.

Abstract

Fiber strength is an important component of cotton fiber quality and is associated with other properties, such as fiber maturity, fineness, and length. Stable QTL qFS07.1, controlling fiber strength, had been identified on chromosome 7 in an upland cotton recombinant inbred line (RIL) population from a cross (CCRI35?×?Yumian1) described in our previous studies. To fine-map qFS07.1, an F₂ population with 2484 individual plants from a cross between recombinant line RIL014 and CCRI35 was established. A total of 1518 SSR primer pairs, including 1062, designed from chromosome 1 of the Gossypium raimondii genome and 456 from chromosome 1 of the G. arboreum genome (corresponding to the QTL region) were used to fine-map qFS07.1, and qFS07.1 was mapped into a 62.6-kb genome region which contained four annotated genes on chromosome A07 of G. hirsutum. RT-qPCR and comparative analysis of candidate genes revealed a leucine-rich repeat protein kinase (LRR RLK) family protein to be a promising candidate gene for qFS07.1. Fine mapping and identification of the candidate gene for qFS07.1 will play a vital role in marker-assisted selection (MAS) and the study of mechanism of cotton fiber development.

     

Identification of candidate genes for key fibre‐related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre‐quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high‐density genetic map by specific‐locus amplified fragment sequencing (SLAF‐seq) to identify QTLs associated with fibre‐quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high‐yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre‐quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF‐based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co‐located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre‐related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre‐length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint‐percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre‐related traits and would be helpful for cloning fibre‐development‐related genes as well as for marker‐assisted genetic improvement in cotton.     

陆地棉主要产量相关性状的SSR标记关联分析

高产优质育种是我国棉花育种的主要目标。寻找与目标性状关联的分子标记,可克服常规育种的盲目性,提高分子标记辅助选择育种的准确性。本研究对118份陆地棉种质资源的衣分、单铃重、单株铃数及子指等4个产量相关性状进行2年2点的表型鉴定,并利用覆盖全基因组的、有多态性的214对SSR标记进行标记与性状的关联分析。结果表明:118份材料的4个产量相关性状表型变异丰富,平均变异系数的变幅在6.1%~19.1%之间,且在各环境中表现较为稳定;基因型分析表明,214对标记共检测到460个等位变异,基因多样性指数平均为0.5151,PIC值平均为0.4587,表明该批标记具有较多的等位变异数和较高的基因多样性;群体结构分析表明该批材料可分为4个亚群,且各类群中材料与地理来源无对应关系;关联分析结果显示,在显著条件下(-log10P1.3,P0.05),共有39个标记位点能够在2个及2个以上的环境中同时检测到,其中有4个标记位点同时与2个以上性状相关联,进一步比较发现,有7个位点与前人研究结果一致,其余32个位点为新发现的位点。研究结果可为陆地棉产量性状遗传改良的分子标记辅助选择提供理论依据。     

Quantitative trait loci mapping and genetic dissection for lint percentage in upland cotton (Gossypium hirsutum)

Lint percentage is an important character of cotton yield components and it is also correlated with cotton fibre development. In this study, we used a high lint percentage variety, Baimian1, and a low lint percentage, TM-1 genetic standard for Gossypium hirsutum, as parents to construct a mapping populations in upland cotton (G. hirsutum). A quantitative trait locus/loci (QTL) analysis of lint percentage was performed by using two mapping procedures; composite interval mapping (CIM), inclusive composite interval mapping (ICIM) and the F_2:3 populations in 2 years. Six main-effect QTL (M-QTL) for lint percentage (four significant and two suggestive) were detected in both years by CIM, and were located on chr. 3, chr. 19, chr. 26 and chr. 5/chr. 19. Of the six QTL, marker intervals and favourable gene sources of the significant M-QTL, qLP-3(2010) and qLP-3(2011) were consistent. These QTL were also detected by ICIM, and therefore, should preferentially be used for marker-assisted selection (MAS) of lint percentage. Another M-QTL, qLP-19(2010), was detected by two mapping procedures, and it could also be a candidate for MAS. We detected the interaction between two M-QTL and environment, and 11 epistatic QTL (E-QTL) and their interaction with environment by using ICIM. The study also found two EST-SSRs, NAU1187 and NAU1255, linked to M-QTL for lint percentage that could be candidate markers affecting cotton fibre development.     

Association between SSR markers and fibre traits in sea island cotton (<Emphasis Type="BoldItalic">Gossypium barbadense</Emphasis>) germplasm resources

Identification of molecular markers associated with fibre traits can accelerate cotton marker-assisted selection (MAS) programmes. In this study, Gossypium barbadense germplasm accessions with diverse origins (\(n = 123\)) were used to perform association analysis of fibre traits with 120 polymorphic simple sequence repeat (SSR) markers. In total, 120 polymorphic primer pairs amplified 258 loci with a mean of 2.15 loci per primer. Population structure analysis identified three main clusters for the accessions, which indicated agreement of genetic and predefined populations. Marker–trait associations (\(n= 58\)) were detected for 10 fibre traits with 26 SSR markers located on 15 chromosomes. The \(R^{2}\) (phenotypic variation explained) ranged from 3.19 to 15.21%. Two markers (NAU5465 and NAU3013) were found to be stably associated with boll number per plant (BNP) and fibre uniformity (UI), respectively. Four markers (BNL252, NAU3424, NAU3324 and CGR5202) associated with fibre quality traits preferentially clustered on the D8 chromosome, which was thus identified as an important candidate region for study molecular mechanisms underlying fibre quality and for use in breeding cotton cultivars for improving fibre quality. This study generated molecular data with a potential for better understanding of the genetic basis of the fibre traits and provided new markers for MAS in G. barbadense breeding programmes.     

Genetic diversity and population structure in the US Upland cotton (Gossypium hirsutum L.)

Key message

Genetic diversity and population structure in the US Upland cotton was established and core sets of allelic richness were identified for developing association mapping populations in cotton.

Abstract

Elite plant breeding programs could likely benefit from the unexploited standing genetic variation of obsolete cultivars without the yield drag typically associated with wild accessions. A set of 381 accessions comprising 378 Upland (Gossypium hirsutum L.) and 3 G. barbadense L. accessions of the United States cotton belt were genotyped using 120 genome-wide SSR markers to establish the genetic diversity and population structure in tetraploid cotton. These accessions represent more than 100 years of Upland cotton breeding in the United States. Genetic diversity analysis identified a total of 546 alleles across 141 marker loci. Twenty-two percent of the alleles in Upland accessions were unique, specific to a single accession. Population structure analysis revealed extensive admixture and identified five subgroups corresponding to Southeastern, Midsouth, Southwest, and Western zones of cotton growing areas in the United States, with the three accessions of G. barbadense forming a separate cluster. Phylogenetic analysis supported the subgroups identified by STRUCTURE. Average genetic distance between G. hirsutum accessions was 0.195 indicating low levels of genetic diversity in Upland cotton germplasm pool. The results from both population structure and phylogenetic analysis were in agreement with pedigree information, although there were a few exceptions. Further, core sets of different sizes representing different levels of allelic richness in Upland cotton were identified. Establishment of genetic diversity, population structure, and identification of core sets from this study could be useful for genetic and genomic analysis and systematic utilization of the standing genetic variation in Upland cotton.     

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

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

Genome-wide linkage mapping of yield-related traits in three Chinese bread wheat populations using high-density SNP markers

Key message

We identified 21 new and stable QTL, and 11 QTL clusters for yield-related traits in three bread wheat populations using the wheat 90 K SNP assay.

Abstract

Identification of quantitative trait loci (QTL) for yield-related traits and closely linked molecular markers is important in order to identify gene/QTL for marker-assisted selection (MAS) in wheat breeding. The objectives of the present study were to identify QTL for yield-related traits and dissect the relationships among different traits in three wheat recombinant inbred line (RIL) populations derived from crosses Doumai?×?Shi 4185 (D?×?S), Gaocheng 8901?×?Zhoumai 16 (G?×?Z) and Linmai 2?×?Zhong 892 (L?×?Z). Using the available high-density linkage maps previously constructed with the wheat 90 K iSelect single nucleotide polymorphism (SNP) array, 65, 46 and 53 QTL for 12 traits were identified in the three RIL populations, respectively. Among them, 34, 23 and 27 were likely to be new QTL. Eighteen common QTL were detected across two or three populations. Eleven QTL clusters harboring multiple QTL were detected in different populations, and the interval 15.5–32.3 cM around the Rht-B1 locus on chromosome 4BS harboring 20 QTL is an important region determining grain yield (GY). Thousand-kernel weight (TKW) is significantly affected by kernel width and plant height (PH), whereas flag leaf width can be used to select lines with large kernel number per spike. Eleven candidate genes were identified, including eight cloned genes for kernel, heading date (HD) and PH-related traits as well as predicted genes for TKW, spike length and HD. The closest SNP markers of stable QTL or QTL clusters can be used for MAS in wheat breeding using kompetitive allele-specific PCR or semi-thermal asymmetric reverse PCR assays for improvement of GY.

     

SSR marker-assisted improvement of fiber qualities in Gossypium hirsutum using G. barbadense introgression lines

Key message

This study demonstrates the first practical use of CSILs for the transfer of fiber quality QTLs into Upland cotton cultivars using SSR markers without detrimentally affecting desirable agronomic characteristics.

Abstract

Gossypium hirsutum is characterized by its high lint production and medium fiber quality compared to extra-long staple cotton G. barbadense. Transferring valuable traits or genes from G. barbadense into G. hirsutum is a promising but challenging approach through a traditional interspecific introgression strategy. We developed one set of chromosome segment introgression lines (CSILs), where TM-1, the genetic standard in G. hirsutum, was used as the recipient parent and the long staple cotton G. barbadense cv. Hai7124 was used as the donor parent by molecular marker-assisted selection (MAS). Among them, four CSILs, IL040-A4-1, IL080-D6-1, IL088-A7-3 and IL019-A2-6, found to be associated with superior fiber qualities including fiber length, strength and fineness QTL in Xinjiang were selected and backcrossed, and transferred these QTLs into three commercial Upland cotton cultivars such as Xinluzao (XLZ) 26, 41 and 42 grown in Xinjiang. By backcrossing and self-pollinating twice, five improved lines (3262-4, 3389-2, 3326-3, 3380-4 and 3426-5) were developed by MAS of background and introgressed segments. In diverse field trials, these QTLs consistently and significantly offered additive effects on the target phenotype. Furthermore, we also pyramided two segments from different CSILs (IL080-D6-1 and IL019-A2-6) into cultivar 0768 to accelerate breeding process purposefully with MAS. The improved lines pyramided by these two introgressed segments showed significant additive epistatic effects in four separate field trials. No significant alteration in yield components was observed in these modified lines. In summary, we first report that these CSILs have great potential to improve fiber qualities in Upland cotton MAS breeding programs.     

Construction of a complete set of alien chromosome addition lines from Gossypium australe in Gossypium hirsutum: morphological,cytological, and genotypic characterization

Key message

We report the first complete set of alien addition lines of G. hirsutum . The characterized lines can be used to introduce valuable traits from G. australe into cultivated cotton.

Abstract

Gossypium australe is a diploid wild cotton species (2n = 26, GG) native to Australia that possesses valuable characteristics unavailable in the cultivated cotton gene pool, such as delayed pigment gland morphogenesis in the seed and resistances to pests and diseases. However, it is very difficult to directly transfer favorable traits into cultivated cotton through conventional gene recombination due to the absence of pairing and crossover between chromosomes of G. australe and Gossypium hirsutum (2n = 52, AADD). To enhance the transfer of favorable genes from wild species into cultivated cotton, we developed a set of hirsutum–australe monosomic alien chromosome addition lines (MAAL) using a combination of morphological survey, microsatellite marker-assisted selection, and molecular cytogenetic analysis. The amphidiploid (2n = 78, AADDGG) of G. australe and G. hirsutum was consecutively backcrossed with upland cotton to develop alien addition lines of individual G. australe chromosomes in G. hirsutum. From these backcross progeny, we generated the first complete set of chromosome addition lines in cotton; 11 of 13 lines are monosomic additions, and chromosomes 7G^a and 13G^a are multiple additions. MAALs of 1G^a and 11G^a were the first to be isolated. The chromosome addition lines can be employed as bridges for the transfer of desired genes from G. australe into G. hirsutum, as well as for gene assignment, isolation of chromosome-specific probes, flow sorting and microdissection of chromosome, development of chromosome-specific ‘‘paints’’ for fluorochrome-labeled DNA fragments, physical mapping, and selective isolation and mapping of cDNAs for a particular G. australe chromosome.