Comprehensive Multi-omics Analysis Reveals Mitochondrial Stress as a Central Biological Hub for Spaceflight Impact.

Spaceflight is known to impose changes on human physiology with unknown molecular etiologies. To reveal these causes, we used a multi-omics, systems biology analytical approach using biomedical profiles from fifty-nine astronauts and data from NASA's GeneLab derived from hundreds of samples flown in space to determine transcriptomic, proteomic, metabolomic, and epigenetic responses to spaceflight. Overall pathway analyses on the multi-omics datasets showed significant enrichment for mitochondrial processes, as well as innate immunity, chronic inflammation, cell cycle, circadian rhythm, and olfactory functions. Importantly, NASA's Twin Study provided a platform to confirm several of our principal findings. Evidence of altered mitochondrial function and DNA damage was also found in the urine and blood metabolic data compiled from the astronaut cohort and NASA Twin Study data, indicating mitochondrial stress as a consistent phenotype of spaceflight.

A Multiscale Map of the Stem Cell State in Pancreatic Adenocarcinoma.

Drug resistance and relapse remain key challenges in pancreatic cancer. Here, we have used RNA sequencing (RNA-seq), chromatin immunoprecipitation (ChIP)-seq, and genome-wide CRISPR analysis to map the molecular dependencies of pancreatic cancer stem cells, highly therapy-resistant cells that preferentially drive tumorigenesis and progression. This integrated genomic approach revealed an unexpected utilization of immuno-regulatory signals by pancreatic cancer epithelial cells. In particular, the nuclear hormone receptor retinoic-acid-receptor-related orphan receptor gamma (RORγ), known to drive inflammation and T cell differentiation, was upregulated during pancreatic cancer progression, and its genetic or pharmacologic inhibition led to a striking defect in pancreatic cancer growth and a marked improvement in survival. Further, a large-scale retrospective analysis in patients revealed that RORγ expression may predict pancreatic cancer aggressiveness, as it positively correlated with advanced disease and metastasis. Collectively, these data identify an orthogonal co-option of immuno-regulatory signals by pancreatic cancer stem cells, suggesting that autoimmune drugs should be evaluated as novel treatment strategies for pancreatic cancer patients.

Non-coding Transcription Instructs Chromatin Folding and Compartmentalization to Dictate Enhancer-Promoter Communication and T Cell Fate.

It is now established that Bcl11b specifies T cell fate. Here, we show that in developing T cells, the Bcl11b enhancer repositioned from the lamina to the nuclear interior. Our search for factors that relocalized the Bcl11b enhancer identified a non-coding RNA named ThymoD (thymocyte differentiation factor). ThymoD-deficient mice displayed a block at the onset of T cell development and developed lymphoid malignancies. We found that ThymoD transcription promoted demethylation at CTCF bound sites and activated cohesin-dependent looping to reposition the Bcl11b enhancer from the lamina to the nuclear interior and to juxtapose the Bcl11b enhancer and promoter into a single-loop domain. These large-scale changes in nuclear architecture were associated with the deposition of activating epigenetic marks across the loop domain, plausibly facilitating phase separation. These data indicate how, during developmental progression and tumor suppression, non-coding transcription orchestrates chromatin folding and compartmentalization to direct with high precision enhancer-promoter communication.

Nuclear morphology is shaped by loop-extrusion programs.

It is well established that neutrophils adopt malleable polymorphonuclear shapes to migrate through narrow interstitial tissue spaces1-3. However, how polymorphonuclear structures are assembled remains unknown4. Here we show that in neutrophil progenitors, halting loop extrusion-a motor-powered process that generates DNA loops by pulling in chromatin5-leads to the assembly of polymorphonuclear genomes. Specifically, we found that in mononuclear neutrophil progenitors, acute depletion of the loop-extrusion loading factor nipped-B-like protein (NIPBL) induced the assembly of horseshoe, banded, ringed and hypersegmented nuclear structures and led to a reduction in nuclear volume, mirroring what is observed during the differentiation of neutrophils. Depletion of NIPBL also induced cell-cycle arrest, activated a neutrophil-specific gene program and conditioned a loss of interactions across topologically associating domains to generate a chromatin architecture that resembled that of differentiated neutrophils. Removing NIPBL resulted in enrichment for mega-loops and interchromosomal hubs that contain genes associated with neutrophil-specific enhancer repertoires and an inflammatory gene program. On the basis of these observations, we propose that in neutrophil progenitors, loop-extrusion programs produce lineage-specific chromatin architectures that permit the packing of chromosomes into geometrically confined lobular structures. Our data also provide a blueprint for the assembly of polymorphonuclear structures, and point to the possibility of engineering de novo nuclear shapes to facilitate the migration of effector cells in densely populated tumorigenic environments.

Variant STAT4 and Response to Ruxolitinib in an Autoinflammatory Syndrome.

BACKGROUND: Disabling pansclerotic morphea (DPM) is a rare systemic inflammatory disorder, characterized by poor wound healing, fibrosis, cytopenias, hypogammaglobulinemia, and squamous-cell carcinoma. The cause is unknown, and mortality is high. METHODS: We evaluated four patients from three unrelated families with an autosomal dominant pattern of inheritance of DPM. Genomic sequencing independently identified three heterozygous variants in a specific region of the gene that encodes signal transducer and activator of transcription 4 (STAT4). Primary skin fibroblast and cell-line assays were used to define the functional nature of the genetic defect. We also assayed gene expression using single-cell RNA sequencing of peripheral-blood mononuclear cells to identify inflammatory pathways that may be affected in DPM and that may respond to therapy. RESULTS: Genome sequencing revealed three novel heterozygous missense gain-of-function variants in STAT4. In vitro, primary skin fibroblasts showed enhanced interleukin-6 secretion, with impaired wound healing, contraction of the collagen matrix, and matrix secretion. Inhibition of Janus kinase (JAK)-STAT signaling with ruxolitinib led to improvement in the hyperinflammatory fibroblast phenotype in vitro and resolution of inflammatory markers and clinical symptoms in treated patients, without adverse effects. Single-cell RNA sequencing revealed expression patterns consistent with an immunodysregulatory phenotype that were appropriately modified through JAK inhibition. CONCLUSIONS: Gain-of-function variants in STAT4 caused DPM in the families that we studied. The JAK inhibitor ruxolitinib attenuated the dermatologic and inflammatory phenotype in vitro and in the affected family members. (Funded by the American Academy of Allergy, Asthma, and Immunology Foundation and others.).

The genomic landscape of Ménière's disease: a path to endolymphatic hydrops.

BACKGROUND: Ménière's disease (MD) is a disorder of the inner ear that causes episodic bouts of severe dizziness, roaring tinnitus, and fluctuating hearing loss. To date, no targeted therapy exists. As such, we have undertaken a large whole genome sequencing study on carefully phenotyped unilateral MD patients with the goal of gene/pathway discovery and a move towards targeted intervention. This study was a retrospective review of patients with a history of Ménière's disease. Genomic DNA, acquired from saliva samples, was purified and subjected to whole genome sequencing. RESULTS: Stringent variant calling, performed on 511 samples passing quality checks, followed by gene-based filtering by recurrence and proximity in molecular interaction networks, led to 481 high priority MD genes. These high priority genes, including MPHOSPH8, MYO18A, TRIOBP, OTOGL, TNC, and MYO6, were previously implicated in hearing loss, balance, and cochlear function, and were significantly enriched in common variant studies of hearing loss. Validation in an independent MD cohort confirmed 82 recurrent genes. Pathway analysis pointed to cell-cell adhesion, extracellular matrix, and cellular energy maintenance as key mediators of MD. Furthermore, the MD-prioritized genes were highly expressed in human inner ear hair cells and dark/vestibular cells, and were differentially expressed in a mouse model of hearing loss. CONCLUSION: By enabling the development of model systems that may lead to targeted therapies and MD screening panels, the genes and variants identified in this study will inform diagnosis and treatment of MD.

Defining the relationship of salivary gland malignancies to novel cell subpopulations in human salivary glands using single nucleus RNA-sequencing.

Salivary glands have essential roles in maintaining oral health, mastication, taste and speech, by secreting saliva. Salivary glands are composed of several types of cells, and each cell type is predicted to be involved in the carcinogenesis of different types of cancers including adenoid cystic carcinoma (ACC), acinic cell carcinoma (AciCC), salivary duct carcinoma (SDC), myoepithelial carcinoma (MECA) and other histology. In our study, we performed single nucleus RNA-seq on three human salivary gland samples to clarify the gene expression profile of each complex cellular component of the salivary glands and related these expression patterns to expression found in salivary gland cancers (SGC) to infer cell of origin. By single nucleus RNA-seq, salivary gland cells were stratified into four clusters: acinar cells, ductal cells 1, ductal cells 2 and myoepithelial cells/stromal cells. The localization of each cell group was verified by IHC of each cluster marker gene, and one group of ductal cells was found to represent intercalated ductal cells labeled with HES1. Furthermore, in comparison with SGC RNA-seq data, acinar cell markers were upregulated in AciCC, but downregulated in ACC and ductal cell markers were upregulated in SDC but downregulated in MECA, suggesting that markers of origin are highly expressed in some SGC. Cell type expressions in specific SGC histology are similar to those found in normal salivary gland populations, indicating a potential etiologic relationship.

Nonalcoholic fatty liver disease risk and histologic severity are associated with genetic polymorphisms in children.

BACKGROUND AND AIMS: NAFLD is the most common chronic liver disease in children. Large pediatric studies identifying single nucleotide polymorphisms (SNPs) associated with risk and histologic severity of NAFLD are limited. Study aims included investigating SNPs associated with risk for NAFLD using family trios and association of candidate alleles with histologic severity. APPROACH AND RESULTS: Children with biopsy-confirmed NAFLD were enrolled from the NASH Clinical Research Network. The Expert Pathology Committee reviewed liver histology. Genotyping was conducted with allele-specific primers for 60 candidate SNPs. Parents were enrolled for trio analysis. To assess risk for NAFLD, the transmission disequilibrium test was conducted in trios. Among cases, regression analysis assessed associations with histologic severity. A total of 822 children with NAFLD had mean age 13.2 years (SD 2.7) and mean ALT 101 U/L (SD 90). PNPLA3 (rs738409) demonstrated the strongest risk ( p = 2.24 × 10 -14 ) for NAFLD. Among children with NAFLD, stratifying by PNPLA3 s738409 genotype, the variant genotype associated with steatosis ( p = 0.005), lobular ( p = 0.03) and portal inflammation ( p = 0.002). Steatosis grade associated with TM6SF2 ( p = 0.0009), GCKR ( p = 0.0032), PNPLA3 rs738409 ( p = 0.0053), and MTTP ( p = 0.0051). Fibrosis stage associated with PARVB rs6006473 ( p = 0.0001), NR1I2 ( p = 0.0021), ADIPOR2 ( p = 0.0038), and OXTR ( p = 0.0065). PNPLA3 rs738409 ( p = 0.0002) associated with borderline zone 1 NASH. CONCLUSIONS: This study demonstrated disease-associated SNPs in children with NAFLD. In particular, rs6006473 was highly associated with severity of fibrosis. These hypothesis-generating results support future mechanistic studies of development of adverse outcomes such as fibrosis and generation of therapeutic targets for NAFLD in children.

Interleukin-18 signaling promotes activation of hepatic stellate cells in mouse liver fibrosis.

BACKGROUND AND AIMS: Nucleotide-binding oligomerization domain-like receptor-family pyrin domain-containing 3 (NLRP3) inflammasome activation has been shown to result in liver fibrosis. Mechanisms and downstream signaling remain incompletely understood. Here, we studied the role of IL-18 in hepatic stellate cells (HSCs), and its impact on liver fibrosis. APPROACH AND RESULTS: We observed significantly increased serum levels of IL-18 (128.4 pg/ml vs. 74.9 pg/ml) and IL-18 binding protein (BP; 46.50 ng/ml vs. 15.35 ng/ml) in patients with liver cirrhosis compared with healthy controls. Single cell RNA sequencing data showed that an immunoregulatory subset of murine HSCs highly expresses Il18 and Il18r1 . Treatment of cultured primary murine HSC with recombinant mouse IL-18 accelerated their transdifferentiation into myofibroblasts. In vivo , IL-18 receptor-deficient mice had reduced liver fibrosis in a model of fibrosis induced by HSC-specific NLRP3 overactivation. Whole liver RNA sequencing analysis from a murine model of severe NASH-induced fibrosis by feeding a choline-deficient, L-amino acid-defined, high fat diet showed that genes related to IL-18 and its downstream signaling were significantly upregulated, and Il18-/- mice receiving this diet for 10 weeks showed protection from fibrotic changes with decreased number of alpha smooth muscle actin-positive cells and collagen deposition. HSC activation triggered by NLRP3 inflammasome activation was abrogated when IL-18 signaling was blocked by its naturally occurring antagonist IL-18BP. Accordingly, we observed that the severe inflammatory phenotype associated with myeloid cell-specific NLRP3 gain-of-function was rescued by IL-18BP. CONCLUSIONS: Our study highlights the role of IL-18 in the development of liver fibrosis by its direct effect on HSC activation identifying IL-18 as a target to treat liver fibrosis.

Comprehensive next-generation sequencing identifies novel putative pathogenic or likely pathogenic germline variants in patients with concurrent tubo-ovarian and endometrial serous and endometrioid carcinomas or precursors.

BACKGROUND: Endometrial serous carcinoma (ESC) and tubo-ovarian high-grade serous carcinoma (HGSC) are characterized by late-stage presentation and high mortality. Current guidelines for prevention recommend risk-reducing salpingo-oophorectomy (RRSO) in patients with hereditary mutations in cancer susceptibility genes. However, HGSC displays extensive genetic heterogeneity with alterations in 168 genes identified in TCGA study, but current germline testing panels are often limited to the handful of recurrently mutated genes, leaving families with rare hereditary gene mutations potentially at-risk. OBJECTIVE: To determine if there are rare germline mutations that may aid in early identification of more patients at-risk for ESC and/or HGSC by evaluating patients with concurrent ESC, HGSC or precursor lesions, and endometrial atypical hyperplasia (CAH) or low-grade endometrial endometrioid adenocarcinoma (LGEEA). METHODS: We performed targeted next-generation sequencing using TSO 500, a 523 gene panel, on formalin-fixed paraffin-embedded tumor and matched benign non-tumor tissue blocks from 5 patients with concurrent ESC, HGSC or precursor lesions, and CAH or LGEEA. RESULTS: We identified germline pathogenic, likely pathogenic or uncertain significance variants in cancer susceptibility genes in 4 of 5 patients - affected genes included GLI1, PIK3R1, FOXP1, FANCD2, INPP4B and H3F3C. Notably, none of these genes were included in the commercially available germline testing panels initially used to evaluate the patients at the time of their diagnoses. CONCLUSION: Comprehensive germline testing of patients with concurrent LGEEA or CAH and ESC, HGSC or precursor lesions may aid in early identification of relatives at-risk for cancer who may be candidates for RRSO with hysterectomy.

Chronic villitis as a distinctive feature of placental injury in maternal SARS-CoV-2 infection.

BACKGROUND: SARS-CoV-2 infection during pregnancy is associated with an increased risk for stillbirth, preeclampsia, and preterm birth. However, this does not seem to be caused by intrauterine fetal infection because vertical transmission is rarely reported. There is a paucity of data regarding the associated placental SARS-CoV-2 histopathology and their relationship with the timing and severity of infection. OBJECTIVE: This study aimed to determine if maternal SARS-CoV-2 infection was associated with specific patterns of placental injury and if these findings differed by gestational age at time of infection or disease severity. STUDY DESIGN: A retrospective cohort study was performed at the University of California San Diego between March 2020 and February 2021. Placentas from pregnancies with a positive SARS-CoV-2 test were matched with 2 sets of controls; 1 set was time-matched by delivery date and sent to pathology for routine clinical indications, and the other was chosen from a cohort of placentas previously collected for research purposes without clinical indications for pathologic examination before the SARS-CoV-2 outbreak. Placental pathologic lesions were defined based on standard criteria and included maternal and fetal vascular malperfusion and acute and chronic inflammatory lesions. A bivariate analysis was performed using the independent Student t test and Pearson chi-square test. A logistic regression was used to control for relevant covariates. Regions of SARS-CoV-2-associated villitis were further investigated using protein-based digital spatial profiling assays on the GeoMx platform, validated by immunohistochemistry, and compared with cases of infectious villitis and villitis of unknown etiology. Differential expression analysis was performed to identify protein expression differences between these groups of villitis. RESULTS: We included 272 SARS-CoV-2 positive cases, 272 time-matched controls, and 272 historic controls. The mean age of SARS-CoV-2 affected subjects was 30.1±5.5 years and the majority were Hispanic (53.7%) and parous (65.7%). SARS-CoV-2 placentas demonstrated a higher frequency of the 4 major patterns of placental injury (all P<.001) than the historic controls. SARS-CoV-2 placentas also showed a higher frequency of chronic villitis and severe chronic villitis (P=.03 for both) than the time-matched controls, which remained significant after controlling for gestational age at delivery (adjusted odds ratio, 1.52; 95% confidence interval, 1.01-2.28; adjusted odds ratio, 2.12; 95% confidence interval, 1.16-3.88, respectively). Digital spatial profiling revealed that programmed death-ligand 1 was increased in villitis-positive regions of the SARS-CoV-2 (logFC, 0.47; adjusted P value =.002) and villitis of unknown etiology (logFC, 0.58; adjusted P value =.003) cases, but it was conversely decreased in villitis-positive regions of the infectious villitis group (log FC, -1.40; adjusted P value <.001). CONCLUSION: Chronic villitis seems to be the most specific histopathologic finding associated with SARS-CoV-2 maternal infection. Chronic villitis involves damage to the vasculosyncytial membrane of the chorionic villi, which are involved in gas and nutrient exchange, suggesting potential mechanisms of placental (and perhaps neonatal) injury, even in the absence of vertical transmission. Surprisingly, changes in protein expression in SARS-CoV-2-associated villitis seem to be more similar to villitis of unknown etiology than to infectious villitis.

The Lung Microenvironment Instructs Gene Transcription in Neonatal and Adult Alveolar Macrophages.

Immaturity of alveolar macrophages (AMs) around birth contributes to the susceptibility of newborns to lung disease. However, the molecular features differentiating neonatal and mature, adult AMs are poorly understood. In this study, we identify the unique transcriptomes and enhancer landscapes of neonatal and adult AMs in mice. Although the core AM signature was similar, murine adult AMs expressed higher levels of genes involved in lipid metabolism, whereas neonatal AMs expressed a largely proinflammatory gene profile. Open enhancer regions identified by an assay for transposase-accessible chromatin followed by high-throughput sequencing (ATAC-seq) contained motifs for nuclear receptors, MITF, and STAT in adult AMs and AP-1 and NF-κB in neonatal AMs. Intranasal LPS activated a similar innate immune response in both neonatal and adult mice, with higher basal expression of inflammatory genes in neonates. The lung microenvironment drove many of the distinguishing gene expression and open chromatin characteristics of neonatal and adult AMs. Neonatal mouse AMs retained high expression of some proinflammatory genes, suggesting that the differences in neonatal AMs result from both inherent cell properties and environmental influences.

Activation of NF-κB and p300/CBP potentiates cancer chemoimmunotherapy through induction of MHC-I antigen presentation.

Many cancers evade immune rejection by suppressing major histocompatibility class I (MHC-I) antigen processing and presentation (AgPP). Such cancers do not respond to immune checkpoint inhibitor therapies (ICIT) such as PD-1/PD-L1 [PD-(L)1] blockade. Certain chemotherapeutic drugs augment tumor control by PD-(L)1 inhibitors through potentiation of T-cell priming but whether and how chemotherapy enhances MHC-I-dependent cancer cell recognition by cytotoxic T cells (CTLs) is not entirely clear. We now show that the lysine acetyl transferases p300/CREB binding protein (CBP) control MHC-I AgPPM expression and neoantigen amounts in human cancers. Moreover, we found that two distinct DNA damaging drugs, the platinoid oxaliplatin and the topoisomerase inhibitor mitoxantrone, strongly up-regulate MHC-I AgPP in a manner dependent on activation of nuclear factor kappa B (NF-κB), p300/CBP, and other transcription factors, but independently of autocrine IFNγ signaling. Accordingly, NF-κB and p300 ablations prevent chemotherapy-induced MHC-I AgPP and abrogate rejection of low MHC-I-expressing tumors by reinvigorated CD8+ CTLs. Drugs like oxaliplatin and mitoxantrone may be used to overcome resistance to PD-(L)1 inhibitors in tumors that had "epigenetically down-regulated," but had not permanently lost MHC-I AgPP activity.

Oxidative stress diverts tRNA synthetase to nucleus for protection against DNA damage.

Tyrosyl-tRNA synthetase (TyrRS) is known for its essential aminoacylation function in protein synthesis. Here we report a function for TyrRS in DNA damage protection. We found that oxidative stress, which often downregulates protein synthesis, induces TyrRS to rapidly translocate from the cytosol to the nucleus. We also found that angiogenin mediates or potentiates this stress-induced translocalization. The nuclear-localized TyrRS activates transcription factor E2F1 to upregulate the expression of DNA damage repair genes such as BRCA1 and RAD51. The activation is achieved through direct interaction of TyrRS with TRIM28 to sequester this vertebrate-specific epigenetic repressor and its associated HDAC1 from deacetylating and suppressing E2F1. Remarkably, overexpression of TyrRS strongly protects against UV-induced DNA double-strand breaks in zebrafish, whereas restricting TyrRS nuclear entry completely abolishes the protection. Therefore, oxidative stress triggers an essential cytoplasmic enzyme used for protein synthesis to translocate to the nucleus to protect against DNA damage.

Efficient population-scale variant analysis and prioritization with VAPr.

Summary: With the growing availability of population-scale whole-exome and whole-genome sequencing, demand for reproducible, scalable variant analysis has spread within genomic research communities. To address this need, we introduce the Python package Variant Analysis and Prioritization (VAPr). VAPr leverages existing annotation tools ANNOVAR and MyVariant.info with MongoDB-based flexible storage and filtering functionality. It offers biologists and bioinformatics generalists easy-to-use and scalable analysis and prioritization of genomic variants from large cohort studies. Availability and implementation: VAPr is developed in Python and is available for free use and extension under the MIT License. An install package is available on PyPi at https://pypi.python.org/pypi/VAPr, while source code and extensive documentation are on GitHub at https://github.com/ucsd-ccbb/VAPr.

Interactive network visualization in Jupyter notebooks: visJS2jupyter.

Motivation: Network biology is widely used to elucidate mechanisms of disease and biological processes. The ability to interact with biological networks is important for hypothesis generation and to give researchers an intuitive understanding of the data. We present visJS2jupyter, a tool designed to embed interactive networks in Jupyter notebooks to streamline network analysis and to promote reproducible research. Results: The tool provides functions for performing and visualizing useful network operations in biology, including network overlap, network propagation around a focal set of genes, and co-localization of two sets of seed genes. visJS2jupyter uses the JavaScript library vis.js to create interactive networks displayed within Jupyter notebook cells with features including drag, click, hover, and zoom. We demonstrate the functionality of visJS2jupyter applied to a biological question, by creating a network propagation visualization to prioritize risk-related genes in autism. Availability and implementation: The visJS2jupyter package is distributed under the MIT License. The source code, documentation and installation instructions are freely available on GitHub at https://github.com/ucsd-ccbb/visJS2jupyter. The package can be downloaded at https://pypi.python.org/pypi/visJS2jupyter. Contact: sbrosenthal@ucsd.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

A Conservation Hatchery Population of Delta Smelt Shows Evidence of Genetic Adaptation to Captivity After 9 Generations.

Genetic adaptation to captivity is a concern for threatened and endangered species held in conservation hatcheries. Here, we present evidence of genetic adaptation to captivity in a conservation hatchery for the endangered delta smelt (Fish Conservation and Culture Laboratory, University of California Davis; FCCL). The FCCL population is genetically managed with parentage analysis and the addition of wild fish each year. Molecular monitoring indicates little loss of genetic variation and low differentiation between the wild and conservation populations. Yet, we found an increase in offspring survival to reproductive maturity during the subsequent spawning season (recovery rate) in crosses that included one or both cultured parents. Crosses with higher levels of hatchery ancestry tend to produce a greater number of offspring that are recovered the following year. The recovery rate of a cross decreases when offspring are raised in a tank with fish of high levels of hatchery ancestry. We suggest changes in fish rearing practices at the FCCL to reduce genetic adaptation to captivity, as delta smelt numbers in the wild continue to decline and the use of FCCL fish for reintroduction becomes more likely.

Omics Pipe: a community-based framework for reproducible multi-omics data analysis.

MOTIVATION: Omics Pipe (http://sulab.scripps.edu/omicspipe) is a computational framework that automates multi-omics data analysis pipelines on high performance compute clusters and in the cloud. It supports best practice published pipelines for RNA-seq, miRNA-seq, Exome-seq, Whole-Genome sequencing, ChIP-seq analyses and automatic processing of data from The Cancer Genome Atlas (TCGA). Omics Pipe provides researchers with a tool for reproducible, open source and extensible next generation sequencing analysis. The goal of Omics Pipe is to democratize next-generation sequencing analysis by dramatically increasing the accessibility and reproducibility of best practice computational pipelines, which will enable researchers to generate biologically meaningful and interpretable results. RESULTS: Using Omics Pipe, we analyzed 100 TCGA breast invasive carcinoma paired tumor-normal datasets based on the latest UCSC hg19 RefSeq annotation. Omics Pipe automatically downloaded and processed the desired TCGA samples on a high throughput compute cluster to produce a results report for each sample. We aggregated the individual sample results and compared them to the analysis in the original publications. This comparison revealed high overlap between the analyses, as well as novel findings due to the use of updated annotations and methods. AVAILABILITY AND IMPLEMENTATION: Source code for Omics Pipe is freely available on the web (https://bitbucket.org/sulab/omics_pipe). Omics Pipe is distributed as a standalone Python package for installation (https://pypi.python.org/pypi/omics_pipe) and as an Amazon Machine Image in Amazon Web Services Elastic Compute Cloud that contains all necessary third-party software dependencies and databases (https://pythonhosted.org/omics_pipe/AWS_installation.html).

Cord Blood Adductomics Reveals Oxidative Stress Exposure Pathways of Bronchopulmonary Dysplasia.

Fetal and neonatal exposures to perinatal oxidative stress (OS) are key mediators of bronchopulmonary dysplasia (BPD). To characterize these exposures, adductomics is an exposure science approach that captures electrophilic addition products (adducts) in blood protein. Adducts are bound to the nucleophilic cysteine loci of human serum albumin (HSA), which has a prolonged half-life. We conducted targeted and untargeted adductomics to test the hypothesis that adducts of OS vary with BPD. We studied 205 preterm infants (≤28 weeks) and 51 full-term infants from an ongoing birth cohort. Infant plasma was collected at birth (cord blood), 1-week, 1-month, and 36-weeks postmenstrual age. HSA was isolated from plasma, trypsin digested, and analyzed using high-performance liquid chromatography-mass spectrometry to quantify previously annotated (known) and unknown adducts. We identified 105 adducts in cord and postnatal blood. A total of 51 known adducts (small thiols, direct oxidation products, and reactive aldehydes) were increased with BPD. Postnatally, serial concentrations of several known OS adducts correlated directly with supplemental oxygen exposure. The application of large-scale adductomics elucidated OS-mediated pathways of BPD. This is the first study to investigate the "neonatal-perinatal exposome" and to identify oxidative stress-related exposure biomarkers that may inform antioxidant strategies to protect the health of future generations of infants.

RNA editing regulates host immune response and T cell homeostasis in SARS-CoV-2 infection.

Adenosine to inosine (A-to-I) RNA editing by ADAR1 has been implicated in maintaining self-tolerance, preventing autoimmunity, and mediating antiviral immunity. Foreign viral double-stranded RNA triggers rapid interferon response and activates ADAR1 in the host immune system. Emerging data points to a role of ADAR1 A-to-I editing in the inflammatory response associated with severe COVID-19 disease. We identify A-to-I editing events within human whole transcriptome data from SARS-CoV-2 infected individuals, non-infected individuals, and individuals with other viral illnesses from nasopharyngeal swabs. High levels of RNA editing in host cells are associated with low SARS-CoV-2 viral load (p = 9.27 E-06), suggesting an inhibitory effect of ADAR1 on viral infection. Additionally, we find differentially expressed genes associated with RNA-modifications and interferon response. Single cell RNA-sequencing analysis of SARS-CoV-2 infected nasopharyngeal swabs reveals that cytotoxic CD8 T cells upregulate ADAR1 in COVID-19 positive samples (p = 0.0269). We further reveal ADAR1 expression increases with CD4 and CD8 T cell activation, and knockdown of ADAR1 leads to apoptosis and aberrant IL-2 secretion. Together, our data suggests A-to-I RNA editing is required to maintain healthy homeostasis of activated T cells to combat SARS-CoV-2 infection.