iCEED: Integrated customized extraction of enzyme data.

Enzymes catalyze diverse biochemical reactions and are building blocks of cellular and metabolic pathways. Data and metadata of enzymes are distributed across databases and are archived in various formats. The enzyme databases provide utilities for efficient searches and downloading enzyme records in batch mode but do not support organism-specific extraction of subsets of data. Users are required to write scripts for parsing entries for customized data extraction prior to downstream analysis. Integrated Customized Extraction of Enzyme Data (iCEED) has been developed to provide organism-specific customized data extraction utilities for seven commonly used enzyme databases and brings these resources under an integrated portal. iCEED provides dropdown menus and search boxes using typehead utility for submission of queries as well as enzyme class-based browsing utility. A utility to facilitate mapping and visualization of functionally important features on the three-dimensional (3D) structures of enzymes is integrated. The customized data extraction utilities provided in iCEED are expected to be useful for biochemists, biotechnologists, computational biologists, and life science researchers to build curated datasets of their choice through an easy to navigate web-based interface. The integrated feature visualization system is useful for a fine-grained understanding of the enzyme structure-function relationship. Desired subsets of data, extracted and curated using iCEED can be subsequently used for downstream processing, analyses, and knowledge discovery. iCEED can also be used for training and teaching purposes.

IDV Typer: An Automated Tool for Lineage Typing of Influenza D Viruses Based on Return Time Distribution.

Influenza D virus (IDV) is the most recent addition to the Orthomyxoviridae family and cattle serve as the primary reservoir. IDV has been implicated in Bovine Respiratory Disease Complex (BRDC), and there is serological evidence of human infection of IDV. Evolutionary changes in the IDV genome have resulted in the expansion of genetic diversity and the emergence of multiple lineages that might expand the host tropism and potentially increase the pathogenicity to animals and humans. Therefore, there is an urgent need for automated, accurate and rapid typing tools for IDV lineage typing. Currently, IDV lineage typing is carried out using BLAST-based searches and alignment-based molecular phylogeny of the hemagglutinin-esterase fusion (HEF) gene sequences, and lineage is assigned to query sequences based on sequence similarity (BLAST search) and proximity to the reference lineages in the tree topology, respectively. To minimize human intervention and lineage typing time, we developed IDV Typer server, implementing alignment-free method based on return time distribution (RTD) of k-mers. Lineages are assigned using HEF gene sequences. The server performs with 100% sensitivity and specificity. The IDV Typer server is the first application of an RTD-based alignment-free method for typing animal viruses.

Neutralizing Antibody Response to Genotypically Diverse Measles Viruses in Clinically Suspected Measles Cases.

The neutralizing antibody (Nt-Ab) response to vaccine and wild-type measles viruses (MeV) was studied in suspected measles cases reported during the years 2012-2016. The neutralization activity against MeV A, D4 and D8 genotypes was studied on sera (Panel A; n = 68 (measles-immunized) and Panel B; n = 50 (unvaccinated)) that were either laboratory confirmed or not confirmed by the presence of IgM antibodies. Additionally, the Nt-Ab response in Panel A was measured against the MeV vaccine and four wild-type viruses. Neutralization results were compared using homology modeling and molecular dynamics simulation (MDS) of MeV-hemagglutinin (H) and fusion (F) proteins. Overall, the Nt-Ab titres for MeV-A were found to be significantly lower than MeV-D4 and MeV-D8 viruses for Panel A. No major difference was noted in Nt-Ab titres between MeV-D8 viruses (Jamnagar and New Delhi), whereas MeV-D4 (Sindhudurg and Bagalkot (BGK) viruses) showed significant differences between Nt-Ab titres for Panel B. Interestingly, the substitutions observed in epitopes of H-protein, L249P and G316A are observed to be unique to MeV-BGK. MDS of H-protein revealed significant fluctuations in neutralizing epitopes due to L249P substitution. The majority of the clinically suspected cases showed Nt-Abs to MeV wild-types. Higher IgG antibody avidity and Nt-Ab titres were noted in IgM-negatives than in IgM-positives cases, indicating reinfection or breakthrough. MDS revealed reduced neutralization due to decreased conformational flexibility in the H-epitope.

Women in the European Virus Bioinformatics Center.

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

Genetic variability in minor capsid protein (L2 gene) of human papillomavirus type 16 among Indian women.

Human papillomavirus type 16 (HPV-16) is the predominant genotype worldwide associated with invasive cervical cancer and hence remains as the focus for diagnostic development and vaccine research. L2, the minor capsid protein forms the packaging unit for the HPV genome along with the L1 protein and is primarily associated with transport of genomic DNA to the nucleus. Unlike L1, L2 is known to elicit cross-neutralizing antibodies and thus becomes a suitable candidate for pan-HPV prophylactic vaccine development. In the present study, a total of 148 cervical HPV-16 isolates from Indian women were analyzed by PCR-directed sequencing, phylogenetic analysis and in silico immunoinformatics tools to determine the L2 variations that may impact the immune response and oncogenesis. Ninety-one SNPs translating to 35 non-synonymous amino acid substitutions were observed, of these 16 substitutions are reported in the Indian isolates for the first time. T245A, L266F, S378V and S384A substitutions were significantly associated with high-grade cervical neoplastic status. Multiple substitutions were observed in samples from high-grade cervical neoplastic status as compared to those from normal cervical status (p = 0.027), specifically from the D3 sub-lineage. It was observed that substitution T85A was part of both, B and T cell epitopes recognized by MHC-I molecules; T245A was common to B and T cell epitopes recognized by MHC-II molecules and S122P/A was common to the region recognized by both MHC-I and MHC-II molecules. These findings reporting L2 protein substitutions have implications on cervical oncogenesis and design of next-generation L2-based HPV vaccines.

Phylogenetic Analysis of Spread of Hepatitis C Virus Identified during HIV Outbreak Investigation, Unnao, India.

An HIV outbreak investigation during 2017-2018 in Unnao District, Uttar Pradesh, India, unearthed high prevalence of hepatitis C virus (HCV) antibodies among the study participants. We investigated these HCV infections by analyzing NS5B and core regions. We observed no correlation between HIV-HCV viral loads and clustering of HCV sequences, regardless of HIV serostatus. All HCV isolates belonged to genotype 3a. Monophyletic clustering of isolates in NS5B phylogeny indicates emergence of the outbreak from a single isolate or its closely related descendants. The nucleotide substitution rate for NS5B was 6 × 10-3 and for core was 2 × 10-3 substitutions/site/year. Estimated time to most recent common ancestor of these isolates was 2012, aligning with the timeline of this outbreak, which might be attributable to unsafe injection practices while seeking healthcare. HIV-HCV co-infection underlines the need for integrated testing, surveillance, strengthening of healthcare systems, community empowerment, and molecular analyses as pragmatic public health tools.

Genetic variations in the long control region of human papillomavirus type 16 isolates from India: implications for cervical carcinogenesis.

Introduction. Infection with high-risk human papillomavirus (HPV) types, specifically HPV type 16 (HPV16), is considered to be the most important risk factor in the development of cervical intraepithelial neoplasia and cancer. The long control region (LCR) is a noncoding region that comprises approximately 10 % of the HPV genome and contains regulatory elements for viral transcription and replication. Sequence variations in LCR may impact on the replication efficiency and oncogenic potential of the virus.Gap statement. Studies documenting variations in LCR of HPV16 isolates pertaining to cervical neoplastic status in India are limited.Aim. The present study was designed to characterize variations in the LCR of Indian isolates of HPV16 and study their association with cervical disease grades.Methodology. The LCR was amplified and sequenced from HPV16 positive cervical samples belonging to different cervical disease grades. Sequences were aligned to identify variations and potential transcription factor binding sites (TFbs) were predicted using the JASPAR database in addition to phylogenetic studies.Results. Among the 163 HPV16 isolates analysed, 47 different nucleotide variations were detected in the LCR, of which 25 are reported for first time in Indian isolates. Point mutations were detected in 35/54 (64.8 %) samples with normal cervical status, 44/50 (88 %) samples with low-grade cervical disease and 53/59 (89.8 %) samples with high-grade cervical disease. Variations T6586C, G6657A and T6850G were significantly associated with high-grade cervical status. Thirteen LCR variations were detected in the binding sites for CEBPB, ETS1, JUN, MYB, NFIL3, PHOX2A and SOX9 transcription factors.Conclusion. The present study helped to identify unique variations in the LCRs of HPV16 Indian isolates. The variations in the A4 sub-lineage were significantly associated with high-grade disease status. The isolates belonging to the A4 and D3 sub-lineages harboured mutations in putative TFbs, implying a potential impact on viral replication and progression to cervical cancer.

Circulation and Evolution of SARS-CoV-2 in India: Let the Data Speak.

The COVID-19 pandemic is a global challenge that impacted 200+ countries. India ranks in the second and third positions in terms of number of reported cases and deaths. Being a populous country with densely packed cities, SARS-CoV-2 spread exponentially. India sequenced ≈0.14% isolates from confirmed cases for pandemic surveillance and contributed ≈1.58% of complete genomes sequenced globally. This study was designed to map the circulating lineage diversity and to understand the evolution of SARS-CoV-2 in India using comparative genomics and population genetics approaches. Despite varied sequencing coverage across Indian States and Union Territories, isolates belonging to variants of concern (VoC) and variants of interest (VoI) circulated, persisted, and diversified during the first seventeen months of the pandemic. Delta and Kappa lineages emerged in India and spread globally. The phylogenetic tree shows lineage-wise monophyletic clusters of VoCs/VoIs and diversified tree topologies for non-VoC/VoI lineages designated as 'Others' in this study. Evolutionary dynamics analyses substantiate a lack of spatio-temporal clustering, which is indicative of multiple global and local introductions. Sites under positive selection and significant variations in spike protein corroborate with the constellation of mutations to be monitored for VoC/VoI as well as substitutions that are characteristic of functions with implications in virus-host interactions, differential glycosylation, immune evasion, and escape from neutralization.

Inhibition of miR-155 Promotes TGF-ß Mediated Suppression of HIV Release in the Cervical Epithelial Cells.

TGF-ß has been shown to play a differential role in either restricting or aiding HIV infection in different cell types, however its role in the cervical cells is hitherto undefined. Among females, more than 80% of infections occur through heterosexual contact where cervicovaginal mucosa plays a critical role, however the early events during the establishment of infection at female genital mucosa are poorly understood. We earlier showed that increased TGF-ß level has been associated with cervical viral shedding in the HIV infected women, however a causal relationship could not be examined. Therefore, here we first established an in vitro cell-associated model of HIV infection in the cervical epithelial cells (ME-180) and demonstrated that TGF-ß plays an important role as a negative regulator of HIV release in the infected cervical epithelial cells. Inhibition of miR-155 upregulated TGF-ß signaling and mRNA expression of host restriction factors such as APOBEC-3G, IFI-16 and IFITM-3, while decreased the HIV release in ME-180 cells. To conclude, this is the first study to decipher the complex interplay between TGF-ß, miR-155 and HIV release in the cervical epithelial cells. Collectively, our data suggest the plausible role of TGF-ß in promoting HIV latency in cervical epithelial cells which needs further investigations.

One year into the pandemic: Short-term evolution of SARS-CoV-2 and emergence of new lineages.

The COVID-19 pandemic was officially declared on March 11th, 2020. Since the very beginning, the spread of the virus has been tracked nearly in real-time by worldwide genome sequencing efforts. As of March 2021, more than 830,000 SARS-CoV-2 genomes have been uploaded in GISAID and this wealth of data allowed researchers to study the evolution of SARS-CoV-2 during this first pandemic year. In parallel, nomenclatures systems, often with poor consistency among each other, have been developed to designate emerging viral lineages. Despite general fears that the virus might mutate to become more virulent or transmissible, SARS-CoV-2 genetic diversity has remained relatively low during the first ~ 8 months of sustained human-to-human transmission. At the end of 2020/beginning of 2021, though, some alarming events started to raise concerns of possible changes in the evolutionary trajectory of the virus. Specifically, three new viral variants associated with extensive transmission have been described as variants of concern (VOC). These variants were first reported in the UK (B.1.1.7), South Africa (B.1.351) and Brazil (P.1). Their designation as VOCs was determined by an increase of local cases and by the high number of amino acid substitutions harboured by these lineages. This latter feature is reminiscent of viral sequences isolated from immunocompromised patients with long-term infection, suggesting a possible causal link. Here we review the events that led to the identification of these lineages, as well as emerging data concerning their possible implications for viral phenotypes, reinfection risk, vaccine efficiency and epidemic potential. Most of the available evidence is, to date, provisional, but still represents a starting point to uncover the potential threat posed by the VOCs. We also stress that genomic surveillance must be strengthened, especially in the wake of the vaccination campaigns.

Genetic and antigenic characterization of wild type rubella viruses isolated from India.

Rubella, is a contagious disease caused by Rubella virus (RuV) that manifests as fever with skin-rashes in children and adults along with complications in pregnant women. WHO-SEAR has set a target for Rubella elimination by 2023. This is the first report of antigenic characterization and genome sequencing of nine RuVs sampled during 1992, 2007-9, and 2015-17 from four Indian states. Comparative analysis of Indian RuVs (2B) with that of global isolates and vaccine strain RA 27/3 (1a) revealed that the observed mutations in structural proteins have no major impact on the 3D structure, function and antigenicity. Indian RuVs formed three major clusters (Pune-1992, Kannur-2009 and Chitradurg-2007) in genome-based phylogeny of global isolates. Neutralizing antibody titers in a panel of serum samples from measles negative cases were significantly higher to the vaccine strain compared to a wild-type 2B isolate (Kannur) with concordance of 91.9%, thereby substantiating the use of current vaccines.

Antibody Class(es) Predictor for Epitopes (AbCPE): A Multi-Label Classification Algorithm.

Development of vaccines and therapeutic antibodies to deal with infectious and other diseases are the most perceptible scientific interventions that have had huge impact on public health including that in the current Covid-19 pandemic. From inactivation methodologies to reverse vaccinology, vaccine development strategies of 21st century have undergone several transformations and are moving towards rational design approaches. These developments are driven by data as the combinatorials involved in antigenic diversity of pathogens and immune repertoire of hosts are enormous. The computational prediction of epitopes is central to these developments and numerous B-cell epitope prediction methods developed over the years in the field of immunoinformatics have contributed enormously. Most of these methods predict epitopes that could potentially bind to an antibody regardless of its type and only a few account for antibody class specific epitope prediction. Recent studies have provided evidence of more than one class of antibodies being associated with a particular disease. Therefore, it is desirable to predict and prioritize 'peptidome' representing B-cell epitopes that can potentially bind to multiple classes of antibodies, as an open problem in immunoinformatics. To address this, AbCPE, a novel algorithm based on multi-label classification approach has been developed for prediction of antibody class(es) to which an epitope can potentially bind. The epitopes binding to one or more antibody classes (IgG, IgE, IgA and IgM) have been used as a knowledgebase to derive features for prediction. Multi-label algorithms, Binary Relevance and Label Powerset were applied along with Random Forest and AdaBoost. Classifier performance was assessed using evaluation measures like Hamming Loss, Precision, Recall and F1 score. The Binary Relevance model based on dipeptide composition, Random Forest and AdaBoost achieved the best results with Hamming Loss of 0.1121 and 0.1074 on training and test sets respectively. The results obtained by AbCPE are promising. To the best of our knowledge, this is the first multi-label method developed for prediction of antibody class(es) for sequential B-cell epitopes and is expected to bring a paradigm shift in the field of immunoinformatics and immunotherapeutic developments in synthetic biology. The AbCPE web server is available at http://bioinfo.unipune.ac.in/AbCPE/Home.html.

Understanding evolution of SARS-CoV-2: A perspective from analysis of genetic diversity of RdRp gene.

Coronavirus disease 2019 emerged as the first example of "Disease X", a hypothetical disease of humans caused by an unknown infectious agent that was named as novel coronavirus and subsequently designated as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The origin of the outbreak at the animal market in Wuhan, China implies it as a case of zoonotic spillover. The study was designed to understand evolution of Betacoronaviruses and in particular diversification of SARS-CoV-2 using RNA dependent RNA polymerase (RdRp) gene, a stable genetic marker. Phylogenetic and population stratification analyses were carried out using maximum likelihood and Bayesian methods, respectively. Molecular phylogeny using RdRp showed that SARS-CoV-2 isolates cluster together. Bat-CoV isolate RaTG13 and Pangolin-CoVs are observed to branch off prior to SARS-CoV-2 cluster. While SARS-CoV form a single cluster, Bat-CoVs form multiple clusters. Population-based analyses revealed that both SARS-CoV-2 and SARS-CoV form separate clusters with no admixture. Bat-CoVs were found to have single and mixed ancestry and clustered as four sub-populations. Population-based analyses of Betacoronaviruses using RdRp revealed that SARS-CoV-2 is a homogeneous population. SARS-CoV-2 appears to have evolved from Bat-CoV isolate RaTG13, which diversified from a common ancestor from which Pangolin-CoVs have also evolved. The admixed Bat-CoV sub-populations indicate that bats serve as reservoirs harboring virus ensembles that are responsible for zoonotic spillovers such as SARS-CoV and SARS-CoV-2. The extent of admixed isolates of Bat-CoVs observed in population diversification studies underline the need for periodic surveillance of bats and other animal reservoirs for potential spillovers as a measure towards preparedness for emergence of zoonosis.

Characterization of diversity of measles viruses in India: Genomic sequencing and comparative genomics studies.

OBJECTIVE: To map genomic diversity of Measles virus (MeV) isolates collected during 2009-2017 from ten states of India. METHODS: Genome sequencing of Indian isolates and comparative genomics with global MeV using phylogeny, population stratification and selection pressure approaches were performed. RESULTS: The first report of complete genome sequences of forty-three Indian MeV isolates belonging to genotypes D4 (eight) and D8 (thirty-five). Three Indian isolates mapped to named strains D4-Enfield, D8-Villupuram and D8-Victoria. Indian D4 isolates deviate from standard genome length due to indels in M-F intergenic region. Estimated nucleotide substitution rates of Indian MeV derived using genome and individual genes are lower than that of global isolates. Phylogeny revealed genotype-based temporal clustering, suggesting existence of two lineages of D4 and three lineages of D8 in India. Absence of spatial clustering suggests role of cross-border travel in MeV transmission. CONCLUSIONS: Evolutionary analyses suggest the need for surveillance of MeV in India, particularly in view of diversified trajectories of D4 and D8 isolates. This study contributes to global measles epidemiology and indicates no major impact on antigenicity in Indian isolates, thereby substantiating the use of current vaccines to meet measles elimination target of 2023 set by World Health Organization for South-East Asia Region.

Characterization of major capsid protein (L1) variants of Human papillomavirus type 16 by cervical neoplastic status in Indian women: Phylogenetic and functional analysis.

The etiological role of infection with Human papillomavirus type 16 (HPV16) in cervical cancer is well established. HPV16 variants, classified based on less than 10% nucleotide variations in the major capsid (L1 ORF) are known to contribute to persistent infection leading to cancer development. L1 protein forms the cornerstone of HPV structure and antigenicity. In the present study, HPV16 L1 variants were characterized by cervical lesion grade and variations in sequences were correlated to structure and function. The L1 gene was analyzed in 152 HPV16 positive cervical samples obtained from Indian women using polymerase chain reaction-directed sequencing. Phylogenetic analysis was carried out for lineage typing. Sixty-one SNPs were detected in L1 genes resulting in 20 nonsynonymous amino acid substitutions of which N56T, N92T, L158F, V178G, N181I, K236T, K443Q, K454T, and K475R are reported in Indian isolates for the first time. The substitutions N181T, T353P, and T389S were significantly associated with high-grade cervical disease. The predominance of lineage A (A1-A4, 84.96%) was observed among the isolates, while the D3 sublineage showed significant association with high-grade cervical lesions. No evidence for recombination and the positive selection was obtained. These substitutions, when mapped on three-dimensional structure, revealed that 11 and 4 substitutions are part of experimentally validated B- and T-cell epitopes, of which T266A and N285T were common to both types of epitopes and may impact HPV vaccine efficacy. The variants identified through this study have the potential to serve as translational leads for designing diagnostic probes and vaccines.

Erratum for Byukusenge et al., "Complete Genome Sequences of Four Bovine Coronavirus Isolates from Pennsylvania".

[This corrects the article DOI: 10.1128/genomeA.00467-18.].

Complete Genome Sequences of Four Bovine Coronavirus Isolates from Pennsylvania.

We report four full-genome sequences of bovine coronavirus (BCoV) isolates from dairy calves in Pennsylvania obtained in 2016 and 2017. BCoV is a pathogen of great importance to cattle health, and this is the first report of full-genome sequences of BCoV from PA cattle.

AllerBase: a comprehensive allergen knowledgebase.

Database URL: http://bioinfo.net.in/AllerBase/Home.html.

Characterization of bacterial community associated with phytoplankton bloom in a eutrophic lake in South Norway using 16S rRNA gene amplicon sequence analysis.

Interactions between different phytoplankton taxa and heterotrophic bacterial communities within aquatic environments can differentially support growth of various heterotrophic bacterial species. In this study, phytoplankton diversity was studied using traditional microscopic techniques and the bacterial communities associated with phytoplankton bloom were studied using High Throughput Sequencing (HTS) analysis of 16S rRNA gene amplicons from the V1-V3 and V3-V4 hypervariable regions. Samples were collected from Lake Akersvannet, a eutrophic lake in South Norway, during the growth season from June to August 2013. Microscopic examination revealed that the phytoplankton community was mostly represented by Cyanobacteria and the dinoflagellate Ceratium hirundinella. The HTS results revealed that Proteobacteria (Alpha, Beta, and Gamma), Bacteriodetes, Cyanobacteria, Actinobacteria and Verrucomicrobia dominated the bacterial community, with varying relative abundances throughout the sampling season. Species level identification of Cyanobacteria showed a mixed population of Aphanizomenon flos-aquae, Microcystis aeruginosa and Woronichinia naegeliana. A significant proportion of the microbial community was composed of unclassified taxa which might represent locally adapted freshwater bacterial groups. Comparison of cyanobacterial species composition from HTS and microscopy revealed quantitative discrepancies, indicating a need for cross validation of results. To our knowledge, this is the first study that uses HTS methods for studying the bacterial community associated with phytoplankton blooms in a Norwegian lake. The study demonstrates the value of considering results from multiple methods when studying bacterial communities.

Genetic diversity and evolution of dengue virus serotype 3: A comparative genomics study.

Dengue virus serotype 3 (DENV-3), one of the four serotypes of Dengue viruses, is geographically diverse. There are five distinct genotypes (I-V) of DENV-3. Emerging strains and lineages of DENV-3 are increasingly being reported. Availability of genomic data for DENV-3 strains provides opportunity to study its population structure. Complete genome sequences are available for 860 strains of four genotypes (I, II, III and V) isolated worldwide and were analyzed using population genetics and evolutionary approaches to map landscape of genomic diversity. DENV-3 population is observed to be stratified into five major subpopulations. Genotype I and II formed independent subpopulations while genotype III is subdivided into three subpopulations (GIII-a, GIII-b and GIII-c) and is therefore heterogeneous. Genotypes I, II and GIII-a subpopulations comprise of Asian strains whereas GIII-c comprises of American strains. GIII-b subpopulation includes mainly of American strains along with a few strains from Sri Lanka. Genetic admixture is predominantly observed in Sri Lankan strains of genotype III and all strains of genotype V. Inter-genotype recombination was observed to occur in non-structural region of several Asian strains whereas extent of recombination was limited in American strains. Significant positive selection was found to be operational on all genes and observed to be the main driving force of genetic diversity. Positive selection was strongly operational on the branches leading to Asian genotypes and helped to delineate the genetic differences between Asian and American lineages. Thus, inter-genotype recombination, migration and adaptive evolution are the major determinants of evolution of DENV-3.