STRIDE-DB: a comprehensive database for exploration of instability and phenotypic relevance of short tandem repeats in the human genome.

Short Tandem Repeats (STRs) are genetic markers made up of repeating DNA sequences. The variations of the STRs are widely studied in forensic analysis, population studies and genetic testing for a variety of neuromuscular disorders. Understanding polymorphic STR variation and its cause is crucial for deciphering genetic information and finding links to various disorders. In this paper, we present STRIDE-DB, a novel and unique platform to explore STR Instability and its Phenotypic Relevance, and a comprehensive database of STRs in the human genome. We utilized RepeatMasker to identify all the STRs in the human genome (hg19) and combined it with frequency data from the 1000 Genomes Project. STRIDE-DB, a user-friendly resource, plays a pivotal role in investigating the relationship between STR variation, instability and phenotype. By harnessing data from genome-wide association studies (GWAS), ClinVar database, Alu loci, Haploblocks in genome and Conservation of the STRs, it serves as an important tool for researchers exploring the variability of STRs in the human genome and its direct impact on phenotypes. STRIDE-DB has its broad applicability and significance in various research domains like forensic sciences and other repeat expansion disorders. Database URL: https://stridedb.igib.res.in.

Investigation of RFC1 tandem nucleotide repeat locus in diverse neurodegenerative outcomes in an Indian cohort.

An intronic bi-allelic pentanucleotide repeat expansion mutation, (AAGGG)400-2000, at AAAAG repeat locus in RFC1 gene, is known as underlying genetic cause in cases with cerebellar ataxia, neuropathy, and vestibular areflexia syndrome (CANVAS) and late-onset sporadic ataxia. Biallelic positive cases carry a common recessive risk haplotype, "AAGA," spanning RFC1 gene. In this study, our aim is to find prevalence of bi-allelic (AAGGG)exp in Indian ataxia and other neurological disorders and investigate the complexity of RFC1 repeat locus and its potential association with neurodegenerative diseases in Indian population-based cohorts. We carried out repeat number and repeat type estimation using flanking PCR and repeat primed PCR (AAAAG/AAAGG/AAGGG) in four Indian disease cohorts and healthy controls. Haplotype assessment of suspected cases was done by genotyping and confirmed by Sanger sequencing. Blood samples and consent of all the cases and detailed clinical details of positive cases were collected in collaboration with A.I.I.M.S. Furthermore, comprehension of RFC1 repeat locus and risk haplotype analysis in Indian background was performed on the NGS data of Indian healthy controls by ExpansionHunter, ExpansionHunter Denovo, and PHASE analysis, respectively. Genetic screening of RFC1-TNR locus in 1998 uncharacterized cases (SCA12: 87; uncharacterized ataxia: 1818, CMT: 93) and 564 heterogenous controls showed that the frequency of subjects with bi-allelic (AAGGG)exp are 1.15%, < 0.05%, 2.15%, and 0% respectively. Two RFC1 positive sporadic late-onset ataxia cases, one bi-allelic (AAGGG)exp and another, (AAAGG)~700/(AAGGG)exp, had recessive risk haplotype and CANVAS symptoms. Long normal alleles, 15-27, are significantly rare in ataxia cohort. In IndiGen control population (IndiGen; N = 1029), long normal repeat range, 15-27, is significantly associated with A3G3 and some rare repeat motifs, AGAGG, AACGG, AAGAG, and AAGGC. Risk-associated "AAGA" haplotype of the original pathogenic expansion of A2G3 was found associated with the A3G3 representing alleles in background population. Apart from bi-allelic (AAGGG)exp, we report cases with a new pathogenic expansion of (AAAGG)exp/(AAGGG)exp in RFC1 and recessive risk haplotype. We found different repeat motifs at RFC1 TNR locus, like AAAAG, AAAGG, AAAGGG, AAAAGG, AAGAG, AACGG, AAGGC, AGAGG, and AAGGG, in Indian background population except ACAGG and (AAAGG)n/(AAGGG)n. Our findings will help in further understanding the role of long normal repeat size and different repeat motifs, specifically AAAGG, AAAGGG, and other rare repeat motifs, at the RFC1 locus.

Sequencing through hyperexpanded Friedreich's ataxia-GAA repeats by nanopore technology: implications in genotype-phenotype correlation.

Friedreich's ataxia, an autosomal recessive disorder, is caused by tandem GAA nucleotide repeat expansions in intron 1 of the frataxin gene. The GAA repeats over 66 in number are considered as pathogenic, and commonly occurring pathogenic repeats are within a range of 600-1200. Clinically, the spectrum of features is confined mainly to neurological tissues; however, cardiomyopathy and diabetes mellitus have been reported in 60 and 30% of the subjects, respectively. The accurate detection of GAA repeat count is of utmost importance for clinical genetic correlation, and no study so far has attempted an approach that is of high-throughput nature and defines the exact sequence of GAA repeats. Largely, the method for detection of GAA repeats so far is either through the conventional polymerase chain reaction-based screening or Southern blot, which remains the gold standard method. We utilized an approach of long-range targeted amplification of FXN-GAA repeats using Oxford Nanopore Technologies MinION platform for accurate estimation of repeat length. We were able to achieve successful amplification of GAA repeats ranging from ∼120 to 1100 at ∼2600× mean coverage. The total throughput achievable through our protocol can allow for screening of up to 96 samples per flow cell in less than 24 h. The proposed method is clinically scalable and deployable for day-to-day diagnostics. In this paper, we demonstrate to resolve the genotype-phenotype correlation of Friedreich's ataxia patients with better accuracy.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

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.

Genetics of Ataxias in Indian Population: A Collative Insight from a Common Genetic Screening Tool.

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.

Respiratory Co-Infections: Modulators of SARS-CoV-2 Patients' Clinical Sub-Phenotype.

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.

Vaccine Breakthrough Infections by SARS-CoV-2 Variants after ChAdOx1 nCoV-19 Vaccination in Healthcare Workers.

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.

Integrated genomic view of SARS-CoV-2 in India.

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.

EasyQC: Tool with Interactive User Interface for Efficient Next-Generation Sequencing Data Quality Control.

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.

A novel mutation in SGSH causing Sanfillipo type 3A Mucopolysaccharidoses in an Indian family.

Mucopolysaccharidoses (MPS) type III also termed as Sanfillipo syndrome, involves defect in enzymes required for degradation of heparan sulphate. We report a clinical case of MPS-III later followed by genetic investigation for MPS-III genes SGSH, NAGLU, HGSNAT and GNS. It allowed us to identify a novel and likely pathogenic variant p. G205R in SGSH. Protein based Inslico prediction and protein modelling suggests aberration of helical structure of SGSH protein and reduced binding affinity for its substrate.

A novel mutation in SLC39A14 causing hypermanganesemia associated with infantile onset dystonia.

BACKGROUND: Mutations in SLC39A14 cause a recessive disorder of manganese (Mn) metabolism that manifests as childhood onset progressive neurodegeneration characterized by parkinsonism and dystonia. METHODS: The present study genetically investigated a case of hypermanganesemia. We describe a family where an affected child with a history of progressive neurodegeneration showed symptoms of dystonia with increased levels of blood Mn and altered signal intensities in globus pallidus and dentate nucleus. Whole exome sequencing was conducted to genetically investigate the pathology in the child, which allowed us to identify a novel homozygous causal mutation in SLC39A14. RESULTS: Insilico modeling of the novel homozygous causal mutation in SLC39A14 predicted that it was deleterious, affecting Mn binding and transportation of metal by transmembrane instability of the protein structure. The clinical features of other reported mutations in SLC39A14 were also reviewed and the clinical spectrum in our case conforms to the described neurological abnormalities. CONCLUSIONS: We conclude that the mutation identified in SLC39A14 in our case is a novel variation linked to recessive disorders of hypermaganesemia and dystonia.