Genetic predisposition, modifiable lifestyles, and their joint effects on human lifespan: evidence from multiple cohort studies.

OBJECTIVE: To investigate the associations across genetic and lifestyle factors with lifespan. DESIGN: A longitudinal cohort study. SETTING: UK Biobank. PARTICIPANTS: 353 742 adults of European ancestry, who were recruited from 2006 to 2010 and were followed up until 2021. EXPOSURES: A polygenic risk score for lifespan with long (highest quintile) risk categories and a weighted healthy lifestyle score, including no current smoking, moderate alcohol consumption, regular physical activity, healthy body shape, adequate sleep duration, and a healthy diet, categorised into favourable, intermediate, and unfavourable lifestyles. MAIN OUTCOME MEASURES: Lifespan defined as the date of death or the censor date minus the date of birth. RESULTS: Of the included 353 742 participants of European ancestry with a median follow-up of 12.86 years, 24 239 death cases were identified. Participants were grouped into three genetically determined lifespan categories including long (20.1%), intermediate (60.1%), and short (19.8%), and into three lifestyle score categories including favourable (23.1%), intermediate (55.6%), and unfavourable (21.3%). The hazard ratio (HR) of death for individuals with a genetic predisposition to a short lifespan was 1.21 (95% CI 1.16 to 1.26) compared to those with a genetic predisposition to a long lifespan. The HR of death for individuals in the unfavourable lifestyle category was 1.78 (95% CI 1.71 to 1.85), compared with those in the favourable lifestyle category. Participants with a genetic predisposition to a short lifespan and an unfavourable lifestyle had 2.04 times (95% CI 1.87 to 2.22) higher rates of death compared with those with a genetic predisposition to a long lifespan and a favourable lifestyle. No multiplicative interaction was detected between the polygenic risk score of lifespan and the weighted healthy lifestyle score (p=0.10). The optimal combination of healthy lifestyles, including never smoking, regular physical activity, adequate sleep duration, and a healthy diet, was derived to decrease risk of premature death (death before 75 years). CONCLUSION: Genetic and lifestyle factors were independently associated with lifespan. Adherence to healthy lifestyles could largely attenuate the genetic risk of a shorter lifespan or premature death. The optimal combination of healthy lifestyles could convey better benefits for a longer lifespan, regardless of genetic background.

Educational attainment, health outcomes and mortality: a within-sibship Mendelian randomization study.

BACKGROUND: Previous Mendelian randomization (MR) studies using population samples (population MR) have provided evidence for beneficial effects of educational attainment on health outcomes in adulthood. However, estimates from these studies may have been susceptible to bias from population stratification, assortative mating and indirect genetic effects due to unadjusted parental genotypes. MR using genetic association estimates derived from within-sibship models (within-sibship MR) can avoid these potential biases because genetic differences between siblings are due to random segregation at meiosis. METHODS: Applying both population and within-sibship MR, we estimated the effects of genetic liability to educational attainment on body mass index (BMI), cigarette smoking, systolic blood pressure (SBP) and all-cause mortality. MR analyses used individual-level data on 72 932 siblings from UK Biobank and the Norwegian HUNT study, and summary-level data from a within-sibship Genome-wide Association Study including >140 000 individuals. RESULTS: Both population and within-sibship MR estimates provided evidence that educational attainment decreased BMI, cigarette smoking and SBP. Genetic variant-outcome associations attenuated in the within-sibship model, but genetic variant-educational attainment associations also attenuated to a similar extent. Thus, within-sibship and population MR estimates were largely consistent. The within-sibship MR estimate of education on mortality was imprecise but consistent with a putative effect. CONCLUSIONS: These results provide evidence of beneficial individual-level effects of education (or liability to education) on adulthood health, independently of potential demographic and family-level confounders.

Unraveling Neuro-Proteogenomic Landscape and Therapeutic Implications for Human Behaviors and Psychiatric Disorders.

Understanding the genetic basis of neuro-related proteins is essential for dissecting the molecular basis of human behavioral traits and the disease etiology of neuropsychiatric disorders. Here, the SCALLOP Consortium conducted a genome-wide association meta-analysis of over 12,500 individuals for 184 neuro-related proteins in human plasma. The analysis identified 117 cis-regulatory protein quantitative trait loci (cis-pQTL) and 166 trans-pQTL. The mapped pQTL capture on average 50% of each protein's heritability. Mendelian randomization analyses revealed multiple proteins showing potential causal effects on neuro-related traits such as sleeping, smoking, feelings, alcohol intake, mental health, and psychiatric disorders. Integrating with established drug information, we validated 13 out of 13 matched combinations of protein targets and diseases or side effects with available drugs, while suggesting hundreds of re-purposing and new therapeutic targets. This consortium effort provides a large-scale proteogenomic resource for biomedical research on human behaviors and other neuro-related phenotypes.

Comparative analysis of transferrin and IgG N-glycosylation in two human populations.

Human plasma transferrin (Tf) N-glycosylation has been mostly studied as a marker for congenital disorders of glycosylation, alcohol abuse, and hepatocellular carcinoma. However, inter-individual variability of Tf N-glycosylation is not known, mainly due to technical limitations of Tf isolation in large-scale studies. Here, we present a highly specific robust high-throughput approach for Tf purification from human blood plasma and detailed characterization of Tf N-glycosylation on the level of released glycans by ultra-high-performance liquid chromatography based on hydrophilic interactions and fluorescence detection (HILIC-UHPLC-FLD), exoglycosidase sequencing, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS). We perform a large-scale comparative study of Tf and immunoglobulin G (IgG) N-glycosylation analysis in two human populations and demonstrate that Tf N-glycosylation is associated with age and sex, along with multiple biochemical and physiological traits. Observed association patterns differ compared to the IgG N-glycome corroborating tissue-specific N-glycosylation and specific N-glycans' role in their distinct physiological functions.

Clinical case study meets population cohort: identification of a BRCA1 pathogenic founder variant in Orcadians.

We multiply ascertained the BRCA1 pathogenic missense variant c.5207T > C; p.Val1736Ala (V1736A) in clinical investigation of breast and ovarian cancer families from Orkney in the Northern Isles of Scotland, UK. We sought to investigate the frequency and clinical relevance of this variant in those of Orcadian ancestry as an exemplar of the value of population cohorts in clinical care, especially in isolated populations. Oral history and birth, marriage and death registrations indicated genealogical linkage of the clinical cases to ancestors from the Isle of Westray, Orkney. Further clinical cases were identified through targeted testing for V1736A in women of Orcadian ancestry attending National Health Service (NHS) genetic clinics for breast and ovarian cancer family risk assessments. The variant segregates with female breast and ovarian cancer in clinically ascertained cases. Separately, exome sequence data from 2088 volunteer participants with three or more Orcadian grandparents, in the ORCADES research cohort, was interrogated to estimate the population prevalence of V1736A in Orcadians. The effects of the variant were assessed using Electronic Health Record (EHR) linkage. Twenty out of 2088 ORCADES research volunteers (~1%) carry V1736A, with a common haplotype around the variant. This allele frequency is ~480-fold higher than in UK Biobank participants. Cost-effectiveness of population screening for BRCA1 founder pathogenic variants has been demonstrated at a carrier frequency below the ~1% observed here. Thus we suggest that Orcadian women should be offered testing for the BRCA1 V1736A founder pathogenic variant, starting with those with known Westray ancestry.

Mendelian randomisation identifies priority groups for prophylactic EBV vaccination.

BACKGROUND: Epstein Barr virus (EBV) infects ~ 95% of the population worldwide and is known to cause adverse health outcomes such as Hodgkin's, non-Hodgkin's lymphomas, and multiple sclerosis. There is substantial interest and investment in developing infection-preventing vaccines for EBV. To effectively deploy such vaccines, it is vital that we understand the risk factors for infection. Why particular individuals do not become infected is currently unknown. The current literature, describes complex, often conflicting webs of intersecting factors-sociodemographic, clinical, genetic, environmental-, rendering causality difficult to decipher. We aimed to use Mendelian randomization (MR) to overcome the issues posed by confounding and reverse causality to determine the causal risk factors for the acquisition of EBV. METHODS: We mapped the complex evidence from the literature prior to this study factors associated with EBV serostatus (as a proxy for infection) into a causal diagram to determine putative risk factors for our study. Using data from the UK Biobank of 8422 individuals genomically deemed to be of white British ancestry between the ages of 40 and 69 at recruitment between the years 2006 and 2010, we performed a genome wide association study (GWAS) of EBV serostatus, followed by a Two Sample MR to determine which putative risk factors were causal. RESULTS: Our GWAS identified two novel loci associated with EBV serostatus. In MR analyses, we confirmed shorter time in education, an increase in number of sexual partners, and a lower age of smoking commencement, to be causal risk factors for EBV serostatus. CONCLUSIONS: Given the current interest and likelihood of a future EBV vaccine, these factors can inform vaccine development and deployment strategies by completing the puzzle of causality. Knowing these risk factors allows identification of those most likely to acquire EBV, giving insight into what age to vaccinate and who to prioritise when a vaccine is introduced.

Cystatin C is associated with adverse COVID-19 outcomes in diverse populations.

COVID-19 has highly variable clinical courses. The search for prognostic host factors for COVID-19 outcome is a priority. We performed logistic regression for ICU admission against a polygenic score (PGS) for Cystatin C (CyC) production in patients with COVID-19. We analyzed the predictive value of longitudinal plasma CyC levels in an independent cohort of patients hospitalized with COVID-19. In four cohorts spanning European and African ancestry populations, we identified a significant association between CyC-production PGS and odds of critical illness (n cases=2,319), with the strongest association captured in the UKB cohort (OR 2.13, 95% CI 1.58-2.87, p=7.12e-7). Plasma proteomics from an independent cohort of hospitalized COVID-19 patients (n cases = 131) demonstrated that CyC production was associated with COVID-specific mortality (p=0.0007). Our findings suggest that CyC may be useful for stratification of patients and it has functional role in the host response to COVID-19.

Investigation of the causal relationships between human IgG N-glycosylation and 12 common diseases associated with changes in the IgG N-glycome.

Changes in the N-glycosylation of immunoglobulin G (IgG) are often observed in pathological states, such as autoimmune, inflammatory, neurodegenerative, cardiovascular diseases and some types of cancer. However, in most cases, it is not clear if the disease onset causes these changes, or if the changes in IgG N-glycosylation are among the risk factors for the diseases. The aim of this study was to investigate the casual relationships between IgG N-glycosylation traits and 12 diseases, in which the alterations of IgG N-glycome were previously reported, using two sample Mendelian randomization (MR) approach. We have performed two sample MR using publicly available summary statistics of genome-wide association studies of IgG N-glycosylation and disease risks. Our results indicate positive causal effect of systemic lupus erythematosus (SLE) on the abundance of N-glycans with bisecting N-acetylglucosamine in the total IgG N-glycome. Therefore, we suggest regarding this IgG glycosylation trait as a biomarker of SLE. We also emphasize the need for more powerful GWAS studies of IgG N-glycosylation to further elucidate the causal effect of IgG N-glycome on the diseases.

Lifestyle and Genetic Factors Modify Parent-of-Origin Effects on the Human Methylome.

BACKGROUND: parent-of-origin effects (POE) play important roles in complex disease and thus understanding their regulation and associated molecular and phenotypic variation are warranted. Previous studies mainly focused on the detection of genomic regions or phenotypes regulated by POE. Understanding whether POE may be modified by environmental or genetic exposures is important for understanding of the source of POE-associated variation, but only a few case studies addressing modifiable POE exist. METHODS: in order to understand this high order of POE regulation, we screened 101 genetic and environmental factors such as 'predicted mRNA expression levels' of DNA methylation/imprinting machinery genes and environmental exposures. POE-mQTL-modifier interaction models were proposed to test the potential of these factors to modify POE at DNA methylation using data from Generation Scotland: The Scottish Family Health Study(N=2315). FINDINGS: a set of vulnerable/modifiable POE-CpGs were identified (modifiable-POE-regulated CpGs, N=3). Four factors, 'lifetime smoking status' and 'predicted mRNA expression levels' of TET2, SIRT1 and KDM1A, were found to significantly modify the POE on the three CpGs in both discovery and replication datasets. We further identified plasma protein and health-related phenotypes associated with the methylation level of one of the identified CpGs. INTERPRETATION: the modifiable POE identified here revealed an important yet indirect path through which genetic background and environmental exposures introduce their effect on DNA methylation, motivating future comprehensive evaluation of the role of these modifiers in complex diseases. FUNDING: NSFC (81971270),H2020-MSCA-ITN(721815), Wellcome (204979/Z/16/Z,104036/Z/14/Z), MRC (MC_UU_00007/10, MC_PC_U127592696), CSO (CZD/16/6,CZB/4/276, CZB/4/710), SFC (HR03006), EUROSPAN (LSHG-CT-2006-018947), BBSRC (BBS/E/D/30002276), SYSU, Arthritis Research UK, NHLBI, NIH.

Genome-wide methylation data improves dissection of the effect of smoking on body mass index.

Variation in obesity-related traits has a genetic basis with heritabilities between 40 and 70%. While the global obesity pandemic is usually associated with environmental changes related to lifestyle and socioeconomic changes, most genetic studies do not include all relevant environmental covariates, so the genetic contribution to variation in obesity-related traits cannot be accurately assessed. Some studies have described interactions between a few individual genes linked to obesity and environmental variables but there is no agreement on their total contribution to differences between individuals. Here we compared self-reported smoking data and a methylation-based proxy to explore the effect of smoking and genome-by-smoking interactions on obesity related traits from a genome-wide perspective to estimate the amount of variance they explain. Our results indicate that exploiting omic measures can improve models for complex traits such as obesity and can be used as a substitute for, or jointly with, environmental records to better understand causes of disease.

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.

Trends in disease incidence and survival and their effect on mortality in Scotland: nationwide cohort study of linked hospital admission and death records 2001-2016.

OBJECTIVES: Identify causes and future trends underpinning Scottish mortality improvements and quantify the relative contributions of disease incidence and survival. DESIGN: Population-based study. SETTING: Linked secondary care and mortality records across Scotland. PARTICIPANTS: 1 967 130 individuals born between 1905 and 1965 and resident in Scotland from 2001 to 2016. MAIN OUTCOME MEASURES: Hospital admission rates and survival within 5 years postadmission for 28 diseases, stratified by sex and socioeconomic status. RESULTS: 'Influenza and pneumonia', 'Symptoms and signs involving circulatory and respiratory systems' and 'Malignant neoplasm of respiratory and intrathoracic organs' were the hospital diagnosis groupings associated with most excess deaths, being both common and linked to high postadmission mortality. Using disease trends, we modelled a mean mortality HR of 0.737 (95% CI 0.730 to 0.745) from one decade of birth to the next, equivalent to a life extension of ~3 years per decade. This improvement was 61% (30%-93%) accounted for by improved disease survival after hospitalisation (principally cancer) with the remainder accounted for by lowered hospitalisation incidence (principally heart disease and cancer). In contrast, deteriorations in infectious disease incidence and survival increased mortality by 9% (~3.3 months per decade). Disease-driven mortality improvements were slightly greater for men than women (due to greater falls in disease incidence), and generally similar across socioeconomic deciles. We project mortality improvements will continue over the next decade but slow by 21% because much progress in disease survival has already been achieved. CONCLUSION: Morbidity improvements broadly explain observed mortality improvements, with progress on prevention and treatment of heart disease and cancer contributing the most. The male-female health gaps are closing, but those between socioeconomic groups are not. Slowing improvements in morbidity may explain recent stalling in improvements of UK period life expectancies. However, these could be offset if we accelerate improvements in the diseases accounting for most deaths and counteract recent deteriorations in infectious disease.

Genomics of 1 million parent lifespans implicates novel pathways and common diseases and distinguishes survival chances.

We use a genome-wide association of 1 million parental lifespans of genotyped subjects and data on mortality risk factors to validate previously unreplicated findings near CDKN2B-AS1, ATXN2/BRAP, FURIN/FES, ZW10, PSORS1C3, and 13q21.31, and identify and replicate novel findings near ABO, ZC3HC1, and IGF2R. We also validate previous findings near 5q33.3/EBF1 and FOXO3, whilst finding contradictory evidence at other loci. Gene set and cell-specific analyses show that expression in foetal brain cells and adult dorsolateral prefrontal cortex is enriched for lifespan variation, as are gene pathways involving lipid proteins and homeostasis, vesicle-mediated transport, and synaptic function. Individual genetic variants that increase dementia, cardiovascular disease, and lung cancer - but not other cancers - explain the most variance. Resulting polygenic scores show a mean lifespan difference of around five years of life across the deciles. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Ageing happens to us all, and as the cabaret singer Maurice Chevalier pointed out, "old age is not that bad when you consider the alternative". Yet, the growing ageing population of most developed countries presents challenges to healthcare systems and government finances. For many older people, long periods of ill health are part of the end of life, and so a better understanding of ageing could offer the opportunity to prolong healthy living into old age. Ageing is complex and takes a long time to study ­ a lifetime in fact. This makes it difficult to discern its causes, among the countless possibilities based on an individual's genes, behaviour or environment. While thousands of regions in an individual's genetic makeup are known to influence their risk of different diseases, those that affect how long they will live have proved harder to disentangle. Timmers et al. sought to pinpoint such regions, and then use this information to predict, based on their DNA, whether someone had a better or worse chance of living longer than average. The DNA of over 500,000 people was read to reveal the specific 'genetic fingerprints' of each participant. Then, after asking each of the participants how long both of their parents had lived, Timmers et al. pinpointed 12 DNA regions that affect lifespan. Five of these regions were new and had not been linked to lifespan before. Across the twelve as a whole several were known to be involved in Alzheimer's disease, smoking-related cancer or heart disease. Looking at the entire genome, Timmers et al. could then predict a lifespan score for each individual, and when they sorted participants into ten groups based on these scores they found that top group lived five years longer than the bottom, on average. Many factors beside genetics influence how long a person will live and our lifespan cannot be read from our DNA alone. Nevertheless, Timmers et al. had hoped to narrow down their search and discover specific genes that directly influence how quickly people age, beyond diseases. If such genes exist, their effects were too small to be detected in this study. The next step will be to expand the study to include more participants, which will hopefully pinpoint further genomic regions and help disentangle the biology of ageing and disease.

Low-frequency variation in TP53 has large effects on head circumference and intracranial volume.

Cranial growth and development is a complex process which affects the closely related traits of head circumference (HC) and intracranial volume (ICV). The underlying genetic influences shaping these traits during the transition from childhood to adulthood are little understood, but might include both age-specific genetic factors and low-frequency genetic variation. Here, we model the developmental genetic architecture of HC, showing this is genetically stable and correlated with genetic determinants of ICV. Investigating up to 46,000 children and adults of European descent, we identify association with final HC and/or final ICV + HC at 9 novel common and low-frequency loci, illustrating that genetic variation from a wide allele frequency spectrum contributes to cranial growth. The largest effects are reported for low-frequency variants within TP53, with 0.5 cm wider heads in increaser-allele carriers versus non-carriers during mid-childhood, suggesting a previously unrecognized role of TP53 transcripts in human cranial development.

GWAS for male-pattern baldness identifies 71 susceptibility loci explaining 38% of the risk.

Male pattern baldness (MPB) or androgenetic alopecia is one of the most common conditions affecting men, reaching a prevalence of ~50% by the age of 50; however, the known genes explain little of the heritability. Here, we present the results of a genome-wide association study including more than 70,000 men, identifying 71 independently replicated loci, of which 30 are novel. These loci explain 38% of the risk, suggesting that MPB is less genetically complex than other complex traits. We show that many of these loci contain genes that are relevant to the pathology and highlight pathways and functions underlying baldness. Finally, despite only showing genome-wide genetic correlation with height, pathway-specific genetic correlations are significant for traits including lifespan and cancer. Our study not only greatly increases the number of MPB loci, illuminating the genetic architecture, but also provides a new approach to disentangling the shared biological pathways underlying complex diseases.

Estrogens regulate glycosylation of IgG in women and men.

The immunologic potency of IgG is modulated by glycosylation, but mechanisms regulating this process are undefined. A role for sex hormones is suggested by differences in IgG glycans between women and men, most prominently with respect to galactose. We therefore assessed IgG galactosylation in 713 healthy adults from 2 cohorts as well as in 159 subjects from 4 randomized controlled studies of endocrine manipulation: postmenopausal women receiving conjugated estrogens, raloxifene, or placebo; premenopausal women deprived of gonadal hormones with leuprolide and treated with estradiol or placebo; men deprived of gonadal hormones with goserelin and given testosterone or placebo; and men deprived of gonadal hormones with goserelin and given testosterone or placebo together with anastrozole to block conversion of testosterone to estradiol. Menopause was associated with an increase in agalactosylated IgG glycans, particularly in the most abundant fucosylated nonbisected (G0F) glycoform. Conjugated estrogens and raloxifene reduced G0F glycans in postmenopausal women, while in premenopausal women leuprolide increased G0F glycans in a manner reversed by estradiol. Among men, goserelin increased G0F glycans, an effect blocked by testosterone through conversion to estradiol. These results establish estrogens as an in vivo modulator of IgG galactosylation in both women and men, defining a pathway by which sex modulates immunity.

Evidence for large-scale gene-by-smoking interaction effects on pulmonary function.

Background: Smoking is the strongest environmental risk factor for reduced pulmonary function. The genetic component of various pulmonary traits has also been demonstrated, and at least 26 loci have been reproducibly associated with either FEV 1 (forced expiratory volume in 1 second) or FEV 1 /FVC (FEV 1 /forced vital capacity). Although the main effects of smoking and genetic loci are well established, the question of potential gene-by-smoking interaction effect remains unanswered. The aim of the present study was to assess, using a genetic risk score approach, whether the effect of these 26 loci on pulmonary function is influenced by smoking. Methods: We evaluated the interaction between smoking exposure, considered as either ever vs never or pack-years, and a 26-single nucleotide polymorphisms (SNPs) genetic risk score in relation to FEV 1 or FEV 1 /FVC in 50 047 participants of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) and SpiroMeta consortia. Results: We identified an interaction ( ßint = -0.036, 95% confidence interval, -0.040 to -0.032, P = 0.00057) between an unweighted 26 SNP genetic risk score and smoking status (ever/never) on the FEV 1 /FVC ratio. In interpreting this interaction, we showed that the genetic risk of falling below the FEV /FVC threshold used to diagnose chronic obstructive pulmonary disease is higher among ever smokers than among never smokers. A replication analysis in two independent datasets, although not statistically significant, showed a similar trend in the interaction effect. Conclusions: This study highlights the benefit of using genetic risk scores for identifying interactions missed when studying individual SNPs and shows, for the first time, that persons with the highest genetic risk for low FEV 1 /FVC may be more susceptible to the deleterious effects of smoking.

Farming, Foreign Holidays, and Vitamin D in Orkney.

Orkney, north of mainland Scotland, has the world's highest prevalence of multiple sclerosis (MS); vitamin D deficiency, a marker of low UV exposure, is also common in Scotland. Strong associations have been identified between vitamin D deficiency and MS, and between UV exposure and MS independent of vitamin D, although causal relationships remain to be confirmed. We aimed to compare plasma 25-hydroxyvitamin D levels in Orkney and mainland Scotland, and establish the determinants of vitamin D status in Orkney. We compared mean vitamin D and prevalence of deficiency in cross-sectional study data from participants in the Orkney Complex Disease Study (ORCADES) and controls in the Scottish Colorectal Cancer Study (SOCCS). We used multivariable regression to identify factors associated with vitamin D levels in Orkney. Mean (standard deviation) vitamin D was significantly higher among ORCADES than SOCCS participants (35.3 (18.0) and 31.7 (21.2), respectively). Prevalence of severe vitamin D deficiency was lower in ORCADES than SOCCS participants (6.6% to 16.2% p = 1.1 x 10(-15)). Older age, farming occupations and foreign holidays were significantly associated with higher vitamin D in Orkney. Although mean vitamin D levels are higher in Orkney than mainland Scotland, this masks variation within the Orkney population which may influence MS risk.

Integrative pathway genomics of lung function and airflow obstruction.

Chronic respiratory disorders are important contributors to the global burden of disease. Genome-wide association studies (GWASs) of lung function measures have identified several trait-associated loci, but explain only a modest portion of the phenotypic variability. We postulated that integrating pathway-based methods with GWASs of pulmonary function and airflow obstruction would identify a broader repertoire of genes and processes influencing these traits. We performed two independent GWASs of lung function and applied gene set enrichment analysis to one of the studies and validated the results using the second GWAS. We identified 131 significantly enriched gene sets associated with lung function and clustered them into larger biological modules involved in diverse processes including development, immunity, cell signaling, proliferation and arachidonic acid. We found that enrichment of gene sets was not driven by GWAS-significant variants or loci, but instead by those with less stringent association P-values. Next, we applied pathway enrichment analysis to a meta-analyzed GWAS of airflow obstruction. We identified several biologic modules that functionally overlapped with those associated with pulmonary function. However, differences were also noted, including enrichment of extracellular matrix (ECM) processes specifically in the airflow obstruction study. Network analysis of the ECM module implicated a candidate gene, matrix metalloproteinase 10 (MMP10), as a putative disease target. We used a knockout mouse model to functionally validate MMP10's role in influencing lung's susceptibility to cigarette smoke-induced emphysema. By integrating pathway analysis with population-based genomics, we unraveled biologic processes underlying pulmonary function traits and identified a candidate gene for obstructive lung disease.

Molecular mechanisms underlying variations in lung function: a systems genetics analysis.

BACKGROUND: Lung function measures reflect the physiological state of the lung, and are essential to the diagnosis of chronic obstructive pulmonary disease (COPD). The SpiroMeta-CHARGE consortium undertook the largest genome-wide association study (GWAS) so far (n=48,201) for forced expiratory volume in 1 s (FEV1) and the ratio of FEV1 to forced vital capacity (FEV1/FVC) in the general population. The lung expression quantitative trait loci (eQTLs) study mapped the genetic architecture of gene expression in lung tissue from 1111 individuals. We used a systems genetics approach to identify single nucleotide polymorphisms (SNPs) associated with lung function that act as eQTLs and change the level of expression of their target genes in lung tissue; termed eSNPs. METHODS: The SpiroMeta-CHARGE GWAS results were integrated with lung eQTLs to map eSNPs and the genes and pathways underlying the associations in lung tissue. For comparison, a similar analysis was done in peripheral blood. The lung mRNA expression levels of the eSNP-regulated genes were tested for associations with lung function measures in 727 individuals. Additional analyses identified the pleiotropic effects of eSNPs from the published GWAS catalogue, and mapped enrichment in regulatory regions from the ENCODE project. Finally, the Connectivity Map database was used to identify potential therapeutics in silico that could reverse the COPD lung tissue gene signature. FINDINGS: SNPs associated with lung function measures were more likely to be eQTLs and vice versa. The integration mapped the specific genes underlying the GWAS signals in lung tissue. The eSNP-regulated genes were enriched for developmental and inflammatory pathways; by comparison, SNPs associated with lung function that were eQTLs in blood, but not in lung, were only involved in inflammatory pathways. Lung function eSNPs were enriched for regulatory elements and were over-represented among genes showing differential expression during fetal lung development. An mRNA gene expression signature for COPD was identified in lung tissue and compared with the Connectivity Map. This in-silico drug repurposing approach suggested several compounds that reverse the COPD gene expression signature, including a nicotine receptor antagonist. These findings represent novel therapeutic pathways for COPD. INTERPRETATION: The system genetics approach identified lung tissue genes driving the variation in lung function and susceptibility to COPD. The identification of these genes and the pathways in which they are enriched is essential to understand the pathophysiology of airway obstruction and to identify novel therapeutic targets and biomarkers for COPD, including drugs that reverse the COPD gene signature in silico. FUNDING: The research reported in this article was not specifically funded by any agency. See Acknowledgments for a full list of funders of the lung eQTL study and the Spiro-Meta CHARGE GWAS.