Genome-wide association and functional follow-up reveals new loci for kidney function

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.     

Genome-wide QTL and bulked transcriptomic analysis reveals new candidate genes for the control of tuber carotenoid content in potato (Solanum tuberosum L.)

Key message

Genome-wide QTL analysis of potato tuber carotenoid content was investigated in populations of Solanum tuberosum Group Phureja that segregate for flesh colour, revealing a novel major QTL on chromosome 9.

Abstract

The carotenoid content of edible plant storage organs is a key nutritional and quality trait. Although the structural genes that encode the biosynthetic enzymes are well characterised, much less is known about the factors that determine overall storage organ content. In this study, genome-wide QTL mapping, in concert with an efficient ‘genetical genomics’ analysis using bulked samples, has been employed to investigate the genetic architecture of potato tuber carotenoid content. Two diploid populations of Solanum tuberosum Group Phureja were genotyped (AFLP, SSR and DArT markers) and analysed for their tuber carotenoid content over two growing seasons. Common to both populations were QTL that explained relatively small proportions of the variation in constituent carotenoids and a major QTL on chromosome 3 explaining up to 71 % of the variation in carotenoid content. In one of the populations (01H15), a second major carotenoid QTL was identified on chromosome 9, explaining up to 20 % of the phenotypic variation. Whereas the major chromosome 3 QTL was likely to be due to an allele of a gene encoding β-carotene hydroxylase, no known carotenoid biosynthetic genes are located in the vicinity of the chromosome 9 QTL. A unique expression profiling strategy using phenotypically distinct bulks comprised individuals with similar carotenoid content provided further support for the QTL mapping to chromosome 9. This study shows the potential of using the potato genome sequence to link genetic maps to data arising from eQTL approaches to enhance the discovery of candidate genes underlying QTLs.     

Genome-wide association analysis reveals putative Alzheimer's disease susceptibility loci in addition to APOE

Alzheimer's disease (AD) is a genetically complex and heterogeneous disorder. To date four genes have been established to either cause early-onset autosomal-dominant AD (APP, PSEN1, and PSEN2(1-4)) or to increase susceptibility for late-onset AD (APOE5). However, the heritability of late-onset AD is as high as 80%, (6) and much of the phenotypic variance remains unexplained to date. We performed a genome-wide association (GWA) analysis using 484,522 single-nucleotide polymorphisms (SNPs) on a large (1,376 samples from 410 families) sample of AD families of self-reported European descent. We identified five SNPs showing either significant or marginally significant genome-wide association with a multivariate phenotype combining affection status and onset age. One of these signals (p = 5.7 x 10(-14)) was elicited by SNP rs4420638 and probably reflects APOE-epsilon4, which maps 11 kb proximal (r2 = 0.78). The other four signals were tested in three additional independent AD family samples composed of nearly 2700 individuals from almost 900 families. Two of these SNPs showed significant association in the replication samples (combined p values 0.007 and 0.00002). The SNP (rs11159647, on chromosome 14q31) with the strongest association signal also showed evidence of association with the same allele in GWA data generated in an independent sample of approximately 1,400 AD cases and controls (p = 0.04). Although the precise identity of the underlying locus(i) remains elusive, our study provides compelling evidence for the existence of at least one previously undescribed AD gene that, like APOE-epsilon4, primarily acts as a modifier of onset age.     

Genetic and physiological analysis of iron biofortification in maize kernels

Background

Maize is a major cereal crop widely consumed in developing countries, which have a high prevalence of iron (Fe) deficiency anemia. The major cause of Fe deficiency in these countries is inadequate intake of bioavailable Fe, where poverty is a major factor. Therefore, biofortification of maize by increasing Fe concentration and or bioavailability has great potential to alleviate this deficiency. Maize is also a model system for genomic research and thus allows the opportunity for gene discovery. Here we describe an integrated genetic and physiological analysis of Fe nutrition in maize kernels, to identify loci that influence grain Fe concentration and bioavailability.

Methodology

Quantitative trait locus (QTL) analysis was used to dissect grain Fe concentration (FeGC) and Fe bioavailability (FeGB) from the Intermated B73 × Mo17 (IBM) recombinant inbred (RI) population. FeGC was determined by ion coupled argon plasma emission spectroscopy (ICP). FeGB was determined by an in vitro digestion/Caco-2 cell line bioassay.

Conclusions

Three modest QTL for FeGC were detected, in spite of high heritability. This suggests that FeGC is controlled by many small QTL, which may make it a challenging trait to improve by marker assisted breeding. Ten QTL for FeGB were identified and explained 54% of the variance observed in samples from a single year/location. Three of the largest FeGB QTL were isolated in sister derived lines and their effect was observed in three subsequent seasons in New York. Single season evaluations were also made at six other sites around North America, suggesting the enhancement of FeGB was not specific to our farm site. FeGB was not correlated with FeGC or phytic acid, suggesting that novel regulators of Fe nutrition are responsible for the differences observed. Our results indicate that iron biofortification of maize grain is achievable using specialized phenotyping tools and conventional plant breeding techniques.     

Genome-wide interaction-based association analysis identified multiple new susceptibility Loci for common diseases

Genome-wide interaction-based association (GWIBA) analysis has the potential to identify novel susceptibility loci. These interaction effects could be missed with the prevailing approaches in genome-wide association studies (GWAS). However, no convincing loci have been discovered exclusively from GWIBA methods, and the intensive computation involved is a major barrier for application. Here, we developed a fast, multi-thread/parallel program named “pair-wise interaction-based association mapping” (PIAM) for exhaustive two-locus searches. With this program, we performed a complete GWIBA analysis on seven diseases with stringent control for false positives, and we validated the results for three of these diseases. We identified one pair-wise interaction between a previously identified locus, C1orf106, and one new locus, TEC, that was specific for Crohn''s disease, with a Bonferroni corrected P<0.05 (P = 0.039). This interaction was replicated with a pair of proxy linked loci (P = 0.013) on an independent dataset. Five other interactions had corrected P<0.5. We identified the allelic effect of a locus close to SLC7A13 for coronary artery disease. This was replicated with a linked locus on an independent dataset (P = 1.09×10⁻⁷). Through a local validation analysis that evaluated association signals, rather than locus-based associations, we found that several other regions showed association/interaction signals with nominal P<0.05. In conclusion, this study demonstrated that the GWIBA approach was successful for identifying novel loci, and the results provide new insights into the genetic architecture of common diseases. In addition, our PIAM program was capable of handling very large GWAS datasets that are likely to be produced in the future.     

Genome-wide association analysis of seedling root development in maize (Zea mays L.)

Background

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

Results

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10^-7). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

Conclusion

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1226-9) contains supplementary material, which is available to authorized users.     

Quantitative trait loci for biofortification traits in maize grain

Detecting genes that influence biofortification traits in cereal grain could help increase the concentrations of bioavailable mineral elements in crops to solve the global mineral malnutrition problem. The aims of this study were to detect the quantitative trait loci (QTLs) for phosphorus (P), iron (Fe), zinc (Zn), and magnesium (Mg) concentrations in maize grain in a mapping population, as well as QTLs for bioavailable Fe, Zn, and Mg, by precalculating their respective ratios with P. Elemental analysis of grain samples was done by coupled plasma-optical emission spectrometry in 294 F(4) lines of a biparental population taken from field trials of over 3 years. The population was mapped using sets of 121 polymorphic markers. QTL analysis revealed 32 significant QTLs detected for 7 traits, of which some were colocalized. The Additive-dominant model revealed highly significant additive effects, suggesting that biofortification traits in maize are generally controlled by numerous small-effect QTLs. Three QTLs for Fe/P, Zn/P, and Mg/P were colocalized on chromosome 3, coinciding with simple sequence repeats marker bnlg1456, which resides in close proximity to previously identified phytase genes (ZM phys1 and phys2). Thus, we recommend the ratios as bioavailability traits in biofortification research.     

Genome-wide association analysis for feed efficiency in Angus cattle

Estimated breeding values for average daily feed intake (AFI; kg/day), residual feed intake (RFI; kg/day) and average daily gain (ADG; kg/day) were generated using a mixed linear model incorporating genomic relationships for 698 Angus steers genotyped with the Illumina BovineSNP50 assay. Association analyses of estimated breeding values (EBVs) were performed for 41,028 single nucleotide polymorphisms (SNPs), and permutation analysis was used to empirically establish the genome-wide significance threshold (P < 0.05) for each trait. SNPs significantly associated with each trait were used in a forward selection algorithm to identify genomic regions putatively harbouring genes with effects on each trait. A total of 53, 66 and 68 SNPs explained 54.12% (24.10%), 62.69% (29.85%) and 55.13% (26.54%) of the additive genetic variation (when accounting for the genomic relationships) in steer breeding values for AFI, RFI and ADG, respectively, within this population. Evaluation by pathway analysis revealed that many of these SNPs are in genomic regions that harbour genes with metabolic functions. The presence of genetic correlations between traits resulted in 13.2% of SNPs selected for AFI and 4.5% of SNPs selected for RFI also being selected for ADG in the analysis of breeding values. While our study identifies panels of SNPs significant for efficiency traits in our population, validation of all SNPs in independent populations will be necessary before commercialization.     

Genome-wide association study reveals genetic architecture of coleoptile length in wheat

Key message

Eight QTL for coleoptile length were identified in a genome-wide association study on a set of 893 wheat accessions, four of which are novel loci.

Abstract

Wheat cultivars with long coleoptiles are preferred in wheat-growing regions where deep planting is practiced. However, the wide use of gibberellic acid (GA)-insensitive dwarfing genes, Rht-B1b and Rht-D1b, makes it challenging to breed dwarf wheat cultivars with long coleoptiles. To understand the genetic basis of coleoptile length, we performed a genome-wide association study on a set of 893 landraces and historical cultivars using 5011 single nucleotide polymorphism (SNP) markers. Structure analysis revealed four subgroups in the association panel. Association analysis results suggested that Rht-B1b and Rht-D1b genes significantly reduced coleoptile length, and eight additional quantitative trait loci (QTL) for coleoptile length were also identified. These QTL explained 1.45–3.18 and 1.36–3.11% of the phenotypic variation in 2015 and 2016, respectively, and their allelic substitution effects ranged from 0.31 to 1.75 cm in 2015, and 0.63–1.55 cm in 2016. Of the eight QTL, QCL.stars-1BS1, QCL.stars-2DS1, QCL.stars-4BS2, and QCL.stars-5BL1 are likely novel loci for coleoptile length. The favorable alleles in each accession ranged from two to eight with an average of 5.8 at eight loci in the panel, and more favorable alleles were significantly associated with longer coleoptile, suggesting that QTL pyramiding is an effective approach to increase wheat coleoptile length.

     

Genome-wide association analysis reveals cryptic alleles as an important factor in heterosis for fatness in chicken F2 population

Genome-wide association studies have become possible in the chicken because of the recent availability of the complete genome sequence, a polymorphism map and high-density single nucleotide polymorphism (SNP) genotyping platforms. We used these tools to study the genetic basis of a very high level of heterosis that was previously observed for fatness in two F₂ populations established by crossing one outbred broiler (meat-type) sire with dams from two unrelated, highly inbred, light-bodied lines (Fayoumi and Leghorn). In each F₂ population, selective genotyping was carried out using phenotypically extreme males for abdominal fat percentage (AF) and about 3000 SNPs. Single-point association analysis of about 500 informative SNPs per cross showed significant association ( P  < 0.01) of 15 and 24 markers with AF in the Broiler × Fayoumi and Broiler × Leghorn crosses respectively. These SNPs were on 10 chromosomes (GGA1, 2, 3, 4, 7, 8, 10, 12, 15 and 27). Interestingly, of the 39 SNPs that were significantly associated with AF, there were about twice as many homozygous genotypes associated with higher AF that traced back to the inbred lines alleles, although the broiler line had on average higher AF. These SNPs are considered to be associated with QTL with cryptic alleles. This study reveals cryptic alleles as an important factor in heterosis for fatness observed in two chicken F₂ populations, and suggests epistasis as the common underlying mechanism for heterosis and cryptic allele expression. The results of this study also demonstrate the power of high marker-density SNP association studies in discovering QTL that were not detected by previous microsatellite-based genotyping studies.     

Genome-wide analysis of lysine catabolism in bacteria reveals new connections with osmotic stress resistance

Lysine is catabolized via the saccharopine pathway in plants and mammals. In this pathway, lysine is converted to α-aminoadipic-δ-semialdehyde (AASA) by lysine-ketoglutarate reductase/saccharopine dehydrogenase (LKR/SDH); thereafter, AASA is converted to aminoadipic acid (AAA) by α-aminoadipic-δ-semialdehyde dehydrogenase (AASADH). Here, we investigate the occurrence, genomic organization and functional role of lysine catabolic pathways among prokaryotes. Surprisingly, only 27 species of the 1478 analyzed contain the lkr and sdh genes, whereas 323 species contain aasadh orthologs. A sdh-related gene, identified in 159 organisms, was frequently found contiguously to an aasadh gene. This gene, annotated as lysine dehydrogenase (lysdh), encodes LYSDH an enzyme that directly converts lysine to AASA. Pipecolate oxidase (PIPOX) and lysine-6-aminotransferase (LAT), that converts lysine to AASA, were also found associated with aasadh. Interestingly, many lysdh–aasadh–containing organisms live under hyperosmotic stress. To test the role of the lysine-to-AASA pathways in the bacterial stress response, we subjected Silicibacter pomeroyi to salt stress. All but lkr, sdh, lysdh and aasadh were upregulated under salt stress conditions. In addition, lysine-supplemented culture medium increased the growth rate of S. pomeroyi under high-salt conditions and induced high-level expression of the lysdh–aasadh operon. Finally, transformation of Escherichia coli with the S. pomeroyi lysdh–aasadh operon resulted in increased salt tolerance. The transformed E. coli accumulated high levels of the compatible solute pipecolate, which may account for the salt resistance. These findings suggest that the lysine-to-AASA pathways identified in this work may have a broad evolutionary importance in osmotic stress resistance.     

Genome-wide association study reveals five nucleotide sequence variants for carcass traits in beef cattle

Genetic associations of nucleotide sequence variants with carcass traits in beef cattle were investigated using a genome-wide single nucleotide polymorphism (SNP) assay. Three hundred and thirteen Korean cattle were genotyped with the Illumina BovineSNP50 BeadChip, and 39,129 SNPs from 311 animals were analysed for each carcass phenotype after filtering by quality assurance. Five sequence markers were associated with one of the meat quantity or quality traits; rs109593638 on chromosome 3 with marbling score, rs109821175 on chromosome 11 and rs110862496 on chromosome 13 with backfat thickness (BFT), and rs110228023 on chromosome 6 and rs110201414 on chromosome 16 with eye muscle area (EMA) (P < 1.27 × 10(-6) , Bonferonni P < 0.05). The ss96319521 SNP, located within a gene with functions of muscle development, dishevelled homolog 1 (DVL1), would be a desirable candidate marker. Individuals with genotype CC at this gene appeared to have increased both EMA and carcass weight. Fine-mapping would be required to refine each of the five association signals shown in the current study for future application in marker-assisted selection for genetic improvement of beef quality and quantity.     

Genome-wide association study of fish oil supplementation on lipid traits in 81,246 individuals reveals new gene-diet interaction loci

Fish oil supplementation is widely used for reducing serum triglycerides (TAGs) but has mixed effects on other circulating cardiovascular biomarkers. Many genetic polymorphisms have been associated with blood lipids, including high- and low-density-lipoprotein cholesterol (HDL-C, LDL-C), total cholesterol, and TAGs. Here, the gene-diet interaction effects of fish oil supplementation on these lipids were analyzed in a discovery cohort of up to 73,962 UK Biobank participants, using a 1-degree-of-freedom (1df) test for interaction effects and a 2-degrees-of-freedom (2df) test to jointly analyze interaction and main effects. Associations with P < 1×10⁻⁶ in either test (26,157; 18,300 unique variants) were advanced to replication in up to 7,284 participants from the Atherosclerosis Risk in Communities (ARIC) Study. Replicated associations reaching 1df P < 0.05 (2,175; 1,763 unique variants) were used in meta-analyses. We found 13 replicated and 159 non-replicated (UK Biobank only) loci with significant 2df joint tests that were predominantly driven by main effects and have been previously reported. Four novel interaction loci were identified with 1df P < 5×10⁻⁸ in meta-analysis. The lead variant in the GJB6-GJB2-GJA3 gene cluster, rs112803755 (A>G; minor allele frequency = 0.041), shows exclusively interaction effects. The minor allele is significantly associated with decreased TAGs in individuals with fish oil supplementation, but with increased TAGs in those without supplementation. This locus is significantly associated with higher GJB2 expression of connexin 26 in adipose tissue; connexin activity is known to change upon exposure to omega-3 fatty acids. Significant interaction effects were also found in three other loci in the genes SLC12A3 (HDL-C), ABCA6 (LDL-C), and MLXIPL (LDL-C), but highly significant main effects are also present. Our study identifies novel gene-diet interaction effects for four genetic loci, whose effects on blood lipids are modified by fish oil supplementation. These findings highlight the need and possibility for personalized nutrition.