Search

Summary: The elucidation of the role of microorganisms in human infections has been hindered by difficulties using conventional culture-based techniques. Here, we present a protocol for the investigation of transcriptionally active microbes (TAMs) using an RNA sequencing (RNA-seq)-based approach. We describe the steps for RNA isolation, viral genome sequencing, RNA-seq library preparation, and metatranscriptomic and transcriptomic analysis. This protocol permits a comprehensive evaluation of TAMs’ contributions to the differential severity of infectious diseases, with a particular focus on diseases such as COVID-19.For complete details on the use and execution of this protocol, please refer to Devi et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Abstract Dengue virus (DENV), known to cause viral infection, belongs to the family Flaviviridae, having four serotypes (DENV1-4) that spreads by the bite of the Aedes aegypti mosquito. India has been suffering from dengue outbreaks annually with widespread epidemics by prevalence of all the four DENV serotypes. The diverse spectrum of clinical manifestations in dengue infection, mild to severe forms, makes the need of timely diagnosis and prompt treatment an essence. The identification of a dengue host response signature in serum can increase the understanding of dengue pathogenesis since most dengue NS1 Ag tests have been developed and evaluated in serum samples. Here, to understand the same, we undertook a dual RNA-sequencing (RNA-Seq) based approach from the serum samples of dengue-infected patients. The results thus yield the early transcriptional signatures that discriminated the high viral reads patients from patients who had low dengue viral reads. We identified a significant upregulation of two sets of genes, key antiviral (IFIT3, RSAD2, SAT1) and vascular dysfunction (TNFS10, CXCL8) related genes in the high viral reads group. Deeper delving of this gene profile revealed a unique two-way response, where the antiviral genes can mediate the disease course to mild, contrarily the increased expression of the other gene set might act as pointers of severe disease course. Further, we explored the hematologic parameters from the complete blood count (CBC), which suggests that lymphocytes (low) and neutrophils (high) might serve as an early predictor of prognosis in dengue infection. Collectively, our findings give insights into the foundation for further investigation of the early host response using the RNA isolated from dengue patients’ serum samples and opens the door for careful monitoring of the early clinical and transcriptome profiles for management of the dengue patients.

IntroductionDengue virus (DENV) is a flavivirus that has emerged as a global health threat, characterized by either asymptomatic or mild self-limiting febrile illness, but a subset of DENV outbreaks have been associated with severe disease. Studies have looked into the host immune response and dengue viral load during infection. However, it remains unknown how the active microbial isolates modulate the dengue viral infection. In this study, we demonstrate the significance of in-depth analysis of microbiota composition in the serum samples of dengue-infected patients.Materials and methodsRNA was extracted from the serum samples collected from 24 dengue positive patients. The human mapped reads generated through RNA-Sequencing (RNA-Seq) were removed, while the unmapped (non-human) reads were employed for microbial taxonomic classification using Kraken2 and Bracken2. Further, we assessed the initial blood parameters analyzing the complete blood count (CBC) profile of the patients.ResultsFindings revealed differential abundance of commensals and pathogenic microbes in the early febrile period of hospitalized dengue patients, segregated into, High Viral Reads (HVR) and Low Viral Reads (LVR). The Campylobacter genus was abundant in the HVR whereas Lactobacillus dominated the LVR patients. At species level, the microbiota of HVR exhibited higher abundance of unique potential opportunistic microbes, compared to the commensal microbes’ enrichment in the LVR patients’. We hypothesize that the DENV might alter the microbiota composition as observed by the increase in preponderance of opportunistic pathogens and an absence of commensals in the HVR. The presence of commensals in the LVR might explain, i) overall lower dengue viral reads compared to the HVR, and ii) shift in lymphocytes (high) and neutrophils (low) counts; resulting in a comparatively milder clinical manifestation in this group. Our findings may help in understanding the co-infection aspect that will be important to develop dengue therapeutics and vaccines.DiscussionThis study highlights the potential of the unexplored roles of the TAMs in modulating the dengue disease severity using the metatranscriptomic sequencing. This study serves to enhance our understanding of the distinctive microbial and hematologic signatures in the early infection stage that differentiate patients with high viral reads patients from those with low dengue viral reads.

IntroductionSeveral efforts have been made to describe the complexity of T cell heterogeneity during the COVID-19 disease; however, there remain gaps in our understanding in terms of the granularity within.MethodsFor this attempt, we performed a single-cell transcriptomic analysis of 33 individuals (4 healthy, 16 COVID-19 positive patients, and 13 COVID-19 recovered individuals).ResultsWe found CD8+ T cell-biased lymphopenia in COVID-19 patients compared to healthy and recovered individuals. We also found an optimal Th1/Th2 ratio, indicating an effective immune response during COVID-19. Expansion of activated CD4+ T and NK T was detected in the COVID-19-positive individuals. Surprisingly, we found cellular and metal ion homeostasis pathways enriched in the COVID-19-positive individuals compared to the healthy and recovered in the CD8+ T cell populations (CD8+ TCM and CD8+ TEM) as well as activated CD4+ T cells.DiscussionIn summary, the COVID-19-positive individuals exhibit a dynamic T cell mediated response. This response may have a possible association with the dysregulation of non-canonical pathways, including housekeeping functions in addition to the conventional antiviral immune response mediated by the T cell subpopulation. These findings considerably extend our insights into the heterogeneity of T cell response during and post-SARS-CoV-2 infection.

Summary: Emergence of new SARS-CoV-2 VOCs jeopardize global vaccine and herd immunity safeguards. VOCs interactions with host microbiota might affect clinical course and outcome. This longitudinal investigation involving Pre-VOC and VOCs (Delta & Omicron) holo-transcriptome based nasopharyngeal microbiome at taxonomic levels followed by metabolic pathway analysis and integrative host-microbiome interaction. VOCs showed enrichment of Proteobacteria with dominance of Pseudomonas. Interestingly, Proteobacteria with superiority of Pseudomonas and Acinetobacter, were highlights of Delta VOC rather than Omicron. Common species comprising the core microbiome across all variants, reiterated the significance of Klebsiella pneumoniae in Delta, and its association with metabolic pathways enhancing inflammation in patients. Microbe-host gene correlation network revealed Acinetobacter baumannii, Pseudomonas stutzeri, and Pseudomonas aeuroginosa modulating immune pathways, which might augment clinical severity in Delta. Importantly, opportunistic species of Acinetobacter, Enterococcus, Prevotella, and Streptococcus were abundant in Delta-mortality. The study establishes a functional association between elevated nasal pathobionts and dysregulated host response, particularly for Delta.

The Omicron variant of SARS-CoV-2 is capable of infecting unvaccinated, vaccinated and previously-infected individuals due to its ability to evade neutralization by antibodies. With multiple sub-lineages of Omicron emerging in the last 12 months, there is inadequate information on the quantitative antibody response generated upon natural infection with Omicron variant and whether these antibodies offer cross-protection against other sub-lineages of Omicron variant. In this study, we characterized the growth kinetics of Kappa, Delta and Omicron variants of SARS-CoV-2 in Calu-3 cells. Relatively higher amounts infectious virus titers, cytopathic effect and disruption of epithelial barrier functions was observed with Delta variant whereas infection with Omicron sub-lineages led to a more robust induction of interferon pathway, lower level of virus replication and mild effect on epithelial barrier. The replication kinetics of BA.1, BA.2 and BA.2.75 sub-lineages of the Omicron variant were comparable in cell culture and natural infection in a subset of individuals led to a significant increase in binding and neutralizing antibodies to the Delta variant and all the three sub-lineages of Omicron but the level of neutralizing antibodies were lowest against the BA.2.75 variant. Finally, we show that Cu2+, Zn2+ and Fe2+ salts inhibited in vitro RdRp activity but only Cu2+ and Fe2+ inhibited both the Delta and Omicron variants in cell culture. Thus, our results suggest that high levels of interferons induced upon infection with Omicron variant may counter virus replication and spread. Waning neutralizing antibody titers rendered subjects susceptible to infection by Omicron variants and natural Omicron infection elicits neutralizing antibodies that can cross-react with other sub-lineages of Omicron and other variants of concern.

The development of COVID 19 vaccines as an effort to mitigate the outbreak, has saved millions of lives globally. However, vaccination breakthroughs have continuously challenged the vaccines' effectiveness and provided incentives to explore facets holding potential to alter vaccination-induced immunity and protection from subsequent infection, especially VOCs (Variants Of Concern). We explored the functional dynamics of nasopharyngeal transcriptionally active microbes (TAMs) between vaccination breakthroughs and unvaccinated SARS-CoV-2 infected individuals. Microbial taxonomic communities were differentially altered with skewed enrichment of bacterial class/genera of Firmicutes and Gammaproteobacteria with grossly reduced phylum Bacteroidetes in vaccination breakthrough individuals. The Bacillus genus was abundant in Firmicutes in vaccination breakthrough whereas Prevotella among Bacteroides dominated the unvaccinated. Also, Pseudomonas and Salmonella of Gammaproteobacteria were overrepresented in vaccination breakthrough, whilst unvaccinated showed presence of several genera, Achromobacter, Bordetella, Burkholderia, Neisseria, Hemophilus, Salmonella and Pseudomonas, belonging to Proteobacteria. At species level, the microbiota of vaccination breakthrough exhibited relatively higher abundance of unique commensals, in comparison to potential opportunistic microbes enrichment in unvaccinated patients' microbiota. Functional metabolic pathways like amino acid biosynthesis, sulphate assimilation, fatty acid and beta oxidation, associated with generation of SCFAs (short chain fatty acids), were enriched in vaccination breakthroughs. Majorly, metabolic pathways of LCFAs biosynthesis (long chain fatty acids; oleate, dodecenoate, palmitoleate, gondoate) were found associated with the unvaccinated. Our research highlights that vaccination decreases the microbial diversity in terms of depleting opportunistic pathogens and increasing the preponderance of commensals with respect to unvaccinated patients. Metabolic pathway analysis substantiates the shift in diversity to functionally modulate immune response generation, which may be related to mild clinical manifestations and faster recovery times during vaccination breakthroughs.

IntroductionDespite numerous efforts to describe COVID-19's immunological landscape, there is still a gap in our understanding of the virus's infections after-effects, especially in the recovered patients. This would be important to understand as we now have huge number of global populations infected by the SARS-CoV-2 as well as variables inclusive of VOCs, reinfections, and vaccination breakthroughs. Furthermore, single-cell transcriptome alone is often insufficient to understand the complex human host immune landscape underlying differential disease severity and clinical outcome.MethodsBy combining single-cell multi-omics (Whole Transcriptome Analysis plus Antibody-seq) and machine learning-based analysis, we aim to better understand the functional aspects of cellular and immunological heterogeneity in the COVID-19 positive, recovered and the healthy individuals.ResultsBased on single-cell transcriptome and surface marker study of 163,197 cells (124,726 cells after data QC) from the 33 individuals (healthy=4, COVID-19 positive=16, and COVID-19 recovered=13), we observed a reduced MHC Class-I-mediated antigen presentation and dysregulated MHC Class-II-mediated antigen presentation in the COVID-19 patients, with restoration of the process in the recovered individuals. B-cell maturation process was also impaired in the positive and the recovered individuals. Importantly, we discovered that a subset of the naive T-cells from the healthy individuals were absent from the recovered individuals, suggesting a post-infection inflammatory stage. Both COVID-19 positive patients and the recovered individuals exhibited a CD40-CD40LG-mediated inflammatory response in the monocytes and T-cell subsets. T-cells, NK-cells, and monocyte-mediated elevation of immunological, stress and antiviral responses were also seen in the COVID-19 positive and the recovered individuals, along with an abnormal T-cell activation, inflammatory response, and faster cellular transition of T cell subtypes in the COVID-19 patients. Importantly, above immune findings were used for a Bayesian network model, which significantly revealed FOS, CXCL8, IL1β, CST3, PSAP, CD45 and CD74 as COVID-19 severity predictors.DiscussionIn conclusion, COVID-19 recovered individuals exhibited a hyper-activated inflammatory response with the loss of B cell maturation, suggesting an impeded post-infection stage, necessitating further research to delineate the dynamic immune response associated with the COVID-19. To our knowledge this is first multi-omic study trying to understand the differential and dynamic immune response underlying the sample subtypes.

ABSTRACT The modulators of severe COVID-19 have emerged as the most intriguing features of SARS-CoV-2 pathogenesis. This is especially true as we are encountering variants of concern (VOC) with increased transmissibility and vaccination breakthroughs. Microbial co-infections are being investigated as one of the crucial factors for exacerbation of disease severity and complications of COVID-19. A key question remains whether early transcriptionally active microbial signature/s in COVID-19 patients can provide a window for future disease severity susceptibility and outcome? Using complementary metagenomics sequencing approaches, respiratory virus oligo panel (RVOP) and Holo-seq, our study highlights the possible functional role of nasopharyngeal early resident transcriptionally active microbes in modulating disease severity, within recovered patients with sub-phenotypes (mild, moderate, severe) and mortality. The integrative analysis combines patients’ clinical parameters, SARS-CoV-2 phylogenetic analysis, microbial differential composition, and their functional role. The clinical sub-phenotypes analysis led to the identification of transcriptionally active bacterial species associated with disease severity. We found significant transcript abundance of Achromobacter xylosoxidans and Bacillus cereus in the mortality, Leptotrichia buccalis in the severe, Veillonella parvula in the moderate, and Actinomyces meyeri and Halomonas sp. in the mild COVID-19 patients. Additionally, the metabolic pathways, distinguishing the microbial functional signatures between the clinical sub-phenotypes, were also identified. We report a plausible mechanism wherein the increased transcriptionally active bacterial isolates might contribute to enhanced inflammatory response and co-infections that could modulate the disease severity in these groups. Current study provides an opportunity for potentially using these bacterial species for screening and identifying COVID-19 patient sub-groups with severe disease outcome and priority medical care. IMPORTANCE COVID-19 is invariably a disease of diverse clinical manifestation, with multiple facets involved in modulating the progression and outcome. In this regard, we investigated the role of transcriptionally active microbial co-infections as possible modulators of disease pathology in hospital admitted SARS-CoV-2 infected patients. Specifically, can there be early nasopharyngeal microbial signatures indicative of prospective disease severity? Based on disease severity symptoms, the patients were segregated into clinical sub-phenotypes: mild, moderate, severe (recovered), and mortality. We identified significant presence of transcriptionally active isolates, Achromobacter xylosoxidans and Bacillus cereus in the mortality patients. Importantly, the bacterial species might contribute toward enhancing the inflammatory responses as well as reported to be resistant to common antibiotic therapy, which together hold potential to alter the disease severity and outcome.

Vaccine development against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been of primary importance to contain the ongoing global pandemic. However, studies have demonstrated that vaccine effectiveness is reduced and the immune response is evaded by variants of concern (VOCs), which include Alpha, Beta, Delta, and, the most recent, Omicron. Subsequently, several vaccine breakthrough (VBT) infections have been reported among healthcare workers (HCWs) due to their prolonged exposure to viruses at healthcare facilities. We conducted a clinico-genomic study of ChAdOx1 (Covishield) VBT cases in HCWs after complete vaccination. Based on the clinical data analysis, most of the cases were categorized as mild, with minimal healthcare support requirements. These patients were divided into two sub-phenotypes based on symptoms: mild and mild plus. Statistical analysis showed a significant correlation of specific clinical parameters with VBT sub-phenotypes. Viral genomic sequence analysis of VBT cases revealed a spectrum of high- and low-frequency mutations. More in-depth analysis revealed the presence of low-frequency mutations within the functionally important regions of SARS-CoV-2 genomes. Emphasizing the potential benefits of surveillance, low-frequency mutations, D144H in the N gene and D138Y in the S gene, were observed to potentially alter the protein secondary structure with possible influence on viral characteristics. Substantiated by the literature, our study highlights the importance of integrative analysis of pathogen genomic and clinical data to offer insights into low-frequency mutations that could be a modulator of VBT infections.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had an enormous burden on the healthcare system worldwide as a consequence of its new emerging variants of concern (VOCs) since late 2019. Elucidating viral genome characteristics and its influence on disease severity and clinical outcome has been one of the crucial aspects toward pandemic management. Genomic surveillance holds the key to identify the spectrum of mutations vis-à-vis disease outcome. Here, in our study, we performed a comprehensive analysis of the mutation distribution among the coronavirus disease 2019 (COVID-19) recovered and mortality patients. In addition to the clinical data analysis, the significant mutations within the two groups were analyzed for their global presence in an effort to understand the temporal dynamics of the mutations globally in comparison with our cohort. Interestingly, we found that all the mutations within the recovered patients showed significantly low global presence, indicating the possibility of regional pool of mutations and the absence of preferential selection by the virus during the course of the pandemic. In addition, we found the mutation S194L to have the most significant occurrence in the mortality group, suggesting its role toward a severe disease progression. Also, we discovered three mutations within the mortality patients with a high cohort and global distribution, which later became a part of variants of interest (VOIs)/VOCs, suggesting its significant role in enhancing viral characteristics. To understand the possible mechanism, we performed molecular dynamics (MD) simulations of nucleocapsid mutations, S194L and S194*, from the mortality and recovered patients, respectively, to examine its impacts on protein structure and stability. Importantly, we observed the mutation S194* within the recovered to be comparatively unstable, hence showing a low global frequency, as we observed. Thus, our study provides integrative insights about the clinical features, mutations significantly associated with the two different clinical outcomes, its global presence, and its possible effects at the structural level to understand the role of mutations in driving the COVID-19 pandemic.

The variability of clinical course and prognosis of COVID-19 highlights the necessity of patient sub-group risk stratification based on clinical data. In this study, clinical data from a cohort of Indian COVID-19 hospitalized patients is used to develop risk stratification and mortality prediction models. We analyzed a set of 70 clinical parameters including physiological and hematological for developing machine learning models to identify biomarkers. We also compared the Indian and Wuhan cohort, and analyzed the role of steroids. A bootstrap averaged ensemble of Bayesian networks was also learned to construct an explainable model for discovering actionable influences on mortality and days to outcome. We discovered blood parameters, diabetes, co-morbidity and SpO2 levels as important risk stratification features, whereas mortality prediction is dependent only on blood parameters. XGboost and logistic regression model yielded the best performance on risk stratification and mortality prediction, respectively (AUC score 0.83, AUC score 0.92). Blood coagulation parameters (ferritin, D-Dimer and INR), immune and inflammation parameters IL6, LDH and Neutrophil (%) are common features for both risk and mortality prediction. Compared with Wuhan patients, Indian patients with extreme blood parameters indicated higher survival rate. Analyses of medications suggest that a higher proportion of survivors and mild patients who were administered steroids had extreme neutrophil and lymphocyte percentages. The ensemble averaged Bayesian network structure revealed serum ferritin to be the most important predictor for mortality and Vitamin D to influence severity independent of days to outcome. The findings are important for effective triage during strains on healthcare infrastructure.

The global coronavirus disease 2019 (COVID-19) pandemic has demonstrated the range of disease severity and pathogen genomic diversity emanating from a singular virus (severe acute respiratory syndrome coronavirus 2, SARS-CoV-2). This diversity in disease manifestations and genomic mutations has challenged healthcare management and resource allocation during the pandemic, especially for countries such as India with a bigger population base. Here, we undertake a combinatorial approach toward scrutinizing the diagnostic and genomic diversity to extract meaningful information from the chaos of COVID-19 in the Indian context. Using methods of statistical correlation, machine learning (ML), and genomic sequencing on a clinically comprehensive patient dataset with corresponding with/without respiratory support samples, we highlight specific significant diagnostic parameters and ML models for assessing the risk of developing severe COVID-19. This information is further contextualized in the backdrop of SARS-CoV-2 genomic features in the cohort for pathogen genomic evolution monitoring. Analysis of the patient demographic features and symptoms revealed that age, breathlessness, and cough were significantly associated with severe disease; at the same time, we found no severe patient reporting absence of physical symptoms. Observing the trends in biochemical/biophysical diagnostic parameters, we noted that the respiratory rate, total leukocyte count (TLC), blood urea levels, and C-reactive protein (CRP) levels were directly correlated with the probability of developing severe disease. Out of five different ML algorithms tested to predict patient severity, the multi-layer perceptron-based model performed the best, with a receiver operating characteristic (ROC) score of 0.96 and an F1 score of 0.791. The SARS-CoV-2 genomic analysis highlighted a set of mutations with global frequency flips and future inculcation into variants of concern (VOCs) and variants of interest (VOIs), which can be further monitored and annotated for functional significance. In summary, our findings highlight the importance of SARS-CoV-2 genomic surveillance and statistical analysis of clinical data to develop a risk assessment ML model.

Disease caused by SARS-CoV-2 coronavirus (COVID-19) led to significant morbidity and mortality worldwide. A systemic hyper-inflammation characterizes severe COVID-19 disease, often associated with acute respiratory distress syndrome (ARDS). Blood biomarkers capable of risk stratification are of great importance in effective triage and critical care of severe COVID-19 patients. Flow cytometry and next-generation sequencing were done on peripheral blood cells and urokinase-type plasminogen activator receptor (suPAR), and cytokines were measured from and mass spectrometry-based proteomics was done on plasma samples from an Indian cohort of COVID-19 patients. Publicly available single-cell RNA sequencing data were analyzed for validation of primary data. Statistical analyses were performed to validate risk stratification. We report here higher plasma abundance of suPAR, expressed by an abnormally expanded myeloid cell population, in severe COVID-19 patients with ARDS. The plasma suPAR level was found to be linked to a characteristic plasma proteome, associated with coagulation disorders and complement activation. Receiver operator characteristic curve analysis to predict mortality identified a cutoff value of suPAR at 1,996.809 pg/ml (odds ratio: 2.9286, 95% confidence interval 1.0427–8.2257). Lower-than-cutoff suPAR levels were associated with a differential expression of the immune transcriptome as well as favorable clinical outcomes, in terms of both survival benefit (hazard ratio: 0.3615, 95% confidence interval 0.1433–0.912) and faster disease remission in our patient cohort. Thus, we identified suPAR as a key pathogenic circulating molecule linking systemic hyperinflammation to the hypercoagulable state and stratifying clinical outcomes in severe COVID-19 patients with ARDS.

Background: Post infection immunity and post vaccination immunity both confer protection against COVID-19. However, there have been many whole genome sequencing proven reinfections and breakthrough infections. Both are most often mild and caused by Variants of Concern (VOC).Methods: The patient in our study underwent serial COVID-19 RT-PCR, blood tests for serology, acute phase reactants, and chest imaging as part of clinical care. We interviewed the patient for clinical history and retrieved reports and case papers. We retrieved stored RT-PCR positive samples for whole genome sequencing (WGS) of SARS-CoV-2 from the patient's breakthrough infections and the presumed index case.Findings: The patient had three RT-PCR confirmed SARS-CoV-2 infections. Two breakthrough infections occurred in quick succession with the first over 3 weeks after complete vaccination with COVISHIELD and despite post-vaccination seroconversion. The first breakthrough infection was due to the Alpha variant and the second due to the Delta variant. The Delta variant infection resulted in hypoxia, hospitalization, and illness lasting seven weeks. Serial serology, acute phase reactants, and chest imaging supported WGS in establishing distinct episodes of infection. WGS established a fully vaccinated family member as the index case.Interpretation: The patient had an Alpha variant breakthrough infection despite past infection, complete vaccination, and seroconversion. Despite boosting after this infection, the patient subsequently had a severe Delta variant breakthrough infection. This was also a WGS proven reinfection and, therefore, a case of breakthrough reinfection. The patient acquired the infection from a fully vaccinated family member.

Human host and pathogen interaction is dynamic in nature and often modulated by co-pathogens with a functional role in delineating the physiological outcome of infection. Co-infection may present either as a pre-existing pathogen which is accentuated by the introduction of a new pathogen or may appear in the form of new infection acquired secondarily due to a compromised immune system. Using diverse examples of co-infecting pathogens such as Human Immunodeficiency Virus, Mycobacterium tuberculosis and Hepatitis C Virus, we have highlighted the role of co-infections in modulating disease severity and clinical outcome. This interaction happens at multiple hierarchies, which are inclusive of stress and immunological responses and together modulate the disease severity. Already published literature provides much evidence in favor of the occurrence of co-infections during SARS-CoV-2 infection, which eventually impacts the Coronavirus disease-19 outcome. The availability of biological models like 3D organoids, mice, cell lines and mathematical models provide us with an opportunity to understand the role and mechanism of specific co-infections. Exploration of multi-omics-based interactions across co-infecting pathogens may provide deeper insights into their role in disease modulation.

Co-infection with ancillary pathogens is a significant modulator of morbidity and mortality in infectious diseases. There have been limited reports of co-infections accompanying SARS-CoV-2 infections, albeit lacking India specific study. The present study has made an effort toward elucidating the prevalence, diversity and characterization of co-infecting respiratory pathogens in the nasopharyngeal tract of SARS-CoV-2 positive patients. Two complementary metagenomics based sequencing approaches, Respiratory Virus Oligo Panel (RVOP) and Holo-seq, were utilized for unbiased detection of co-infecting viruses and bacteria. The limited SARS-CoV-2 clade diversity along with differential clinical phenotype seems to be partially explained by the observed spectrum of co-infections. We found a total of 43 bacteria and 29 viruses amongst the patients, with 18 viruses commonly captured by both the approaches. In addition to SARS-CoV-2, Human Mastadenovirus, known to cause respiratory distress, was present in a majority of the samples. We also found significant differences of bacterial reads based on clinical phenotype. Of all the bacterial species identified, ∼60% have been known to be involved in respiratory distress. Among the co-pathogens present in our sample cohort, anaerobic bacteria accounted for a preponderance of bacterial diversity with possible role in respiratory distress. Clostridium botulinum, Bacillus cereus and Halomonas sp. are anaerobes found abundantly across the samples. Our findings highlight the significance of metagenomics based diagnosis and detection of SARS-CoV-2 and other respiratory co-infections in the current pandemic to enable efficient treatment administration and better clinical management. To our knowledge this is the first study from India with a focus on the role of co-infections in SARS-CoV-2 clinical sub-phenotype.

This study elucidated the clinical, humoral immune response and genomic analysis of vaccine breakthrough (VBT) infections after ChAdOx1 nCoV-19/Covishield vaccine in healthcare workers (HCWs). Amongst 1858 HCWs, 1639 had received either two doses (1346) or a single dose (293) of ChAdOx1 nCoV-19 vaccine. SARS-CoV-2 IgG antibodies and neutralizing antibodies were measured in the vaccinated group and the development of SARS-CoV-2 infection was monitored.Forty-six RT-PCR positive samples from the 203 positive samples were subjected to whole genome sequencing (WGS). Of the 203 (10.92%) infected HCWs, 21.46% (47/219) were non-vaccinated, which was significantly more than 9.52% (156/1639) who were vaccinated and infection was higher in doctors and nurses. Unvaccinated HCWs had 1.57 times higher risk compared to partially vaccinated HCWs and 2.49 times higher risk than those who were fully vaccinated.The partially vaccinated were at higher risk than the fully vaccinated (RR 1.58). Antibody non-response was seen in 3.44% (4/116), low antibody levels in 15.51% (18/116) and medium levels were found in 81.03% (94/116). Fully vaccinated HCWs had a higher antibody response at day 42 than those who were partially vaccinated (8.96 + 4.00 vs. 7.17 + 3.82). Whole genome sequencing of 46 samples revealed that the Delta variant (B.1.617.2) was predominant (69.5%). HCWs who had received two doses of vaccine showed better protection from mild, moderate, or severe infection, with a higher humoral immune response than those who had received a single dose. The genomic analysis revealed the predominance of the Delta variant (B.1.617.2) in the VBT infections.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) manifests a broad spectrum of clinical presentations, varying in severity from asymptomatic to mortality. As the viral infection spread, it evolved and developed into many variants of concern. Understanding the impact of mutations in the SARS-CoV-2 genome on the clinical phenotype and associated co-morbidities is important for treatment and preventionas the pandemic progresses. Based on the mild, moderate, and severe clinical phenotypes, we analyzed the possible association between both, the clinical sub-phenotypes and genomic mutations with respect to the severity and outcome of the patients. We found a significant association between the requirement of respiratory support and co-morbidities. We also identified six SARS-CoV-2 genome mutations that were significantly correlated with severity and mortality in our cohort. We examined structural alterations at the RNA and protein levels as a result of three of these mutations: A26194T, T28854T, and C25611A, present in the Orf3a and N protein. The RNA secondary structure change due to the above mutations can be one of the modulators of the disease outcome. Our findings highlight the importance of integrative analysis in which clinical and genetic components of the disease are co-analyzed. In combination with genomic surveillance, the clinical outcome-associated mutations could help identify individuals for priority medical support.