COMBREX: a project to accelerate the functional annotation of prokaryotic genomes.

COMBREX (http://combrex.bu.edu) is a project to increase the speed of the functional annotation of new bacterial and archaeal genomes. It consists of a database of functional predictions produced by computational biologists and a mechanism for experimental biochemists to bid for the validation of those predictions. Small grants are available to support successful bids.

Antibody titer levels and the effect on subsequent SARS-CoV-2 infection in a large US-based cohort.

Despite a growing amount of data around the kinetics and durability of the antibody response induced by vaccination and previous infection, there is little understanding of whether or not a given quantitative level of antibodies correlates to protection against SARS-CoV-2 infection or reinfection. In this study, we examine SARS-CoV-2 anti-spike receptor binding domain (RBD) antibody titers and subsequent SARS-CoV-2 reverse transcription polymerase chain reaction (RT-PCR) tests in a large cohort of US-based patients. We analyzed antibody test results in a cohort of 22,204 individuals, 6.8% (n = 1,509) of whom eventually tested positive for SARS-CoV-2 RNA, suggesting infection or reinfection. Kaplan-Meier curves were plotted to understand the effect of various levels of anti-spike RBD antibody titers (classified into discrete ranges) on subsequent RT-PCR positivity rates. Statistical analyses included fitting a Cox proportional hazards model to estimate the age-, sex- and exposure-adjusted hazard ratios for S antibody titer, using zip-code positivity rates by week as a proxy for COVID-19 exposure. It was found that the best models of the temporally associated infection risk were those based on log antibody titer level (HR = 0.836 (p < 0.05)). When titers were binned, the hazard ratio associated with antibody titer >250 Binding Antibody Units (BAU) was 0.27 (p < 0.05, 95% CI [0.18, 0.41]), while the hazard ratio associated with previous infection was 0.20 (p < 0.05, 95% CI [0.10, 0.39]). Fisher exact odds ratio (OR) for Ab titers <250 BAU showed OR = 2.84 (p < 0.05; 95% CI: [2.30, 3.53]) for predicting the outcome of a subsequent PCR test. Antibody titer levels correlate with protection against subsequent SARS-CoV-2 infection or reinfection when examining a cohort of real-world patients who had the spike RBD antibody assay performed.

Follow-Up SARS-CoV-2 PCR Testing Outcomes From a Large Reference Lab in the US.

By analyzing COVID-19 sequential COVID-19 test results of patients across the United States, we herein attempt to quantify some of the observations we've made around long-term infection (and false-positive rates), as well as provide observations on the uncertainty of sampling variability and other dynamics of COVID-19 infection in the United States. Retrospective cohort study of a registry of RT-PCR testing results for all patients tested at any of the reference labs operated by Labcorp® including both positive, negative, and inconclusive results, from March 1, 2020 to January 28, 2021, including patients from all 50 states and outlying US territories. The study included 22 million patients with RT-PCR qualitative test results for SARS-CoV-2, of which 3.9 million had more than one test at Labcorp. We observed a minuscule <0.1% basal positive rate for follow up tests >115 days, which could account for false positives, long-haulers, and/or reinfection but is indistinguishable in the data. In observing repeat-testing, for patients who have a second test after a first RT-PCR, 30% across the cohort tested negative on the second test. For patients who test positive first and subsequently negative within 96 h (40% of positive test results), 18% of tests will subsequently test positive within another 96-h span. For those who first test negative and then positive within 96 h (2.3% of negative tests), 56% will test negative after a third and subsequent 96-h period. The sudden changes in RT-PCR test results for SARS-CoV-2 from this large cohort study suggest that negative test results during active infection or exposure can change rapidly within just days or hours. We also demonstrate that there does not appear to be a basal false positive rate among patients who test positive >115 days after their first RT-PCR positive test while failing to observe any evidence of widespread reinfection.

Protocol dependence of sequencing-based gene expression measurements.

RNA Seq provides unparalleled levels of information about the transcriptome including precise expression levels over a wide dynamic range. It is essential to understand how technical variation impacts the quality and interpretability of results, how potential errors could be introduced by the protocol, how the source of RNA affects transcript detection, and how all of these variations can impact the conclusions drawn. Multiple human RNA samples were used to assess RNA fragmentation, RNA fractionation, cDNA synthesis, and single versus multiple tag counting. Though protocols employing polyA RNA selection generate the highest number of non-ribosomal reads and the most precise measurements for coding transcripts, such protocols were found to detect only a fraction of the non-ribosomal RNA in human cells. PolyA RNA excludes thousands of annotated and even more unannotated transcripts, resulting in an incomplete view of the transcriptome. Ribosomal-depleted RNA provides a more cost-effective method for generating complete transcriptome coverage. Expression measurements using single tag counting provided advantages for assessing gene expression and for detecting short RNAs relative to multi-read protocols. Detection of short RNAs was also hampered by RNA fragmentation. Thus, this work will help researchers choose from among a range of options when analyzing gene expression, each with its own advantages and disadvantages.

RNA sequencing and quantitation using the Helicos Genetic Analysis System.

The recent transition in gene expression analysis technology to ultra high-throughput cDNA sequencing provides a means for higher quantitation sensitivity across a wider dynamic range than previously possible. Sensitivity of detection is mostly a function of the sheer number of sequence reads generated. Typically, RNA is converted to cDNA using random hexamers and the cDNA is subsequently sequenced (RNA-Seq). With this approach, higher read numbers are generated for long transcripts as compared to short ones. This length bias necessitates the generation of very high read numbers to achieve sensitive quantitation of short, low-expressed genes. To eliminate this length bias, we have developed an ultra high-throughput sequencing approach where only a single read is generated for each transcript molecule (single-molecule sequencing Digital Gene Expression (smsDGE)). So, for example, equivalent quantitation accuracy of the yeast transcriptome can be achieved by smsDGE using only 25% of the reads that would be required using RNA-Seq. For sample preparation, RNA is first reverse-transcribed into single-stranded cDNA using oligo-dT as a primer. A poly-A tail is then added to the 3' ends of cDNA to facilitate the hybridization of the sample to the Helicos(®) single-molecule sequencing Flow-Cell to which a poly dT oligo serves as the substrate for subsequent sequencing by synthesis. No PCR, sample-size selection, or ligation steps are required, thus avoiding possible biases that may be introduced by such manipulations. Each tailed cDNA sample is injected into one of 50 flow-cell channels and sequenced on the Helicos(®) Genetic Analysis System. Thus, 50 samples are sequenced simultaneously generating 10-20 million sequence reads on average for each sample channel. The sequence reads can then be aligned to the reference of choice such as the transcriptome, for quantitation of known transcripts, or the genome for novel transcript discovery. This chapter provides a summary of the methods required for smsDGE.

Aberrant overexpression of satellite repeats in pancreatic and other epithelial cancers.

Satellite repeats in heterochromatin are transcribed into noncoding RNAs that have been linked to gene silencing and maintenance of chromosomal integrity. Using digital gene expression analysis, we showed that these transcripts are greatly overexpressed in mouse and human epithelial cancers. In 8 of 10 mouse pancreatic ductal adenocarcinomas (PDACs), pericentromeric satellites accounted for a mean 12% (range 1 to 50%) of all cellular transcripts, a mean 40-fold increase over that in normal tissue. In 15 of 15 human PDACs, alpha satellite transcripts were most abundant and HSATII transcripts were highly specific for cancer. Similar patterns were observed in cancers of the lung, kidney, ovary, colon, and prostate. Derepression of satellite transcripts correlated with overexpression of the long interspersed nuclear element 1 (LINE-1) retrotransposon and with aberrant expression of neuroendocrine-associated genes proximal to LINE-1 insertions. The overexpression of satellite transcripts in cancer may reflect global alterations in heterochromatin silencing and could potentially be useful as a biomarker for cancer detection.

Error tolerant indexing and alignment of short reads with covering template families.

The rapid adoption of high-throughput next generation sequence data in biological research is presenting a major challenge for sequence alignment tools­specifically, the efficient alignment of vast amounts of short reads to large references in the presence of differences arising from sequencing errors and biological sequence variations. To address this challenge, we developed a short read aligner for high-throughput sequencer data that is tolerant of errors or mutations of all types­namely, substitutions, deletions, and insertions. The aligner utilizes a multi-stage approach in which template-based indexing is used to identify candidate regions for alignment with dynamic programming. A template is a pair of gapped seeds, with one used with the read and one used with the reference. In this article, we focus on the development of template families that yield error-tolerant indexing up to a given error-budget. A general algorithm for finding those families is presented, and a recursive construction that creates families with higher error tolerance from ones with a lower error tolerance is developed.

Biological process linkage networks.

BACKGROUND: The traditional approach to studying complex biological networks is based on the identification of interactions between internal components of signaling or metabolic pathways. By comparison, little is known about interactions between higher order biological systems, such as biological pathways and processes. We propose a methodology for gleaning patterns of interactions between biological processes by analyzing protein-protein interactions, transcriptional co-expression and genetic interactions. At the heart of the methodology are the concept of Linked Processes and the resultant network of biological processes, the Process Linkage Network (PLN). RESULTS: We construct, catalogue, and analyze different types of PLNs derived from different data sources and different species. When applied to the Gene Ontology, many of the resulting links connect processes that are distant from each other in the hierarchy, even though the connection makes eminent sense biologically. Some others, however, carry an element of surprise and may reflect mechanisms that are unique to the organism under investigation. In this aspect our method complements the link structure between processes inherent in the Gene Ontology, which by its very nature is species-independent. As a practical application of the linkage of processes we demonstrate that it can be effectively used in protein function prediction, having the power to increase both the coverage and the accuracy of predictions, when carefully integrated into prediction methods. CONCLUSIONS: Our approach constitutes a promising new direction towards understanding the higher levels of organization of the cell as a system which should help current efforts to re-engineer ontologies and improve our ability to predict which proteins are involved in specific biological processes.

Quantification of the yeast transcriptome by single-molecule sequencing.

We present single-molecule sequencing digital gene expression (smsDGE), a high-throughput, amplification-free method for accurate quantification of the full range of cellular polyadenylated RNA transcripts using a Helicos Genetic Analysis system. smsDGE involves a reverse-transcription and polyA-tailing sample preparation procedure followed by sequencing that generates a single read per transcript. We applied smsDGE to the transcriptome of Saccharomyces cerevisiae strain DBY746, using 6 of the available 50 channels in a single sequencing run, yielding on average 12 million aligned reads per channel. Using spiked-in RNA, accurate quantitative measurements were obtained over four orders of magnitude. High correlation was demonstrated across independent flow-cell channels, instrument runs and sample preparations. Transcript counting in smsDGE is highly efficient due to the representation of each transcript molecule by a single read. This efficiency, coupled with the high throughput enabled by the single-molecule sequencing platform, provides an alternative method for expression profiling.

Comparative genome sequencing of Drosophila pseudoobscura: chromosomal, gene, and cis-element evolution.

We have sequenced the genome of a second Drosophila species, Drosophila pseudoobscura, and compared this to the genome sequence of Drosophila melanogaster, a primary model organism. Throughout evolution the vast majority of Drosophila genes have remained on the same chromosome arm, but within each arm gene order has been extensively reshuffled, leading to a minimum of 921 syntenic blocks shared between the species. A repetitive sequence is found in the D. pseudoobscura genome at many junctions between adjacent syntenic blocks. Analysis of this novel repetitive element family suggests that recombination between offset elements may have given rise to many paracentric inversions, thereby contributing to the shuffling of gene order in the D. pseudoobscura lineage. Based on sequence similarity and synteny, 10,516 putative orthologs have been identified as a core gene set conserved over 25-55 million years (Myr) since the pseudoobscura/melanogaster divergence. Genes expressed in the testes had higher amino acid sequence divergence than the genome-wide average, consistent with the rapid evolution of sex-specific proteins. Cis-regulatory sequences are more conserved than random and nearby sequences between the species--but the difference is slight, suggesting that the evolution of cis-regulatory elements is flexible. Overall, a pattern of repeat-mediated chromosomal rearrangement, and high coadaptation of both male genes and cis-regulatory sequences emerges as important themes of genome divergence between these species of Drosophila.

Whole-genome annotation by using evidence integration in functional-linkage networks.

The advent of high-throughput biology has catalyzed a remarkable improvement in our ability to identify new genes. A large fraction of newly discovered genes have an unknown functional role, particularly when they are specific to a particular lineage or organism. These genes, currently labeled "hypothetical," might support important biological cell functions and could potentially serve as targets for medical, diagnostic, or pharmacogenomic studies. An important challenge to the scientific community is to associate these newly predicted genes with a biological function that can be validated by experimental screens. In the absence of sequence or structural homology to known genes, we must rely on advanced biotechnological methods, such as DNA chips and protein-protein interaction screens as well as computational techniques to assign putative functions to these genes. In this article, we propose an effective methodology for combining biological evidence obtained in several high-throughput experimental screens and integrating this evidence in a way that provides consistent functional assignments to hypothetical genes. We use the visualization method of propagation diagrams to illustrate the flow of functional evidence that supports the functional assignments produced by the algorithm. Our results contain a number of predictions and furnish strong evidence that integration of functional information is indeed a promising direction for improving the accuracy and robustness of functional genomics.