Article Methods Cell lines Cell lines were purchased from ATCC and were not formally authenticated, but confirmation of expected gene expression patterns were performed for RNA-seq and eCLIP experiments. Cell lines were routinely tested for mycoplasma contamination (MycoAlert, Lonza).

The highly transmissible B.1.1.7 variant of SARS-CoV-2, first identified in the United Kingdom, has gained a foothold across the world. Using S gene target failure (SGTF) and SARS-CoV-2 genomic sequencing, we investigated the prevalence and dynamics of this variant in the United States (U.S.), tracking it back to its early emergence. We found that while the fraction of B.1.1.7 varied by state, the variant increased at a logistic rate with a roughly weekly doubling rate and an increased transmission of 40-50%. We revealed several independent introductions of B.1.1.7 into the U.S. as early as late November 2020, with community transmission spreading it to most states within months. We show that the U.S. is on a similar trajectory as other countries where B.1.1.7 became dominant, requiring immediate and decisive action to minimize COVID-19 morbidity and mortality.

As SARS-CoV-2 continues to spread and evolve, detecting emerging variants early is critical for public health interventions. Inferring lineage prevalence by clinical testing is infeasible at scale, especially in areas with limited resources, participation, or testing and/or sequencing capacity, which can also introduce biases1–3. SARS-CoV-2 RNA concentration in wastewater successfully tracks regional infection dynamics and provides less biased abundance estimates than clinical testing4,5. Tracking virus genomic sequences in wastewater would improve community prevalence estimates and detect emerging variants. However, two factors limit wastewater-based genomic surveillance: low-quality sequence data and inability to estimate relative lineage abundance in mixed samples. Here we resolve these critical issues to perform a high-resolution, 295-day wastewater and clinical sequencing effort, in the controlled environment of a large university campus and the broader context of the surrounding county. We developed and deployed improved virus concentration protocols and deconvolution software that fully resolve multiple virus strains from wastewater. We detected emerging variants of concern up to 14 days earlier in wastewater samples, and identified multiple instances of virus spread not captured by clinical genomic surveillance. Our study provides a scalable solution for wastewater genomic surveillance that allows early detection of SARS-CoV-2 variants and identification of cryptic transmission.

Direct RNA sequencing holds great promise for the de novo identification of RNA modifications at single-coordinate resolution; however, interpretation of raw sequencing output to discover modified bases remains a challenge. Using Oxford Nanopore's direct RNA sequencing technology, we developed a random forest classifier trained using experimentally detected N 6-methyladenosine (m 6 A) sites within DRACH motifs. Our software MINES (m 6 A Identification using Nanopore Sequencing) assigned m 6 A methylation status to more than 13,000 previously unannotated DRACH sites in endogenous HEK293T transcripts and identified more than 40,000 sites with isoform-level resolution in a human mammary epithelial cell line. These sites displayed sensitivity to the m 6 A writer, METTL3, and eraser, ALKBH5, respectively. MINES (https:// github.com/YeoLab/MINES.git) enables m 6 A annotation at single coordinate-level resolution from direct RNA nanopore sequencing.

SUMMARY Multiple division cycles without growth are a characteristic feature of early embryogenesis. The female germline loads proteins and RNAs into oocytes to support these divisions, which lack many quality control mechanisms operating in somatic cells undergoing growth. Here we describe a small RNA-Argonaute pathway that ensures early embryonic divisions in C. elegans by employing catalytic slicing activity to broadly tune, instead of silence, germline gene expression. Misregulation of one target, a kinesin-13 microtubule depolymerase, underlies a major phenotype associated with pathway loss. Tuning of target transcript levels is guided by density of homologous small RNAs, whose generation must ultimately be related to target sequence. Thus, the tuning action of a small RNA-catalytic Argonaute pathway generates oocytes capable of supporting embryogenesis. We speculate that the specialized nature of germline chromatin led to emergence of small RNA-catalytic Argonaute pathways in the female germline as a post-transcriptional control layer to optimize oocyte composition.

The N6-methyladenosine (m6A) modification is the most prevalent post-transcriptional mRNA modification, regulating mRNA decay and splicing. It plays a major role during normal development, differentiation, and disease progression. The modification is regulated by a set of writer, eraser, and reader proteins. The YTH domain family of proteins consists of three homologous m6A-binding proteins, Ythdf1, Ythdf2, and Ythdf3, which were suggested to have different cellular functions. However, their sequence similarity and their tendency to bind the same targets suggest that they may have overlapping roles. We systematically knocked out (KO) the Mettl3 writer, each of the Ythdf readers, and the three readers together (triple-KO). We then estimated the effect in vivo in mouse gametogenesis, postnatal viability, and in vitro in mouse embryonic stem cells (mESCs). In gametogenesis, Mettl3-KO severity is increased as the deletion occurs earlier in the process, and Ythdf2 has a dominant role that cannot be compensated by Ythdf1 or Ythdf3, due to differences in readers’ expression pattern across different cell types, both in quantity and in spatial location. Knocking out the three readers together and systematically testing viable offspring genotypes revealed a redundancy in the readers’ role during early development that is Ythdf1/2/3 gene dosage-dependent. Finally, in mESCs there is compensation between the three Ythdf reader proteins, since the resistance to differentiate and the significant effect on mRNA decay occur only in the triple-KO cells and not in the single KOs. Thus, we suggest a new model for the Ythdf readers function, in which there is profound dosage-dependent redundancy when all three readers are equivalently coexpressed in the same cell types.

Deadenylases are best known for degrading the poly(A) tail during mRNA decay. The deadenylase family has expanded throughout evolution and, in mammals, consists of 12 Mg2+-dependent 3’ end ribonucleases with mostly unknown substrate specificity1. Pontocerebellar hypoplasia type 7 (PCH7) is a unique recessive syndrome characterized by neurodegeneration with ambiguous genitalia2 (MIM%614969). We studied 12 human families with PCH7, uncovering biallelic, loss of function mutations in TOE1 (NC_000001.11), which encodes an unconventional deadenylase3,4. Toe1-morphant zebrafish displayed mid- and hind-brain degeneration, modeling PCH-like structural defects in vivo. Surprisingly, we found TOE1 associated with incompletely processed small nuclear (sn)RNAs of the spliceosome, which is responsible for pre-mRNA splicing. These pre-snRNAs contained 3’ genome-encoded tails often followed by post-transcriptionally added adenosines. Human cells with reduced levels of TOE1 accumulated 3’ end-extended pre-snRNAs, and immuno-isolated TOE1 complex was sufficient for 3’ end maturation of snRNAs. Our findings reveal the cause of a neurodegenerative syndrome linked to snRNA maturation and uncover a key factor involved in processing of snRNA 3’ ends.

A metagenomic analysis was performed on a soil profile from a wet tundra site in northern Alaska. The goal was to link existing biogeochemical knowledge of the system with the organisms and genes responsible for the relevant metabolic pathways. We specifically investigated how the importance of iron (Fe) oxides and humic substances (HS) as terminal electron acceptors in this ecosystem is expressed genetically, and how respiratory and fermentative processes varied with soil depth into the active layer and into the upper permafrost. Overall, the metagenomes reflected a microbial community enriched in a diverse range of anaerobic pathways, with a preponderance of known Fe reducing species at all depths in the profile. The abundance of sequences associated with anaerobic metabolic processes generally increased with depth, while aerobic cytochrome c oxidases decreased. Methanogenesis genes and methanogen genomes followed the pattern of CH4 fluxes : they increased steeply with depth into the active layer, but declined somewhat over the transition zone between the lower active layer and the upper permafrost. The latter was relatively enriched in fermentative and anaerobic respiratory pathways. A survey of decaheme cytochromes (MtrA, MtrC and their homologs) revealed that this is a promising approach to identifying potential reducers of Fe(III) or HS, and indicated a possible role for Acidobacteria as Fe reducers in these soils. Methanogens appear to coexist in the same layers, though in lower abundance, with Fe reducing bacteria and other potential competitors, including acetogens. These observations provide a rich set of hypotheses for further targeted study.