Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

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 uncovers novel genomic regions and candidate genes of yield-related traits in upland cotton

Key message

A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield.

Abstract

Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

     

Genomic analyses reveal the genetic basis of early maturity and identification of loci and candidate genes in upland cotton (Gossypium hirsutum L.)

Although upland cotton (Gossypium hirsutism L.) originated in the tropics, this early maturity cotton can be planted as far north as 46°N in China due to the accumulation of numerous phenotypic and physiological adaptations during domestication. However, how the genome of early maturity cotton has been altered by strong human selection remains largely unknown. Herein, we report a cotton genome variation map generated by the resequencing of 436 cotton accessions. Whole‐genome scans for sweep regions identified 357 putative selection sweeps covering 4.94% (112 Mb) of the upland cotton genome, including 5184 genes. These genes were functionally related to flowering time control, hormone catabolism, ageing and defence response adaptations to environmental changes. A genome‐wide association study (GWAS) for seven early maturity traits identified 307 significant loci, 22.48% (69) of which overlapped with putative selection sweeps that occurred during the artificial selection of early maturity cotton. Several previously undescribed candidate genes associated with early maturity were identified by GWAS. This study provides insights into the genetic basis of early maturity in upland cotton as well as breeding resources for cotton improvement.     

Identification of loci governing eight agronomic traits using a GBS‐GWAS approach and validation by QTL mapping in soya bean

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome‐wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping‐by‐sequencing (GBS) approach was used to provide dense genome‐wide marker coverage (>47 000 SNPs) for a panel of 304 short‐season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.     

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.     

Genetic mapping and quantitative trait locus analysis of fiber quality traits using a three-parent composite population in upland cotton (Gossypium hirsutum L.)

Composite cross populations (CP) developed from three or more cultivars/lines are frequently used to improve agronomic and economic traits in crop cultivar development programs. Employing CP in linkage map construction and quantitative trait locus (QTL) mapping may increase the marker density of upland cotton (Gossypium hirsutum L.) genetic maps, exploit more adequate gene resources and facilitate marker-assisted selection (MAS). To construct a relatively high-density map and identify QTL associated with fiber quality traits in upland cotton, three elite upland cultivars/lines, Yumian 1, CRI 35 and 7,235, were used to obtain the segregating population, Yumian 1/CRI 35//Yumian 1/7,235. A genetic map containing 978 simple sequence repeat (SSR) loci and 69 linkage groups was constructed; the map spanned 4,184.4 cM, covering approximately 94.1% of the entire tetraploid cotton genome. A total of 63 QTL were detected, explaining 8.1–55.8% of the total phenotypic variance: 11 QTL for fiber elongation, 16 QTL for fiber length, 9 QTL for fiber micronaire reading, 10 QTL for fiber strength and 17 QTL for fiber length uniformity. The genetic map and QTL detected for fiber quality traits are promising for further breeding programs of upland cotton with improved fiber quality.     

Resequencing core accessions of a pedigree identifies derivation of genomic segments and key agronomic trait loci during cotton improvement

Upland cotton (Gossypium hirsutum) is the world's largest source of natural fibre and dominates the global textile industry. Hybrid cotton varieties exhibit strong heterosis that confers high fibre yields, yet the genome‐wide effects of artificial selection that have influenced Upland cotton during its breeding history are poorly understood. Here, we resequenced Upland cotton genomes and constructed a variation map of an intact breeding pedigree comprising seven elite and 19 backbone parents. Compared to wild accessions, the 26 pedigree accessions underwent strong artificial selection during domestication that has resulted in reduced genetic diversity but stronger linkage disequilibrium and higher extents of selective sweeps. In contrast to the backbone parents, the elite parents have acquired significantly improved agronomic traits, with an especially pronounced increase in the lint percentage. Notably, identify by descent (IBD) tracking revealed that the elite parents inherited abundant beneficial trait segments and loci from the backbone parents and our combined analyses led to the identification of a core genomic segment which was inherited in the elite lines from the parents Zhong 7263 and Ejing 1 and that was strongly associated with lint percentage. Additionally, SNP correlation analysis of this core segment showed that a non‐synonymous SNP (A‐to‐G) site in a gene encoding the cell wall‐associated receptor‐like kinase 3 (GhWAKL3) protein was highly correlated with increased lint percentage. Our results substantially increase the valuable genomics resources available for future genetic and functional genomics studies of cotton and reveal insights that will facilitate yield increases in the molecular breeding of cotton.     

A genome‐wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro‐ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome‐wide association study was conducted for these traits using 52 303 DArT‐seq markers to identify marker‐trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty‐five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN‐B1 and vrn‐B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker‐assisted selection in wheat.     

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.

     

Genome‐wide association analyses reveal the genetic basis of combining ability in rice

Combining ability is a measure for selecting elite parents and predicting hybrid performance in plant breeding. However, the genetic basis of combining ability remains unclear and a global view of combining ability from diverse mating designs is lacking. We developed a North Carolina II (NCII) population of 96 Oryza sativa and four male sterile lines to identify parents of greatest value for hybrid rice production. Statistical analyses indicated that general combining ability (GCA) and specific combining ability (SCA) contributed variously to different agronomic traits. In a genome‐wide association study (GWAS) of agronomic traits, GCA and SCA, we identified 34 significant associations (P < 2.39 × 10^?7). The superior alleles of GCA loci (Ghd8, GS3 and qSSR4) accumulated in parental lines with high GCA and explained 30.03% of GCA variance in grain yield, indicating that molecular breeding of high GCA parental lines is feasible. The distinct distributions of these QTLs contributed to the differentiation of parental GCA in subpopulations. GWAS of SCA identified 12 more loci that showed dominance on corresponding agronomic traits. We conclude that the accumulation of superior GCA and SCA alleles is an important contributor to heterosis and QTLs that greatly contributed to combining ability in our study would accelerate the identification of elite inbred lines and breeding of super hybrids.     

Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

Genome‐wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield‐limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome‐wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific‐locus amplified fragment sequencing (SLAF‐seq) approach. Trait–SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10^–5 (P = 1/85 630, –log₁₀P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P‐value (–log₁₀P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372‐kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR‐NBS‐LRR domains that may be involved in disease resistance response. Real‐time quantitative PCR and virus‐induced gene silencing (VIGS) analysis showed that CG02 was specific to up‐regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.     

Genome-wide association study-based identification genes influencing agronomic traits in rice (Oryza sativa L.)

Rice is one of the most important cereal crops, providing the daily dietary intake for approximately 50% of the global human population. Here, we re-sequenced 259 rice accessions, generating 1371.65 Gb of raw data. Furthermore, we performed genome-wide association studies (GWAS) on 13 agronomic traits using 2.8 million single nucleotide polymorphisms (SNPs) characterized in 259 rice accessions. Phenotypic data and best linear unbiased prediction (BLUP) values of each of the 13 traits over two years of each trait were used for the GWAS. The results showed that 816 SNP signals were significantly associated with the 13 agronomic traits. Then we detected candidate genes related to target traits within 200 kb upstream and downstream of the associated SNP loci, based on linkage disequilibrium (LD) blocks in the whole rice genome. These candidate genes were further identified through haplotype block constructions. This comprehensive study provides a timely and important genomic resource for breeding high yielding rice cultivars.     

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin–streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N₅₀ of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F₂ population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome (‘A’ or ‘D’) each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.     

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 detection of genetic loci and candidate genes for teat number and body conformation traits at birth in Chinese Sushan pigs

Body conformation at birth and teat number are economically important traits in the pig industry, as these traits are usually explored to evaluate the growth and reproductive potential of piglets. To detect genetic loci and candidate genes for these traits, we performed a GWAS on 269 pigs from a recently developed Chinese breed (Sushan) using 38  128 informative SNPs on the Affymetrix Porcine SNP 55K Array. In total, we detected one genome‐wide significant (P = 1.31e‐6) SNP for teat number on chromosome X and 15 chromosome‐wide significant SNPs for teat number, body weight, body length, chest circumference and cannon circumference at birth on chromosomes 1, 3, 4, 6, 7, 9, 10, 13, 14, 15, 17 and 18. The most significant SNP had an additive effect of 0.74 × total teat number, explaining 20% of phenotypic variance. Five significant SNPs resided in the previously reported quantitative trait loci for these traits and seven significant SNPs had a pleiotropic effect on multiple traits. Intriguingly, 12 of the genes nearest to the significant SNPs are functionally related to body conformation and teat number traits, including SPRED2, MKX, TMSB4X and ESR1. GO analysis revealed that candidate genes proximal to the significant SNPs were enriched in the G‐protein coupled receptor and steroid hormone‐mediated signaling pathway. Our findings shed light on the genetic basis of the measured traits and provide molecular markers especially for the genetic improvement of teat number in Sushan and related pigs.