Search

RNA-sequencing (RNA-seq) is a gold-standard method to profile genome-wide changes in gene expression. RNA-seq uses high-throughput sequencing technology to quantify the amount of RNA in a biological sample. With the increasing popularity of RNA-seq, many variations on the protocol have been proposed to extract unique and relevant information from biological samples. 3' Tag-Seq (also called TagSeq, 3' Tag-RNA-Seq, and Quant-Seq 3' mRNA-Seq) is one RNA-seq variation where the 3' end of the transcript is selected and amplified to yield one copy of cDNA from each transcript in the biological sample. We present a simple, easy-to-use, and publicly available computational workflow to analyze 3' Tag-Seq data. The workflow begins by trimming sequence adapters from raw FASTQ files. The trimmed sequence reads are checked for quality using FastQC and aligned to the reference genome, and then read counts are obtained using STAR. Differential gene expression analysis is performed using DESeq2, based on differential analysis of gene count data. The outputs of this workflow are MA plots, tables of differentially expressed genes, and UpSet plots. This protocol is intended for users specifically interested in analyzing 3' Tag-Seq data, and thus normalizations based on transcript length are not performed within the workflow. Future updates to this workflow could include custom analyses based on the gene counts table as well as data visualization enhancements. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Running the 3' Tag-Seq workflow Support Protocol: Generating genome indices.

In response to the COVID-19 pandemic, many higher educational institutions moved their courses on-line in hopes of slowing disease spread. The advent of multiple highly-effective vaccines offers the promise of a return to "normal" in-person operations, but it is not clear if-or for how long-campuses should employ non-pharmaceutical interventions such as requiring masks or capping the size of in-person courses. In this study, we develop and fine-tune a model of COVID-19 spread to UC Merced's student and faculty population. We perform a global sensitivity analysis to consider how both pharmaceutical and non-pharmaceutical interventions impact disease spread. Our work reveals that vaccines alone may not be sufficient to eradicate disease dynamics and that significant contact with an infectious surrounding community will maintain infections on-campus. Our work provides a foundation for higher-education planning allowing campuses to balance the benefits of in-person instruction with the ability to quarantine/isolate infectious individuals.

Genetic regulatory networks (GRNs) regulate the flow of genetic information from the genome to expressed messenger RNAs (mRNAs) and thus are critical to controlling the phenotypic characteristics of cells. Numerous methods exist for profiling mRNA transcript levels and identifying protein-DNA binding interactions at the genome-wide scale. These enable researchers to determine the structure and output of transcriptional regulatory networks, but uncovering the complete structure and regulatory logic of GRNs remains a challenge. The field of GRN inference aims to meet this challenge using computational modeling to derive the structure and logic of GRNs from experimental data and to encode this knowledge in Boolean networks, Bayesian networks, ordinary differential equation (ODE) models, or other modeling frameworks. However, most existing models do not incorporate dynamic transcriptional data since it has historically been less widely available in comparison to "static" transcriptional data. We report the development of an evolutionary algorithm-based ODE modeling approach (named EA) that integrates kinetic transcription data and the theory of attractor matching to infer GRN architecture and regulatory logic. Our method outperformed six leading GRN inference methods, none of which incorporate kinetic transcriptional data, in predicting regulatory connections among TFs when applied to a small-scale engineered synthetic GRN in Saccharomyces cerevisiae. Moreover, we demonstrate the potential of our method to predict unknown transcriptional profiles that would be produced upon genetic perturbation of the GRN governing a two-state cellular phenotypic switch in Candida albicans. We established an iterative refinement strategy to facilitate candidate selection for experimentation; the experimental results in turn provide validation or improvement for the model. In this way, our GRN inference approach can expedite the development of a sophisticated mathematical model that can accurately describe the structure and dynamics of the in vivo GRN.

Regulatory transcription factors control many important biological processes, including cellular differentiation, responses to environmental perturbations and stresses, and host-pathogen interactions. Determining the genome-wide binding of regulatory transcription factors to DNA is essential to understanding the function of transcription factors in these often complex biological processes. Cleavage under targets and release using nuclease (CUT&RUN) is a modern method for genome-wide mapping of in vivo protein-DNA binding interactions that is an attractive alternative to the traditional and widely used chromatin immunoprecipitation followed by sequencing (ChIP-seq) method. CUT&RUN is amenable to a higher-throughput experimental setup and has a substantially higher dynamic range with lower per-sample sequencing costs than ChIP-seq. Here, a comprehensive CUT&RUN protocol and accompanying data analysis workflow tailored for genome-wide analysis of transcription factor-DNA binding interactions in the human fungal pathogen Candida albicans are described. This detailed protocol includes all necessary experimental procedures, from epitope tagging of transcription factor-coding genes to library preparation for sequencing; additionally, it includes a customized computational workflow for CUT&RUN data analysis.

CRISPR/Cas-induced genome editing is a powerful tool for genetic engineering, however, targeting constraints limit which loci are editable with this method. Since the length of a DNA sequence impacts the likelihood it overlaps a unique target site, precision editing of small genomic features with CRISPR/Cas remains an obstacle. We introduce a two-step genome editing strategy that virtually eliminates CRISPR/Cas targeting constraints and facilitates precision genome editing of elements as short as a single base-pair at virtually any locus in any organism that supports CRISPR/Cas-induced genome editing. Our two-step approach first replaces the locus of interest with an "AddTag" sequence, which is subsequently replaced with any engineered sequence, and thus circumvents the need for direct overlap with a unique CRISPR/Cas target site. In this study, we demonstrate the feasibility of our approach by editing transcription factor binding sites within Candida albicans that could not be targeted directly using the traditional gene-editing approach. We also demonstrate the utility of the AddTag approach for combinatorial genome editing and gene complementation analysis, and we present a software package that automates the design of AddTag editing.

Amyloid appearance is a rare event that is promoted in the presence of other aggregated proteins. These aggregates were thought to act by templating the formation of an assembly-competent nucleation seed, but we find an unanticipated role for them in enhancing the persistence of amyloid after it arises. Specifically, Saccharomyces cerevisiae Rnq1 amyloid reduces chaperone-mediated disassembly of Sup35 amyloid, promoting its persistence in yeast. Mathematical modeling and corresponding in vivo experiments link amyloid persistence to the conformationally defined size of the Sup35 nucleation seed and suggest that amyloid is actively cleared by disassembly below this threshold to suppress appearance of the [PSI+] prion in vivo. Remarkably, this framework resolves multiple known inconsistencies in the appearance and curing of yeast prions. Thus, our observations establish the size of the nucleation seed as a previously unappreciated characteristic of prion variants that is key to understanding transitions between prion states.

Prions are proteins capable of adopting misfolded conformations and transmitting these conformations to other normally folded proteins. Prions are most commonly known for causing fatal neurodegenerative diseases in mammals but are also associated with several harmless phenotypes in yeast. A distinct feature of prion propagation is the existence of different phenotypical variants, called strains. It is widely accepted that these strains correspond to different conformational states of the protein, but the mechanisms driving their interactions remain poorly understood. This study uses mathematical modeling to provide insight into this problem. We show that the classical model of prion dynamics allows at most one conformational strain to stably propagate. In order to conform to biological observations of strain coexistence and co-stability, we develop an extension of the classical model by introducing a novel prion species consistent with biological studies. Qualitative analysis of this model reveals a new variety of behavior. Indeed, it allows for stable coexistence of different strains in a wide parameter range, and it also introduces intricate initial condition dependency. These new behaviors are consistent with experimental observations of prions in both mammals and yeast. As such, our model provides a valuable tool for investigating the underlying mechanisms of prion propagation and the link between prion strains and strain specific phenotypes. The consideration of a novel prion species brings a change in perspective on prion biology and we use our model to generate hypotheses about prion infectivity.

CD8 T cells can play both a protective and pathogenic role in inflammation and autoimmune development. Recent studies have highlighted the ability of CD8 T cells to function as T follicular helper (Tfh) cells in the germinal center in the context of infection. However, whether this phenomenon occurs in autoimmunity and contributes to autoimmune pathogenesis is largely unexplored. In this study, we show that CD8 T cells acquire a CD4 Tfh profile in the absence of functional regulatory T cells in both the IL-2-deficient and scurfy mouse models. Depletion of CD8 T cells mitigates autoimmune pathogenesis in IL-2-deficient mice. CD8 T cells express the B cell follicle-localizing chemokine receptor CXCR5, a principal Tfh transcription factor Bcl6, and the Tfh effector cytokine IL-21. CD8 T cells localize to the B cell follicle, express B cell costimulatory proteins, and promote B cell differentiation and Ab isotype class switching. These data reveal a novel contribution of autoreactive CD8 T cells to autoimmune disease, in part, through CD4 follicular-like differentiation and functionality.

Prions adopt alternative, self-replicating protein conformations and thereby determine novel phenotypes that are often irreversible. Nevertheless, dominant-negative prion mutants can revert phenotypes associated with some conformations. These observations suggest that, while intervention is possible, distinct inhibitors must be developed to overcome the conformational plasticity of prions. To understand the basis of this specificity, we determined the impact of the G58D mutant of the Sup35 prion on three of its conformational variants, which form amyloids in S. cerevisiae. G58D had been previously proposed to have unique effects on these variants, but our studies suggest a common mechanism. All variants, including those reported to be resistant, are inhibited by G58D but at distinct doses. G58D lowers the kinetic stability of the associated amyloid, enhancing its fragmentation by molecular chaperones, promoting Sup35 resolubilization, and leading to amyloid clearance particularly in daughter cells. Reducing the availability or activity of the chaperone Hsp104, even transiently, reverses curing. Thus, the specificity of inhibition is determined by the sensitivity of variants to the mutant dosage rather than mode of action, challenging the view that a unique inhibitor must be developed to combat each variant.

Prions are proteins most commonly associated with fatal neurodegenerative diseases in mammals but are also responsible for a number of harmless heritable phenotypes in yeast. These states arise when a misfolded form of a protein appears and, rather than be removed by cellular quality control mechanisms, persists. The misfolded prion protein forms aggregates and is capable of converting normally folded protein to the misfolded state through direct interaction between the two forms. The dominant mathematical model for prion aggregate dynamics has been the nucleated polymerization model (NPM) which considers the dynamics of only the normal protein and the aggregates. However, for yeast prions the molecular chaperone Hsp104 is essential for prion propagation. Further, although mammals do not express Hsp104, experimental assays have shown Hsp104 also interacts with mammalian prion aggregates. In this study, we generalize the NPM to account for molecular chaperones and develop what we call the enzyme-limited nucleated polymerization model (ELNPM). We discuss existence, uniqueness and stability of solutions to our model and demonstrate that the NPM represents a quasi-steady-state reduction of our model. We validate the ELNPM by demonstrating agreement with experimental results on the yeast prion PSI(+) that could not be supported by the NPM. Finally, we demonstrate that, in contrast to the NPM, the ELNPM permits the coexistence of multiple prion strains.

The nucleated polymerization model is a mathematical framework that has been applied to aggregation and fragmentation processes in both the discrete and continuous setting. In particular, this model has been the canonical framework for analyzing the dynamics of protein aggregates arising in prion and amyloid diseases such as as Alzheimer's and Parkinson's disease. We present an explicit steady-state solution to the aggregate size distribution governed by the discrete nucleated polymerization equations. Steady-state solutions have been previously obtained under the assumption of continuous aggregate sizes; however, the discrete solution allows for direct computation and parameter inference, as well as facilitates estimates on the accuracy of the continuous approximation.

The self-assembly of alternative conformations of normal proteins into amyloid aggregates has been implicated in both the acquisition of new functions and in the appearance and progression of disease. However, while these amyloidogenic pathways are linked to the emergence of new phenotypes, numerous studies have uncoupled the accumulation of aggregates from their biological consequences, revealing currently underappreciated complexity in the determination of these traits. Here, to explore the molecular basis of protein-only phenotypes, we focused on the Saccharomyces cerevisiae Sup35/[PSI(+)] prion, which confers a translation termination defect and expression level-dependent toxicity in its amyloid form. Our studies reveal that aggregated Sup35 retains its normal function as a translation release factor. However, fluctuations in the composition and size of these complexes specifically alter the level of this aggregate-associated activity and thereby the severity and toxicity of the amyloid state. Thus, amyloid heterogeneity is a crucial contributor to protein-only phenotypes.

Abstract Background Polymorphic inversions are a source of genetic variability with a direct impact on recombination frequencies. Given the difficulty of their experimental study, computational methods have been developed to infer their existence in a large number of individuals using genome-wide data of nucleotide variation. Methods based on haplotype tagging of known inversions attempt to classify individuals as having a normal or inverted allele. Other methods that measure differences between linkage disequilibrium attempt to identify regions with inversions but unable to classify subjects accurately, an essential requirement for association studies. Results We present a novel method to both identify polymorphic inversions from genome-wide genotype data and classify individuals as containing a normal or inverted allele. Our method, a generalization of a published method for haplotype data 1, utilizes linkage between groups of SNPs to partition a set of individuals into normal and inverted subpopulations. We employ a sliding window scan to identify regions likely to have an inversion, and accumulation of evidence from neighboring SNPs is used to accurately determine the inversion status of each subject. Further, our approach detects inversions directly from genotype data, thus increasing its usability to current genome-wide association studies (GWAS). Conclusions We demonstrate the accuracy of our method to detect inversions and classify individuals on principled-simulated genotypes, produced by the evolution of an inversion event within a coalescent model 2. We applied our method to real genotype data from HapMap Phase III to characterize the inversion status of two known inversions within the regions 17q21 and 8p23 across 1184 individuals. Finally, we scan the full genomes of the European Origin (CEU) and Yoruba (YRI) HapMap samples. We find population-based evidence for 9 out of 15 well-established autosomic inversions, and for 52 regions previously predicted by independent experimental methods in ten (9+1) individuals 34. We provide efficient implementations of both genotype and haplotype methods as a unified R package inveRsion.