Chronic Obstructive Pulmonary Disease Exacerbations Increase the Risk of Subsequent Cardiovascular Events: A Longitudinal Analysis of the COPDGene Study.

BACKGROUND: Cardiovascular disease (CVD) is the most important comorbidity in patients with chronic obstructive pulmonary disease (COPD). COPD exacerbations not only contribute to COPD progression but may also elevate the risk of CVD. This study aimed to determine whether COPD exacerbations increase the risk of subsequent CVD events using up to 15 years of prospective longitudinal follow-up data from the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) study. METHODS AND RESULTS: The COPDGene study is a large, multicenter, longitudinal investigation of COPD, including subjects at enrollment aged 45 to 80 years with a minimum of 10 pack-years of smoking history. Cox proportional hazards models and Kaplan-Meier survival curves were used to assess the risk of a composite end point of CVD based on the COPD exacerbation rate. Frequent exacerbators exhibited a higher cumulative incidence of composite CVD end points than infrequent exacerbators, irrespective of the presence of CVD at baseline. After adjusting for covariates, frequent exacerbators still maintained higher hazard ratios (HRs) than the infrequent exacerbator group (without CVD: HR, 1.81 [95% CI, 1.47-2.22]; with CVD: HR, 1.92 [95% CI, 1.51-2.44]). This observation remained consistently significant in moderate to severe COPD subjects and the preserved ratio impaired spirometry population. In the mild COPD population, frequent exacerbators showed a trend toward more CVD events. CONCLUSIONS: COPD exacerbations are associated with an increased risk of subsequent cardiovascular events in subjects with and without preexisting CVD. Patients with COPD experiencing frequent exacerbations may necessitate careful monitoring and additional management for subsequent potential CVD. REGISTRATION: URL: https://www.clinicaltrials.gov; Unique identifier: NCT00608764.

RaggedExperiment: the missing link between genomic ranges and matrices in Bioconductor.

SUMMARY: The RaggedExperiment R / Bioconductor package provides lossless representation of disparate genomic ranges across multiple specimens or cells, in conjunction with efficient and flexible calculations of rectangular-shaped summaries for downstream analysis. Applications include statistical analysis of somatic mutations, copy number, methylation, and open chromatin data. RaggedExperiment is compatible with multimodal data analysis as a component of MultiAssayExperiment data objects, and simplifies data representation and transformation for software developers and analysts. MOTIVATION AND RESULTS: Measurement of copy number, mutation, single nucleotide polymorphism, and other genomic attributes that may be stored as VCF files produce "ragged" genomic ranges data: i.e. across different genomic coordinates in each sample. Ragged data are not rectangular or matrix-like, presenting informatics challenges for downstream statistical analyses. We present the RaggedExperiment R/Bioconductor data structure for lossless representation of ragged genomic data, with associated reshaping tools for flexible and efficient calculation of tabular representations to support a wide range of downstream statistical analyses. We demonstrate its applicability to copy number and somatic mutation data across 33 TCGA cancer datasets.

Maternal vitamin D status modifies the effects of early life tobacco exposure on child lung function.

BACKGROUND: Prior studies suggest that vitamin D may modify the effects of environmental exposures; however, none have investigated gestational vitamin D and cumulative tobacco smoke exposure (TSE) throughout pregnancy and early life. OBJECTIVES: This study investigated the effects of early life TSE on child lung function and the modulatory effects of gestational vitamin D on this association. METHODS: The VDAART (Vitamin D Antenatal Asthma Reduction Trial) recruited nonsmoking pregnant women and followed the mother-child pairs to age 6 years. TSE was assessed with questionnaires and plasma cotinine measurements in the mothers (10-18 and 32-38 gestational weeks) and children (1, 3, and 6 years). Cumulative TSE was calculated from the repeated cotinine measurements. 25-hydroxyvitamin D (25[OH]D) levels were measured at 10-18 and 32-38 gestational weeks. Lung function was assessed at 6 years with spirometry and impulse oscillometry. RESULTS: Of the 476 mother-child pairs, 205 (43%) had increased cotinine levels at ≥1 time point. Cumulative TSE was associated with decreased FEV1 (ß = -0.043 L, P = .018) and increased respiratory resistance at 5 Hz (R5; ß = 0.060 kPa/L/s, P = .002). This association persisted in subjects with insufficient (<30 ng/mL) 25(OH)D levels throughout pregnancy (ß = 0.077 kPa/L/s, P = .016 for R5) but not among those with sufficient levels throughout pregnancy. CONCLUSIONS: Cumulative TSE from pregnancy to childhood is associated with dose- and duration-dependent decreases in child lung function at 6 years even in the absence of reported maternal smoking. Gestational vitamin D may modulate this effect and have therapeutic potential for minimizing the adverse effect of TSE on lung throughout early life. RANDOMIZED TRIAL: Maternal Vitamin D Supplementation to Prevent Childhood Asthma (VDAART); clinicaltrials.gov identifier: NCT00920621.

Higher circulating α-carotene was associated with better cognitive function: an evaluation among the MIND trial participants.

There is emerging evidence linking fruit and vegetable consumption and cognitive function. However, studies focusing on the nutrients underlying this relationship are lacking. We aim to examine the association between plasma nutrients and cognition in a population at risk for cognitive decline with a suboptimal diet. The Mediterranean-DASH Intervention for Neurodegenerative Delay (MIND) trial is a randomized controlled intervention that examines the effects of the MIND diet to prevent cognitive decline. The primary outcome is global cognition. A multivariate linear model was used to investigate the association between blood nutrients and global and/or domain-specific cognition. The model was adjusted for age, sex, education, study site, smoking status, cognitive activities and physical activities. High plasma α-carotene was associated with better global cognition. Participants in the highest tertile of plasma α-carotene had a higher global cognition z score of 0⋅17 when compared with individuals in the lowest tertile (P 0⋅002). Circulating α-carotene levels were also associated with higher semantic memory scores (P for trend 0⋅007). Lutein and zeaxanthin (combined) was positively associated with higher semantic memory scores (P for trend 0⋅009). Our study demonstrated that higher α-carotene levels in blood were associated with higher global cognition scores in a US population at risk for cognitive decline. The higher α-carotene levels in blood reflected greater intakes of fruits, other types of vegetables and lesser intakes of butter and margarine and meat. The higher circulating levels of lutein plus zeaxanthin reflected a dietary pattern with high intakes of fruits, green leafy, other vegetables and cheese, and low consumption of fried foods. Objective nutrient markers in the blood can better characterize dietary intake, which may facilitate the implementation of a tailored dietary intervention for the prevention of cognitive decline.

Low gestational vitamin D level and childhood asthma are related to impaired lung function in high-risk children.

BACKGROUND: Lung function impairment in early life often persists into adulthood. Therefore, identifying risk factors for low childhood lung function is crucial. OBJECTIVE: We examined the effect of 25-hydroxyvitamin D (25[OH]D) level and childhood asthma phenotype on childhood lung function in the Vitamin D Antenatal Asthma Reduction Trial (VDAART). METHODS: The 25(OH)D level was measured at set time points in mothers during pregnancy and in children during early life. On the basis of parental reports, children were categorized into 3 clinical phenotypes: asymptomatic/infrequent wheeze, early transient wheeze, and asthma at age 6 years. Lung function was assessed with impulse oscillometry at ages 4, 5, and 6 years and with spirometry at ages 5 and 6 years. RESULTS: A total of 570 mother-child pairs were included in this post hoc analysis. Mean gestational 25(OH)D-level quartiles were negatively associated with child respiratory resistance at 5 Hz (R5) from age 4 to 6 years (ß, -0.021 kPa/L/s; 95% CI, -0.035 to -0.007; P = .003) and positively associated with FEV1 (ß, 0.018 L; 95% CI, 0.005-0.031; P = .008) and forced vital capacity (ß, 0.022 L; 95% CI, 0.009-0.036; P = .002) from age 5 to 6 years. Children with asthma at age 6 years had lower lung function from age 4 to 6 years than the asymptomatic/infrequent wheeze group (ß, 0.065 kPa/L/s; 95% CI, 0.028 to 0.102; P < .001 for R5 and ß, -0.063 L; 95% CI, -0.099 to -0.028; P < .001 for FEV1). CONCLUSIONS: Low gestational 25(OH)D level and childhood asthma are important risk factors for decreased lung function in early childhood.

Multiomic Integration of Public Oncology Databases in Bioconductor.

PURPOSE: Investigations of the molecular basis for the development, progression, and treatment of cancer increasingly use complementary genomic assays to gather multiomic data, but management and analysis of such data remain complex. The cBioPortal for cancer genomics currently provides multiomic data from > 260 public studies, including The Cancer Genome Atlas (TCGA) data sets, but integration of different data types remains challenging and error prone for computational methods and tools using these resources. Recent advances in data infrastructure within the Bioconductor project enable a novel and powerful approach to creating fully integrated representations of these multiomic, pan-cancer databases. METHODS: We provide a set of R/Bioconductor packages for working with TCGA legacy data and cBioPortal data, with special considerations for loading time; efficient representations in and out of memory; analysis platform; and an integrative framework, such as MultiAssayExperiment. Large methylation data sets are provided through out-of-memory data representation to provide responsive loading times and analysis capabilities on machines with limited memory. RESULTS: We developed the curatedTCGAData and cBioPortalData R/Bioconductor packages to provide integrated multiomic data sets from the TCGA legacy database and the cBioPortal web application programming interface using the MultiAssayExperiment data structure. This suite of tools provides coordination of diverse experimental assays with clinicopathological data with minimal data management burden, as demonstrated through several greatly simplified multiomic and pan-cancer analyses. CONCLUSION: These integrated representations enable analysts and tool developers to apply general statistical and plotting methods to extensive multiomic data through user-friendly commands and documented examples.

Imaging-AMARETTO: An Imaging Genomics Software Tool to Interrogate Multiomics Networks for Relevance to Radiography and Histopathology Imaging Biomarkers of Clinical Outcomes.

PURPOSE: The availability of increasing volumes of multiomics, imaging, and clinical data in complex diseases such as cancer opens opportunities for the formulation and development of computational imaging genomics methods that can link multiomics, imaging, and clinical data. METHODS: Here, we present the Imaging-AMARETTO algorithms and software tools to systematically interrogate regulatory networks derived from multiomics data within and across related patient studies for their relevance to radiography and histopathology imaging features predicting clinical outcomes. RESULTS: To demonstrate its utility, we applied Imaging-AMARETTO to integrate three patient studies of brain tumors, specifically, multiomics with radiography imaging data from The Cancer Genome Atlas (TCGA) glioblastoma multiforme (GBM) and low-grade glioma (LGG) cohorts and transcriptomics with histopathology imaging data from the Ivy Glioblastoma Atlas Project (IvyGAP) GBM cohort. Our results show that Imaging-AMARETTO recapitulates known key drivers of tumor-associated microglia and macrophage mechanisms, mediated by STAT3, AHR, and CCR2, and neurodevelopmental and stemness mechanisms, mediated by OLIG2. Imaging-AMARETTO provides interpretation of their underlying molecular mechanisms in light of imaging biomarkers of clinical outcomes and uncovers novel master drivers, THBS1 and MAP2, that establish relationships across these distinct mechanisms. CONCLUSION: Our network-based imaging genomics tools serve as hypothesis generators that facilitate the interrogation of known and uncovering of novel hypotheses for follow-up with experimental validation studies. We anticipate that our Imaging-AMARETTO imaging genomics tools will be useful to the community of biomedical researchers for applications to similar studies of cancer and other complex diseases with available multiomics, imaging, and clinical data.

Global Alliance for Genomics and Health Meets Bioconductor: Toward Reproducible and Agile Cancer Genomics at Cloud Scale.

PURPOSE: Institutional efforts toward the democratization of cloud-scale data and analysis methods for cancer genomics are proceeding rapidly. As part of this effort, we bridge two major bioinformatic initiatives: the Global Alliance for Genomics and Health (GA4GH) and Bioconductor. METHODS: We describe in detail a use case in pancancer transcriptomics conducted by blending implementations of the GA4GH Workflow Execution Services and Tool Registry Service concepts with the Bioconductor curatedTCGAData and BiocOncoTK packages. RESULTS: We carried out the analysis with a formally archived workflow and container at dockstore.org and a workspace and notebook at app.terra.bio. The analysis identified relationships between microsatellite instability and biomarkers of immune dysregulation at a finer level of granularity than previously reported. Our use of standard approaches to containerization and workflow programming allows this analysis to be replicated and extended. CONCLUSION: Experimental use of dockstore.org and app.terra.bio in concert with Bioconductor enabled novel statistical analysis of large genomic projects without the need for local supercomputing resources but involved challenges related to container design, script archiving, and unit testing. Best practices and cost/benefit metrics for the management and analysis of globally federated genomic data and annotation are evolving. The creation and execution of use cases like the one reported here will be helpful in the development and comparison of approaches to federated data/analysis systems in cancer genomics.

Subtyping COPD by Using Visual and Quantitative CT Imaging Features.

BACKGROUND: Multiple studies have identified COPD subtypes by using visual or quantitative evaluation of CT images. However, there has been no systematic assessment of a combined visual and quantitative CT imaging classification. We integrated visually defined patterns of emphysema with quantitative imaging features and spirometry data to produce a set of 10 nonoverlapping CT imaging subtypes, and we assessed differences between subtypes in demographic features, physiological characteristics, longitudinal disease progression, and mortality. METHODS: We evaluated 9,080 current and former smokers in the COPDGene study who had available volumetric inspiratory and expiratory CT images obtained using a standardized imaging protocol. We defined 10 discrete, nonoverlapping CT imaging subtypes: no CT imaging abnormality, paraseptal emphysema (PSE), bronchial disease, small airway disease, mild emphysema, upper lobe predominant centrilobular emphysema (CLE), lower lobe predominant CLE, diffuse CLE, visual without quantitative emphysema, and quantitative without visual emphysema. Baseline and 5-year longitudinal characteristics and mortality were compared across these CT imaging subtypes. RESULTS: The overall mortality differed significantly between groups (P < .01) and was highest in the 3 moderate to severe CLE groups. Subjects having quantitative but not visual emphysema and subjects with visual but not quantitative emphysema were unique groups with mild COPD, at risk for progression, and with likely different underlying mechanisms. Subjects with PSE and/or moderate to severe CLE had substantial progression of emphysema over 5 years compared with findings in subjects with no CT imaging abnormality (P < .01). CONCLUSIONS: The combination of visual and quantitative CT imaging features reflects different underlying pathological processes in the heterogeneous COPD syndrome and provides a useful approach to reclassify types of COPD. TRIAL REGISTRY: ClinicalTrials.gov; No.: NCT00608764; URL: www.clinicaltrials.gov.

restfulSE: A semantically rich interface for cloud-scale genomics with Bioconductor.

Bioconductor's SummarizedExperiment class unites numerical assay quantifications with sample- and experiment-level metadata.  SummarizedExperiment is the standard Bioconductor class for assays that produce matrix-like data, used by over 200 packages.  We describe the restfulSE package, a deployment of  this data model that supports remote storage.  We illustrate use of SummarizedExperiment with remote HDF5 and Google BigQuery back ends, with two applications in cancer genomics.  Our intent is to allow the use of familiar and semantically meaningful programmatic idioms to query genomic data, while abstracting the remote interface from end users and developers.

Sex differences in gene expression in response to ischemia in the human left ventricular myocardium.

Sex differences exist in the prevalence, presentation and outcomes of ischemic heart disease (IHD). Females have higher risk of heart failure post-myocardial infarction relative to males and are two to three times more likely to die after coronary artery bypass grafting surgery. We examined sex differences in human myocardial gene expression in response to ischemia. Left ventricular biopsies from 68 male/46 female patients undergoing aortic valve replacement surgery were obtained at baseline and after a median 74 min of cold cardioplegic arrest/ischemia. Transcriptomes were quantified by RNA-sequencing. Cell-type enrichment analysis was used to estimate the identity and relative proportions of different cell types in each sample. A sex-specific response to ischemia was observed for 271 genes. Notably, the expression FAM5C, PLA2G4E and CYP1A1 showed an increased expression in females compared to males due to ischemia and DIO3, MT1G and CMA1 showed a decreased expression in females compared to males due to ischemia. Functional annotation analysis revealed sex-specific modulation of the oxytocin signaling pathway and common pathway of fibrin clot formation. Expression quantitative trait locus (eQTL) analysis identified variant-by-sex interaction eQTLs, indicative of sex differences in the genotypic effects on gene expression. Cell-type enrichment analysis showed sex-bias in proportion of specific cell types. Common lymphoid progenitor cells and M2 macrophages were found to increase in female samples from pre- to post-ischemia, but no change was observed in male samples. These differences in response to myocardial ischemia provide insight into the sexual dimorphism of IHD and may aid in the development of sex-specific therapies that reduce myocardial injury.

Alcohol consumption and breast tumor gene expression.

BACKGROUND: Alcohol consumption is an established risk factor for breast cancer and the association generally appears stronger among estrogen receptor (ER)-positive tumors. However, the biological mechanisms underlying this association are not completely understood. METHODS: We analyzed messenger RNA (mRNA) microarray data from both invasive breast tumors (N = 602) and tumor-adjacent normal tissues (N = 508) from participants diagnosed with breast cancer in the Nurses' Health Study (NHS) and NHSII. Multivariable linear regression, controlling for other known breast cancer risk factors, was used to identify differentially expressed genes by pre-diagnostic alcohol intake. For pathway analysis, we performed gene set enrichment analysis (GSEA). Differentially expressed genes or enriched pathway-defined gene sets with false discovery rate (FDR) <0.1 identified in tumors were validated in RNA sequencing data of invasive breast tumors (N = 166) from The Cancer Genome Atlas. RESULTS: No individual genes were significantly differentially expressed by alcohol consumption in the NHS/NHSII. However, GSEA identified 33 and 68 pathway-defined gene sets at FDR <0.1 among 471 ER+ and 127 ER- tumors, respectively, all of which were validated. Among ER+ tumors, consuming 10+ grams of alcohol per day (vs. 0) was associated with upregulation in RNA metabolism and transport, cell cycle regulation, and DNA repair, and downregulation in lipid metabolism. Among ER- tumors, in addition to upregulation in RNA processing and cell cycle, alcohol intake was linked to overexpression of genes involved in cytokine signaling, including interferon and transforming growth factor (TGF)-ß signaling pathways, and translation and post-translational modifications. Lower lipid metabolism was observed in both ER+ tumors and ER+ tumor-adjacent normal samples. Most of the significantly enriched gene sets identified in ER- tumors showed a similar enrichment pattern among ER- tumor-adjacent normal tissues. CONCLUSIONS: Our data suggest that moderate alcohol consumption (i.e. 10+ grams/day, equivalent to one or more drinks/day) is associated with several specific and reproducible biological processes and pathways, which adds potential new insight into alcohol-related breast carcinogenesis.

Expression Quantitative Trait loci (QTL) in tumor adjacent normal breast tissue and breast tumor tissue.

We investigate 71 single nucleotide polymorphisms (SNPs) identified in meta-analytic studies of genome-wide association studies (GWAS) of breast cancer, the majority of which are located in intergenic or intronic regions. To explore regulatory impacts of these variants we conducted expression quantitative loci (eQTL) analyses on tissue samples from 376 invasive postmenopausal breast cancer cases in the Nurses' Health Study (NHS) diagnosed from 1990-2004. Expression analysis was conducted on all formalin-fixed paraffin-embedded (FFPE) tissue samples (and on 264 adjacent normal samples) using the Affymetrix Human Transcriptome Array. Significance and ranking of associations between tumor receptor status and expression variation was preserved between NHS FFPE and TCGA fresh-frozen sample sets (Spearman r = 0.85, p<10^-10 for 17 of the 21 Oncotype DX recurrence signature genes). At an FDR threshold of 10%, we identified 27 trans-eQTLs associated with expression variation in 217 distinct genes. SNP-gene associations can be explored using an open-source interactive browser distributed in a Bioconductor package. Using a new a procedure for testing hypotheses relating SNP content to expression patterns in gene sets, defined as molecular function pathways, we find that loci on 6q14 and 6q25 affect various gene sets and molecular pathways (FDR < 10%). Although the ultimate biological interpretation of the GWAS-identified variants remains to be uncovered, this study validates the utility of expression analysis of this FFPE expression set for more detailed integrative analyses.

Early pregnancy vitamin D status and risk of preeclampsia.

BACKGROUND: Low vitamin D status in pregnancy was proposed as a risk factor of preeclampsia. METHODS: We assessed the effect of vitamin D supplementation (4,400 vs. 400 IU/day), initiated early in pregnancy (10-18 weeks), on the development of preeclampsia. The effects of serum vitamin D (25-hydroxyvitamin D [25OHD]) levels on preeclampsia incidence at trial entry and in the third trimester (32-38 weeks) were studied. We also conducted a nested case-control study of 157 women to investigate peripheral blood vitamin D-associated gene expression profiles at 10 to 18 weeks in 47 participants who developed preeclampsia. RESULTS: Of 881 women randomized, outcome data were available for 816, with 67 (8.2%) developing preeclampsia. There was no significant difference between treatment (N = 408) or control (N = 408) groups in the incidence of preeclampsia (8.08% vs. 8.33%, respectively; relative risk: 0.97; 95% CI, 0.61-1.53). However, in a cohort analysis and after adjustment for confounders, a significant effect of sufficient vitamin D status (25OHD ≥30 ng/ml) was observed in both early and late pregnancy compared with insufficient levels (25OHD <30 ng/ml) (adjusted odds ratio, 0.28; 95% CI, 0.10-0.96). Differential expression of 348 vitamin D-associated genes (158 upregulated) was found in peripheral blood of women who developed preeclampsia (FDR <0.05 in the Vitamin D Antenatal Asthma Reduction Trial [VDAART]; P < 0.05 in a replication cohort). Functional enrichment and network analyses of this vitamin D-associated gene set suggests several highly functional modules related to systematic inflammatory and immune responses, including some nodes with a high degree of connectivity. CONCLUSIONS: Vitamin D supplementation initiated in weeks 10-18 of pregnancy did not reduce preeclampsia incidence in the intention-to-treat paradigm. However, vitamin D levels of 30 ng/ml or higher at trial entry and in late pregnancy were associated with a lower risk of preeclampsia. Differentially expressed vitamin D-associated transcriptomes implicated the emergence of an early pregnancy, distinctive immune response in women who went on to develop preeclampsia. TRIAL REGISTRATION: ClinicalTrials.gov NCT00920621. FUNDING: Quebec Breast Cancer Foundation and Genome Canada Innovation Network. This trial was funded by the National Heart, Lung, and Blood Institute. For details see Acknowledgments.

Lung VITAL: Rationale, design, and baseline characteristics of an ancillary study evaluating the effects of vitamin D and/or marine omega-3 fatty acid supplements on acute exacerbations of chronic respiratory disease, asthma control, pneumonia and lung function in adults.

Laboratory and observational research studies suggest that vitamin D and marine omega-3 fatty acids may reduce risk for pneumonia, acute exacerbations of respiratory diseases including chronic obstructive lung disease (COPD) or asthma, and decline of lung function, but prevention trials with adequate dosing, adequate power, and adequate time to follow-up are lacking. The ongoing Lung VITAL study is taking advantage of a large clinical trial-the VITamin D and OmegA-3 TriaL (VITAL)--to conduct the first major evaluation of the influences of vitamin D and marine omega-3 fatty acid supplementation on pneumonia risk, respiratory exacerbation episodes, asthma control and lung function in adults. VITAL is a 5-year U.S.-wide randomized, double-blind, placebo-controlled, 2 × 2 factorial trial of supplementation with vitamin D3 ([cholecalciferol], 2000 IU/day) and marine omega-3 FA (Omacor® fish oil, eicosapentaenoic acid [EPA]+docosahexaenoic acid [DHA], 1g/day) for primary prevention of CVD and cancer among men and women, at baseline aged ≥50 and ≥55, respectively, with 5107 African Americans. In a subset of 1973 participants from 11 urban U.S. centers, lung function is measured before and two years after randomization. Yearly follow-up questionnaires assess incident pneumonia in the entire randomized population, and exacerbations of respiratory disease, asthma control and dyspnea in a subpopulation of 4314 randomized participants enriched, as shown in presentation of baseline characteristics, for respiratory disease, respiratory symptoms, and history of cigarette smoking. Self-reported pneumonia hospitalization will be confirmed by medical record review, and exacerbations will be confirmed by Center for Medicare and Medicaid Services data review.

Age, Sexual Dimorphism, and Disease Associations in the Developing Human Fetal Lung Transcriptome.

The fetal origins of disease hypothesis suggests that variations in the course of prenatal lung development may affect life-long pulmonary function growth, decline, and pathobiology. Many studies support the existence of differences in the developing lung trajectory in males and females, and sex-specific differences in the prevalence of chronic lung diseases, such as asthma and bronchopulmonary dysplasia. The objectives of this study were to investigate the early developing fetal lung for transcriptomic correlates of postconception age (maturity) and sex, and their associations with chronic lung diseases. We analyzed whole-lung transcriptome profiles of 61 females and 78 males at 54-127 days postconception (dpc) from nonsmoking mothers using unsupervised principal component analysis and supervised linear regression models. We identified dominant transcriptomic correlates for postconception age and sex with corresponding gene sets that were enriched for developing lung structural and functional ontologies. We observed that the transcriptomic sex difference was not a uniform global time shift/lag, rather, lungs of males appear to be more mature than those of females before 96 dpc, and females appear to be more mature than males after 96 dpc. The age correlate gene set was consistently enriched for asthma and bronchopulmonary dysplasia genes, but the sex correlate gene sets were not. Despite sex differences in the developing fetal lung transcriptome, postconception age appears to be more dominant than sex in the effect of early fetal lung developments on disease risk during this early pseudoglandular phase of development.

Public data and open source tools for multi-assay genomic investigation of disease.

Molecular interrogation of a biological sample through DNA sequencing, RNA and microRNA profiling, proteomics and other assays, has the potential to provide a systems level approach to predicting treatment response and disease progression, and to developing precision therapies. Large publicly funded projects have generated extensive and freely available multi-assay data resources; however, bioinformatic and statistical methods for the analysis of such experiments are still nascent. We review multi-assay genomic data resources in the areas of clinical oncology, pharmacogenomics and other perturbation experiments, population genomics and regulatory genomics and other areas, and tools for data acquisition. Finally, we review bioinformatic tools that are explicitly geared toward integrative genomic data visualization and analysis. This review provides starting points for accessing publicly available data and tools to support development of needed integrative methods.

Genetic control of gene expression at novel and established chronic obstructive pulmonary disease loci.

Genetic risk loci have been identified for a wide range of diseases through genome-wide association studies (GWAS), but the relevant functional mechanisms have been identified for only a small proportion of these GWAS-identified loci. By integrating results from the largest current GWAS of chronic obstructive disease (COPD) with expression quantitative trait locus (eQTL) analysis in whole blood and sputum from 121 subjects with COPD from the ECLIPSE Study, this analysis identifies loci that are simultaneously associated with COPD and the expression of nearby genes (COPD eQTLs). After integrative analysis, 19 COPD eQTLs were identified, including all four previously identified genome-wide significant loci near HHIP, FAM13A, and the 15q25 and 19q13 loci. For each COPD eQTL, fine mapping and colocalization analysis to identify causal shared eQTL and GWAS variants identified a subset of sites with moderate-to-strong evidence of harboring at least one shared variant responsible for both the eQTL and GWAS signals. Transcription factor binding site (TFBS) analysis confirms that multiple COPD eQTL lead SNPs disrupt TFBS, and enhancer enrichment analysis for loci with the strongest colocalization signals showed enrichment for blood-related cell types (CD3 and CD4+ T cells, lymphoblastoid cell lines). In summary, integrative eQTL and GWAS analysis confirms that genetic control of gene expression plays a key role in the genetic architecture of COPD and identifies specific blood-related cell types as likely participants in the functional pathway from GWAS-associated variant to disease phenotype.

Smoking-Associated Site-Specific Differential Methylation in Buccal Mucosa in the COPDGene Study.

DNA methylation is a complex, tissue-specific phenomenon that can reflect both endogenous factors and exogenous exposures. Buccal brushings represent an easily accessible source of DNA, which may be an appropriate surrogate tissue in the study of environmental exposures and chronic respiratory diseases. Buccal brushings were obtained from a subset of current and former smokers from the COPDGene study. Genome-wide DNA methylation data were obtained in the discovery cohort (n = 82) using the Illumina HumanMethylation450K array. Empirical Bayes methods were used to test for differential methylation by current smoking status at 468,219 autosomal CpG sites using linear models adjusted for age, sex, and race. Pyrosequencing was performed in a nonoverlapping replication cohort (n = 130). Current smokers were significantly younger than former smokers in both the discovery and replication cohorts. Seven CpG sites were associated with current smoking at a false discovery rate less than 0.05 in the discovery cohort. Six of the seven significant sites were pyrosequenced in the replication cohort; five CpG sites, including sites annotated to CYP1B1 and PARVA, were replicated. Correlations between cumulative smoke exposure and time since smoking cessation were observed in a subset of the significantly associated CpG sites. A significant correlation between reduced lung function and increased radiographic emphysema with methylation at cg02162897 (CYP1B1) was observed among female subjects. Site-specific methylation of DNA isolated from buccal mucosa is associated with exposure to cigarette smoke, and may provide insights into the mechanisms underlying differential susceptibility toward the development of smoking-related chronic respiratory diseases.