Model-based analysis of sample index hopping reveals its widespread artifacts in multiplexed single-cell RNA-sequencing

Single-cell immune repertoire analysis

Single-cell T cell and B cell antigen receptor-sequencing data analysis can potentially perform in-depth assessments of adaptive immune cells that inform on understanding immune cell development to tracking clonal expansion in disease and therapy. However, it has been extremely challenging to analyze and interpret T cells and B cells and their adaptive immune receptor repertoires at the single-cell level due to not only the complexity of the data but also the underlying biology. In this Review, we delve into the computational breakthroughs that have transformed the analysis of single-cell T cell and B cell antigen receptor-sequencing data.

Characterization of the brain virome in human immunodeficiency virus infection and substance use disorder

Viruses can infect the brain in individuals with and without HIV-infection: however, the brain virome is poorly characterized. Metabolic alterations have been identified which predispose people to substance use disorder (SUD), but whether these could be triggered by viral infection of the brain is unknown. We used a target-enrichment, deep sequencing platform and bioinformatic pipeline named "ViroFind", for the unbiased characterization of DNA and RNA viruses in brain samples obtained from the National Neuro-AIDS Tissue Consortium. We analyzed fresh frozen post-mortem prefrontal cortex from 72 individuals without known viral infection of the brain, including 16 HIV+/SUD+, 20 HIV+/SUD-, 16 HIV-/SUD+, and 20 HIV-/SUD-. The average age was 52.3 y and 62.5% were males. We identified sequences from 26 viruses belonging to 11 viral taxa. These included viruses with and without known pathogenic potential or tropism to the nervous system, with sequence coverage ranging from 0.03 to 99.73% of the viral genomes. In SUD+ people, HIV-infection was associated with a higher total number of viruses, and HIV+/SUD+ compared to HIV-/SUD+ individuals had an increased frequency of Adenovirus (68.8 vs 0%; p<0.001) and Epstein-Barr virus (EBV) (43.8 vs 6.3%; p=0.037) as well as an increase in Torque Teno virus (TTV) burden. Conversely, in HIV+ people, SUD was associated with an increase in frequency of Hepatitis C virus, (25 in HIV+/SUD+ vs 0% in HIV+/SUD-; p=0.031). Finally, HIV+/SUD- compared to HIV-/SUD- individuals had an increased frequency of EBV (50 vs 0%; p<0.001) and an increase in TTV viral burden, but a decreased Adenovirus viral burden. These data demonstrate an unexpectedly high variety in the human brain virome, identifying targets for future research into the impact of these taxa on the central nervous system. ViroFind could become a valuable tool for monitoring viral dynamics in various compartments, monitoring outbreaks, and informing vaccine development.

Transcriptomic diversity of cell types across the adult human brain

The human brain directs complex behaviors, ranging from fine motor skills to abstract intelligence, but the diversity of cell types that support these skills has not been fully described. In this work, we used single-nucleus RNA sequencing to systematically survey cells across the entire adult human brain. We sampled more than three million nuclei from approximately 100 dissections across the forebrain, midbrain, and hindbrain in three postmortem donors. Our analysis identified 461 clusters and 3313 subclusters organized largely according to developmental origins and revealing high diversity in midbrain and hindbrain neurons. Astrocytes and oligodendrocyte-lineage cells also exhibited regional diversity at multiple scales. The transcriptomic census of the entire human brain presented in this work provides a resource for understanding the molecular diversity of the human brain in health and disease.

Single-cell multiomic understanding of HIV-1 reservoir at epigenetic, transcriptional, and protein levels

PURPOSE OF REVIEW

The success of HIV-1 eradication strategies relies on in-depth understanding of HIV-1-infected cells. However, HIV-1-infected cells are extremely heterogeneous and rare. Single-cell multiomic approaches resolve the heterogeneity and rarity of HIV-1-infected cells.

RECENT FINDINGS

Advancement in single-cell multiomic approaches enabled HIV-1 reservoir profiling across the epigenetic (ATAC-seq), transcriptional (RNA-seq), and protein levels (CITE-seq). Using HIV-1 RNA as a surrogate, ECCITE-seq identified enrichment of HIV-1-infected cells in clonally expanded cytotoxic CD4+ T cells. Using HIV-1 DNA PCR-activated microfluidic sorting, FIND-seq captured the bulk transcriptome of HIV-1 DNA+ cells. Using targeted HIV-1 DNA amplification, PheP-seq identified surface protein expression of intact versus defective HIV-1-infected cells. Using ATAC-seq to identify HIV-1 DNA, ASAP-seq captured transcription factor activity and surface protein expression of HIV-1 DNA+ cells. Combining HIV-1 mapping by ATAC-seq and HIV-1 RNA mapping by RNA-seq, DOGMA-seq captured the epigenetic, transcriptional, and surface protein expression of latent and transcriptionally active HIV-1-infected cells. To identify reproducible biological insights and authentic HIV-1-infected cells and avoid false-positive discovery of artifacts, we reviewed current practices of single-cell multiomic experimental design and bioinformatic analysis.

SUMMARY

Single-cell multiomic approaches may identify innovative mechanisms of HIV-1 persistence, nominate therapeutic strategies, and accelerate discoveries.

Decoding the impact of nuclear organization on antigenic variation in parasites

Antigenic variation as a strategy to evade the host adaptive immune response has evolved in divergent pathogens. Antigenic variation involves restricted, and often mutually exclusive, expression of dominant antigens and a periodic switch in antigen expression during infection. In eukaryotes, nuclear compartmentalization, including three-dimensional folding of the genome and physical separation of proteins in compartments or condensates, regulates mutually exclusive gene expression and chromosomal translocations. In this Review, we discuss the impact of nuclear organization on antigenic variation in the protozoan pathogens Trypanosoma brucei and Plasmodium falciparum. In particular, we highlight the relevance of nuclear organization in both mutually exclusive antigen expression and genome stability, which underlie antigenic variation.

A Refined, Controlled 16S rRNA Gene Sequencing Approach Reveals Limited Detection of Cerebrospinal Fluid Microbiota in Children with Bacterial Meningitis

Advances in both laboratory and computational components of high-throughput 16S amplicon sequencing (16S HTS) have markedly increased its sensitivity and specificity. Additionally, these refinements have better delineated the limits of sensitivity, and contributions of contamination to these limits, for 16S HTS that are particularly relevant for samples with low bacterial loads, such as human cerebrospinal fluid (CSF). The objectives of this work were to (i) optimize the performance of 16S HTS in CSF samples with low bacterial loads by defining and addressing potential sources of error, and (ii) perform refined 16S HTS on CSF samples from children diagnosed with bacterial meningitis and compare results with those from microbiological cultures. Several bench and computational approaches were taken to address potential sources of error for low bacterial load samples. We compared DNA yields and sequencing results after applying three different DNA extraction approaches to an artificially constructed mock-bacterial community. We also compared two postsequencing computational contaminant removal strategies, decontam R and full contaminant sequence removal. All three extraction techniques followed by decontam R yielded similar results for the mock community. We then applied these methods to 22 CSF samples from children diagnosed with meningitis, which has low bacterial loads relative to other clinical infection samples. The refined 16S HTS pipelines identified the cultured bacterial genus as the dominant organism for only 3 of these samples. We found that all three DNA extraction techniques followed by decontam R generated similar DNA yields for mock communities at the low bacterial loads representative of CSF samples. However, the limits of detection imposed by reagent contaminants and methodologic bias precluded the accurate detection of bacteria in CSF from children with culture-confirmed meningitis using these approaches, despite rigorous controls and sophisticated computational approaches. Although we did not find current DNA-based diagnostics to be useful for pediatric meningitis samples, the utility of these methods for CSF shunt infection remains undefined. Future advances in sample processing methods to minimize or eliminate contamination will be required to improve the sensitivity and specificity of these methods for pediatric meningitis. IMPORTANCE Advances in both laboratory and computational components of high-throughput 16S amplicon sequencing (16S HTS) have markedly increased its sensitivity and specificity. These refinements have better delineated the limits of sensitivity, and contributions of contamination to these limits, for 16S HTS that are particularly relevant for samples with low bacterial loads such as human cerebrospinal fluid (CSF). The objectives of this work were to (i) optimize the performance of 16S HTS in CSF samples by defining and addressing potential sources of error, and (ii) perform refined 16S HTS on CSF samples from children diagnosed with bacterial meningitis and compare results with those from microbiological cultures. We found that the limits of detection imposed by reagent contaminants and methodologic bias precluded the accurate detection of bacteria in CSF from children with culture-confirmed meningitis using these approaches, despite rigorous controls and sophisticated computational approaches.

DNA Capture and Enrichment: A Culture-Independent Approach for Characterizing the Genomic Diversity of Pathogenic Leptospira Species

Because they are difficult to culture, obtaining genomic information from Leptospira spp. is challenging, hindering the overall understanding of leptospirosis. We designed and validated a culture-independent DNA capture and enrichment system for obtaining Leptospira genomic information from complex human and animal samples. It can be utilized with a variety of complex sample types and diverse species as it was designed using the pan-genome of all known pathogenic Leptospira spp. This system significantly increases the proportion of Leptospira DNA contained within DNA extracts obtained from complex samples, oftentimes reaching >95% even when some estimated starting proportions were <1%. Sequencing enriched extracts results in genomic coverage similar to sequenced isolates, thereby enabling enriched complex extracts to be analyzed together with whole genome sequences from isolates, which facilitates robust species identification and high-resolution genotyping. The system is flexible and can be readily updated when new genomic information becomes available. Implementation of this DNA capture and enrichment system will improve efforts to obtain genomic data from unculturable Leptospira-positive human and animal samples. This, in turn, will lead to a better understanding of the overall genomic diversity and gene content of Leptospira spp. that cause leptospirosis, aiding epidemiology and the development of improved diagnostics and vaccines.

Metagenomic analysis reveals novel dietary-related viruses in the gut virome of marmosets hybrids (Callithrix jacchus x Callithrix penicillata), Brazil

Viral metagenomics has contributed enormously to the characterization of a wide range of viruses infecting animals of all phyla in the last decades. Among Neotropical primates, especially those introduced, knowledge about viral diversity remains poorly studied. Therefore, using metagenomics based on virus enrichment, we explored the viral microbiota present in the feces of introduced common marmosets (Callithrix sp.) in three locations from the Silva Jardim region in the State of Rio de Janeiro, Brazil. Fecal samples were collected from nine marmosets, pooled into three sample pools, and sequenced on Illumina MiSeq platform. Sequence reads were analyzed using a viral metagenomic analysis pipeline and two novel insect viruses belonging to the Parvoviridae and Baculoviridae families were identified. The complete genome of a densovirus (Parvoviridae family) of 5,309 nucleotides (nt) was obtained. The NS1 and VP1 proteins share lower than 32% sequence identity with the corresponding proteins of known members of the subfamily Densovirinae. Phylogenetic analysis suggests that this virus represents a new genus, provisionally named Afoambidensovirus due to its discovery in the Brazilian Atlantic Forest. The novel species received the name Afoambidensovirus incertum 1. The complete circular genome of a baculovirus of 107,191 nt was also obtained, showing 60.8% sequence identity with the most closely related member of the Baculoviridae family. Phylogenetic analysis suggests that this virus represents a new species in the Betabaculovirus genus, provisionally named Betabaculovirus incertum 1. In addition, sequences from several families of arthropods in the three pools evaluated were characterized (contigs ranging from 244 to 6,750 nt), corroborating the presence of possible insect hosts with which these new viruses may be associated. Our study expands the knowledge about two viral families known to infect insects, an important component of the marmosets' diet. This identification in hosts' feces samples demonstrates one of the many uses of this type of data and could serve as a basis for future research characterizing viruses in wildlife using noninvasive samples.

Forensic Analysis of Novel SARS2r-CoV Identified in Game Animal Datasets in China Shows Evolutionary Relationship to Pangolin GX CoV Clade and Apparent Genetic Experimentation

Pangolins are the only animals other than bats proposed to have been infected with SARS-CoV-2 related coronaviruses (SARS2r-CoVs) prior to the COVID-19 pandemic. Here, we examine the novel SARS2r-CoV we previously identified in game animal metatranscriptomic datasets sequenced by the Nanjing Agricultural University in 2022, and find that sections of the partial genome phylogenetically group with Guangxi pangolin CoVs (GX PCoVs), while the full RdRp sequence groups with bat-SL-CoVZC45. While the novel SARS2r-CoV is found in 6 pangolin datasets, it is also found in 10 additional NGS datasets from 5 separate mammalian species and is likely related to contamination by a laboratory researched virus. Absence of bat mitochondrial sequences from the datasets, the fragmentary nature of the virus sequence and the presence of a partial sequence of a cloning vector attached to a SARS2r-CoV read suggests that it has been cloned. We find that NGS datasets containing the novel SARS2r-CoV are contaminated with significant Homo sapiens genetic material, and numerous viruses not associated with the host animals sampled. We further identify the dominant human haplogroup of the contaminating H. sapiens genetic material to be F1c1a1, which is of East Asian provenance. The association of this novel SARS2r-CoV with both bat CoV and the GX PCoV clades is an important step towards identifying the origin of the GX PCoVs.

Sequencing introduced false positive rare taxa lead to biased microbial community diversity, assembly, and interaction interpretation in amplicon studies

BACKGROUND

Increasing studies have demonstrated potential disproportionate functional and ecological contributions of rare taxa in a microbial community. However, the study of the microbial rare biosphere is hampered by their inherent scarcity and the deficiency of currently available techniques. Sample-wise cross contaminations might be introduced by sample index misassignment in the most widely used metabarcoding amplicon sequencing approach. Although downstream bioinformatic quality control and clustering or denoising algorithms could remove sequencing errors and non-biological artifact reads, no algorithm could eliminate high quality reads from sample-wise cross contaminations introduced by index misassignment, making it difficult to distinguish between bona fide rare taxa and potential false positives in metabarcoding studies.

RESULTS

We thoroughly evaluated the rate of index misassignment of the widely used NovaSeq 6000 and DNBSEQ-G400 sequencing platforms using both commercial and customized mock communities, and observed significant lower (0.08% vs. 5.68%) fraction of potential false positive reads for DNBSEQ-G400 as compared to NovaSeq 6000. Significant batch effects could be caused by stochastically introduced false positive or false negative rare taxa. These false detections could also lead to inflated alpha diversity of relatively simple microbial communities and underestimated that of complex ones. Further test using a set of cow rumen samples reported differential rare taxa by different sequencing platforms. Correlation analysis of the rare taxa detected by each sequencing platform demonstrated that the rare taxa identified by DNBSEQ-G400 platform had a much higher possibility to be correlated with the physiochemical properties of rumen fluid as compared to NovaSeq 6000 platform. Community assembly mechanism and microbial network correlation analysis indicated that false positive or negative rare taxa detection could lead to biased community assembly mechanism and identification of fake keystone species of the community.

CONCLUSIONS

We highly suggest proper positive/negative/blank controls, technical replicate settings, and proper sequencing platform selection in future amplicon studies, especially when the microbial rare biosphere would be focused.

GigaAssay - An adaptable high-throughput saturation mutagenesis assay platform

High-throughput assay systems have had a large impact on understanding the mechanisms of basic cell functions. However, high-throughput assays that directly assess molecular functions are limited. Herein, we describe the "GigaAssay", a modular high-throughput one-pot assay system for measuring molecular functions of thousands of genetic variants at once. In this system, each cell was infected with one virus from a library encoding thousands of Tat mutant proteins, with each viral particle encoding a random unique molecular identifier (UMI). We demonstrate proof of concept by measuring transcription of a GFP reporter in an engineered reporter cell line driven by binding of the HIV Tat transcription factor to the HIV long terminal repeat. Infected cells were flow-sorted into 3 bins based on their GFP fluorescence readout. The transcriptional activity of each Tat mutant was calculated from the ratio of signals from each bin. The use of UMIs in the GigaAssay produced a high average accuracy (95%) and positive predictive value (98%) determined by comparison to literature benchmark data, known C-terminal truncations, and blinded independent mutant tests. Including the substitution tolerance with structure/function analysis shows restricted substitution types spatially concentrated in the Cys-rich region. Tat has abundant intragenic epistasis (10%) when single and double mutants are compared.

Clinical and Genetic Re-Evaluation of Inherited Retinal Degeneration Pedigrees following Initial Negative Findings on Panel-Based Next Generation Sequencing

Although rare, inherited retinal degenerations (IRDs) are the most common reason for blind registration in the working age population. They are highly genetically heterogeneous (>300 known genetic loci), and confirmation of a molecular diagnosis is a prerequisite for many therapeutic clinical trials and approved treatments. First-tier genetic testing of IRDs with panel-based next-generation sequencing (pNGS) has a diagnostic yield of ≈70-80%, leaving the remaining more challenging cases to be resolved by second-tier testing methods. This study describes the phenotypic reassessment of patients with a negative result from first-tier pNGS and the rationale, outcomes, and cost of second-tier genetic testing approaches. Removing non-IRD cases from consideration and utilizing case-appropriate second-tier genetic testing techniques, we genetically resolved 56% of previously unresolved pedigrees, bringing the overall resolve rate to 92% (388/423). At present, pNGS remains the most cost-effective first-tier approach for the molecular assessment of diverse IRD populations Second-tier genetic testing should be guided by clinical (i.e., reassessment, multimodal imaging, electrophysiology), and genetic (i.e., single alleles in autosomal recessive disease) indications to achieve a genetic diagnosis in the most cost-effective manner.

Early sample tagging and pooling enables simultaneous SARS-CoV-2 detection and variant sequencing

Most severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) diagnostic tests have relied on RNA extraction followed by reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays. Whereas automation improved logistics and different pooling strategies increased testing capacity, highly multiplexed next-generation sequencing (NGS) diagnostics remain a largely untapped resource. NGS tests have the potential to markedly increase throughput while providing crucial SARS-CoV-2 variant information. Current NGS-based detection and genotyping assays for SARS-CoV-2 are costly, mostly due to parallel sample processing through multiple steps. Here, we have established ApharSeq, in which samples are barcoded in the lysis buffer and pooled before reverse transcription. We validated this assay by applying ApharSeq to more than 500 clinical samples from the Clinical Virology Laboratory at Hadassah hospital in a robotic workflow. The assay was linear across five orders of magnitude, and the limit of detection was Ct 33 (~1000 copies/ml, 95% sensitivity) with >99.5% specificity. ApharSeq provided targeted high-confidence genotype information due to unique molecular identifiers incorporated into this method. Because of early pooling, we were able to estimate a 10- to 100-fold reduction in labor, automated liquid handling, and reagent requirements in high-throughput settings compared to current testing methods. The protocol can be tailored to assay other host or pathogen RNA targets simultaneously. These results suggest that ApharSeq can be a promising tool for current and future mass diagnostic challenges.

Genetic demultiplexing of pooled single-cell RNA-sequencing samples in cancer facilitates effective experimental design

BACKGROUND

Pooling cells from multiple biological samples prior to library preparation within the same single-cell RNA sequencing experiment provides several advantages, including lower library preparation costs and reduced unwanted technological variation, such as batch effects. Computational demultiplexing tools based on natural genetic variation between individuals provide a simple approach to demultiplex samples, which does not require complex additional experimental procedures. However, to our knowledge these tools have not been evaluated in cancer, where somatic variants, which could differ between cells from the same sample, may obscure the signal in natural genetic variation.

RESULTS

Here, we performed in silico benchmark evaluations by combining raw sequencing reads from multiple single-cell samples in high-grade serous ovarian cancer, which has a high copy number burden, and lung adenocarcinoma, which has a high tumor mutational burden. Our results confirm that genetic demultiplexing tools can be effectively deployed on cancer tissue using a pooled experimental design, although high proportions of ambient RNA from cell debris reduce performance.

CONCLUSIONS

This strategy provides significant cost savings through pooled library preparation. To facilitate similar analyses at the experimental design phase, we provide freely accessible code and a reproducible Snakemake workflow built around the best-performing tools found in our in silico benchmark evaluations, available at https://github.com/lmweber/snp-dmx-cancer.

Macrobdella decora: Old World Leech Gut Microbial Community Structure Conserved in a New World Leech

Leeches are found in terrestrial, aquatic, and marine habitats on all continents. Sanguivorous leeches have been used in medicine for millennia. Modern scientific uses include studies of neurons, anticoagulants, and gut microbial symbioses. Hirudo verbana, the European medicinal leech, maintains a gut community dominated by two bacterial symbionts, Aeromonas veronii and Mucinivorans hirudinis, which sometimes account for as much as 97% of the total crop microbiota. The highly simplified gut anatomy and microbiome of H. verbana make it an excellent model organism for studying gut microbial dynamics. The North American medicinal leech, Macrobdella decora, is a hirudinid leech native to Canada and the northern United States. In this study, we show that M. decora symbiont communities are very similar to those in H. verbana. We performed an extensive study using field-caught M. decora and purchased H. verbana from two suppliers. Deep sequencing of the V4 region of the 16S rRNA gene allowed us to determine that the core microbiome of M. decora consists of Bacteroides, Aeromonas, Proteocatella, and Butyricicoccus. The analysis revealed that the compositions of the gut microbiomes of the two leech species were significantly different at all taxonomic levels. The R ² value was highest at the genus and amplicon sequence variant (ASV) levels and much lower at the phylum, class, and order levels. The gut and bladder microbial communities were distinct. We propose that M. decora is an alternative to H. verbana for studies of wild-caught animals and provide evidence for the conservation of digestive-tract and bladder symbionts in annelid models.IMPORTANCE Building evidence implicates the gut microbiome in critical animal functions such as regulating digestion, nutrition, immune regulation, and development. Simplified, phylogenetically diverse models for hypothesis testing are necessary because of the difficulty of assigning causative relationships in complex gut microbiomes. Previous research used Hirudo verbana as a tractable animal model of digestive-tract symbioses. Our data show that Macrobdella decora may work just as well without the drawback of being an endangered organism and with the added advantage of easy access to field-caught specimens. The similarity of the microbial community structures of species from two different continents reveals the highly conserved nature of the microbial symbionts in sanguivorous leeches.

A multiplexed, next generation sequencing platform for high-throughput detection of SARS-CoV-2

Population scale sweeps of viral pathogens, such as SARS-CoV-2, require high intensity testing for effective management. Here, we describe "Systematic Parallel Analysis of RNA coupled to Sequencing for Covid-19 screening" (C19-SPAR-Seq), a multiplexed, scalable, readily automated platform for SARS-CoV-2 detection that is capable of analyzing tens of thousands of patient samples in a single run. To address strict requirements for control of assay parameters and output demanded by clinical diagnostics, we employ a control-based Precision-Recall and Receiver Operator Characteristics (coPR) analysis to assign run-specific quality control metrics. C19-SPAR-Seq coupled to coPR on a trial cohort of several hundred patients performs with a specificity of 100% and sensitivity of 91% on samples with low viral loads, and a sensitivity of >95% on high viral loads associated with disease onset and peak transmissibility. This study establishes the feasibility of employing C19-SPAR-Seq for the large-scale monitoring of SARS-CoV-2 and other pathogens.

Dual indexed library design enables compatibility of in-Drop single-cell RNA-sequencing with exAMP chemistry sequencing platforms

Background

The increasing demand of single-cell RNA-sequencing (scRNA-seq) experiments, such as the number of experiments and cells queried per experiment, necessitates higher sequencing depth coupled to high data quality. New high-throughput sequencers, such as the Illumina NovaSeq 6000, enables this demand to be filled in a cost-effective manner. However, current scRNA-seq library designs present compatibility challenges with newer sequencing technologies, such as index-hopping, and their ability to generate high quality data has yet to be systematically evaluated.

Results

Here, we engineered a dual-indexed library structure, called TruDrop, on top of the inDrop scRNA-seq platform to solve these compatibility challenges, such that TruDrop libraries and standard Illumina libraries can be sequenced alongside each other on the NovaSeq. On scRNA-seq libraries, we implemented a previously-documented countermeasure to the well-described problem of index-hopping, demonstrated significant improvements in base-calling accuracy on the NovaSeq, and provided an example of multiplexing twenty-four scRNA-seq libraries simultaneously. We showed favorable comparisons in transcriptional diversity of TruDrop compared with prior inDrop libraries.

Conclusions

Our approach enables cost-effective, high throughput generation of sequencing data with high quality, which should enable more routine use of scRNA-seq technologies.