Transcriptome variation in human tissues revealed by long-read sequencing.

Regulation of transcript structure generates transcript diversity and plays an important role in human disease1-7. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure8-16. In this Article, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from Genotype-Tissue Expression (GTEx) tissues and cell lines, complementing the GTEx resource. We identified just over 70,000 novel transcripts for annotated genes, and validated the protein expression of 10% of novel transcripts. We developed a new computational package, LORALS, to analyse the genetic effects of rare and common variants on the transcriptome by allele-specific analysis of long reads. We characterized allele-specific expression and transcript structure events, providing new insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb the transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we used this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.

An autoimmune disease risk variant: A trans master regulatory effect mediated by IRF1 under immune stimulation?

Functional mechanisms remain unknown for most genetic loci associated to complex human traits and diseases. In this study, we first mapped trans-eQTLs in a data set of primary monocytes stimulated with LPS, and discovered that a risk variant for autoimmune disease, rs17622517 in an intron of C5ORF56, affects the expression of the transcription factor IRF1 20 kb away. The cis-regulatory effect specific to IRF1 is active under early immune stimulus, with a large number of trans-eQTL effects across the genome under late LPS response. Using CRISPRi silencing, we showed that perturbation of the SNP locus downregulates IRF1 and causes widespread transcriptional effects. Genome editing by CRISPR had suggestive recapitulation of the LPS-specific trans-eQTL signal and lent support for the rs17622517 site being functional. Our results suggest that this common genetic variant affects inter-individual response to immune stimuli via regulation of IRF1. For this autoimmune GWAS locus, our work provides evidence of the functional variant, demonstrates a condition-specific enhancer effect, identifies IRF1 as the likely causal gene in cis, and indicates that overactivation of the downstream immune-related pathway may be the cellular mechanism increasing disease risk. This work not only provides rare experimental validation of a master-regulatory trans-eQTL, but also demonstrates the power of eQTL mapping to build mechanistic hypotheses amenable for experimental follow-up using the CRISPR toolkit.

Evaluation of noninvasive biospecimens for transcriptome studies.

Transcriptome studies disentangle functional mechanisms of gene expression regulation and may elucidate the underlying biology of disease processes. However, the types of tissues currently collected typically assay a single post-mortem timepoint or are limited to investigating cell types found in blood. Noninvasive tissues may improve disease-relevant discovery by enabling more complex longitudinal study designs, by capturing different and potentially more applicable cell types, and by increasing sample sizes due to reduced collection costs and possible higher enrollment from vulnerable populations. Here, we develop methods for sampling noninvasive biospecimens, investigate their performance across commercial and in-house library preparations, characterize their biology, and assess the feasibility of using noninvasive tissues in a multitude of transcriptomic applications. We collected buccal swabs, hair follicles, saliva, and urine cell pellets from 19 individuals over three to four timepoints, for a total of 300 unique biological samples, which we then prepared with replicates across three library preparations, for a final tally of 472 transcriptomes. Of the four tissues we studied, we found hair follicles and urine cell pellets to be most promising due to the consistency of sample quality, the cell types and expression profiles we observed, and their performance in disease-relevant applications. This is the first study to thoroughly delineate biological and technical features of noninvasive samples and demonstrate their use in a wide array of transcriptomic and clinical analyses. We anticipate future use of these biospecimens will facilitate discovery and development of clinical applications.

SARS-CoV-2 Sequence Analysis during COVID-19 Case Surge, Liberia, 2021.

In June 2021, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases surged in Liberia. SARS-CoV-2 sequences from patients hospitalized during March-July 2021 revealed the Delta variant was in Liberia in early March and was dominant in June, irrespective of geography. Mutations and deletions suggest multiple SARS-CoV-2 Delta variant introductions.

Capture Sequencing Enables Sensitive Detection of Tick-Borne Agents in Human Blood.

Assay sensitivity can be a limiting factor in the use of PCR as a tool for the detection of tick-borne pathogens in blood. We evaluated the performance of Tick-borne disease Capture Sequencing Assay (TBDCapSeq), a capture sequencing assay targeting tick-borne agents, to test 158 whole blood specimens obtained from the Lyme Disease Biobank. These included samples from 98 individuals with signs and symptoms of acute Lyme disease, 25 healthy individuals residing in Lyme disease endemic areas, and 35 samples collected from patients admitted to the Massachusetts General Hospital or referred to the infectious disease clinic. Compared to PCR, TBDCapSeq had better sensitivity and could identify infections with a wider range of tick-borne agents. TBDCapSeq identified a higher rate of samples positive for Borrelia burgdorferi (8 vs. 1 by PCR) and Babesia microti (26 vs. 15 by PCR). TBDCapSeq also identified previously unknown infections with Borrelia miyamotoi, Ehrlichia, and Rickettsia species. Overall, TBDCapSeq identified a pathogen in 43 samples vs. 23 using PCR, with four co-infections detected versus zero by PCR. We conclude that capture sequencing enables superior detection of tick-borne agents relative to PCR.

Tension promotes kinetochore-microtubule release by Aurora B kinase.

To ensure accurate chromosome segregation, interactions between kinetochores and microtubules are regulated by a combination of mechanics and biochemistry. Tension provides a signal to discriminate attachment errors from bi-oriented kinetochores with sisters correctly attached to opposite spindle poles. Biochemically, Aurora B kinase phosphorylates kinetochores to destabilize interactions with microtubules. To link mechanics and biochemistry, current models regard tension as an input signal to locally regulate Aurora B activity. Here, we show that the outcome of kinetochore phosphorylation depends on tension. Using optogenetics to manipulate Aurora B at individual kinetochores, we find that kinase activity promotes microtubule release when tension is high. Conversely, when tension is low, Aurora B activity promotes depolymerization of kinetochore-microtubules while maintaining attachment. Thus, phosphorylation converts a catch-bond, in which tension stabilizes attachments, to a slip-bond, which releases microtubules under tension. We propose that tension is a signal inducing distinct error-correction pathways, with release or depolymerization being advantageous for typical errors characterized by high or low tension, respectively.

Development of a capture sequencing assay for enhanced detection and genotyping of tick-borne pathogens.

Inadequate sensitivity has been the primary limitation for implementing high-throughput sequencing for studies of tick-borne agents. Here we describe the development of TBDCapSeq, a sequencing assay that uses hybridization capture probes that cover the complete genomes of the eleven most common tick-borne agents found in the United States. The probes are used for solution-based capture and enrichment of pathogen nucleic acid followed by high-throughput sequencing. We evaluated the performance of TBDCapSeq to surveil samples that included human whole blood, mouse tissues, and field-collected ticks. For Borrelia burgdorferi and Babesia microti, the sensitivity of TBDCapSeq was comparable and occasionally exceeded the performance of agent-specific quantitative PCR and resulted in 25 to > 10,000-fold increase in pathogen reads when compared to standard unbiased sequencing. TBDCapSeq also enabled genome analyses directly within vertebrate and tick hosts. The implementation of TBDCapSeq could have major impact in studies of tick-borne pathogens by improving detection and facilitating genomic research that was previously unachievable with standard sequencing approaches.

A polyclonal allelic expression assay for detecting regulatory effects of transcript variants.

We present an assay to experimentally test the regulatory effects of genetic variants within transcripts using CRISPR/Cas9 followed by targeted sequencing. We applied the assay to 32 premature stop-gained variants across the genome and in two Mendelian disease genes, 33 putative causal variants of eQTLs, and 62 control variants in HEK293T cells, replicating a subset of variants in HeLa cells. We detected significant effects in the expected direction (in 60% of variants), demonstrating the ability of the assay to capture regulatory effects of eQTL variants and nonsense-mediated decay triggered by premature stop-gained variants. The results suggest a utility for validating transcript-level effects of genetic variants.