Protocol for high-throughput DNA methylation profiling in rat tissues using automated reduced representation bisulfite sequencing.

Although reduced representation bisulfite sequencing (RRBS) measures DNA methylation (DNAme) across CpG-rich genomic regions with high sensitivity, the assay can be time-consuming and prone to batch effects. Here, we present a high-throughput, automated RRBS protocol starting with DNA extraction from frozen rat tissues. We describe steps for RRBS library preparation, library quality control, and sequencing. We also detail an optimized pipeline for sequencing data processing. This protocol has been applied successfully to DNAme profiling across multiple rat tissues. For complete details on the use and execution of this protocol, please refer to Nair et al.1.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

The mitochondrial multi-omic response to exercise training across rat tissues.

Mitochondria have diverse functions critical to whole-body metabolic homeostasis. Endurance training alters mitochondrial activity, but systematic characterization of these adaptations is lacking. Here, the Molecular Transducers of Physical Activity Consortium mapped the temporal, multi-omic changes in mitochondrial analytes across 19 tissues in male and female rats trained for 1, 2, 4, or 8 weeks. Training elicited substantial changes in the adrenal gland, brown adipose, colon, heart, and skeletal muscle. The colon showed non-linear response dynamics, whereas mitochondrial pathways were downregulated in brown adipose and adrenal tissues. Protein acetylation increased in the liver, with a shift in lipid metabolism, whereas oxidative proteins increased in striated muscles. Exercise-upregulated networks were downregulated in human diabetes and cirrhosis. Knockdown of the central network protein 17-beta-hydroxysteroid dehydrogenase 10 (HSD17B10) elevated oxygen consumption, indicative of metabolic stress. We provide a multi-omic, multi-tissue, temporal atlas of the mitochondrial response to exercise training and identify candidates linked to mitochondrial dysfunction.

Class IIa HDAC4 and HDAC7 cooperatively regulate gene transcription in Th17 cell differentiation.

Class II histone deacetylases (HDACs) are important in regulation of gene transcription during T cell development. However, our understanding of their cell-specific functions is limited. In this study, we reveal that class IIa Hdac4 and Hdac7 (Hdac4/7) are selectively induced in transcription, guiding the lineage-specific differentiation of mouse T-helper 17 (Th17) cells from naive CD4+ T cells. Importantly, Hdac4/7 are functionally dispensable in other Th subtypes. Mechanistically, Hdac4 interacts with the transcription factor (TF) JunB, facilitating the transcriptional activation of Th17 signature genes such as Il17a/f. Conversely, Hdac7 collaborates with the TF Aiolos and Smrt/Ncor1-Hdac3 corepressors to repress transcription of Th17 negative regulators, including Il2, in Th17 cell differentiation. Inhibiting Hdac4/7 through pharmacological or genetic methods effectively mitigates Th17 cell-mediated intestinal inflammation in a colitis mouse model. Our study uncovers molecular mechanisms where HDAC4 and HDAC7 function distinctively yet cooperatively in regulating ordered gene transcription during Th17 cell differentiation. These findings suggest a potential therapeutic strategy of targeting HDAC4/7 for treating Th17-related inflammatory diseases, such as ulcerative colitis.

Comprehensive interrogation of human skeletal muscle reveals a dissociation between insulin resistance and mitochondrial capacity.

Insulin resistance and blunted mitochondrial capacity in skeletal muscle are often synonymous, however, this association remains controversial. The aim of this study was to perform an in-depth multifactorial comparison of skeletal muscle mitochondrial capacity between individuals who were lean and active (Active, n = 9), individuals with obesity (Obese, n = 9), and individuals with obesity, insulin resistance, and type 2 diabetes (T2D, n = 22). Mitochondrial capacity was assessed by ex vivo mitochondrial respiration with fatty-acid and glycolytic-supported protocols adjusted for mitochondrial content (mtDNA and citrate synthase activity). Supercomplex assembly was measured by Blue Native (BN)-PAGE and immunoblot. Tricarboxylic (TCA) cycle intermediates were assessed with targeted metabolomics. Exploratory transcriptomics and DNA methylation analyses were performed to uncover molecular differences affecting mitochondrial function among the three groups. We reveal no discernable differences in skeletal muscle mitochondrial content, mitochondrial capacity, supercomplex assembly, TCA cycle intermediates, and mitochondrial molecular profiles between obese individuals with and without T2D that had comparable levels of confounding factors (body mass index, age, and aerobic capacity). We highlight that lean, active individuals have greater mitochondrial content, mitochondrial capacity, supercomplex assembly, and TCA cycle intermediates. These phenotypical changes are reflected at the level of DNA methylation and gene transcription. The collective observation of comparable muscle mitochondrial capacity in individuals with obesity and T2D (vs. individuals without T2D) underscores a dissociation from skeletal muscle insulin resistance. Clinical trial number: NCT01911104.NEW & NOTEWORTHY Whether impaired mitochondrial capacity contributes to skeletal muscle insulin resistance is debated. Our multifactorial analysis shows no differences in skeletal muscle mitochondrial content, mitochondrial capacity, and mitochondrial molecular profiles between obese individuals with and without T2D that had comparable levels of confounding factors (BMI, age, aerobic capacity). We highlight that lean, active individuals have enhanced skeletal muscle mitochondrial capacity that is also reflected at the level of DNA methylation and gene transcription.

Development of a Nickel-Catalyzed N-N Coupling for the Synthesis of Hydrazides.

A nickel-catalyzed N-N cross-coupling for the synthesis of hydrazides is reported. O-Benzoylated hydroxamates were efficiently coupled with a broad range of aryl and aliphatic amines via nickel catalysis to form hydrazides in an up to 81% yield. Experimental evidence implicates the intermediacy of electrophilic Ni-stabilized acyl nitrenoids and the formation of a Ni(I) catalyst via silane-mediated reduction. This report constitutes the first example of an intermolecular N-N coupling compatible with secondary aliphatic amines.

Dysfunction of ubiquitin protein ligase MYCBP2 leads to cell resilience in human breast cancers.

Breast cancer is the most common type of cancer among women worldwide, and it is estimated that 294 000 new diagnoses and 37 000 deaths will occur each year in the United States alone by 2030. Large-scale genomic studies have identified a number of genetic loci with alterations in breast cancer. However, identification of the genes that are critical for tumorgenicity still remains a challenge. Here, we perform a comprehensive functional multi-omics analysis of somatic mutations in breast cancer and identify previously unknown key regulators of breast cancer tumorgenicity. We identify dysregulation of MYCBP2, an E3 ubiquitin ligase and an upstream regulator of mTOR signaling, is accompanied with decreased disease-free survival. We validate MYCBP2 as a key target through depletion siRNA using in vitro apoptosis assays in MCF10A, MCF7 and T47D cells. We demonstrate that MYCBP2 loss is associated with resistance to apoptosis from cisplatin-induced DNA damage and cell cycle changes, and that CHEK1 inhibition can modulate MYCBP2 activity and caspase cleavage. Furthermore, we show that MYCBP2 knockdown is associated with transcriptomic responses in TSC2 and in apoptosis genes and interleukins. Therefore, we show that MYCBP2 is an important genetic target that represents a key node regulating multiple molecular pathways in breast cancer corresponding with apparent drug resistance in our study.

Histone H3 lysine 27 crotonylation mediates gene transcriptional repression in chromatin.

Histone lysine acylation, including acetylation and crotonylation, plays a pivotal role in gene transcription in health and diseases. However, our understanding of histone lysine acylation has been limited to gene transcriptional activation. Here, we report that histone H3 lysine 27 crotonylation (H3K27cr) directs gene transcriptional repression rather than activation. Specifically, H3K27cr in chromatin is selectively recognized by the YEATS domain of GAS41 in complex with SIN3A-HDAC1 co-repressors. Proto-oncogenic transcription factor MYC recruits GAS41/SIN3A-HDAC1 complex to repress genes in chromatin, including cell-cycle inhibitor p21. GAS41 knockout or H3K27cr-binding depletion results in p21 de-repression, cell-cycle arrest, and tumor growth inhibition in mice, explaining a causal relationship between GAS41 and MYC gene amplification and p21 downregulation in colorectal cancer. Our study suggests that H3K27 crotonylation signifies a previously unrecognized, distinct chromatin state for gene transcriptional repression in contrast to H3K27 trimethylation for transcriptional silencing and H3K27 acetylation for transcriptional activation.

Dissecting aortic aneurysm in Marfan syndrome is associated with losartan-sensitive transcriptomic modulation of aortic cells.

To improve our limited understanding of the pathogenesis of thoracic aortic aneurysm (TAA) that leads to acute aortic dissection, single-cell RNA sequencing (scRNA-seq) was employed to profile disease-relevant transcriptomic changes of aortic cell populations in a well-characterized mouse model of the most commonly diagnosed form of Marfan syndrome (MFS). As result, 2 discrete subpopulations of aortic cells (SMC3 and EC4) were identified only in the aorta of Fbn1mgR/mgR mice. SMC3 cells highly express genes related to extracellular matrix formation and nitric oxide signaling, whereas the EC4 transcriptional profile is enriched in smooth muscle cell (SMC), fibroblast, and immune cell-related genes. Trajectory analysis predicted close phenotypic modulation between SMC3 and EC4, which were therefore analyzed together as a discrete MFS-modulated (MFSmod) subpopulation. In situ hybridization of diagnostic transcripts located MFSmod cells at the intima of Fbn1mgR/mgR aortas. Reference-based data set integration revealed transcriptomic similarity between MFSmod- and SMC-derived cell clusters modulated in human TAA. Consistent with the angiotensin II type I receptor (At1r) contribution to TAA development, MFSmod cells were absent in the aorta of Fbn1mgR/mgR mice treated with the At1r antagonist losartan. Altogether, our findings indicate that a discrete dynamic alteration of aortic cell identity is associated with dissecting TAA in MFS mice and increased risk of aortic dissection in MFS patients.

Isolation of nuclei from frozen human subcutaneous adipose tissue for full-length single-nuclei transcriptional profiling.

Automated single-cell dispensing is incompatible with white adipose tissue (WAT) due to lipid-laden adipocytes. Single-nuclei RNA-Seq permits transcriptional profiling of all cells from WAT. Human WAT faces unique technical challenges in isolating nuclei compared to rodent tissue due to greater extra-cellular matrix content and larger lipid droplets. In this protocol, we detail how to isolate nuclei from frozen subcutaneous human WAT for single-nuclei RNA-Seq. For complete information on the generation and use of this protocol, please refer to Whytock et al. (2022).1.

Multi-omics-based analysis of high grade serous ovarian cancer subtypes reveals distinct molecular processes linked to patient prognosis.

Despite advancements in treatment, high-grade serous ovarian cancer (HGSOC) is still characterized by poor patient outcomes. To understand the molecular heterogeneity of this disease, which underlies the challenge in selecting optimal treatments for HGSOC patients, we have integrated genomic, transcriptomic, and epigenetic information to identify seven new HGSOC subtypes using a multiscale clustering method. These subtypes not only have significantly distinct overall survival, but also exhibit unique patterns of gene expression, microRNA expression, DNA methylation, and copy number alterations. As determined by our analysis, patients with similar clinical outcomes have distinct profiles of activated or repressed cellular processes, including cell cycle, epithelial-to-mesenchymal transition, immune activation, interferon response, and cilium organization. Furthermore, we performed a multiscale gene co-expression network analysis to identify subtype-specific key regulators and predicted optimal targeted therapies based on subtype-specific gene expression. In summary, this study provides new insights into the cellular heterogeneity of the HGSOC genomic, epigenetic, and transcriptomic landscapes and provides a basis for future studies into precision medicine for HGSOC patients.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training.

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

The mitochondrial multi-omic response to exercise training across tissues.

Mitochondria are adaptable organelles with diverse cellular functions critical to whole-body metabolic homeostasis. While chronic endurance exercise training is known to alter mitochondrial activity, these adaptations have not yet been systematically characterized. Here, the Molecular Transducers of Physical Activity Consortium (MoTrPAC) mapped the longitudinal, multi-omic changes in mitochondrial analytes across 19 tissues in male and female rats endurance trained for 1, 2, 4 or 8 weeks. Training elicited substantial changes in the adrenal gland, brown adipose, colon, heart and skeletal muscle, while we detected mild responses in the brain, lung, small intestine and testes. The colon response was characterized by non-linear dynamics that resulted in upregulation of mitochondrial function that was more prominent in females. Brown adipose and adrenal tissues were characterized by substantial downregulation of mitochondrial pathways. Training induced a previously unrecognized robust upregulation of mitochondrial protein abundance and acetylation in the liver, and a concomitant shift in lipid metabolism. The striated muscles demonstrated a highly coordinated response to increase oxidative capacity, with the majority of changes occurring in protein abundance and post-translational modifications. We identified exercise upregulated networks that are downregulated in human type 2 diabetes and liver cirrhosis. In both cases HSD17B10, a central dehydrogenase in multiple metabolic pathways and mitochondrial tRNA maturation, was the main hub. In summary, we provide a multi-omic, cross-tissue atlas of the mitochondrial response to training and identify candidates for prevention of disease-associated mitochondrial dysfunction.

BRD4-PRC2 represses transcription of T-helper 2-specific negative regulators during T-cell differentiation.

BRD4 is a well-recognized transcriptional activator, but how it regulates gene transcriptional repression in a cell type-specific manner has remained elusive. In this study, we report that BRD4 works with Polycomb repressive complex 2 (PRC2) to repress transcriptional expression of the T-helper 2 (Th2)-negative regulators Foxp3 and E3-ubiqutin ligase Fbxw7 during lineage-specific differentiation of Th2 cells from mouse primary naïve CD4+ T cells. Brd4 binds to the lysine-acetylated-EED subunit of the PRC2 complex via its second bromodomain (BD2) to facilitate histone H3 lysine 27 trimethylation (H3K27me3) at target gene loci and thereby transcriptional repression. We found that Foxp3 represses transcription of Th2-specific transcription factor Gata3, while Fbxw7 promotes its ubiquitination-directed protein degradation. BRD4-mediated repression of Foxp3 and Fbxw7 in turn promotes BRD4- and Gata3-mediated transcriptional activation of Th2 cytokines including Il4, Il5, and Il13. Chemical inhibition of the BRD4 BD2 induces transcriptional de-repression of Foxp3 and Fbxw7, and thus transcriptional downregulation of Il4, Il5, and Il13, resulting in inhibition of Th2 cell lineage differentiation. Our study presents a previously unappreciated mechanism of BRD4's role in orchestrating a Th2-specific transcriptional program that coordinates gene repression and activation, and safeguards cell lineage differentiation.

Single-Cell RNA Sequencing Reveals New Basic and Translational Insights in the Cystic Fibrosis Lung.

Cystic fibrosis (CF) is a multisystemic, autosomal recessive disorder caused by mutations in the CFTR (cystic fibrosis transmembrane conductance regulator) gene, with the majority of morbidity and mortality extending from lung disease. Single-cell RNA sequencing (scRNA-seq) has been leveraged in the lung and elsewhere in the body to articulate discrete cell populations, describing cell types, states, and lineages as well as their roles in health and disease. In this translational review, we provide an overview of the current applications of scRNA-seq to the study of the normal and CF lungs, allowing the beginning of a new cellular and molecular narrative of CF lung disease, and we highlight some of the future opportunities to further leverage scRNA-seq and complementary single-cell technologies in the study of CF as we bridge from scientific understanding to clinical application.

Selective Inhibitors of Autophagy Reveal New Link between the Cell Cycle and Autophagy and Lead to Discovery of Novel Synergistic Drug Combinations.

Autophagy is a conserved metabolic pathway that is central to many diseases. Recently, there has been a lot of interest in targeting autophagy with small molecule inhibitors as a possible therapeutic strategy. However, many of the compounds used for autophagy are nonselective. Here, we explored the inhibition of autophagy in pancreatic cancer cells using established selective small molecule inhibitors and discovered an unexpected link between the autophagy pathway and progression through the cell cycle. Our findings revealed that treatments with inhibitors that have different autophagy pathway targets block cell replication and activate other metabolic pathways to compensate for the blockade in autophagy. An unbiased screen looking for known drugs that might synergize with autophagy inhibition revealed new combination treatments that might provide a blueprint for therapeutic approaches to pancreatic cancer. The drugs quizartinib and THZ1 showed a strong synergistic effect in pancreatic cells with autophagy inhibition.

Single cell full-length transcriptome of human subcutaneous adipose tissue reveals unique and heterogeneous cell populations.

White adipose tissue (WAT) is a complex mixture of adipocytes and non-adipogenic cells. Characterizing the cellular composition of WAT is critical for identifying where potential alterations occur that impact metabolism. Most single-cell (sc) RNA-Seq studies focused on the stromal vascular fraction (SVF) which does not contain adipocytes and have used technology that has a 3' or 5' bias. Using full-length sc/single-nuclei (sn) RNA-Seq technology, we interrogated the transcriptional composition of WAT using: snRNA-Seq of whole tissue, snRNA-Seq of isolated adipocytes, and scRNA-Seq of SVF. Whole WAT snRNA-Seq provided coverage of major cell types, identified three distinct adipocyte clusters, and was capable of tracking adipocyte differentiation with pseudotime. Compared to WAT, adipocyte snRNA-Seq was unable to match adipocyte heterogeneity. SVF scRNA-Seq provided greater resolution of non-adipogenic cells. These findings provide critical evidence for the utility of sc full-length transcriptomics in WAT and SVF in humans.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer.

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Sedentary and Trained Older Men Have Distinct Circulating Exosomal microRNA Profiles at Baseline and in Response to Acute Exercise.

Exercise has multi-systemic benefits and attenuates the physiological impairments associated with aging. Emerging evidence suggests that circulating exosomes mediate some of the beneficial effects of exercise via the transfer of microRNAs between tissues. However, the impact of regular exercise and acute exercise on circulating exosomal microRNAs (exomiRs) in older populations remains unknown. In the present study, we analyzed circulating exomiR expression in endurance-trained elderly men (n = 5) and age-matched sedentary males (n = 5) at baseline (Pre), immediately after a forty minute bout of aerobic exercise on a cycle ergometer (Post), and three hours after this acute exercise (3hPost). Following the isolation and enrichment of exosomes from plasma, exosome-enriched preparations were characterized and exomiR levels were determined by sequencing. The effect of regular exercise on circulating exomiRs was assessed by comparing the baseline expression levels in the trained and sedentary groups. The effect of acute exercise was determined by comparing baseline and post-training expression levels in each group. Regular exercise resulted in significantly increased baseline expression of three exomiRs (miR-486-5p, miR-215-5p, miR-941) and decreased expression of one exomiR (miR-151b). Acute exercise altered circulating exomiR expression in both groups. However, exomiRs regulated by acute exercise in the trained group (7 miRNAs at Post and 8 at 3hPost) were distinct from those in the sedentary group (9 at Post and 4 at 3hPost). Pathway analysis prediction and reported target validation experiments revealed that the majority of exercise-regulated exomiRs are targeting genes that are related to IGF-1 signaling, a pathway involved in exercise-induced muscle and cardiac hypertrophy. The immediately post-acute exercise exomiR signature in the trained group correlates with activation of IGF-1 signaling, whereas in the sedentary group it is associated with inhibition of IGF-1 signaling. While further validation is needed, including measurements of IGF-1/IGF-1 signaling in blood or skeletal muscle, our results suggest that training status may counteract age-related anabolic resistance by modulating circulating exomiR profiles both at baseline and in response to acute exercise.

Molecular Transducers of Physical Activity Consortium (MoTrPAC): Mapping the Dynamic Responses to Exercise.

Exercise provides a robust physiological stimulus that evokes cross-talk among multiple tissues that when repeated regularly (i.e., training) improves physiological capacity, benefits numerous organ systems, and decreases the risk for premature mortality. However, a gap remains in identifying the detailed molecular signals induced by exercise that benefits health and prevents disease. The Molecular Transducers of Physical Activity Consortium (MoTrPAC) was established to address this gap and generate a molecular map of exercise. Preclinical and clinical studies will examine the systemic effects of endurance and resistance exercise across a range of ages and fitness levels by molecular probing of multiple tissues before and after acute and chronic exercise. From this multi-omic and bioinformatic analysis, a molecular map of exercise will be established. Altogether, MoTrPAC will provide a public database that is expected to enhance our understanding of the health benefits of exercise and to provide insight into how physical activity mitigates disease.