Genetic Landscape of the ACE2 Coronavirus Receptor.

BACKGROUND: SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood. METHODS: We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data. RESULTS: We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells. CONCLUSIONS: Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

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.

Insights into the genetic basis of retinal detachment.

Retinal detachment (RD) is a serious and common condition, but genetic studies to date have been hampered by the small size of the assembled cohorts. In the UK Biobank data set, where RD was ascertained by self-report or hospital records, genetic correlations between RD and high myopia or cataract operation were, respectively, 0.46 (SE = 0.08) and 0.44 (SE = 0.07). These correlations are consistent with known epidemiological associations. Through meta-analysis of genome-wide association studies using UK Biobank RD cases (N = 3 977) and two cohorts, each comprising ~1 000 clinically ascertained rhegmatogenous RD patients, we uncovered 11 genome-wide significant association signals. These are near or within ZC3H11B, BMP3, COL22A1, DLG5, PLCE1, EFEMP2, TYR, FAT3, TRIM29, COL2A1 and LOXL1. Replication in the 23andMe data set, where RD is self-reported by participants, firmly establishes six RD risk loci: FAT3, COL22A1, TYR, BMP3, ZC3H11B and PLCE1. Based on the genetic associations with eye traits described to date, the first two specifically impact risk of a RD, whereas the last four point to shared aetiologies with macular condition, myopia and glaucoma. Fine-mapping prioritized the lead common missense variant (TYR S192Y) as causal variant at the TYR locus and a small set of credible causal variants at the FAT3 locus. The larger study size presented here, enabled by resources linked to health records or self-report, provides novel insights into RD aetiology and underlying pathological pathways.

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 (p<5×10-8). 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.

Genome-wide meta-analysis identifies novel gender specific loci associated with thyroid antibodies level in Croatians.

Autoimmune thyroid diseases (AITD) are multifactorial endocrine diseases most frequently accompanied by Tg and TPO autoantibodies. Both antibodies have a higher prevalence in females and act under a strong genetic influence. To identify novel variants underlying thyroid antibody levels, we performed GWAS meta-analysis on the plasma levels of TgAb and TPOAb in three Croatian cohorts, as well as gender specific GWAS and a bivariate analysis. No significant association was detected with the level of TgAb and TPOAb in the meta-analysis of GWAS or bivariate results for all individuals. The bivariate analysis in females only revealed a genome-wide significant association for the locus near GRIN3A (rs4457391, P = 7.76 × 10-9). The same locus had borderline association with TPOAb levels in females (rs1935377, P = 8.58 × 10-8). In conclusion, we identified a novel gender specific locus associated with TgAb and TPOAb levels. Our findings provide a novel insight into genetic and gender differences associated with thyroid antibodies.

New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics.

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a heritable disease. Siblings of POAG cases have a ten-fold increased risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup-disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed so far. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG risk.

Genome-wide association meta-analysis for total thyroid hormone levels in Croatian population.

Thyroid hormones (THs) are key regulators of cellular growth, development, and metabolism. The thyroid gland secretes two THs, thyroxine (T4) and triiodothyronine (T3), into the plasma where they are almost all bound reversibly to plasma proteins. Free forms of THs are metabolically active, however, they represent a very small fraction of total TH levels. No genome-wide studies have been performed to date on total TH levels, comprising of protein-bound and free forms of THs. To detect genetic variants associated with total TH levels, we carried out the first GWAS meta-analysis of total T4 levels in 1121 individuals from two Croatian cohorts (Split and Korcula). We also performed GWAS analyses of total T3 levels in 577 individuals and T3/T4 ratio in 571 individuals from the Split cohort. The top association in GWAS meta-analysis of total T4 was detected for an intronic variant within SLC22A9 gene (rs12282281, P = 4.00 × 10-7). Within the same region, a genome-wide significant variant (rs11822642, P = 2.50 × 10-8) for the T3/T4 ratio was identified. SLC22A9 encodes for an organic anion transporter protein expressed predominantly in the liver and belongs to the superfamily of solute carriers (SLC), a large group of transport membrane proteins. The transport of THs across the plasma membrane in peripheral tissues is facilitated by the membrane proteins, and all TH transport proteins known to date belong to the same SLC superfamily as SLC22A9. These results suggest a potential role for SLC22A9 as a novel transporter protein of THs.

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

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.

When do myopia genes have their effect? Comparison of genetic risks between children and adults.

Previous studies have identified many genetic loci for refractive error and myopia. We aimed to investigate the effect of these loci on ocular biometry as a function of age in children, adolescents, and adults. The study population consisted of three age groups identified from the international CREAM consortium: 5,490 individuals aged <10 years; 5,000 aged 10-25 years; and 16,274 aged >25 years. All participants had undergone standard ophthalmic examination including measurements of axial length (AL) and corneal radius (CR). We examined the lead SNP at all 39 currently known genetic loci for refractive error identified from genome-wide association studies (GWAS), as well as a combined genetic risk score (GRS). The beta coefficient for association between SNP genotype or GRS versus AL/CR was compared across the three age groups, adjusting for age, sex, and principal components. Analyses were Bonferroni-corrected. In the age group <10 years, three loci (GJD2, CHRNG, ZIC2) were associated with AL/CR. In the age group 10-25 years, four loci (BMP2, KCNQ5, A2BP1, CACNA1D) were associated; and in adults 20 loci were associated. Association with GRS increased with age; ß = 0.0016 per risk allele (P = 2 × 10-8 ) in <10 years, 0.0033 (P = 5 × 10-15 ) in 10- to 25-year-olds, and 0.0048 (P = 1 × 10-72 ) in adults. Genes with strongest effects (LAMA2, GJD2) had an early effect that increased with age. Our results provide insights on the age span during which myopia genes exert their effect. These insights form the basis for understanding the mechanisms underlying high and pathological myopia.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Colocalization of corneal resistance factor GWAS loci with GTEx e/sQTLs highlights plausible candidate causal genes for keratoconus postnatal corneal stroma weakening.

Background: Genome-wide association studies (GWAS) for corneal resistance factor (CRF) have identified 100s of loci and proved useful to uncover genetic determinants for keratoconus, a corneal ectasia of early-adulthood onset and common indication of corneal transplantation. In the current absence of studies to probe the impact of candidate causal variants in the cornea, we aimed to fill some of this knowledge gap by leveraging tissue-shared genetic effects. Methods: 181 CRF signals were examined for evidence of colocalization with genetic signals affecting steady-state gene transcription and splicing in adult, non-eye, tissues of the Genotype-Tissue Expression (GTEx) project. Expression of candidate causal genes thus nominated was evaluated in single cell transcriptomes from adult cornea, limbus and conjunctiva. Fine-mapping and colocalization of CRF and keratoconus GWAS signals was also deployed to support their sharing causal variants. Results and discussion: 26.5% of CRF causal signals colocalized with GTEx v8 signals and nominated genes enriched in genes with high and specific expression in corneal stromal cells amongst tissues examined. Enrichment analyses carried out with nearest genes to all 181 CRF GWAS signals indicated that stromal cells of the limbus could be susceptible to signals that did not colocalize with GTEx's. These cells might not be well represented in GTEx and/or the genetic associations might have context specific effects. The causal signals shared with GTEx provide new insights into mediation of CRF genetic effects, including modulation of splicing events. Functionally relevant roles for several implicated genes' products in providing tensile strength, mechano-sensing and signaling make the corresponding genes and regulatory variants prime candidates to be validated and their roles and effects across tissues elucidated. Colocalization of CRF and keratoconus GWAS signals strengthened support for shared causal variants but also highlighted many ways into which likely true shared signals could be missed when using readily available GWAS summary statistics.

Typing myalgic encephalomyelitis by infection at onset: A DecodeME study.

Background: People with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) experience core symptoms of post-exertional malaise, unrefreshing sleep, and cognitive impairment. Despite numbering 0.2-0.4% of the population, no laboratory test is available for their diagnosis, no effective therapy exists for their treatment, and no scientific breakthrough regarding pathogenesis has been made. It remains unknown, despite decades of small-scale studies, whether individuals experience different types of ME/CFS separated by onset-type, sex or age. Methods: DecodeME is a large population-based study of ME/CFS that recruited 17,074 participants in the first 3 months following full launch. Detailed questionnaire responses from UK-based participants who all reported being diagnosed with ME/CFS by a health professional provided an unparalleled opportunity to investigate, using logistic regression, whether ME/CFS severity or onset type is significantly associated with sex, age, illness duration, comorbid conditions or symptoms. Results: The well-established sex-bias among ME/CFS patients is evident in the initial DecodeME cohort: 83.5% of participants were females. What was not known previously was that females tend to have more comorbidities than males. Moreover, being female, being older and being over 10 years from ME/CFS onset are significantly associated with greater severity. Five different ME/CFS onset types were examined in the self-reported data: those with ME/CFS onset (i) after glandular fever (infectious mononucleosis); (ii) after COVID-19 infection; (iii) after other infections; (iv) without an infection at onset; and, (v) where the occurrence of an infection at or preceding onset is not known. Among other findings, ME/CFS onset with unknown infection status was significantly associated with active fibromyalgia. Conclusions: DecodeME participants differ in symptoms, comorbid conditions and/or illness severity when stratified by their sex-at-birth and/or infection around the time of ME/CFS onset.

Myalgic Encephalomyelitis / Chronic Fatigue Syndrome (ME/CFS) is a chronic disease that affects an estimated 250,000 people in the UK. Its defining symptom is post-exertional malaise, an excessive delayed worsening of symptoms following even minor physical or mental exertion. For those with it, ME/CFS means disability and poor quality of life. DecodeME is a research study which is looking for DNA differences between people with ME/CFS and people without any health problems. People with ME/CFS who take part in DecodeME complete a questionnaire that assesses their symptoms and whether they will then be invited to donate a DNA sample. This paper analyses the answers to this questionnaire; we will publish results of the DNA analysis separately. So far, more than 17 thousand people with ME/CFS have completed the DecodeME questionnaire. Their answers help us to address the question: "Are there different types of ME/CFS linked to different causes and how severe it becomes?" Results show that people with ME/CFS do not form a single group reporting similar symptoms and additional medical conditions. Instead, participants who had an infection at the start of their ME/CFS reported a different pattern of symptoms and conditions compared to those without an infection. It is well known that most people with ME/CFS are females. What was not clear previously was that females tend to have more additional health conditions. Also, being female, being older and being over 10 years from ME/CFS onset all make it more likely that someone is more severely affected by their ME/CFS. These findings could indicate that by studying people with different ME/CFS onset-types separately ­ rather than analysing all people with ME/CFS together ­ it will be easier to understand what is going wrong.

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.

The trans-ancestral genomic architecture of glycemic traits.

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Long-term evolution of transposable elements.

Transposable elements are often considered parasitic DNA sequences, able to invade the genome of their host thanks to their self-replicating ability. This colonization process has been extensively studied, both theoretically and experimentally, but their long-term coevolution with the genomes is still poorly understood. In this work, we aim to challenge previous population genetics models by considering features of transposable elements as quantitative, rather than discrete, variables. We also describe more realistic transposable element dynamics by accounting for the variability of the insertion effect, from deleterious to adaptive, as well as mutations leading to a loss of transposition activity and to nonautonomous copies. Individual-based simulations of the behavior of a transposable-element family over several thousand generations show different ways in which active or inactive copies can be maintained for a very long time. Results reveal an unexpected impact of genetic drift on the "junk DNA" content of the genome and strongly question the likelihood of the sustainable long-term stable transposition-selection equilibrium on which numerous previous works were based.

Fine-mapping and cell-specific enrichment at corneal resistance factor loci prioritize candidate causal regulatory variants.

Corneal resistance factor (CRF) is altered during corneal diseases progression. Genome-wide-association studies (GWAS) indicated potential CRF and disease genetics overlap. Here, we characterise 135 CRF loci following GWAS in 76029 UK Biobank participants. Enrichment of extra-cellular matrix gene-sets, genetic correlation with corneal thickness (70% (SE = 5%)), reported keratoconus risk variants at 13 loci, all support relevance to corneal stroma biology. Fine-mapping identifies a subset of 55 highly likely causal variants, 91% of which are non-coding. Genomic features enrichments, using all associated variants, also indicate prominent regulatory causal role. We newly established open chromatin landscapes in two widely-used human cornea immortalised cell lines using ATAC-seq. Variants associated with CRF were significantly enriched in regulatory regions from the corneal stroma-derived cell line and enrichment increases to over 5 fold for variants prioritised by fine-mapping-including at GAS7, SMAD3 and COL6A1 loci. Our analysis generates many hypotheses for future functional validation of aetiological mechanisms.

Genome-Wide Analysis Identifies Two Susceptibility Loci for Positive Thyroid Peroxidase and Thyroglobulin Antibodies.

INTRODUCTION: Thyroid peroxidase (TPO) and thyroglobulin (Tg) are main components of the thyroid gland and play an essential role in thyroid hormone synthesis. The development of antibodies to thyroid peroxidase (TPOAb) and thyroglobulin (TgAb) is the major diagnostic hallmark and early indicator of autoimmune thyroid disease. TPOAb and TgAb are under strong genetic influence; however, genetic factors that determine thyroid antibody positivity are largely unknown. MATERIALS AND METHODS: To identify novel loci associated with TPOAb and/or TgAb positivity, we performed a genome-wide meta-analysis in a total of 2613 individuals from Croatia. Participants with elevated plasma TPOAb and/or TgAb were defined as cases (N = 619) and those with TPOAb and TgAb within reference values were defined as controls (N = 1994). RESULTS: We identified 2 novel loci, of which 1 is located within the YES1 gene (rs77284350, P = 1.50 × 10-8), and the other resides within the IRF8 gene (rs16939945, P = 5.04 × 10-8). CONCLUSIONS: Although the observed variants were associated with TPOAb and TgAb positivity for the first time, both YES1 and IRF8 were previously linked to susceptibility to other autoimmune diseases, and represent plausible biological candidates. This study adds to the knowledge of genetics underlying thyroid antibodies and provides a good basis for further research.

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.

Genomic and drug target evaluation of 90 cardiovascular proteins in 30,931 individuals.

Circulating proteins are vital in human health and disease and are frequently used as biomarkers for clinical decision-making or as targets for pharmacological intervention. Here, we map and replicate protein quantitative trait loci (pQTL) for 90 cardiovascular proteins in over 30,000 individuals, resulting in 451 pQTLs for 85 proteins. For each protein, we further perform pathway mapping to obtain trans-pQTL gene and regulatory designations. We substantiate these regulatory findings with orthogonal evidence for trans-pQTLs using mouse knockdown experiments (ABCA1 and TRIB1) and clinical trial results (chemokine receptors CCR2 and CCR5), with consistent regulation. Finally, we evaluate known drug targets, and suggest new target candidates or repositioning opportunities using Mendelian randomization. This identifies 11 proteins with causal evidence of involvement in human disease that have not previously been targeted, including EGF, IL-16, PAPPA, SPON1, F3, ADM, CASP-8, CHI3L1, CXCL16, GDF15 and MMP-12. Taken together, these findings demonstrate the utility of large-scale mapping of the genetics of the proteome and provide a resource for future precision studies of circulating proteins in human health.