A methylation clock model of mild SARS-CoV-2 infection provides insight into immune dysregulation.

DNA methylation comprises a cumulative record of lifetime exposures superimposed on genetically determined markers. Little is known about methylation dynamics in humans following an acute perturbation, such as infection. We characterized the temporal trajectory of blood epigenetic remodeling in 133 participants in a prospective study of young adults before, during, and after asymptomatic and mildly symptomatic SARS-CoV-2 infection. The differential methylation caused by asymptomatic or mildly symptomatic infections was indistinguishable. While differential gene expression largely returned to baseline levels after the virus became undetectable, some differentially methylated sites persisted for months of follow-up, with a pattern resembling autoimmune or inflammatory disease. We leveraged these responses to construct methylation-based machine learning models that distinguished samples from pre-, during-, and postinfection time periods, and quantitatively predicted the time since infection. The clinical trajectory in the young adults and in a diverse cohort with more severe outcomes was predicted by the similarity of methylation before or early after SARS-CoV-2 infection to the model-defined postinfection state. Unlike the phenomenon of trained immunity, the postacute SARS-CoV-2 epigenetic landscape we identify is antiprotective.

Transcription factor GATA2 may potentiate follicle-stimulating hormone production in mice via induction of the BMP antagonist gremlin in gonadotrope cells.

Mammalian reproduction depends on the gonadotropins, follicle-stimulating hormone (FSH), and luteinizing hormone, which are secreted by pituitary gonadotrope cells. The zinc-finger transcription factor GATA2 was previously implicated in FSH production in male mice; however, its mechanisms of action and role in females were not determined. To directly address GATA2 function in gonadotropes, we generated and analyzed gonadotrope-specific Gata2 KO mice using the Cre-lox system. We found that while conditional KO (cKO) males exhibited ∼50% reductions in serum FSH levels and pituitary FSHß subunit (Fshb) expression relative to controls, FSH production was apparently normal in cKO females. In addition, RNA-seq analysis of purified gonadotropes from control and cKO males revealed a profound decrease in expression of gremlin (Grem1), a bone morphogenetic protein (BMP) antagonist. We show Grem1 was expressed in gonadotropes, but not other cell lineages, in the adult male mouse pituitary. Furthermore, Gata2, Grem1, and Fshb mRNA levels were significantly higher in the pituitaries of WT males relative to females but decreased in males treated with estradiol and increased following ovariectomy in control but not cKO females. Finally, we found that recombinant gremlin stimulated Fshb expression in pituitary cultures from WT mice. Collectively, the data suggest that GATA2 promotes Grem1 expression in gonadotropes and that the gremlin protein potentiates FSH production. The mechanisms of gremlin action have not yet been established but may involve attenuation of BMP binding to activin type II receptors in gonadotropes, facilitating induction of Fshb transcription by activins or related ligands.

SARS-CoV-2 Transmission among Marine Recruits during Quarantine.

BACKGROUND: The efficacy of public health measures to control the transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has not been well studied in young adults. METHODS: We investigated SARS-CoV-2 infections among U.S. Marine Corps recruits who underwent a 2-week quarantine at home followed by a second supervised 2-week quarantine at a closed college campus that involved mask wearing, social distancing, and daily temperature and symptom monitoring. Study volunteers were tested for SARS-CoV-2 by means of quantitative polymerase-chain-reaction (qPCR) assay of nares swab specimens obtained between the time of arrival and the second day of supervised quarantine and on days 7 and 14. Recruits who did not volunteer for the study underwent qPCR testing only on day 14, at the end of the quarantine period. We performed phylogenetic analysis of viral genomes obtained from infected study volunteers to identify clusters and to assess the epidemiologic features of infections. RESULTS: A total of 1848 recruits volunteered to participate in the study; within 2 days after arrival on campus, 16 (0.9%) tested positive for SARS-CoV-2, 15 of whom were asymptomatic. An additional 35 participants (1.9%) tested positive on day 7 or on day 14. Five of the 51 participants (9.8%) who tested positive at any time had symptoms in the week before a positive qPCR test. Of the recruits who declined to participate in the study, 26 (1.7%) of the 1554 recruits with available qPCR results tested positive on day 14. No SARS-CoV-2 infections were identified through clinical qPCR testing performed as a result of daily symptom monitoring. Analysis of 36 SARS-CoV-2 genomes obtained from 32 participants revealed six transmission clusters among 18 participants. Epidemiologic analysis supported multiple local transmission events, including transmission between roommates and among recruits within the same platoon. CONCLUSIONS: Among Marine Corps recruits, approximately 2% who had previously had negative results for SARS-CoV-2 at the beginning of supervised quarantine, and less than 2% of recruits with unknown previous status, tested positive by day 14. Most recruits who tested positive were asymptomatic, and no infections were detected through daily symptom monitoring. Transmission clusters occurred within platoons. (Funded by the Defense Health Agency and others.).

Molecular Signature Predictive of Long-Term Liver Fibrosis Progression to Inform Antifibrotic Drug Development.

BACKGROUND & AIMS: There is a major unmet need to assess the prognostic impact of antifibrotics in clinical trials because of the slow rate of liver fibrosis progression. We aimed to develop a surrogate biomarker to predict future fibrosis progression. METHODS: A fibrosis progression signature (FPS) was defined to predict fibrosis progression within 5 years in patients with hepatitis C virus and nonalcoholic fatty liver disease (NAFLD) with no to minimal fibrosis at baseline (n = 421) and was validated in an independent NAFLD cohort (n = 78). The FPS was used to assess response to 13 candidate antifibrotics in organotypic ex vivo cultures of clinical fibrotic liver tissues (n = 78) and cenicriviroc in patients with nonalcoholic steatohepatitis enrolled in a clinical trial (n = 19, NCT02217475). A serum protein-based surrogate FPS was developed and tested in a cohort of compensated cirrhosis patients (n = 122). RESULTS: A 20-gene FPS was defined and validated in an independent NAFLD cohort (adjusted odds ratio, 10.93; area under the receiver operating characteristic curve, 0.86). Among computationally inferred fibrosis-driving FPS genes, BCL2 was confirmed as a potential pharmacologic target using clinical liver tissues. Systematic ex vivo evaluation of 13 candidate antifibrotics identified rational combination therapies based on epigallocatechin gallate, which were validated for enhanced antifibrotic effect in ex vivo culture of clinical liver tissues. In patients with nonalcoholic steatohepatitis treated with cenicriviroc, FPS modulation was associated with 1-year fibrosis improvement accompanied by suppression of the E2F pathway. Induction of the PPARα pathway was absent in patients without fibrosis improvement, suggesting a benefit of combining PPARα agonism to improve the antifibrotic efficacy of cenicriviroc. A 7-protein serum protein-based surrogate FPS was associated with the development of decompensation in cirrhosis patients. CONCLUSION: The FPS predicts long-term fibrosis progression in an etiology-agnostic manner, which can inform antifibrotic drug development.

Heterogeneous origins and functions of mouse skeletal muscle-resident macrophages.

Tissue-resident macrophages can originate from embryonic or adult hematopoiesis. They play important roles in a wide range of biological processes including tissue remodeling during organogenesis, organ homeostasis, repair following injury, and immune response to pathogens. Although the origins and tissue-specific functions of resident macrophages have been extensively studied in many other tissues, they are not well characterized in skeletal muscle. In the present study, we have characterized the ontogeny of skeletal muscle-resident macrophages by lineage tracing and bone marrow transplant experiments. We demonstrate that skeletal muscle-resident macrophages originate from both embryonic hematopoietic progenitors located within the yolk sac and fetal liver as well as definitive hematopoietic stem cells located within the bone marrow of adult mice. Single-cell-based transcriptome analyses revealed that skeletal muscle-resident macrophages are distinctive from resident macrophages in other tissues as they express a distinct complement of transcription factors and are composed of functionally diverse subsets correlating to their origins. Functionally, skeletal muscle-resident macrophages appear to maintain tissue homeostasis and promote muscle growth and regeneration.

Skeletal muscle transcriptome response to a bout of endurance exercise in physically active and sedentary older adults.

Age-related declines in cardiorespiratory fitness and physical function are mitigated by regular endurance exercise in older adults. This may be due, in part, to changes in the transcriptional program of skeletal muscle following repeated bouts of exercise. However, the impact of chronic exercise training on the transcriptional response to an acute bout of endurance exercise has not been clearly determined. Here, we characterized baseline differences in muscle transcriptome and exercise-induced response in older adults who were active/endurance trained or sedentary. RNA-sequencing was performed on vastus lateralis biopsy specimens obtained before, immediately after, and 3 h following a bout of endurance exercise (40 min of cycling at 60%-70% of heart rate reserve). Using a recently developed bioinformatics approach, we found that transcript signatures related to type I myofibers, mitochondria, and endothelial cells were higher in active/endurance-trained adults and were associated with key phenotypic features including V̇o2peak, ATPmax, and muscle fiber proportion. Immune cell signatures were elevated in the sedentary group and linked to visceral and intermuscular adipose tissue mass. Following acute exercise, we observed distinct temporal transcriptional signatures that were largely similar among groups. Enrichment analysis revealed catabolic processes were uniquely enriched in the sedentary group at the 3-h postexercise timepoint. In summary, this study revealed key transcriptional signatures that distinguished active and sedentary adults, which were associated with difference in oxidative capacity and depot-specific adiposity. The acute response signatures were consistent with beneficial effects of endurance exercise to improve muscle health in older adults irrespective of exercise history and adiposity.NEW & NOTEWORTHY Muscle transcript signatures associated with oxidative capacity and immune cells underlie important phenotypic and clinical characteristics of older adults who are endurance trained or sedentary. Despite divergent phenotypes, the temporal transcriptional signatures in response to an acute bout of endurance exercise were largely similar among groups. These data provide new insight into the transcriptional programs of aging muscle and the beneficial effects of endurance exercise to promote healthy aging in older adults.

SARS-CoV-2 Outbreak Dynamics in an Isolated US Military Recruit Training Center With Rigorous Prevention Measures.

BACKGROUND: Marine recruits training at Parris Island experienced an unexpectedly high rate of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, despite preventive measures including a supervised, 2-week, pre-entry quarantine. We characterize SARS-CoV-2 transmission in this cohort. METHODS: Between May and November 2020, we monitored 2,469 unvaccinated, mostly male, Marine recruits prospectively during basic training. If participants tested negative for SARS-CoV-2 by quantitative polymerase chain reaction (qPCR) at the end of quarantine, they were transferred to the training site in segregated companies and underwent biweekly testing for 6 weeks. We assessed the effects of coronavirus disease 2019 (COVID-19) prevention measures on other respiratory infections with passive surveillance data, performed phylogenetic analysis, and modeled transmission dynamics and testing regimens. RESULTS: Preventive measures were associated with drastically lower rates of other respiratory illnesses. However, among the trainees, 1,107 (44.8%) tested SARS-CoV-2-positive, with either mild or no symptoms. Phylogenetic analysis of viral genomes from 580 participants revealed that all cases but one were linked to five independent introductions, each characterized by accumulation of mutations across and within companies, and similar viral isolates in individuals from the same company. Variation in company transmission rates (mean reproduction number R 0 ; 5.5 [95% confidence interval [CI], 5.0, 6.1]) could be accounted for by multiple initial cases within a company and superspreader events. Simulations indicate that frequent rapid-report testing with case isolation may minimize outbreaks. CONCLUSIONS: Transmission of wild-type SARS-CoV-2 among Marine recruits was approximately twice that seen in the community. Insights from SARS-CoV-2 outbreak dynamics and mutations spread in a remote, congregate setting may inform effective mitigation strategies.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle.

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Innate Immune Response to Influenza Virus at Single-Cell Resolution in Human Epithelial Cells Revealed Paracrine Induction of Interferon Lambda 1.

Early interactions of influenza A virus (IAV) with respiratory epithelium might determine the outcome of infection. The study of global cellular innate immune responses often masks multiple aspects of the mechanisms by which populations of cells work as organized and heterogeneous systems to defeat virus infection, and how the virus counteracts these systems. In this study, we experimentally dissected the dynamics of IAV and human epithelial respiratory cell interaction during early infection at the single-cell level. We found that the number of viruses infecting a cell (multiplicity of infection [MOI]) influences the magnitude of virus antagonism of the host innate antiviral response. Infections performed at high MOIs resulted in increased viral gene expression per cell and stronger antagonist effect than infections at low MOIs. In addition, single-cell patterns of expression of interferons (IFN) and IFN-stimulated genes (ISGs) provided important insights into the contributions of the infected and bystander cells to the innate immune responses during infection. Specifically, the expression of multiple ISGs was lower in infected than in bystander cells. In contrast with other IFNs, IFN lambda 1 (IFNL1) showed a widespread pattern of expression, suggesting a different cell-to-cell propagation mechanism more reliant on paracrine signaling. Finally, we measured the dynamics of the antiviral response in primary human epithelial cells, which highlighted the importance of early innate immune responses at inhibiting virus spread.IMPORTANCE Influenza A virus (IAV) is a respiratory pathogen of high importance to public health. Annual epidemics of seasonal IAV infections in humans are a significant public health and economic burden. IAV also causes sporadic pandemics, which can have devastating effects. The main target cells for IAV replication are epithelial cells in the respiratory epithelium. The cellular innate immune responses induced in these cells upon infection are critical for defense against the virus, and therefore, it is important to understand the complex interactions between the virus and the host cells. In this study, we investigated the innate immune response to IAV in the respiratory epithelium at the single-cell level, providing a better understanding on how a population of epithelial cells functions as a complex system to orchestrate the response to virus infection and how the virus counteracts this system.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Combination of Gene Expression Signature and Model for End-Stage Liver Disease Score Predicts Survival of Patients With Severe Alcoholic Hepatitis.

BACKGROUND & AIMS: Patients with severe alcoholic hepatitis (AH) have a high risk of death within 90 days. Corticosteroids, which can cause severe adverse events, are the only treatment that increases short-term survival. It is a challenge to predict outcomes of patients with severe AH. Therefore, we developed a scoring system to predict patient survival, integrating baseline molecular and clinical variables. METHODS: We obtained fixed liver biopsy samples from 71 consecutive patients diagnosed with severe AH and treated with corticosteroids from July 2006 through December 2013 in Brussels, Belgium (derivation cohort). Gene expression patterns were analyzed by microarrays and clinical data were collected for 180 days. We identified gene expression signatures and clinical data that are associated with survival without liver transplantation at 90 and 180 days after initiation of corticosteroid therapy. Findings were validated using liver biopsies from 48 consecutive patients with severe AH treated with corticosteroids, collected from March 2010 through February 2015 at hospitals in Belgium and Switzerland (validation cohort 1) and in liver biopsies from 20 patients (9 received corticosteroid treatment), collected from January 2012 through May 2015 in the United States (validation cohort 2). RESULTS: We integrated data on expression patterns of 123 genes and the model for end-stage liver disease (MELD) scores to assign patients to groups with poor survival (29% survived 90 days and 26% survived 180 days) and good survival (76% survived 90 days and 65% survived 180 days) (P < .001) in the derivation cohort. We named this assignment system the gene signature-MELD (gs-MELD) score. In validation cohort 1, the gs-MELD score discriminated patients with poor survival (43% survived 90 days) from those with good survival (96% survived 90 days) (P < .001). The gs-MELD score also discriminated between patients with a poor survival at 180 days (34% survived) and a good survival at 180 days (84% survived) (P < .001). The time-dependent area under the receiver operator characteristic curve for the score was 0.86 (95% confidence interval 0.73-0.99) for survival at 90 days, and 0.83 (95% confidence interval 0.71-0.96) for survival at 180 days. This score outperformed other clinical models to predict survival of patients with severe AH in validation cohort 1. In validation cohort 2, the gs-MELD discriminated patients with a poor survival at 90 days (12% survived) from those with a good survival at 90 days (100%) (P < .001). CONCLUSIONS: We integrated data on baseline liver gene expression pattern and the MELD score to create the gs-MELD scoring system, which identifies patients with severe AH, treated or not with corticosteroids, most and least likely to survive for 90 and 180 days.

Pharmacokinetics and relative bioavailability of tiamulin in broiler chicken as influenced by different routes of administration.

The bioavailability and pharmacokinetic disposition of tiamulin in broiler chicken were investigated after administration through the crop, drinking water, and feed at 40 mg/kg body weight. Residues of tiamulin in tissues of broiler chicken were also assessed. Plasma and tissue concentrations of tiamulin were analyzed by reverse-phase high-performance liquid chromatography (HPLC) method. Plasma concentration-time data were described by the non-compartmental model for all three routes, and pharmacokinetic parameters were calculated. There were no significant differences (p > 0.05) in pharmacokinetic parameters and mean plasma concentrations of tiamulin between three routes tested (crop, water, and feed), indicating equal efficacy. Tiamulin residues in edible tissues (muscles, skin, and fat) were lower than the advocated maximum residue limit (MRL of 0.1 µg/g and that of liver was 1 µg/g) on the 3rd day. No traces were found on the 5th day after drug administration. This indicated that the withdrawal period (less than 5 days) is very short, which makes it safer. This study shows that tiamulin can be used with equal efficacy through all routes of administration in broiler chicken (crop, water, and feed).

Clinicopathological indices to predict hepatocellular carcinoma molecular classification.

BACKGROUND & AIMS: Hepatocellular carcinoma (HCC) is the second most lethal cancer caused by lack of effective therapies. Although promising, HCC molecular classification, which enriches potential responders to specific therapies, has not yet been assessed in clinical trials of anti-HCC drugs. We aimed to overcome these challenges by developing clinicopathological surrogate indices of HCC molecular classification. METHODS: Hepatocellular carcinoma classification defined in our previous transcriptome meta-analysis (S1, S2 and S3 subclasses) was implemented in an FDA-approved diagnostic platform (Elements assay, NanoString). Ninety-six HCC tumours (training set) were assayed to develop molecular subclass-predictive indices based on clinicopathological features, which were independently validated in 99 HCC tumours (validation set). Molecular deregulations associated with the histopathological features were determined by pathway analysis. Sample sizes for HCC clinical trials enriched with specific molecular subclasses were determined. RESULTS: Hepatocellular carcinoma subclass-predictive indices were steatohepatitic (SH)-HCC variant and immune cell infiltrate for S1 subclass, macrotrabecular/compact pattern, lack of pseudoglandular pattern, and high serum alpha-foetoprotein (>400 ng/ml) for S2 subclass, and microtrabecular pattern, lack of SH-HCC and clear cell variants, and lower histological grade for S3 subclass. Macrotrabecular/compact pattern, a predictor of S2 subclass, was associated with the activation of therapeutically targetable oncogene YAP and stemness markers EPCAM/KRT19. BMP4 was associated with pseudoglandular pattern. Subclass-predictive indices-based patient enrichment reduced clinical trial sample sizes from 121, 184 and 53 to 30, 43 and 22 for S1, S2 and S3 subclass-targeting therapies respectively. CONCLUSIONS: Hepatocellular carcinoma molecular subclasses can be enriched by clinicopathological indices tightly associated with deregulation of therapeutically targetable molecular pathways.

A genomic and clinical prognostic index for hepatitis C-related early-stage cirrhosis that predicts clinical deterioration.

OBJECTIVE: The number of patients with HCV-related cirrhosis is increasing, leading to a rising risk of complications and death. Prognostic stratification in patients with early-stage cirrhosis is still challenging. We aimed to develop and validate a clinically useful prognostic index based on genomic and clinical variables to identify patients at high risk of disease progression. DESIGN: We developed a prognostic index, comprised of a 186-gene signature validated in our previous genome-wide profiling study, bilirubin (>1 mg/dL) and platelet count (<100,000/mm(3)), in an Italian HCV cirrhosis cohort (training cohort, n=216, median follow-up 10 years). The gene signature test was implemented using a digital transcript counting (nCounter) assay specifically developed for clinical use and the prognostic index was evaluated using archived specimens from an independent cohort of HCV-related cirrhosis in the USA (validation cohort, n=145, median follow-up 8 years). RESULTS: In the training cohort, the prognostic index was associated with hepatic decompensation (HR=2.71, p=0.003), overall death (HR=6.00, p<0.001), hepatocellular carcinoma (HR=3.31, p=0.001) and progression of Child-Turcotte-Pugh class (HR=6.70, p<0.001). The patients in the validation cohort were stratified into high-risk (16%), intermediate-risk (42%) or low-risk (42%) groups by the prognostic index. The high-risk group had a significantly increased risk of hepatic decompensation (HR=7.36, p<0.001), overall death (HR=3.57, p=0.002), liver-related death (HR=6.49, p<0.001) and all liver-related adverse events (HR=4.98, p<0.001). CONCLUSIONS: A genomic and clinical prognostic index readily available for clinical use was successfully validated, warranting further clinical evaluation for prognostic prediction and clinical trial stratification and enrichment for preventive interventions.

Low-variance RNAs identify Parkinson's disease molecular signature in blood.

The diagnosis of Parkinson's disease (PD) is usually not established until advanced neurodegeneration leads to clinically detectable symptoms. Previous blood PD transcriptome studies show low concordance, possibly resulting from the use of microarray technology, which has high measurement variation. The Leucine-rich repeat kinase 2 (LRRK2) G2019S mutation predisposes to PD. Using preclinical and clinical studies, we sought to develop a novel statistically motivated transcriptomic-based approach to identify a molecular signature in the blood of Ashkenazi Jewish PD patients, including LRRK2 mutation carriers. Using a digital gene expression platform to quantify 175 messenger RNA (mRNA) markers with low coefficients of variation (CV), we first compared whole-blood transcript levels in mouse models (1) overexpressing wild-type (WT) LRRK2, (2) overexpressing G2019S LRRK2, (3) lacking LRRK2 (knockout), and (4) and in WT controls. We then studied an Ashkenazi Jewish cohort of 34 symptomatic PD patients (both WT LRRK2 and G2019S LRRK2) and 32 asymptomatic controls. The expression profiles distinguished the four mouse groups with different genetic background. In patients, we detected significant differences in blood transcript levels both between individuals differing in LRRK2 genotype and between PD patients and controls. Discriminatory PD markers included genes associated with innate and adaptive immunity and inflammatory disease. Notably, gene expression patterns in levodopa-treated PD patients were significantly closer to those of healthy controls in a dose-dependent manner. We identify whole-blood mRNA signatures correlating with LRRK2 genotype and with PD disease state. This approach may provide insight into pathogenesis and a route to early disease detection.

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.

Localization of BRCA1 protein in breast cancer tissue and cell lines with mutations.

BACKGROUND: The breast and ovarian cancer susceptibility gene (BRCA1) encodes a tumor suppressor. The BRCA1 protein is found primarily in cell nuclei and plays an important role in the DNA damage response and transcriptional regulation. Deficiencies in DNA repair capabilities have been associated with higher histopathological grade and worse prognosis in breast cancer. METHODS: In order to investigate the subcellular distribution of BRCA1 in tumor tissue we randomly selected 22 breast carcinomas and tested BRCA1 protein localization in frozen and contiguous formalin-fixed, paraffin embedded (FFPE) tissue, using pressure cooker antigen-retrieval and the MS110 antibody staining. To assess the impact of BRCA1 germline mutations on protein localization, we retrospectively tested 16 of the tumor specimens to determine whether they contained the common Ashkenazi Jewish founder mutations in BRCA1 (185delAG, 5382insC), and BRCA2 (6174delT). We also compared co-localization of BRCA1 and nucleolin in MCF7 cells (wild type) and a mutant BRCA1 cell line, HCC1937 (5382insC). RESULTS: In FFPE tissue, with MS110 antibody staining, we frequently found reduced BRCA1 nuclear staining in breast tumor tissue compared to normal tissue, and less BRCA1 staining with higher histological grade in the tumors. However, in the frozen sections, BRCA1 antibody staining showed punctate, intra-nuclear granules in varying numbers of tumor, lactating, and normal cells. Two mutation carriers were identified and were confirmed by gene sequencing. We have also compared co-localization of BRCA1 and nucleolin in MCF7 cells (wild type) and a mutant BRCA1 cell line, HCC1937 (5382insC) and found altered sub-nuclear and nucleolar localization patterns consistent with a functional impact of the mutation on protein localization. CONCLUSIONS: The data presented here support a role for BRCA1 in the pathogenesis of sporadic and inherited breast cancers. The use of well-characterized reagents may lead to further insights into the function of BRCA1 and possibly the further development of targeted therapeutics.

Activating Transcription Factor 3 Stimulates Follicle-Stimulating Hormone-ß Expression In Vitro But Is Dispensable for Follicle-Stimulating Hormone Production in Murine Gonadotropes In Vivo.

Follicle-stimulating hormone (FSH), a dimeric glycoprotein produced by pituitary gonadotrope cells, regulates spermatogenesis in males and ovarian follicle growth in females. Hypothalamic gonadotropin-releasing hormone (GnRH) stimulates FSHß subunit gene (Fshb) transcription, though the underlying mechanisms are poorly understood. To address this gap in knowledge, we examined changes in pituitary gene expression in GnRH-deficient mice (hpg) treated with a regimen of exogenous GnRH that increases pituitary Fshb but not luteinizing hormone ß (Lhb) messenger RNA levels. Activating transcription factor 3 (Atf3) was among the most upregulated genes. Activating transcription factor 3 (ATF3) can heterodimerize with members of the activator protein 1 family to regulate gene transcription. Co-expression of ATF3 with JunB stimulated murine Fshb, but not Lhb, promoter-reporter activity in homologous LßT2b cells. ATF3 also synergized with a constitutively active activin type I receptor to increase endogenous Fshb expression in these cells. Nevertheless, FSH production was intact in gonadotrope-specific Atf3 knockout [conditional knockout (cKO)] mice. Ovarian follicle development, ovulation, and litter sizes were equivalent between cKOs and controls. Testis weights and sperm counts did not differ between genotypes. Following gonadectomy, increases in LH secretion were enhanced in cKO animals. Though FSH levels did not differ between genotypes, post-gonadectomy increases in pituitary Fshb and gonadotropin α subunit expression were more pronounced in cKO than control mice. These data indicate that ATF3 can selectively stimulate Fshb expression in vitro but is not required for FSH production in vivo.

Integrated single-cell multiome analysis reveals muscle fiber-type gene regulatory circuitry modulated by endurance exercise.

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.