Search

Associations between mitochondrial genome, OXPHOS pathway, genes and prostate cancer risk overall and by self-reported race/ethnicity

Background: Mitochondria are involved in many processes that are central to the life and death of a cell. Oxidative phosphorylation (OXPHOS), in particular, is known to be altered in carcinogenesis, leading to an increase in the production of reactive oxidative species and glycolysis, one of the hallmarks of cancer cells. Because of this, genetic variation in the mitochondrial genome, which encodes for part of the OXPHOS pathway, has been suggested to play a role in many cancers, including prostate cancer.Methods: We comprehensively examined the role of the mitochondrial genome and prostate cancer risk in 4,086 prostate cancer cases and 3,698 controls from the Multiethnic Cohort (MEC), testing 350 mitochondrial SNPs (mtSNPs) in five racial/ethnic populations—Africans, Asian Americans, Europeans, Latinos, and Native Hawaiians. Logistic regression was conducted to examine single mitochondrial SNP and haplogroup associations. The sequence kernel association test was conducted for gene and pathway analysis.Results: Eleven mtSNPs and haplogroup N were nominally associated with overall prostate cancer risk at P < 0.05. The mitochondrial DNA-encoded OXPHOS pathway, complexes, and genes were not associated with prostate cancer risk. No significant associations were identified after multiple testing corrections (all FDR q > 0.20).Conclusions: The mitochondrial genome was not associated with prostate cancer risk in our study of 7,784 subjects from the MEC.Impact: Our comprehensive study does not support the role of the mitochondrial genome in the risk of prostate cancer. Cancer Epidemiol Biomarkers Prev; 25(6); 1001–3. ©2016 AACR.

Background: In the USA, genetically admixed populations have the highest asthma prevalence and severe asthma exacerbations rates. This could be explained not only by environmental factors but also by genetic variants that exert ethnic-specific effects. However, no admixture mapping has been performed for severe asthma exacerbations. Objective: We sought to identify genetic variants associated with severe asthma exacerbations in Hispanic/Latino subgroups by means of admixture mapping analyses and fine mapping, and to assess their transferability to other populations and potential functional roles. Methods: We performed an admixture mapping in 1124 Puerto Rican and 625 Mexican American children with asthma. Fine-mapping of the significant peaks was performed via allelic testing of common and rare variants. We performed replication across Hispanic/Latino subgroups, and the transferability to non-Hispanic/Latino populations was assessed in 1001 African Americans, 1250 Singaporeans and 941 Europeans with asthma. The effects of the variants on gene expression and DNA methylation from whole blood were also evaluated in participants with asthma and in silico with data obtained through public databases. Results: Genomewide significant associations of Indigenous American ancestry with severe asthma exacerbations were found at 5q32 in Mexican Americans as well as at 13q13-q13.2 and 3p13 in Puerto Ricans. The single nucleotide polymorphism (SNP) rs1144986 (C5orf46) showed consistent effects for severe asthma exacerbations across Hispanic/Latino subgroups, but it was not validated in non-Hispanics/Latinos. This SNP was associated with DPYSL3 DNA methylation and SCGB3A2 gene expression levels. Conclusions: Admixture mapping study of asthma exacerbations revealed a novel locus that exhibited Hispanic/Latino-specific effects and regulated DPYSL3 and SCGB3A2. This work was funded by the Spanish Ministry of Science and Innovation MCIN/AEI/10.13039/501100011033 and the European Regional Development Fund 'ERDF A way of making Europe' by the European Union grant [SAF2017-83417R], and by MCIN/AEI/10.13039/501100011033 [grant PID2020-116274RB-I00]. Whole-genome sequencing for the TOPMed (Trans-Omics in Precision Medicine) program was supported by the NHLBI. Whole-genome sequencing for GALA II (NHLBI TOPMed: Genes-environments and Admixture in Latino Americans) Study (phs000920) was performed at the New York Genome Center [3R01HL117004-01S3]. Centralised read mapping and genotype calling, along with variant quality metrics and filtering, were provided by the TOPMed Informatics Research Center [3R01HL-117626-02S1]. Phenotype harmonization, data management, sample-identity quality control, and general study coordination were provided by the TOPMed Data Coordinating Center [3R01HL-120393-02S1]. WGS of part of GALA II was performed by New York Genome Center under The Centers for Common Disease Genomics of the Genome Sequencing Program (GSP) Grant [UM1 HG008901]. The GSP Coordinating Center [U24 HG008956] contributed to cross-program scientific initiatives and provided logistical and general study coordination. GSP is funded by the National Human Genome Research Institute, the National Heart, Lung, and Blood Institute, and the National Eye Institute [NA]. This work was also supported by the Sandler Family Foundation[NA], the American Asthma Foundation [NA], the RWJF Amos Medical Faculty Development Program [NA], Harry Wm. and Diana V. Hind Distinguished Professor in Pharmaceutical Sciences II [NA], the National Heart, Lung, and Blood Institute of the National Institutes of Health [R01HL117004, R01HL128439, R01HL135156, X01HL134589, R01HL155024, R01AI079139, R01HL141845, R01HL118267 and R01DK113003], the National Institute of Health and Environmental Health Sciences [R01ES015794 and R21ES24844], the National Institute on Minority Health and Health Disparities [R01MD010443 and R56MD013312], the Tobacco-Related Disease Research Program [Award Number 24RT-0025 and 27IR-0030]. The BAMSE study was funded by the Swedish Heart-Lung Foundation and the Swedish Research Council and Region Stockholm (ALF project and database maintenance). The INMA study was funded by the Instituto de Salud Carlos III [Red INMA G03/176, CB06/02/0041], Spanish Ministry of Health [FIS-FEDER PI16/1288, FIS-FEDER PI19/1338; Miguel Servet FEDER 15/0025 and 20/0006], Generalitat Valenciana [BEST/2020/ 059]. This work was also partially funded by Fundación Canaria Instituto de Investigación Sanitaria de Canarias [PIFIISC19/17]. JV was funded by ISCIII and the European Regional Development Fund 'ERDF A way of making Europe' by the European Union [PI16/00049 and PI19/00141]. MP-Y was funded by the Ramón y Cajal Program [RYC-2015-17205] by MCIN/AEI/10.13039/501100011033 and by the European Social Fund 'ESF Investing in your future' and by the EAACI Allergopharma Award 2021. JV and M.P.-Y were funded Instituto de Salud Carlos III (ISCIII) [CB06/06/1088]. EH-L was supported by a fellowship awarded by MCIN/AEI/10.13039/501100011033 and by 'ESF Investing in your future' [PRE2018-083837]. JP-G was supported by a fellowship awarded by the Spanish Ministry of Universities [FPU19/02175].

Introduction DNA methylation studies have associated methylation levels at different CpG sites or genomic regions with lung function. Moreover, genetic ancestry has been associated with lung function in Latinos. However, no epigenome-wide association study (EWAS) of lung function has been performed in this population. Here, we aimed to identify DNA methylation patterns associated with lung function in pediatric asthma among Latinos. Results We conducted an EWAS in whole blood from 250 Puerto Rican and 148 Mexican American children and young adults with asthma. A total of five CpGs exceeded the genome-wide significance threshold of p = 1.17 × 10−7 in the combined analyses from Puerto Ricans and Mexican Americans: cg06035600 (MAP3K6, p = 6.13 × 10−8) showed significant association with pre-bronchodilator Tiffeneau–Pinelli index, the probes cg00914963 (TBC1D16, p = 1.04 × 10−7), cg16405908 (MRGPRE, p = 2.05 × 10−8), and cg07428101 (MUC2, p = 5.02 × 10−9) were associated with post-bronchodilator forced vital capacity (FVC), and cg20515679 (KCNJ6) with post-bronchodilator Tiffeneau–Pinelli index (p = 1.13 × 10−8). However, these markers did not show significant associations in publicly available data from Europeans (p > 0.05). A methylation quantitative trait loci analysis revealed that methylation levels at these CpG sites were regulated by genetic variation in Latinos and the Biobank-based Integrative Omics Studies (BIOS) consortium. Additionally, two differentially methylated regions in REXOC and AURKC were associated with pre-bronchodilator Tiffeneau–Pinelli index (adjusted p Conclusions We replicated previous associations of epigenetic markers with lung function in whole blood and identified novel population-specific associations shared among Latino subgroups.

Telomeres play a key role in chromosomal maintenance and stability. To date, few studies have investigated the association of leukocyte telomere length with risk of cancer incidence and all-cause mortality in a large prospective cohort, particularly of the Asian population. Relative telomere lengths in genomic DNA from peripheral blood samples were quantified using a validated quantitative real-time PCR among 26 540 middle-aged or older Chinese adults. Hazard ratios (HRs) and 95% confidence intervals (CIs) of cancer and deaths by quintiles of telomere length were calculated using the Cox proportional hazards regression method with adjustment for age, sex and other potential confounders. After baseline blood collection, 4353 persons developed cancer and 7609 died. Participants with the longest decile of telomeres had a 26% (95% CI: 11%-44%) higher risk of total cancer incidence compared to the shortest decile after controlling for age, sex and other potential founders (Ptrend < .0001). In contrast, longer telomeres were associated with lower risk of all-cause mortality (HR = 0.93; 95% CI: 0.84-1.03), noncancer death (HR = 0.81; 95% CI: 0.71-0.92), specifically, death from chronic obstructive pulmonary disease and pneumonia (HR = 0.79, 95% CI: 0.70-0.89) and digestive diseases (HR = 0.60, 95% CI: 0.42-0.88). Our findings demonstrated that longer telomeres are associated with increased risk of cancer development overall and several common cancer types including breast, rectal, prostate, pancreatic cancer and lung adenocarcinoma. Our study also confirmed that longer telomeres are associated with a reduced risk of noncancer related death.

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

Although nasopharyngeal samples have been considered the gold standard for COVID-19 testing, variability in viral load across different anatomical sites could cause nasopharyngeal samples to be less sensitive than saliva or nasal samples in certain cases. Self-collected samples have logistical advantages over nasopharyngeal samples, making them amenable to population-scale screening. To evaluate sampling alternatives for population screening, we collected nasopharyngeal, saliva, and nasal samples from two cohorts with varied levels and types of symptoms. In a mixed cohort of 60 symptomatic and asymptomatic participants, we found that saliva had 88% concordance with nasopharyngeal samples when tested in the same testing lab (n = 41) and 68% concordance when tested in different testing labs (n = 19). In a second cohort of 20 participants hospitalized for COVID-19, saliva had 74% concordance with nasopharyngeal samples tested in the same testing lab but detected virus in two participants that tested negative with nasopharyngeal samples on the same day. Medical record review showed that the saliva-based testing sensitivity was related to the timing of symptom onset and disease stage. We find that no sample site will be perfectly sensitive for COVID-19 testing in all situations, and the significance of negative results will always need to be determined in the context of clinical signs and symptoms. Saliva retained high clinical sensitivity for early-stage and presymptomatic COVID-19 while allowing easier collection, minimizing the exposure of health care workers, and need for personal protective equipment and making it a viable option for population-scale testing. IMPORTANCE Methods for COVID-19 detection are necessary for public health efforts to monitor the spread of disease. Nasopharyngeal samples have been considered the best approach for COVID-19 testing. However, alternative samples like self-collected saliva offer advantages for population-scale screening. Meta-analyses of recent studies suggest that saliva is useful for detecting SARS-CoV-2; however, differences in disease prevalence, sample collection, and analysis methods still confound strong conclusions on the utility of saliva compared to nasopharyngeal samples. Here, we find that the sensitivity of saliva testing is related to both the timing of the sample collection relative to symptom onset and the disease stage. Importantly, several clinical vignettes in our cohorts highlight the challenges of medical decision making with limited knowledge of the associations between laboratory test data and the natural biology of infection.

BackgroundSARS-CoV-2 viral entry may disrupt angiotensin II (Ang II) homeostasis in part via ACE2 downregulation, potentially contributing to COVID-19 induced lung injury. Preclinical models of viral pneumonias that utilize ACE2 demonstrate Ang II type 1 receptor (AT1R) blockade mitigates lung injury, though observational COVID-19 data addressing the effect of AT1R blockade remain mixed.MethodsMulticenter, blinded, placebo-controlled randomized trial of losartan (50 mg PO twice daily for 10 days) versus placebo. Hospitalized patients with COVID-19 and a respiratory sequential organ failure assessment score of at least 1 and not already taking a renin-angiotensin-aldosterone system (RAAS) inhibitor were eligible. The primary outcome was the imputed partial pressure of oxygen to fraction of inspired oxygen (PaO2/FiO2) ratio at 7 days. Secondary outcomes included ordinal COVID-19 severity, oxygen, ventilator, and vasopressor-free days, and mortality. Losartan pharmacokinetics (PK) and RAAS components [Ang II, angiotensin-(1–7) (Ang-(1–7)), ACE, ACE2] were measured in a subgroup of participants.FindingsFrom April 2020 - February 2021, 205 participants were randomized, 101 to losartan and 104 to placebo. Compared to placebo, losartan did not significantly affect PaO2/FiO2 ratio at 7 days [difference of -24.8 (95% -55.6 to 6.1; p=0.12)]. Losartan did not improve any secondary clinical outcome, but worsened vasopressor-free days. PK data were consistent with appropriate steady-state concentrations, but we observed no significant effect of losartan on RAAS components.InterpretationInitiation of orally administered losartan to hospitalized patients with COVID-19 and acute lung injury does not improve PaO2 / FiO2 ratio at 7 days. These data may have implications for ongoing clinical trials.Trial RegistrationLosartan for Patients With COVID-19 Requiring Hospitalization (NCT04312009), https://clinicaltrials.gov/ct2/show/NCT04312009

Background The SARS-CoV-2 virus enters cells via Angiotensin-converting enzyme 2 (ACE2), disrupting the renin-angiotensin-aldosterone axis, potentially contributing to lung injury. Treatment with angiotensin receptor blockers (ARBs), such as losartan, may mitigate these effects, though induction of ACE2 could increase viral entry, replication, and worsen disease. Methods This study represents a placebo-controlled blinded randomized clinical trial (RCT) to test the efficacy of losartan on outpatients with COVID-19 across three hospital systems with numerous community sites in Minnesota, U.S. Participants included symptomatic outpatients with COVID-19 not already taking ACE-inhibitors or ARBs, enrolled within 7 days of symptom onset. Patients were randomized to 1:1 losartan (25 mg orally twice daily unless estimated glomerular filtration rate, eGFR, was reduced, when dosing was reduced to once daily) versus placebo for 10 days, and all patients and outcome assesors were blinded. The primary outcome was all-cause hospitalization within 15 days. Secondary outcomes included functional status, dyspnea, temperature, and viral load. (clinicatrials.gov, NCT04311177, closed to new participants) Findings From April to November 2020, 117 participants were randomized 58 to losartan and 59 to placebo, and all were analyzed under intent to treat principles. The primary outcome did not differ significantly between the two arms based on Barnard's test [losartan arm: 3 events (5.2% 95% CI 1.1, 14.4%) versus placebo arm: 1 event (1.7%; 95% CI 0.0, 9.1%)]; proportion difference -3.5% (95% CI -13.2, 4.8%); p = 0.32]. Viral loads were not statistically different between treatment groups at any time point. Adverse events per 10 patient days did not differ signifcantly [0.33 (95% CI 0.22–0.49) for losartan vs. 0.37 (95% CI 0.25–0.55) for placebo]. Due to a lower than expected hospitalization rate and low likelihood of a clinically important treatment effect, the trial was terminated early. Interpretation In this multicenter blinded RCT for outpatients with mild symptomatic COVID-19 disease, losartan did not reduce hospitalizations, though assessment was limited by low event rate. Importantly, viral load was not statistically affected by treatment. This study does not support initiation of losartan for low-risk outpatients. Funding This study was supported by Minnesota Partnership for Biotechnology and Medical Genomics (CON000000076883).

The role of low frequency variants associated with telomere length homeostasis in chronic diseases and mortalities is relatively understudied in the East-Asian population. Here we evaluated low frequency variants, including 1,915,154 Asian specific variants, for leukocyte telomere length (LTL) associations among 25,533 Singapore Chinese samples. Three East Asian specific variants in/near POT1, TERF1 and STN1 genes are associated with LTL (Meta-analysis P 2.49×10−14–6.94×10−10). Rs79314063, a missense variant (p.Asp410His) at POT1, shows effect 5.3 fold higher and independent of a previous common index SNP. TERF1 (rs79617270) and STN1 (rs139620151) are linked to LTL-associated common index SNPs at these loci. Rs79617270 is associated with cancer mortality [HR95%CI = 1.544 (1.173, 2.032), PAdj = 0.018] and 4.76% of the association between the rs79617270 and colon cancer is mediated through LTL. Overall, genetically determined LTL is particularly associated with lung adenocarcinoma [HR95%CI = 1.123 (1.051, 1.201), Padj = 0.007]. Ethnicity-specific low frequency variants may affect LTL homeostasis and associate with certain cancers., Xuling Chang et al. study whether low frequency variants are associated with leukocyte telomere length and whether these variants predispose to incident cancers and mortalities among East Asians. They find that three East Asian specific variants at POT1, TERF1 and STN1 loci are associated with leukocyte telomere length and report a mechanistic pathway linking TERF1, leukocyte telomere length and incident colon cancer.

Genetic factors underlying leukocyte telomere length (LTL) may provide insights into telomere homeostasis, with direct links to disease susceptibility. Genetic evaluation of 23,096 Singaporean Chinese samples identifies 10 genome-wide loci (P, Shortening of leukocyte telomere length (LTL) is associated with age and increased risk for various chronic diseases. Here, the authors report genome-wide association studies for LTL in Singaporean Chinese populations and find that longer LTL associates with less severe outcomes of respiratory disease phenotypes.

Background Age and disease prevalence are the 2 biggest risk factors for Coronavirus disease 2019 (COVID-19) symptom severity and death. We therefore hypothesized that increased biological age, beyond chronological age, may be driving disease-related trends in COVID-19 severity. Methods Using the UK Biobank England data, we tested whether a biological age estimate (PhenoAge) measured more than a decade prior to the COVID-19 pandemic was predictive of 2 COVID-19 severity outcomes (inpatient test positivity and COVID-19-related mortality with inpatient test-confirmed COVID-19). Logistic regression models were used with adjustment for age at the pandemic, sex, ethnicity, baseline assessment centers, and preexisting diseases/conditions. Results Six hundred and thirteen participants tested positive at inpatient settings between March 16 and April 27, 2020, 154 of whom succumbed to COVID-19. PhenoAge was associated with increased risks of inpatient test positivity and COVID-19-related mortality (ORMortality = 1.63 per 5 years, 95% CI: 1.43–1.86, p = 4.7 × 10−13) adjusting for demographics including age at the pandemic. Further adjustment for preexisting diseases/conditions at baseline (ORM = 1.50, 95% CI: 1.30–1.73 per 5 years, p = 3.1 × 10−8) and at the early pandemic (ORM = 1.21, 95% CI: 1.04–1.40 per 5 years, p = .011) decreased the association. Conclusions PhenoAge measured in 2006–2010 was associated with COVID-19 severity outcomes more than 10 years later. These associations were partly accounted for by prevalent chronic diseases proximate to COVID-19 infection. Overall, our results suggest that aging biomarkers, like PhenoAge may capture long-term vulnerability to diseases like COVID-19, even before the accumulation of age-related comorbid conditions.

Background: The SARS-CoV-2 virus enters cells via the ACE2 receptor, disrupting the renin-angiotensin-aldosterone axis, potentially contributing to lung injury. Treatment with angiotensin receptor blockers(ARBs) may mitigate these effects, though induction of ACE2 expression could increase viral entry, replication, and worsen disease. Methods: This study represents a placebo-controlled blinded randomized clinical trial (RCT) to test the efficacy of losartan on outpatients with COVID-19 across three hospital systems in Minnesota, U.S. Participants included symptomatic outpatients with COVID-19 not already taking ACE-inhibitors or ARBs, enrolled within 7 days of symptom onset. Patients were randomized to losartan versus placebo for 10 days. The primary outcome was all-cause hospitalization within 15 days. Secondary outcomes included functional status, dyspnea, temperature, and viral load. Findings: From April-November 2020, 117 participants were randomized 58 to losartan and 59 to placebo. The primary outcome did not differ significantly between the two arms based on Barnard’s test [losartan arm: 3 events (5.2%; 95% CI 1.8-14.1%) versus placebo arm: 1 event (1.7% 95% CI 0.01-9.0%); proportion difference -3.5% (95% CI -12.7-4.5%); p=0.32]. Functional status, self-reported dyspnea, maximum daily temperatures, and viral load were not significantly different between treatment groups. Due to a lower than expected hospitalization rate and low likelihood of a clinically important treatment effect, the trial was terminated early. Interpretation: In this multicenter blinded RCT for outpatients with mild symptomatic COVID-19 disease, losartan did not significantly reduce hospitalizations, affect functional status or assessment of dyspnea. Viral load was not significantly affected by treatment. This study does not support initiation of losartan for low-risk outpatients. Trial Registration: The trial was conducted under the authority of the Food and Drug Administration (IND 148365), was registered on clinicaltrials.gov prior to study initialization (NCT04311177). Funding Statement: This study was supported by Minnesota Partnership for Biotechnology and Medical Genomics (CON000000076883) Declaration of Interests: The authors have no financial conflicts of interest to disclose. MAP, MHN, TWS, and CJT have received additional funding from the Bill and Melinda Gates Foundation to conduct a clinical trial of losartan in inpatients with COVID-19. Ethics Approval Statement: The protocol was approved by a central institutional review board and underwent local context review. The study was conducted following good clinical practice guidelines under the oversight of an independent data safety monitoring board (DSMB).

Widespread high-throughput testing for identification of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection by RT-PCR has been a foundation in the response to the coronavirus disease 2019 (COVID-19) pandemic. Quality assurance metrics for these RT-PCR tests are still evolving as testing is widely implemented. As testing increases, it is important to understand performance characteristics and the errors associated with these tests. Herein, we investigate a high-throughput, laboratory-developed SARS-CoV-2 RT-PCR assay to determine whether modeling can generate quality control metrics that identify false-positive (FP) results due to contamination. This study reviewed repeated clinical samples focusing on positive samples that test negative on re-extraction and PCR, likely representing false positives. To identify and predict false-positive samples, we constructed machine learning–derived models based on the extraction method used. These models identified variables associated with false-positive results across all methods, with sensitivities for predicting FP results ranging between 67% and 100%. Application of the models to all results predicted a total FP rate of 0.08% across all samples, or 2.3% of positive results, similar to reports for other RT-PCR tests for RNA viruses. These models can predict quality control parameters, enabling laboratories to generate decision trees that reduce interpretation errors, allow for automated reflex testing of samples with a high FP probability, improve workflow efficiency, and increase diagnostic accuracy for patient care.

With no known treatments or vaccine, COVID-19 presents a major threat, particularly to older adults, who account for the majority of severe illness and deaths. The age-related susceptibility is partly explained by increased comorbidities including dementia and type II diabetes [1]. While it is unclear why these diseases predispose risk, we hypothesize that increased biological age, rather than chronological age, may be driving disease-related trends in COVID-19 severity with age. To test this hypothesis, we applied our previously validated biological age measure (PhenoAge) [2] composed of chronological age and nine clinical chemistry biomarkers to data of 347,751 participants from a large community cohort in the United Kingdom (UK Biobank), recruited between 2006 and 2010. Other data included disease diagnoses (to 2017), mortality data (to 2020), and the UK national COVID-19 test results (to May 31, 2020) [3]. Accelerated aging 10-14 years prior to the start of the COVID-19 pandemic was associated with test positivity (OR=1.15 per 5-year acceleration, 95% CI: 1.08 to 1.21, p=3.2×10−6) and all-cause mortality with test-confirmed COVID-19 (OR=1.25, per 5-year acceleration, 95% CI: 1.09 to 1.44, p=0.002) after adjustment for demographics including current chronological age and pre-existing diseases or conditions. The corresponding areas under the curves were 0.669 and 0.803, respectively. Biological aging, as captured by PhenoAge, is a better predictor of COVID-19 severity than chronological age, and may inform risk stratification initiatives, while also elucidating possible underlying mechanisms, particularly those related to inflammaging.

BackgroundThe global COVID-19 pandemic has led to an urgent need for scalable methods for clinical diagnostics and viral tracking. Next generation sequencing technologies have enabled large-scale genomic surveillance of SARS-CoV-2 as thousands of isolates are being sequenced around the world and deposited in public data repositories. A number of methods using both short- and long-read technologies are currently being applied for SARS-CoV-2 sequencing, including amplicon approaches, metagenomic methods, and sequence capture or enrichment methods. Given the small genome size, the ability to sequence SARS-CoV-2 at scale is limited by the cost and labor associated with making sequencing libraries.ResultsHere we describe a low-cost, streamlined, all amplicon-based method for sequencing SARS-CoV-2, which bypasses costly and time-consuming library preparation steps. We benchmark this tailed amplicon method against both the ARTIC amplicon protocol and sequence capture approaches and show that an optimized tailed amplicon approach achieves comparable amplicon balance, coverage metrics, and variant calls to the ARTIC v3 approach.ConclusionsThe tailed amplicon method we describe represents a cost-effective and highly scalable method for SARS-CoV-2 sequencing.

Osteosarcoma is an aggressive tumor of the bone that primarily affects young adults and adolescents. Osteosarcoma is characterized by genomic chaos and heterogeneity. While inactivation of tumor suppressor p53TP53is nearly universal other high frequency mutations or structural variations have not been identified. Despite this genomic heterogeneity, key conserved transcriptional programs associated with survival have been identified across human, canine and induced murine osteosarcoma. The epigenomic landscape, including DNA methylation, plays a key role in establishing transcriptional programs in all cell types. The role of epigenetic dysregulation has been studied in a variety of cancers but has yet to be explored at scale in osteosarcoma. Here we examined genome-wide DNA methylation patterns in 24 human and 44 canine osteosarcoma samples identifying groups of highly correlated DNA methylation marks in human and canine osteosarcoma samples. We also link specific DNA methylation patterns to key transcriptional programs in both human and canine osteosarcoma. Building on previous work, we built a DNA methylation-based measure for the presence and abundance of various immune cell types in osteosarcoma. Finally, we determined that the underlying state of the tumor, and not changes in cell composition, were the main driver of differences in DNA methylation across the human and canine samples.SignificanceThis is the first large scale study of DNA methylation in osteosarcoma and lays the ground work for the exploration of DNA methylation programs that help establish conserved transcriptional programs in the context of different genomic landscapes.

Telomeres are crucial in the maintenance of chromosome integrity and genomic stability. Critically short telomeres can trigger programed cell death while cells with longer telomeres may have increased likelihood of replicative errors, resulting in genetic mutations and chromosomal alterations, and ultimately promoting oncogenesis. Data on telomere length and lung cancer risk from large prospective cohort studies are spare. Relative telomere length in peripheral blood leukocytes was quantified using a validated monochrome multiplex quantitative polymerase chain reaction (qPCR) method in 26,540 participants of the Singapore Chinese Health Study. After a follow-up of 12 years, 654 participants developed lung cancer including 288 adenocarcinoma, 113 squamous cell carcinoma and 253 other/unknown histological type. The Cox proportional hazard regression was used to estimate hazard ratio (HR) and 95% confidence interval (CI). HR of lung adenocarcinoma for individuals in the highest comparing the lowest 20 percentile of telomere length was 2.84 (95% CI 1.94-4.14, ptrend

BACKGROUND: Graft-versus-host-disease (GVHD) is a multistep process which involves T cell recognition and priming toward alloantigen, expansion, acquisition of effector function, and repeated tissue injury, resulting in clinical manifestations of the disease. All of these processes have considerable metabolic demands and understanding the key role of mitochondria in cellular metabolism as it relates to GVHD has increased significantly. Mitochondrial DNA (mtDNA) haplotypes have been linked to functional differences in vitro, suggesting they have functional differences at an organismal level. We previously used mtDNA typing to assess the impact of mtDNA haplotypes on outcomes of ~400 allo-HCT patients. This pilot study identified uncommon mtDNA haplotypes potentially associated with inferior outcomes. OBJECTIVE: We sought to validate pilot findings of associations between donor and recipient mitochondrial haplotypes and transplant outcome. STUDY DESIGN: We examined a cohort of 4,143 donor-recipient pairs obtained from the Center for International Blood and Marrow Transplant Research. mtDNA was extracted from whole blood or peripheral blood mononuclear cells from donors and recipients and sequenced to discern haplotype. We used multiple regression analysis to examine the independent association of mtDNA haplotype with overall survival and grade III-IV acute GVHD (aGVHD) adjusting for known risk factors for poor transplant outcome. RESULTS: Neither recipient nor donor mtDNA haplotype reached groupwise significance for overall survival (p-values = 0.26 and 0.39, respectively) or grade III-IV aGVHD (p-values = 0.68 and 0.57, respectively). Adjustment for genomically-determined ancestry in the subset of donor-recipient pairs for which this was available did not materially change results. CONCLUSIONS: We conclude that our original finding was due to chance in a small sample size and that there is essentially no evidence that mtDNA haplotype or haplotype mismatch contribute to risk of serious outcomes following allogeneic transplant.

Proofreading polymerases have 3′ to 5′ exonuclease activity that allows the excision and correction of mis-incorporated bases during DNA replication. In a previous study, we demonstrated that in addition to correcting substitution errors and lowering the error rate of DNA amplification, proofreading polymerases can also edit PCR primers to match template sequences. Primer editing is a feature that can be advantageous in certain experimental contexts, such as amplicon-based microbiome profiling. Here we develop a set of synthetic DNA standards to report on primer editing activity and use these standards to dissect this phenomenon. The primer editing standards allow next-generation sequencing-based enzymological measurements, reveal the extent of editing, and allow the comparison of different polymerases and cycling conditions. We demonstrate that proofreading polymerases edit PCR primers in a concentration-dependent manner, and we examine whether primer editing exhibits any sequence specificity. In addition, we use these standards to show that primer editing is tunable through the incorporation of phosphorothioate linkages. Finally, we demonstrate the ability of primer editing to robustly rescue the drop-out of taxa with 16S rRNA gene-targeting primer mismatches using mock communities and human skin microbiome samples.

A long-standing question of human immunodeficiency virus (HIV) genetic variation and evolution has been whether differences exist in mutation rate and/or mutation spectra among HIV types (i.e., HIV-1 versus HIV-2) and among HIV groups (i.e., HIV-1 groups M–P and HIV-2 groups A–H) and HIV-1 Group M subtypes (i.e., subtypes A–D, F–H, and J–K). To address this, we developed a new single-strand consensus sequencing assay for the determination of HIV mutation frequencies and spectra using the Illumina sequencing platform. This assay enables parallel and standardized comparison of HIV mutagenesis among various viral vectors with lower background error than traditional methods of Illumina library preparation. We found significant differences in viral mutagenesis between HIV types but intriguingly no significant differences among HIV-1 Group M subtypes. More specifically, HIV-1 exhibited higher transition frequencies than HIV-2, due mostly to single G-to-A mutations and (to a lesser extent) G-to-A hypermutation. These data suggest that HIV-2 RT exhibits higher fidelity during viral replication, and taken together, these findings demonstrate that HIV type but not subtype significantly affects viral mutation frequencies and spectra. These differences may inform antiviral and vaccine strategies.

Objectives The purpose of this study was to determine the level of heterogeneity in high grade serous ovarian cancer (HGSOC) by analyzing RNA expression in single epithelial and cancer associated stromal cells. In addition, we explored the possibility of identifying subgroups based on pathway activation and pre-defined signatures from cancer stem cells and chemo-resistant cells.

The zebra mussel,Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources forDreissenaor related bivalves, with nearly 450 million years of divergence between zebra mussels and its closest sequenced relative. Although theD. polymorphagenome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate the highest quality molluscan assembly to date. Through comparative analysis and transcriptomics experiments we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact Steamer-Like Elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found thatD. polymorphahave an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

EgeUn###, Rationale: Several common and rare genetic variants have been associated with idiopathic pulmonary fibrosis, a progressive fibrotic condition that is localized to the lung. Objectives: To develop an integrated understanding of the rare and common variants located in multiple loci that have been reported to contribute to the risk of disease. Methods: We performed deep targeted resequencing (3.69 Mb of DNA) in cases (n = 3,624) and control subjects (n = 4,442) across genes and regions previously associated with disease. We tested for associations between disease and 1) individual common variants via logistic regression and 2) groups of rare variants via sequence kernel association tests. Measurements and Main Results: Statistically significant common variant association signals occurred in all 10 of the regions chosen based on genome-wide association studies. The strongest risk variant is the MUC5B promoter variant rs35705950, with an odds ratio of 5.45 (95% confidence interval, 4.91-6.06) for one copy of the risk allele and 18.68 (95% confidence interval, 13.34-26.17) for two copies of the risk allele (P = 9.60 3 102295). In addition to identifying for the first time that rare variation in FAM13A is associated with disease, we confirmed the role of rare variation in the TERT and RTEL1 gene regions in the risk of IPF, and found that the FAM13A and TERT regions have independent common and rare variant signals. Conclusions: A limited number of common and rare variants contribute to the risk of idiopathic pulmonary fibrosis in each of the resequencing regions, and these genetic variants focus on biological mechanisms of host defense and cell senescence. Copyright © 2019 by the American Thoracic Society, Háskóli Íslands, HI University of Chicago University of Ulsan, UOU Monash University, MU University of Nottingham University of Edinburgh University of Colorado Denver Universitat de Barcelona University of Pittsburgh, Pitt Harvard School of Dental Medicine Massachusetts General Hospital, MGH University of Pennsylvania, Penn University College Cork, UCC City, University of London, City University of California, San Francisco, UCSF School of Medicine, Vanderbilt University Ege Üniversitesi University of Virginia, UV Aarhus Universitetshospital Central Manchester University Hospitals NHS Foundation Trust Columbia University Novartis Sunovion COPD Foundation Pfizer National Heart, Lung, and Blood Institute, NHLBI: R01-HL097163, UH3-HL123442, P01-HL092870 U.S. Department of Defense, DOD: W81XWH-17-1-0597, U01 HL089897, U01 HL089856 AstraZeneca Siemens Foundation North Carolina GlaxoSmithKline Foundation, 1National Jewish Health, Denver, Colorado; 2School of Public Health; 3Department of Medicine, and 54Department of Immunology, University of Colorado Denver, Denver, Colorado; 4Department of Medicine, Vanderbilt University School of Medicine, Nashville, Tennessee; 5Brigham and Women’s Hospital, Harvard School of Medicine, Boston, Massachusetts; 6Oklahoma Medical Research Foundation, Oklahoma City, Oklahoma;7Asan Medical Center, University of Ulsan College of Medicine, Seoul, Korea; 8MRC Centre for Inflammation Research, University of Edinburgh, Edinburgh, United Kingdom; 9Respiratory Medicine Unit, Royal Infirmary of Edinburgh, Edinburgh, United Kingdom; 10Biomedical Research Centre, University of Nottingham, Nottingham, United Kingdom; 11Royal Brompton Hospital and Imperial College, London, United Kingdom; 12Simmons Center for Interstitial Lung Disease, University of Pittsburgh, Pittsburgh, Pennsylvania; 13Department of Medicine, University of California, San Francisco, San Francisco, California; 14Gilead Sciences, Foster City, California; 15Department of Medicine, University of Chicago, Chicago, Illinois; 16Department of Medicine, University of Virginia, Charlottesville, Virginia; 17Perelman School of Medicine at the University of Pennsylvania, Philadelphia, Pennsylvania; 18National Hospital Organization Kinki-Chuo Chest Medical Center, Osaka, Japan; 19National Hospital Organization Tokyo National Hospital, Tokyo, Japan; 20Helmholtz Zentrum München, Neuherberg, Germany; 21University Hospital Munich, Munich, Germany; 22Department of Pulmonology, Ege University Hospital, Bornova, Izmir, Turkey; 23Royal Prince Alfred Hospital and University of Sydney, Sydney, Australia; 24Alfred Hospital and Monash University, Melbourne, Australia; 25Pulmonary Medicine, GB Morgagni Hospital, Forlì, Italy; 26Department of Diseases of the Thorax, Ospedale GB Morgagni, Forlì, Italy; 27Université Paris Diderot and Hôpital Bichat, Paris, France; 28Royal Papworth Hospital and Cambridge University Hospitals NHS Foundation Trust, Cambridge, United Kingdom; 29Respiratory Department, University Hospital of Bellvitge, University of Barcelona, Barcelona, Spain; 30National University Hospital of Iceland, University of Iceland, Reykjavik, Iceland; 31Department of Medicine, Columbia University Irving Medical Center, New York, New York; 32Division of Pulmonary and Critical Care Medicine, Massachusetts General Hospital, Boston, Massachusetts; 33Department of Respiratory Diseases and Allergy, Aarhus University Hospital, Aarhus, Denmark; 34Instituto Nacional de Enfermedades Respiratorias “Ismael Cosio Villegas,” México City, México; 35Universidad Nacional Autónoma de México, México City, México; 36Cork University Hospital and University College Cork, Cork, Ireland; 37St. Vincent’s University Hospital, Dublin, Ireland; 38School of Medicine, University College Dublin, Dublin, Ireland; 39Department of Respiratory Medicine, First Faculty of Medicine Charles University and Thomayer Hospital, Prague, Czech Republic; 40University of British Columbia, Vancouver, Canada; 41Tokyo Medical and Dental University, Tokyo, Japan; 42Institute for Respiratory Health and; 43Centre for Cell Therapy and Regenerative Medicine, School of Biomedical Sciences, The University of Western Australia, Perth, Australia; 44Department of Medicine, Warren Alpert Medical School of Brown University, Providence, Rhode Island; 45Advanced Lung Disease and Transplant Program, Inova Fairfax Hospital, Falls Church, Virginia; 46Department of Pulmonary Medicine, Gazi University School of Medicine, Ankara, Turkey; 47Department of Medicine, University of Alabama at Birmingham, Birmingham, Alabama; 48Ruhrlandklinik, University Hospital, University of Duisburg-Essen, Essen, Germany; 49Department of Medicine, Tallaght University Hospital, Trinity College Dublin, Dublin, Ireland; 50CardioPulmonary Reserach Center, Alliance Pulmonary Group, Guaynabo, Puerto Rico; 51Department of Internal Medicine, University of Genoa, Genoa, Italy; 52Biomedical Genomics Center, University of Minnesota; Minneapolis, Minnesota; and 53Northwest Genomics Center, University of Washington, Seattle, Washington -- This research was supported by grants from the NHLBI (R01-HL097163, P01-HL092870, and UH3-HL123442) and the Department of Defense (W81XWH-17-1-0597). The COPDGene project was supported by Award Number U01 HL089897 and Award Number U01 HL089856 from the NHLBI. The COPDGene project was also supported by the COPD Foundation through contributions made to an industry advisory board comprised of AstraZeneca, Boehringer Ingelheim, GlaxoSmithKline, Novartis, Pfizer, Siemens, and Sunovion.

Clinical laboratories have adopted next generation sequencing (NGS) as a gold standard for the diagnosis of hereditary disorders because of its analytic accuracy, high throughput, and potential for cost-effectiveness. We describe the implementation of a single broad-based NGS sequencing assay to meet the genetic testing needs at the University of Minnesota. A single hybrid capture library preparation was used for each test ordered, data was informatically blinded to clinically-ordered genes, and identified variants were reviewed and classified by genetic counselors and molecular pathologists. We performed 2509 sequencing tests from August 2012 till December 2017. The diagnostic yield has remained steady at 25%, but the number of variants of uncertain significance (VUS) included in a patient report decreased over time with 50% of the patient reports including at least one VUS in 2012 and only 22% of the patient reports reporting a VUS in 2017 (p = .002). Among the various clinical specialties, the diagnostic yield was highest in dermatology (60% diagnostic yield) and ophthalmology (42% diagnostic yield) while the diagnostic yield was lowest in gastrointestinal diseases and pulmonary diseases (10% detection yield in both specialties). Deletion/duplication analysis was also implemented in a subset of panels ordered, with 9% of samples having a diagnostic finding using the deletion/duplication analysis. We have demonstrated the feasibility of this broad-based NGS platform to meet the needs of our academic institution by aggregating a sufficient sample volume from many individually rare tests and providing a flexible ordering for custom, patient-specific panels. Keywords: Next generation sequencing, Molecular diagnostics, Panel testing, Diagnostic yield, Variants of uncertain significance, Copy number variation

Quantification of DNA sequence tags from engineered constructs such as plasmids, transposons, or other transgenes underlies many functional genomics measurements. Typically, such measurements rely on PCR followed by next-generation sequencing. However, PCR amplification can introduce significant quantitative error. We describe REcount, a novel PCR-free direct counting method. Comparing measurements of defined plasmid pools to droplet digital PCR data demonstrates that REcount is highly accurate and reproducible. We use REcount to provide new insights into clustering biases due to molecule length across different Illumina sequencers and illustrate the impacts on interpretation of next-generation sequencing data and the economics of data generation. Electronic supplementary material The online version of this article (10.1186/s13059-019-1691-6) contains supplementary material, which is available to authorized users.

BackgroundThe mitochondrial genome encodes for thirty-seven proteins, among them thirteen are essential for the oxidative phosphorylation (OXPHOS) system. Inherited variation in mitochondrial genes may influence cancer development through changes in mitochondrial proteins, altering the OXPHOS process and promoting the production of reactive oxidative species.MethodsTo investigate the association between mitochondrial genetic variation and breast cancer risk, we tested 314 mitochondrial SNPs (mtSNPs), capturing four complexes of the mitochondrial OXPHOS pathway and mtSNP groupings for rRNA and tRNA, in 2,723 breast cancer cases and 3,260 controls from the Multiethnic Cohort Study.ResultsWe examined the collective set of 314 mtSNPs as well as subsets of mtSNPs grouped by mitochondrial OXPHOS pathway, complexes, and genes, using the sequence kernel association test and adjusting for age, sex, and principal components of global ancestry. We also tested haplogroup associations using unconditional logistic regression and adjusting for the same covariates. Stratified analyses were conducted by self-reported maternal race/ethnicity. No significant mitochondrial OXPHOS pathway, gene, and haplogroup associations were observed in African Americans, Asian Americans, Latinos, and Native Hawaiians. In European Americans, a global test of all genetic variants of the mitochondrial genome identified an association with breast cancer risk (P = 0.017, q = 0.102). In mtSNP-subset analysis, the gene MT-CO2 (P = 0.001, q = 0.09) in Complex IV (cytochrome c oxidase) and MT-ND2 (P = 0.004, q = 0.19) in Complex I (NADH dehydrogenase (ubiquinone)) were significantly associated with breast cancer risk.ConclusionsIn summary, our findings suggest that collective mitochondrial genetic variation and particularly in the MT-CO2 and MT-ND2 may play a role in breast cancer risk among European Americans. Further replication is warranted in larger populations and future studies should evaluate the contribution of mitochondrial proteins encoded by both the nuclear and mitochondrial genomes to breast cancer risk.

Quantification of DNA sequence tags associated with engineered genetic constructs underlies many genomics measurements. Typically, such measurements are done using PCR to enrich sequence tags and add adapters, followed by next-generation sequencing (NGS). However, PCR amplification can introduce significant quantitative error into these measurements. Here we describe REcount, a novel PCR-free direct counting method for NGS-based quantification of engineered genetic constructs. By comparing measurements of defined plasmid pools to droplet digital PCR data, we demonstrate that this method is highly accurate and reproducible. We further demonstrate that the REcount approach is amenable to multiplexing through the use of orthogonal restriction enzymes. Finally, we use REcount to provide new insights into clustering biases due to molecule length across different Illumina sequencing platforms.

The purpose of this study was to determine the level of heterogeneity in high grade serous ovarian cancer (HGSOC) by analyzing RNA expression in single epithelial and cancer associated stromal cells. In addition, we explored the possibility of identifying subgroups based on pathway activation and pre-defined signatures from cancer stem cells and chemo-resistant cells.A fresh, HGSOC tumor specimen derived from ovary was enzymatically digested and depleted of immune infiltrating cells. RNA sequencing was performed on 92 single cells and 66 of these single cell datasets passed quality control checks. Sequences were analyzed using multiple bioinformatics tools, including clustering, principle components analysis, and geneset enrichment analysis to identify subgroups and activated pathways. Immunohistochemistry for ovarian cancer, stem cell and stromal markers was performed on adjacent tumor sections.Analysis of the gene expression patterns identified two major subsets of cells characterized by epithelial and stromal gene expression patterns. The epithelial group was characterized by proliferative genes including genes associated with oxidative phosphorylation and MYC activity, while the stromal group was characterized by increased expression of extracellular matrix (ECM) genes and genes associated with epithelial-to-mesenchymal transition (EMT). Neither group expressed a signature correlating with published chemo-resistant gene signatures, but many cells, predominantly in the stromal subgroup, expressed markers associated with cancer stem cells.Single cell sequencing provides a means of identifying subpopulations of cancer cells within a single patient. Single cell sequence analysis may prove to be critical for understanding the etiology, progression and drug resistance in ovarian cancer.

Background: Massively parallel, or next-generation, sequencing is a powerful technique for the assessment of somatic genomic alterations in cancer samples. Numerous gene targets can be interrogated simultaneously with a high degree of sensitivity. The clinical standard of care for many advanced solid and hematologic malignancies currently requires mutation analysis of several genes in the front-line setting, making focused next generation sequencing (NGS) assays an effective tool for clinical molecular diagnostic laboratories. Methods: We have utilized an integrated microfluidics circuit (IFC) technology for multiplex PCR-based library preparation coupled with a bioinformatic method designed to enhance indel detection. A parallel low input PCR-based library preparation method was developed for challenging specimens with low DNA yield. Computational data filters were written to optimize analytic sensitivity and specificity for clinically relevant variants. Results: Minimum sequencing coverage and precision of variant calls were the two primary criteria used to establish minimum DNA mass input onto the IFC. Wet-bench and bioinformatics protocols were modified based on data from the optimization and familiarization process to improve assay performance. The NGS platform was then clinically validated for single nucleotide and indel (up to 93 base pair) variant detection with overall analytic accuracy of 98% (97% sensitivity; 100% specificity) using as little as 3 ng of formalin-fixed, paraffin-embedded DNA or 0.3 ng of unfixed DNA. Conclusions: We created a targeted clinical NGS assay for common solid and hematologic cancers with high sensitivity, high specificity, and the flexibility to test very limited tissue samples often encountered in routine clinical practice.

Although next-generation sequencing (NGS) can revolutionize molecular diagnostics, several hurdles remain in the implementation of this technology in clinical laboratories.To validate and implement an NGS panel for genetic diagnosis of more than 100 inherited diseases, such as neurologic conditions, congenital hearing loss and eye disorders, developmental disorders, nonmalignant diseases treated by hematopoietic cell transplantation, familial cancers, connective tissue disorders, metabolic disorders, disorders of sexual development, and cardiac disorders. The diagnostic gene panels ranged from 1 to 54 genes with most of panels containing 10 genes or fewer.We used a liquid hybridization-based, target-enrichment strategy to enrich 10 067 exons in 568 genes, followed by NGS with a HiSeq 2000 sequencing system (Illumina, San Diego, California).We successfully sequenced 97.6% (9825 of 10 067) of the targeted exons to obtain a minimum coverage of 20× at all bases. We demonstrated 100% concordance in detecting 19 pathogenic single-nucleotide variations and 11 pathogenic insertion-deletion mutations ranging in size from 1 to 18 base pairs across 18 samples that were previously characterized by Sanger sequencing. Using 4 pairs of blinded, duplicate samples, we demonstrated a high degree of concordance (99%) among the blinded, duplicate pairs.We have successfully demonstrated the feasibility of using the NGS platform to multiplex genetic tests for several rare diseases and the use of cloud computing for bioinformatics analysis as a relatively low-cost solution for implementing NGS in clinical laboratories.

A study of genome-wide gene expression in major depressive disorder (MDD) was undertaken in a large population-based sample to determine whether altered expression levels of genes and pathways could provide insights into biological mechanisms that are relevant to this disorder. Gene expression studies have the potential to detect changes that may be because of differences in common or rare genomic sequence variation, environmental factors or their interaction. We recruited a European ancestry sample of 463 individuals with recurrent MDD and 459 controls, obtained self-report and semi-structured interview data about psychiatric and medical history and other environmental variables, sequenced RNA from whole blood and genotyped a genome-wide panel of common single-nucleotide polymorphisms. We used analytical methods to identify MDD-related genes and pathways using all of these sources of information. In analyses of association between MDD and expression levels of 13 857 single autosomal genes, accounting for multiple technical, physiological and environmental covariates, a significant excess of low P-values was observed, but there was no significant single-gene association after genome-wide correction. Pathway-based analyses of expression data detected significant association of MDD with increased expression of genes in the interferon α/β signaling pathway. This finding could not be explained by potentially confounding diseases and medications (including antidepressants) or by computationally estimated proportions of white blood cell types. Although cause-effect relationships cannot be determined from these data, the results support the hypothesis that altered immune signaling has a role in the pathogenesis, manifestation, and/or the persistence and progression of MDD.