Single cell transcriptomic analyses reveal diverse and dynamic changes of distinct populations of lung interstitial macrophages in hypoxia-induced pulmonary hypertension.

Introduction: Hypoxia is a common pathological driver contributing to various forms of pulmonary vascular diseases leading to pulmonary hypertension (PH). Pulmonary interstitial macrophages (IMs) play pivotal roles in immune and vascular dysfunction, leading to inflammation, abnormal remodeling, and fibrosis in PH. However, IMs' response to hypoxia and their role in PH progression remain largely unknown. We utilized a murine model of hypoxia-induced PH to investigate the repertoire and functional profiles of IMs in response to acute and prolonged hypoxia, aiming to elucidate their contributions to PH development. Methods: We conducted single-cell transcriptomic analyses to characterize the repertoire and functional profiles of murine pulmonary IMs following exposure to hypobaric hypoxia for varying durations (0, 1, 3, 7, and 21 days). Hallmark pathways from the mouse Molecular Signatures Database were utilized to characterize the molecular function of the IM subpopulation in response to hypoxia. Results: Our analysis revealed an early acute inflammatory phase during acute hypoxia exposure (Days 1-3), which was resolved by Day 7, followed by a pro-remodeling phase during prolonged hypoxia (Days 7-21). These phases were marked by distinct subpopulations of IMs: MHCIIhiCCR2+EAR2+ cells characterized the acute inflammatory phase, while TLF+VCAM1hi cells dominated the pro-remodeling phase. The acute inflammatory phase exhibited enrichment in interferon-gamma, IL-2, and IL-6 pathways, while the pro-remodeling phase showed dysregulated chemokine production, hemoglobin clearance, and tissue repair profiles, along with activation of distinct complement pathways. Discussion: Our findings demonstrate the existence of distinct populations of pulmonary interstitial macrophages corresponding to acute and prolonged hypoxia exposure, pivotal in regulating the inflammatory and remodeling phases of PH pathogenesis. This understanding offers potential avenues for targeted interventions, tailored to specific populations and distinct phases of the disease. Moreover, further identification of triggers for pro-remodeling IMs holds promise in unveiling novel therapeutic strategies for pulmonary hypertension.

Single-cell RNA sequencing and binary hierarchical clustering define lung interstitial macrophage heterogeneity in response to hypoxia.

Few studies have examined lung interstitial macrophage (IM) molecular phenotypes after being exposed to hypoxia in vivo at the single-cell level, even though macrophages contribute to hypoxic pulmonary hypertension (PH). We aimed to determine IM diversity and its association with hypoxia-induced PH. We hypothesized that integrating single-cell RNA sequencing (scRNAseq) and binary hierarchal clustering (BHC) could resolve IM heterogeneity under normal homeostatic conditions and changes induced by hypoxia exposure. Cx3cr1GFP/+ reporter mice were exposed to normoxic conditions (∼21% [Formula: see text]) or exposed to 1 day (D1) or 7 days (D7) of hypoxia (∼10% [Formula: see text]). We used flow cytometry to isolate Cx3cr1+ IMs and the 10X Genomics platform for scRNAseq, Cell Ranger, Seurat, ClusterMap, monocle, ingenuity pathway analysis, and Fisher's exact test (q value < 0.05) for functional investigations. n = 374 (normoxia), n = 2,526 (D1), and n = 1,211 (D7) IMs were included in the analyses. We identified three normoxia-related cell types, five hypoxia-associated cell types that emerged at D1, and three that appeared at D7. We describe the existence of a putative resident trained innate IM, which is present in normoxia, transiently depleted at D1, and recovered after 7 days of sustained hypoxia. We also define a rare putative pathogenic population associated with transcripts implicated in PH development that emerges at D7. In closing, we describe the successful integration of BHC with scRNAseq to determine IM heterogeneity and its association with PH. These results shed light on how resident-trained innate IMs become more heterogeneous but ultimately accustomed to hypoxia.

"Omics" data integration and functional analyses link Enoyl-CoA hydratase, short chain 1 to drug refractory dilated cardiomyopathy.

BACKGROUND: Large-scale "omics" datasets have not been leveraged and integrated with functional analyses to discover potential drivers of cardiomyopathy. This study addresses the knowledge gap. METHODS: We coupled RNA sequence (RNA-Seq) variant detection and transcriptome profiling with pathway analysis to model drug refractory dilated cardiomyopathy (drDCM) using the BaseSpace sequencing hub and Ingenuity Pathway Analysis. We used RNA-Seq case-control datasets (n = 6 cases, n = 4 controls), exome sequence familial DCM datasets (n = 3 Italians, n = 5 Italians, n = 5 Chinese), and controls from the HapMap project (n = 5 Caucasians, and n = 5 Asians) for disease modeling and putative mutation discovery. Variant replication datasets: n = 128 cases and n = 15 controls. Source of datasets: NCBI Sequence Read Archive. STATISTICS: Pairwise differential expression analyses to determine differentially expressed genes and t-tests to calculate p-values. We adjusted for false discovery rates and reported q-values. We used chi-square tests to assess independence among variables, the Fisher's Exact Tests and overlap p-values for the pathways and p-scores to rank network. RESULTS: Data revealed that ECHS1(enoyl-CoA hydratase, short chain 1(log2(foldchange) = 1.63329) hosts a mirtron, MIR3944 expressed in drDCM (FPKM = 5.2857) and not in controls (FPKM = 0). Has-miR3944-3p is a putative target of BAG1 (BCL2 associated athanogene 1(log2(foldchange) = 1.31978) and has-miR3944-5p of ITGAV (integrin subunit alpha V(log2(foldchange) = 1.46107) and RHOD (ras homolog family member D(log2(foldchange) = 1.28851). There is an association between ECHS1:11 V/A(rs10466126) and drDCM (p = 0.02496). The interaction (p = 2.82E-07) between ECHS1:75 T/I(rs1049951) and ECHS1:rs10466126 is associated with drDCM (p < 2.2e-16). ECHS1:rs10466126 and ECHS1:rs1049951 are in linkage disequilibrium (D' = 1). The interaction (p = 7.84E-08) between ECHS1:rs1049951 and the novel ECHS1:c.41insT variant is associated with drDCM (p < 2.2e-16). The interaction (p = 0.001096) between DBT (Dihydrolipoamide branched chain transacylase E2):384G/S(rs12021720) and ECHS1:rs10466126 is associated with drDCM (p < 2.2e-16). At the mRNA level, there is an association between ECHS1 (log2(foldchange) = 1.63329; q = 0.013927) and DBT (log2(foldchange) = 0.955072; q = 0.0368792) with drDCM. ECHS1 is involved in valine (-log (p = 3.39E00)), isoleucine degradation (p = 0.00457), fatty acid ß-oxidation (-log(p) = 2.83E00), and drug metabolism:cytochrome P450 (z-score = 2.07985196) pathways. The mitochondria (-log(p) = 8.73E00), oxidative phosphorylation (-log(p) = 5.35E00) and TCA-cycle II (-log(p) = 2.70E00) are dysfunctional. CONCLUSIONS: We introduce an integrative data strategy that considers the interplay between the DNA, mRNA, and associated pathways, which represents a possible diagnostic, prognostic, biomarker, and personalized treatment discovery approach in genomically heterogeneous diseases.

Whole exome sequencing identifies a troponin T mutation hot spot in familial dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) commonly causes heart failure and shows extensive genetic heterogeneity that may be amenable to newly developed next-generation DNA sequencing of the exome. In this study we report the successful use of exome sequencing to identify a pathogenic variant in the TNNT2 gene using segregation analysis in a large DCM family. Exome sequencing was performed on three distant relatives from a large family with a clear DCM phenotype. Missense, nonsense, and splice variants were analyzed for segregation among the three affected family members and confirmed in other relatives by direct sequencing. A c.517T C>T, Arg173Trp TNNT2 variant segregated with all affected family members and was also detected in one additional DCM family in our registry. The inclusion of segregation analysis using distant family members markedly improved the bioinformatics filtering process by removing from consideration variants that were not shared by all affected subjects. Haplotype analysis confirmed that the variant found in both DCM families was located on two distinct haplotypes, supporting the notion of independent mutational events in each family. In conclusion, an exome sequencing strategy that includes segregation analysis using distant affected relatives within a family represents a viable diagnostic strategy in a genetically heterogeneous disease like DCM.

Genes and environment: novel, functional polymorphism in the human cathepsin L (CTSL1) promoter disrupts a xenobiotic response element (XRE) to alter transcription and blood pressure.

BACKGROUND: Cathepsin L (CTSL1) catalyzes the formation of peptides that influence blood pressure (BP). Naturally occurring genetic variation or targeted ablation of the Ctsl1 locus in mice yield cardiovascular pathology. Here, we searched for genetic variation across the human CTSL1 locus and probed its functional effects, especially in the proximal promoter. METHODS AND RESULTS: Systematic polymorphism discovery by re-sequencing across CTSL1 in 81 patients uncovered 38 genetic variants, five of which were relatively common (MAF >5%), creating a single linkage disequilibrium block in multiple biogeographic ancestries. One of these five common variants lay in a functional domain of the gene: promoter C-171A (rs3118869), which disrupts a predicted xenobiotic response element (XRE; match C>A). In transfected CTSL1 promoter/luciferase reporter plasmids, C-171A allele influenced transcription (C>A, P = 3.36E-6), and transcription was also augmented by co-exposure to the aryl hydrocarbon receptor (AHR) complex (AHR:ARNT) in the presence of their ligand dioxin (P = 6.81E-8); allele (C vs. A) and AHR:ARNT/dioxin stimulus interacted to control gene expression (interaction P = 0.033). Endogenous Ctsl1, Ahr, and Arnt transcripts were present in chromaffin cells. Promoter functional C-171A genotype also predicted hypertension (P = 1.0E-3), SBP (P = 4.0E-4), and DBP (P = 3.0E-3), in an additive pattern for diploid genotypes (A/A > C/A > C/C) in 868 patients, and the results were extended by validation analysis into an independent population sample of 986 patients. CONCLUSION: We conclude that common genetic variation in the proximal CTSL1 promoter, especially at position C-171A, is functional in cells, and alters transcription so as to explain the association of CTSL1 with BP in vivo. At the XRE, endogenous genetic variation plus exogenous aryl hydrocarbon stimulation interact to control CTSL1 gene expression. These results unveil a novel control point whereby heredity and environment can intersect to control a complex trait, and point to new transcriptional strategies for intervention into transmitter biosynthesis and its cardiovascular consequences.

Heritability of nonalcoholic fatty liver disease.

BACKGROUND & AIMS: Nonalcoholic fatty liver disease (NAFLD) is the most common chronic liver disease in the United States. The etiology is believed to be multifactorial with a substantial genetic component; however, the heritability of NAFLD is undetermined. Therefore, a familial aggregation study was performed to test the hypothesis that NAFLD is highly heritable. METHODS: Overweight children with biopsy-proven NAFLD and overweight children without NAFLD served as probands. Family members were studied, including the use of magnetic resonance imaging to quantify liver fat fraction. Fatty liver was defined as a liver fat fraction of 5% or higher. Etiologies for fatty liver other than NAFLD were excluded. Narrow-sense heritability estimates for fatty liver (dichotomous) and fat fraction (continuous) were calculated using variance components analysis adjusted for covariate effects. RESULTS: Fatty liver was present in 17% of siblings and 37% of parents of overweight children without NAFLD. Fatty liver was significantly more common in siblings (59%) and parents (78%) of children with NAFLD. Liver fat fraction was correlated with body mass index, although the correlation was significantly stronger for families of children with NAFLD than those without NAFLD. Adjusted for age, sex, race, and body mass index, the heritability of fatty liver was 1.000 and of liver fat fraction was 0.386. CONCLUSIONS: Family members of children with NAFLD should be considered at high risk for NAFLD. These data suggest that familial factors are a major determinant of whether an individual has NAFLD. Studies examining the complex relations between genes and environment in the development and progression of NAFLD are warranted.