Faculty Bibliography

BACKGROUND:Neurologic complications of the human immunodeficiency virus (HIV) are common in treated individuals and toxicity of certain antiretroviral therapies (ART) may contribute to cognitive impairment. We investigated exposures to specific ART and cognition among women living with HIV (WLWH). METHODS:Virologically suppressed (viral load <200 copies/mL during at least two semi-annual visits) WLWH and age/race matched HIV-seronegative controls enrolled in the Women's Interagency HIV Study (WIHS) who completed at least two biennial cognitive assessments were included. Analysis of WLWH was restricted to those with exposure to the drug class of interest and a nucleoside reverse transcriptase inhibitor backbone (NRTI). Generalized estimating equations were used to evaluate repeated measures of cognition over time in association with ART class exposure. RESULTS:Among 1,242 eligible WLWH, 20% (n=247) had isolated drug exposure to non-nucleoside reverse transcriptase inhibitors (NNRTI), 18% (n=219) to protease inhibitors (PI), and 6% (n=79) to integrase inhibitors (INSTI) with a NRTI backbone. Cognitive assessments were performed at a median of 3 biennial visits {IQR 2-4 visits}. At the index assessment, 21% of WLWH demonstrated global cognitive impairment versus 29% at their last cognitive assessment. In multivariable analyses adjusted, WLWH exposed to NNRTIs demonstrated verbal learning improvements (mean T-score change 1.3, p=0.020) as compared to other treated women. Compared to HIV-seronegative women, WLWH exposed to PIs had worse verbal learning (mean T-score difference -2.62, p=0.002) and verbal memory performance (mean T-score difference -1.74, p=0.032) at baseline. Compared to HIV-seronegative women, WLWH exposed to PIs had improvements in verbal learning (mean T-score slope difference 0.36, p=0.025) and verbal memory (mean T-score slope difference 0.32, p=0.042). CONCLUSIONS:We noted emerging trends in cognition in WLWH exposed to specific ART. Ongoing study of this relatively young group is important to characterize long term cognitive outcomes and effect of antiretrovirals as treatment guidelines evolve.

BACKGROUND:Alterations in gut microbiota and blood metabolomic profiles have been implicated in HIV infection and cardiovascular disease. However, it remains unclear whether alterations in gut microbiota may contribute to disrupted host blood metabolomic profiles in relation to atherosclerosis, especially in the context of HIV infection. METHODS:We analyzed cross-sectional associations between gut microbiota features and carotid artery plaque in 361 women with or at high risk of HIV (67% HIV+), and further integrated plaque-associated microbial features with plasma lipidomic/metabolomic profiles. Furthermore, in 737 women and men, we examined prospective associations of baseline gut bacteria-associated lipidomic and metabolomic profiles with incident carotid artery plaque over 7-year follow-up. RESULTS:were associated with multiple lipids/metabolites which were clustered into 8 modules in network. A module comprised of 9 lysophosphatidylcholines and lysophosphatidylethanolamines and a module comprised of 9 diglycerides were associated with increased risk of carotid artery plaque (risk ratio [95% CI], 1.34 [1.09-1.64] and 1.24 [1.02-1.51] per SD increment, respectively). Functional analyses identified bacterial enzymes in lipid metabolism associated with these plasma lipids. In particular, phospholipase A1 and A2 are the key enzymes in the reactions producing lysophosphatidylcholines and lysophosphatidylethanolamines. CONCLUSIONS:Among individuals with or at high risk of HIV infection, we identified altered gut microbiota and related functional capacities in the lipid metabolism associated with disrupted plasma lipidomic profiles and carotid artery atherosclerosis.

MicroRNAs (miRNAs) may be considered important regulators of risk for type 2 diabetes (T2D). The aim of the present study was to identify novel sets of miRNAs associated with T2D risk, as well as their gene and pathway targets. Circulating miRNAs (n=59) were measured in plasma from participants in a previously completed clinical trial (n=82). An agnostic statistical approach was applied to identify novel sets of miRNAs with optimal co‑expression patterns. In silico analyses were used to identify the messenger RNA and biological pathway targets of the miRNAs within each factor. A total of three factors of miRNAs were identified, containing 18, seven and two miRNAs each. Eight biological pathways were revealed to contain genes targeted by the miRNAs in all three factors, 38 pathways contained genes targeted by the miRNAs in two factors, and 55, 18 and two pathways were targeted by the miRNAs in a single factor, respectively (all q<0.05). The pathways containing genes targeted by miRNAs in the largest factor shared a common theme of biological processes related to metabolism and inflammation. By contrast, the pathways containing genes targeted by miRNAs in the second largest factor were related to endocrine function and hormone activity. The present study focused on the pathways uniquely targeted by each factor of miRNAs in order to identify unique mechanisms that may be associated with a subset of individuals. Further exploration of the genes and pathways related to these biological themes may provide insights about the subtypes of T2D and lead to the identification of novel therapeutic targets.

BACKGROUND:Assessing the effect of alcohol consumption on biological age is essential for understanding alcohol use-related comorbidities and mortality. Previously developed epigenetic clocks are mainly based on DNA methylation in heterogeneous cell types, which provide limited knowledge on the impacts of alcohol consumption at the individual cellular level. Evidence shows that monocytes play an important role in both alcohol-induced pathophysiology and the aging process. In this study, we developed a novel monocyte-based DNA methylation clock (MonoDNAmAge) to assess the impact of alcohol consumption on monocyte age. METHODS: = 2242). Using the MonoDNAmAge clock and four established clocks (i.e., HorvathDNAmAge, HannumDNAmAge, PhenoDNAmAge, GrimDNAmAge), we then evaluated the effect of alcohol consumption on epigenetic aging in the three cohorts [i.e., Yale Stress Center Community Study (YSCCS), Veteran Aging Cohort Study (VACS), Women's Interagency HIV Study (WIHS)] using linear and quadratic models. RESULTS: CONCLUSIONS:The results suggest a nonlinear relationship between alcohol consumption and its effects on epigenetic age. Considering adverse effects of alcohol consumption on health, nonlinear effects of alcohol use should be interpreted with caution. The findings, for the first time, highlight the complex effects of alcohol consumption on biological aging.

Here we report three epigenome-wide association studies (EWAS) of DNA methylation on self-reported race, global genetic ancestry, and local genetic ancestry in admixed Americans from three sets of samples, including internal and external replications (N_{total }= 1224). Our EWAS on local ancestry (LA) identified the largest number of ancestry-associated DNA methylation sites and also featured the highest replication rate. Furthermore, by incorporating ancestry origins of genetic variations, we identified 36 methylation quantitative trait loci (meQTL) clumps for LA-associated CpGs that cannot be captured by a model that assumes identical genetic effects across ancestry origins. Lead SNPs at 152 meQTL clumps had significantly different genetic effects in the context of an African or European ancestry background. Local ancestry information enables superior capture of ancestry-associated methylation signatures and identification of ancestry-specific genetic effects on DNA methylation. These findings highlight the importance of incorporating local ancestry for EWAS in admixed samples from multi-ancestry cohorts.

BACKGROUND:Gene expression is regulated by transcription factors, cofactors, and epigenetic mechanisms. Coexpressed genes indicate similar functional categories and gene networks. Detecting gene-gene coexpression is important for understanding the underlying mechanisms of cellular function and human diseases. A common practice of identifying coexpressed genes is to test the correlation of expression in a set of genes. In single-cell RNA-seq data, an important challenge is the abundance of zero values, so-called "dropout", which results in biased estimation of gene-gene correlations for downstream analyses. In recent years, efforts have been made to recover coexpressed genes in scRNA-seq data. Here, our goal is to detect coexpressed gene pairs to reduce the "dropout" effect in scRNA-seq data using a novel graph-based k-partitioning method by merging transcriptomically similar cells. RESULTS:We observed that the number of zero values was reduced among the merged transcriptomically similar cell clusters. Motivated by this observation, we leveraged a graph-based algorithm and develop an R package, scCorr, to recover the missing gene-gene correlation in scRNA-seq data that enables the reliable acquisition of cluster-based gene-gene correlations in three independent scRNA-seq datasets. The graphically partitioned cell clusters did not change the local cell community. For example, in scRNA-seq data from peripheral blood mononuclear cells (PBMCs), the gene-gene correlation estimated by scCorr outperformed the correlation estimated by the nonclustering method. Among 85 correlated gene pairs in a set of 100 clusters, scCorr detected 71 gene pairs, while the nonclustering method detected only 4 pairs of a dataset from PBMCs. The performance of scCorr was comparable to those of three previously published methods. As an example of downstream analysis using scCorr, we show that scCorr accurately identified a known cell type (i.e., CD4+ T cells) in PBMCs with a receiver operating characteristic area under the curve of 0.96. CONCLUSIONS:Our results demonstrate that scCorr is a robust and reliable graph-based method for identifying correlated gene pairs, which is fundamental to network construction, gene-gene interaction, and cellular omic analyses. scCorr can be quickly and easily implemented to minimize zero values in scRNA-seq analysis and is freely available at https://github.com/CBIIT-CGBB/scCorr .

Prior studies focused on circulating microRNAs and the risk for complex diseases have shown inconsistent findings. The majority of studies focused on European and East Asian racial or ethnic groups, however, ancestry was not typically reported. We evaluated the risk for type 2 diabetes as an exemplar to show that race and ethnic group may contribute to inconsistent validation of previous findings of associations with microRNAs.

PURPOSE:Gestational diabetes (GDM) is associated with increased risk for preterm birth and related complications for both the pregnant person and newborn. Changes in gene expression have the potential to characterize complex interactions between genetic and behavioral/environmental risk factors for GDM. Our goal was to summarize the state of the science about changes in gene expression and GDM. DESIGN:The systematic review was conducted using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. METHODS:PubMed articles about humans, in English, from any date were included if they described mRNA transcriptome or microRNA findings from blood samples in adults with GDM compared with adults without GDM. RESULTS:Sixteen articles were found representing 1355 adults (n=674 with GDM, n=681 controls) from 12 countries. Three studies reported transcriptome results and thirteen reported microRNA findings. Identified pathways described various aspects of diabetes pathogenesis, including glucose and insulin signaling, regulation, and transport; natural killer cell mediated cytotoxicity; and fatty acid biosynthesis and metabolism. Studies described 135 unique miRNAs that were associated with GDM, of which eight (miR-16-5p, miR-17-5p, miR-20a-5p, miR-29a-3p, miR-195-5p, miR-222-3p, miR-210-3p, and miR-342-3p) were described in 2 or more studies. Findings suggest that miRNA levels vary based on the time in pregnancy when GDM develops, the time point at which they were measured, sex assigned at birth of the offspring, and both the pre-pregnancy and gestational body mass index of the pregnant person. CONCLUSIONS:The mRNA, miRNA, gene targets, and pathways identified in this review contribute to our understanding of GDM pathogenesis; however, further research is warranted to validate previous findings. In particular, longitudinal repeated-measures designs are needed that control for participant characteristics (e.g., weight), use standardized data collection methods and analysis tools, and are sufficiently powered to detect differences between subgroups. Findings may be used to improve early diagnosis, prevention, medication choice and/or clinical treatment of patients with GDM.