Antiviral responses induced by Tdap-IPV vaccination are associated with persistent humoral immunity to Bordetella pertussis.

Many countries continue to experience pertussis epidemics despite widespread vaccination. Waning protection after booster vaccination has highlighted the need for a better understanding of the immunological factors that promote durable protection. Here we apply systems vaccinology to investigate antibody responses in adolescents in the Netherlands (N = 14; NL) and the United Kingdom (N = 12; UK) receiving a tetanus-diphtheria-acellular pertussis-inactivated poliovirus (Tdap-IPV) vaccine. We report that early antiviral and interferon gene expression signatures in blood correlate to persistence of pertussis-specific antibody responses. Single-cell analyses of the innate response identified monocytes and myeloid dendritic cells (MoDC) as principal responders that upregulate antiviral gene expression and type-I interferon cytokine production. With public data, we show that Tdap vaccination stimulates significantly lower antiviral/type-I interferon responses than Tdap-IPV, suggesting that IPV may promote antiviral gene expression. Subsequent in vitro stimulation experiments demonstrate TLR-dependent, IPV-specific activation of the pro-inflammatory p38 MAP kinase pathway in MoDCs. Together, our data provide insights into the molecular host response to pertussis booster vaccination and demonstrate that IPV enhances innate immune activity associated with persistent, pertussis-specific antibody responses.

Primary Exposure to SARS-CoV-2 via Infection or Vaccination Determines Mucosal Antibody-Dependent ACE2 Binding Inhibition.

BACKGROUND: Mucosal antibodies play a critical role in preventing SARS-CoV-2 infections or reinfections by blocking the interaction of the receptor-binding domain (RBD) with the angiotensin-converting enzyme 2 (ACE2) receptor on the cell surface. In this study, we investigated the difference between the mucosal antibody response after primary infection and vaccination. METHODS: We assessed longitudinal changes in the quantity and capacity of nasal antibodies to neutralize the interaction of RBD with the ACE2 receptor using the spike protein and RBD from ancestral SARS-CoV-2 (Wuhan-Hu-1), as well as the RBD from the Delta and Omicron variants. RESULTS: Significantly higher mucosal IgA concentrations were detected postinfection vs postvaccination, while vaccination induced higher IgG concentrations. However, ACE2-inhibiting activity did not differ between the cohorts. Regarding whether IgA or IgG drove ACE2 inhibition, infection-induced binding inhibition was driven by both isotypes, while postvaccination binding inhibition was mainly driven by IgG. CONCLUSIONS: Our study provides new insights into the relationship between antibody isotypes and neutralization by using a sensitive and high-throughput ACE2 binding inhibition assay. Key differences are highlighted between vaccination and infection at the mucosal level, showing that despite differences in the response quantity, postinfection and postvaccination ACE2 binding inhibition capacity did not differ.

CD7 activation regulates cytotoxicity-driven pathology in systemic sclerosis, yielding a target for selective cell depletion.

OBJECTIVES: Cytotoxic T cells and natural killer (NK) cells are central effector cells in cancer and infections. Their effector response is regulated by activating and inhibitory receptors. The regulation of these cells in systemic autoimmune diseases such as systemic sclerosis (SSc) is less defined. METHODS: We conducted ex vivo analysis of affected skin and blood samples from 4 SSc patient cohorts (a total of 165 SSc vs 80 healthy individuals) using single-cell transcriptomics, flow cytometry and multiplex immunofluorescence staining. We further analysed the effects of costimulatory modulation in functional assays, and in a severely affected SSc patient who was treated on compassionate use with a novel anti-CD3/CD7 immunotoxin treatment. RESULTS: Here, we show that SSc-affected skin contains elevated numbers of proliferating T cells, cytotoxic T cells and NK cells. These cells selectively express the costimulatory molecule CD7 in association with cytotoxic, proinflammatory and profibrotic genes, especially in recent-onset and severe disease. We demonstrate that CD7 regulates the cytolytic activity of T cells and NK cells and that selective depletion of CD7+ cells prevents cytotoxic cell-induced fibroblast contraction and inhibits their profibrotic phenotype. Finally, anti-CD3/CD7 directed depletive treatment eliminated CD7+ skin cells and stabilised disease manifestations in a severely affected SSc patient. CONCLUSION: Together, the findings imply costimulatory molecules as key regulators of cytotoxicity-driven pathology in systemic autoimmune disease, yielding CD7 as a novel target for selective depletion of pathogenic cells.

Quantifying the deformability of malaria-infected red blood cells using deep learning trained on synthetic cells.

Several hematologic diseases, including malaria, diabetes, and sickle cell anemia, result in a reduced red blood cell deformability. This deformability can be measured using a microfluidic device with channels of varying width. Nevertheless, it is challenging to algorithmically recognize large numbers of red blood cells and quantify their deformability from image data. Deep learning has become the method of choice to handle noisy and complex image data. However, it requires a significant amount of labeled data to train the neural networks. By creating images of cells and mimicking noise and plasticity in those images, we generate synthetic data to train a network to detect and segment red blood cells from video-recordings, without the need for manually annotated labels. Using this new method, we uncover significant differences between the deformability of RBCs infected with different strains of Plasmodium falciparum, providing clues to the variation in virulence of these strains.

Effect of age and season on respiratory mucosal immune marker profiles.

BACKGROUND: The upper respiratory tract is continuously exposed to microorganisms and noxious elements, leading to local immune responses and the secretion of immune markers. While several studies describe immune marker profiles in respiratory mucosal samples in defined patient cohorts, mucosal immune profiles from the general population during the different seasons are lacking. Such baseline profiles are essential to understand the effect of various exposures to the mucosal immune system throughout life. OBJECTIVE: We sought to establish baseline local upper respiratory mucosal immune profiles in the general population and assess these profiles with regard to age, sex, seasonality, and basic health and lifestyle factors. METHODS: We measured the concentrations of 35 immune markers involved in a broad range of immunological processes at the mucosa in nasopharyngeal swab samples from 951 individuals, aged 0 to 86 years, from a nationwide study. RESULTS: Clustering analysis showed that immune marker profiles clearly reflected immunological functions, such as tissue regeneration and antiviral responses. Immune marker concentrations changed strongly with seasonality and age, with the most profound changes occurring in the first 25 years of life; they were also associated with sex, body mass index, smoking, mild symptoms of airway infection, and chronic asthma and hay fever. CONCLUSION: Immunological analyses of noninvasive mucosal samples provide insight into mucosal immune responses to microbial and noxious element exposure in the general population. These data provide a baseline for future studies on respiratory mucosal immune responses and for the development of mucosal immune-based diagnostics.

Longitudinal analysis on the ecological dynamics of the cervicovaginal microbiome in hrHPV infection.

The cervicovaginal microbiome (CVM) is a dynamic continuous microenvironment that can be clustered in microbial community state types (CSTs) and is associated with women's cervical health. Lactobacillus-depleted communities particularly associate with an increased susceptibility for persistence of high-risk human papillomavirus (hrHPV) infections and progression of disease, but the long-term ecological dynamics of CSTs after hrHPV infection diagnosis remain poorly understood. To determine such dynamics, we examined the CVM of our longitudinal cohort of 141 women diagnosed with hrHPV infection at baseline with collected cervical smears at two timepoints six-months apart. Here we describe that the long-term microbiome dissimilarity has a positive correlation with microbial diversity at both visits and that women with high abundance and dominance for Lactobacillus iners at baseline exhibit more similar microbiome composition at second visit than women with Lactobacillus-depleted communities at baseline. We further show that the species Lactobacillus acidophilus and Megasphaera genomosp type 1 associate with CST changes between both visits. Lastly, we also observe that Gardnerella vaginalis is associated with the stability of Lactobacillus-depleted communities while L. iners is associated with the instability of Megasphaera genomosp type 1-dominated communities. Our data suggest dynamic patterns of cervicovaginal CSTs during hrHPV infection, which could be potentially used to develop microbiome-based therapies against infection progression towards disease.

Comparative Clustering (CompaCt) of eukaryote complexomes identifies novel interactions and sheds light on protein complex evolution.

Complexome profiling allows large-scale, untargeted, and comprehensive characterization of protein complexes in a biological sample using a combined approach of separating intact protein complexes e.g., by native gel electrophoresis, followed by mass spectrometric analysis of the proteins in the resulting fractions. Over the last decade, its application has resulted in a large collection of complexome profiling datasets. While computational methods have been developed for the analysis of individual datasets, methods for large-scale comparative analysis of complexomes from multiple species are lacking. Here, we present Comparative Clustering (CompaCt), that performs fully automated integrative analysis of complexome profiling data from multiple species, enabling systematic characterization and comparison of complexomes. CompaCt implements a novel method for leveraging orthology in comparative analysis to allow systematic identification of conserved as well as taxon-specific elements of the analyzed complexomes. We applied this method to a collection of 53 complexome profiles spanning the major branches of the eukaryotes. We demonstrate the ability of CompaCt to robustly identify the composition of protein complexes, and show that integrated analysis of multiple datasets improves characterization of complexes from specific complexome profiles when compared to separate analyses. We identified novel candidate interactors and complexes in a number of species from previously analyzed datasets, like the emp24, the V-ATPase and mitochondrial ATP synthase complexes. Lastly, we demonstrate the utility of CompaCt for the automated large-scale characterization of the complexome of the mosquito Anopheles stephensi shedding light on the evolution of metazoan protein complexes. CompaCt is available from https://github.com/cmbi/compact-bio.

Uncharacterized protein C17orf80 - a novel interactor of human mitochondrial nucleoids.

Molecular functions of many human proteins remain unstudied, despite the demonstrated association with diseases or pivotal molecular structures, such as mitochondrial DNA (mtDNA). This small genome is crucial for the proper functioning of mitochondria, the energy-converting organelles. In mammals, mtDNA is arranged into macromolecular complexes called nucleoids that serve as functional stations for its maintenance and expression. Here, we aimed to explore an uncharacterized protein C17orf80, which was previously detected close to the nucleoid components by proximity labelling mass spectrometry. To investigate the subcellular localization and function of C17orf80, we took advantage of immunofluorescence microscopy, interaction proteomics and several biochemical assays. We demonstrate that C17orf80 is a mitochondrial membrane-associated protein that interacts with nucleoids even when mtDNA replication is inhibited. In addition, we show that C17orf80 is not essential for mtDNA maintenance and mitochondrial gene expression in cultured human cells. These results provide a basis for uncovering the molecular function of C17orf80 and the nature of its association with nucleoids, possibly leading to new insights about mtDNA and its expression.

Understanding structure-guided variant effect predictions using 3D convolutional neural networks.

Predicting pathogenicity of missense variants in molecular diagnostics remains a challenge despite the available wealth of data, such as evolutionary information, and the wealth of tools to integrate that data. We describe DeepRank-Mut, a configurable framework designed to extract and learn from physicochemically relevant features of amino acids surrounding missense variants in 3D space. For each variant, various atomic and residue-level features are extracted from its structural environment, including sequence conservation scores of the surrounding amino acids, and stored in multi-channel 3D voxel grids which are then used to train a 3D convolutional neural network (3D-CNN). The resultant model gives a probabilistic estimate of whether a given input variant is disease-causing or benign. We find that the performance of our 3D-CNN model, on independent test datasets, is comparable to other widely used resources which also combine sequence and structural features. Based on the 10-fold cross-validation experiments, we achieve an average accuracy of 0.77 on the independent test datasets. We discuss the contribution of the variant neighborhood in the model's predictive power, in addition to the impact of individual features on the model's performance. Two key features: evolutionary information of residues in the variant neighborhood and their solvent accessibilities were observed to influence the predictions. We also highlight how predictions are impacted by the underlying disease mechanisms of missense mutations and offer insights into understanding these to improve pathogenicity predictions. Our study presents aspects to take into consideration when adopting deep learning approaches for protein structure-guided pathogenicity predictions.

Parallel evolution and enhanced virulence upon in vivo passage of an RNA virus in Drosophila melanogaster.

Virus evolution is strongly affected by antagonistic co-evolution of virus and host. Host immunity positively selects for viruses that evade the immune response, which in turn may drive counter-adaptations in host immune genes. We investigated how host immune pressure shapes virus populations, using the fruit fly Drosophila melanogaster and its natural pathogen Drosophila C virus (DCV), as a model. We performed an experimental evolution study in which DCV was serially passaged for ten generations in three fly genotypes differing in their antiviral RNAi response: wild-type flies and flies in which the endonuclease gene Dicer-2 was either overexpressed or inactivated. All evolved virus populations replicated more efficiently in vivo and were more virulent than the parental stock. The number of polymorphisms increased in all three host genotypes with passage number, which was most pronounced in Dicer-2 knockout flies. Mutational analysis showed strong parallel evolution, as mutations accumulated in a specific region of the VP3 capsid protein in every lineage in a host genotype-independent manner. The parental tyrosine at position ninety-five of VP3 was substituted with either one of five different amino acids in fourteen out of fifteen lineages. However, no consistent amino acid changes were observed in the viral RNAi suppressor gene 1A, nor elsewhere in the genome in any of the host backgrounds. Our study indicates that the RNAi response restricts the sequence space that can be explored by viral populations. Moreover, our study illustrates how evolution towards higher virulence can be a highly reproducible, yet unpredictable process.

Prior exposure to B. pertussis shapes the mucosal antibody response to acellular pertussis booster vaccination.

Bordetella pertussis (Bp), the causative agent of pertussis, continues to circulate despite widespread vaccination programs. An important question is whether and how (sub)clinical infections shape immune memory to Bp, particularly in populations primed with acellular pertussis vaccines (aP). Here, we examine the prevalence of mucosal antibodies against non-vaccine antigens in aP-primed children and adolescents of the BERT study (NCT03697798), using antibody binding to a Bp mutant strain lacking aP antigens (Bp_mut). Our study identifies increased levels of mucosal IgG and IgA binding to Bp_mut in older aP-primed individuals, suggesting different Bp exposure between aP-primed birth cohorts, in line with pertussis disease incidence data. To examine whether Bp exposure influences vaccination responses, we measured mucosal antibody responses to aP booster vaccination as a secondary study outcome. Although booster vaccination induces significant increases in mucosal antibodies to Bp in both cohorts, the older age group that had higher baseline antibodies to Bp_ mut shows increased persistence of antibodies after vaccination.

In-depth insights into cervicovaginal microbial communities and hrHPV infections using high-resolution microbiome profiling.

The cervicovaginal microbiome (CVM) correlates with women's cervical health, and variations in its composition are associated with high-risk human papillomavirus (hrHPV) infection outcomes. Cervicovaginal microbes have been grouped into five community state types (CSTs) based on microbial community composition and abundance. However, studying the impact of CSTs in health and disease is challenging because the current sequencing technologies have limited confident discrimination between closely related and yet functionally different bacterial species. Circular probe-based RNA sequencing (ciRNAseq) achieves high-resolution microbiome profiling and therefore provides in-depth and unambiguous knowledge about the composition of the CVM. Based on ciRNAseq profiling of a large cohort of cervical smears (n = 541), we here define subgroups of CSTs I, III, and IV based on intra-CST differences with respect to abundances of Lactobacillus acidophilus (CSTs I-A vs. I-B and CSTs III-A vs. III-B), Lactobacillus iners (CSTs I-A vs. I-B and CSTs III-A vs. III-B), and Megasphaera genomosp type 1 (CSTs IV-A vs. IV-B). Our results further support the existence of subgroups of CST IV-C that are dominant for non-Lactobacillus species and have intermediate microbial diversity. We also show that CST V is associated with uninfected conditions, and CST IV-A associates with hrHPV-induced cervical disease. In conclusion, we characterized new subdivisions of cervicovaginal CSTs, which may further advance our understanding of women's cervical health and hrHPV-related progression to disease.

Assessing the Cervicovaginal Microbiota in the Context of hrHPV Infections: Temporal Dynamics and Therapeutic Strategies.

Cervical cancer is the third leading cause of female cancers globally, resulting in more than 300,000 deaths every year. The majority of all cervical cancers are caused by persistent infections with high-risk human papillomaviruses (hrHPV) that can progress to cancer via a series of premalignant lesions. Most women, however, clear this infection within a year, concomitant with disease regression. Both hrHPV clearance and disease regression have been associated with the composition of the cervicovaginal microenvironment, which is defined by the host immune system and the cervicovaginal microbiome (CVM). A healthy microbiome is generally characterized by a high abundance of Lactobacillus species, and a change in the composition may cause bacterial vaginosis (BV), which is associated with an increased susceptibility to persistent hrHPV infections and disease. In this review, the composition of the CVM is discussed, with emphasis on the possible causes that drive changes in the cervicovaginal microbiota in relation to hrHPV infections, disease progression, and disease regression. The literature search focused on the composition of the CVM and its correlation with hrHPV infections and neoplastic lesions as well as the current efforts to adjust the microbiome against adverse viral outcomes.

Nutrition or nature: using elementary flux modes to disentangle the complex forces shaping prokaryote pan-genomes.

BACKGROUND: Microbial pan-genomes are shaped by a complex combination of stochastic and deterministic forces. Even closely related genomes exhibit extensive variation in their gene content. Understanding what drives this variation requires exploring the interactions of gene products with each other and with the organism's external environment. However, to date, conceptual models of pan-genome dynamics often represent genes as independent units and provide limited information about their mechanistic interactions. RESULTS: We simulated the stochastic process of gene-loss using the pooled genome-scale metabolic reaction networks of 46 taxonomically diverse bacterial and archaeal families as proxies for their pan-genomes. The frequency by which reactions are retained in functional networks when stochastic gene loss is simulated in diverse environments allowed us to disentangle the metabolic reactions whose presence depends on the metabolite composition of the external environment (constrained by "nutrition") from those that are independent of the environment (constrained by "nature"). By comparing the frequency of reactions from the first group with their observed frequencies in bacterial and archaeal families, we predicted the metabolic niches that shaped the genomic composition of these lineages. Moreover, we found that the lineages that were shaped by a more diverse metabolic niche also occur in more diverse biomes as assessed by global environmental sequencing datasets. CONCLUSION: We introduce a computational framework for analyzing and interpreting pan-reactomes that provides novel insights into the ecological and evolutionary drivers of pan-genome dynamics.

Targeted RNA next generation sequencing analysis of cervical smears can predict the presence of hrHPV-induced cervical lesions.

BACKGROUND: Because most cervical cancers are caused by high-risk human papillomaviruses (hrHPVs), cervical cancer prevention programs increasingly employ hrHPV testing as a primary test. The high sensitivity of HPV tests is accompanied by low specificity, resulting in high rates of overdiagnosis and overtreatment. Targeted circular probe-based RNA next generation sequencing (ciRNAseq) allows for the quantitative detection of RNAs of interest with high sequencing depth. Here, we examined the potential of ciRNAseq-testing on cervical scrapes to identify hrHPV-positive women at risk of having or developing high-grade cervical intraepithelial neoplasia (CIN). METHODS: We performed ciRNAseq on 610 cervical scrapes from the Dutch cervical cancer screening program to detect gene expression from 15 hrHPV genotypes and from 429 human genes. Differentially expressed hrHPV- and host genes in scrapes from women with outcome "no CIN" or "CIN2+" were identified and a model was built to distinguish these groups. RESULTS: Apart from increasing percentages of hrHPV oncogene expression from "no CIN" to high-grade cytology/histology, we identified genes involved in cell cycle regulation, tyrosine kinase signaling pathways, immune suppression, and DNA repair being expressed at significantly higher levels in scrapes with high-grade cytology and histology. Machine learning using random forest on all the expression data resulted in a model that detected 'no CIN' versus CIN2+ in an independent data set with sensitivity and specificity of respectively 85 ± 8% and 72 ± 13%. CONCLUSIONS: CiRNAseq on exfoliated cells in cervical scrapes measures hrHPV-(onco)gene expression and host gene expression in one single assay and in the process identifies HPV genotype. By combining these data and applying machine learning protocols, the risk of CIN can be calculated. Because ciRNAseq can be performed in high-throughput, making it cost-effective, it can be a promising screening technology to stratify women at risk of CIN2+. Further increasing specificity by model improvement in larger cohorts is warranted.

CG7630 is the Drosophila melanogaster homolog of the cytochrome c oxidase subunit COX7B.

The mitochondrial respiratory chain (MRC) is composed of four multiheteromeric enzyme complexes. According to the endosymbiotic origin of mitochondria, eukaryotic MRC derives from ancestral proteobacterial respiratory structures consisting of a minimal set of complexes formed by a few subunits associated with redox prosthetic groups. These enzymes, which are the "core" redox centers of respiration, acquired additional subunits, and increased their complexity throughout evolution. Cytochrome c oxidase (COX), the terminal component of MRC, has a highly interspecific heterogeneous composition. Mammalian COX consists of 14 different polypeptides, of which COX7B is considered the evolutionarily youngest subunit. We applied proteomic, biochemical, and genetic approaches to investigate the COX composition in the invertebrate model Drosophila melanogaster. We identified and characterized a novel subunit which is widely different in amino acid sequence, but similar in secondary and tertiary structures to COX7B, and provided evidence that this object is in fact replacing the latter subunit in virtually all protostome invertebrates. These results demonstrate that although individual structures may differ the composition of COX is functionally conserved between vertebrate and invertebrate species.

Mitochondrial complex complexification.

Variation in complex composition provides clues about the function of individual subunits.

BCG-induced trained immunity enhances acellular pertussis vaccination responses in an explorative randomized clinical trial.

Acellular pertussis (aP) booster vaccines are central to pertussis immunization programs, although their effectiveness varies. The Bacille Calmette-Guérin (BCG) vaccine is a prototype inducer of trained immunity, which enhances immune responses to subsequent infections or vaccinations. While previous clinical studies have demonstrated that trained immunity can protect against heterologous infections, its effect on aP vaccines in humans is unknown. We conducted a clinical study in order to determine the immunological effects of trained immunity on pertussis vaccination. Healthy female volunteers were randomly assigned to either receive BCG followed by a booster dose of tetanus-diphteria-pertussis inactivated polio vaccine (Tdap-IPV) 3 months later (BCG-trained), BCG + Tdap-IPV concurrently, or Tdap-IPV followed by BCG 3 months later. Primary outcomes were pertussis-specific humoral, T- and B-cell responses and were quantified at baseline of Tdap-IPV vaccination and 2 weeks thereafter. As a secondary outcome in the BCG-trained cohort, ex vivo leukocyte responses were measured in response to unrelated stimuli before and after BCG vaccination. BCG vaccination 3 months prior to, but not concurrent with, Tdap-IPV improves pertussis-specific Th1-cell and humoral responses, and also increases total memory B cell responses. These responses were correlated with enhanced IL-6 and IL-1ß production at the baseline of Tdap-IPV vaccination in the BCG-trained cohort. Our study demonstrates that prior BCG vaccination potentiates immune responses to pertussis vaccines and that biomarkers of trained immunity are the most reliable correlates of those responses.

Mitochondrial RNA processing defect caused by a SUPV3L1 mutation in two siblings with a novel neurodegenerative syndrome.

SUPV3L1 encodes a helicase that is mainly localized in the mitochondria. It has been shown in vitro to possess both double-stranded RNA and DNA unwinding activity that is ATP-dependent. Here we report the first two patients for this gene who presented with a homozygous preliminary stop codon resulting in a C-terminal truncation of the SUPV3L1 protein. They presented with a characteristic phenotype of neurodegenerative nature with progressive spastic paraparesis, growth restriction, hypopigmentation, and predisposition to autoimmune disease. Ophthalmological examination showed severe photophobia with corneal erosions, optic atrophy, and pigmentary retinopathy, while neuroimaging showed atrophy of the optic chiasm and the pons with calcification of putamina, with intermittent and mild elevation of lactate. We show that the amino acids that are eliminated by the preliminary stop codon are highly conserved and are predicted to form an amphipathic helix. To investigate if the mutation causes mitochondrial dysfunction, we examined fibroblasts of the proband. We observed very low expression of the truncated protein, a reduction in the mature ND6 mRNA species as well as the accumulation of double-stranded RNA. Lentiviral complementation with the full-length SUPV3L1 cDNA partly restored the observed RNA phenotypes, supporting that the SUPV3L1 mutation in these patients is pathogenic and the cause of the disease.

Novel high-resolution targeted sequencing of the cervicovaginal microbiome.

BACKGROUND: The cervicovaginal microbiome (CVM) plays a significant role in women's cervical health and disease. Microbial alterations at the species level and characteristic community state types (CST) have been associated with acquisition and persistence of high-risk human papillomavirus (hrHPV) infections that may result in progression of cervical lesions to malignancy. Current sequencing methods, especially most commonly used multiplex 16S rRNA gene sequencing, struggle to fully clarify these changes because they generally fail to provide sufficient taxonomic resolution to adequately perform species-level associative studies. To improve CVM species designation, we designed a novel sequencing tool targeting microbes at the species taxonomic rank and examined its potential for profiling the CVM. RESULTS: We introduce an accessible and practical circular probe-based RNA sequencing (CiRNAseq) technology with the potential to profile and quantify the CVM. In vitro and in silico validations demonstrate that CiRNAseq can distinctively detect species in a mock mixed microbial environment, with the output data reflecting its ability to estimate microbes' abundance. Moreover, compared to 16S rRNA gene sequencing, CiRNAseq provides equivalent results but with improved sequencing sensitivity. Analyses of a cohort of cervical smears from hrHPV-negative women versus hrHPV-positive women with high-grade cervical intraepithelial neoplasia confirmed known differences in CST occurring in the CVM of women with hrHPV-induced lesions. The technique also revealed variations in microbial diversity and abundance in the CVM of hrHPV-positive women when compared to hrHPV-negative women. CONCLUSIONS: CiRNAseq is a promising tool for studying the interplay between the CVM and hrHPV in cervical carcinogenesis. This technology could provide a better understanding of cervicovaginal CST and microbial species during health and disease, prompting the discovery of biomarkers, additional to hrHPV, that can help detect high-grade cervical lesions.