Background: India first detected SARS-CoV-2, causal agent of COVID-19 in late January 2020, imported from Wuhan, China. From March 2020 onwards, the importation of cases from countries in the rest of the world followed by seeding of local transmission triggered further outbreaks in India. Methods: We used ARTIC protocol-based tiling amplicon sequencing of SARS-CoV-2 (n=104) from different states of India using a combination of MinION and MinIT sequencing from Oxford Nanopore Technology to understand how introduction and local transmission occurred. Results: The analyses revealed multiple introductions of SARS-CoV-2 genomes, including the A2a cluster from Europe and the USA, A3 cluster from Middle East and A4 cluster (haplotype redefined) from Southeast Asia (Indonesia, Thailand and Malaysia) and Central Asia (Kyrgyzstan). The local transmission and persistence of genomes A4, A2a and A3 was also observed in the studied locations. The most prevalent genomes with patterns of variance (confined in a cluster) remain unclassified, and are here proposed as A4-clade based on its divergence within the A cluster. Conclusions: The viral haplotypes may link their persistence to geo-climatic conditions and host response. Multipronged strategies including molecular surveillance based on real-time viral genomic data is of paramount importance for a timely management of the pandemic.

We conclude that the mutation identified in SLC39A14 in our case is a novel variation linked to recessive disorders of hypermaganesemia and dystonia.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

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.

Since its zoonotic transmission in the human host, the SARS-CoV-2 virus has infected millions and has diversified extensively. A hallmark feature of viral system survival is their continuous evolution and adaptation within the host. RNA editing via APOBEC and ADAR family of enzymes has been recently implicated as the major driver of intra-host variability of the SARS-CoV-2 genomes. Analysis of the intra-host single-nucleotide variations (iSNVs) in SARS-CoV-2 genomes at spatio-temporal scales can provide insights on the consequence of RNA editing on the establishment, spread and functional outcomes of the virus. In this study, using 1,347 transcriptomes of COVID-19 infected patients across various populations, we find variable prevalence of iSNVs with distinctly higher levels in Indian population. Our results also suggest that iSNVs can likely establish variants in a population. These iSNVs may also contribute to key structural and functional changes in the Spike protein that confer antibody resistance.

Cerebellar ataxias (CAs) represent a group of autosomal dominant and recessive neurodegenerative disorders affecting cerebellum with or without spinal cord. Overall, CAs have preponderance for tandem nucleotide repeat expansions as an etiological factor (10 TREs explain nearly 30–40% of ataxia cohort globally). The experience of 10 years of common genetic ataxia subtypes for ≈5600 patients’ referrals (Pan‐India) received at a single center is shared herein. Frequencies (in %, n) of SCA types and FRDA in the sample cohort are observed as follows: SCA12 (8.6%, 490); SCA2 (8.5%, 482); SCA1 (4.8%, 272); SCA3 (2%, 113); SCA7 (0.5%, 28); SCA6 (0.1%, 05); SCA17 (0.1%, 05), and FRDA (2.2%, 127). A significant amount of variability in TRE lengths at each locus is observed, we noted presence of biallelic expansion, co‐occurrence of SCA‐subtypes, and the presence of premutable normal alleles. The frequency of mutated GAA‐FRDA allele in healthy controls is 1/158 (0.63%), thus an expected FRDA prevalence of 1:100 000 persons. The data of this study are relevant not only for clinical decision making but also for guidance in direction of genetic investigations, transancestral comparison of genotypes, and lastly provide insight for policy decision for the consideration of SCAs under rare disease category.

The advent of next-generation sequencing (NGS) technologies has revolutionized the world of genomic research. Millions of sequences are generated in a short period of time and they provide intriguing insights to the researcher. Many NGS platforms have evolved over a period of time and their efficiency has been ever increasing. Still, primarily because of the chemistry, glitch in the sequencing machine and human handling errors, some artifacts tend to exist in the final sequence data set. These sequence errors have a profound impact on the downstream analyses and may provide misleading information. Hence, filtering of these erroneous reads has become inevitable and myriad of tools are available for this purpose. However, many of them are accessible as a command line interface that requires the user to enter each command manually. Here, we report EasyQC, a tool for NGS data quality control (QC) with a graphical user interface providing options to carry out trimming of NGS reads based on quality, length, homopolymer, and ambiguous bases. EasyQC also possesses features such as format converter, paired end merger, adapter trimmer, and a graph generator that generates quality distribution, length distribution, GC content, and base composition graphs. Comparison of raw and processed sequence data sets using EasyQC suggested significant increase in overall quality of the sequences. Testing of EasyQC using NGS data sets on a standalone desktop proved to be relatively faster. EasyQC is developed using PERL modules and can be executed in Windows and Linux platforms. With the various QC features, easy interface for end users, and cross-platform compatibility, EasyQC would be a valuable addition to the already existing tools facilitating better downstream analyses.

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.