Grant W. Montgomery

Genome-wide meta-analysis increases to 71 the number of confirmed Crohn's disease susceptibility loci

We undertook a meta-analysis of six Crohns disease genome-wide association studies (GWAS) comprising 6,333 affected individuals (cases) and 15,056 controls and followed up the top association signals in 15,694 cases, 14,026 controls and 414 parent-offspring trios. We identified 30 new susceptibility loci meeting genome-wide significance (P < 5 × 10−8). A series of in silico analyses highlighted particular genes within these loci and, together with manual curation, implicated functionally interesting candidate genes including SMAD3, ERAP2, IL10, IL2RA, TYK2, FUT2, DNMT3A, DENND1B, BACH2 and TAGAP. Combined with previously confirmed loci, these results identify 71 distinct loci with genome-wide significant evidence for association with Crohns disease.

Mutations in an oocyte-derived growth factor gene ( BMP15 ) cause increased ovulation rate and infertility in a dosage-sensitive manner

Multiple ovulations are uncommon in humans, cattle and many breeds of sheep. Pituitary gonadotrophins and as yet unidentified ovarian factors precisely regulate follicular development so that, normally, only one follicle is selected to ovulate. The Inverdale (FecXI) sheep, however, carries a naturally occurring X-linked mutation that causes increased ovulation rate and twin and triplet births in heterozygotes (FecXI/FecX+; ref. 1), but primary ovarian failure in homozygotes (FecXI/FecXI; ref. 2). Germ-cell development, formation of the follicle and the earliest stages of follicular growth are normal in FecXI/FecXI sheep, but follicular development beyond the primary stage is impaired. A second family unrelated to the Inverdale sheep also has the same X-linked phenotype (Hanna, FecXH). Crossing FecXI with FecXH animals produces FecXI/FecXH infertile females phenotypically indistinguishable from FecXI/FecXI females. We report here that the FecXI locus maps to an orthologous chromosomal region syntenic to human Xp11.2–11.4, which contains BMP15, encoding bone morphogenetic protein 15 (also known as growth differentiation factor 9B (GDF9B)). Whereas BMP15 is a member of the transforming growth factor β (TGFβ) superfamily and is specifically expressed in oocytes, its function is unknown. We show that independent germline point mutations exist in FecXI and FecXH carriers. These findings establish that BMP15 is essential for female fertility and that natural mutations in an ovary-derived factor can cause both increased ovulation rate and infertility phenotypes in a dosage-sensitive manner.

Systematic identification of trans eQTLs as putative drivers of known disease associations

Identifying the downstream effects of disease-associated SNPs is challenging. To help overcome this problem, we performed expression quantitative trait locus (eQTL) meta-analysis in non-transformed peripheral blood samples from 5,311 individuals with replication in 2,775 individuals. We identified and replicated trans eQTLs for 233 SNPs (reflecting 103 independent loci) that were previously associated with complex traits at genome-wide significance. Some of these SNPs affect multiple genes in trans that are known to be altered in individuals with disease: rs4917014, previously associated with systemic lupus erythematosus (SLE), altered gene expression of C1QB and five type I interferon response genes, both hallmarks of SLE. DeepSAGE RNA sequencing showed that rs4917014 strongly alters the 3′ UTR levels of IKZF1 in cis, and chromatin immunoprecipitation and sequencing analysis of the trans-regulated genes implicated IKZF1 as the causal gene. Variants associated with cholesterol metabolism and type 1 diabetes showed similar phenomena, indicating that large-scale eQTL mapping provides insight into the downstream effects of many trait-associated variants.

Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis (vol 5, 4926, 2014)

Corrigendum: Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis

Loci influencing lipid levels and coronary heart disease risk in 16 European population cohorts

Recent genome-wide association (GWA) studies of lipids have been conducted in samples ascertained for other phenotypes, particularly diabetes. Here we report the first GWA analysis of loci affecting total cholesterol (TC), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol and triglycerides sampled randomly from 16 population-based cohorts and genotyped using mainly the Illumina HumanHap300-Duo platform. Our study included a total of 17,797–22,562 persons, aged 18–104 years and from geographic regions spanning from the Nordic countries to Southern Europe. We established 22 loci associated with serum lipid levels at a genome-wide significance level (P < 5 × 10−8), including 16 loci that were identified by previous GWA studies. The six newly identified loci in our cohort samples are ABCG5 (TC, P = 1.5 × 10−11; LDL, P = 2.6 × 10−10), TMEM57 (TC, P = 5.4 × 10−10), CTCF-PRMT8 region (HDL, P = 8.3 × 10−16), DNAH11 (LDL, P = 6.1 × 10−9), FADS3-FADS2 (TC, P = 1.5 × 10−10; LDL, P = 4.4 × 10−13) and MADD-FOLH1 region (HDL, P = 6 × 10−11). For three loci, effect sizes differed significantly by sex. Genetic risk scores based on lipid loci explain up to 4.8% of variation in lipids and were also associated with increased intima media thickness (P = 0.001) and coronary heart disease incidence (P = 0.04). The genetic risk score improves the screening of high-risk groups of dyslipidemia over classical risk factors.

A Versatile Gene-Based Test for Genome-wide Association Studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

DNA methylation profiles in monozygotic and dizygotic twins

Twin studies have provided the basis for genetic and epidemiological studies in human complex traits. As epigenetic factors can contribute to phenotypic outcomes, we conducted a DNA methylation analysis in white blood cells (WBC), buccal epithelial cells and gut biopsies of 114 monozygotic (MZ) twins as well as WBC and buccal epithelial cells of 80 dizygotic (DZ) twins using 12K CpG island microarrays. Here we provide the first annotation of epigenetic metastability of ∼6,000 unique genomic regions in MZ twins. An intraclass correlation (ICC)-based comparison of matched MZ and DZ twins showed significantly higher epigenetic difference in buccal cells of DZ co-twins (P = 1.2 × 10−294). Although such higher epigenetic discordance in DZ twins can result from DNA sequence differences, our in silico SNP analyses and animal studies favor the hypothesis that it is due to epigenomic differences in the zygotes, suggesting that molecular mechanisms of heritability may not be limited to DNA sequence differences.

Genetic loci for retinal arteriolar microcirculation

Narrow arterioles in the retina have been shown to predict hypertension as well as other vascular diseases, likely through an increase in the peripheral resistance of the microcirculatory flow. In this study, we performed a genome-wide association study in 18,722 unrelated individuals of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium and the Blue Mountain Eye Study, to identify genetic determinants associated with variations in retinal arteriolar caliber. Retinal vascular calibers were measured on digitized retinal photographs using a standardized protocol. One variant (rs2194025 on chromosome 5q14 near the myocyte enhancer factor 2C MEF2C gene) was associated with retinal arteriolar caliber in the meta-analysis of the discovery cohorts at genome-wide significance of P-value <5×10−8. This variant was replicated in an additional 3,939 individuals of European ancestry from the Australian Twins Study and Multi-Ethnic Study of Atherosclerosis (rs2194025, P-value = 2.11×10−12 in combined meta-analysis of discovery and replication cohorts). In independent studies of modest sample sizes, no significant association was found between this variant and clinical outcomes including coronary artery disease, stroke, myocardial infarction or hypertension. In conclusion, we found one novel loci which underlie genetic variation in microvasculature which may be relevant to vascular disease. The relevance of these findings to clinical outcomes remains to be determined.

Highly Prolific Booroola Sheep Have a Mutation in the Intracellular Kinase Domain of Bone Morphogenetic Protein IB Receptor (ALK-6) That Is Expressed in Both Oocytes and Granulosa Cells

Abstract The Booroola fecundity gene (FecB) increases ovulation rate and litter size in sheep and is inherited as a single autosomal locus. The effect of FecB is additive for ovulation rate (increasing by about 1.6 corpora lutea per cycle for each copy) and has been mapped to sheep chromosome 6q23–31, which is syntenic to human chromosome 4q21–25. Bone morphogenetic protein IB (BMP-IB) receptor (also known as ALK-6), which binds members of the transforming growth factor-β (TGF-β) superfamily, is located in the region containing the FecB locus. Booroola sheep have a mutation (Q249R) in the highly conserved intracellular kinase signaling domain of the BMP-IB receptor. The mutation segregated with the FecB phenotype in the Booroola backcross and half-sib flocks of sheep with no recombinants. The mutation was not found in individuals from a number of sheep breeds not derived from the Booroola strain. BMPR-IB was expressed in the ovary and in situ hybridization revealed its specific location to the oocyte and the granulosa cell. Expression of mRNA encoding the BMP type II receptor was widespread throughout the ovary. The mutation in BMPR-IB found in Booroola sheep is the second reported defect in a gene from the TGF-β pathway affecting fertility in sheep following the recent discovery of mutations in the growth factor, GDF9b/BMP15.

Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits

We present an approximate conditional and joint association analysis that can use summary-level statistics from a meta-analysis of genome-wide association studies (GWAS) and estimated linkage disequilibrium (LD) from a reference sample with individual-level genotype data. Using this method, we analyzed meta-analysis summary data from the GIANT Consortium for height and body mass index (BMI), with the LD structure estimated from genotype data in two independent cohorts. We identified 36 loci with multiple associated variants for height (38 leading and 49 additional SNPs, 87 in total) via a genome-wide SNP selection procedure. The 49 new SNPs explain approximately 1.3% of variance, nearly doubling the heritability explained at the 36 loci. We did not find any locus showing multiple associated SNPs for BMI. The method we present is computationally fast and is also applicable to case-control data, which we demonstrate in an example from meta-analysis of type 2 diabetes by the DIAGRAM Consortium.