Maternal Filaggrin Mutations Increase the Risk of Atopic Dermatitis in Children: An Effect Independent of Mutation Inheritance

Epidemiological studies suggest that allergy risk is preferentially transmitted through mothers. This can be due to genomic imprinting, where the phenotype effect of an allele depends on its parental origin, or due to maternal effects reflecting the maternal genome''s influence on the child during prenatal development. Loss-of-function mutations in the filaggrin gene (FLG) cause skin barrier deficiency and strongly predispose to atopic dermatitis (AD). We investigated the 4 most prevalent European FLG mutations (c.2282del4, p.R501X, p.R2447X, and p.S3247X) in two samples including 759 and 450 AD families. We used the multinomial and maximum-likelihood approach implemented in the PREMIM/EMIM tool to model parent-of-origin effects. Beyond the known role of FLG inheritance in AD (R1_{meta-analysis} = 2.4, P = 1.0 x 10⁻³⁶), we observed a strong maternal FLG genotype effect that was consistent in both independent family sets and for all 4 mutations analysed. Overall, children of FLG-carrier mothers had a 1.5-fold increased AD risk (S1 = 1.50, P_{meta-analysis} = 8.4 x 10⁻⁸). Our data point to two independent and additive effects of FLG mutations: i) carrying a mutation and ii) having a mutation carrier mother. The maternal genotype effect was independent of mutation inheritance and can be seen as a non-genetic transmission of a genetic effect. The FLG maternal effect was observed only when mothers had allergic sensitization (elevated allergen-specific IgE antibody plasma levels), suggesting that FLG mutation-induced systemic immune responses in the mother may influence AD risk in the child. Notably, the maternal effect reported here was stronger than most common genetic risk factors for AD recently identified through genome-wide association studies (GWAS). Our study highlights the power of family-based studies in the identification of new etiological mechanisms and reveals, for the first time, a direct influence of the maternal genotype on the offspring’s susceptibility to a common human disease.     

Assessing the Causal Relationship of Maternal Height on Birth Size and Gestational Age at Birth: A Mendelian Randomization Analysis

Background

Observational epidemiological studies indicate that maternal height is associated with gestational age at birth and fetal growth measures (i.e., shorter mothers deliver infants at earlier gestational ages with lower birth weight and birth length). Different mechanisms have been postulated to explain these associations. This study aimed to investigate the casual relationships behind the strong association of maternal height with fetal growth measures (i.e., birth length and birth weight) and gestational age by a Mendelian randomization approach.

Methods and Findings

We conducted a Mendelian randomization analysis using phenotype and genome-wide single nucleotide polymorphism (SNP) data of 3,485 mother/infant pairs from birth cohorts collected from three Nordic countries (Finland, Denmark, and Norway). We constructed a genetic score based on 697 SNPs known to be associated with adult height to index maternal height. To avoid confounding due to genetic sharing between mother and infant, we inferred parental transmission of the height-associated SNPs and utilized the haplotype genetic score derived from nontransmitted alleles as a valid genetic instrument for maternal height. In observational analysis, maternal height was significantly associated with birth length (p = 6.31 × 10⁻⁹), birth weight (p = 2.19 × 10⁻¹⁵), and gestational age (p = 1.51 × 10⁻⁷). Our parental-specific haplotype score association analysis revealed that birth length and birth weight were significantly associated with the maternal transmitted haplotype score as well as the paternal transmitted haplotype score. Their association with the maternal nontransmitted haplotype score was far less significant, indicating a major fetal genetic influence on these fetal growth measures. In contrast, gestational age was significantly associated with the nontransmitted haplotype score (p = 0.0424) and demonstrated a significant (p = 0.0234) causal effect of every 1 cm increase in maternal height resulting in ~0.4 more gestational d. Limitations of this study include potential influences in causal inference by biological pleiotropy, assortative mating, and the nonrandom sampling of study subjects.

Conclusions

Our results demonstrate that the observed association between maternal height and fetal growth measures (i.e., birth length and birth weight) is mainly defined by fetal genetics. In contrast, the association between maternal height and gestational age is more likely to be causal. In addition, our approach that utilizes the genetic score derived from the nontransmitted maternal haplotype as a genetic instrument is a novel extension to the Mendelian randomization methodology in casual inference between parental phenotype (or exposure) and outcomes in offspring.     

How Genome-Wide SNP-SNP Interactions Relate to Nasopharyngeal Carcinoma Susceptibility

This study is the first to use genome-wide association study (GWAS) data to evaluate the multidimensional genetic architecture underlying nasopharyngeal cancer. Since analysis of data from GWAS confirms a close and consistent association between elevated risk for nasopharyngeal carcinoma (NPC) and major histocompatibility complex class 1 genes, our goal here was to explore lesser effects of gene-gene interactions. We conducted an exhaustive genome-wide analysis of GWAS data of NPC, revealing two-locus interactions occurring between single nucleotide polymorphisms (SNPs), and identified a number of suggestive interaction loci which were missed by traditional GWAS analyses. Although none of the interaction pairs we identified passed the genome-wide Bonferroni-adjusted threshold for significance, using independent GWAS data from the same population (Stage 2), we selected 66 SNP pairs in 39 clusters with P<0.01. We identified that in several chromosome regions, multiple suggestive interactions group to form a block-like signal, effectively reducing the rate of false discovery. The strongest cluster of interactions involved the CREB5 gene and a SNP rs1607979 on chromosome 17q22 (P = 9.86×10⁻¹¹) which also show trans-expression quantitative loci (eQTL) association in Chinese population. We then detected a complicated cis-interaction pattern around the NPC-associated HLA-B locus, which is immediately adjacent to copy-number variations implicated in male susceptibility for NPC. While it remains to be seen exactly how and to what degree SNP-SNP interactions such as these affect susceptibility for nasopharyngeal cancer, future research on these questions holds great promise for increasing our understanding of this disease’s genetic etiology, and possibly also that of other gene-related cancers.     

Genome-Wide Association Study of Lp-PLA2 Activity and Mass in the Framingham Heart Study

Lipoprotein-associated phospholipase A₂ (Lp-PLA₂) is an emerging risk factor and therapeutic target for cardiovascular disease. The activity and mass of this enzyme are heritable traits, but major genetic determinants have not been explored in a systematic, genome-wide fashion. We carried out a genome-wide association study of Lp-PLA₂ activity and mass in 6,668 Caucasian subjects from the population-based Framingham Heart Study. Clinical data and genotypes from the Affymetrix 550K SNP array were obtained from the open-access Framingham SHARe project. Each polymorphism that passed quality control was tested for associations with Lp-PLA₂ activity and mass using linear mixed models implemented in the R statistical package, accounting for familial correlations, and controlling for age, sex, smoking, lipid-lowering-medication use, and cohort. For Lp-PLA₂ activity, polymorphisms at four independent loci reached genome-wide significance, including the APOE/APOC1 region on chromosome 19 (p = 6×10⁻²⁴); CELSR2/PSRC1 on chromosome 1 (p = 3×10⁻¹⁵); SCARB1 on chromosome 12 (p = 1×10⁻⁸) and ZNF259/BUD13 in the APOA5/APOA1 gene region on chromosome 11 (p = 4×10⁻⁸). All of these remained significant after accounting for associations with LDL cholesterol, HDL cholesterol, or triglycerides. For Lp-PLA₂ mass, 12 SNPs achieved genome-wide significance, all clustering in a region on chromosome 6p12.3 near the PLA2G7 gene. Our analyses demonstrate that genetic polymorphisms may contribute to inter-individual variation in Lp-PLA₂ activity and mass.     

Genome-Wide Joint Meta-Analysis of SNP and SNP-by-Smoking Interaction Identifies Novel Loci for Pulmonary Function

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV₁), and its ratio to forced vital capacity (FEV₁/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV₁ and FEV₁/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P _JMA =5.00×10⁻¹¹), HLA-DQB1 and HLA-DQA2 (smallest P _JMA =4.35×10⁻⁹), and KCNJ2 and SOX9 (smallest P _JMA =1.28×10⁻⁸) were associated with FEV₁/FVC or FEV₁ in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.     

Genome-Wide Association Studies of Asthma in Population-Based Cohorts Confirm Known and Suggested Loci and Identify an Additional Association near HLA

Rationale

Asthma has substantial morbidity and mortality and a strong genetic component, but identification of genetic risk factors is limited by availability of suitable studies.

Objectives

To test if population-based cohorts with self-reported physician-diagnosed asthma and genome-wide association (GWA) data could be used to validate known associations with asthma and identify novel associations.

Methods

The APCAT (Analysis in Population-based Cohorts of Asthma Traits) consortium consists of 1,716 individuals with asthma and 16,888 healthy controls from six European-descent population-based cohorts. We examined associations in APCAT of thirteen variants previously reported as genome-wide significant (P<5x10⁻⁸) and three variants reported as suggestive (P<5×10⁻⁷). We also searched for novel associations in APCAT (Stage 1) and followed-up the most promising variants in 4,035 asthmatics and 11,251 healthy controls (Stage 2). Finally, we conducted the first genome-wide screen for interactions with smoking or hay fever.

Main Results

We observed association in the same direction for all thirteen previously reported variants and nominally replicated ten of them. One variant that was previously suggestive, rs11071559 in RORA, now reaches genome-wide significance when combined with our data (P = 2.4×10⁻⁹). We also identified two genome-wide significant associations: rs13408661 near IL1RL1/IL18R1 (P _{Stage1+Stage2} = 1.1x10⁻⁹), which is correlated with a variant recently shown to be associated with asthma (rs3771180), and rs9268516 in the HLA region (P _{Stage1+Stage2} = 1.1x10⁻⁸), which appears to be independent of previously reported associations in this locus. Finally, we found no strong evidence for gene-environment interactions with smoking or hay fever status.

Conclusions

Population-based cohorts with simple asthma phenotypes represent a valuable and largely untapped resource for genetic studies of asthma.     

Genome-wide association identifies a common variant in the reelin gene that increases the risk of schizophrenia only in women

Sex differences in schizophrenia are well known, but their genetic basis has not been identified. We performed a genome-wide association scan for schizophrenia in an Ashkenazi Jewish population using DNA pooling. We found a female-specific association with rs7341475, a SNP in the fourth intron of the reelin (RELN) gene (p = 2.9 × 10⁻⁵ in women), with a significant gene-sex effect (p = 1.8 × 10⁻⁴). We studied rs7341475 in four additional populations, totaling 2,274 cases and 4,401 controls. A significant effect was observed only in women, replicating the initial result (p = 2.1 × 10⁻³ in women; p = 4.2 × 10⁻³ for gene-sex interaction). Based on all populations the estimated relative risk of women carrying the common genotype is 1.58 (p = 8.8 × 10⁻⁷; p = 1.6 × 10⁻⁵ for gene-sex interaction). The female-specific association between RELN and schizophrenia is one of the few examples of a replicated sex-specific genetic association in any disease.     

Pleiotropy method reveals genetic overlap between orofacial clefts at multiple novel loci from GWAS of multi-ethnic trios

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts– the most common craniofacial birth defects in humans– are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10⁻⁸, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10⁻⁶, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes.     

Resequencing and Association Analysis of PTPRA,a Possible Susceptibility Gene for Schizophrenia and Autism Spectrum Disorders

Background

The PTPRA gene, which encodes the protein RPTP-α, is critical to neurodevelopment. Previous linkage studies, genome-wide association studies, controlled expression analyses and animal models support an association with both schizophrenia and autism spectrum disorders, both of which share a substantial portion of genetic risks.

Methods

We sequenced the protein-encoding areas of the PTPRA gene for single nucleotide polymorphisms or small insertions/deletions (InDel) in 382 schizophrenia patients. To validate their association with the disorders, rare (minor allele frequency <1%), missense mutations as well as one InDel in the 3′UTR region were then genotyped in another independent sample set comprising 944 schizophrenia patients, 336 autism spectrum disorders patients, and 912 healthy controls.

Results

Eight rare mutations, including 3 novel variants, were identified during the mutation-screening phase. In the following association analysis, L59P, one of the two missense mutations, was only observed among patients of schizophrenia. Additionally, a novel duplication in the 3′UTR region, 174620_174623dupTGAT, was predicted to be located within a Musashi Binding Element.

Major Conclusions

No evidence was seen for the association of rare, missense mutations in the PTPRA gene with schizophrenia or autism spectrum disorders; however, we did find some rare variants with possibly damaging effects that may increase the susceptibility of carriers to the disorders.     

Stratification by Smoking Status Reveals an Association of CHRNA5-A3-B4 Genotype with Body Mass Index in Never Smokers

We previously used a single nucleotide polymorphism (SNP) in the CHRNA5-A3-B4 gene cluster associated with heaviness of smoking within smokers to confirm the causal effect of smoking in reducing body mass index (BMI) in a Mendelian randomisation analysis. While seeking to extend these findings in a larger sample we found that this SNP is associated with 0.74% lower body mass index (BMI) per minor allele in current smokers (95% CI -0.97 to -0.51, P = 2.00×10⁻¹⁰), but also unexpectedly found that it was associated with 0.35% higher BMI in never smokers (95% CI +0.18 to +0.52, P = 6.38×10⁻⁵). An interaction test confirmed that these estimates differed from each other (P = 4.95×10⁻¹³). This difference in effects suggests the variant influences BMI both via pathways unrelated to smoking, and via the weight-reducing effects of smoking. It would therefore be essentially undetectable in an unstratified genome-wide association study of BMI, given the opposite association with BMI in never and current smokers. This demonstrates that novel associations may be obscured by hidden population sub-structure. Stratification on well-characterized environmental factors known to impact on health outcomes may therefore reveal novel genetic associations.     

Results of Genome-Wide Analyses on Neurodevelopmental Phenotypes at Four-Year Follow-Up following Cardiac Surgery in Infancy

Background

Adverse neurodevelopmental sequelae are reported among children who undergo early cardiac surgery to repair congenital heart defects (CHD). APOE genotype has previously been determined to contribute to the prediction of these outcomes. Understanding further genetic causes for the development of poor neurobehavioral outcomes should enhance patient risk stratification and improve both prevention and treatment strategies.

Methods

We performed a prospective observational study of children who underwent cardiac surgery before six months of age; this included a neurodevelopmental evaluation between their fourth and fifth birthdays. Attention and behavioral skills were assessed through parental report utilizing the Attention Deficit-Hyperactivity Disorder-IV scale preschool edition (ADHD-IV), and Child Behavior Checklist (CBCL/1.5-5), respectively. Of the seven investigated, three neurodevelopmental phenotypes met genomic quality control criteria. Linear regression was performed to determine the effect of genome-wide genetic variation on these three neurodevelopmental measures in 316 subjects.

Results

This genome-wide association study identified single nucleotide polymorphisms (SNPs) associated with three neurobehavioral phenotypes in the postoperative children ADHD-IV Impulsivity/Hyperactivity, CBCL/1.5-5 PDPs, and CBCL/1.5-5 Total Problems. The most predictive SNPs for each phenotype were: a LGALS8 intronic SNP, rs4659682, associated with ADHD-IV Impulsivity (P = 1.03×10⁻⁶); a PCSK5 intronic SNP, rs2261722, associated with CBCL/1.5-5 PDPs (P = 1.11×10⁻⁶); and an intergenic SNP, rs11617488, 50 kb from FGF9, associated with CBCL/1.5-5 Total Problems (P = 3.47×10⁻⁷). 10 SNPs (3 for ADHD-IV Impulsivity, 5 for CBCL/1.5-5 PDPs, and 2 for CBCL/1.5-5 Total Problems) had p<10⁻⁵.

Conclusions

No SNPs met genome-wide significance for our three neurobehavioral phenotypes; however, 10 SNPs reached a threshold for suggestive significance (p<10⁻⁵). Given the unique nature of this cohort, larger studies and/or replication are not possible. Studies to further investigate the mechanisms through which these newly identified genes may influence neurodevelopment dysfunction are warranted.     

Longitudinal Genome-Wide Association of Cardiovascular Disease Risk Factors in the Bogalusa Heart Study

Cardiovascular disease (CVD) is the leading cause of death worldwide. Recent genome-wide association (GWA) studies have pinpointed many loci associated with CVD risk factors in adults. It is unclear, however, if these loci predict trait levels at all ages, if they are associated with how a trait develops over time, or if they could be used to screen individuals who are pre-symptomatic to provide the opportunity for preventive measures before disease onset. We completed a genome-wide association study on participants in the longitudinal Bogalusa Heart Study (BHS) and have characterized the association between genetic factors and the development of CVD risk factors from childhood to adulthood. We report 7 genome-wide significant associations involving CVD risk factors, two of which have been previously reported. Top regions were tested for replication in the Young Finns Study (YF) and two associations strongly replicated: rs247616 in CETP with HDL levels (combined P = 9.7×10⁻²⁴), and rs445925 at APOE with LDL levels (combined P = 8.7×10⁻¹⁹). We show that SNPs previously identified in adult cross-sectional studies tend to show age-independent effects in the BHS with effect sizes consistent with previous reports. Previously identified variants were associated with adult trait levels above and beyond those seen in childhood; however, variants with time-dependent effects were also promising predictors. This is the first GWA study to evaluate the role of common genetic variants in the development of CVD risk factors in children as they advance through adulthood and highlights the utility of using longitudinal studies to identify genetic predictors of adult traits in children.     

A Genome-Wide Association Study Uncovers a Genetic Locus Associated with Thoracic-to-Hip Ratio in Koreans

The thoracic-to-hip circumference ratio (THR) is an anthropometric marker recently described as a predictor of type 2 diabetes. In this study, we performed a genome-wide association study (GWAS) followed by confirmatory analyses to identify genetic markers associated with THR. A total of 7,240 Korean subjects (4,988 for the discovery stage and 2,252 for the confirmatory analyses) were recruited for this study, and genome-wide single nucleotide polymorphism (SNP) genotyping of the initial 4,988 individuals was performed using Affymetrix Human SNP array 5.0. Linear regression analysis was then performed to adjust for the effects of age, sex, and current diabetes medication status on the THR of the study subjects. In the initial discovery stage, there was a statistically nominal association between minor alleles of SNP markers on chromosomes 4, 8, 10, and 12, and THR changes (p < 5.0 × 10⁻⁶). The subsequent confirmatory analyses of these markers, however, only detected a significant association between two SNPs in the HECTD4 gene and decreased THRs. Notably, this association was detected in male (rs11066280: p = 1.14 × 10⁻²; rs2074356: p = 1.10 × 10⁻²), but not in female subjects. Meanwhile, the combined results from the two analyses (initial and confirmatory) indicated that minor alleles of these two intronic variants exhibited a significant genome-wide association with decreased THR in the male subjects (n = 3,155; rs11066280: effect size = −0.008624, p = 6.19 × 10⁻⁹; rs2074356: effect size = −0.008762, p = 1.89 × 10⁻⁸). Furthermore, minor alleles of these two SNPs exhibited protective effects on patients’ risks for developing type 2 diabetes. In conclusion, we have identified two genetic variations in HECTD4 that are associated with THR, particularly in men.     

Nonadditive Effects of Genes in Human Metabolomics

Genome-wide association studies (GWAS) are widely applied to analyze the genetic effects on phenotypes. With the availability of high-throughput technologies for metabolite measurements, GWAS successfully identified loci that affect metabolite concentrations and underlying pathways. In most GWAS, the effect of each SNP on the phenotype is assumed to be additive. Other genetic models such as recessive, dominant, or overdominant were considered only by very few studies. In contrast to this, there are theories that emphasize the relevance of nonadditive effects as a consequence of physiologic mechanisms. This might be especially important for metabolites because these intermediate phenotypes are closer to the underlying pathways than other traits or diseases. In this study we analyzed systematically nonadditive effects on a large panel of serum metabolites and all possible ratios (22,801 total) in a population-based study [Cooperative Health Research in the Region of Augsburg (KORA) F4, N = 1,785]. We applied four different 1-degree-of-freedom (1-df) tests corresponding to an additive, dominant, recessive, and overdominant trait model as well as a genotypic model with two degree-of-freedom (2-df) that allows a more general consideration of genetic effects. Twenty-three loci were found to be genome-wide significantly associated (Bonferroni corrected P ≤ 2.19 × 10⁻¹²) with at least one metabolite or ratio. For five of them, we show the evidence of nonadditive effects. We replicated 17 loci, including 3 loci with nonadditive effects, in an independent study (TwinsUK, N = 846). In conclusion, we found that most genetic effects on metabolite concentrations and ratios were indeed additive, which verifies the practice of using the additive model for analyzing SNP effects on metabolites.     

Haplotypes at LBX1 Have Distinct Inheritance Patterns with Opposite Effects in Adolescent Idiopathic Scoliosis

Adolescent idiopathic scoliosis (AIS) is a clinically significant disorder with high heritability that affects 2–4% of the population. Genome-wide association studies have identified LBX1 as a strong susceptibility locus for AIS in Asian and Caucasian populations. Here we further dissect the genetic association with AIS in a Caucasian population. To identify genetic markers associated with AIS we employed a genome-wide association study (GWAS) design comparing 620 female Caucasian patients who developed idiopathic scoliosis during adolescence with 1,287 ethnically matched females who had normal spinal curves by skeletal maturity. The genomic region around LBX1 was imputed and haplotypes investigated for genetic signals under different inheritance models. The strongest signal was identified upstream of LBX1 (rs11190878, P_trend = 4.18×10^-9, OR = 0.63[0.54–0.74]). None of the remaining SNPs pass the genome-wide significance threshold. We found rs11190870, downstream of LBX1 and previously associated with AIS in Asian populations, to be in modest linkage disequilibrium (LD) with rs11190878 (r² = 0.40, D'' = 0.81). Haplotype analysis shows that rs11190870 and rs11190878 track a single risk factor that resides on the ancestral haplotype and is shared across ethnic groups. We identify six haplotypes at the LBX1 locus including two strongly associated haplotypes; a recessive risk haplotype (TTA, Control_freq = 0.52, P = 1.25×10^-9, OR = 1.56), and a co-dominant protective haplotype (CCG, Control_freq = 0.28, P = 2.75×10^-7, OR = 0.65). Together the association signals from LBX1 explain 1.4% of phenotypic variance. Our results identify two clinically relevant haplotypes in the LBX1-region with opposite effects on AIS risk. The study demonstrates the utility of haplotypes over un-phased SNPs for individualized risk assessment by more strongly delineating individuals at risk for AIS without compromising the effect size.     

Parent-Of-Origin Effects in Autism Identified through Genome-Wide Linkage Analysis of 16,000 SNPs

Background

Autism is a common heritable neurodevelopmental disorder with complex etiology. Several genome-wide linkage and association scans have been carried out to identify regions harboring genes related to autism or autism spectrum disorders, with mixed results. Given the overlap in autism features with genetic abnormalities known to be associated with imprinting, one possible reason for lack of consistency would be the influence of parent-of-origin effects that may mask the ability to detect linkage and association.

Methods and Findings

We have performed a genome-wide linkage scan that accounts for potential parent-of-origin effects using 16,311 SNPs among families from the Autism Genetic Resource Exchange (AGRE) and the National Institute of Mental Health (NIMH) autism repository. We report parametric (GH, Genehunter) and allele-sharing linkage (Aspex) results using a broad spectrum disorder case definition. Paternal-origin genome-wide statistically significant linkage was observed on chromosomes 4 (LOD_GH = 3.79, empirical p<0.005 and LOD_Aspex = 2.96, p = 0.008), 15 (LOD_GH = 3.09, empirical p<0.005 and LOD_Aspex = 3.62, empirical p = 0.003) and 20 (LOD_GH = 3.36, empirical p<0.005 and LOD_Aspex = 3.38, empirical p = 0.006).

Conclusions

These regions may harbor imprinted sites associated with the development of autism and offer fruitful domains for molecular investigation into the role of epigenetic mechanisms in autism.     

Adiposity significantly modifies genetic risk for dyslipidemia

Recent genome-wide association studies have identified multiple loci robustly associated with plasma lipids, which also contribute to extreme lipid phenotypes. However, these common genetic variants explain <12% of variation in lipid traits. Adiposity is also an important determinant of plasma lipoproteins, particularly plasma TGs and HDL cholesterol (HDLc) concentrations. Thus, interactions between genes and clinical phenotypes may contribute to this unexplained heritability. We have applied a weighted genetic risk score (GRS) for both plasma TGs and HDLc in two large cohorts at the extremes of BMI. Both BMI and GRS were strongly associated with these lipid traits. A significant interaction between obese/lean status and GRS was noted for each of TG (P_Interaction = 2.87 × 10⁻⁴) and HDLc (P_Interaction = 1.05 × 10⁻³). These interactions were largely driven by SNPs tagging APOA5, glucokinase receptor (GCKR), and LPL for TG, and cholesteryl ester transfer protein (CETP), GalNAc-transferase (GALNT2), endothelial lipase (LIPG), and phospholipid transfer protein (PLTP) for HDLc. In contrast, the GRS_{LDL cholesterol} × adiposity interaction was not significant. Sexual dimorphism was evident for the GRS_HDL on HDLc in obese (P_Interaction = 0.016) but not lean subjects. SNP by BMI interactions may provide biological insight into specific genetic associations and missing heritability.     

Sex-stratified Genome-wide Association Studies Including 270,000 Individuals Show Sexual Dimorphism in Genetic Loci for Anthropometric Traits

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10⁻⁸), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.     

Pleiotropy among Common Genetic Loci Identified for Cardiometabolic Disorders and C-Reactive Protein

Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10^-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes.     

Variations in the FRA10AC1 Fragile Site and 15q21 Are Associated with Cerebrospinal Fluid Aβ1-42 Level

Proteolytic fragments of amyloid and post-translational modification of tau species in Cerebrospinal fluid (CSF) as well as cerebral amyloid deposition are important biomarkers for Alzheimer’s Disease. We conducted genome-wide association study to identify genetic factors influencing CSF biomarker level, cerebral amyloid deposition, and disease progression. The genome-wide association study was performed via a meta-analysis of two non-overlapping discovery sample sets to identify genetic variants other than APOE ε4 predictive of the CSF biomarker level (Aβ_1–42, t-Tau, p-Tau_181P, t-Tau:Aβ_1–42 ratio, and p-Tau_181P:Aβ_1–42 ratio) in patients enrolled in the Alzheimer’s Disease Neuroimaging Initiative (ADNI) study. Loci passing a genome-wide significance threshold of P < 5 x 10⁻⁸ were followed-up for replication in an independent sample set. We also performed joint meta-analysis of both discovery sample sets together with the replication sample set. In the discovery phase, we identified variants in FRA10AC1 associated with CSF Aβ_1–42 level passing the genome-wide significance threshold (directly genotyped SNV rs10509663 P _FE = 1.1 x 10⁻⁹, imputed SNV rs116953792 P _FE = 3.5 x 10⁻¹⁰), rs116953792 (P _one-sided = 0.04) achieved replication. This association became stronger in the joint meta-analysis (directly genotyped SNV rs10509663 P _FE = 1.7 x 10⁻⁹, imputed SNV rs116953792 P _FE = 7.6 x 10⁻¹¹). Additionally, we identified locus 15q21 (imputed SNV rs1503351 P _FE = 4.0 x 10⁻⁸) associated with CSF Aβ_1–42 level. No other variants passed the genome-wide significance threshold for other CSF biomarkers in either the discovery sample sets or joint analysis. Gene set enrichment analyses suggested that targeted genes mediated by miR-33, miR-146, and miR-193 were enriched in various GWAS analyses. This finding is particularly important because CSF biomarkers confer disease susceptibility and may be predictive of the likelihood of disease progression in Alzheimer’s Disease.