The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

With the rapid advancement of sequencing technologies, next generation sequencing (NGS) analysis has been widely applied in cancer genomics research. More recently, NGS has been adopted in clinical oncology to advance personalized medicine. Clinical applications of precision oncology require accurate tests that can distinguish tumor-specific mutations from artifacts introduced during NGS processes or data analysis. Therefore, there is an urgent need to develop best practices in cancer mutation detection using NGS and the need for standard reference data sets for systematically measuring accuracy and reproducibility across platforms and methods. Within the SEQC2 consortium context, we established paired tumor-normal reference samples and generated whole-genome (WGS) and whole-exome sequencing (WES) data using sixteen library protocols, seven sequencing platforms at six different centers. We systematically interrogated somatic mutations in the reference samples to identify factors affecting detection reproducibility and accuracy in cancer genomes. These large cross-platform/site WGS and WES datasets using well-characterized reference samples will represent a powerful resource for benchmarking NGS technologies, bioinformatics pipelines, and for the cancer genomics studies.

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.

Crewed National Aeronautics and Space Administration (NASA) missions to other solar system bodies are currently being planned. One high-profile scientific focus during such expeditions would be life detection, specifically the discovery of past or present microbial life, if they exist. However, both humans and associated objects typically carry a high microbial burden. Thus, it is essential to distinguish between microbes brought with the expedition and those present on the exploring planets. Modern spacesuits are unique, customized spacecraft which provide protection, mobility and life support to crew during spacewalks, yet they vent, and the mobility of microbes through spacesuits has not been studied.

The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples.

Here, we performed shotgun metagenome sequencing of swab samples collected on floors at a train station in Narita City, Chiba, Japan. The taxonomic analysis revealed that Actinobacteria and Proteobacteria were the dominant phyla. The data will contribute to insight into the microbiome community on the surfaces of urban built environments.

Multiple sclerosis (MS) is a complex disease of the CNS thought to require an environmental trigger. Gut dysbiosis is common in MS, but specific causative species are unknown. To address this knowledge gap, we used sensitive and quantitative PCR detection to show that people with MS were more likely to harbor and show a greater abundance of epsilon toxin-producing (ETX-producing) strains of C. perfringens within their gut microbiomes compared with individuals who are healthy controls (HCs). Isolates derived from patients with MS produced functional ETX and had a genetic architecture typical of highly conjugative plasmids. In the active immunization model of experimental autoimmune encephalomyelitis (EAE), where pertussis toxin (PTX) is used to overcome CNS immune privilege, ETX can substitute for PTX. In contrast to PTX-induced EAE, where inflammatory demyelination is largely restricted to the spinal cord, ETX-induced EAE caused demyelination in the corpus callosum, thalamus, cerebellum, brainstem, and spinal cord, more akin to the neuroanatomical lesion distribution seen in MS. CNS endothelial cell transcriptional profiles revealed ETX-induced genes that are known to play a role in overcoming CNS immune privilege. Together, these findings suggest that ETX-producing C. perfringens strains are biologically plausible pathogens in MS that trigger inflammatory demyelination in the context of circulating myelin autoreactive lymphocytes.

Strains of Agrobacterium genomospecies 3 (i.e., genomovar G3 of the Agrobacterium tumefaciens species complex) have been previously isolated from diverse environments, including in association with plant roots, with algae, as part of a lignocellulose degrading community, from a hospital environment, as a human opportunistic pathogen, or as reported in this study, from a surface within the International Space Station. Polyphasic taxonomic methods revealed the relationship of Agrobacterium G3 strains to other Agrobacterium spp., which supports the description of a novel species. The G3 strains tested (n = 9) were phenotypically distinguishable among the strains from other genomospecies of the genus Agrobacterium. Phylogenetic analyses of the 16S rRNA gene, gyrB gene, multi-locus sequence analysis, and 1,089-gene core-genome gene concatenate concur that tested G3 strains belong to the Agrobacterium genus and they form a clade distinct from other validly described Agrobacterium species. The distinctiveness of this clade was confirmed by average nucleotide identity (ANI) and in silico digital DNA-DNA hybridization (dDDH) comparisons between the G3 tested strains and all known Agrobacterium species type strains, since obtained values were considerably below the 95% (ANI) and 70% (dDDH) thresholds used for the species delineation. According to the core-genome phylogeny and ANI comparisons, the closest relatives of G3 strains were Agrobacterium sp. strains UGM030330-04 and K599, members of a novel genomospecies we propose to call genomovar G21. Using this polyphasic approach, we characterized the phenotypic and genotypic synapomorphies of Agrobacterium G3, showing it is a bona fide bacterial species, well separated from previously named Agrobacterium species or other recognized genomic species. We thus propose the name Agrobacterium tomkonis for this species previously referred to as Agrobacterium genomospecies 3. The type strain of A. tomkonis is IIF1SW-B1T (= LMG 32164 = NRRL B-65602). Comparative genomic analysis show A. tomkonis strains have species-specific genes associated with secretion of secondary metabolites, including an exopolysaccharide and putative adhesins and resistance to copper. A. tomkonis specific gene functions notably relate to surface adhesion and could be involved to colonize nutrient-poor and harsh habitats. The A. tomkonis strains from the ISS showed presence of a 40-kbp plasmid and several other potential mobile genetic elements detected that could also be part of conjugative elements or integrated prophages.

In addition to mutations, epigenetic alterations are important contributors to malignant transformation and tumor progression. The aim of this work was to identify epigenetic events in which promoter or gene body DNA methylation induces gene expression changes that drive melanocyte malignant transformation and metastasis. We previously developed a linear mouse model of melanoma progression consisting of spontaneously immortalized melanocytes, premalignant melanocytes, a nonmetastatic tumorigenic, and a metastatic cell line. Here, through the integrative analysis of methylome and transcriptome data, we identified the relationship between promoter and/or gene body DNA methylation alterations and gene expression in early, intermediate, and late stages of melanoma progression. We identified adenylate cyclase type 3 (Adcy3) and inositol polyphosphate 4-phosphatase type II (Inpp4b), which affect tumor growth and metastatic potential, respectively. Importantly, the gene expression and DNA methylation profiles found in this murine model of melanoma progression were correlated with available clinical data from large population-based primary melanoma cohorts, revealing potential prognostic markers.

The epigenetic changes associated with melanoma progression to advanced and metastatic stages are still poorly understood. To shed light on the CpG methylation dynamics during melanoma development, we analyzed the methylome profiles of a four-stage cell line model of melanoma progression: non-tumorigenic melanocytes (melan-a), premalignant melanocytes (4C), non-metastatic melanoma cells (4C11-), and metastatic melanoma cells (4C11+). We identified 540 hypo- and 37 hypermethylated gene promoters that together characterized a malignancy signature, and 646 hypo- and 520 hypermethylated promoters that distinguished a metastasis signature. Differentially methylated genes from these signatures were correlated with overall survival using TCGA-SKCM methylation data. Moreover, multivariate Cox analyses with LASSO regularization identified panels of 33 and 31 CpGs, respectively, from the malignancy and metastasis signatures that predicted poor survival. We found a concordant relationship between DNA methylation and transcriptional levels for genes from the malignancy (Pyroxd2 and Ptgfrn) and metastasis (Arnt2, Igfbp4 and Ptprf) signatures, which were both also correlated with melanoma prognosis. Altogether, this study reveals novel CpGs methylation markers associated with malignancy and metastasis that collectively could improve the survival prediction of melanoma patients.

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.

Molecular methods have been used to detect human pathogens in wastewater with sampling typically performed at wastewater treatment plants (WWTP) and upstream locations within the sewer system. A wastewater-based surveillance (WBS) program was established at the University of Miami (UM) in 2020, which included measurements of SARS-CoV-2 levels in wastewater from its hospital and within the regional WWTP. In addition to the development of a SARS-CoV-2 quantitative PCR (qPCR) assay, qPCR assays to detect other human pathogens of interest were also developed at UM. Here we report on the use of a modified set of reagents published by the CDC to detect nucleic acids of Monkeypox virus (MPXV) which emerged during May of 2022 to become a concern worldwide. Samples collected from the University hospital and from the regional WWTP were processed through DNA and RNA workflows and analyzed by qPCR to detect a segment of the MPXV CrmB gene. Results show positive detections of MPXV nucleic acids in the hospital and wastewater treatment plant wastewater which coincided with clinical cases in the community and mirrored the overall trend of nationwide MPXV cases reported to the CDC. We recommend the expansion of current WBS programs' methods to detect a broader range of pathogens of concern in wastewater and present evidence that viral RNA in human cells infected by a DNA virus can be detected in wastewater.

With the advent of long-term human habitation in space and on the moon, understanding how the built environment microbiome of space habitats differs from Earth habitats, and how microbes survive, proliferate and spread in space conditions, is becoming more important. The microbial tracking mission series has been monitoring the microbiome of the International Space Station (ISS) for almost a decade. During this mission series, six unique strains of Gram-stain-positive bacteria, including two spore-forming and three non-spore-forming species, were isolated from the environmental surfaces of the ISS. The analysis of their 16S rRNA gene sequences revealed > 99% similarities with previously described bacterial species. To further explore their phylogenetic affiliation, whole genome sequencing was undertaken. For all strains, the gyrB gene exhibited < 93% similarity with closely related species, which proved effective in categorizing these ISS strains as novel species. Average nucleotide identity and digital DNA-DNA hybridization values, when compared to any known bacterial species, were < 94% and <50% respectively for all species described here. Traditional biochemical tests, fatty acid profiling, polar lipid, and cell wall composition analyses were performed to generate phenotypic characterization of these ISS strains. A study of the shotgun metagenomic reads from the ISS samples, from which the novel species were isolated, showed that only 0.1% of the total reads mapped to the novel species, supporting the idea that these novel species are rare in the ISS environments. In-depth annotation of the genomes unveiled a variety of genes linked to amino acid and derivative synthesis, carbohydrate metabolism, cofactors, vitamins, prosthetic groups, pigments, and protein metabolism. Further analysis of these ISS-isolated organisms revealed that, on average, they contain 46 genes associated with virulence, disease, and defense. The main predicted functions of these genes are: conferring resistance to antibiotics and toxic compounds, and enabling invasion and intracellular resistance. After conducting antiSMASH analysis, it was found that there are roughly 16 cluster types across the six strains, including β-lactone and type III polyketide synthase (T3PKS) clusters. Based on these multi-faceted taxonomic methods, it was concluded that these six ISS strains represent five novel species, which we propose to name as follows: Arthrobacter burdickii IIF3SC-B10T (= NRRL B-65660T = DSM 115933T), Leifsonia virtsii F6_8S_P_1AT (= NRRL B-65661T = DSM 115931T), Leifsonia williamsii F6_8S_P_1BT (= NRRL B-65662T = DSM 115932T), Paenibacillus vandeheii F6_3S_P_1CT (= NRRL B-65663T = DSM 115940T), and Sporosarcina highlanderae F6_3S_P_2T (= NRRL B-65664T = DSM 115943T). Identifying and characterizing the genomes and phenotypes of novel microbes found in space habitats, like those explored in this study, is integral for expanding our genomic databases of space-relevant microbes. This approach offers the only reliable method to determine species composition, track microbial dispersion, and anticipate potential threats to human health from monitoring microbes on the surfaces and equipment within space habitats. By unraveling these microbial mysteries, we take a crucial step towards ensuring the safety and success of future space missions.

An essential step in the discovery of molecular mechanisms contributing to disease phenotypes and efficient experimental planning is the development of weighted hypotheses that estimate the functional effects of sequence variants discovered by high-throughput genomics. With the increasing specialization of the bioinformatics resources, creating analytical workflows that seamlessly integrate data and bioinformatics tools developed by multiple groups becomes inevitable. Here we present a case study of a use of the distributed analytical environment integrating four complementary specialized resources, namely the Lynx platform, VISTA RViewer, the Developmental Brain Disorders Database (DBDB), and the RaptorX server, for the identification of high-confidence candidate genes contributing to pathogenesis of spina bifida. The analysis resulted in prediction and validation of deleterious mutations in the SLC19A placental transporter in mothers of the affected children that causes narrowing of the outlet channel and therefore leads to the reduced folate permeation rate. The described approach also enabled correct identification of several genes, previously shown to contribute to pathogenesis of spina bifida, and suggestion of additional genes for experimental validations. The study demonstrates that the seamless integration of bioinformatics resources enables fast and efficient prioritization and characterization of genomic factors and molecular networks contributing to the phenotypes of interest.

Obtaining RNA-seq measurements involves a complex data analytical process with a large number of competing algorithms as options. There is much debate about which of these methods provides the best approach. Unfortunately, it is currently difficult to evaluate their performance due in part to a lack of sensitive assessment metrics. We present a series of statistical summaries and plots to evaluate the performance in terms of specificity and sensitivity, available as a R/Bioconductor package ( http://bioconductor.org/packages/rnaseqcomp ). Using two independent datasets, we assessed seven competing pipelines. Performance was generally poor, with two methods clearly underperforming and RSEM slightly outperforming the rest.

Genetic heterogeneity contributes to clinical outcome and progression of most tumors, but little is known about allelic diversity for epigenetic compartments, and almost no data exist for acute myeloid leukemia (AML). We examined epigenetic heterogeneity as assessed by cytosine methylation within defined genomic loci with four CpGs (epialleles), somatic mutations, and transcriptomes of AML patient samples at serial time points. We observed that epigenetic allele burden is linked to inferior outcome and varies considerably during disease progression. Epigenetic and genetic allelic burden and patterning followed different patterns and kinetics during disease progression. We observed a subset of AMLs with high epiallele and low somatic mutation burden at diagnosis, a subset with high somatic mutation and lower epiallele burdens at diagnosis, and a subset with a mixed profile, suggesting distinct modes of tumor heterogeneity. Genes linked to promoter-associated epiallele shifts during tumor progression showed increased single-cell transcriptional variance and differential expression, suggesting functional impact on gene regulation. Thus, genetic and epigenetic heterogeneity can occur with distinct kinetics likely to affect the biological and clinical features of tumors.

Cohesin complex members have recently been identified as putative tumor suppressors in hematologic and epithelial malignancies. The cohesin complex guides chromosome segregation; however, cohesin mutant leukemias do not show genomic instability. We hypothesized that reduced cohesin function alters chromatin structure and disrupts cis-regulatory architecture of hematopoietic progenitors. We investigated the consequences of Smc3 deletion in normal and malignant hematopoiesis. Biallelic Smc3 loss induced bone marrow aplasia with premature sister chromatid separation and revealed an absolute requirement for cohesin in hematopoietic stem cell (HSC) function. In contrast, Smc3 haploinsufficiency increased self-renewal in vitro and in vivo, including competitive transplantation. Smc3 haploinsufficiency reduced coordinated transcriptional output, including reduced expression of transcription factors and other genes associated with lineage commitment. Smc3 haploinsufficiency cooperated with Flt3-ITD to induce acute leukemia in vivo, with potentiated Stat5 signaling and altered nucleolar topology. These data establish a dose dependency for cohesin in regulating chromatin structure and HSC function.

Somatic Addition of Sex Combs Like 1 (ASXL1) mutations occur in 10-30% of patients with myeloid malignancies, most commonly in myelodysplastic syndromes (MDSs), and are associated with adverse outcome. Germline ASXL1 mutations occur in patients with Bohring-Opitz syndrome. Here, we show that constitutive loss of Asxl1 results in developmental abnormalities, including anophthalmia, microcephaly, cleft palates, and mandibular malformations. In contrast, hematopoietic-specific deletion of Asxl1 results in progressive, multilineage cytopenias and dysplasia in the context of increased numbers of hematopoietic stem/progenitor cells, characteristic features of human MDS. Serial transplantation of Asxl1-null hematopoietic cells results in a lethal myeloid disorder at a shorter latency than primary Asxl1 knockout (KO) mice. Asxl1 deletion reduces hematopoietic stem cell self-renewal, which is restored by concomitant deletion of Tet2, a gene commonly co-mutated with ASXL1 in MDS patients. Moreover, compound Asxl1/Tet2 deletion results in an MDS phenotype with hastened death compared with single-gene KO mice. Asxl1 loss results in a global reduction of H3K27 trimethylation and dysregulated expression of known regulators of hematopoiesis. RNA-Seq/ChIP-Seq analyses of Asxl1 in hematopoietic cells identify a subset of differentially expressed genes as direct targets of Asxl1. These findings underscore the importance of Asxl1 in Polycomb group function, development, and hematopoiesis.

Phosphatidylinositol-3-kinase (PI3K) signaling is constitutive in most human cancers. Selective inhibition of PI3Kδ (p110δ) by GS-1101 has emerged as a promising therapy in chronic lymphocytic leukemia and indolent lymphomas. In aggressive non-Hodgkin lymphomas such as mantle cell lymphoma (MCL), however, efficacy has been observed, but the extent and duration of tumor control is modest. To determine if tumor killing by GS-1101 is cell cycle-dependent, we show in primary MCL cells by whole-transcriptome sequencing that, despite aberrant expression and recurrent mutations in Cyclin D1, mutations are rare in coding regions of CDK4, RB1 and other genes that control G1-S cell cycle progression or PI3K/AKT signaling. PI3Kδ is the predominant PI3K catalytic subunit expressed, and inhibition by GS-1101 transiently inhibits AKT phosphorylation but not proliferation in MCL cells. Induction of prolonged early G1-arrest (pG1) by selective inhibition of CDK4/CDK6 with PD 0332991 amplifies and sustains PI3Kδ inhibition, which leads to robust apoptosis. Accordingly, inhibition of PI3Kδ induces apoptosis of primary MCL tumor cells once they have ceased to cycle ex vivo, and this killing is enhanced by PD 0332991 inhibition of CDK4/CDK6. PIK3IP1, a negative PI3K regulator, appears to mediate pG1 sensitization to PI3K inhibition; it is markedly reduced in MCL tumor cells compared with normal peripheral B cells, profoundly induced in pG1 and required for pG1 sensitization to GS-1101. Thus, the magnitude and duration of PI3K inhibition and tumor killing by GS-1101 is pG1-dependent, suggesting induction of pG1 by CDK4/CDK6 inhibition as a strategy to sensitize proliferating lymphoma cells to PI3K inhibition.

We have undertaken a genome-wide analysis of rare copy-number variation (CNV) in 1124 autism spectrum disorder (ASD) families, each comprised of a single proband, unaffected parents, and, in most kindreds, an unaffected sibling. We find significant association of ASD with de novo duplications of 7q11.23, where the reciprocal deletion causes Williams-Beuren syndrome, characterized by a highly social personality. We identify rare recurrent de novo CNVs at five additional regions, including 16p13.2 (encompassing genes USP7 and C16orf72) and Cadherin 13, and implement a rigorous approach to evaluating the statistical significance of these observations. Overall, large de novo CNVs, particularly those encompassing multiple genes, confer substantial risks (OR = 5.6; CI = 2.6-12.0, p = 2.4 × 10(-7)). We estimate there are 130-234 ASD-related CNV regions in the human genome and present compelling evidence, based on cumulative data, for association of rare de novo events at 7q11.23, 15q11.2-13.1, 16p11.2, and Neurexin 1.