Simulated lunar microgravity transiently arrests growth and induces osteocyte-chondrocyte lineage differentiation in human Wharton's jelly stem cells.

Human Wharton's jelly stem cells (hWJSCs) are multipotent stem cells that are extensively employed in biotechnology applications. However, the impact of simulated lunar microgravity (sµG) on the growth, differentiation, and viability of this cell population is incompletely characterized. We aimed to determine whether acute (72 h) exposure to sµG elicited changes in growth and lineage differentiation in hWJSCs and if putative changes were maintained once exposure to terrestrial gravity (1.0 G) was restored. hWJSCs were cultured under standard 1.0 G conditions prior to being passaged and cultured under sµG (0.16 G) using a random positioning machine. Relative to control, hWJSCs cultured under sµG exhibited marked reductions in growth but not viability. Cell population expression of characteristic stemness markers (CD 73, 90, 105) was significantly reduced under sµG conditions. hWJSCs had 308 significantly upregulated and 328 significantly downregulated genes when compared to 1.0 G culture conditions. Key markers of cell replication, including MKI67, were inhibited. Significant upregulation of osteocyte-chondrocyte lineage markers, including SERPINI1, MSX2, TFPI2, BMP6, COMP, TMEM119, LUM, HGF, CHI3L1 and SPP1, and downregulation of cell fate regulators, including DNMT1 and EZH2, were detected in sµG-exposed hWJSCs. When returned to 1.0 G for 3 days, sµG-exposed hWJSCs had accelerated growth, and expression of stemness markers increased, approaching normal (i.e. 95%) levels. Our data support earlier findings that acute sµG significantly reduces the cell division potential of hWJSCs and suggest that acute sµG-exposure induces reversible changes in cell growth accompanied by osteocyte-chondrocyte changes in lineage differentiation.

Deep transcriptome profiling reveals limited conservation of A-to-I RNA editing in Xenopus.

BACKGROUND: Xenopus has served as a valuable model system for biomedical research over the past decades. Notably, ADAR was first detected in frog oocytes and embryos as an activity that unwinds RNA duplexes. However, the scope of A-to-I RNA editing by the ADAR enzymes in Xenopus remains underexplored. RESULTS: Here, we identify millions of editing events in Xenopus with high accuracy and systematically map the editome across developmental stages, adult organs, and species. We report diverse spatiotemporal patterns of editing with deamination activity highest in early embryogenesis before zygotic genome activation and in the ovary. Strikingly, editing events are poorly conserved across different Xenopus species. Even sites that are detected in both X. laevis and X. tropicalis show largely divergent editing levels or developmental profiles. In protein-coding regions, only a small subset of sites that are found mostly in the brain are well conserved between frogs and mammals. CONCLUSIONS: Collectively, our work provides fresh insights into ADAR activity in vertebrates and suggest that species-specific editing may play a role in each animal's unique physiology or environmental adaptation.

Time-dependent specific molecular signatures of inflammation and remodelling are associated with trimethylamine-N-oxide (TMAO)-induced endothelial cell dysfunction.

Endothelial dysfunction is a critical initiating factor contributing to cardiovascular diseases, involving the gut microbiome-derived metabolite trimethylamine N-oxide (TMAO). This study aims to clarify the time-dependent molecular pathways by which TMAO mediates endothelial dysfunction through transcriptomics and metabolomics analyses in human microvascular endothelial cells (HMEC-1). Cell viability and reactive oxygen species (ROS) generation were also evaluated. TMAO treatment for either 24H or 48H induces reduced cell viability and enhanced oxidative stress. Interestingly, the molecular signatures were distinct between the two time-points. Specifically, few Gene Ontology biological processes (BPs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were modulated after a short (24H) compared to a long (48H) treatment. However, the KEGG signalling pathways namely "tumour necrosis factor (TNF)" and "cytokine-cytokine receptor interaction" were downregulated at 24H but activated at 48H. In addition, at 48H, BPs linked to inflammatory phenotypes were activated (confirming KEGG results), while BPs linked to extracellular matrix (ECM) structural organisation, endothelial cell proliferation, and collagen metabolism were repressed. Lastly, metabolic profiling showed that arachidonic acid, prostaglandins, and palmitic acid were enriched at 48H. This study demonstrates that TMAO induces distinct time-dependent molecular signatures involving inflammation and remodelling pathways, while pathways such as oxidative stress are also modulated, but in a non-time-dependent manner.

Multi-region sampling with paired sample sequencing analyses reveals sub-groups of patients with novel patient-specific dysregulation in Hepatocellular Carcinoma.

BACKGROUND: Conventional differential expression (DE) testing compares the grouped mean value of tumour samples to the grouped mean value of the normal samples, and may miss out dysregulated genes in small subgroup of patients. This is especially so for highly heterogeneous cancer like Hepatocellular Carcinoma (HCC). METHODS: Using multi-region sampled RNA-seq data of 90 patients, we performed patient-specific differential expression testing, together with the patients' matched adjacent normal samples. RESULTS: Comparing the results from conventional DE analysis and patient-specific DE analyses, we show that the conventional DE analysis omits some genes due to high inter-individual variability present in both tumour and normal tissues. Dysregulated genes shared in small subgroup of patients were useful in stratifying patients, and presented differential prognosis. We also showed that the target genes of some of the current targeted agents used in HCC exhibited highly individualistic dysregulation pattern, which may explain the poor response rate. DISCUSSION/CONCLUSION: Our results highlight the importance of identifying patient-specific DE genes, with its potential to provide clinically valuable insights into patient subgroups for applications in precision medicine.

Hypoxia induces HIF1α-dependent epigenetic vulnerability in triple negative breast cancer to confer immune effector dysfunction and resistance to anti-PD-1 immunotherapy.

The hypoxic tumor microenvironment has been implicated in immune escape, but the underlying mechanism remains elusive. Using an in vitro culture system modeling human T cell dysfunction and exhaustion in triple-negative breast cancer (TNBC), we find that hypoxia suppresses immune effector gene expression, including in T and NK cells, resulting in immune effector cell dysfunction and resistance to immunotherapy. We demonstrate that hypoxia-induced factor 1α (HIF1α) interaction with HDAC1 and concurrent PRC2 dependency causes chromatin remolding resulting in epigenetic suppression of effector genes and subsequent immune dysfunction. Targeting HIF1α and the associated epigenetic machinery can reverse the immune effector dysfunction and overcome resistance to PD-1 blockade, as demonstrated both in vitro and in vivo using syngeneic and humanized mice models. These findings identify a HIF1α-mediated epigenetic mechanism in immune dysfunction and provide a potential strategy to overcome immune resistance in TNBC.

IFI16-dependent STING signaling is a crucial regulator of anti-HER2 immune response in HER2+ breast cancer.

Relapse to anti-HER2 monoclonal antibody (mAb) therapies, such as trastuzumab in HER2+ breast cancer (BC), is associated with residual disease progression due to resistance to therapy. Here, we identify interferon-γ inducible protein 16 (IFI16)-dependent STING signaling as a significant determinant of trastuzumab responses in HER2+ BC. We show that down-regulation of immune-regulated genes (IRG) is specifically associated with poor survival of HER2+, but not other BC subtypes. Among IRG, IFI16 is identified as a direct target of EZH2, the underexpression of which leads to deficient STING activation and downstream CXCL10/11 expression in response to trastuzumab treatment. Dual inhibition of EZH2 and histone deacetylase (HDAC) significantly activates IFI16-dependent immune responses to trastuzumab. Notably, a combination of a novel histone methylation inhibitor with an HDAC inhibitor induces complete tumor eradication and long-term T cell memory in a HER2+ BC mouse model. Our findings demonstrate an epigenetic regulatory mechanism suppressing the expression of the IFI16-CXCL10/11 signaling pathway that provides a survival advantage to HER2+ BC to confer resistance to trastuzumab treatment.

Population-based plasma lipidomics reveals developmental changes in metabolism and signatures of obesity risk: a mother-offspring cohort study.

BACKGROUND: Lipids play a vital role in health and disease, but changes to their circulating levels and the link with obesity remain poorly characterized in expecting mothers and their offspring in early childhood. METHODS: LC-MS/MS-based quantitation of 480 lipid species was performed on 2491 plasma samples collected at 4 time points in the mother-offspring Asian cohort GUSTO (Growing Up in Singapore Towards healthy Outcomes). These 4 time points constituted samples collected from mothers at 26-28 weeks of gestation (n=752) and 4-5 years postpartum (n=650), and their offspring at birth (n=751) and 6 years of age (n=338). Linear regression models were used to identify the pregnancy and developmental age-specific variations in the plasma lipidomic profiles, and their association with obesity risk. An independent birth cohort (n=1935), the Barwon Infant Study (BIS), comprising mother-offspring dyads of Caucasian origin was used for validation. RESULTS: Levels of 36% of the profiled lipids were significantly higher (absolute fold change > 1.5 and Padj < 0.05) in antenatal maternal circulation as compared to the postnatal phase, with phosphatidylethanolamine levels changing the most. Compared to antenatal maternal lipids, cord blood showed lower concentrations of most lipid species (79%) except lysophospholipids and acylcarnitines. Changes in lipid concentrations from birth to 6 years of age were much higher in magnitude (log2FC=-2.10 to 6.25) than the changes observed between a 6-year-old child and an adult (postnatal mother) (log2FC=-0.68 to 1.18). Associations of cord blood lipidomic profiles with birth weight displayed distinct trends compared to the lipidomic profiles associated with child BMI at 6 years. Comparison of the results between the child and adult BMI identified similarities in association with consistent trends (R2=0.75). However, large number of lipids were associated with BMI in adults (67%) compared to the children (29%). Pre-pregnancy BMI was specifically associated with decrease in the levels of phospholipids, sphingomyelin, and several triacylglycerol species in pregnancy. CONCLUSIONS: In summary, our study provides a detailed landscape of the in utero lipid environment provided by the gestating mother to the growing fetus, and the magnitude of changes in plasma lipidomic profiles from birth to early childhood. We identified the effects of adiposity on the circulating lipid levels in pregnant and non-pregnant women as well as offspring at birth and at 6 years of age. Additionally, the pediatric vs maternal overlap of the circulating lipid phenotype of obesity risk provides intergenerational insights and early opportunities to track and intervene the onset of metabolic adversities. CLINICAL TRIAL REGISTRATION: This birth cohort is a prospective observational study, which was registered on 1 July 2010 under the identifier NCT01174875 .

Genetic Link Determining the Maternal-Fetal Circulation of Vitamin D.

Vitamin D is an essential micronutrient whose demand is heightened during pregnancy to support the growth of the fetus. Furthermore, the fetus does not produce vitamin D and hence relies exclusively on the supply of maternal vitamin D through the placenta. Vitamin D inadequacy is linked with pregnancy complications and adverse infant outcomes. Hence, early predictive markers of vitamin D inadequacy such as genetic vulnerability are important to both mother and offspring. In this multi-ethnic Asian birth cohort study, we report the first genome-wide association analysis (GWAS) of maternal and fetal vitamin D in circulation. For this, 25-hydroxyvitamin D (25OHD) was measured in the antenatal blood of mothers during mid gestation (n=942), and the cord blood of their offspring at birth (n=812). Around ~7 million single nucleotide polymorphisms (SNPs) were regressed against 25OHD concentrations to identify genetic risk variants. About 41% of mothers had inadequate 25OHD (≤75nmol/L) during pregnancy. Antenatal 25OHD was associated with ethnicity [Malay (Β=-22.32nmol/L, p=2.3×10-26); Indian (Β=-21.85, p=3.1×10-21); reference Chinese], age (Β=0.47/year, p=0.0058), and supplement intake (Β=16.47, p=2.4×10-13). Cord blood 25OHD highly correlated with antenatal vitamin D (r=0.75) and was associated with ethnicity [Malay (Β=-4.44, p=2.2×10-7); Indian (Β=-1.99, p=0.038); reference Chinese]. GWAS analysis identified rs4588, a missense variant in the group-specific component (GC) gene encoding vitamin D binding protein (VDBP), and its defining haplotype, as a risk factor for low antenatal (Β=-8.56/T-allele, p=1.0×10-9) and cord blood vitamin D (Β=-3.22/T-allele, p=1.0×10-8) in all three ethnicities. We also discovered a novel association in a SNP downstream of CYP2J2 (rs10789082), a gene involved in 25-hydroxylation of vitamin D, with vitamin D in pregnant women (Β=-7.68/G-allele, p=1.5×10-8), but not their offspring. As the prevention and early detection of suboptimal vitamin D levels are of profound importance to both mother and offspring's health, the genetic risk variants identified in this study allow risk assessment and precision in early intervention of vitamin D deficiency.

Maternal height, gestational diabetes mellitus and pregnancy complications.

AIMS: To explore the glucose-overload hypothesis of artefactual gestational diabetes (GDM) diagnosis in shorter women during oral glucose tolerance testing (OGTT), by investigating associations between height and maternal glycemia; and GDM and pregnancy complications in height-groups. METHODS: Women from GUSTO (n = 1100, 2009-2010) and NUH (n = 4068, 2017-2018) cohorts underwent a mid-gestation two and three time-point 75 g 2-hour OGTT, respectively. GDM-related complications (hypertensive disorders of pregnancy, preterm delivery, emergency cesarean section, neonatal intensive care unit admission, macrosomia, birthweight) were compared within shorter and taller groups, dichotomized by ethnic-specific median height. RESULTS: Using WHO-1999 criteria, 18.8% (GUSTO) to 22.9% (NUH) of women were diagnosed with GDM-1999; and by WHO-2013 criteria, 21.9% (NUH) had GDM-2013. Each 5-cm height increment was inversely associated with GDM-1999 (adjusted odds ratio [aOR, 95% CI] = 0.81 [0.76-0.87], 2-h glycemia (adjusted ß [aß, 95% CI] = -0.171 mmol/L [-0.208, -0.135]) and 1-h glycemia (aß = -0.160 mmol/L [-0.207, -0.112]). The inverse association between height and 2-h glycemia was most marked in "Other" ethnicities (Eurasians/Caucasians/mixed/other Asians) and Indians, followed by Chinese, then Malays. Compared with non-GDM, GDM-1999 was associated with preterm delivery (aOR = 1.76 [1.19-2.61]) and higher birthweight (aß = 57.16 g [20.95, 93.38]) only among taller but not shorter women. CONCLUSIONS: Only taller women had an increased odds of GDM-related pregnancy complications. An artefactual GDM diagnosis due to glucose-overload among shorter women is plausible.

Molecular docking-aided identification of small molecule inhibitors targeting ß-catenin-TCF4 interaction.

Here we report a molecular docking-based approach to identify small molecules that can target the ß-catenin (ß-cat)-TCF4 protein-protein interaction (PPI), a key effector complex for nuclear Wnt signaling activity. Specifically, we developed and optimized a computational model of ß-cat using publicly available ß-cat protein crystal structures, and existing ß-cat-TCF4 interaction inhibitors as the training set. Using our computational model to an in silico screen predicted 27 compounds as good binders to ß-cat, of which 3 were identified to be effective against a Wnt-responsive luciferase reporter. In vitro functional validation experiments revealed GB1874 as an inhibitor of the Wnt pathway that targets the ß-cat-TCF4 PPI. GB1874 also affected the proliferation and stemness of Wnt-addicted colorectal cancer (CRC) cells in vitro. Encouragingly, GB1874 inhibited the growth of CRC tumor xenografts in vivo, thus demonstrating its potential for further development into therapeutics against Wnt-associated cancer indications.

Oxygen Glucose Deprivation Induced Prosurvival Autophagy Is Insufficient to Rescue Endothelial Function.

Endothelial dysfunction, referring to a disturbance in the vascular homeostasis, has been implicated in many disease conditions including ischemic/reperfusion injury and atherosclerosis. Endothelial mitochondria have been increasingly recognized as a regulator of calcium homeostasis which has implications in the execution of diverse cellular events and energy production. The mitochondrial calcium uniporter complex through which calcium enters the mitochondria is composed of several proteins, including the pore-forming subunit MCU and its regulators MCUR1, MICU1, and MICU2. Mitochondrial calcium overload leads to opening of MPTP (mitochondrial permeability transition pore) and results in apoptotic cell death. Whereas, blockage of calcium entry into the mitochondria results in reduced ATP production thereby activates AMPK-mediated pro-survival autophagy. Here, we investigated the expression of mitochondrial calcium uniporter complex components (MCU, MCUR1, MICU1, and MICU2), induction of autophagy and apoptotic cell death in endothelial cells in response to oxygen-glucose deprivation. Human pulmonary microvascular endothelial cells (HPMVECs) were subjected to oxygen-glucose deprivation (OGD) at 3-h timepoints up to 12 h. Interestingly, except MCUR1 which was significantly downregulated, all other components of the uniporter (MCU, MICU1, and MICU2) remained unchanged. MCUR1 downregulation has been shown to activate AMPK mediated pro-survival autophagy. Similarly, MCUR1 downregulation in response to OGD resulted in AMPK phosphorylation and LC3 processing indicating the activation of pro-survival autophagy. Despite the activation of autophagy, OGD induced Caspase-mediated apoptotic cell death. Blockade of autophagy did not reduce OGD-induced apoptotic cell death whereas serum starvation conferred enough cellular and functional protection. In conclusion, the autophagic flux induced by MCUR1 downregulation in response to OGD is insufficient in protecting endothelial cells from undergoing apoptotic cell death and requires enhancement of autophagic flux by additional means such as serum starvation.

Regulatory sites for splicing in human basal ganglia are enriched for disease-relevant information.

Genome-wide association studies have generated an increasing number of common genetic variants associated with neurological and psychiatric disease risk. An improved understanding of the genetic control of gene expression in human brain is vital considering this is the likely modus operandum for many causal variants. However, human brain sampling complexities limit the explanatory power of brain-related expression quantitative trait loci (eQTL) and allele-specific expression (ASE) signals. We address this, using paired genomic and transcriptomic data from putamen and substantia nigra from 117 human brains, interrogating regulation at different RNA processing stages and uncovering novel transcripts. We identify disease-relevant regulatory loci, find that splicing eQTLs are enriched for regulatory information of neuron-specific genes, that ASEs provide cell-specific regulatory information with evidence for cellular specificity, and that incomplete annotation of the brain transcriptome limits interpretation of risk loci for neuropsychiatric disease. This resource of regulatory data is accessible through our web server, http://braineacv2.inf.um.es/.

Transcriptomic and genetic analyses reveal potential causal drivers for intractable partial epilepsy.

Mesial temporal lobe epilepsy with hippocampal sclerosis represents the most common epilepsy syndrome in adult patients with medically intractable partial epilepsy. Mesial temporal lobe epilepsy is usually regarded as a polygenic and complex disorder, still poorly understood but probably caused and perpetuated by dysregulation of numerous biological networks and cellular functions. The study of gene expression changes by single nucleotide polymorphisms in regulatory elements (expression quantitative trait loci, eQTLs) has been shown to be a powerful complementary approach to the detection and understanding of risk loci by genome-wide association studies. We performed a whole (gene and exon-level) transcriptome analysis on cortical tissue samples (Brodmann areas 20 and 21) from 86 patients with mesial temporal lobe epilepsy with hippocampal sclerosis and 75 neurologically healthy controls. Genome-wide genotyping data from the same individuals (patients and controls) were analysed and paired with the transcriptome data. We report potential epilepsy-risk eQTLs, some of which are specific to tissue from patients with mesial temporal lobe epilepsy with hippocampal sclerosis. We also found large transcriptional and splicing deregulation in mesial temporal lobe epilepsy with hippocampal sclerosis tissue as well as gene networks involving neuronal and glial mechanisms that provide new insights into the cause and maintenance of the seizures. These data (available via the 'Seizubraineac' web-tool resource, www.seizubraineac.org) will facilitate the identification of new therapeutic targets and biomarkers as well as genetic risk variants that could influence epilepsy and pharmacoresistance.

Correction: Synergistic activation of pro-inflammatory type-2 CD8 + T lymphocytes by lipid mediators in severe eosinophilic asthma.

This article was originally published under standard licence, but has now been made available under a [CC BY 4.0] license. The PDF and HTML versions of the paper have been modified accordingly.

Synergistic activation of pro-inflammatory type-2 CD8+ T lymphocytes by lipid mediators in severe eosinophilic asthma.

Human type-2 CD8+ T cells are a cell population with potentially important roles in allergic disease. We investigated this in the context of severe asthma with persistent airway eosinophilia-a phenotype associated with high exacerbation risk and responsiveness to type-2 cytokine-targeted therapies. In two independent cohorts we show that, in contrast to Th2 cells, type-2 cytokine-secreting CD8+CRTH2+ (Tc2) cells are enriched in blood and airways in severe eosinophilic asthma. Concentrations of prostaglandin D2 (PGD2) and cysteinyl leukotriene E4 (LTE4) are also increased in the airways of the same group of patients. In vitro PGD2 and LTE4 function synergistically to trigger Tc2 cell recruitment and activation in a TCR-independent manner. These lipids regulate diverse genes in Tc2 cells inducing type-2 cytokines and many other pro-inflammatory cytokines and chemokines, which could contribute to eosinophilia. These findings are consistent with an important innate-like role for human Tc2 cells in severe eosinophilic asthma and suggest a potential target for therapeutic intervention in this and other diseases.

Development of an objective gene expression panel as an alternative to self-reported symptom scores in human influenza challenge trials.

BACKGROUND: Influenza challenge trials are important for vaccine efficacy testing. Currently, disease severity is determined by self-reported scores to a list of symptoms which can be highly subjective. A more objective measure would allow for improved data analysis. METHODS: Twenty-one volunteers participated in an influenza challenge trial. We calculated the daily sum of scores (DSS) for a list of 16 influenza symptoms. Whole blood collected at baseline and 24, 48, 72 and 96 h post challenge was profiled on Illumina HT12v4 microarrays. Changes in gene expression most strongly correlated with DSS were selected to train a Random Forest model and tested on two independent test sets consisting of 41 individuals profiled on a different microarray platform and 33 volunteers assayed by qRT-PCR. RESULTS: 1456 probes are significantly associated with DSS at 1% false discovery rate. We selected 19 genes with the largest fold change to train a random forest model. We observed good concordance between predicted and actual scores in the first test set (r = 0.57; RMSE = -16.1%) with the greatest agreement achieved on samples collected approximately 72 h post challenge. Therefore, we assayed samples collected at baseline and 72 h post challenge in the second test set by qRT-PCR and observed good concordance (r = 0.81; RMSE = -36.1%). CONCLUSIONS: We developed a 19-gene qRT-PCR panel to predict DSS, validated on two independent datasets. A transcriptomics based panel could provide a more objective measure of symptom scoring in future influenza challenge studies. Trial registration Samples were obtained from a clinical trial with the ClinicalTrials.gov Identifier: NCT02014870, first registered on December 5, 2013.

Novel insights into the genetics of smoking behaviour, lung function, and chronic obstructive pulmonary disease (UK BiLEVE): a genetic association study in UK Biobank.

BACKGROUND: Understanding the genetic basis of airflow obstruction and smoking behaviour is key to determining the pathophysiology of chronic obstructive pulmonary disease (COPD). We used UK Biobank data to study the genetic causes of smoking behaviour and lung health. METHODS: We sampled individuals of European ancestry from UK Biobank, from the middle and extremes of the forced expiratory volume in 1 s (FEV1) distribution among heavy smokers (mean 35 pack-years) and never smokers. We developed a custom array for UK Biobank to provide optimum genome-wide coverage of common and low-frequency variants, dense coverage of genomic regions already implicated in lung health and disease, and to assay rare coding variants relevant to the UK population. We investigated whether there were shared genetic causes between different phenotypes defined by extremes of FEV1. We also looked for novel variants associated with extremes of FEV1 and smoking behaviour and assessed regions of the genome that had already shown evidence for a role in lung health and disease. We set genome-wide significance at p<5 × 10(-8). FINDINGS: UK Biobank participants were recruited from March 15, 2006, to July 7, 2010. Sample selection for the UK BiLEVE study started on Nov 22, 2012, and was completed on Dec 20, 2012. We selected 50,008 unique samples: 10,002 individuals with low FEV1, 10,000 with average FEV1, and 5002 with high FEV1 from each of the heavy smoker and never smoker groups. We noted a substantial sharing of genetic causes of low FEV1 between heavy smokers and never smokers (p=2.29 × 10(-16)) and between individuals with and without doctor-diagnosed asthma (p=6.06 × 10(-11)). We discovered six novel genome-wide significant signals of association with extremes of FEV1, including signals at four novel loci (KANSL1, TSEN54, TET2, and RBM19/TBX5) and independent signals at two previously reported loci (NPNT and HLA-DQB1/HLA-DQA2). These variants also showed association with COPD, including in individuals with no history of smoking. The number of copies of a 150 kb region containing the 5' end of KANSL1, a gene that is important for epigenetic gene regulation, was associated with extremes of FEV1. We also discovered five new genome-wide significant signals for smoking behaviour, including a variant in NCAM1 (chromosome 11) and a variant on chromosome 2 (between TEX41 and PABPC1P2) that has a trans effect on expression of NCAM1 in brain tissue. INTERPRETATION: By sampling from the extremes of the lung function distribution in UK Biobank, we identified novel genetic causes of lung function and smoking behaviour. These results provide new insight into the specific mechanisms underlying airflow obstruction, COPD, and tobacco addiction, and show substantial shared genetic architecture underlying airflow obstruction across individuals, irrespective of smoking behaviour and other airway disease. FUNDING: Medical Research Council.

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.

Farm environment during infancy and lung function at the age of 31: a prospective birth cohort study in Finland.

OBJECTIVES: Farming as an occupation is considered a risk factor for asthma and reduced lung function. By contrast, living on a farm during infancy has been reported to be associated with lower risk of asthma in adulthood. However, little is known about the association between farming environment during infancy and lung function in adulthood. We aimed to study the prospective longitudinal association between farming environment during infancy and lung function in adulthood. DESIGN: A prospective birth cohort study. SETTING: Northern Finland. PARTICIPANTS: 5666 participants born in 1966 were followed up at the age of 31 years. PRIMARY OUTCOME MEASURES: Spirometry at the age of 31 years. RESULTS: To be born into a farmer's family was associated with higher forced expiratory volume in 1 s (FEV1) (36 mL; 95% CI 6 to 67 mL) and forced vital capacity (FVC) (40 mL; 95% CI 5 to 75 mL) at the age of 31 years. Contact with farm animals during infancy was associated with higher FEV1. No associations were seen with FEV1/FVC (FEV1/FVC ratio). Having dogs in childhood revealed similar associations. There was a suggestive dose-dependent association with the number of animal species during childhood and higher FEV1 and FVC at adulthood, especially among women. CONCLUSIONS: Farming environment in early life may have a positive impact on lung function in adulthood.