Genome and clonal hematopoiesis stability contrasts with immune, cfDNA, mitochondrial, and telomere length changes during short duration spaceflight

Background

The Inspiration4 (I4) mission, the first all-civilian orbital flight mission, investigated the physiological effects of short-duration spaceflight through a multi-omic approach. Despite advances, there remains much to learn about human adaptation to spaceflight's unique challenges, including microgravity, immune system perturbations, and radiation exposure.

Methods

To provide a detailed genetics analysis of the mission, we collected dried blood spots pre-, during, and post-flight for DNA extraction. Telomere length was measured by quantitative PCR, while whole genome and cfDNA sequencing provided insight into genomic stability and immune adaptations. A robust bioinformatic pipeline was used for data analysis, including variant calling to assess mutational burden.

Result

Telomere elongation occurred during spaceflight and shortened after return to Earth. Cell-free DNA analysis revealed increased immune cell signatures post-flight. No significant clonal hematopoiesis of indeterminate potential (CHIP) or whole-genome instability was observed. The long-term gene expression changes across immune cells suggested cellular adaptations to the space environment persisting months post-flight.

Conclusion

Our findings provide valuable insights into the physiological consequences of short-duration spaceflight, with telomere dynamics and immune cell gene expression adapting to spaceflight and persisting after return to Earth. CHIP sequencing data will serve as a reference point for studying the early development of CHIP in astronauts, an understudied phenomenon as previous studies have focused on career astronauts. This study will serve as a reference point for future commercial and non-commercial spaceflight, low Earth orbit (LEO) missions, and deep-space exploration.

Resistance to PRMT5-targeted therapy in mantle cell lymphoma

ABSTRACT

Mantle cell lymphoma (MCL) is an incurable B-cell non-Hodgkin lymphoma, and patients who relapse on targeted therapies have poor prognosis. Protein arginine methyltransferase 5 (PRMT5), an enzyme essential for B-cell transformation, drives multiple oncogenic pathways and is overexpressed in MCL. Despite the antitumor activity of PRMT5 inhibition (PRT-382/PRT-808), drug resistance was observed in a patient-derived xenograft (PDX) MCL model. Decreased survival of mice engrafted with these PRMT5 inhibitor-resistant cells vs treatment-naive cells was observed (P = .005). MCL cell lines showed variable sensitivity to PRMT5 inhibition. Using PRT-382, cell lines were classified as sensitive (n = 4; 50% inhibitory concentration [IC50], 20-140 nM) or primary resistant (n = 4; 340-1650 nM). Prolonged culture of sensitive MCL lines with drug escalation produced PRMT5 inhibitor-resistant cell lines (n = 4; 200-500 nM). This resistant phenotype persisted after prolonged culture in the absence of drug and was observed with PRT-808. In the resistant PDX and cell line models, symmetric dimethylarginine reduction was achieved at the original PRMT5 inhibitor IC50, suggesting activation of alternative resistance pathways. Bulk RNA sequencing of resistant cell lines and PDX relative to sensitive or short-term-treated cells, respectively, highlighted shared upregulation of multiple pathways including mechanistic target of rapamycin kinase [mTOR] signaling (P < 10-5 and z score > 0.3 or < 0.3). Single-cell RNA sequencing analysis demonstrated a strong shift in global gene expression, with upregulation of mTOR signaling in resistant PDX MCL samples. Targeted blockade of mTORC1 with temsirolimus overcame the PRMT5 inhibitor-resistant phenotype, displayed therapeutic synergy in resistant MCL cell lines, and improved survival of a resistant PDX.

Elevated A-to-I RNA editing in COVID-19 infected individuals

Given the current status of coronavirus disease 2019 (COVID-19) as a global pandemic, it is of high priority to gain a deeper understanding of the disease's development and how the virus impacts its host. Adenosine (A)-to-Inosine (I) RNA editing is a post-transcriptional modification, catalyzed by the ADAR family of enzymes, that can be considered part of the inherent cellular defense mechanism as it affects the innate immune response in a complex manner. It was previously reported that various viruses could interact with the host's ADAR enzymes, resulting in epigenetic changes both to the virus and the host. Here, we analyze RNA-seq of nasopharyngeal swab specimens as well as whole-blood samples of COVID-19 infected individuals and show a significant elevation in the global RNA editing activity in COVID-19 compared to healthy controls. We also detect specific coding sites that exhibit higher editing activity. We further show that the increment in editing activity during the disease is temporary and returns to baseline shortly after the symptomatic period. These significant epigenetic changes may contribute to the immune system response and affect adverse outcomes seen in post-viral cases.

Viral activation and ecological restructuring characterize a microbiome axis of spaceflight-associated immune activation

Maintenance of astronaut health during spaceflight will require monitoring and potentially modulating their microbiomes, which play a role in some space-derived health disorders. However, documenting the response of microbiota to spaceflight has been difficult thus far due to mission constraints that lead to limited sampling. Here, we executed a six-month longitudinal study centered on a three-day flight to quantify the high-resolution microbiome response to spaceflight. Via paired metagenomics and metatranscriptomics alongside single immune profiling, we resolved a microbiome "architecture" of spaceflight characterized by time-dependent and taxonomically divergent microbiome alterations across 750 samples and ten body sites. We observed pan-phyletic viral activation and signs of persistent changes that, in the oral microbiome, yielded plaque-associated pathobionts with strong associations to immune cell gene expression. Further, we found enrichments of microbial genes associated with antibiotic production, toxin-antitoxin systems, and stress response enriched universally across the body sites. We also used strain-level tracking to measure the potential propagation of microbial species from the crew members to each other and the environment, identifying microbes that were prone to seed the capsule surface and move between the crew. Finally, we identified associations between microbiome and host immune cell shifts, proposing both a microbiome axis of immune changes during flight as well as the sources of some of those changes. In summary, these datasets and methods reveal connections between crew immunology, the microbiome, and their likely drivers and lay the groundwork for future microbiome studies of spaceflight.

PRMT5 supports multiple oncogenic pathways in mantle cell lymphoma

Mantle cell lymphoma (MCL) is an incurable B-cell malignancy with an overall poor prognosis, particularly for patients that progress on targeted therapies. Novel, more durable treatment options are needed for patients with MCL. Protein arginine methyltransferase 5 (PRMT5) is overexpressed in MCL and plays an important oncogenic role in this disease via epigenetic and posttranslational modification of cell cycle regulators, DNA repair genes, components of prosurvival pathways, and RNA splicing regulators. The mechanism of targeting PRMT5 in MCL remains incompletely characterized. Here, we report on the antitumor activity of PRMT5 inhibition in MCL using integrated transcriptomics of in vitro and in vivo models of MCL. Treatment with a selective small-molecule inhibitor of PRMT5, PRT-382, led to growth arrest and cell death and provided a therapeutic benefit in xenografts derived from patients with MCL. Transcriptional reprograming upon PRMT5 inhibition led to restored regulatory activity of the cell cycle (p-RB/E2F), apoptotic cell death (p53-dependent/p53-independent), and activation of negative regulators of B-cell receptor-PI3K/AKT signaling (PHLDA3, PTPROt, and PIK3IP1). We propose pharmacologic inhibition of PRMT5 for patients with relapsed/refractory MCL and identify MTAP/CDKN2A deletion and wild-type TP53 as biomarkers that predict a favorable response. Selective targeting of PRMT5 has significant activity in preclinical models of MCL and warrants further investigation in clinical trials.

Core mitochondrial genes are down-regulated during SARS-CoV-2 infection of rodent and human hosts

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral proteins bind to host mitochondrial proteins, likely inhibiting oxidative phosphorylation (OXPHOS) and stimulating glycolysis. We analyzed mitochondrial gene expression in nasopharyngeal and autopsy tissues from patients with coronavirus disease 2019 (COVID-19). In nasopharyngeal samples with declining viral titers, the virus blocked the transcription of a subset of nuclear DNA (nDNA)-encoded mitochondrial OXPHOS genes, induced the expression of microRNA 2392, activated HIF-1α to induce glycolysis, and activated host immune defenses including the integrated stress response. In autopsy tissues from patients with COVID-19, SARS-CoV-2 was no longer present, and mitochondrial gene transcription had recovered in the lungs. However, nDNA mitochondrial gene expression remained suppressed in autopsy tissue from the heart and, to a lesser extent, kidney, and liver, whereas mitochondrial DNA transcription was induced and host-immune defense pathways were activated. During early SARS-CoV-2 infection of hamsters with peak lung viral load, mitochondrial gene expression in the lung was minimally perturbed but was down-regulated in the cerebellum and up-regulated in the striatum even though no SARS-CoV-2 was detected in the brain. During the mid-phase SARS-CoV-2 infection of mice, mitochondrial gene expression was starting to recover in mouse lungs. These data suggest that when the viral titer first peaks, there is a systemic host response followed by viral suppression of mitochondrial gene transcription and induction of glycolysis leading to the deployment of antiviral immune defenses. Even when the virus was cleared and lung mitochondrial function had recovered, mitochondrial function in the heart, kidney, liver, and lymph nodes remained impaired, potentially leading to severe COVID-19 pathology.

Geospatially-resolved public-health surveillance via wastewater sequencing

Wastewater, which contains everything from pathogens to pollutants, is a geospatially-and temporally-linked microbial fingerprint of a given population. As a result, it can be leveraged for monitoring multiple dimensions of public health across locales and time. Here, we integrate targeted and bulk RNA sequencing (n=1,419 samples) to track the viral, bacterial, and functional content over geospatially distinct areas within Miami Dade County from 2020-2022. First, we used targeted amplicon sequencing (n=966) to track diverse SARS-CoV-2 variants across space and time, and we found a tight correspondence with clinical caseloads from University students (N = 1,503) and Miami-Dade County hospital patients (N = 3,939 patients), as well as an 8-day earlier detection of the Delta variant in wastewater vs. in patients. Additionally, in 453 metatranscriptomic samples, we demonstrate that different wastewater sampling locations have clinically and public-health-relevant microbiota that vary as a function of the size of the human population they represent. Through assembly, alignment-based, and phylogenetic approaches, we also detect multiple clinically important viruses (e.g., norovirus ) and describe geospatial and temporal variation in microbial functional genes that indicate the presence of pollutants. Moreover, we found distinct profiles of antimicrobial resistance (AMR) genes and virulence factors across campus buildings, dorms, and hospitals, with hospital wastewater containing a significant increase in AMR abundance. Overall, this effort lays the groundwork for systematic characterization of wastewater to improve public health decision making and a broad platform to detect emerging pathogens.

IDH2 and TET2 mutations synergize to modulate T Follicular Helper cell functional interaction with the AITL microenvironment

Angioimmunoblastic T cell lymphoma (AITL) is a peripheral T cell lymphoma that originates from T follicular helper (Tfh) cells and exhibits a prominent tumor microenvironment (TME). IDH2 and TET2 mutations co-occur frequently in AITL, but their contribution to tumorigenesis is poorly understood. We developed an AITL mouse model that is driven by Idh2 and Tet2 mutations. Malignant Tfh cells display aberrant transcriptomic and epigenetic programs that impair TCR signaling. Neoplastic Tfh cells bearing combined Idh2 and Tet2 mutations show altered cross-talk with germinal center B cells that promotes B cell clonal expansion while decreasing Fas-FasL interaction and reducing B cell apoptosis. The plasma cell count and angiogenesis are also increased in the Idh2-mutated tumors, implying a major relationship between Idh2 mutation and the characteristic AITL TME. Our mouse model recapitulates several features of human IDH2-mutated AITL and provides a rationale for exploring therapeutic targeting of Tfh-TME cross-talk for AITL patients.

Distinct and opposite effects of leukemogenic Idh and Tet2 mutations in hematopoietic stem and progenitor cells

Mutations in IDH1, IDH2, and TET2 are recurrently observed in myeloid neoplasms. IDH1 and IDH2 encode isocitrate dehydrogenase isoforms, which normally catalyze the conversion of isocitrate to α-ketoglutarate (α-KG). Oncogenic IDH1/2 mutations confer neomorphic activity, leading to the production of D-2-hydroxyglutarate (D-2-HG), a potent inhibitor of α-KG-dependent enzymes which include the TET methylcytosine dioxygenases. Given their mutual exclusivity in myeloid neoplasms, IDH1, IDH2, and TET2 mutations may converge on a common oncogenic mechanism. Contrary to this expectation, we observed that they have distinct, and even opposite, effects on hematopoietic stem and progenitor cells in genetically engineered mice. Epigenetic and single-cell transcriptomic analyses revealed that Idh2^R172K and Tet2 loss-of-function have divergent consequences on the expression and activity of key hematopoietic and leukemogenic regulators. Notably, chromatin accessibility and transcriptional deregulation in Idh2^R172K cells were partially disconnected from DNA methylation alterations. These results highlight unanticipated divergent effects of IDH1/2 and TET2 mutations, providing support for the optimization of genotype-specific therapies.

Microbiome and metagenomic analysis of Lake Hillier Australia reveals pigment-rich polyextremophiles and wide-ranging metabolic adaptations

Lake Hillier is a hypersaline lake known for its distinctive bright pink color. The cause of this phenomenon in other hypersaline sites has been attributed to halophiles, Dunaliella, and Salinibacter, however, a systematic analysis of the microbial communities, their functional features, and the prevalence of pigment-producing-metabolisms has not been previously studied. Through metagenomic sequencing and culture-based approaches, our results evidence that Lake Hillier is composed of a diverse set of microorganisms including archaea, bacteria, algae, and viruses. Our data indicate that the microbiome in Lake Hillier is composed of multiple pigment-producer microbes, including Dunaliella, Salinibacter, Halobacillus, Psychroflexus, Halorubrum, many of which are cataloged as polyextremophiles. Additionally, we estimated the diversity of metabolic pathways in the lake and determined that many of these are related to pigment production. We reconstructed complete or partial genomes for 21 discrete bacteria (N = 14) and archaea (N = 7), only 2 of which could be taxonomically annotated to previously observed species. Our findings provide the first metagenomic study to decipher the source of the pink color of Australia's Lake Hillier. The study of this pink hypersaline environment is evidence of a microbial consortium of pigment producers, a repertoire of polyextremophiles, a core microbiome and potentially novel species.

Relationships between SARS-CoV-2 in Wastewater and COVID-19 Clinical Cases and Hospitalizations, with and without Normalization against Indicators of Human Waste

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) in wastewater has been used to track community infections of coronavirus disease-2019 (COVID-19), providing critical information for public health interventions. Since levels in wastewater are dependent upon human inputs, we hypothesize that tracking infections can be improved by normalizing wastewater concentrations against indicators of human waste [Pepper Mild Mottle Virus (PMMoV), β-2 Microglobulin (B2M), and fecal coliform]. In this study, we analyzed SARS-CoV-2 and indicators of human waste in wastewater from two sewersheds of different scales: a University campus and a wastewater treatment plant. Wastewater data were combined with complementary COVID-19 case tracking to evaluate the efficiency of wastewater surveillance for forecasting new COVID-19 cases and, for the larger scale, hospitalizations. Results show that the normalization of SARS-CoV-2 levels by PMMoV and B2M resulted in improved correlations with COVID-19 cases for campus data using volcano second generation (V2G)-qPCR chemistry (r s = 0.69 without normalization, r s = 0.73 with normalization). Mixed results were obtained for normalization by PMMoV for samples collected at the community scale. Overall benefits from normalizing with measures of human waste depend upon qPCR chemistry and improves with smaller sewershed scale. We recommend further studies that evaluate the efficacy of additional normalization targets.

TARGETED DOWN REGULATION OF CORE MITOCHONDRIAL GENES DURING SARS-COV-2 INFECTION

Defects in mitochondrial oxidative phosphorylation (OXPHOS) have been reported in COVID-19 patients, but the timing and organs affected vary among reports. Here, we reveal the dynamics of COVID-19 through transcription profiles in nasopharyngeal and autopsy samples from patients and infected rodent models. While mitochondrial bioenergetics is repressed in the viral nasopharyngeal portal of entry, it is up regulated in autopsy lung tissues from deceased patients. In most disease stages and organs, discrete OXPHOS functions are blocked by the virus, and this is countered by the host broadly up regulating unblocked OXPHOS functions. No such rebound is seen in autopsy heart, results in severe repression of genes across all OXPHOS modules. Hence, targeted enhancement of mitochondrial gene expression may mitigate the pathogenesis of COVID-19.

System-wide transcriptome damage and tissue identity loss in COVID-19 patients

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections.

Genomic Characterization of the Titan-like Cell Producing Naganishia tulchinskyi, the First Novel Eukaryote Isolated from the International Space Station

Multiple strains of a novel yeast belonging to genus Naganishia were isolated from environmental surfaces aboard the International Space Station (ISS). These strains exhibited a phenotype similar to Titan cell (~10 µm diameter) morphology when grown under a combination of simulated microgravity and 5% CO₂ conditions. Confocal, scanning, and transmission electron microscopy revealed distinct morphological differences between the microgravity-grown cells and the standard Earth gravity-grown cells, including larger cells and thicker cell walls, altered intracellular morphology, modifications to extracellular fimbriae, budding, and the shedding of bud scars. Phylogenetic analyses via multi-locus sequence typing indicated that these ISS strains represented a single species in the genus Naganishia and were clustered with Naganishia diffluens. The name Naganishia tulchinskyi is proposed to accommodate these strains, with IF6SW-B1^T as the holotype. The gene ontologies were assigned to the cell morphogenesis, microtubule-based response, and response to UV light, suggesting a variety of phenotypes that are well suited to respond to microgravity and radiation. Genomic analyses also indicated that the extracellular region, outer membrane, and cell wall were among the highest cellular component results, thus implying a set of genes associated with Titan-like cell plasticity. Finally, the highest molecular function matches included cytoskeletal motor activity, microtubule motor activity, and nuclear export signal receptor activity.

Assessing reproducibility of inherited variants detected with short-read whole genome sequencing

BACKGROUND

Reproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS.

RESULTS

To dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×.

CONCLUSIONS

Our findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

Betacoronavirus-specific alternate splicing

Viruses can subvert a number of cellular processes including splicing in order to block innate antiviral responses, and many viruses interact with cellular splicing machinery. SARS-CoV-2 infection was shown to suppress global mRNA splicing, and at least 10 SARS-CoV-2 proteins bind specifically to one or more human RNAs. Here, we investigate 17 published experimental and clinical datasets related to SARS-CoV-2 infection, datasets from the betacoronaviruses SARS-CoV and MERS, as well as Streptococcus pneumonia, HCV, Zika virus, Dengue virus, influenza H3N2, and RSV. We show that genes showing differential alternative splicing in SARS-CoV-2 have a similar functional profile to those of SARS-CoV and MERS and affect a diverse set of genes and biological functions, including many closely related to virus biology. Additionally, the differentially spliced transcripts of cells infected by coronaviruses were more likely to undergo intron-retention, contain a pseudouridine modification, and have a smaller number of exons as compared with differentially spliced transcripts in the control groups. Viral load in clinical COVID-19 samples was correlated with isoform distribution of differentially spliced genes. A significantly higher number of ribosomal genes are affected by differential alternative splicing and gene expression in betacoronavirus samples, and the betacoronavirus differentially spliced genes are depleted for binding sites of RNA-binding proteins. Our results demonstrate characteristic patterns of differential splicing in cells infected by SARS-CoV-2, SARS-CoV, and MERS. The alternative splicing changes observed in betacoronaviruses infection potentially modify a broad range of cellular functions, via changes in the functions of the products of a diverse set of genes involved in different biological processes.

Comparative Analysis of Single-Cell RNA Sequencing Platforms and Methods

Single-cell RNA sequencing (scRNA-seq) offers great new opportunities for increasing our understanding of complex biological processes. In particular, development of an accurate Human Cell Atlas is largely dependent on the rapidly advancing technologies and molecular chemistries employed in scRNA-seq. These advances have already allowed an increase in throughput for scRNA-seq from 96 to 80,000 cells on a single instrument run by capturing cells within nanoliter droplets. Although this increase in throughput is critical for many experimental questions, a thorough comparison between microfluidic-based, plate-based, and droplet-based technologies or between multiple available platforms utilizing these technologies is largely lacking. Here, we report scRNA-seq data from SUM149PT cells treated with the histone deacetylase inhibitor trichostatin A versus untreated controls across several scRNA-seq platforms (Fluidigm C1, WaferGen iCell8, 10x Genomics Chromium Controller, and Illumina/BioRad ddSEQ). The primary goal of this project was to demonstrate RNA sequencing methods for profiling the ultra-low amounts of RNA present in individual cells, and this report discusses the results of the study, as well as technical challenges and lessons learned and present general guidelines for best practices in sample preparation and analysis.

The SEQC2 epigenomics quality control (EpiQC) study

BACKGROUND

Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify types of cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and cross-validated resource for epigenetics research from the FDA's Epigenomics Quality Control Group.

RESULTS

Each sample is processed in multiple replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), enzymatic deamination method (EMSeq), targeted methylation sequencing (Illumina Methyl Capture EPIC), single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies, and 850k Illumina methylation arrays. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, bwa-meth, and BitMapperBS), we find overall high concordance between assays, but also differences in efficiency of read mapping, CpG capture, coverage, and platform performance, and variable performance across 26 microarray normalization algorithms.

CONCLUSIONS

The data provided herein can guide the use of these DNA reference materials in epigenomics research, as well as provide best practices for experimental design in future studies. By leveraging seven human cell lines that are designated as publicly available reference materials, these data can be used as a baseline to advance epigenomics research.

Investigation of product-derived lymphoma following infusion of piggyBac-modified CD19 chimeric antigen receptor T cells

We performed a phase 1 clinical trial to evaluate outcomes in patients receiving donor-derived CD19-specific chimeric antigen receptor (CAR) T cells for B-cell malignancy that relapsed or persisted after matched related allogeneic hemopoietic stem cell transplant. To overcome the cost and transgene-capacity limitations of traditional viral vectors, CAR T cells were produced using the piggyBac transposon system of genetic modification. Following CAR T-cell infusion, 1 patient developed a gradually enlarging retroperitoneal tumor due to a CAR-expressing CD4+ T-cell lymphoma. Screening of other patients led to the detection, in an asymptomatic patient, of a second CAR T-cell tumor in thoracic para-aortic lymph nodes. Analysis of the first lymphoma showed a high transgene copy number, but no insertion into typical oncogenes. There were also structural changes such as altered genomic copy number and point mutations unrelated to the insertion sites. Transcriptome analysis showed transgene promoter-driven upregulation of transcription of surrounding regions despite insulator sequences surrounding the transgene. However, marked global changes in transcription predominantly correlated with gene copy number rather than insertion sites. In both patients, the CAR T-cell-derived lymphoma progressed and 1 patient died. We describe the first 2 cases of malignant lymphoma derived from CAR gene-modified T cells. Although CAR T cells have an enviable record of safety to date, our results emphasize the need for caution and regular follow-up of CAR T recipients, especially when novel methods of gene transfer are used to create genetically modified immune therapies. This trial was registered at www.anzctr.org.au as ACTRN12617001579381.

Role of miR-2392 in driving SARS-CoV-2 infection

MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional gene regulation that have a major impact on many diseases and provide an exciting avenue toward antiviral therapeutics. From patient transcriptomic data, we determined that a circulating miRNA, miR-2392, is directly involved with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) machinery during host infection. Specifically, we show that miR-2392 is key in driving downstream suppression of mitochondrial gene expression, increasing inflammation, glycolysis, and hypoxia, as well as promoting many symptoms associated with coronavirus disease 2019 (COVID-19) infection. We demonstrate that miR-2392 is present in the blood and urine of patients positive for COVID-19 but is not present in patients negative for COVID-19. These findings indicate the potential for developing a minimally invasive COVID-19 detection method. Lastly, using in vitro human and in vivo hamster models, we design a miRNA-based antiviral therapeutic that targets miR-2392, significantly reduces SARS-CoV-2 viability in hamsters, and may potentially inhibit a COVID-19 disease state in humans.

DNA methylation-calling tools for Oxford Nanopore sequencing: a survey and human epigenome-wide evaluation

BACKGROUND

Nanopore long-read sequencing technology greatly expands the capacity of long-range, single-molecule DNA-modification detection. A growing number of analytical tools have been developed to detect DNA methylation from nanopore sequencing reads. Here, we assess the performance of different methylation-calling tools to provide a systematic evaluation to guide researchers performing human epigenome-wide studies.

RESULTS

We compare seven analytic tools for detecting DNA methylation from nanopore long-read sequencing data generated from human natural DNA at a whole-genome scale. We evaluate the per-read and per-site performance of CpG methylation prediction across different genomic contexts, CpG site coverage, and computational resources consumed by each tool. The seven tools exhibit different performances across the evaluation criteria. We show that the methylation prediction at regions with discordant DNA methylation patterns, intergenic regions, low CG density regions, and repetitive regions show room for improvement across all tools. Furthermore, we demonstrate that 5hmC levels at least partly contribute to the discrepancy between bisulfite and nanopore sequencing. Lastly, we provide an online DNA methylation database ( https://nanome.jax.org ) to display the DNA methylation levels detected by nanopore sequencing and bisulfite sequencing data across different genomic contexts.

CONCLUSIONS

Our study is the first systematic benchmark of computational methods for detection of mammalian whole-genome DNA modifications in nanopore sequencing. We provide a broad foundation for cross-platform standardization and an evaluation of analytical tools designed for genome-scale modified base detection using nanopore sequencing.

Epigenetic Forensics for Suspect Identification and Age Prediction

Background: Genetic testing at crime scenes is an instrumental molecular technique to identify or eliminate suspects, as well as to overturn wrongful convictions. Yet, genotyping alone cannot reveal the age of a sample, which could help advance the utility of crime scene samples for suspect identification. The distribution of cytosine methylation within a DNA sample can be leveraged to determine the epigenetic age of someone's blood. Methodology: We sought to demonstrate the ability of DNA methylation markers to accurately discern the age of blood spots from an actual crime scene, a "mock" crime scene, and also from a tube of blood stored in ethylenediaminetetraacetic acid for >20 years. This was achieved by quantifying methylation within known age-associated genetic loci across each DNA sample. We observed a strong linear coefficient (0.91) and high overall correlation (R ² = 0.963) between the known age of a sample and the predicted age. Conclusion: We show that novel methods for targeted methylation and low-input whole-genome bisulfite sequencing can enable a novel and improved forensic profile of a crime scene that discerns not only who was present at the crime, but also their age. Finally, we use this model to discern the age and provenance of a blood sample that was used in a criminal investigation.

Performance assessment of DNA sequencing platforms in the ABRF Next-Generation Sequencing Study

Assessing the reproducibility, accuracy and utility of massively parallel DNA sequencing platforms remains an ongoing challenge. Here the Association of Biomolecular Resource Facilities (ABRF) Next-Generation Sequencing Study benchmarks the performance of a set of sequencing instruments (HiSeq/NovaSeq/paired-end 2 × 250-bp chemistry, Ion S5/Proton, PacBio circular consensus sequencing (CCS), Oxford Nanopore Technologies PromethION/MinION, BGISEQ-500/MGISEQ-2000 and GS111) on human and bacterial reference DNA samples. Among short-read instruments, HiSeq 4000 and X10 provided the most consistent, highest genome coverage, while BGI/MGISEQ provided the lowest sequencing error rates. The long-read instrument PacBio CCS had the highest reference-based mapping rate and lowest non-mapping rate. The two long-read platforms PacBio CCS and PromethION/MinION showed the best sequence mapping in repeat-rich areas and across homopolymers. NovaSeq 6000 using 2 × 250-bp read chemistry was the most robust instrument for capturing known insertion/deletion events. This study serves as a benchmark for current genomics technologies, as well as a resource to inform experimental design and next-generation sequencing variant calling.

Molecular Characterization of Increased Amplicon Lengths in SARS-CoV-2 Reverse Transcription Loop-Mediated Isothermal Amplification Assays

Loop-mediated isothermal amplification (LAMP) is a power tool for the amplification of specific RNA and DNA targets. Much like PCR, LAMP requires primers that surround a target amplification region and generates exponential product through a unique highly specific daisy-chain single-temperature amplification reaction. However, until recently, attempts to amplify targets of greater than 200 base pairs (bp) have been mostly unsuccessful and limited to short amplicon targets of less than 150 bp. Although short amplicons have the benefit of a rapid detection (<40 min), they do not allow for the prediction of RNA integrity based on RNA length and possible intactness. In this study, 8 primer sets were developed using 2 LAMP primer-specific software packages against the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nucleocapsid gene with insert lengths ranging from 262 to 945 bp in order to amplify and infer the integrity of viral RNA. Because these amplification lengths are greater than the current methods that use an insert length of 130 or less, they require a longer incubation, modified primer and temperature strategies, and the addition of specific adjuncts to prevent nonspecific amplification. This proof of concept study resulted in successful reverse transcription LAMP reactions for amplicon targets of 262, 687, 693, and 945 bp using a clinical nasopharyngeal NP sample, purified SARS-CoV-2 RNA, and crude lysate containing inactivated virus.

Establishing community reference samples, data and call sets for benchmarking cancer mutation detection using whole-genome sequencing

The lack of samples for generating standardized DNA datasets for setting up a sequencing pipeline or benchmarking the performance of different algorithms limits the implementation and uptake of cancer genomics. Here, we describe reference call sets obtained from paired tumor-normal genomic DNA (gDNA) samples derived from a breast cancer cell line-which is highly heterogeneous, with an aneuploid genome, and enriched in somatic alterations-and a matched lymphoblastoid cell line. We partially validated both somatic mutations and germline variants in these call sets via whole-exome sequencing (WES) with different sequencing platforms and targeted sequencing with >2,000-fold coverage, spanning 82% of genomic regions with high confidence. Although the gDNA reference samples are not representative of primary cancer cells from a clinical sample, when setting up a sequencing pipeline, they not only minimize potential biases from technologies, assays and informatics but also provide a unique resource for benchmarking 'tumor-only' or 'matched tumor-normal' analyses.

The Great Deceiver: miR-2392's Hidden Role in Driving SARS-CoV-2 Infection

MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional gene regulation that have a major impact on many diseases and provides an exciting avenue towards antiviral therapeutics. From patient transcriptomic data, we have discovered a circulating miRNA, miR-2392, that is directly involved with SARS-CoV-2 machinery during host infection. Specifically, we show that miR-2392 is key in driving downstream suppression of mitochondrial gene expression, increasing inflammation, glycolysis, and hypoxia as well as promoting many symptoms associated with COVID-19 infection. We demonstrate miR-2392 is present in the blood and urine of COVID-19 positive patients, but not detected in COVID-19 negative patients. These findings indicate the potential for developing a novel, minimally invasive, COVID-19 detection method. Lastly, using in vitro human and in vivo hamster models, we have developed a novel miRNA-based antiviral therapeutic that targets miR-2392, significantly reduces SARS-CoV-2 viability in hamsters and may potentially inhibit a COVID-19 disease state in humans.

Betacoronavirus-specific alternate splicing

Viruses can subvert a number of cellular processes in order to block innate antiviral responses, and many viruses interact with cellular splicing machinery. SARS-CoV-2 infection was shown to suppress global mRNA splicing, and at least 10 SARS-CoV-2 proteins bind specifically to one or more human RNAs. Here, we investigate 17 published experimental and clinical datasets related to SARS-CoV-2 infection as well as datasets from the betacoronaviruses SARS-CoV and MERS as well as Streptococcus pneumonia, HCV, Zika virus, Dengue virus, influenza H3N2, and RSV. We show that genes showing differential alternative splicing in SARS-CoV-2 have a similar functional profile to those of SARS-CoV and MERS and affect a diverse set of genes and biological functions, including many closely related to virus biology. Additionally, the differentially spliced transcripts of cells infected by coronaviruses were more likely to undergo intron-retention, contain a pseudouridine modification and a smaller number of exons than differentially spliced transcripts in the control groups. Viral load in clinical COVID-19 samples was correlated with isoform distribution of differentially spliced genes. A significantly higher number of ribosomal genes are affected by DAS and DGE in betacoronavirus samples, and the betacoronavirus differentially spliced genes are depleted for binding sites of RNA-binding proteins. Our results demonstrate characteristic patterns of differential splicing in cells infected by SARS-CoV-2, SARS-CoV, and MERS, potentially modifying a broad range of cellular functions and affecting a diverse set of genes and biological functions.

Haplotype diversity and sequence heterogeneity of human telomeres

Telomeres are regions of repetitive nucleotide sequences capping the ends of eukaryotic chromosomes that protect against deterioration, and whose lengths can be correlated with age and adverse health risk factors. Yet, given their length and repetitive nature, telomeric regions are not easily reconstructed from short-read sequencing, thus making telomere sequencing, mapping, and variant resolution challenging problems. Recently, long-read sequencing, with read lengths measuring in hundreds of kilobase pairs, has made it possible to routinely read into telomeric regions and inspect their sequence structure. Here, we describe a framework for extracting telomeric reads from whole-genome single-molecule sequencing experiments, including de novo identification of telomere repeat motifs and repeat types, and also describe their sequence variation. We find that long, complex telomeric stretches and repeats can be accurately captured with long-read sequencing, observe extensive sequence heterogeneity of human telomeres, discover and localize noncanonical telomere sequence motifs (both previously reported, as well as novel), and validate them in short-read sequence data. These data reveal extensive intra- and inter-population diversity of repeats in telomeric haplotypes, reveal higher paternal inheritance of telomeric variants, and represent the first motif composition maps of multi-kilobase-pair human telomeric haplotypes across three distinct ancestries (Ashkenazi, Chinese, and Utah), which can aid in future studies of genetic variation, aging, and genome biology.

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.

Fundamental Biological Features of Spaceflight: Advancing the Field to Enable Deep-Space Exploration

Research on astronaut health and model organisms have revealed six features of spaceflight biology that guide our current understanding of fundamental molecular changes that occur during space travel. The features include oxidative stress, DNA damage, mitochondrial dysregulation, epigenetic changes (including gene regulation), telomere length alterations, and microbiome shifts. Here we review the known hazards of human spaceflight, how spaceflight affects living systems through these six fundamental features, and the associated health risks of space exploration. We also discuss the essential issues related to the health and safety of astronauts involved in future missions, especially planned long-duration and Martian missions.

The spatial landscape of lung pathology during COVID-19 progression

Recent studies have provided insights into the pathology of and immune response to COVID-191-8. However, a thorough investigation of the interplay between infected cells and the immune system at sites of infection has been lacking. Here we use high-parameter imaging mass cytometry9 that targets the expression of 36 proteins to investigate the cellular composition and spatial architecture of acute lung injury in humans (including injuries derived from SARS-CoV-2 infection) at single-cell resolution. These spatially resolved single-cell data unravel the disordered structure of the infected and injured lung, alongside the distribution of extensive immune infiltration. Neutrophil and macrophage infiltration are hallmarks of bacterial pneumonia and COVID-19, respectively. We provide evidence that SARS-CoV-2 infects predominantly alveolar epithelial cells and induces a localized hyperinflammatory cell state that is associated with lung damage. We leverage the temporal range of fatal outcomes of COVID-19 in relation to the onset of symptoms, which reveals increased macrophage extravasation and increased numbers of mesenchymal cells and fibroblasts concomitant with increased proximity between these cell types as the disease progresses-possibly as a result of attempts to repair the damaged lung tissue. Our data enable us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. We use this landscape to characterize the pathophysiology of the human lung from its macroscopic presentation to the single-cell level, which provides an important basis for understanding COVID-19 and lung pathology in general.

Cross-oncopanel study reveals high sensitivity and accuracy with overall analytical performance depending on genomic regions

BACKGROUND

Targeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing.

RESULTS

All panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5-20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden.

CONCLUSION

This comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.

A verified genomic reference sample for assessing performance of cancer panels detecting small variants of low allele frequency

BACKGROUND

Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance.

RESULTS

In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5-100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels.

CONCLUSION

These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

SARS-CoV-2 genomic diversity and the implications for qRT-PCR diagnostics and transmission

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 lineages that have spread throughout the world. In this study, we investigated 129 RNA-seq data sets and 6928 consensus genomes to contrast the intra-host and inter-host diversity of SARS-CoV-2. Our analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights intra-host single nucleotide variant (iSNV) and SNP similarity, albeit with differences in C > U changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of insertions and deletions contribute to the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

Shotgun transcriptome, spatial omics, and isothermal profiling of SARS-CoV-2 infection reveals unique host responses, viral diversification, and drug interactions

In less than nine months, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) killed over a million people, including >25,000 in New York City (NYC) alone. The COVID-19 pandemic caused by SARS-CoV-2 highlights clinical needs to detect infection, track strain evolution, and identify biomarkers of disease course. To address these challenges, we designed a fast (30-minute) colorimetric test (LAMP) for SARS-CoV-2 infection from naso/oropharyngeal swabs and a large-scale shotgun metatranscriptomics platform (total-RNA-seq) for host, viral, and microbial profiling. We applied these methods to clinical specimens gathered from 669 patients in New York City during the first two months of the outbreak, yielding a broad molecular portrait of the emerging COVID-19 disease. We find significant enrichment of a NYC-distinctive clade of the virus (20C), as well as host responses in interferon, ACE, hematological, and olfaction pathways. In addition, we use 50,821 patient records to find that renin-angiotensin-aldosterone system inhibitors have a protective effect for severe COVID-19 outcomes, unlike similar drugs. Finally, spatial transcriptomic data from COVID-19 patient autopsy tissues reveal distinct ACE2 expression loci, with macrophage and neutrophil infiltration in the lungs. These findings can inform public health and may help develop and drive SARS-CoV-2 diagnostic, prevention, and treatment strategies.

Systemic Tissue and Cellular Disruption from SARS-CoV-2 Infection revealed in COVID-19 Autopsies and Spatial Omics Tissue Maps

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus has infected over 115 million people and caused over 2.5 million deaths worldwide. Yet, the molecular mechanisms underlying the clinical manifestations of COVID-19, as well as what distinguishes them from common seasonal influenza virus and other lung injury states such as Acute Respiratory Distress Syndrome (ARDS), remains poorly understood. To address these challenges, we combined transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues, matched with spatial protein and expression profiling (GeoMx) across 357 tissue sections. These results define both body-wide and tissue-specific (heart, liver, lung, kidney, and lymph nodes) damage wrought by the SARS-CoV-2 infection, evident as a function of varying viral load (high vs. low) during the course of infection and specific, transcriptional dysregulation in splicing isoforms, T cell receptor expression, and cellular expression states. In particular, cardiac and lung tissues revealed the largest degree of splicing isoform switching and cell expression state loss. Overall, these findings reveal a systemic disruption of cellular and transcriptional pathways from COVID-19 across all tissues, which can inform subsequent studies to combat the mortality of COVID-19, as well to better understand the molecular dynamics of lethal SARS-CoV-2 infection and other viruses.

A diagnostic host response biosignature for COVID-19 from RNA profiling of nasal swabs and blood

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes coronavirus disease-19 (COVID-19), has emerged as the cause of a global pandemic. We used RNA sequencing to analyze 286 nasopharyngeal (NP) swab and 53 whole-blood (WB) samples from 333 patients with COVID-19 and controls. Overall, a muted immune response was observed in COVID-19 relative to other infections (influenza, other seasonal coronaviruses, and bacterial sepsis), with paradoxical down-regulation of several key differentially expressed genes. Hospitalized patients and outpatients exhibited up-regulation of interferon-associated pathways, although heightened and more robust inflammatory responses were observed in hospitalized patients with more clinically severe illness. Two-layer machine learning-based host classifiers consisting of complete (>1000 genes), medium (<100), and small (<20) gene biomarker panels identified COVID-19 disease with 85.1-86.5% accuracy when benchmarked using an independent test set. SARS-CoV-2 infection has a distinct biosignature that differs between NP swabs and WB and can be leveraged for COVID-19 diagnosis.

Transcriptional response modules characterize IL-1β and IL-6 activity in COVID-19

Dysregulated IL-1β and IL-6 responses have been implicated in the pathogenesis of severe Coronavirus Disease 2019 (COVID-19). Innovative approaches for evaluating the biological activity of these cytokines in vivo are urgently needed to complement clinical trials of therapeutic targeting of IL-1β and IL-6 in COVID-19. We show that the expression of IL-1β or IL-6 inducible transcriptional signatures (modules) reflects the bioactivity of these cytokines in immunopathology modelled by juvenile idiopathic arthritis (JIA) and rheumatoid arthritis. In COVID-19, elevated expression of IL-1β and IL-6 response modules, but not the cytokine transcripts themselves, is a feature of infection in the nasopharynx and blood but is not associated with severity of COVID-19 disease, length of stay, or mortality. We propose that IL-1β and IL-6 transcriptional response modules provide a dynamic readout of functional cytokine activity in vivo, aiding quantification of the biological effects of immunomodulatory therapies in COVID-19.

Cell-free DNA (cfDNA) and Exosome Profiling from a Year-Long Human Spaceflight Reveals Circulating Biomarkers

Liquid biopsies based on cell-free DNA (cfDNA) or exosomes provide a noninvasive approach to monitor human health and disease but have not been utilized for astronauts. Here, we profile cfDNA characteristics, including fragment size, cellular deconvolution, and nucleosome positioning, in an astronaut during a year-long mission on the International Space Station, compared to his identical twin on Earth and healthy donors. We observed a significant increase in the proportion of cell-free mitochondrial DNA (cf-mtDNA) inflight, and analysis of post-flight exosomes in plasma revealed a 30-fold increase in circulating exosomes and patient-specific protein cargo (including brain-derived peptides) after the year-long mission. This longitudinal analysis of astronaut cfDNA during spaceflight and the exosome profiles highlights their utility for astronaut health monitoring, as well as cf-mtDNA levels as a potential biomarker for physiological stress or immune system responses related to microgravity, radiation exposure, and the other unique environmental conditions of spaceflight.

Transcriptional response modules characterise IL-1β and IL-6 activity in COVID-19

Dysregulated IL-1β and IL-6 responses have been implicated in the pathogenesis of severe Coronavirus Disease 2019 (COVID-19). Innovative approaches for evaluating the biological activity of these cytokines in vivo are urgently needed to complement clinical trials of therapeutic targeting of IL-1β and IL-6 in COVID-19. We show that the expression of IL-1β or IL-6 inducible transcriptional signatures (modules) reflects the bioactivity of these cytokines in immunopathology modelled by juvenile idiopathic arthritis (JIA) and rheumatoid arthritis. In COVID-19, elevated expression of IL-1β and IL-6 response modules, but not the cytokine transcripts themselves, is a feature of infection in the nasopharynx and blood, but is not associated with severity of COVID-19 disease, length of stay or mortality. We propose that IL-1β and IL-6 transcriptional response modules provide a dynamic readout of functional cytokine activity in vivo, aiding quantification of the biological effects of immunomodulatory therapies in COVID-19.

Clonal Hematopoiesis Before, During, and After Human Spaceflight

Clonal hematopoiesis (CH) occurs when blood cells harboring an advantageous mutation propagate faster than others. These mutations confer a risk for hematological cancers and cardiovascular disease. Here, we analyze CH in blood samples from a pair of twin astronauts over 4 years in bulk and fractionated cell populations using a targeted CH panel, linked-read whole-genome sequencing, and deep RNA sequencing. We show CH with distinct mutational profiles and increasing allelic fraction that includes a high-risk, TET2 clone in one subject and two DNMT3A mutations on distinct alleles in the other twin. These astronauts exhibit CH almost two decades prior to the mean age at which it is typically detected and show larger shifts in clone size than age-matched controls or radiotherapy patients, based on a longitudinal cohort of 157 cancer patients. As such, longitudinal monitoring of CH may serve as an important metric for overall cancer and cardiovascular risk in astronauts.

Trinchant et al. examined twin astronauts for clonal hematopoiesis (CH). Some high-risk CH clones (TET2 and DNMT3A) were observed two decades before expected, with TET2 decreasing in spaceflight and elevating later post flight. Thus, CH is an important metric for overall cancer and cardiovascular risk in astronauts.

Temporal Telomere and DNA Damage Responses in the Space Radiation Environment

Telomeres, repetitive terminal features of chromosomes essential for maintaining genome integrity, shorten with cell division, lifestyle factors and stresses, and environmental exposures, and so they provide a robust biomarker of health, aging, and age-related diseases. We assessed telomere length dynamics (changes over time) in three unrelated astronauts before, during, and after 1-year or 6-month missions aboard the International Space Station (ISS). Similar to our results for National Aeronautics and Space Administration's (NASA's) One-Year Mission twin astronaut (Garrett-Bakelman et al., 2019), significantly longer telomeres were observed during spaceflight for two 6-month mission astronauts. Furthermore, telomere length shortened rapidly after return to Earth for all three crewmembers and, overall, telomere length tended to be shorter after spaceflight than before spaceflight. Consistent with chronic exposure to the space radiation environment, signatures of persistent DNA damage responses were also detected, including mitochondrial and oxidative stress, inflammation, and telomeric and chromosomal aberrations, which together provide potential mechanistic insight into spaceflight-specific telomere elongation.

Immune complement and coagulation dysfunction in adverse outcomes of SARS-CoV-2 infection

Understanding the pathophysiology of SARS-CoV-2 infection is critical for therapeutic and public health strategies. Viral-host interactions can guide discovery of disease regulators, and protein structure function analysis points to several immune pathways, including complement and coagulation, as targets of coronaviruses. To determine whether conditions associated with dysregulated complement or coagulation systems impact disease, we performed a retrospective observational study and found that history of macular degeneration (a proxy for complement-activation disorders) and history of coagulation disorders (thrombocytopenia, thrombosis and hemorrhage) are risk factors for SARS-CoV-2-associated morbidity and mortality-effects that are independent of age, sex or history of smoking. Transcriptional profiling of nasopharyngeal swabs demonstrated that in addition to type-I interferon and interleukin-6-dependent inflammatory responses, infection results in robust engagement of the complement and coagulation pathways. Finally, in a candidate-driven genetic association study of severe SARS-CoV-2 disease, we identified putative complement and coagulation-associated loci including missense, eQTL and sQTL variants of critical complement and coagulation regulators. In addition to providing evidence that complement function modulates SARS-CoV-2 infection outcome, the data point to putative transcriptional genetic markers of susceptibility. The results highlight the value of using a multimodal analytical approach to reveal determinants and predictors of immunity, susceptibility and clinical outcome associated with infection.

Hidden genomic diversity of SARS-CoV-2: implications for qRT-PCR diagnostics and transmission

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 clades that have spread throughout the world. In this study, we investigated over 7,000 SARS-CoV-2 datasets to unveil both intrahost and interhost diversity. Our intrahost and interhost diversity analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights iSNV and SNP similarity, albeit with high variability in C>T changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of small indels fuel the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

Identification of Radioresponsive Genes in Esophageal Cancer from Longitudinal and Single Cell Exome Sequencing

PURPOSE

The majority (70%) of the esophageal squamous cell carcinoma (ESCC) cases in the world occur in China, where radiation therapy is the most common treatment. Yet the majority of ESCC patients still relapse.

METHODS AND MATERIALS

To better understand the genetic basis of radiation therapy resistance for ESCC, we performed longitudinal, whole-exome sequencing throughout radiation therapy on 42 patient tumor samples, including single-cell whole-exome sequencing for 147 cells for 2 patients.

RESULTS

Significant allelic changes were observed during clinical irradiation, with 42 recurrent radioresponsive genes (sensitive and resistant) identified in multiple patients, including NOTCH1, MAML3, CDKN2A, NFE2L2, GAS2L2, OBSCN and TP53, with the last 3 genes implicated as radioresponsive in both bulk and single-cell whole-exome sequencing. Most (37/42) radioresponsive genes showed regional variegation in both radioresistant and radiosensitive mutations, with a paucity of resistant-only mutations (2.5%). A subset of sensitive mutations in 10 genes and resistant mutations in 18 genes defined a significantly improved prognosis and the shortest time for locoregional recurrence, respectively, indicating possible clinical utility. We also confirmed these significant mutational signatures in orthogonal Cancer Genome Atlas ESCC cohorts.

CONCLUSIONS

Overall, our results quantify the allelic shifts underlying radioresponse in bulk and single-cell ESCC exomes for the first time, provide a temporal resolution to such mutational dynamics, and offer new therapeutic target genes and loci for esophageal and potentially other cancers.

Identification of Immune complement function as a determinant of adverse SARS-CoV-2 infection outcome

Understanding the pathophysiology of SARS-CoV-2 infection is critical for therapeutics and public health intervention strategies. Viral-host interactions can guide discovery of regulators of disease outcomes, and protein structure function analysis points to several immune pathways, including complement and coagulation, as targets of the coronavirus proteome. To determine if conditions associated with dysregulation of the complement or coagulation systems impact adverse clinical outcomes, we performed a retrospective observational study of 11,116 patients who presented with suspected SARS-CoV-2 infection. We found that history of macular degeneration (a proxy for complement activation disorders) and history of coagulation disorders (thrombocytopenia, thrombosis, and hemorrhage) are risk factors for morbidity and mortality in SARS-CoV-2 infected patients - effects that could not be explained by age, sex, or history of smoking. Further, transcriptional profiling of nasopharyngeal (NP) swabs from 650 control and SARS-CoV-2 infected patients demonstrated that in addition to innate Type-I interferon and IL-6 dependent inflammatory immune responses, infection results in robust engagement and activation of the complement and coagulation pathways. Finally, we conducted a candidate driven genetic association study of severe SARS-CoV-2 disease. Among the findings, our scan identified putative complement and coagulation associated loci including missense, eQTL and sQTL variants of critical regulators of the complement and coagulation cascades. In addition to providing evidence that complement function modulates SARS-CoV-2 infection outcome, the data point to putative transcriptional genetic markers of susceptibility. The results highlight the value of using a multi-modal analytical approach, combining molecular information from virus protein structure-function analysis with clinical informatics, transcriptomics, and genomics to reveal determinants and predictors of immunity, susceptibility, and clinical outcome associated with infection.

Illuminating Genetic Mysteries of the Dead Sea Scrolls

The discovery of the 2,000-year-old Dead Sea Scrolls had an incomparable impact on the historical understanding of Judaism and Christianity. "Piecing together" scroll fragments is like solving jigsaw puzzles with an unknown number of missing parts. We used the fact that most scrolls are made from animal skins to "fingerprint" pieces based on DNA sequences. Genetic sorting of the scrolls illuminates their textual relationship and historical significance. Disambiguating the contested relationship between Jeremiah fragments supplies evidence that some scrolls were brought to the Qumran caves from elsewhere; significantly, they demonstrate that divergent versions of Jeremiah circulated in parallel throughout Israel (ancient Judea). Similarly, patterns discovered in non-biblical scrolls, particularly the Songs of the Sabbath Sacrifice, suggest that the Qumran scrolls represent the broader cultural milieu of the period. Finally, genetic analysis divorces debated fragments from the Qumran scrolls. Our study demonstrates that interdisciplinary approaches enrich the scholar's toolkit.