Karol Estrada

Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus

To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1×10−11 observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ±500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6×10−31 and an effect size explaining between 0.6%–1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42×10−10) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9×10−16) and rs7801723 (P = 8.9×10−28), also mapping to C7orf58 (r2 = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life.

Evidence of Inbreeding Depression on Human Height

Stature is a classical and highly heritable complex trait, with 80%–90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ2 = 83.89, df = 1; p = 5.2×10−20). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

GRIMP: A web- and grid-based tool for high-speed analysis of large-scale genome-wide association using imputed data

Summary: The current fast growth of genome-wide association studies (GWAS) combined with now common computationally expensive imputation requires the online access of large user groups to high-performance computing resources capable of analyzing rapidly and efficiently millions of genetic markers for ten thousands of individuals. Here, we present a web-based interface—called GRIMP—to run publicly available genetic software for extremely large GWAS on scalable super-computing grid infrastructures. This is of major importance for the enlargement of GWAS with the availability of whole-genome sequence data from the 1000 Genomes Project and for future whole-population efforts. Contact: ta.knoch@taknoch.org; f.rivadeneira@erasmusmc.nl

Mapping of Gene Expression Reveals CYP27A1 as a Susceptibility Gene for Sporadic ALS

Amyotrophic lateral sclerosis (ALS) is a progressive, neurodegenerative disease characterized by loss of upper and lower motor neurons. ALS is considered to be a complex trait and genome-wide association studies (GWAS) have implicated a few susceptibility loci. However, many more causal loci remain to be discovered. Since it has been shown that genetic variants associated with complex traits are more likely to be eQTLs than frequency-matched variants from GWAS platforms, we conducted a two-stage genome-wide screening for eQTLs associated with ALS. In addition, we applied an eQTL analysis to finemap association loci. Expression profiles using peripheral blood of 323 sporadic ALS patients and 413 controls were mapped to genome-wide genotyping data. Subsequently, data from a two-stage GWAS (3,568 patients and 10,163 controls) were used to prioritize eQTLs identified in the first stage (162 ALS, 207 controls). These prioritized eQTLs were carried forward to the second sample with both gene-expression and genotyping data (161 ALS, 206 controls). Replicated eQTL SNPs were then tested for association in the second-stage GWAS data to find SNPs associated with disease, that survived correction for multiple testing. We thus identified twelve cis eQTLs with nominally significant associations in the second-stage GWAS data. Eight SNP-transcript pairs of highest significance (lowest p = 1.27×10−51) withstood multiple-testing correction in the second stage and modulated CYP27A1 gene expression. Additionally, we show that C9orf72 appears to be the only gene in the 9p21.2 locus that is regulated in cis, showing the potential of this approach in identifying causative genes in association loci in ALS. This study has identified candidate genes for sporadic ALS, most notably CYP27A1. Mutations in CYP27A1 are causal to cerebrotendinous xanthomatosis which can present as a clinical mimic of ALS with progressive upper motor neuron loss, making it a plausible susceptibility gene for ALS.

Assessment of Gene-by-Sex Interaction Effect on Bone Mineral Density

Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene‐by‐sex autosomal interactions genome‐wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome‐wide meta‐analysis of gene‐by‐sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single‐nucleotide polymorphisms (SNPs; p < 1 × 10−5) in an additional set of 24,763 individuals. Gene‐by‐sex interaction and sex‐specific effects were examined in these 12 SNPs. We detected one novel genome‐wide significant interaction associated with LS‐BMD at the Chr3p26.1‐p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10−5; female effect = −0.007 and p = 3.3 × 10−2), and 11 suggestive loci associated with either FN‐ or LS‐BMD in discovery cohorts. However, there was no evidence for genome‐wide significant (p < 5 × 10−8) gene‐by‐sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome‐wide significant evidence for gene‐by‐sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene‐by‐sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP.

C9orf72 and UNC13A are shared risk loci for amyotrophic lateral sclerosis and frontotemporal dementia: A genome-wide meta-analysis

Substantial clinical, pathological, and genetic overlap exists between amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP‐43 inclusions have been found in both ALS and FTD cases (FTD‐TDP). Recently, a repeat expansion in C9orf72 was identified as the causal variant in a proportion of ALS and FTD cases. We sought to identify additional evidence for a common genetic basis for the spectrum of ALS‐FTD.

Phenome-wide association studies across large population cohorts support drug target validation

Phenome-wide association studies (PheWAS) have been proposed as a possible aid in drug development through elucidating mechanisms of action, identifying alternative indications, or predicting adverse drug events (ADEs). Here, we select 25 single nucleotide polymorphisms (SNPs) linked through genome-wide association studies (GWAS) to 19 candidate drug targets for common disease indications. We interrogate these SNPs by PheWAS in four large cohorts with extensive health information (23andMe, UK Biobank, FINRISK, CHOP) for association with 1683 binary endpoints in up to 697,815 individuals and conduct meta-analyses for 145 mapped disease endpoints. Our analyses replicate 75% of known GWAS associations (P < 0.05) and identify nine study-wide significant novel associations (of 71 with FDR < 0.1). We describe associations that may predict ADEs, e.g., acne, high cholesterol, gout, and gallstones with rs738409 (p.I148M) in PNPLA3 and asthma with rs1990760 (p.T946A) in IFIH1. Our results demonstrate PheWAS as a powerful addition to the toolkit for drug discovery.Testing the association between genetic variants and a range of phenotypes can assist drug development. Here, in a phenome-wide association study in up to 697,815 individuals, Diogo et al. identify genotype–phenotype associations predicting efficacy, alternative indications or adverse drug effects.

A whole-genome sequence study identifies genetic risk factors for neuromyelitis optica

Neuromyelitis optica (NMO) is a rare autoimmune disease that affects the optic nerve and spinal cord. Most NMO patients ( > 70%) are seropositive for circulating autoantibodies against aquaporin 4 (NMO-IgG+). Here, we meta-analyze whole-genome sequences from 86 NMO cases and 460 controls with genome-wide SNP array from 129 NMO cases and 784 controls to test for association with SNPs and copy number variation (total N = 215 NMO cases, 1244 controls). We identify two independent signals in the major histocompatibility complex (MHC) region associated with NMO-IgG+, one of which may be explained by structural variation in the complement component 4 genes. Mendelian Randomization analysis reveals a significant causal effect of known systemic lupus erythematosus (SLE), but not multiple sclerosis (MS), risk variants in NMO-IgG+. Our results suggest that genetic variants in the MHC region contribute to the etiology of NMO-IgG+ and that NMO-IgG+ is genetically more similar to SLE than MS.Neuromyelitis optica (NMO) is a rare autoimmune condition characterized by inflammation and demyelination of the optic nerve and the spinal cord. Here, Estrada et al. identify NMO susceptibility variants in the MHC region and find that autoantibody-positive NMO genetically overlaps with lupus.