Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Primary care electronic health records provide an easily accessible approach ...to investigate TRD clinical and genetic characteristics. MDD defined from primary care records in UK Biobank (UKB) and EXCEED studies was compared with other measures of depression and tested for association with MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least 6 weeks. Clinical-demographic characteristics, SNP-based heritability (h
) and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases. In 230,096 and 8926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD, respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 71-88% with other MDD definitions. In UKB, TRD vs healthy controls and non-TRD vs healthy controls h
was comparable (0.25 SE = 0.04 and 0.19 SE = 0.02, respectively). TRD vs non-TRD was positively associated with the PRS of attention deficit hyperactivity disorder, with lower socio-economic status, obesity, higher neuroticism and other unfavourable clinical characteristics. This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.
Display omitted
•Gene-air pollution interaction effects on lung function are poorly understood.•Seven new potential gene-air pollution interaction signals identified.•Interaction effects identified ...are large enough to be clinically relevant.•Observed interactions include one previously identified lung function signal.•High-risk genetic subgroups potentially more susceptible to outdoor air pollution.
Impaired lung function is predictive of mortality and is a key component of chronic obstructive pulmonary disease. Lung function has a strong genetic component but is also affected by environmental factors such as increased exposure to air pollution, but the effect of their interactions is not well understood.
To identify interactions between genetic variants and air pollution measures which affect COPD risk and lung function. Additionally, to determine whether previously identified lung function genetic association signals showed evidence of interaction with air pollution, considering both individual effects and combined effects using a genetic risk score (GRS).
We conducted a genome-wide gene-air pollution interaction analysis of spirometry measures with three measures of air pollution at home address: particulate matter (PM2.5 & PM10) and nitrogen dioxide (NO2), in approximately 300,000 unrelated European individuals from UK Biobank. We explored air pollution interactions with previously identified lung function signals and determined their combined interaction effect using a GRS.
We identified seven new genome-wide interaction signals (P<5×10-8), and a further ten suggestive interaction signals (P<5×10-7). Additionally, we found statistical evidence of interaction for FEV1/FVC between PM2.5 and previously identified lung function signal, rs10841302, near AEBP2, suggesting increased susceptibility as copies of the G allele increased (but size of the impact was small - interaction beta: -0.363 percentage points, 95% CI: -0.523, -0.203 per 5 µg/m3). There was no observed interaction between air pollutants and the weighted GRS.
We carried out the largest genome-wide gene-air pollution interaction study of lung function and identified potential effects of clinically relevant size and significance. We observed up to 440 ml lower lung function for certain genotypes when exposed to mean levels of outdoor air pollution, which is approximately equivalent to nine years of average normal loss of lung function in adults.
Smoking is a leading risk factor for many of the top ten causes of death worldwide. Of the 1.3 billion smokers globally, 80% live in low- and middle-income countries, where the number of deaths due ...to tobacco use is expected to double in the next decade according to the World Health Organization. Genetic studies have helped to identify biological pathways for smoking behaviours, but have mostly focussed on individuals of European ancestry or living in either North America or Europe. We performed a genome-wide association study of two smoking behaviour traits in 10,558 men of African ancestry living in five African countries and the UK. Eight independent variants were associated with either smoking initiation or cessation at P-value < 5 × 10
, four being monomorphic or rare in European populations. Gene prioritisation strategy highlighted five genes, including SEMA6D, previously described as associated with several smoking behaviour traits. These results confirm the importance of analysing underrepresented populations in genetic epidemiology, and the urgent need for larger genomic studies to boost discovery power to better understand smoking behaviours, as well as many other traits.
Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 ...independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.
Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to ...undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.
The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect ...associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma.
We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping.
Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11).
At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.
Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide; smoking is the main risk factor for COPD, but genetic factors are also relevant contributors. Genome-wide ...association studies (GWAS) of the lung function measures used in the diagnosis of COPD have identified a number of loci, however association signals are often broad and collectively these loci only explain a small proportion of the heritability. In order to examine the association with COPD risk of genetic variants down to low allele frequencies, to aid fine-mapping of association signals and to explain more of the missing heritability, we undertook a targeted sequencing study in 300 COPD cases and 300 smoking controls for 26 loci previously reported to be associated with lung function. We used a pooled sequencing approach, with 12 pools of 25 individuals each, enabling high depth (30x) coverage per sample to be achieved. This pooled design maximised sample size and therefore power, but led to challenges during variant-calling since sequencing error rates and minor allele frequencies for rare variants can be very similar. For this reason we employed a rigorous quality control pipeline for variant detection which included the use of 3 independent calling algorithms. In order to avoid false positive associations we also developed tests to detect variants with potential batch effects and removed them before undertaking association testing. We tested for the effects of single variants and the combined effect of rare variants within a locus. We followed up the top signals with data available (only 67% of collapsing methods signals) in 4,249 COPD cases and 11,916 smoking controls from UK Biobank. We provide suggestive evidence for the combined effect of rare variants on COPD risk in TNXB and in sliding windows within MECOM and upstream of HHIP. These findings can lead to an improved understanding of the molecular pathways involved in the development of COPD.
Homozygosity for the
Z allele causes α
-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, ...and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (β=1.02 cm, p=3.91×10
), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV
) (β=19.36 mL, p=9.21×10
) and FEV
/forced vital capacity (β=0.0031, p=1.22×10
) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.
Non-alcoholic fatty liver disease (NAFLD) is a complex trait with an estimated prevalence of 25% globally. We aimed to identify the genetic variant underlying a four-generation family with ...progressive NAFLD leading to cirrhosis, decompensation, and development of hepatocellular carcinoma in the absence of common risk factors such as obesity and type 2 diabetes.
Exome sequencing and genome comparisons were used to identify the likely causal variant. We extensively characterised the clinical phenotype and post-prandial metabolic responses of family members with the identified novel variant in comparison with healthy non-carriers and wild-type patients with NAFLD. Variant-expressing hepatocyte-like cells (HLCs) were derived from human-induced pluripotent stem cells generated from homozygous donor skin fibroblasts and restored to wild-type using CRISPR-Cas9. The phenotype was assessed using imaging, targeted RNA analysis, and molecular expression arrays.
We identified a rare causal variant c.1691T>C p.I564T (rs745447480) in MTTP, encoding microsomal triglyceride transfer protein (MTP), associated with progressive NAFLD, unrelated to metabolic syndrome and without characteristic features of abetalipoproteinaemia. HLCs derived from a homozygote donor had significantly lower MTP activity and lower lipoprotein ApoB secretion than wild-type cells, while having similar levels of MTP mRNA and protein. Cytoplasmic triglyceride accumulation in HLCs triggered endoplasmic reticulum stress, secretion of pro-inflammatory mediators, and production of reactive oxygen species.
We have identified and characterised a rare causal variant in MTTP, and homozygosity for MTTP p.I564T is associated with progressive NAFLD without any other manifestations of abetalipoproteinaemia. Our findings provide insights into mechanisms driving progressive NAFLD.
A rare genetic variant in the gene MTTP has been identified as responsible for the development of severe non-alcoholic fatty liver disease in a four-generation family with no typical disease risk factors. A cell line culture created harbouring this variant gene was characterised to understand how this genetic variation leads to a defect in liver cells, which results in accumulation of fat and processes that promote disease. This is now a useful model for studying the disease pathways and to discover new ways to treat common types of fatty liver disease.
Display omitted
•A rare MTTP variant, p.I564T, leads to progressive fatty liver disease and cirrhosis.•Hepatocyte-like cells expressing MTP 564T have reduced liver lipoprotein secretion capability.•Reactive oxygen species, NF-κB pathway, and extracellular matrix remodelling drive lipotoxicity.
The UK Biobank (UKB) has recently released genotypes on 152,328 individuals together with extensive phenotypic and lifestyle information. We present a new phasing method, SHAPEIT3, that can handle ...such biobank-scale data sets and results in switch error rates as low as ∼0.3%. The method exhibits O(NlogN) scaling with sample size N, enabling fast and accurate phasing of even larger cohorts.
PDF
