Comparative Methods for Demystifying Spatial Transcriptomics.

Spatial Transcriptomics (ST), coined as the term for parallel RNA-Seq on cell populations ordered spatially on a histological tissue section, has recently become increasingly popular, especially in experiments where microfluidics-based single-cell sequencing fails, such as assays on neurons. ST platforms, like the 10x Visium technology investigated herein, therefore produce in a single experiment simultaneously thousands of RNA readouts, captured by an array of micrometer scale spots under the histological section. Therefore, a central challenge of analyzing ST experiments consists of analyzing the gene expression morphology of all spots to delineate clusters of similar cell mixtures, which are then compared to each other to identify up- or down-regulated marker genes. Moreover, another level of complexity in ST experiments, compared to traditional RNA-Seq, is imposed by staining the tissue section with protein markers of cells or cell components to identify spots providing relevant information afterward. The corresponding microscopy images need to be analyzed in addition to the RNA-Seq read mappings on the reference genome and transcriptome sequences. Focusing on the software suite provided by the Visium platform manufacturer, we break down the ST analysis pipeline into its four essential steps-the image analysis, the read alignment, the gene quantification, and the spot clustering-and compare results obtained when using reads from different subsets of spots and/or when employing alternative genome or transcriptome references. Our comparative analyses demonstrate the impact of spot selection and the choice of genome/transcriptome references on the analysis results when employing the manufacturer's pipeline.

Reliability of continuous vital sign monitoring in post-operative patients employing consumer-grade fitness trackers: A randomised pilot trial.

Introduction: Fitness trackers can provide continuous monitoring of vital signs and thus have the potential to become a complementary, mobile and effective tool for early detection of patient deterioration and post-operative complications. Methods: To evaluate potential implementations in acute care setting, we included 36 patients after moderate to major surgery in a recent randomised pilot trial to compare the performance of vital sign monitoring by three different fitness trackers (Apple Watch 7, Garmin Fenix 6pro and Withings ScanWatch) with established standard clinical monitors in post-anaesthesia care units and monitoring wards. Results: During a cumulative period of 56 days, a total of 53,197 heart rate (HR) measurements, as well as 12,219 measurements of the peripheral blood oxygen saturation (SpO2) and 28,954 respiratory rate (RR) measurements were collected by fitness trackers. Under real-world conditions, HR monitoring was accurate and reliable across all benchmarked devices (r = [0.95;0.98], p < 0.001; Bias = [-0.74 bpm;-0.01 bpm]; MAPE∼2%). However, the performance of SpO2 (r = [0.21;0.68]; p < 0.001; Bias = [-0.46%;-2.29%]; root-mean-square error = [2.82%;4.1%]) monitoring was substantially inferior. RR measurements could not be obtained for two of the devices, therefore exclusively the accuracy of the Garmin tracker could be evaluated (r = 0.28, p < 0.001; Bias = -1.46/min). Moreover, the time resolution of the vital sign measurements highly depends on the tracking device, ranging from 0.7 to 117.94 data points per hour. Conclusion: According to the results of the present study, tracker devices are generally reliable and accurate for HR monitoring, whereas SpO2 and RR measurements should be interpreted carefully, considering the clinical context of the respective patients.

Evaluating blood oxygen saturation measurements by popular fitness trackers in postoperative patients: A prospective clinical trial.

Blood oxygen saturation is an important clinical parameter, especially in postoperative hospitalized patients, monitored in clinical practice by arterial blood gas (ABG) and/or pulse oximetry that both are not suitable for a long-term continuous monitoring of patients during the entire hospital stay, or beyond. Technological advances developed recently for consumer-grade fitness trackers could-at least in theory-help to fill in this gap, but benchmarks on the applicability and accuracy of these technologies in hospitalized patients are currently lacking. We therefore conducted at the postanaesthesia care unit under controlled settings a prospective clinical trial with 201 patients, comparing in total >1,000 oxygen blood saturation measurements by fitness trackers of three brands with the ABG gold standard and with pulse oximetry. Our results suggest that, despite of an overall still tolerable measuring accuracy, comparatively high dropout rates severely limit the possibilities of employing fitness trackers, particularly during the immediate postoperative period of hospitalized patients.

Accuracy and Systematic Biases of Heart Rate Measurements by Consumer-Grade Fitness Trackers in Postoperative Patients: Prospective Clinical Trial.

BACKGROUND: Over the recent years, technological advances of wrist-worn fitness trackers heralded a new era in the continuous monitoring of vital signs. So far, these devices have primarily been used for sports. OBJECTIVE: However, for using these technologies in health care, further validations of the measurement accuracy in hospitalized patients are essential but lacking to date. METHODS: We conducted a prospective validation study with 201 patients after moderate to major surgery in a controlled setting to benchmark the accuracy of heart rate measurements in 4 consumer-grade fitness trackers (Apple Watch 7, Garmin Fenix 6 Pro, Withings ScanWatch, and Fitbit Sense) against the clinical gold standard (electrocardiography). RESULTS: All devices exhibited high correlation (r≥0.95; P<.001) and concordance (rc≥0.94) coefficients, with a relative error as low as mean absolute percentage error <5% based on 1630 valid measurements. We identified confounders significantly biasing the measurement accuracy, although not at clinically relevant levels (mean absolute error<5 beats per minute). CONCLUSIONS: Consumer-grade fitness trackers appear promising in hospitalized patients for monitoring heart rate. TRIAL REGISTRATION: ClinicalTrials.gov NCT05418881; https://www.clinicaltrials.gov/ct2/show/NCT05418881.

The Use of Non-Invasive Continuous Blood Pressure Measuring (ClearSight®) during Central Neuraxial Anaesthesia for Caesarean Section-A Retrospective Validation Study.

The close monitoring of blood pressure during a caesarean section performed under central neuraxial anaesthesia should be the standard of safe anaesthesia. As classical oscillometric and invasive blood pressure measuring have intrinsic disadvantages, we investigated a novel, non-invasive technique for continuous blood pressure measuring. Methods: In this monocentric, retrospective data analysis, the reliability of continuous non-invasive blood pressure measuring using ClearSight® (Edwards Lifesciences Corporation) is validated in 31 women undergoing central neuraxial anaesthesia for caesarean section. In addition, patients and professionals evaluated ClearSight® through questioning. Results: 139 measurements from 11 patients were included in the final analysis. Employing Bland-Altman analyses, we identified a bias of -10.8 mmHg for systolic, of -0.45 mmHg for diastolic and of +0.68 mmHg for mean arterial blood pressure measurements. Pooling all paired measurements resulted in a Pearson correlation coefficient of 0.7 for systolic, of 0.67 for diastolic and of 0.75 for mean arterial blood pressure. Compensating the interindividual differences in linear regressions of the paired measurements provided improved correlation coefficients of 0.73 for systolic, of 0.9 for diastolic and of 0.89 for mean arterial blood pressure measurements. Discussion: Diastolic and mean arterial blood pressure are within an acceptable range of deviation from the reference method, according to the Association for the Advancement of Medical Instrumentation (AAMI) in the patient collective under study. Both patients and professionals prefer ClearSight® to oscillometric blood pressure measurement in regard of comfort and handling.

Analysis of ovarian transcriptomes reveals thousands of novel genes in the insect vector Rhodnius prolixus.

Rhodnius prolixus is a Triatominae insect species and a primary vector of Chagas disease. The genome of R. prolixus has been recently sequenced and partially assembled, but few transcriptome analyses have been performed to date. In this study, we describe the stage-specific transcriptomes obtained from previtellogenic stages of oogenesis and from mature eggs. By analyzing ~ 228 million paired-end RNA-Seq reads, we significantly improved the current genome annotations for 9206 genes. We provide extended 5' and 3' UTRs, complete Open Reading Frames, and alternative transcript variants. Strikingly, using a combination of genome-guided and de novo transcriptome assembly we found more than two thousand novel genes, thus increasing the number of genes in R. prolixus from 15,738 to 17,864. We used the improved transcriptome to investigate stage-specific gene expression profiles during R. prolixus oogenesis. Our data reveal that 11,127 genes are expressed in the early previtellogenic stage of oogenesis and their transcripts are deposited in the developing egg including key factors regulating germline development, genome integrity, and the maternal-zygotic transition. In addition, GO term analyses show that transcripts encoding components of the steroid hormone receptor pathway, cytoskeleton, and intracellular signaling are abundant in the mature eggs, where they likely control early embryonic development upon fertilization. Our results significantly improve the R. prolixus genome and transcriptome and provide novel insight into oogenesis and early embryogenesis in this medically relevant insect.

Echocardiographic Measurements in a Preclinical Model of Chronic Chagasic Cardiomyopathy in Dogs: Validation and Reproducibility.

Background: The failure to translate preclinical results to the clinical setting is the rule, not the exception. One reason that is frequently overlooked is whether the animal model reproduces distinctive features of human disease. Another is the reproducibility of the method used to measure treatment effects in preclinical studies. Left ventricular (LV) function improvement is the most common endpoint in preclinical cardiovascular disease studies, while echocardiography is the most frequently used method to evaluate LV function. In this work, we conducted a robust echocardiographic evaluation of LV size and function in dogs chronically infected by Trypanosoma cruzi. Methods and Results: Echocardiography was performed blindly by two distinct observers in mongrel dogs before and between 6 and 9 months post infection. Parameters analyzed included end-systolic volume (ESV), end-diastolic volume (EDV), ejection fraction (EF), and fractional shortening (FS). We observed a significant LVEF and FS reduction in infected animals compared to controls, with no significant variation in volumes. However, the effect of chronic infection in systolic function was quite variable, with EF ranging from 17 to 66%. Using the cut-off value of EF ≤ 40%, established for dilated cardiomyopathy (DCM) in dogs, only 28% of the infected dogs were affected by the chronic infection. Conclusions: The canine model of CCC mimics human disease, reproducing the percentage of individuals that develop heart failure during the chronic infection. It is thus mandatory to establish inclusion criteria in the experimental design of canine preclinical studies to account for the variable effect that chronic infection has on systolic function.

Characterization of HPV integration, viral gene expression and E6E7 alternative transcripts by RNA-Seq: A descriptive study in invasive cervical cancer.

Scarce data are available on the expression of papillomavirus genome and the frequency of alternatively spliced E6E7 mRNAs in invasive cervical cancer. We carried out a comprehensive characterization of HPV expression by RNA-Seq analysis in 22 invasive cervical cancer with HPV16 or HPV18, characterizing the presence of integrated/episomal viral DNA, the integration sites in human genome and the proportion of alternative splicing products of E6 and E7 genes. The expression patterns suggested the presence of episomal and/or integrated viral DNA, with integration detected in most tumors, frequently occurring within human genes in HPV18+ and in intergenic regions in HPV16+ tumors. Alternative splicing of E6E7 transcripts showed E6*I as the most frequent isoform for both viral types, followed by E6*II and E6/E7 (unspliced) transcripts in HPV16+, and by E6/E7 in HPV18+ tumors. Previously described E6*VI and E6*V transcript isoforms for HPV16, and E6*X for HPV18, were rare or not detected.

Transcriptomic and functional analyses of the piRNA pathway in the Chagas disease vector Rhodnius prolixus.

The piRNA pathway is a surveillance system that guarantees oogenesis and adult fertility in a range of animal species. The pathway is centered on PIWI clade Argonaute proteins and the associated small non-coding RNAs termed piRNAs. In this study, we set to investigate the evolutionary conservation of the piRNA pathway in the hemimetabolous insect Rhodnius prolixus. Our transcriptome profiling reveals that core components of the pathway are expressed during previtellogenic stages of oogenesis. Rhodnius' genome harbors four putative piwi orthologs. We show that Rp-piwi2, Rp-piwi3 and Rp-ago3, but not Rp-piwi1 transcripts are produced in the germline tissues and maternally deposited in the mature eggs. Consistent with a role in Rhodnius oogenesis, parental RNAi against the Rp-piwi2, Rp-piwi3 and Rp-ago3 results in severe egg laying and female adult fertility defects. Furthermore, we show that the reduction of the Rp-piwi2 levels by parental RNAi disrupts oogenesis by causing a dramatic loss of trophocytes, egg chamber degeneration and oogenesis arrest. Intriguingly, the putative Rp-Piwi2 protein features a polyglutamine tract at its N-terminal region, which is conserved in PIWI proteins encoded in the genome of other Triatomine species. Together with R. prolixus, these hematophagous insects are primary vectors of the Chagas disease. Thus, our data shed more light on the evolution of the piRNA pathway and provide a framework for the development of new control strategies for Chagas disease insect vectors.

An automated method for detecting alternatively spliced protein domains.

Motivation: Alternative splicing (AS) has been demonstrated to play a role in shaping eukaryotic gene diversity at the transcriptional level. However, the impact of AS on the proteome is still controversial. Studies that seek to explore the effect of AS at the proteomic level are hampered by technical difficulties in the cumbersome process of casting forth and back between genome, transcriptome and proteome space coordinates, and the naïve prediction of protein domains in the presence of AS suffers many redundant sequence scans that emerge from constitutively spliced regions that are shared between alternative products of a gene. Results: We developed the AstaFunk pipeline that computes for every generic transcriptome all domains that are altered by AS events in a systematic and efficient manner. In a nutshell, our method employs Viterbi dynamic programming, which guarantees to find all score-optimal hits of the domains under consideration, while complementary optimizations at different levels avoid redundant and other irrelevant computations. We evaluate AstaFunk qualitatively and quantitatively using RNAseq in well-studied genes with AS, and on large-scale employing entire transcriptomes. Our study confirms complementary reports that the effect of most AS events on the proteome seems to be rather limited, but our results also pinpoint several cases where AS could have a major impact on the function of a protein domain. Availability and implementation: The JAVA implementation of AstaFunk is available as an open source project on http://astafunk.sammeth.net. Supplementary information: Supplementary data are available at Bioinformatics online.

Landscape of the spliced leader trans-splicing mechanism in Schistosoma mansoni.

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

The effects of death and post-mortem cold ischemia on human tissue transcriptomes.

Post-mortem tissues samples are a key resource for investigating patterns of gene expression. However, the processes triggered by death and the post-mortem interval (PMI) can significantly alter physiologically normal RNA levels. We investigate the impact of PMI on gene expression using data from multiple tissues of post-mortem donors obtained from the GTEx project. We find that many genes change expression over relatively short PMIs in a tissue-specific manner, but this potentially confounding effect in a biological analysis can be minimized by taking into account appropriate covariates. By comparing ante- and post-mortem blood samples, we identify the cascade of transcriptional events triggered by death of the organism. These events do not appear to simply reflect stochastic variation resulting from mRNA degradation, but active and ongoing regulation of transcription. Finally, we develop a model to predict the time since death from the analysis of the transcriptome of a few readily accessible tissues.

Comparative Genomics in Drosophila.

Since the pioneering studies of Thomas Hunt Morgan and coworkers at the dawn of the twentieth century, Drosophila melanogaster and its sister species have tremendously contributed to unveil the rules underlying animal genetics, development, behavior, evolution, and human disease. Recent advances in DNA sequencing technologies launched Drosophila into the post-genomic era and paved the way for unprecedented comparative genomics investigations. The complete sequencing and systematic comparison of the genomes from 12 Drosophila species represents a milestone achievement in modern biology, which allowed a plethora of different studies ranging from the annotation of known and novel genomic features to the evolution of chromosomes and, ultimately, of entire genomes. Despite the efforts of countless laboratories worldwide, the vast amount of data that were produced over the past 15 years is far from being fully explored.In this chapter, we will review some of the bioinformatic approaches that were developed to interrogate the genomes of the 12 Drosophila species. Setting off from alignments of the entire genomic sequences, the degree of conservation can be separately evaluated for every region of the genome, providing already first hints about elements that are under purifying selection and therefore likely functional. Furthermore, the careful analysis of repeated sequences sheds light on the evolutionary dynamics of transposons, an enigmatic and fascinating class of mobile elements housed in the genomes of animals and plants. Comparative genomics also aids in the computational identification of the transcriptionally active part of the genome, first and foremost of protein-coding loci, but also of transcribed nevertheless apparently noncoding regions, which were once considered "junk" DNA. Eventually, the synergy between functional and comparative genomics also facilitates in silico and in vivo studies on cis-acting regulatory elements, like transcription factor binding sites, that due to the high degree of sequence variability usually impose increased challenges for bioinformatics approaches.

Comparative Genomics in Homo sapiens.

Genomes can be compared at different levels of divergence, either between species or within species. Within species genomes can be compared between different subpopulations, such as human subpopulations from different continents. Investigating the genomic differences between different human subpopulations is important when studying complex diseases that are affected by many genetic variants, as the variants involved can differ between populations. The 1000 Genomes Project collected genome-scale variation data for 2504 human individuals from 26 different populations, enabling a systematic comparison of variation between human subpopulations. In this chapter, we present step-by-step a basic protocol for the identification of population-specific variants employing the 1000 Genomes data. These variants are subsequently further investigated for those that affect the proteome or RNA splice sites, to investigate potentially biologically relevant differences between the populations.

Correction: Gene Expansion Shapes Genome Architecture in the Human Pathogen Lichtheimia corymbifera: An Evolutionary Genomics Analysis in the Ancient Terrestrial Mucorales (Mucoromycotina).

[This corrects the article DOI: 10.1371/journal.pgen.1004496.].

Sequence variation between 462 human individuals fine-tunes functional sites of RNA processing.

Recent advances in the cost-efficiency of sequencing technologies enabled the combined DNA- and RNA-sequencing of human individuals at the population-scale, making genome-wide investigations of the inter-individual genetic impact on gene expression viable. Employing mRNA-sequencing data from the Geuvadis Project and genome sequencing data from the 1000 Genomes Project we show that the computational analysis of DNA sequences around splice sites and poly-A signals is able to explain several observations in the phenotype data. In contrast to widespread assessments of statistically significant associations between DNA polymorphisms and quantitative traits, we developed a computational tool to pinpoint the molecular mechanisms by which genetic markers drive variation in RNA-processing, cataloguing and classifying alleles that change the affinity of core RNA elements to their recognizing factors. The in silico models we employ further suggest RNA editing can moonlight as a splicing-modulator, albeit less frequently than genomic sequence diversity. Beyond existing annotations, we demonstrate that the ultra-high resolution of RNA-Seq combined from 462 individuals also provides evidence for thousands of bona fide novel elements of RNA processing-alternative splice sites, introns, and cleavage sites-which are often rare and lowly expressed but in other characteristics similar to their annotated counterparts.

Human genomics. The human transcriptome across tissues and individuals.

Transcriptional regulation and posttranscriptional processing underlie many cellular and organismal phenotypes. We used RNA sequence data generated by Genotype-Tissue Expression (GTEx) project to investigate the patterns of transcriptome variation across individuals and tissues. Tissues exhibit characteristic transcriptional signatures that show stability in postmortem samples. These signatures are dominated by a relatively small number of genes­which is most clearly seen in blood­though few are exclusive to a particular tissue and vary more across tissues than individuals. Genes exhibiting high interindividual expression variation include disease candidates associated with sex, ethnicity, and age. Primary transcription is the major driver of cellular specificity, with splicing playing mostly a complementary role; except for the brain, which exhibits a more divergent splicing program. Variation in splicing, despite its stochasticity, may play in contrast a comparatively greater role in defining individual phenotypes.

Human genomics. Effect of predicted protein-truncating genetic variants on the human transcriptome.

Accurate prediction of the functional effect of genetic variation is critical for clinical genome interpretation. We systematically characterized the transcriptome effects of protein-truncating variants, a class of variants expected to have profound effects on gene function, using data from the Genotype-Tissue Expression (GTEx) and Geuvadis projects. We quantitated tissue-specific and positional effects on nonsense-mediated transcript decay and present an improved predictive model for this decay. We directly measured the effect of variants both proximal and distal to splice junctions. Furthermore, we found that robustness to heterozygous gene inactivation is not due to dosage compensation. Our results illustrate the value of transcriptome data in the functional interpretation of genetic variants.

Analysis of alternative splicing events in custom gene datasets by AStalavista.

Alternative splicing (AS) is a eukaryotic principle to derive more than one RNA product from transcribed genes by removing distinct subsets of introns from a premature polymer. We know today that this process is highly regulated and makes up a large part of the differences between species, cell types, and states. The key to compare AS across different genes or organisms is to tokenize the AS phenomenon into atomary units, so-called AS events. These events then usually are grouped by common patterns to investigate the underlying molecular mechanisms that drive their regulation. However, attempts to decompose loci with AS observations into events are often hampered by applying a limited set of a priori defined event patterns which are not capable to describe all AS configurations and therefore cannot decompose the phenomenon exhaustively. In this chapter, we describe working scenarios of AStalavista, a computational method that reports all AS events reflected by transcript annotations. We show how to practically employ AStalavista to study AS variation in complex transcriptomes, as characterized by the human GENCODE annotation. Our examples demonstrate how the inherent and universal AStalavista paradigm allows for an automatic delineation of AS events in custom gene datasets. Additionally, we sketch an example of an AStalavista use case including next-generation sequencing data (RNA-Seq) to enrich the landscape of discovered AS events.

Nova1 is a master regulator of alternative splicing in pancreatic beta cells.

Alternative splicing (AS) is a fundamental mechanism for the regulation of gene expression. It affects more than 90% of human genes but its role in the regulation of pancreatic beta cells, the producers of insulin, remains unknown. Our recently published data indicated that the 'neuron-specific' Nova1 splicing factor is expressed in pancreatic beta cells. We have presently coupled specific knockdown (KD) of Nova1 with RNA-sequencing to determine all splice variants and downstream pathways regulated by this protein in beta cells. Nova1 KD altered the splicing of nearly 5000 transcripts. Pathway analysis indicated that these genes are involved in exocytosis, apoptosis, insulin receptor signaling, splicing and transcription. In line with these findings, Nova1 silencing inhibited insulin secretion and induced apoptosis basally and after cytokine treatment in rodent and human beta cells. These observations identify a novel layer of regulation of beta cell function, namely AS controlled by key splicing regulators such as Nova1.