Vascular dysfunction is a hallmark of chronic diseases in elderly. The contribution of the vasculature to lung repair and fibrosis is not fully understood. Here, we performed an epigenetic and ...transcriptional analysis of lung endothelial cells (ECs) from young and aged mice during the resolution or progression of bleomycin-induced lung fibrosis. We identified the transcription factor ETS-related gene (ERG) as putative orchestrator of lung capillary homeostasis and repair, and whose function is dysregulated in aging. ERG dysregulation is associated with reduced chromatin accessibility and maladaptive transcriptional responses to injury. Loss of endothelial ERG enhances paracrine fibroblast activation in vitro, and impairs lung fibrosis resolution in young mice in vivo. scRNA-seq of ERG deficient mouse lungs reveales transcriptional and fibrogenic abnormalities resembling those associated with aging and human lung fibrosis, including reduced number of general capillary (gCap) ECs. Our findings demonstrate that lung endothelial chromatin remodeling deteriorates with aging leading to abnormal transcription, vascular dysrepair, and persistent fibrosis following injury.
The mammalian circadian clock relies on the transcription factor CLOCK:BMAL1 to coordinate the rhythmic expression of thousands of genes. Consistent with the various biological functions under clock ...control, rhythmic gene expression is tissue-specific despite an identical clockwork mechanism in every cell. Here we show that BMAL1 DNA binding is largely tissue-specific, likely because of differences in chromatin accessibility between tissues and cobinding of tissue-specific transcription factors. Our results also indicate that BMAL1 ability to drive tissue-specific rhythmic transcription is associated with not only the activity of BMAL1-bound enhancers but also the activity of neighboring enhancers. Characterization of physical interactions between BMAL1 enhancers and other
-regulatory regions by RNA polymerase II chromatin interaction analysis by paired-end tag (ChIA-PET) reveals that rhythmic BMAL1 target gene expression correlates with rhythmic chromatin interactions. These data thus support that much of BMAL1 target gene transcription depends on BMAL1 capacity to rhythmically regulate a network of enhancers.
Background
COVID-19 is caused by the SARS-CoV-2 virus and has strikingly heterogeneous clinical manifestations, with most individuals contracting mild disease but a substantial minority experiencing ...fulminant cardiopulmonary symptoms or death. The clinical covariates and the laboratory tests performed on a patient provide robust statistics to guide clinical treatment. Deep learning approaches on a data set of this nature enable patient stratification and provide methods to guide clinical treatment.
Objective
Here, we report on the development and prospective validation of a state-of-the-art machine learning model to provide mortality prediction shortly after confirmation of SARS-CoV-2 infection in the Mayo Clinic patient population.
Methods
We retrospectively constructed one of the largest reported and most geographically diverse laboratory information system and electronic health record of COVID-19 data sets in the published literature, which included 11,807 patients residing in 41 states of the United States of America and treated at medical sites across 5 states in 3 time zones. Traditional machine learning models were evaluated independently as well as in a stacked learner approach by using AutoGluon, and various recurrent neural network architectures were considered. The traditional machine learning models were implemented using the AutoGluon-Tabular framework, whereas the recurrent neural networks utilized the TensorFlow Keras framework. We trained these models to operate solely using routine laboratory measurements and clinical covariates available within 72 hours of a patient’s first positive COVID-19 nucleic acid test result.
Results
The GRU-D recurrent neural network achieved peak cross-validation performance with 0.938 (SE 0.004) as the area under the receiver operating characteristic (AUROC) curve. This model retained strong performance by reducing the follow-up time to 12 hours (0.916 SE 0.005 AUROC), and the leave-one-out feature importance analysis indicated that the most independently valuable features were age, Charlson comorbidity index, minimum oxygen saturation, fibrinogen level, and serum iron level. In the prospective testing cohort, this model provided an AUROC of 0.901 and a statistically significant difference in survival (P<.001, hazard ratio for those predicted to survive, 95% CI 0.043-0.106).
Conclusions
Our deep learning approach using GRU-D provides an alert system to flag mortality for COVID-19–positive patients by using clinical covariates and laboratory values within a 72-hour window after the first positive nucleic acid test result.
The differentiation of B cells into antibody secreting plasma cells (PCs) is governed by a strict regulatory network that results in expression of specific transcriptomes along the activation ...continuum. In vitro models yielding significant numbers of PCs phenotypically identical to the in vivo state enable investigation of pathways, metabolomes, and non-coding (ncRNAs) not previously identified. The objective of our study was to characterize ncRNA expression during human B cell activation and differentiation. To achieve this, we used an in vitro system and performed RNA-seq on resting and activated B cells and PCs. Characterization of coding gene transcripts, including immunoglobulin (Ig), validated our system and also demonstrated that memory B cells preferentially differentiated into PCs. Importantly, we identified more than 980 ncRNA transcripts that are differentially expressed across the stages of activation and differentiation, some of which are known to target transcription, proliferation, cytoskeletal, autophagy and proteasome pathways. Interestingly, ncRNAs located within Ig loci may be targeting both Ig and non-Ig-related transcripts. ncRNAs associated with B cell malignancies were also identified. Taken together, this system provides a platform to study the role of specific ncRNAs in B cell differentiation and altered expression of those ncRNAs involved in B cell malignancies.
Antibody secreting plasma cells (PCs) play an important role in effective humoral immune responses. The low frequency of bone marrow PCs in humans makes it challenging to obtain sufficient numbers of ...PCs for biologic studies. Previous studies have employed in vitro model systems to generate cells that morphologically, phenotypically, and functionally resemble normal polyclonal PCs. Gene expression profiles of in vitro generated PCs (IVPCs) mirror their normal counterparts, however to date extensive immunoglobulin (Ig) repertoire analysis of IVPCs is lacking. Here, we used a modified 3-step protocol to generate IVPCs and used RNA-seq to explore the transcriptome with emphasis on the Ig repertoire of plasmablasts and PCs.
Total B cells were isolated from 3 normal donors and cultured with various cytokines and the B cell activators CpG ODN and CD40L. RNA was obtained from freshly isolated B cells (Day 0; D0) as well as from Day 4 (D4) plasmablasts, and Day 10 (D10) IVPCs. Morphologically, D10 cells exhibited typical PC morphology, including an eccentric nucleus and perinuclear hof. RNA-seq was performed on total RNA from all 3 donors and time points using the Standard TRuSeq v2 library prep and with paired end sequencing on the Illumina HiSeq 4000 platform.
Principle component analysis of gene expression data showed that D0, D4 and D10 cells could be clearly segregated across all 3 normal donors. Of importance, transcripts previously described as distinguishing B cells from PCs were found to be differentially expressed including overexpression of CXCR5, CD19, EBF, CD83, PAX5, IRF8 in D0 B cells and overexpression of IRF4, Blimp-1, XBP1, BCMA, SLAMF7, Syndecan-1, CD38 and CD27 in IVPCs, thus validating our in vitro model for generating PCs. Furthermore, expression of cell cycle related transcripts such as CKS1, CDK1, and CCDN2 followed the pattern of low expression in resting B cells, increased expression in plasmablasts, and decreased expression in IVPCs confirming the cells are actively cycling in a manner comparable to cells in vivo. D10 IVPCs also overexpressed transcripts known to be upregulated during the unfolded protein response.
As expected from Ig secreting cells, D10 IVPCs had an over-representation of Ig transcripts. At D0, resting B cells had high levels of IgD and IgM heavy chain (HC) transcripts. At D10, IgM transcripts modestly increased with Log2 fold change (FC) = 3 and as expected, IgD levels decreased significantly (Log2 FC = -2.2). IgA and IgG isotype transcripts significantly increased at D10 (Log2 FC > 6.0) with the IgG4 subtype having the greatest Log2 FC at 8.4.
Next we focused on the Ig repertoire of D0, D4, and D10 cells. By aligning to known germline Ig sequences in IMGT/V-Quest (www.imgt.org) and then assembling the paired ends of D0, D4 and D10 Ig transcripts, we were able to analyze the Ig repertoire. Since the Ig HC variable (V) region is encoded by V, diversity (D) and joining (J) segments, only fragments that could be confidently determined were considered. All but 3 IGHV transcripts (IGHV3-35, IGHV3-47 and IGHV7-8) and 2 IGHD transcripts (IGHD4-4 and IGHD5-5) were found and all IGHJ segments were represented across the differentiation spectrum. In D0 cells, the number of unique VDJ combinations ranged from 643 to 863 across all 3 normal samples and increased to a range of 2524 to 2867 in D10 IVPCs. When looking at the differential expression of each VDJ combination from D0 to D10, a pairwise t-test for relative frequency showed that there was no significant change greater than 1%, suggesting the repertoire diversity was not skewed, thus proving the conditions for stimulation were not targeting any one starting B cell.
Our data also allowed us to track clonal expansions during differentiation as defined by the increasing frequency of sequences with identical nucleotide sequence in the V region and CDR3 (including D and J regions). Hence, a single sequence could be tracked from D0 to D10. Of interest, in a small sampling of the total available sequences, only those B cells with a mutated IGHV region, characteristic of a memory B cell, went on to expand in this system whereas B cells with an unmutated IGHV did not.
Our analysis of the Ig repertoire of IVPCs suggests this system provides a functional model to study Ig repertoire along the B cell differentiation process and further delineate the conditions that may result in a clonal expansion, a hallmark of many hematologic malignancies including multiple myeloma.
No relevant conflicts of interest to declare.
FOXE3 is a lens-specific transcription factor with a highly conserved forkhead domain previously implicated in congenital primary aphakia and anterior segment dysgenesis. Here, we identify new ...recessive FOXE3 mutations causative for microphthalmia, sclerocornea, primary aphakia, and glaucoma in two extended consanguineous families by SNP array genotyping followed by a candidate gene approach. Following an additional screen of 236 subjects with developmental eye anomalies, we report two further novel heterozygous mutations segregating in a dominant fashion in two different families. Although the dominant mutations were penetrant, they gave rise to highly variable phenotypes including iris and chorioretinal colobomas, Peters' anomaly, and isolated cataract (cerulean type and early onset adult nuclear and cortical cataract). Using in situ hybridization in human embryos, we demonstrate expression of FOXE3 restricted to lens tissue, predominantly in the anterior epithelium, suggesting that the extralenticular phenotypes caused by FOXE3 mutations are most likely to be secondary to abnormal lens formation. Our findings suggest that mutations in FOXE3 can give rise to a broad spectrum of eye anomalies, largely, but not exclusively related to lens development, and that both dominant and recessive inheritance patterns can be represented. We suggest including FOXE3 in the diagnostic genetic screening for these anomalies. Hum Mutat 30:1-9, 2009.