Genome-Wide Association of Stem Water Soluble Carbohydrates in Bread Wheat

Water soluble carbohydrates (WSC) in stems play an important role in buffering grain yield in wheat against biotic and abiotic stresses; however, knowledge of genes controlling WSC is very limited. We conducted a genome-wide association study (GWAS) using a high-density 90K SNP array to better understand the genetic basis underlying WSC, and to explore marker-based breeding approaches. WSC was evaluated in an association panel comprising 166 Chinese bread wheat cultivars planted in four environments. Fifty two marker-trait associations (MTAs) distributed across 23 loci were identified for phenotypic best linear unbiased estimates (BLUEs), and 11 MTAs were identified in two or more environments. Liner regression showed a clear dependence of WSC BLUE scores on numbers of favorable (increasing WSC content) and unfavorable alleles (decreasing WSC), indicating that genotypes with higher numbers of favorable or lower numbers of unfavorable alleles had higher WSC content. In silico analysis of flanking sequences of trait-associated SNPs revealed eight candidate genes related to WSC content grouped into two categories based on the type of encoding proteins, namely, defense response proteins and proteins triggered by environmental stresses. The identified SNPs and candidate genes related to WSC provide opportunities for breeding higher WSC wheat cultivars.     

Association of a CONSTANS-LIKE gene to flowering and height in autotetraploid alfalfa

In alfalfa (Medicago sativa), an autotetraploid forage legume, stem length is a major component of forage yield, quality and competing ability. In this species, flowering date is not a breeding criterion. Association mapping based on a candidate gene approach has given good results in plants, including autotetraploid species for which genetic analyses are complex. The role of a CONSTANS-LIKE gene, identified as a candidate for stem elongation and flowering date in the model legume M. truncatula, was tested for association with the same traits in alfalfa. Four hundred genotypes from ten cultivars were evaluated for stem height and flowering date in two locations during 4 years. They were genotyped with simple sequence repeat markers and a low structuration was noticed. Primers were designed to amplify and sequence two regions of the alfalfa gene homologous to CONSTANS-LIKE. Single nucleotide polymorphisms (SNPs) were detected and their allelic dose in each genotype was scored. Linkage disequilibrium within CONSTANS-LIKE rapidly decreased as expected. Eight SNPs with a frequency above 10% were detected over 1,010 bp (one SNP every 126 bp on average) in the 400 genotypes. This number was lower than observed in a neutral gene (a SNP every 31 bp on average). Highly significant associations of three SNPs to flowering date and stem height were identified. Each SNP explained up to 4.2% of the genetic variance. Thus, as in the model species, the CONSTANS-LIKE gene was shown to be involved in flowering date and stem height in alfalfa.     

Confirmation and fine‐mapping of clinical mastitis and somatic cell score QTL in Nordic Holstein cattle

A genome‐wide association study of 2098 progeny‐tested Nordic Holstein bulls genotyped for 36 387 SNPs on 29 autosomes was conducted to confirm and fine‐map quantitative trait loci (QTL) for mastitis traits identified earlier using linkage analysis with sparse microsatellite markers in the same population. We used linear mixed model analysis where a polygenic genetic effect was fitted as a random effect and single SNPs were successively included as fixed effects in the model. We detected 143 SNP‐by‐trait significant associations (P < 0.0001) on 20 chromosomes affecting mastitis‐related traits. Among them, 21 SNP‐by‐trait combinations exceeded the genome‐wide significant threshold. For 12 chromosomes, both the present association study and the previous linkage study detected QTL, and of these, six were in the same chromosomal locations. Strong associations of SNPs with mastitis traits were observed on bovine autosomes 6, 13, 14 and 20. Possible candidate genes for these QTL were identified. Identification of SNPs in linkage disequilibrium with QTL will enable marker‐based selection for mastitis resistance. The candidate genes identified should be further studied to detect candidate polymorphisms underlying these QTL.     

Allelic variation in the perennial ryegrass FLOWERING LOCUS T gene is associated with changes in flowering time across a range of populations

The Arabidopsis (Arabidopsis thaliana) FLOWERING LOCUS T (FT) gene and its orthologs in other plant species (e.g. rice [Oryza sativa] OsFTL2/Hd3a) have an established role in the photoperiodic induction of flowering response. The genomic and phenotypic variations associated with the perennial ryegrass (Lolium perenne) ortholog of FT, designated LpFT3, was assessed in a diverse collection of nine European germplasm populations, which together constituted an association panel of 864 plants. Sequencing and genotyping of a series of amplicons derived from the nine populations, containing the complete exon and intron sequences as well as 5' and 3' noncoding sequences of LpFT3, identified a total of seven haplotypes. Genotyping assays designed to detect the genomic variation showed that three haplotypes were present in approximately equal proportions and represented 84% of the total, with a fourth representing a further 11%. Of the three major haplotypes, two were predicted to code for identical protein products and the third contained two amino acid substitutions. Association analysis using either a mixed model with a relationship matrix to correct for population structure and relatedness or structured association with further correction using genomic control indicated significant associations between LpFT3 and variation in flowering time. These associations were corroborated in a validation population segregating for the same major alleles. The most "diagnostic" region of genomic variation was situated 5' of the coding sequence. Analysis of this region identified that the interhaplotype variation was closely associated with sequence motifs that were apparently conserved in the 5' region of orthologs of LpFT3 from other plant species. These may represent cis-regulatory elements involved in influencing the expression of this gene.     

Validation of growth-related quantitative trait loci markers in different Exopalaemon carinicauda families for marker-assisted selection

We detected growth-related QTL and associated markers from the backcross population of Exopalaemon carinicauda in the previous study. Based on our previous study, the 47 SNP markers associated with candidate growth trait QTL were selected to analyze the association between these markers and body weight (BW), body length and abdominal segment length traits in four different populations including wild population, a full-sib family, a half-sib family and a backcross population for evaluating their potential application of marker-assisted selection in E. carinicauda. The general linear model (GLM) and mixed linear model were applied and the associations between SNP loci and three growth-related traits verified. The results showed that the Marker79268 and Marker100644 were significantly associated with the BW trait in more than three populations by the GLM method. The Marker100644 was significantly associated with BW in the full-sib family, half-sib family and backcross populations by the GLM and mixed linear model methods. Our findings will provide useful SNP markers to go forward to improve growth performance through more refined marker-assisted selection in E. carinicauda.     

Effects of Favorable Alleles for Water-Soluble Carbohydrates at Grain Filling on Grain Weight under Drought and Heat Stresses in Wheat

Drought, heat and other abiotic stresses during grain filling can result in reductions in grain weight. Conserved water-soluble carbohydrates (WSC) at early grain filling play an important role in partial compensation of reduced carbon supply. A diverse population of 262 historical winter wheat accessions was used in the present study. There were significant correlations between 1000-grain weight (TGW) and four types of WSC, viz. (1) total WSC at the mid-grain filling stage (14 days after flowering) produced by leaves and non-leaf organs; (2) WSC contributed by current leaf assimilation during the mid-grain filling; (3) WSC in non-leaf organs at the mid-grain filling, excluding the current leaf assimilation; and (4) WSC used for respiration and remobilization during the mid-grain filling. Association and favorable allele analyses of 209 genome-wide SSR markers and the four types of WSC were conducted using a mixed linear model. Seven novel favorable WSC alleles exhibited positive individual contributions to TGW, which were verified under 16 environments. Dosage effects of pyramided favorable WSC alleles and significantly linear correlations between the number of favorable WSC alleles and TGW were observed. Our results suggested that pyramiding more favorable WSC alleles was effective for improving both WSC and grain weight in future wheat breeding programs.     

Breeding for Bio-ethanol Production in Lolium perenne L.: Association of Allelic Variation with High Water-Soluble Carbohydrate Content

Increasing the extractable sugar yield from perennial crops is one strategy to generate renewable fuels such as bio-ethanol. Lolium perenne L. (perennial ryegrass) can contain significant (>30% dry matter) water-soluble sugars in the form of polymeric fructan which is readily extracted, broken down and fermented to bio-ethanol. A population of L. perenne generated from four parents which differed in water-soluble carbohydrate (WSC) content was subjected to multiple rounds of selection and recombination on the basis of early spring WSC content to produce a high WSC, and a low WSC population. A control population was generated by selecting the same number of plants at random. The alleles present at six candidate gene loci were analysed before and after selection and correlated to WSC content. Significant differences in the allele frequency of L. perenne soluble-acid invertase1:4 were observed between the three populations with one haplotype significantly associated with the high WSC C2^S+ population (after three rounds of selection and two rounds of recombination). Moreover, WSC content was also associated with biomass accumulation. Thus, in addition to a 2.84-fold increase in WSC yield, the C2^S+ population also had 1.48-fold more biomass per plant, resulting in 3.9-fold higher WSC yield per plant than the control population.     

Favorable Alleles for Stem Water-Soluble Carbohydrates Identified by Association Analysis Contribute to Grain Weight under Drought Stress Conditions in Wheat

Drought is a major environmental constraint to crop distribution and productivity. Stem water-soluble carbohydrates (WSC) buffer wheat grain yield against conditions unfavorable for photosynthesis during the grain filling stage. In this study, 262 winter wheat accessions and 209 genome-wide SSR markers were collected and used to undertake association analysis based on a mixed linear model (MLM). The WSC in different internodes at three growth stages and 1000-grain weight (TGW) were investigated under four environmental regimes (well-watered, drought stress during the whole growth period, and two levels of terminal drought stress imposed by chemical desiccation under the well-watered and drought stress during the whole growth period conditions). Under diverse drought stress conditions, WSC in lower internodes showed significant positive correlations with TGW, especially at the flowering stage under well-watered conditions and at grain filling under drought stress. Sixteen novel WSC-favorable alleles were identified, and five of them contributed to significantly higher TGW. In addition, pyramiding WSC favorable alleles was not only effective for obtaining accessions with higher WSC, but also for enhancing TGW under different water regimes. During the past fifty years of wheat breeding, WSC was selected incidentally. The average number of favorable WSC alleles increased from 1.13 in the pre-1960 period to 4.41 in the post-2000 period. The results indicate a high potential for using marker-assisted selection to pyramid WSC favorable alleles in improving WSC and TGW in wheat.     

Extremely high cytoplasmic diversity in natural and breeding populations of Lolium (Poaceae)

Ten chloroplast microsatellite markers were used to characterise chloroplast genetic diversity at allelic and haplotypic level in 104 accessions of Lolium perenne, other Lolium species, Festuca species and x Festulolium cultivars. Furthermore, genetic relationships between the accessions and biogeographic distribution of haplotypes were investigated using a range of Nei's population genetic diversity measures and analysis of molecular variance (AMOVA). An extremely high number (511) of haplotypes was detected in 1575 individuals. Nei's gene diversity values among L. perenne accessions ranged between 0 and 0.333. Much of the L. perenne European ecotype diversity (61%), as calculated using AMOVA, could be attributed to within-population variance and this is likely caused by, and maintained by, high levels of natural and anthropogenic seed dispersal. Plastid gene pools and maternal lineages for L. perenne could be clearly identified. Evidence was found, using AMOVA, to show a likely migration route of L. perenne from Southern regions of Europe northwards.     

Identifying genetic components controlling fertility in the outcrossing grass species perennial ryegrass (Lolium perenne) by quantitative trait loci analysis and comparative genetics

Mutational load and resource allocation factors and their effects on limiting seed set were investigated in ryegrass by comparative mapping genomics and quantitative trait loci (QTL) analysis in two perennial ryegrass (Lolium perenne) mapping families sharing common genetic markers. Quantitative trait loci for seed-set were identified on chromosome (LG) 7 in both families and on LG4 of the F2/WSC family. On LG7, seed-set and heading date QTLs colocalized in both families and cannot be unequivocally resolved. Comparative genomics suggests that the LG7 region is syntenous to a region of rice LG6 which contains both fertility (S5(n)) and heading date (Hd1, Hd3a) candidate genes. The LG4 region is syntenous to a region of rice LG3 which contains a fertility (S33) candidate gene. QTL maxima for seed-set and heading date on LG4 in the F2/WSC family are separated by c. 8 cm, indicating distinct genetic control. Low seed set is under the control of recessive genes at both LG4 and LG7 locations. The identification of QTLs associated with seed set, a major component of seed yield in perennial ryegrass, indicates that mutational load associated with these genomic regions can be mitigated through marker-assisted selection.     

Genetic association mapping identifies single nucleotide polymorphisms in genes that affect abscisic acid levels in maize floral tissues during drought

In maize, water stress at flowering causes loss of kernel set and productivity. While changes in the levels of sugars and abscisic acid (ABA) are thought to play a role in this stress response, the mechanistic basis and genes involved are not known. A candidate gene approach was used with association mapping to identify loci involved in accumulation of carbohydrates and ABA metabolites during stress. A panel of single nucleotide polymorphisms (SNPs) in genes from these metabolic pathways and in genes for reproductive development and stress response was used to genotype 350 tropical and subtropical maize inbred lines that were well watered or water stressed at flowering. Pre-pollination ears, silks, and leaves were analysed for sugars, starch, proline, ABA, ABA-glucose ester, and phaseic acid. ABA and sugar levels in silks and ears were negatively correlated with their growth. Association mapping with 1229 SNPs in 540 candidate genes identified an SNP in the maize homologue of the Arabidopsis MADS-box gene, PISTILLATA, which was significantly associated with phaseic acid in ears of well-watered plants, and an SNP in pyruvate dehydrogenase kinase, a key regulator of carbon flux into respiration, that was associated with silk sugar concentration. An SNP in an aldehyde oxidase gene was significantly associated with ABA levels in silks of water-stressed plants. Given the short range over which decay of linkage disequilibrium occurs in maize, the results indicate that allelic variation in these genes affects ABA and carbohydrate metabolism in floral tissues during drought.     

Association between Expression Quantitative Trait Loci and Metabolic Traits in Two Korean Populations

Most genome-wide association studies consider genes that are located closest to single nucleotide polymorphisms (SNPs) that are highly significant for those studies. However, the significance of the associations between SNPs and candidate genes has not been fully determined. An alternative approach that used SNPs in expression quantitative trait loci (eQTL) was reported previously for Crohn’s disease; it was shown that eQTL-based preselection for follow-up studies was a useful approach for identifying risk loci from the results of moderately sized GWAS. In this study, we propose an approach that uses eQTL SNPs to support the functional relationships between an SNP and a candidate gene in a genome-wide association study. The genome-wide SNP genotypes and 10 biochemical measures (fasting glucose levels, BUN, serum albumin levels, AST, ALT, gamma GTP, total cholesterol, HDL cholesterol, triglycerides, and LDL cholesterol) were obtained from the Korean Association Resource (KARE) consortium. The eQTL SNPs were isolated from the SNP dataset based on the RegulomeDB eQTL-SNP data from the ENCODE projects and two recent eQTL reports. A total of 25,658 eQTL SNPs were tested for their association with the 10 metabolic traits in 2 Korean populations (Ansung and Ansan). The proportion of phenotypic variance explained by eQTL and non-eQTL SNPs showed that eQTL SNPs were more likely to be associated with the metabolic traits genetically compared with non-eQTL SNPs. Finally, via a meta-analysis of the two Korean populations, we identified 14 eQTL SNPs that were significantly associated with metabolic traits. These results suggest that our approach can be expanded to other genome-wide association studies.     

Haplotype analysis improved evidence for candidate genes for intramuscular fat percentage from a genome wide association study of cattle

In genome wide association studies (GWAS), haplotype analyses of SNP data are neglected in favour of single point analysis of associations. In a recent GWAS, we found that none of the known candidate genes for intramuscular fat (IMF) had been identified. In this study, data from the GWAS for these candidate genes were re-analysed as haplotypes. First, we confirmed that the methodology would find evidence for association between haplotypes in candidate genes of the calpain-calpastatin complex and musculus longissimus lumborum peak force (LLPF), because these genes had been confirmed through single point analysis in the GWAS. Then, for intramuscular fat percent (IMF), we found significant partial haplotype substitution effects for the genes ADIPOQ and CXCR4, as well as suggestive associations to the genes CEBPA, FASN, and CAPN1. Haplotypes for these genes explained 80% more of the phenotypic variance compared to the best single SNP. For some genes the analyses suggested that there was more than one causative mutation in some genes, or confirmed that some causative mutations are limited to particular subgroups of a species. Fitting the SNPs and their interactions simultaneously explained a similar amount of the phenotypic variance compared to haplotype analyses. Haplotype analysis is a neglected part of the suite of tools used to analyse GWAS data, would be a useful method to extract more information from these data sets, and may contribute to reducing the missing heritability problem.     

The SNP c.1326T>G in the non-SMC condensin I complex, subunit G (NCAPG) gene encoding a p.Ile442Met variant is associated with an increase in body frame size at puberty in cattle

Recently, we had located a bovine carcass weight QTL, CW-2, to a 591-kb interval on BTA6 and have identified the SNP c.1326T>G in the NCAPG (non-SMC condensin I complex, subunit G) gene that leads to the amino acid change p.Ile442Met in the NCAPG protein, which is a candidate causative variation. Here, we examined the association of the NCAPG:c.1326T>G locus with linear skeletal measurements of growth-associated traits during adolescence, which is a period of intensive growth, using two historically and geographically distant cattle populations: 792 Japanese Black steers and 161 F(2) bulls of an experimental cross from Charolais and German Holstein. In both populations, the SNP NCAPG:c.1326T>G was associated with each component of body frame size: height, length and width at puberty. The associations of CW-2 with height- and length-associated traits were observed at an earlier growth period compared to the associations with thickness- and width-associated traits, indicating that the primary effect of the CW-2 QTL may possibly be exerted on skeletal growth. The significant associations of the NCAPG:c.1326T>G locus with growth-associated skeletal measurements are similar to the effects of the syntenic region on human chromosome 4 that are associated with adult height in humans, supporting the hypothesis that CW-2 is analogous to the human locus and pointing to a conserved growth-associated locus or chromosomal region present in both species.     

Identification of favorable SNP alleles and candidate genes responsible for inflorescence-related traits via GWAS in chrysanthemum

Key message

81 SNPs were identified for three inflorescence-related traits, in which 15 were highly favorable. Two dCAPS markers were developed for future MAS breeding, and six candidate genes were predicted.

Abstract

Chrysanthemum is a leading ornamental species worldwide and demonstrates a wealth of morphological variation. Knowledge about the genetic basis of its phenotypic variation for key horticultural traits can contribute to its effective management and genetic improvement. In this study, we conducted a genome-wide association study (GWAS) based on two years of phenotype data and a set of 92,617 single nucleotide polymorphisms (SNPs) using a panel of 107 diverse cut chrysanthemums to dissect the genetic control of three inflorescence-related traits. A total of 81 SNPs were significantly associated with the three inflorescence-related traits (capitulum diameter, number of ray florets and flowering time) in at least one environment, with an individual allele explaining 22.72–38.67% of the phenotypic variation. Fifteen highly favorable alleles were identified for the three target traits by computing the phenotypic effect values for the stable associations detected in 2 year-long trials at each locus. Dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between highly favorable allele numbers and corresponding phenotypic performance were observed. Two highly favorable SNP alleles correlating to flowering time and capitulum diameter were converted to derived cleaved amplified polymorphic sequence (dCAPS) markers to facilitate future breeding. Finally, six putative candidate genes were identified that contribute to flowering time and capitulum diameter. These results serve as a foundation for analyzing the genetic mechanisms underlying important horticultural traits and provide valuable insights into molecular marker-assisted selection (MAS) in chrysanthemum breeding programs.

     

QTL-seq identifies an early flowering QTL located near Flowering Locus T in cucumber

Key message

Next-generation sequencing enabled a fast discovery of a major QTL controlling early flowering in cucumber, corresponding to the FT gene conditioning flowering time in Arabidopsis.

Abstract

Next-generation sequencing technologies are making it faster and more efficient to establish the association of agronomic traits with molecular markers or candidate genes, which is the requirement for marker-assisted selection in molecular breeding. Early flowering is an important agronomic trait in cucumber (Cucumis sativus L.), but the underlying genetic mechanism is unknown. In this study, we identified a candidate gene for early flowering QTL, Ef1.1 through QTL-seq. Segregation analysis in F₂ and BC₁ populations derived from a cross between two inbred lines “Muromskij” (early flowering) and “9930” (late flowering) suggested quantitative nature of flowering time in cucumber. Genome-wide comparison of SNP profiles between the early and late-flowering bulks constructed from F₂ plants identified a major QTL, designated Ef1.1 on cucumber chromosome 1 for early flowering in Muromskij, which was confirmed by microsatellite marker-based classical QTL mapping in the F₂ population. Joint QTL-seq and traditional QTL analysis delimited Ef1.1 to an 890 kb genomic region. A cucumber gene, Csa1G651710, was identified in this region, which is a homolog of the FLOWERING LOCUS T (FT), the main flowering switch gene in Arabidopsis. Quantitative RT-PCR study of the expression level of Csa1G651710 revealed significantly higher expression in early flowering genotypes. Data presented here provide support for Csa1G651710 as a possible candidate gene for early flowering in the cucumber line Muromskij.     

Evaluation of SNPs in the chicken HMGA2 gene as markers for body weight gain

A QTL affecting body weight in chickens has been mapped to GGA1, between the markers GCT0006 and MCW0106. The gene HMGA2, which was previously identified as a candidate gene for determining body height in humans and mice, is also conspicuously close to the MCW0106 marker in chickens. Subsequently, 14 SNP markers of HMGA2 were genotyped in CAU chicken resource populations, and the associations between body weight and those SNP markers that displayed polymorphisms were analysed. Three SNPs (rs13849241, rs15231472 and rs13849381) were found to be significantly correlated with body weight in chickens (P < 0.05). Furthermore, haplotypes constructed based on these three SNPs were also discovered to be associated with body weight in chickens at the ages of 6, 7, 9 and 12 weeks. These results suggest that the chicken HMGA2 gene is indeed involved in body weight gain.     

Tree scanning: a method for using haplotype trees in phenotype/genotype association studies

We use evolutionary trees of haplotypes to study phenotypic associations by exhaustively examining all possible biallelic partitions of the tree, a technique we call tree scanning. If the first scan detects significant associations, additional rounds of tree scanning are used to partition the tree into three or more allelic classes. Two worked examples are presented. The first is a reanalysis of associations between haplotypes at the Alcohol Dehydrogenase locus in Drosophila melanogaster that was previously analyzed using a nested clade analysis, a more complicated technique for using haplotype trees to detect phenotypic associations. Tree scanning and the nested clade analysis yield the same inferences when permutation testing is used with both approaches. The second example is an analysis of associations between variation in various lipid traits and genetic variation at the Apolipoprotein E (APOE) gene in three human populations. Tree scanning successfully identified phenotypic associations expected from previous analyses. Tree scanning for the most part detected more associations and provided a better biological interpretative framework than single SNP analyses. We also show how prior information can be incorporated into the tree scan by starting with the traditional three electrophoretic alleles at APOE. Tree scanning detected genetically determined phenotypic heterogeneity within all three electrophoretic allelic classes. Overall, tree scanning is a simple, powerful, and flexible method for using haplotype trees to detect phenotype/genotype associations at candidate loci.     

Identification of perennial ryegrass (Lolium perenne (L.)) and meadow fescue (Festuca pratensis (Huds.)) candidate orthologous sequences to the rice Hd1(Se1) and barley HvCO1 CONSTANS-like genes through comparative mapping and microsynteny

Microsynteny with rice and comparative genetic mapping were used to identify candidate orthologous sequences to the rice Hd1(Se1) gene in Lolium perenne and Festuca pratensis. A F. pratensis bacterial artificial chromosome (BAC) library was screened with a marker (S2539) physically close to Hd1 in rice to identify the equivalent genomic region in F. pratensis. The BAC sequence was used to identify and map the same region in L. perenne. Predicted protein sequences for L. perenne and F. pratensis Hd1 candidates (LpHd1 and FpHd1) indicated they were CONSTANS-like zinc finger proteins with 61-62% sequence identity with rice Hd1 and 72% identity with barley HvCO1. LpHd1 and FpHd1 were physically linked in their respective genomes (< 4 kb) to marker S2539, which was mapped to L. perenne chromosome 7. The identified candidate orthologues of rice Hd1 and barley HvCO1 in L. perenne and F. pratensis map to chromosome 7, a region of the L. perenne genome which has a degree of conserved genetic synteny both with rice chromosome 6, which contains Hd1, and barley chromosome 7H, which contains HvCO1.