Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects.

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

Genetic interactions drive heterogeneity in causal variant effect sizes for gene expression and complex traits.

Despite the growing number of genome-wide association studies (GWASs), it remains unclear to what extent gene-by-gene and gene-by-environment interactions influence complex traits in humans. The magnitude of genetic interactions in complex traits has been difficult to quantify because GWASs are generally underpowered to detect individual interactions of small effect. Here, we develop a method to test for genetic interactions that aggregates information across all trait-associated loci. Specifically, we test whether SNPs in regions of European ancestry shared between European American and admixed African American individuals have the same causal effect sizes. We hypothesize that in African Americans, the presence of genetic interactions will drive the causal effect sizes of SNPs in regions of European ancestry to be more similar to those of SNPs in regions of African ancestry. We apply our method to two traits: gene expression in 296 African Americans and 482 European Americans in the Multi-Ethnic Study of Atherosclerosis (MESA) and low-density lipoprotein cholesterol (LDL-C) in 74K African Americans and 296K European Americans in the Million Veteran Program (MVP). We find significant evidence for genetic interactions in our analysis of gene expression; for LDL-C, we observe a similar point estimate, although this is not significant, most likely due to lower statistical power. These results suggest that gene-by-gene or gene-by-environment interactions modify the effect sizes of causal variants in human complex traits.

DNA demethylase activity maintains intestinal cells in an undifferentiated state following loss of APC.

Although genome-wide hypomethylation is a hallmark of many cancers, roles for active DNA demethylation during tumorigenesis are unknown. Here, loss of the APC tumor suppressor gene causes upregulation of a DNA demethylase system and the concomitant hypomethylation of key intestinal cell fating genes. Notably, this hypomethylation maintained zebrafish intestinal cells in an undifferentiated state that was released upon knockdown of demethylase components. Mechanistically, the demethylase genes are directly activated by Pou5f1 and Cebpß and are indirectly repressed by retinoic acid, which antagonizes Pou5f1 and Cebpß. Apc mutants lack retinoic acid as a result of the transcriptional repression of retinol dehydrogenase l1 via a complex that includes Lef1, Groucho2, Ctbp1, Lsd1, and Corest. Our findings imply a model wherein APC controls intestinal cell fating through a switch in DNA methylation dynamics. Wild-type APC and retinoic acid downregulate demethylase components, thereby promoting DNA methylation of key genes and helping progenitors commit to differentiation.

Overcoming the impacts of two-step batch effect correction on gene expression estimation and inference.

Nonignorable technical variation is commonly observed across data from multiple experimental runs, platforms, or studies. These so-called batch effects can lead to difficulty in merging data from multiple sources, as they can severely bias the outcome of the analysis. Many groups have developed approaches for removing batch effects from data, usually by accommodating batch variables into the analysis (one-step correction) or by preprocessing the data prior to the formal or final analysis (two-step correction). One-step correction is often desirable due it its simplicity, but its flexibility is limited and it can be difficult to include batch variables uniformly when an analysis has multiple stages. Two-step correction allows for richer models of batch mean and variance. However, prior investigation has indicated that two-step correction can lead to incorrect statistical inference in downstream analysis. Generally speaking, two-step approaches introduce a correlation structure in the corrected data, which, if ignored, may lead to either exaggerated or diminished significance in downstream applications such as differential expression analysis. Here, we provide more intuitive and more formal evaluations of the impacts of two-step batch correction compared to existing literature. We demonstrate that the undesired impacts of two-step correction (exaggerated or diminished significance) depend on both the nature of the study design and the batch effects. We also provide strategies for overcoming these negative impacts in downstream analyses using the estimated correlation matrix of the corrected data. We compare the results of our proposed workflow with the results from other published one-step and two-step methods and show that our methods lead to more consistent false discovery controls and power of detection across a variety of batch effect scenarios. Software for our method is available through GitHub (https://github.com/jtleek/sva-devel) and will be available in future versions of the $\texttt{sva}$ R package in the Bioconductor project (https://bioconductor.org/packages/release/bioc/html/sva.html).

Psychosocial stressors and breast cancer gene expression in the Black Women's Health Study.

PURPOSE: Prior studies indicate that the physiologic response to stress can affect gene expression. We evaluated differential gene expression in breast cancers collected from Black women with high versus low exposure to psychosocial stressors. METHODS: We analyzed tumor RNA sequencing data from 417 Black Women's Health Study breast cancer cases with data on early life trauma and neighborhood disadvantage. We conducted age-adjusted differential gene expression analyses and pathway analyses. We also evaluated Conserved Transcriptional Response to Adversity (CTRA) contrast scores, relative fractions of immune cell types, T cell exhaustion, and adrenergic signaling. Analyses were run separately for estrogen receptor positive (ER+; n = 299) and ER- (n = 118) cases. RESULTS: Among ER+ cases, the top differentially expressed pathways by stress exposure were related to RNA and protein metabolism. Among ER- cases, they were related to developmental biology, signal transduction, metabolism, and the immune system. Targeted analyses indicated greater immune pathway enrichment with stress exposure for ER- cases, and possible relevance of adrenergic signaling for ER+ cases. CTRA contrast scores did not differ by stress exposure, but in analyses of the CTRA components, ER- breast cancer cases with high neighborhood disadvantage had higher pro-inflammatory gene expression (p = 0.039) and higher antibody gene expression (p = 0.006) compared to those with low neighborhood disadvantage. CONCLUSION: There are multiple pathways through which psychosocial stress exposure may influence breast tumor biology. Given the present findings on inflammation and immune response in ER- tumors, further research to identify stress-induced changes in the etiology and progression of ER- breast cancer is warranted.

Differences in Metabolomic Profiles Between Black and White Women and Risk of Coronary Heart Disease: an Observational Study of Women From Four US Cohorts.

BACKGROUND: Racial differences in metabolomic profiles may reflect underlying differences in social determinants of health by self-reported race and may be related to racial disparities in coronary heart disease (CHD) among women in the United States. However, the magnitude of differences in metabolomic profiles between Black and White women in the United States has not been well-described. It also remains unknown whether such differences are related to differences in CHD risk. METHODS: Plasma metabolomic profiles were analyzed using liquid chromatography-tandem mass spectrometry in the WHI-OS (Women's Health Initiative-Observational Study; 138 Black and 696 White women), WHI-HT trials (WHI-Hormone Therapy; 156 Black and 1138 White women), MESA (Multi-Ethnic Study of Atherosclerosis; 114 Black and 219 White women), JHS (Jackson Heart Study; 1465 Black women with 107 incident CHD cases), and NHS (Nurses' Health Study; 2506 White women with 136 incident CHD cases). First, linear regression models were used to estimate associations between self-reported race and 472 metabolites in WHI-OS (discovery); findings were replicated in WHI-HT and validated in MESA. Second, we used elastic net regression to construct a racial difference metabolomic pattern (RDMP) representing differences in the metabolomic patterns between Black and White women in the WHI-OS; the RDMP was validated in the WHI-HT and MESA. Third, using conditional logistic regressions in the WHI (717 CHD cases and 719 matched controls), we examined associations of metabolites with large differences in levels by race and the RDMP with risk of CHD, and the results were replicated in Black women from the JHS and White women from the NHS. RESULTS: Of the 472 tested metabolites, levels of 259 (54.9%) metabolites, mostly lipid metabolites and amino acids, significantly differed between Black and White women in both WHI-OS and WHI-HT after adjusting for baseline characteristics, socioeconomic status, lifestyle factors, baseline health conditions, and medication use (false discovery rate <0.05); similar trends were observed in MESA. The RDMP, composed of 152 metabolites, was identified in the WHI-OS and showed significantly different distributions between Black and White women in the WHI-HT and MESA. Higher RDMP quartiles were associated with an increased risk of incident CHD (odds ratio=1.51 [0.97-2.37] for the highest quartile comparing to the lowest; Ptrend=0.02), independent of self-reported race and known CHD risk factors. In race-stratified analyses, the RDMP-CHD associations were more pronounced in White women. Similar patterns were observed in Black women from the JHS and White women from the NHS. CONCLUSIONS: Metabolomic profiles significantly and substantially differ between Black and White women and may be associated with CHD risk and racial disparities in US women.

Arsenic Exposure, Blood DNA Methylation, and Cardiovascular Disease.

BACKGROUND: Epigenetic dysregulation has been proposed as a key mechanism for arsenic-related cardiovascular disease (CVD). We evaluated differentially methylated positions (DMPs) as potential mediators on the association between arsenic and CVD. METHODS: Blood DNA methylation was measured in 2321 participants (mean age 56.2, 58.6% women) of the Strong Heart Study, a prospective cohort of American Indians. Urinary arsenic species were measured using high-performance liquid chromatography coupled to inductively coupled plasma mass spectrometry. We identified DMPs that are potential mediators between arsenic and CVD. In a cross-species analysis, we compared those DMPs with differential liver DNA methylation following early-life arsenic exposure in the apoE knockout (apoE-/-) mouse model of atherosclerosis. RESULTS: A total of 20 and 13 DMPs were potential mediators for CVD incidence and mortality, respectively, several of them annotated to genes related to diabetes. Eleven of these DMPs were similarly associated with incident CVD in 3 diverse prospective cohorts (Framingham Heart Study, Women's Health Initiative, and Multi-Ethnic Study of Atherosclerosis). In the mouse model, differentially methylated regions in 20 of those genes and DMPs in 10 genes were associated with arsenic. CONCLUSIONS: Differential DNA methylation might be part of the biological link between arsenic and CVD. The gene functions suggest that diabetes might represent a relevant mechanism for arsenic-related cardiovascular risk in populations with a high burden of diabetes.

Comparison of gene set scoring methods for reproducible evaluation of tuberculosis gene signatures.

BACKGROUND: Blood-based transcriptional gene signatures for tuberculosis (TB) have been developed with potential use to diagnose disease. However, an unresolved issue is whether gene set enrichment analysis of the signature transcripts alone is sufficient for prediction and differentiation or whether it is necessary to use the original model created when the signature was derived. Intra-method comparison is complicated by the unavailability of original training data and missing details about the original trained model. To facilitate the utilization of these signatures in TB research, comparisons between gene set scoring methods cross-data validation of original model implementations are needed. METHODS: We compared the performance of 19 TB gene signatures across 24 transcriptomic datasets using both rrebuilt original models and gene set scoring methods. Existing gene set scoring methods, including ssGSEA, GSVA, PLAGE, Singscore, and Zscore, were used as alternative approaches to obtain the profile scores. The area under the ROC curve (AUC) value was computed to measure performance. Correlation analysis and Wilcoxon paired tests were used to compare the performance of enrichment methods with the original models. RESULTS: For many signatures, the predictions from gene set scoring methods were highly correlated and statistically equivalent to the results given by the original models. In some cases, PLAGE outperformed the original models when considering signatures' weighted mean AUC values and the AUC results within individual studies. CONCLUSION: Gene set enrichment scoring of existing gene sets can distinguish patients with active TB disease from other clinical conditions with equivalent or improved accuracy compared to the original methods and models. These data justify using gene set scoring methods of published TB gene signatures for predicting TB risk and treatment outcomes, especially when original models are difficult to apply or implement.

Hepatotoxicity of silver nanoparticles: Benchmark concentration modeling of an in vitro transcriptomics study in human iPSC-derived hepatocytes.

Despite two decades of research on silver nanoparticle (AgNP) toxicity, a safe threshold for exposure has not yet been established, albeit being critically needed for risk assessment and regulatory decision-making. Traditionally, a point-of-departure (PoD) value is derived from dose response of apical endpoints in animal studies using either the no-observed-adverse-effect level (NOAEL) approach, or benchmark dose (BMD) modeling. To develop new approach methodologies (NAMs) to inform human risk assessment of AgNPs, we conducted a concentration response modeling of the transcriptomic changes in hepatocytes derived from human induced pluripotent stem cells (iPSCs) after being exposed to a wide range concentration (0.01-25 µg/ml) of AgNPs for 24 h. A plausible transcriptomic PoD of 0.21 µg/ml was derived for a pathway related to the mode-of-action (MOA) of AgNPs, and a more conservative PoD of 0.10 µg/ml for a gene ontology (GO) term not apparently associated with the MOA of AgNPs. A reference dose (RfD) could be calculated from either of the PoDs as a safe threshold for AgNP exposure. The current study illustrates the usefulness of in vitro transcriptomic concentration response study using human cells as a NAM for toxicity study of chemicals that lack adequate toxicity data to inform human risk assessment.

Trends in childhood body mass index between 1936 and 2011 showed that underweight remained more common than obesity among 398 970 Danish school children.

AIM: To examine trends in all body mass index (BMI) groups in children from 1936 to 2011. METHODS: We included 197 694 girls and 201 276 boys from the Copenhagen School Health Records Register, born between 1930 and 1996, with longitudinal weight and height measurements (6-14 years). Using International Obesity Task Force criteria, BMI was classified as underweight, normal-weight, overweight and obesity. Sex- and age-specific prevalences were calculated. RESULTS: From the 1930s, the prevalence of underweight was stable until a small increase occurred from 1950 to 1970s, and thereafter it declined into the early 2000s. Using 7-year-olds as an example, underweight changed from 10% to 7% in girls and from 9% to 6% in boys during the study period. The prevalence of overweight plateaued from 1950 to 1970s and then steeply increased from 1970s onwards and in 1990-2000s 15% girls and 11% boys at 7 years had overweight. The prevalence of obesity particularly increased from 1980s onwards and in 1990-2000s 5% girls and 4% boys at 7 years had obesity. These trends slightly differed by age. CONCLUSION: Among Danish schoolchildren, the prevalence of underweight was greater than overweight until the 1980s and greater than obesity throughout the period. Thus, monitoring the prevalence of childhood underweight remains an important public health issue.

Non-parental Childcare During Early Childhood and Problem Behaviour Trajectories from Ages 5 to 14 Years.

Using data from the nationally representative Millennium Cohort Study, this study examined the association between age of starting and weekly hours in formal childcare between birth and 5 years with internalising and externalising behaviour trajectories from ages 5 to 14 years in England (N = 6194 children). Associations were analysed using multilevel general linear regression models, with adjustment for socio-economic position, maternal mental health, demographics, and child temperament. Later entry was associated with more internalising behaviours at age 14 years. Children who spent > 40 h per week in childcare between birth and 3 years displayed more externalising behaviour at 5 years than children who did not attend childcare. Controlling for socio-economic position and parental mental health attenuated findings.

AtheroSpectrum Reveals Novel Macrophage Foam Cell Gene Signatures Associated With Atherosclerotic Cardiovascular Disease Risk.

BACKGROUND: Whereas several interventions can effectively lower lipid levels in people at risk for atherosclerotic cardiovascular disease (ASCVD), cardiovascular event risks remain, suggesting an unmet medical need to identify factors contributing to cardiovascular event risk. Monocytes and macrophages play central roles in atherosclerosis, but studies have yet to provide a detailed view of macrophage populations involved in increased ASCVD risk. METHODS: A novel macrophage foaming analytics tool, AtheroSpectrum, was developed using 2 quantitative indices depicting lipid metabolism and the inflammatory status of macrophages. A machine learning algorithm was developed to analyze gene expression patterns in the peripheral monocyte transcriptome of MESA participants (Multi-Ethnic Study of Atherosclerosis; set 1; n=911). A list of 30 genes was generated and integrated with traditional risk factors to create an ASCVD risk prediction model (30-gene cardiovascular disease risk score [CR-30]), which was subsequently validated in the remaining MESA participants (set 2; n=228); performance of CR-30 was also tested in 2 independent human atherosclerotic tissue transcriptome data sets (GTEx [Genotype-Tissue Expression] and GSE43292). RESULTS: Using single-cell transcriptomic profiles (GSE97310, GSE116240, GSE97941, and FR-FCM-Z23S), AtheroSpectrum detected 2 distinct programs in plaque macrophages-homeostatic foaming and inflammatory pathogenic foaming-the latter of which was positively associated with severity of atherosclerosis in multiple studies. A pool of 2209 pathogenic foaming genes was extracted and screened to select a subset of 30 genes correlated with cardiovascular event in MESA set 1. A cardiovascular disease risk score model (CR-30) was then developed by incorporating this gene set with traditional variables sensitive to cardiovascular event in MESA set 1 after cross-validation generalizability analysis. The performance of CR-30 was then tested in MESA set 2 (P=2.60×10-4; area under the receiver operating characteristic curve, 0.742) and 2 independent data sets (GTEx: P=7.32×10-17; area under the receiver operating characteristic curve, 0.664; GSE43292: P=7.04×10-2; area under the receiver operating characteristic curve, 0.633). Model sensitivity tests confirmed the contribution of the 30-gene panel to the prediction model (likelihood ratio test; df=31, P=0.03). CONCLUSIONS: Our novel computational program (AtheroSpectrum) identified a specific gene expression profile associated with inflammatory macrophage foam cells. A subset of 30 genes expressed in circulating monocytes jointly contributed to prediction of symptomatic atherosclerotic vascular disease. Incorporating a pathogenic foaming gene set with known risk factors can significantly strengthen the power to predict ASCVD risk. Our programs may facilitate both mechanistic investigations and development of therapeutic and prognostic strategies for ASCVD risk.

Racial and Ethnic Differences in All-Cause and Cardiovascular Disease Mortality: The MESA Study.

BACKGROUND: Despite improvements in population health, marked racial and ethnic disparities in longevity and cardiovascular disease (CVD) mortality persist. This study aimed to describe risks for all-cause and CVD mortality by race and ethnicity, before and after accounting for socioeconomic status (SES) and other factors, in the MESA study (Multi-Ethnic Study of Atherosclerosis). METHODS: MESA recruited 6814 US adults, 45 to 84 years of age, free of clinical CVD at baseline, including Black, White, Hispanic, and Chinese individuals (2000-2002). Using Cox proportional hazards modeling with time-updated covariates, we evaluated the association of self-reported race and ethnicity with all-cause and adjudicated CVD mortality, with progressive adjustments for age and sex, SES (neighborhood SES, income, education, and health insurance), lifestyle and psychosocial risk factors, clinical risk factors, and immigration history. RESULTS: During a median of 15.8 years of follow-up, 22.8% of participants (n=1552) died, of which 5.3% (n=364) died of CVD. After adjusting for age and sex, Black participants had a 34% higher mortality hazard (hazard ratio [HR], 1.34 [95% CI, 1.19-1.51]), Chinese participants had a 21% lower mortality hazard (HR, 0.79 [95% CI, 0.66-0.95]), and there was no mortality difference in Hispanic participants (HR, 0.99 [95% CI, 0.86-1.14]) compared with White participants. After adjusting for SES, the mortality HR for Black participants compared with White participants was reduced (HR, 1.16 [95% CI, 1.01-1.34]) but still statistically significant. With adjustment for SES, the mortality hazards for Chinese and Hispanic participants also decreased in comparison with White participants. After further adjustment for additional risk factors and immigration history, Hispanic participants (HR, 0.77 [95% CI, 0.63-0.94]) had a lower mortality risk than White participants, and hazard ratios for Black participants (HR, 1.08 [95% CI, 0.92-1.26]) and Chinese participants (HR, 0.81 [95% CI, 0.60-1.08]) were not significantly different from those of White participants. Similar trends were seen for CVD mortality, although the age- and sex-adjusted HR for CVD mortality for Black participants compared with White participants was greater than all-cause mortality (HR, 1.72 [95% CI, 1.34-2.21] compared with HR, 1.34 [95% CI, 1.19-1.51]). CONCLUSIONS: These results highlight persistent racial and ethnic differences in overall and CVD mortality, largely attributable to social determinants of health, and support the need to identify and act on systemic factors that shape differences in health across racial and ethnic groups.

Whole Genome Sequence Analysis of the Plasma Proteome in Black Adults Provides Novel Insights Into Cardiovascular Disease.

BACKGROUND: Plasma proteins are critical mediators of cardiovascular processes and are the targets of many drugs. Previous efforts to characterize the genetic architecture of the plasma proteome have been limited by a focus on individuals of European descent and leveraged genotyping arrays and imputation. Here we describe whole genome sequence analysis of the plasma proteome in individuals with greater African ancestry, increasing our power to identify novel genetic determinants. METHODS: Proteomic profiling of 1301 proteins was performed in 1852 Black adults from the Jackson Heart Study using aptamer-based proteomics (SomaScan). Whole genome sequencing association analysis was ascertained for all variants with minor allele count ≥5. Results were validated using an alternative, antibody-based, proteomic platform (Olink) as well as replicated in the Multi-Ethnic Study of Atherosclerosis and the HERITAGE Family Study (Health, Risk Factors, Exercise Training and Genetics). RESULTS: We identify 569 genetic associations between 479 proteins and 438 unique genetic regions at a Bonferroni-adjusted significance level of 3.8×10-11. These associations include 114 novel locus-protein relationships and an additional 217 novel sentinel variant-protein relationships. Novel cardiovascular findings include new protein associations at the APOE gene locus including ZAP70 (sentinel single nucleotide polymorphism [SNP] rs7412-T, ß=0.61±0.05, P=3.27×10-30) and MMP-3 (ß=-0.60±0.05, P=1.67×10-32), as well as a completely novel pleiotropic locus at the HPX gene, associated with 9 proteins. Further, the associations suggest new mechanisms of genetically mediated cardiovascular disease linked to African ancestry; we identify a novel association between variants linked to APOL1-associated chronic kidney and heart disease and the protein CKAP2 (rs73885319-G, ß=0.34±0.04, P=1.34×10-17) as well as an association between ATTR amyloidosis and RBP4 levels in community-dwelling individuals without heart failure. CONCLUSIONS: Taken together, these results provide evidence for the functional importance of variants in non-European populations, and suggest new biological mechanisms for ancestry-specific determinants of lipids, coagulation, and myocardial function.

Risk factor associations with individual myocardial infarction subtypes and acute non-ischemic myocardial injury in the Multi-Ethnic Study of Atherosclerosis (MESA): Design and rationale.

BACKGROUND: Despite different prevalence, pathobiology, and prognosis between etiologically distinct myocardial infarction (MI) subtypes, prospective study of risk factor for MI in large NHLBI-sponsored cardiovascular cohorts is limited to acute MI as a singular entity. Therefore, we sought to utilize the Multi-Ethnic Study of Atherosclerosis (MESA), a large prospective primary prevention cardiovascular study, to define the incidence and risk factor profile of individual myocardial injury subtypes. METHODS: We describe the rationale and design of re-adjudicating 4,080 events that occurred over the first 14 years of follow-up in MESA for the presence and subtype of myocardial injury as defined by the Fourth Universal Definition of MI: MI type 1 to 5, acute non-ischemic myocardial injury, and chronic myocardial injury. The project utilizes a 2-physician adjudication process via examination of medical records, abstracted data collection forms, cardiac biomarker results, and electrocardiograms of all relevant clinical events. Comparison of the magnitude and direction of associations between baseline traditional and novel cardiovascular risk factors with incident and recurrent acute MI subtypes and acute non-ischemic myocardial injury events will be made. CONCLUSIONS: This project will result in one of the first large prospective cardiovascular cohort with modern classification of acute MI subtypes, as well as a full accounting of non-ischemic myocardial injury events, with implications for numerous ongoing and future studies in MESA. By creating precise MI phenotypes, and defining their epidemiology, this project will allow for discovery of novel pathobiology-specific risk factors, allow for development of more accurate risk prediction, and suggest more targeted preventive strategies.

Gene expression associations with body mass index in the Multi-Ethnic Study of Atherosclerosis.

BACKGROUND/OBJECTIVES: Obesity, defined as excessive fat accumulation that represents a health risk, is increasing in adults and children, reaching global epidemic proportions. Body mass index (BMI) correlates with body fat and future health risk, yet differs in prediction by fat distribution, across populations and by age. Nonetheless, few genetic studies of BMI have been conducted in ancestrally diverse populations. Gene expression association with BMI was assessed in the Multi-Ethnic Study of Atherosclerosis (MESA) in four self-identified race and ethnicity (SIRE) groups to identify genes associated with obesity. SUBJECTS/METHODS: RNA-sequencing was performed on 1096 MESA participants (37.8% white, 24.3% Hispanic, 28.4% African American, and 9.5% Chinese American) and linear models were used to assess the association of expression from each gene for its effect on BMI, adjusting for age, sex, sequencing center, study site, five expression and four genetic principal components in each self-identified race group. Sample-size-weighted meta-analysis was performed to identify genes with BMI-associated expression across ancestry groups. RESULTS: Within individual SIRE groups, there were zero to three genes whose expression is significantly (p < 1.97 × 10-6) associated with BMI. Across all groups, 45 genes were identified by meta-analysis whose expression was significantly associated with BMI, explaining 29.7% of BMI variation. The 45 genes are expressed in a variety of tissues and cell types and are enriched for obesity-related processes including erythrocyte function, oxygen binding and transport, and JAK-STAT signaling. CONCLUSIONS: We have identified genes whose expression is significantly associated with obesity in a multi-ethnic cohort. We have identified novel genes associated with BMI as well as confirmed previously identified genes from earlier genetic analyses. These novel genes and their biological pathways represent new targets for understanding the biology of obesity as well as new therapeutic intervention to reduce obesity and improve global public health.

Respiratory syncytial virus M2-1 protein associates non-specifically with viral messenger RNA and with specific cellular messenger RNA transcripts.

Respiratory syncytial virus (RSV) is a major cause of respiratory disease in infants and the elderly. RSV is a non-segmented negative strand RNA virus. The viral M2-1 protein plays a key role in viral transcription, serving as an elongation factor to enable synthesis of full-length mRNAs. M2-1 contains an unusual CCCH zinc-finger motif that is conserved in the related human metapneumovirus M2-1 protein and filovirus VP30 proteins. Previous biochemical studies have suggested that RSV M2-1 might bind to specific virus RNA sequences, such as the transcription gene end signals or poly A tails, but there was no clear consensus on what RSV sequences it binds. To determine if M2-1 binds to specific RSV RNA sequences during infection, we mapped points of M2-1:RNA interactions in RSV-infected cells at 8 and 18 hours post infection using crosslinking immunoprecipitation with RNA sequencing (CLIP-Seq). This analysis revealed that M2-1 interacts specifically with positive sense RSV RNA, but not negative sense genome RNA. It also showed that M2-1 makes contacts along the length of each viral mRNA, indicating that M2-1 functions as a component of the transcriptase complex, transiently associating with nascent mRNA being extruded from the polymerase. In addition, we found that M2-1 binds specific cellular mRNAs. In contrast to the situation with RSV mRNA, M2-1 binds discrete sites within cellular mRNAs, with a preference for A/U rich sequences. These results suggest that in addition to its previously described role in transcription elongation, M2-1 might have an additional role involving cellular RNA interactions.

Distinct COPD subtypes in former smokers revealed by gene network perturbation analysis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) varies significantly in symptomatic and physiologic presentation. Identifying disease subtypes from molecular data, collected from easily accessible blood samples, can help stratify patients and guide disease management and treatment. METHODS: Blood gene expression measured by RNA-sequencing in the COPDGene Study was analyzed using a network perturbation analysis method. Each COPD sample was compared against a learned reference gene network to determine the part that is deregulated. Gene deregulation values were used to cluster the disease samples. RESULTS: The discovery set included 617 former smokers from COPDGene. Four distinct gene network subtypes are identified with significant differences in symptoms, exercise capacity and mortality. These clusters do not necessarily correspond with the levels of lung function impairment and are independently validated in two external cohorts: 769 former smokers from COPDGene and 431 former smokers in the Multi-Ethnic Study of Atherosclerosis (MESA). Additionally, we identify several genes that are significantly deregulated across these subtypes, including DSP and GSTM1, which have been previously associated with COPD through genome-wide association study (GWAS). CONCLUSIONS: The identified subtypes differ in mortality and in their clinical and functional characteristics, underlining the need for multi-dimensional assessment potentially supplemented by selected markers of gene expression. The subtypes were consistent across cohorts and could be used for new patient stratification and disease prognosis.

Untargeted Metabolomics Analyses and Contaminant Chemistry of Dreissenid Mussels at the Maumee River Area of Concern in the Great Lakes.

Bivalves serve as an ideal ecological indicator; hence, their use by the NOAA Mussel Watch Program to monitor environmental health. This study aimed to expand the baseline knowledge of using metabolic end points in environmental monitoring by investigating the dreissenid mussel metabolome in the field. Dreissenids were caged at four locations along the Maumee River for 30 days. The mussel metabolome was measured using nuclear magnetic resonance spectroscopy, and mussel tissue chemical contaminants were analyzed using gas or liquid chromatography coupled with mass spectrometry. All Maumee River sites had a distinct mussel metabolome compared to the reference site and revealed changes in the energy metabolism and amino acids. Data also highlighted the importance of considering seasonality or handling effects on the metabolome at the time of sampling. The furthest upstream site presented a specific mussel tissue chemical signature of pesticides (atrazine and metolachlor), while a downstream site, located at Toledo's wastewater treatment plant, was characterized by polycyclic aromatic hydrocarbons and other organic contaminants. Further research into the dreissenid mussel's natural metabolic cycle and metabolic response to specific anthropogenic stressors is necessary before successful implementation of metabolomics in a biomonitoring program.

Development of prognostic scoring system for predicting 1-year mortality among pulmonary tuberculosis patients in South India.

BACKGROUND: Development of a prediction model using baseline characteristics of tuberculosis (TB) patients at the time of diagnosis will aid us in early identification of the high-risk groups and devise pertinent strategies accordingly. Hence, we did this study to develop a prognostic-scoring model for predicting the death among newly diagnosed drug sensitive pulmonary TB patients in South India. METHODS: We undertook a longitudinal analysis of cohort data under the Regional Prospective Observational Research for Tuberculosis India consortium. Multivariable cox regression using the stepwise backward elimination procedure was used to select variables for the model building and the nomogram-scoring system was developed with the final selected model. RESULTS: In total, 54 (4.6%) out of the 1181 patients had died during the 1-year follow-up period. The TB mortality rate was 0.20 per 1000 person-days. Eight variables (age, gender, functional limitation, anemia, leukopenia, thrombocytopenia, diabetes, neutrophil-lymphocyte ratio) were selected and a nomogram was built using these variables. The discriminatory power was 0.81 (95% confidence interval: 0.75-0.86) and this model was well-calibrated. Decision curve analysis showed that the model is beneficial at a threshold probability ~15-65%. CONCLUSIONS: This scoring system could help the clinicians and policy makers to devise targeted interventions and in turn reduce the TB mortality in India.