Genetic Variants Associated with Serum Thyroid Stimulating Hormone (TSH) Levels in European Americans and African Americans from the eMERGE Network

Thyroid stimulating hormone (TSH) hormone levels are normally tightly regulated within an individual; thus, relatively small variations may indicate thyroid disease. Genome-wide association studies (GWAS) have identified variants in PDE8B and FOXE1 that are associated with TSH levels. However, prior studies lacked racial/ethnic diversity, limiting the generalization of these findings to individuals of non-European ethnicities. The Electronic Medical Records and Genomics (eMERGE) Network is a collaboration across institutions with biobanks linked to electronic medical records (EMRs). The eMERGE Network uses EMR-derived phenotypes to perform GWAS in diverse populations for a variety of phenotypes. In this report, we identified serum TSH levels from 4,501 European American and 351 African American euthyroid individuals in the eMERGE Network with existing GWAS data. Tests of association were performed using linear regression and adjusted for age, sex, body mass index (BMI), and principal components, assuming an additive genetic model. Our results replicate the known association of PDE8B with serum TSH levels in European Americans (rs2046045 p = 1.85×10⁻¹⁷, β = 0.09). FOXE1 variants, associated with hypothyroidism, were not genome-wide significant (rs10759944: p = 1.08×10⁻⁶, β = −0.05). No SNPs reached genome-wide significance in African Americans. However, multiple known associations with TSH levels in European ancestry were nominally significant in African Americans, including PDE8B (rs2046045 p = 0.03, β = −0.09), VEGFA (rs11755845 p = 0.01, β = −0.13), and NFIA (rs334699 p = 1.50×10⁻³, β = −0.17). We found little evidence that SNPs previously associated with other thyroid-related disorders were associated with serum TSH levels in this study. These results support the previously reported association between PDE8B and serum TSH levels in European Americans and emphasize the need for additional genetic studies in more diverse populations.     

The p53 Codon 72 PRO/PRO Genotype May Be Associated with Initial Central Visual Field Defects in Caucasians with Primary Open Angle Glaucoma

Background

Loss of vision in glaucoma is due to apoptotic retinal ganglion cell loss. While p53 modulates apoptosis, gene association studies between p53 variants and glaucoma have been inconsistent. In this study we evaluate the association between a p53 variant functionally known to influence apoptosis (codon 72 Pro/Arg) and the subset of primary open angle glaucoma (POAG) patients with early loss of central visual field.

Methods

Genotypes for the p53 codon 72 polymorphism (Pro/Arg) were obtained for 264 POAG patients and 400 controls from the U.S. and in replication studies for 308 POAG patients and 178 controls from Australia (GIST). The glaucoma patients were divided into two groups according to location of initial visual field defect (either paracentral or peripheral). All cases and controls were Caucasian with European ancestry.

Results

The p53-PRO/PRO genotype was more frequent in the U.S. POAG patients with early visual field defects in the paracentral regions compared with those in the peripheral regions or control group (p = 2.7×10⁻⁵). We replicated this finding in the GIST cohort (p  = 7.3×10⁻³, and in the pooled sample (p = 6.6×10⁻⁷) and in a meta-analysis of both the US and GIST datasets (1.3×10⁻⁶, OR 2.17 (1.58–2.98 for the PRO allele).

Conclusions

These results suggest that the p53 codon 72 PRO/PRO genotype is potentially associated with early paracentral visual field defects in primary open-angle glaucoma patients.     

Mitochondrial Genome Analysis of Primary Open Angle Glaucoma Patients

Primary open angle glaucoma (POAG) is a multi-factorial optic disc neuropathy characterized by accelerating damage of the retinal ganglion cells and atrophy of the optic nerve head. The vulnerability of the optic nerve damage leading to POAG has been postulated to result from oxidative stress and mitochondrial dysfunction. In this study, we investigated the possible involvement of the mitochondrial genomic variants in 101 patients and 71 controls by direct sequencing of the entire mitochondrial genome. The number of variable positions in the mtDNA with respect to the revised Cambridge Reference Sequence (rCRS), have been designated “Segregating Sites”. The segregating sites present only in the patients or controls have been designated “Unique Segregating Sites (USS)”. The population mutation rate (θ = 4N_eμ) as estimated by Watterson’s θ (θ_w), considering only the USS, was significantly higher among the patients (p = 9.8×10⁻¹⁵) compared to controls. The difference in θ_w and the number of USS were more pronounced when restricted to the coding region (p<1.31×10⁻²¹ and p = 0.006607, respectively). Further analysis of the region revealed non-synonymous variations were significantly higher in Complex I among the patients (p = 0.0053). Similar trends were retained when USS was considered only within complex I (frequency 0.49 vs 0.31 with p<0.0001 and mutation rate p-value <1.49×10⁻⁴³) and ND5 within its gene cluster (frequency 0.47 vs 0.23 with p<0.0001 and mutation rate p-value <4.42×10⁻⁴⁷). ND5 is involved in the proton pumping mechanism. Incidentally, glaucomatous trabecular meshwork cells have been reported to be more sensitive to inhibition of complex I activity. Thus mutations in ND5, expected to inhibit complex I activity, could lead to generation of oxidative stress and favor glaucomatous condition.     

The PSEN1, p.E318G Variant Increases the Risk of Alzheimer's Disease in APOE-ε4 Carriers

The primary constituents of plaques (Aβ42/Aβ40) and neurofibrillary tangles (tau and phosphorylated forms of tau [ptau]) are the current leading diagnostic and prognostic cerebrospinal fluid (CSF) biomarkers for AD. In this study, we performed deep sequencing of APP, PSEN1, PSEN2, GRN, APOE and MAPT genes in individuals with extreme CSF Aβ42, tau, or ptau levels. One known pathogenic mutation (PSEN1 p.A426P), four high-risk variants for AD (APOE p.L46P, MAPT p.A152T, PSEN2 p.R62H and p.R71W) and nine novel variants were identified. Surprisingly, a coding variant in PSEN1, p.E318G (rs17125721-G) exhibited a significant association with high CSF tau (p = 9.2×10⁻⁴) and ptau (p = 1.8×10⁻³) levels. The association of the p.E318G variant with Aβ deposition was observed in APOE-ε4 allele carriers. Furthermore, we found that in a large case-control series (n = 5,161) individuals who are APOE-ε4 carriers and carry the p.E318G variant are at a risk of developing AD (OR = 10.7, 95% CI = 4.7–24.6) that is similar to APOE-ε4 homozygous (OR = 9.9, 95% CI = 7.2.9–13.6), and double the risk for APOE-ε4 carriers that do not carry p.E318G (OR = 3.9, 95% CI = 3.4–4.4). The p.E318G variant is present in 5.3% (n = 30) of the families from a large clinical series of LOAD families (n = 565) and exhibited a higher frequency in familial LOAD (MAF = 2.5%) than in sporadic LOAD (MAF = 1.6%) (p = 0.02). Additionally, we found that in the presence of at least one APOE-ε4 allele, p.E318G is associated with more Aβ plaques and faster cognitive decline. We demonstrate that the effect of PSEN1, p.E318G on AD susceptibility is largely dependent on an interaction with APOE-ε4 and mediated by an increased burden of Aβ deposition.     

Common Variants of FTO Are Associated with Childhood Obesity in a Cross-Sectional Study of 3,126 Urban Indian Children

Background

FTO variants are robustly associated with obesity and related traits in many population and shown to have variable impact during life course. Although studies have shown association of FTO variants with adiposity in adult Indian, its association in Indian children is yet to be confirmed.

Methods

Here we examined association of FTO variants (rs9939609 and rs8050136) with obesity and related anthropometric and biochemical traits in 3,126 Indian children (aged 11–17 years) including 2,230 normal-weight and 896 over-weight/obese children. We also compared effects observed in the present study with that observed in previous studies on South Asian adults and children of other ethnic groups.

Results

The variant rs9939609 showed significant association with risk of obesity [OR = 1.21, P = 2.5×10⁻³] and its measures BMI, weight, waist circumference and hip circumference [β range = 0.11 to 0.14 Z-score units; P range = 1.3×10⁻⁴ to 1.6×10⁻⁷] in children. The observed effect sizes in Indian children were similar to those reported for European children. Variant rs9939609 explained 0.88% of BMI variance in Indian children. The effect sizes of rs9939609 on BMI and WC were ∼2 fold higher in children than adults. Interestingly rs9939609 was also associated with serum levels of thyroid stimulating hormone (TSH) [β = 0.10 Z-score, P = 5.8×10⁻³]. The other variant rs8050136 was in strong linkage disequilibrium with rs9939609 (r² = 0.97) and provided similar association results.

Conclusion

The study provides first report of association of FTO variants with obesity and related anthropometric traits in Indian children with higher impact in children compared to adults. We also demonstrated association of FTO variant with serum levels of TSH, indicating putative influence of FTO in hypothalamic-pituitary-thyroid axis.     

A Genome-Wide Association Study of Optic Disc Parameters

The optic nerve head is involved in many ophthalmic disorders, including common diseases such as myopia and open-angle glaucoma. Two of the most important parameters are the size of the optic disc area and the vertical cup-disc ratio (VCDR). Both are highly heritable but genetically largely undetermined. We performed a meta-analysis of genome-wide association (GWA) data to identify genetic variants associated with optic disc area and VCDR. The gene discovery included 7,360 unrelated individuals from the population-based Rotterdam Study I and Rotterdam Study II cohorts. These cohorts revealed two genome-wide significant loci for optic disc area, rs1192415 on chromosome 1p22 (p = 6.72×10⁻¹⁹) within 117 kb of the CDC7 gene and rs1900004 on chromosome 10q21.3-q22.1 (p = 2.67×10⁻³³) within 10 kb of the ATOH7 gene. They revealed two genome-wide significant loci for VCDR, rs1063192 on chromosome 9p21 (p = 6.15×10⁻¹¹) in the CDKN2B gene and rs10483727 on chromosome 14q22.3-q23 (p = 2.93×10⁻¹⁰) within 40 kbp of the SIX1 gene. Findings were replicated in two independent Dutch cohorts (Rotterdam Study III and Erasmus Rucphen Family study; N = 3,612), and the TwinsUK cohort (N = 843). Meta-analysis with the replication cohorts confirmed the four loci and revealed a third locus at 16q12.1 associated with optic disc area, and four other loci at 11q13, 13q13, 17q23 (borderline significant), and 22q12.1 for VCDR. ATOH7 was also associated with VCDR independent of optic disc area. Three of the loci were marginally associated with open-angle glaucoma. The protein pathways in which the loci of optic disc area are involved overlap with those identified for VCDR, suggesting a common genetic origin.     

A Novel MMP12 Locus Is Associated with Large Artery Atherosclerotic Stroke Using a Genome-Wide Age-at-Onset Informed Approach

Genome-wide association studies (GWAS) have begun to identify the common genetic component to ischaemic stroke (IS). However, IS has considerable phenotypic heterogeneity. Where clinical covariates explain a large fraction of disease risk, covariate informed designs can increase power to detect associations. As prevalence rates in IS are markedly affected by age, and younger onset cases may have higher genetic predisposition, we investigated whether an age-at-onset informed approach could detect novel associations with IS and its subtypes; cardioembolic (CE), large artery atherosclerosis (LAA) and small vessel disease (SVD) in 6,778 cases of European ancestry and 12,095 ancestry-matched controls. Regression analysis to identify SNP associations was performed on posterior liabilities after conditioning on age-at-onset and affection status. We sought further evidence of an association with LAA in 1,881 cases and 50,817 controls, and examined mRNA expression levels of the nearby genes in atherosclerotic carotid artery plaques. Secondly, we performed permutation analyses to evaluate the extent to which age-at-onset informed analysis improves significance for novel loci. We identified a novel association with an MMP12 locus in LAA (rs660599; p = 2.5×10⁻⁷), with independent replication in a second population (p = 0.0048, OR(95% CI) = 1.18(1.05–1.32); meta-analysis p = 2.6×10⁻⁸). The nearby gene, MMP12, was significantly overexpressed in carotid plaques compared to atherosclerosis-free control arteries (p = 1.2×10⁻¹⁵; fold change = 335.6). Permutation analyses demonstrated improved significance for associations when accounting for age-at-onset in all four stroke phenotypes (p<0.001). Our results show that a covariate-informed design, by adjusting for age-at-onset of stroke, can detect variants not identified by conventional GWAS.     

t-10, c-12 CLA Dietary Supplementation Inhibits Atherosclerotic Lesion Development Despite Adverse Cardiovascular and Hepatic Metabolic Marker Profiles

Animal and human studies have indicated that fatty acids such as the conjugated linoleic acids (CLA) found in milk could potentially alter the risk of developing metabolic disorders including diabetes and cardiovascular disease (CVD). Using susceptible rodent models (apoE^−/− and LDLr^−/− mice) we investigated the interrelationship between mouse strain, dietary conjugated linoleic acids and metabolic markers of CVD. Despite an adverse metabolic risk profile, atherosclerosis (measured directly by lesion area), was significantly reduced with t-10, c-12 CLA and mixed isomer CLA (Mix) supplementation in both apoE^−/− (p<0.05, n = 11) and LDLr^−/− mice (p<0.01, n = 10). Principal component analysis was utilized to delineate the influence of multiple plasma and tissue metabolites on the development of atherosclerosis. Group clustering by dietary supplementation was evident, with the t-10, c-12 CLA supplemented animals having distinct patterns, suggestive of hepatic insulin resistance, regardless of mouse strain. The effect of CLA supplementation on hepatic lipid and fatty acid composition was explored in the LDLr^−/− strain. Dietary supplementation with t-10, c-12 CLA significantly increased liver weight (p<0.05, n = 10), triglyceride (p<0.01, n = 10) and cholesterol ester content (p<0.01, n = 10). Furthermore, t-10, c-12 CLA also increased the ratio of 18∶1 to 18∶0 fatty acid in the liver suggesting an increase in the activity of stearoyl-CoA desaturase. Changes in plasma adiponectin and liver weight with t-10, c-12 CLA supplementation were evident within 3 weeks of initiation of the diet. These observations provide evidence that the individual CLA isomers have divergent mechanisms of action and that t-10, c-12 CLA rapidly changes plasma and liver markers of metabolic syndrome, despite evidence of reduction in atherosclerosis.     

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.     

The NOD2 p.Leu1007fsX1008 Mutation (rs2066847) Is a Stronger Predictor of the Clinical Course of Crohn's Disease than the FOXO3A Intron Variant rs12212067

Background

Very recently, a sub-analysis of genome-wide association scans revealed that the non-coding single nucleotide polymorphism (SNP) rs12212067 in the FOXO3A gene is associated with a milder course of Crohn''s disease (CD) (Cell 2013;155:57–69). The aim of our study was to evaluate the clinical value of the SNP rs12212067 in predicting the severity of CD by correlating CD patient genotype status with the most relevant complications of CD such as stenoses, fistulas, and CD-related surgery.

Methodology/Principal Findings

We genotyped 550 CD patients for rs12212067 (FOXO3A) and the three common CD-associated NOD2 mutations rs2066844, rs2066847, and rs2066847 and performed genotype-phenotype analyses.

Results

No significant phenotypic differences were found between the wild-type genotype TT of the FOXO3A SNP rs12212067 and the minor genotypes TG and GG independently from NOD2 variants. The allele frequency of the minor G allele was 12.7%. Age at diagnosis, disease duration, body mass index, surgery rate, stenoses, fistula, need for immunosuppressive therapy, and disease course were not significantly different. In contrast, the NOD2 mutant p.Leu1007fsX1008 (rs2066847) was highly associated with penetrating CD (p = 0.01), the development of fistulas (p = 0.01) and stenoses (p = 0.01), and ileal disease localization (p = 0.03). Importantly, the NOD2 SNP rs2066847 was a strong separator between an aggressive and a mild course of CD (p = 2.99×10⁻⁵), while the FOXO3A SNP rs12212067 did not separate between mild and aggressive CD behavior in our cohort (p = 0.35). 96.2% of the homozygous NOD2 p.Leu1007fsX1008 carriers had an aggressive disease behavior compared to 69.3% of the patients with the NOD2 wild-type genotype (p = 0.007).

Conclusion/Significance

In clinical practice, the NOD2 variant p.Leu1007fsX1008 (rs2066847), in particular in homozygous form, is a much stronger marker for a severe clinical phenotype than the FOXO3A rs12212067 SNP for a mild disease course on an individual patient level despite its important impact on the inflammatory response of monocytes.     

Distribution and Medical Impact of Loss-of-Function Variants in the Finnish Founder Population

Exome sequencing studies in complex diseases are challenged by the allelic heterogeneity, large number and modest effect sizes of associated variants on disease risk and the presence of large numbers of neutral variants, even in phenotypically relevant genes. Isolated populations with recent bottlenecks offer advantages for studying rare variants in complex diseases as they have deleterious variants that are present at higher frequencies as well as a substantial reduction in rare neutral variation. To explore the potential of the Finnish founder population for studying low-frequency (0.5–5%) variants in complex diseases, we compared exome sequence data on 3,000 Finns to the same number of non-Finnish Europeans and discovered that, despite having fewer variable sites overall, the average Finn has more low-frequency loss-of-function variants and complete gene knockouts. We then used several well-characterized Finnish population cohorts to study the phenotypic effects of 83 enriched loss-of-function variants across 60 phenotypes in 36,262 Finns. Using a deep set of quantitative traits collected on these cohorts, we show 5 associations (p<5×10⁻⁸) including splice variants in LPA that lowered plasma lipoprotein(a) levels (P = 1.5×10⁻¹¹⁷). Through accessing the national medical records of these participants, we evaluate the LPA finding via Mendelian randomization and confirm that these splice variants confer protection from cardiovascular disease (OR = 0.84, P = 3×10⁻⁴), demonstrating for the first time the correlation between very low levels of LPA in humans with potential therapeutic implications for cardiovascular diseases. More generally, this study articulates substantial advantages for studying the role of rare variation in complex phenotypes in founder populations like the Finns and by combining a unique population genetic history with data from large population cohorts and centralized research access to National Health Registers.     

Association of HK2 and NCK2 with Normal Tension Glaucoma in the Japanese Population

Although family studies and genome-wide association studies have shown that genetic factors play a role in glaucoma, it has been difficult to identify the specific genetic variants involved. We tested 669 single nucleotide polymorphisms (SNPs) from the region of chromosome 2 that includes the GLC1B glaucoma locus for association with primary open-angle glaucoma (POAG) and normal tension glaucoma (NTG) in the Japanese population. We performed a two-stage case-control study. The first cohort consisted of 123 POAG cases, 121 NTG cases and 120 controls: the second cohort consisted of 187 POAG cases, 286 NTG cases, and 271 controls. Out of six SNPs showing significant association with POAG in the first round screening, seven SNPs were tested in the second round. Rs678350 in the HK2 gene coding sequence showed significant allelic (p = 0.0027 in Stage Two, 2.7XE-4 in meta-analysis) association with POAG, and significant allelic (p = 4.7XE-4 in Stage Two, 1.0XE-5 in meta-analysis) association with NTG. Although alleles in the TMEM182 gene did not show significant association with glaucoma in the second round, subjects with the A/A allele in TMEM182 rs869833 showed worse visual field mean deviation (p = 0.01). Even though rs2033008 in the NCK2 gene coding sequence did not show significant association in the first round, it had previously shown association with NTG so it was tested for association with NTG in round 2 (p = 0.0053 in Stage Two). Immunohistochemistry showed that both HK2 and NCK2 are expressed in the retinal ganglion cell layer. Once multi-testing was taken into account, only HK2 showed significant association with POAG and NTG in Stage Two. Our data also support previous reports of NCK2 association with NTG, and raise questions about what role TMEM182 might play in phenotypic variability. Our data suggest that HK2 may play an important role in NTG in the Japanese population.     

Genome-Wide Association Study of Maternal and Inherited Loci for Conotruncal Heart Defects

Conotruncal and related heart defects (CTDs) are a group of serious and relatively common birth defects. Although both maternal and inherited genotypes are thought to play a role in the etiology of CTDs, few specific genetic risk factors have been identified. To determine whether common variants acting through the genotype of the mother (e.g. via an in utero effect) or the case are associated with CTDs, we conducted a genome-wide association study of 750 CTD case-parent triads, with follow-up analyses in 358 independent triads. Log-linear analyses were used to assess the association of CTDs with the genotypes of both the mother and case. No association achieved genomewide significance in either the discovery or combined (discovery+follow-up) samples. However, three loci with p-values suggestive of association (p<10⁻⁵) in the discovery sample had p-values <0.05 in the follow-up sample and p-values in the combined data that were lower than in the discovery sample. These included suggestive association with an inherited intergenic variant at 20p12.3 (rs6140038, combined p = 1.0×10⁻⁵) and an inherited intronic variant in KCNJ4 at 22q13.1 (rs2267386, combined p = 9.8×10⁻⁶), as well as with a maternal variant in SLC22A24 at 11q12.3 (rs11231379, combined p = 4.2×10⁻⁶). These observations suggest novel candidate loci for CTDs, including loci that appear to be associated with the risk of CTDs via the maternal genotype, but further studies are needed to confirm these associations.     

Polymorphisms of CYP51A1 from Cholesterol Synthesis: Associations with Birth Weight and Maternal Lipid Levels and Impact on CYP51 Protein Structure

We investigated the housekeeping cytochrome P450 CYP51A1 encoding lanosterol 14α-demethylase from cholesterol synthesis that was so far not directly linked to human disorders. By direct sequencing of CYP51A1 in 188 women with spontaneous preterm delivery and 188 unrelated preterm infants (gestational age <37 weeks) we identified 22 variants where 10 are novel and rare. In infants there were two novel CYP51A1 variants where damaging effects of p.Tyr145Asp from the substrate recognition region, but not p.Asn193Asp, were predicted by PolyPhen2 and SIFT. This was confirmed by molecular modeling showing that Tyr145Asp substitution results in changed electrostatic potential of the CYP51 protein surface and lengthened distance to the heme which prevents hydrogen bonding. The CYP51 Tyr145Asp mutation is rare and thus very interesting for further structure/function relationship studies. From the 12 identified known variants rs6465348 was chosen for family based association studies due to its high minor allele frequency. Interestingly, this CYP51A1 common variant associates with small for gestational age weight in newborns (p = 0.028) and lower blood total cholesterol and low density lipoprotein cholesterol levels in mothers in 2nd trimester of pregnancy (p = 0.042 and p = 0.046 respectively). Our results indicate a new link between a cholesterol synthesis gene CYP51A1 and pregnancy pathologies.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

Genome-Wide Association Studies Identify Two Novel BMP15 Mutations Responsible for an Atypical Hyperprolificacy Phenotype in Sheep

Some sheep breeds are naturally prolific, and they are very informative for the studies of reproductive genetics and physiology. Major genes increasing litter size (LS) and ovulation rate (OR) were suspected in the French Grivette and the Polish Olkuska sheep populations, respectively. To identify genetic variants responsible for the highly prolific phenotype in these two breeds, genome-wide association studies (GWAS) followed by complementary genetic and functional analyses were performed. Highly prolific ewes (cases) and normal prolific ewes (controls) from each breed were genotyped using the Illumina OvineSNP50 Genotyping Beadchip. In both populations, an X chromosome region, close to the BMP15 gene, harbored clusters of markers with suggestive evidence of association at significance levels between 1E⁻⁰⁵ and 1E⁻⁰⁷. The BMP15 candidate gene was then sequenced, and two novel non-conservative mutations called FecX^Gr and FecX^O were identified in the Grivette and Olkuska breeds, respectively. The two mutations were associated with the highly prolific phenotype (p_FecX^Gr = 5.98E⁻⁰⁶ and p_FecX^O = 2.55E⁻⁰⁸). Homozygous ewes for the mutated allele showed a significantly increased prolificacy (FecX^Gr/FecX^Gr, LS = 2.50±0.65 versus FecX⁺/FecX^Gr, LS = 1.93±0.42, p<1E⁻⁰³ and FecX^O/FecX^O, OR = 3.28±0.85 versus FecX⁺/FecX^O, OR = 2.02±0.47, p<1E⁻⁰³). Both mutations are located in very well conserved motifs of the protein and altered the BMP15 signaling activity in vitro using a BMP-responsive luciferase test in COV434 granulosa cells. Thus, we have identified two novel mutations in the BMP15 gene associated with increased LS and OR. Notably, homozygous FecX^Gr/FecX^Gr Grivette and homozygous FecX^O/FecX^O Olkuska ewes are hyperprolific in striking contrast with the sterility exhibited by all other known homozygous BMP15 mutations. Our results bring new insights into the key role played by the BMP15 protein in ovarian function and could contribute to a better understanding of the pathogenesis of women′s fertility disorders.     

Pulmonary and Cardiac Function in Asymptomatic Obese Subjects and Changes following a Structured Weight Reduction Program: A Prospective Observational Study

Background

The prevalence of obesity is rising. Obesity can lead to cardiovascular and ventilatory complications through multiple mechanisms. Cardiac and pulmonary function in asymptomatic subjects and the effect of structured dietary programs on cardiac and pulmonary function is unclear.

Objective

To determine lung and cardiac function in asymptomatic obese adults and to evaluate whether weight loss positively affects functional parameters.

Methods

We prospectively evaluated bodyplethysmographic and echocardiographic data in asymptomatic subjects undergoing a structured one-year weight reduction program.

Results

74 subjects (32 male, 42 female; mean age 42±12 years) with an average BMI 42.5±7.9, body weight 123.7±24.9 kg were enrolled. Body weight correlated negatively with vital capacity (R = −0.42, p<0.001), FEV1 (R = −0.497, p<0.001) and positively with P 0.1 (R = 0.32, p = 0.02) and myocardial mass (R = 0.419, p = 0.002). After 4 months the study subjects had significantly reduced their body weight (−26.0±11.8 kg) and BMI (−8.9±3.8) associated with a significant improvement of lung function (absolute changes: vital capacity +5.5±7.5% pred., p<0.001; FEV1+9.8±8.3% pred., p<0.001, ITGV+16.4±16.0% pred., p<0.001, SR tot −17.4±41.5% pred., p<0.01). Moreover, P0.1/Pimax decreased to 47.7% (p<0.01) indicating a decreased respiratory load. The change of FEV1 correlated significantly with the change of body weight (R = −0.31, p = 0.03). Echocardiography demonstrated reduced myocardial wall thickness (−0.08±0.2 cm, p = 0.02) and improved left ventricular myocardial performance index (−0.16±0.35, p = 0.02). Mitral annular plane systolic excursion (+0.14, p = 0.03) and pulmonary outflow acceleration time (AT +26.65±41.3 ms, p = 0.001) increased.

Conclusion

Even in asymptomatic individuals obesity is associated with abnormalities in pulmonary and cardiac function and increased myocardial mass. All the abnormalities can be reversed by a weight reduction program.     

CX3CR1 Is a Modifying Gene of Survival and Progression in Amyotrophic Lateral Sclerosis

The objective of this study was to investigate the association of functional variants of the human CX3CR1 gene (Fractalkine receptor) with the risk of Amyotrophic Lateral Sclerosis (ALS), the survival and the progression rate of the disease symptoms in a Spanish ALS cohort. 187 ALS patients (142 sporadic [sALS] and 45 familial) and 378 controls were recruited. We investigated CX3CR1 V249I (rs3732379) and T280M (rs3732378) genotypes and their haplotypes as predictors of survival, the progression rate of the symptoms (as measured by ALSFRS-R and FVC decline) and the risk of suffering ALS disease. The results indicated that sALS patients with CX3CR1 249^I/I or 249^V/I genotypes presented a shorter survival time (42.27±4.90) than patients with 249^V/V genotype (67.65±7.42; diff −25.49 months 95%CI [−42.79,−8.18]; p = 0.004; adj-p = 0.018). The survival time was shorter in sALS patients with spinal topography and CX3CR1 249^I alleles (diff = −29.78 months; 95%CI [−49.42,−10.14]; p = 0.003). The same effects were also observed in the spinal sALS patients with 249^I–280^M haplotype (diff = −27.02 months; 95%CI [−49.57, −4.48]; p = 0.019). In the sALS group, the CX3CR1 249^I variant was associated with a faster progression of the disease symptoms (OR = 2.58; 95IC% [1.32, 5.07]; p = 0.006; adj-p = 0.027). There was no evidence for association of these two CX3CR1 variants with ALS disease risk. The association evidenced herein is clinically relevant and indicates that CX3CR1 could be a disease-modifying gene in sALS. The progression rate of the disease''s symptoms and the survival time is affected in patients with one or two copies of the CX3CR1 249I allele. The CX3CR1 is the most potent ALS survival genetic factor reported to date. These results reinforce the role of the immune system in ALS pathogenesis.