Systematic assessment of long-read RNA-seq methods for transcript identification and quantification.

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

Measuring the buffering capacity of gene silencing in Saccharomyces cerevisiae.

Gene silencing in budding yeast is mediated by Sir protein binding to unacetylated nucleosomes to form a chromatin structure that inhibits transcription. Transcriptional silencing is characterized by the high-fidelity transmission of the silent state. Despite its relative stability, the constituent parts of the silent state are in constant flux, giving rise to a model that silent loci can tolerate such fluctuations without functional consequences. However, the level of tolerance is unknown, and we developed methods to measure the threshold of histone acetylation that causes the silent chromatin state to switch to the active state as well as to measure the levels of the enzymes and structural proteins necessary for silencing. We show that loss of silencing required 50 to 75% acetyl-mimic histones, though the precise levels were influenced by silencer strength and upstream activating sequence (UAS) enhancer/promoter strength. Measurements of repressor protein levels necessary for silencing showed that reducing SIR4 gene dosage two- to threefold significantly weakened silencing, though reducing the gene copy numbers for Sir2 or Sir3 to the same extent did not significantly affect silencing suggesting that Sir4 was a limiting component in gene silencing. Calculations suggest that a mere twofold reduction in the ability of acetyltransferases to acetylate nucleosomes across a large array of nucleosomes may be sufficient to generate a transcriptionally silent domain.

Transcriptional silencing in Saccharomyces cerevisiae: known unknowns.

Transcriptional silencing in Saccharomyces cerevisiae is a persistent and highly stable form of gene repression. It involves DNA silencers and repressor proteins that bind nucleosomes. The silenced state is influenced by numerous factors including the concentration of repressors, nature of activators, architecture of regulatory elements, modifying enzymes and the dynamics of chromatin.Silencers function to increase the residence time of repressor Sir proteins at silenced domains while clustering of silenced domains enables increased concentrations of repressors and helps facilitate long-range interactions. The presence of an accessible NDR at the regulatory regions of silenced genes, the cycling of chromatin configurations at regulatory sites, the mobility of Sir proteins, and the non-uniform distribution of the Sir proteins across the silenced domain, all result in silenced chromatin that only stably silences weak promoters and enhancers via changes in transcription burst duration and frequency.These data collectively suggest that silencing is probabilistic and the robustness of silencing is achieved through sub-optimization of many different nodes of action such that a stable expression state is generated and maintained even though individual constituents are in constant flux.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration.

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Probing chromatin accessibility with small molecule DNA intercalation and nanopore sequencing.

Genome-wide identification of chromatin organization and structure has been generally probed by measuring accessibility of the underlying DNA to nucleases or methyltransferases. These methods either only observe the positioning of a single nucleosome or rely on large enzymes to modify or cleave the DNA. We developed adduct sequencing (Add-seq), a method to probe chromatin accessibility by treating chromatin with the small molecule angelicin, which preferentially intercalates into DNA not bound to core nucleosomes. We show that Nanopore sequencing of the angelicin-modified DNA is possible and allows visualization and analysis of long single molecules with distinct chromatin structure. The angelicin modification can be detected from the Nanopore current signal data using a neural network model trained on unmodified and modified chromatin-free DNA. Applying Add-seq to Saccharomyces cerevisiae nuclei, we identified expected patterns of accessibility around annotated gene loci in yeast. We also identify individual clusters of single molecule reads displaying different chromatin structure at specific yeast loci, which demonstrates heterogeneity in the chromatin structure of the yeast population. Thus, using Add-seq, we are able to profile DNA accessibility in the yeast genome across long molecules.

DNA polymerase epsilon, acetylases and remodellers cooperate to form a specialized chromatin structure at a tRNA insulator.

Insulators bind transcription factors and use chromatin remodellers and modifiers to mediate insulation. In this report, we identified proteins required for the efficient formation and maintenance of a specialized chromatin structure at the yeast tRNA insulator. The histone acetylases, SAS-I and NuA4, functioned in insulation, independently of tRNA and did not participate in the formation of the hypersensitive site at the tRNA. In contrast, DNA polymerase epsilon, functioned with the chromatin remodeller, Rsc, and the histone acetylase, Rtt109, to generate a histone-depleted region at the tRNA insulator. Rsc and Rtt109 were required for efficient binding of TFIIIB to the tRNA insulator, and the bound transcription factor and Rtt109 in turn were required for the binding of Rsc to tRNA. Robust insulation during growth and cell division involves the formation of a hypersensitive site at the insulator during chromatin maturation together with competition between acetylases and deacetylases.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration.

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification.

The Long-read RNA-Seq Genome Annotation Assessment Project (LRGASP) Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. The consortium generated over 427 million long-read sequences from cDNA and direct RNA datasets, encompassing human, mouse, and manatee species, using different protocols and sequencing platforms. These data were utilized by developers to address challenges in transcript isoform detection and quantification, as well as de novo transcript isoform identification. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. When aiming to detect rare and novel transcripts or when using reference-free approaches, incorporating additional orthogonal data and replicate samples are advised. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

A new SARS-CoV-2 lineage that shares mutations with known Variants of Concern is rejected by automated sequence repository quality control.

We report a SARS-CoV-2 lineage that shares N501Y, P681H, and other mutations with known variants of concern, such as B.1.1.7. This lineage, which we refer to as B.1.x (COG-UK sometimes references similar samples as B.1.324.1), is present in at least 20 states across the USA and in at least six countries. However, a large deletion causes the sequence to be automatically rejected from repositories, suggesting that the frequency of this new lineage is underestimated using public data. Recent dynamics based on 339 samples obtained in Santa Cruz County, CA, USA suggest that B.1.x may be increasing in frequency at a rate similar to that of B.1.1.7 in Southern California. At present the functional differences between this variant B.1.x and other circulating SARS-CoV-2 variants are unknown, and further studies on secondary attack rates, viral loads, immune evasion and/or disease severity are needed to determine if it poses a public health concern. Nonetheless, given what is known from well-studied circulating variants of concern, it seems unlikely that the lineage could pose larger concerns for human health than many already globally distributed lineages. Our work highlights a need for rapid turnaround time from sequence generation to submission and improved sequence quality control that removes submission bias. We identify promising paths toward this goal.

A method for campus-wide SARS-CoV-2 surveillance at a large public university.

The systematic screening of asymptomatic and pre-symptomatic individuals is a powerful tool for controlling community transmission of infectious disease on college campuses. Faced with a paucity of testing in the beginning of the COVID-19 pandemic, many universities developed molecular diagnostic laboratories focused on SARS-CoV-2 diagnostic testing on campus and in their broader communities. We established the UC Santa Cruz Molecular Diagnostic Lab in early April 2020 and began testing clinical samples just five weeks later. Using a clinically-validated laboratory developed test (LDT) that avoided supply chain constraints, an automated sample pooling and processing workflow, and a custom laboratory information management system (LIMS), we expanded testing from a handful of clinical samples per day to thousands per day with the testing capacity to screen our entire campus population twice per week. In this report we describe the technical, logistical, and regulatory processes that enabled our pop-up lab to scale testing and reporting capacity to thousands of tests per day.

Permutational analysis of Saccharomyces cerevisiae regulatory elements.

Gene expression in Saccharomyces cerevisiae is regulated at multiple levels. Genomic and epigenomic mapping of transcription factors and chromatin factors has led to the delineation of various modular regulatory elements-enhancers (upstream activating sequences), core promoters, 5' untranslated regions (5' UTRs) and transcription terminators/3' untranslated regions (3' UTRs). However, only a few of these elements have been tested in combinations with other elements and the functional interactions between the different modular regulatory elements remain under explored. We describe a simple and rapid approach to build a combinatorial library of regulatory elements and have used this library to study 26 different enhancers, core promoters, 5' UTRs and transcription terminators/3' UTRs to estimate the contribution of individual regulatory parts in gene expression. Our combinatorial analysis shows that while enhancers initiate gene expression, core promoters modulate the levels of enhancer-mediated expression and can positively or negatively affect expression from even the strongest enhancers. Principal component analysis (PCA) indicates that enhancer and promoter function can be explained by a single principal component while UTR function involves multiple functional components. The PCA also highlights outliers and suggest differences in mechanisms of regulation by individual elements. Our data also identify numerous regulatory cassettes composed of different individual regulatory elements that exhibit equivalent gene expression levels. These data thus provide a catalog of elements that could in future be used in the design of synthetic regulatory circuits.

H2A.Z functions to regulate progression through the cell cycle.

Histone H2A variants are highly conserved proteins found ubiquitously in nature and thought to perform specialized functions in the cell. Studies in yeast on the histone H2A variant H2A.Z have shown a role for this protein in transcription as well as chromosome segregation. Our studies have focused on understanding the role of H2A.Z during cell cycle progression. We found that htz1delta cells were delayed in DNA replication and progression through the cell cycle. Furthermore, cells lacking H2A.Z required the S-phase checkpoint pathway for survival. We also found that H2A.Z localized to the promoters of cyclin genes, and cells lacking H2A.Z were delayed in the induction of these cyclin genes. Several different models are proposed to explain these observations.

tRNA Genes Affect Chromosome Structure and Function via Local Effects.

The genome is packaged and organized in an ordered, nonrandom manner, and specific chromatin segments contact nuclear substructures to mediate this organization. tRNA genes (tDNAs) are binding sites for transcription factors and architectural proteins and are thought to play an important role in the organization of the genome. In this study, we investigate the roles of tDNAs in genomic organization and chromosome function by editing a chromosome so that it lacked any tDNAs. Surprisingly our analyses of this tDNA-less chromosome show that loss of tDNAs does not grossly affect chromatin architecture or chromosome tethering and mobility. However, loss of tDNAs affects local nucleosome positioning and the binding of SMC proteins at these loci. The absence of tDNAs also leads to changes in centromere clustering and a reduction in the frequency of long-range HML-HMR heterochromatin clustering with concomitant effects on gene silencing. We propose that the tDNAs primarily affect local chromatin structure, which results in effects on long-range chromosome architecture.

Breaking through to the other side: silencers and barriers.

The establishment and restriction of transcriptionally inactive regions in the nucleus is mediated by silencer and barrier elements. Silencer-bound proteins recruit additional factors to establish the silenced domain during the S-phase of the cell cycle but, contrary to previous models, DNA replication is not a pre-requisite for the establishment. Characteristically, silenced domains contain hypoacetylated histones and recent data have identified residue-specific methylation of histone H3 as an additional signature that distinguishes active regions from inactive ones. Peaks of acetylated histones demarcate the boundaries between these regions and recruitment of HAT activities provides a mechanism to restrict the spread of heterochromatin.

Heterochromatin formation via recruitment of DNA repair proteins.

Heterochromatin formation and nuclear organization are important in gene regulation and genome fidelity. Proteins involved in gene silencing localize to sites of damage and some DNA repair proteins localize to heterochromatin, but the biological importance of these correlations remains unclear. In this study, we examined the role of double-strand-break repair proteins in gene silencing and nuclear organization. We find that the ATM kinase Tel1 and the proteins Mre11 and Esc2 can silence a reporter gene dependent on the Sir, as well as on other repair proteins. Furthermore, these proteins aid in the localization of silenced domains to specific compartments in the nucleus. We identify two distinct mechanisms for repair protein-mediated silencing-via direct and indirect interactions with Sir proteins, as well as by tethering loci to the nuclear periphery. This study reveals previously unknown interactions between repair proteins and silencing proteins and suggests insights into the mechanism underlying genome integrity.

Transcription independent insulation at TFIIIC-dependent insulators.

Chromatin insulators separate active from repressed chromatin domains. In yeast the RNA pol III transcription machinery bound to tRNA genes function with histone acetylases and chromatin remodelers to restrict the spread of heterochromatin. Our results collectively demonstrate that binding of TFIIIC is necessary for insulation but binding of TFIIIB along with TFIIIC likely improves the probability of complex formation at an insulator. Insulation by this transcription factor occurs in the absence of RNA polymerase III or polymerase II but requires specific histone acetylases and chromatin remodelers. This analysis identifies a minimal set of factors required for insulation.

Long-range communication between the silencers of HMR.

Gene regulation involves long-range communication between silencers, enhancers, and promoters. In Saccharomyces cerevisiae, silencers flank transcriptionally repressed genes to mediate regional silencing. Silencers recruit the Sir proteins, which then spread along chromatin to encompass the entire silenced domain. In this report we have employed a boundary trap assay, an enhancer activity assay, chromatin immunoprecipitations, and chromosome conformation capture analyses to demonstrate that the two HMR silencer elements are in close proximity and functionally communicate with one another in vivo. We further show that silencing is necessary for these long-range interactions, and we present models for Sir-mediated silencing based upon these results.