High-resolution within-sewer SARS-CoV-2 surveillance facilitates informed intervention.

To assist in the COVID-19 public health guidance on a college campus, daily composite wastewater samples were withdrawn at 20 manhole locations across the University of Colorado Boulder campus. Low-cost autosamplers were fabricated in-house to enable an economical approach to this distributed study. These sample stations operated from August 25th until November 23rd during the fall 2020 semester, with 1512 samples collected. The concentration of SARS-CoV-2 in each sample was quantified through two comparative reverse transcription quantitative polymerase chain reactions (RT-qPCRs). These methods were distinct in the utilization of technical replicates and normalization to an endogenous control. (1) Higher temporal resolution compensates for supply chain or other constraints that prevent technical or biological replicates. (2) The data normalized by an endogenous control agreed with the raw concentration data, minimizing the utility of normalization. The raw wastewater concentration values reflected SARS-CoV-2 prevalence on campus as detected by clinical services. Overall, combining the low-cost composite sampler with a method that quantifies the SARS-CoV-2 signal within six hours enabled actionable and time-responsive data delivered to key stakeholders. With daily reporting of the findings, wastewater surveillance assisted in decision making during critical phases of the pandemic on campus, from detecting individual cases within populations ranging from 109 to 2048 individuals to monitoring the success of on-campus interventions.

Higher Viral Load Drives Infrequent Severe Acute Respiratory Syndrome Coronavirus 2 Transmission Between Asymptomatic Residence Hall Roommates.

BACKGROUND: The coronavirus disease 2019 pandemic spread to >200 countries in <6 months. To understand coronavirus spread, determining transmission rate and defining factors that increase transmission risk are essential. Most cases are asymptomatic, but people with asymptomatic infection have viral loads indistinguishable from those in symptomatic people, and they do transmit severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, asymptomatic cases are often undetected. METHODS: Given high residence hall student density, the University of Colorado Boulder established a mandatory weekly screening test program. We analyzed longitudinal data from 6408 students and identified 116 likely transmission events in which a second roommate tested positive within 14 days of the index roommate. RESULTS: Although the infection rate was lower in single-occupancy rooms (10%) than in multiple-occupancy rooms (19%), interroommate transmission occurred only about 20% of the time. Cases were usually asymptomatic at the time of detection. Notably, individuals who likely transmitted had an average viral load approximately 6.5-fold higher than individuals who did not (mean quantification cycle [Cq], 26.2 vs 28.9). Although students with diagnosed SARS-CoV-2 infection moved to isolation rooms, there was no difference in time to isolation between cases with or without interroommate transmission. CONCLUSIONS: This analysis argues that interroommate transmission occurs infrequently in residence halls and provides strong correlative evidence that viral load is proportional to transmission probability.

Just 2% of SARS-CoV-2-positive individuals carry 90% of the virus circulating in communities.

We analyze data from the fall 2020 pandemic response efforts at the University of Colorado Boulder, where more than 72,500 saliva samples were tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using qRT-PCR. All samples were collected from individuals who reported no symptoms associated with COVID-19 on the day of collection. From these, 1,405 positive cases were identified. The distribution of viral loads within these asymptomatic individuals was indistinguishable from what has been previously observed in symptomatic individuals. Regardless of symptomatic status, ∼50% of individuals who test positive for SARS-CoV-2 seem to be in noninfectious phases of the disease, based on having low viral loads in a range from which live virus has rarely been isolated. We find that, at any given time, just 2% of individuals carry 90% of the virions circulating within communities, serving as viral "supercarriers" and possibly also superspreaders.

Saliva TwoStep for rapid detection of asymptomatic SARS-CoV-2 carriers.

Here, we develop a simple molecular test for SARS-CoV-2 in saliva based on reverse transcription loop-mediated isothermal amplification. The test has two steps: (1) heat saliva with a stabilization solution and (2) detect virus by incubating with a primer/enzyme mix. After incubation, saliva samples containing the SARS-CoV-2 genome turn bright yellow. Because this test is pH dependent, it can react falsely to some naturally acidic saliva samples. We report unique saliva stabilization protocols that rendered 295 healthy saliva samples compatible with the test, producing zero false positives. We also evaluated the test on 278 saliva samples from individuals who were infected with SARS-CoV-2 but had no symptoms at the time of saliva collection, and from 54 matched pairs of saliva and anterior nasal samples from infected individuals. The Saliva TwoStep test described herein identified infections with 94% sensitivity and >99% specificity in individuals with sub-clinical (asymptomatic or pre-symptomatic) infections.

Just 2% of SARS-CoV-2-positive individuals carry 90% of the virus circulating in communities.

We analyze data from the Fall 2020 pandemic response efforts at the University of Colorado Boulder (USA), where more than 72,500 saliva samples were tested for SARS-CoV-2 using quantitative RT-PCR. All samples were collected from individuals who reported no symptoms associated with COVID-19 on the day of collection. From these, 1,405 positive cases were identified. The distribution of viral loads within these asymptomatic individuals was indistinguishable from what has been previously reported in symptomatic individuals. Regardless of symptomatic status, approximately 50% of individuals who test positive for SARS-CoV-2 seem to be in non-infectious phases of the disease, based on having low viral loads in a range from which live virus has rarely been isolated. We find that, at any given time, just 2% of individuals carry 90% of the virions circulating within communities, serving as viral "super-carriers" and possibly also super-spreaders.

Saliva TwoStep for rapid detection of asymptomatic SARS-CoV-2 carriers.

Here, we develop a simple molecular test for SARS-CoV-2 in saliva based on reverse transcription loop-mediated isothermal amplification (RT-LAMP). The test has two steps: 1) heat saliva with a stabilization solution, and 2) detect virus by incubating with a primer/enzyme mix. After incubation, saliva samples containing the SARS-CoV-2 genome turn bright yellow. Because this test is pH dependent, it can react falsely to some naturally acidic saliva samples. We report unique saliva stabilization protocols that rendered 295 healthy saliva samples compatible with the test, producing zero false positives. We also evaluated the test on 278 saliva samples from individuals who were infected with SARS-CoV-2 but had no symptoms at the time of saliva collection, and from 54 matched pairs of saliva and anterior nasal samples from infected individuals. The Saliva TwoStep test described herein identified infections with 94% sensitivity and >99% specificity in individuals with sub-clinical (asymptomatic or pre-symptomatic) infections.

Alternative RNA structures formed during transcription depend on elongation rate and modify RNA processing.

Most RNA processing occurs co-transcriptionally. We interrogated nascent pol II transcripts by chemical and enzymatic probing and determined how the "nascent RNA structureome" relates to splicing, A-I editing and transcription speed. RNA folding within introns and steep structural transitions at splice sites are associated with efficient co-transcriptional splicing. A slow pol II mutant elicits extensive remodeling into more folded conformations with increased A-I editing. Introns that become more structured at their 3' splice sites get co-transcriptionally excised more efficiently. Slow pol II altered folding of intronic Alu elements where cryptic splicing and intron retention are stimulated, an outcome mimicked by UV, which decelerates transcription. Slow transcription also remodeled RNA folding around alternative exons in distinct ways that predict whether skipping or inclusion is favored, even though it occurs post-transcriptionally. Hence, co-transcriptional RNA folding modulates post-transcriptional alternative splicing. In summary, the plasticity of nascent transcripts has widespread effects on RNA processing.

Neurodegeneration, Heterochromatin, and Double-Stranded RNA.

Changes in chromatin and epigenetic modifications have been associated with aging and aging-associated neurodegenerative diseases, although the causal relationship between these changes and disease-related pathology has been unclear. Recent studies have now made direct connections between neurodegeneration-associated proteins and derepression of repetitive element transcription due to changes in heterochromatin. We suggest that this derepression leads to an increased accumulation of intracellular double-stranded RNA (dsRNA), with an attendant induction of innate immune responses that contribute to the neuroinflammation found in essentially all age-associated neurodegenerative diseases.

The Caenorhabditis elegans Ortholog of TDP-43 Regulates the Chromatin Localization of the Heterochromatin Protein 1 Homolog HPL-2.

TDP-1 is the Caenorhabditis elegans ortholog of mammalian TDP-43, which is strongly implicated in the etiology of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS). We discovered that deletion of the tdp-1 gene results in enhanced nuclear RNA interference (RNAi). As nuclear RNAi in C. elegans involves chromatin changes moderated by HPL-2, a homolog of heterochromatin protein 1 (HP1), we investigated the interaction of TDP-1 and HPL-2. We found that TDP-1 and HPL-2 interact directly and that loss of TDP-1 dramatically alters the chromatin association of HPL-2. We showed previously that deletion of the tdp-1 gene results in transcriptional alterations and the accumulation of double-stranded RNA (dsRNA). These molecular changes are replicated in an hpl-2 deletion strain, consistent with HPL-2 acting in consort with TDP-1 to modulate these aspects of RNA metabolism. Our observations identify novel mechanisms by which HP1 homologs can be recruited to chromatin and by which nuclear depletion of human TDP-43 may lead to changes in RNA metabolism that are relevant to disease.

Transcription elongation rate affects nascent histone pre-mRNA folding and 3' end processing.

Transcription elongation rate influences cotranscriptional pre-mRNA maturation, but how such kinetic coupling works is poorly understood. The formation of nonadenylated histone mRNA 3' ends requires recognition of an RNA structure by stem-loop-binding protein (SLBP). We report that slow transcription by mutant RNA polymerase II (Pol II) caused accumulation of polyadenylated histone mRNAs that extend past the stem-loop processing site. UV irradiation, which decelerates Pol II elongation, also induced long poly(A)+ histone transcripts. Inhibition of 3' processing by slow Pol II correlates with failure to recruit SLBP to histone genes. Chemical probing of nascent RNA structure showed that the stem-loop fails to fold in transcripts made by slow Pol II, thereby explaining the absence of SLBP and failure to process 3' ends. These results show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

RNA Pol II Dynamics Modulate Co-transcriptional Chromatin Modification, CTD Phosphorylation, and Transcriptional Direction.

Eukaryotic genes are marked by conserved post-translational modifications on the RNA pol II C-terminal domain (CTD) and the chromatin template. How the 5'-3' profiles of these marks are established is poorly understood. Using pol II mutants in human cells, we found that slow transcription repositioned specific co-transcriptionally deposited chromatin modifications; histone H3 lysine 36 trimethyl (H3K36me3) shifted within genes toward 5' ends, and histone H3 lysine 4 dimethyl (H3K4me2) extended farther upstream of start sites. Slow transcription also evoked a hyperphosphorylation of CTD Ser2 residues at 5' ends of genes that is conserved in yeast. We propose a "dwell time in the target zone" model to explain the effects of transcriptional dynamics on the establishment of co-transcriptionally deposited protein modifications. Promoter-proximal Ser2 phosphorylation is associated with a longer pol II dwell time at start sites and reduced transcriptional polarity because of strongly enhanced divergent antisense transcription at promoters. These results demonstrate that pol II dynamics help govern the decision between sense and divergent antisense transcription.

Coupling of RNA Polymerase II Transcription Elongation with Pre-mRNA Splicing.

Pre-mRNA maturation frequently occurs at the same time and place as transcription by RNA polymerase II. The co-transcriptionality of mRNA processing has permitted the evolution of mechanisms that functionally couple transcription elongation with diverse events that occur on the nascent RNA. This review summarizes the current understanding of the relationship between transcriptional elongation through a chromatin template and co-transcriptional splicing including alternative splicing decisions that affect the expression of most human genes.

Localization of RNAPII and 3' end formation factor CstF subunits on C. elegans genes and operons.

Transcription termination is mechanistically coupled to pre-mRNA 3' end formation to prevent transcription much beyond the gene 3' end. C. elegans, however, engages in polycistronic transcription of operons in which 3' end formation between genes is not accompanied by termination. We have performed RNA polymerase II (RNAPII) and CstF ChIP-seq experiments to investigate at a genome-wide level how RNAPII can transcribe through multiple poly-A signals without causing termination. Our data shows that transcription proceeds in some ways as if operons were composed of multiple adjacent single genes. Total RNAPII shows a small peak at the promoter of the gene cluster and a much larger peak at 3' ends. These 3' peaks coincide with maximal phosphorylation of Ser2 within the C-terminal domain (CTD) of RNAPII and maximal localization of the 3' end formation factor CstF. This pattern occurs at all 3' ends including those at internal sites in operons where termination does not occur. Thus the normal mechanism of 3' end formation does not always result in transcription termination. Furthermore, reduction of CstF50 by RNAi did not substantially alter the pattern of CstF64, total RNAPII, or Ser2 phosphorylation at either internal or terminal 3' ends. However, CstF50 RNAi did result in a subtle reduction of CstF64 binding upstream of the site of 3' cleavage, suggesting that the CstF50/CTD interaction may facilitate bringing the 3' end machinery to the transcription complex.

Pre-mRNA splicing is facilitated by an optimal RNA polymerase II elongation rate.

Alternative splicing modulates expression of most human genes. The kinetic model of cotranscriptional splicing suggests that slow elongation expands and that fast elongation compresses the "window of opportunity" for recognition of upstream splice sites, thereby increasing or decreasing inclusion of alternative exons. We tested the model using RNA polymerase II mutants that change average elongation rates genome-wide. Slow and fast elongation affected constitutive and alternative splicing, frequently altering exon inclusion and intron retention in ways not predicted by the model. Cassette exons included by slow and excluded by fast elongation (type I) have weaker splice sites, shorter flanking introns, and distinct sequence motifs relative to "slow-excluded" and "fast-included" exons (type II). Many rate-sensitive exons are misspliced in tumors. Unexpectedly, slow and fast elongation often both increased or both decreased inclusion of a particular exon or retained intron. These results suggest that an optimal rate of transcriptional elongation is required for normal cotranscriptional pre-mRNA splicing.

TDP-1, the Caenorhabditis elegans ortholog of TDP-43, limits the accumulation of double-stranded RNA.

Caenorhabditis elegans mutants deleted for TDP-1, an ortholog of the neurodegeneration-associated RNA-binding protein TDP-43, display only mild phenotypes. Nevertheless, transcriptome sequencing revealed that many RNAs were altered in accumulation and/or processing in the mutant. Analysis of these transcriptional abnormalities demonstrates that a primary function of TDP-1 is to limit formation or stability of double-stranded RNA. Specifically, we found that deletion of tdp-1: (1) preferentially alters the accumulation of RNAs with inherent double-stranded structure (dsRNA); (2) increases the accumulation of nuclear dsRNA foci; (3) enhances the frequency of adenosine-to-inosine RNA editing; and (4) dramatically increases the amount of transcripts immunoprecipitable with a dsRNA-specific antibody, including intronic sequences, RNAs with antisense overlap to another transcript, and transposons. We also show that TDP-43 knockdown in human cells results in accumulation of dsRNA, indicating that suppression of dsRNA is a conserved function of TDP-43 in mammals. Altered accumulation of structured RNA may account for some of the previously described molecular phenotypes (e.g., altered splicing) resulting from reduction of TDP-43 function.

Neurotoxic effects of TDP-43 overexpression in C. elegans.

RNA-binding protein TDP-43 has been associated with multiple neurodegenerative diseases, including amyotrophic lateral sclerosis and frontotemporal lobar dementia. We have engineered pan-neuronal expression of human TDP-43 protein in Caenorhabditis elegans, with the goal of generating a convenient in vivo model of TDP-43 function and neurotoxicity. Transgenic worms with the neuronal expression of human TDP-43 exhibit an 'uncoordinated' phenotype and have abnormal motorneuron synapses. Caenorhabditis elegans contains a single putative ortholog of TDP-43, designated TDP-1, which we show can support alternative splicing of CFTR in a cell-based assay. Neuronal overexpression of TDP-1 also results in an uncoordinated phenotype, while genetic deletion of the tdp-1 gene does not affect movement or alter motorneuron synapses. By using the uncoordinated phenotype as a read-out of TDP-43 overexpression neurotoxicty, we have investigated the contribution of specific TDP-43 domains and subcellular localization to toxicity. Full-length (wild-type) human TDP-43 expressed in C. elegans is localized to the nucleus. Deletion of either RNA recognition domain (RRM1 or RRM2) completely blocks neurotoxicity, as does deletion of the C-terminal region. These deleted TDP-43 variants still accumulate in the nucleus, although their subnuclear distribution is altered. Interestingly, fusion of TDP-1 C-terminal sequences to TDP-43 missing its C-terminal domain restores normal subnuclear localization and toxicity in C. elegans and CFTR splicing in cell-based assays. Overexpression of wild-type, full-length TDP-43 in mammalian cells (differentiated M17 cells) can also result in cell toxicity. Our results demonstrate that in vivo TDP-43 neurotoxicity can result from nuclear activity of overexpressed full-length protein.

Functional redundancy of worm spliceosomal proteins U1A and U2B''.

In Caenorhabditis elegans, the small nuclear ribonucleoprotein (snRNP)-associated proteins U1A and U2B'' are approximately 50% identical to each other, and neither bears signature characteristics of mammalian U1A or U2B'' or the single Drosophila homolog, SNF. We show here that the genes that encode these proteins (rnp-2 and rnp-3) are cotranscribed in an operon, and that ribonucleoprotein RNP-2 is U1 snRNP-associated (U1A) whereas RNP-3 is U2 snRNP-associated (U2B''). U2B'' interacts with U2 even in the absence of another U2 snRNP protein, U2A'. Like U1A and U2B'' from yeast, plants, and vertebrates, worm U1A and U2B'' are more similar to each other than they are to other U1A or U2B'' proteins, respectively. Even though U1A and U2B'' interact with different snRNPs, they are functionally redundant; knockout of both is required for a lethal phenotype. Interestingly, U1A associates with U2 RNA when U2B'' is deleted. Thus, the two members of this gene family normally function as components of different snRNPs but apparently remain capable of performing the function of the other. Redundancy results from the fact that one protein can substitute for the other, even though it normally does not.

Roles of the Homothorax/Meis/Prep homolog UNC-62 and the Exd/Pbx homologs CEH-20 and CEH-40 in C. elegans embryogenesis.

Co-factor homeodomain proteins such as Drosophila Homothorax (Hth) and Extradenticle (Exd) and their respective vertebrate homologs, the Meis/Prep and Pbx proteins, can increase the DNA-binding specificity of Hox protein transcription factors and appear to be required for many of their developmental functions. We show that the unc-62 gene encodes the C. elegans ortholog of Hth, and that maternal-effect unc-62 mutations can cause severe posterior disorganization during embryogenesis (Nob phenotype), superficially similar to that seen in embryos lacking function of either the two posterior-group Hox genes nob-1 and php-3 or the caudal homolog pal-1. Other zygotically acting unc-62 alleles cause earlier embryonic arrest or incompletely penetrant larval lethality with variable morphogenetic defects among the survivors, suggesting that unc-62 functions are required at several stages of development. The differential accumulation of four unc-62 transcripts is consistent with multiple functions. The C. elegans exd homologs ceh-20 and ceh-40 interact genetically with unc-62 and may have overlapping roles in embryogenesis: neither CEH-20 nor CEH-40 appears to be required when the other is present, but loss of both functions causes incompletely penetrant embryonic lethality in the presence of unc-62(+) and complete embryonic lethality in the presence of an unc-62 hypomorphic allele.