Genome-Wide Association Analysis and Genomic Prediction of Thyroglobulin Plasma Levels

Thyroglobulin (Tg) is an iodoglycoprotein produced by thyroid follicular cells which acts as an essential substrate for thyroid hormone synthesis. To date, only one genome-wide association study (GWAS) of plasma Tg levels has been performed by our research group. Utilizing recent advancements in computation and modeling, we apply a Bayesian approach to the probabilistic inference of the genetic architecture of Tg. We fitted a Bayesian sparse linear mixed model (BSLMM) and a frequentist linear mixed model (LMM) of 7,289,083 variants in 1096 healthy European-ancestry participants of the Croatian Biobank. Meta-analysis with two independent cohorts (total n = 2109) identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) within the ST6GAL1 gene ([Formula: see text]). BSLMM revealed additional association signals on chromosomes 1, 8, 10, and 14. For ST6GAL1 and the newly uncovered genes, we provide physiological and pathophysiological explanations of how their expression could be associated with variations in plasma Tg levels. We found that the SNP-heritability of Tg is 17% and that 52% of this variation is due to a small number of 16 variants that have a major effect on Tg levels. Our results suggest that the genetic architecture of plasma Tg is not polygenic, but influenced by a few genes with major effects

SNP and Haplotype Regional Heritability Mapping (SNHap-RHM):Joint Mapping of Common and Rare Variation Affecting Complex Traits

We describe a genome-wide analytical approach, SNP and Haplotype Regional Heritability Mapping (SNHap-RHM), that provides regional estimates of the heritability across locally defined regions in the genome. This approach utilises relationship matrices that are based on sharing of SNP and haplotype alleles at local haplotype blocks delimited by recombination boundaries in the genome. We implemented the approach on simulated data and show that the haplotype-based regional GRMs capture variation that is complementary to that captured by SNP-based regional GRMs, and thus justifying the fitting of the two GRMs jointly in a single analysis (SNHap-RHM). SNHap-RHM captures regions in the genome contributing to the phenotypic variation that existing genome-wide analysis methods may fail to capture. We further demonstrate that there are real benefits to be gained from this approach by applying it to real data from about 20,000 individuals from the Generation Scotland: Scottish Family Health Study. We analysed height and major depressive disorder (MDD). We identified seven genomic regions that are genome-wide significant for height, and three regions significant at a suggestive threshold (p-value < 1 × 10(−5)) for MDD. These significant regions have genes mapped to within 400 kb of them. The genes mapped for height have been reported to be associated with height in humans. Similarly, those mapped for MDD have been reported to be associated with major depressive disorder and other psychiatry phenotypes. The results show that SNHap-RHM presents an exciting new opportunity to analyse complex traits by allowing the joint mapping of novel genomic regions tagged by either SNPs or haplotypes, potentially leading to the recovery of some of the “missing” heritability

Characterisation of an inflammation-related epigenetic score and its association with cognitive ability

BACKGROUND: Chronic systemic inflammation has been associated with incident dementia, but its association with age-related cognitive decline is less clear. The acute responses of many inflammatory biomarkers mean they may provide an unreliable picture of the chronicity of inflammation. Recently, a large-scale epigenome-wide association study identified DNA methylation correlates of C-reactive protein (CRP)-a widely used acute-phase inflammatory biomarker. DNA methylation is thought to be relatively stable in the short term, marking it as a potentially useful signature of exposure. METHODS: We utilise a DNA methylation-based score for CRP and investigate its trajectories with age, and associations with cognitive ability in comparison with serum CRP and a genetic CRP score in a longitudinal study of older adults (n = 889) and a large, cross-sectional cohort (n = 7028). RESULTS: We identified no homogeneous trajectories of serum CRP with age across the cohorts, whereas the epigenetic CRP score was consistently found to increase with age (standardised β = 0.07 and 0.01) and to do so more rapidly in males compared to females. Additionally, the epigenetic CRP score had higher test-retest reliability compared to serum CRP, indicating its enhanced temporal stability. Higher serum CRP was not found to be associated with poorer cognitive ability (standardised β = - 0.08 and - 0.05); however, a consistent negative association was identified between cognitive ability and the epigenetic CRP score in both cohorts (standardised β = - 0.15 and - 0.08). CONCLUSIONS: An epigenetic proxy of CRP may provide a more reliable signature of chronic inflammation, allowing for more accurate stratification of individuals, and thus clearer inference of associations with incident health outcomes

Linking protein to phenotype with Mendelian Randomization detects 38 proteins with causal roles in human diseases and traits

To efficiently transform genetic associations into drug targets requires evidence that a particular gene, and its encoded protein, contribute causally to a disease. To achieve this, we employ a three-step proteome-by-phenome Mendelian Randomization (MR) approach. In step one, 154 protein quantitative trait loci (pQTLs) were identified and independently replicated. From these pQTLs, 64 replicated locally-acting variants were used as instrumental variables for proteome-by-phenome MR across 846 traits (step two). When its assumptions are met, proteome-by-phenome MR, is equivalent to simultaneously running many randomized controlled trials. Step 2 yielded 38 proteins that significantly predicted variation in traits and diseases in 509 instances. Step 3 revealed that amongst the 271 instances from GeneAtlas (UK Biobank), 77 showed little evidence of pleiotropy (HEIDI), and 92 evidence of colocalization (eCAVIAR). Results were wide ranging: including, for example, new evidence for a causal role of tyrosine-protein phosphatase non-receptor type substrate 1 (SHPS1; SIRPA) in schizophrenia, and a new finding that intestinal fatty acid binding protein (FABP2) abundance contributes to the pathogenesis of cardiovascular disease. We also demonstrated confirmatory evidence for the causal role of four further proteins (FGF5, IL6R, LPL, LTA) in cardiovascular disease risk

Genome-wide meta-analysis identifies novel loci associated with parathyroid hormone level

Abstract Background Parathyroid hormone (PTH) is one of the principal regulators of calcium homeostasis. Although serum PTH level is mostly accounted by genetic factors, genetic background underlying PTH level is insufficiently known. Therefore, the aim of this study was to identify novel genetic variants associated with PTH levels. Methods We performed GWAS meta-analysis within two genetically isolated Croatian populations followed by replication analysis in a Croatian mainland population and we also combined results across all three analyzed populations. The analyses included 2596 individuals. A total of 7,411,206 variants, imputed using the 1000 Genomes reference panel, were analysed for the association. In addition, a sex-specific GWAS meta-analyses were performed. Results Polymorphisms with the lowest P-values were located on chromosome 4 approximately 84 kb of the 5′ of RASGEF1B gene. The most significant SNP was rs11099476 (P = 1.15 × 10−8). Sex-specific analysis identified genome-wide significant association of the variant rs77178854, located within DPP10 gene in females only (P = 2.21 × 10− 9). There were no genome-wide significant findings in the meta-analysis of males. Conclusions We identified two biologically plausible novel loci associated with PTH levels, providing us with further insights into the genetics of this complex trait

An epigenome-wide association study of sex-specific chronological ageing

Background Advanced age is associated with cognitive and physical decline and is a major risk factor for a multitude of disorders. There is also a gap in life expectancy between males and females. DNA methylation differences have been shown to be associated with both age and sex. Here, we investigate age-by-sex differences in blood-based DNA methylation in an unrelated cohort of 2586 individuals between the ages of 18 and 87 years, with replication in a further 4450 individuals between the ages of 18 and 93 years. Methods Linear regression models were applied, with stringent genome-wide significance thresholds (p < 3.6 x 10(-8)) used in both the discovery and replication data. A second, highly conservative mixed linear model method that better controls the false-positive rate was also applied, using the same genome-wide significance thresholds. Results Using the linear regression method, 52 autosomal and 597 X-linked CpG sites, mapping to 251 unique genes, replicated with concordant effect size directions in the age-by-sex interaction analysis. The site with the greatest difference mapped to GAGE10, an X-linked gene. Here, DNA methylation levels remained stable across the male adult age range (DNA methylation by age r = 0.02) but decreased across female adult age range (DNA methylation by age r = - 0.61). One site (cg23722529) with a significant age-by-sex interaction also had a quantitative trait locus (rs17321482) that is a genome-wide significant variant for prostate cancer. The mixed linear model method identified 11 CpG sites associated with the age-by-sex interaction. Conclusion The majority of differences in age-associated DNA methylation trajectories between sexes are present on the X chromosome. Several of these differences occur within genes that have been implicated in sexually dimorphic traits

Author Correction: Cross-ancestry genome-wide association analysis of corneal thickness strengthens link between complex and Mendelian eye diseases.

Emmanuelle Souzeau, who contributed to analysis of data, was inadvertently omitted from the author list in the originally published version of this Article. This has now been corrected in both the PDF and HTML versions of the Article

An actionable KCNH2 Long QT Syndrome variant detected by sequence and haplotype analysis in a population research cohort

The Viking Health Study Shetland is a population-based research cohort of 2,122 volunteer participants with ancestry from the Shetland Isles in northern Scotland. The high kinship and detailed phenotype data support a range of approaches for associating rare genetic variants, enriched in this isolate population, with quantitative traits and diseases. As an exemplar, the c.1750G > A; p.Gly584Ser variant within the coding sequence of the KCNH2 gene implicated in Long QT Syndrome (LQTS), which occurred once in 500 whole genome sequences from this population, was investigated. Targeted sequencing of the KCNH2 gene in family members of the initial participant confirmed the presence of the sequence variant and identified two further members of the same family pedigree who shared the variant. Investigation of these three related participants for whom single nucleotide polymorphism (SNP) array genotypes were available allowed a unique shared haplotype of 1.22 Mb to be defined around this locus. Searching across the full cohort for this haplotype uncovered two additional apparently unrelated individuals with no known genealogical connection to the original kindred. All five participants with the defined haplotype were shown to share the rare variant by targeted Sanger sequencing. If this result were verified in a healthcare setting, it would be considered clinically actionable, and has been actioned in relatives ascertained independently through clinical presentation. The General Practitioners of four study participants with the rare variant were alerted to the research findings by letters outlining the phenotype (prolonged electrocardiographic QTc interval). A lack of detectable haplotype sharing between c.1750G > A; p.Gly584Ser chromosomes from previously reported individuals from Finland and those in this study from Shetland suggests that this mutation has arisen more than once in human history. This study showcases the potential value of isolate population-based research resources for genomic medicine. It also illustrates some challenges around communication of actionable findings in research participants in this context.Peer reviewe