Twenty-five years of Genomes OnLine Database (GOLD): data updates and new features in v.9.

The Genomes OnLine Database (GOLD) (https://gold.jgi.doe.gov/) at the Department of Energy Joint Genome Institute (DOE-JGI) continues to maintain its role as one of the flagship genomic metadata repositories of the world. The ever-increasing number of projects and metadata are freely available to the user community world-wide. GOLD's metadata is consumed by scientists and remains an important source for large-scale comparative genomics analysis initiatives. Encouraged by this active user engagement and growth, GOLD has continued to add new components and capabilities. The new features such as a public Application Programming Interface (API) and Ecosystem landing page as well as the growth of different entities in this current GOLD v.9 edition are described in detail in this manuscript.

An expectation-maximization framework for comprehensive prediction of isoform-specific functions.

MOTIVATION: Advances in RNA sequencing technologies have achieved an unprecedented accuracy in the quantification of mRNA isoforms, but our knowledge of isoform-specific functions has lagged behind. There is a need to understand the functional consequences of differential splicing, which could be supported by the generation of accurate and comprehensive isoform-specific gene ontology annotations. RESULTS: We present isoform interpretation, a method that uses expectation-maximization to infer isoform-specific functions based on the relationship between sequence and functional isoform similarity. We predicted isoform-specific functional annotations for 85 617 isoforms of 17 900 protein-coding human genes spanning a range of 17 430 distinct gene ontology terms. Comparison with a gold-standard corpus of manually annotated human isoform functions showed that isoform interpretation significantly outperforms state-of-the-art competing methods. We provide experimental evidence that functionally related isoforms predicted by isoform interpretation show a higher degree of domain sharing and expression correlation than functionally related genes. We also show that isoform sequence similarity correlates better with inferred isoform function than with gene-level function. AVAILABILITY AND IMPLEMENTATION: Source code, documentation, and resource files are freely available under a GNU3 license at https://github.com/TheJacksonLaboratory/isopretEM and https://zenodo.org/record/7594321.

Genomes OnLine Database (GOLD) v.8: overview and updates.

The Genomes OnLine Database (GOLD) (https://gold.jgi.doe.gov/) is a manually curated, daily updated collection of genome projects and their metadata accumulated from around the world. The current version of the database includes over 1.17 million entries organized broadly into Studies (45 770), Organisms (387 382) or Biosamples (101 207), Sequencing Projects (355 364) and Analysis Projects (283 481). These four levels contain over 600 metadata fields, which includes 76 controlled vocabulary (CV) tables containing 3873 terms. GOLD provides an interactive web user interface for browsing and searching by a wide range of project and metadata fields. Users can enter details about their own projects in GOLD, which acts as a gatekeeper to ensure that metadata is accurately documented before submitting sequence information to the Integrated Microbial Genomes (IMG) system for analysis. In order to maintain a reference dataset for use by members of the scientific community, GOLD also imports projects from public repositories such as GenBank and SRA. The current status of the database, along with recent updates and improvements are described in this manuscript.

Prenatal phenotyping: A community effort to enhance the Human Phenotype Ontology.

Technological advances in both genome sequencing and prenatal imaging are increasing our ability to accurately recognize and diagnose Mendelian conditions prenatally. Phenotype-driven early genetic diagnosis of fetal genetic disease can help to strategize treatment options and clinical preventive measures during the perinatal period, to plan in utero therapies, and to inform parental decision-making. Fetal phenotypes of genetic diseases are often unique and at present are not well understood; more comprehensive knowledge about prenatal phenotypes and computational resources have an enormous potential to improve diagnostics and translational research. The Human Phenotype Ontology (HPO) has been widely used to support diagnostics and translational research in human genetics. To better support prenatal usage, the HPO consortium conducted a series of workshops with a group of domain experts in a variety of medical specialties, diagnostic techniques, as well as diseases and phenotypes related to prenatal medicine, including perinatal pathology, musculoskeletal anomalies, neurology, medical genetics, hydrops fetalis, craniofacial malformations, cardiology, neonatal-perinatal medicine, fetal medicine, placental pathology, prenatal imaging, and bioinformatics. We expanded the representation of prenatal phenotypes in HPO by adding 95 new phenotype terms under the Abnormality of prenatal development or birth (HP:0001197) grouping term, and revised definitions, synonyms, and disease annotations for most of the 152 terms that existed before the beginning of this effort. The expansion of prenatal phenotypes in HPO will support phenotype-driven prenatal exome and genome sequencing for precision genetic diagnostics of rare diseases to support prenatal care.

Antibacterial activity of flavonoid isolated from Trianthema decandra against Pseudomonas aeruginosa and molecular docking study of FabZ.

The natural product flavonoid demonstrates an extensive sort of pharmacological properties including antimicrobial activity. Although its Pseudomonas aeruginosa inhibition has been discovered, no target for action against flavonoid has been revealed to date. The anti - P. aeruginosa activity of the 2 - (3', 4' dihydroxy-phenyl) - 3, 5, 7-trihydroxy-chromen-4-one isolated from T. decandra was evaluated by disc diffusion and minimum inhibitory concentration methods. The molecular docking of the flavonoid isolated from T. decandra was carried out using CDOCKER (Discovery Studio 2.0). The flavonoid isolated from T. decandra was found to inhibit the growth of P. aeruginosa and the zone of inhibition was found to be 22 ±â€¯0.04 mm at 20 µg/ml while chloramphenicol showed 23 ±â€¯0.05 mm at 30 µg/ml. P. aeruginosa was found to be the most sensitive to both isolated flavonoid and standard control chloramphenicol with MIC values 39.05 µg/ml and 25 µg/ml respectively. Further, the FAS II ß-hydroxyacyl-ACP (FabZ) of P. aeruginosa was found to be a potential target of the flavonoid as it docked in silico effectively. Our work has demonstrated the anti - P. aeruginosa activity of flavonoid isolated from T. decandra and also resulted in the elucidation of a plausible mechanism of action of the isolated flavonoid by inhibiting the FabZ using in silico analysis.

Mycobacterium tuberculosis Whole Genome Sequences From Southern India Suggest Novel Resistance Mechanisms and the Need for Region-Specific Diagnostics.

BACKGROUND.: India is home to 25% of all tuberculosis cases and the second highest number of multidrug resistant cases worldwide. However, little is known about the genetic diversity and resistance determinants of Indian Mycobacterium tuberculosis, particularly for the primary lineages found in India, lineages 1 and 3. METHODS.: We whole genome sequenced 223 randomly selected M. tuberculosis strains from 196 patients within the Tiruvallur and Madurai districts of Tamil Nadu in Southern India. Using comparative genomics, we examined genetic diversity, transmission patterns, and evolution of resistance. RESULTS.: Genomic analyses revealed (11) prevalence of strains from lineages 1 and 3, (11) recent transmission of strains among patients from the same treatment centers, (11) emergence of drug resistance within patients over time, (11) resistance gained in an order typical of strains from different lineages and geographies, (11) underperformance of known resistance-conferring mutations to explain phenotypic resistance in Indian strains relative to studies focused on other geographies, and (11) the possibility that resistance arose through mutations not previously implicated in resistance, or through infections with multiple strains that confound genotype-based prediction of resistance. CONCLUSIONS.: In addition to substantially expanding the genomic perspectives of lineages 1 and 3, sequencing and analysis of M. tuberculosis whole genomes from Southern India highlight challenges of infection control and rapid diagnosis of resistant tuberculosis using current technologies. Further studies are needed to fully explore the complement of diversity and resistance determinants within endemic M. tuberculosis populations.

A corpus of GA4GH Phenopackets: case-level phenotyping for genomic diagnostics and discovery.

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present phenopacket-store. Version 0.1.12 of phenopacket-store includes 4916 phenopackets representing 277 Mendelian and chromosomal diseases associated with 236 genes, and 2872 unique pathogenic alleles curated from 605 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

De novo TRPM3 missense variant associated with neurodevelopmental delay and manifestations of cerebral palsy.

We identified a de novo heterozygous transient receptor potential cation channel subfamily M (melastatin) member 3 (TRPM3) missense variant, p.(Asn1126Asp), in a patient with developmental delay and manifestations of cerebral palsy (CP) using phenotype-driven prioritization analysis of whole-genome sequencing data with Exomiser. The variant is localized in the functionally important ion transport domain of the TRPM3 protein and predicted to impact the protein structure. Our report adds TRPM3 to the list of Mendelian disease-associated genes that can be associated with CP and provides further evidence for the pathogenicity of the variant p.(Asn1126Asp).

Resistance-associated epitopes of HIV-1C-highly probable candidates for a multi-epitope vaccine.

Earlier studies have identified a large number of immunogenic epitopes in HIV-1. Efforts are required to prioritize these epitopes in order to identify the best candidates for formulating an effective multi-epitope vaccine for HIV. We modeled 155 known cytotoxic T lymphocyte epitopes of HIV-1 subtype C on the 3D structure of HLA-A*0201, HLA-B*2705, and HLA-B*5101 using MODPROPEP, as these alleles are known to be associated with resistance to HIV/slow progression to AIDS. Thirty-six epitopes were identified to bind to all the three HLA alleles with better binding affinity than the control peptides complexed with each HLA allele but not to any of the HLA alleles reported to be associated with susceptibility to HIV infection/rapid progression to disease. As increase in stability of the epitope-HLA complex results in increased immunogenicity, the short-listed epitopes could be suitable candidates for vaccine development. Twenty of the 36 epitopes were polyfunctional in nature adding to their immunological relevance for vaccine design. Further, 9 of the 20 polyfunctional epitopes were found to bind to all three resistance-associated HLA alleles using an additional method, adding worth to their potential as candidates for a vaccine formulation for HIV-1C.

Microbiome Metadata Standards: Report of the National Microbiome Data Collaborative's Workshop and Follow-On Activities.

Microbiome samples are inherently defined by the environment in which they are found. Therefore, data that provide context and enable interpretation of measurements produced from biological samples, often referred to as metadata, are critical. Important contributions have been made in the development of community-driven metadata standards; however, these standards have not been uniformly embraced by the microbiome research community. To understand how these standards are being adopted, or the barriers to adoption, across research domains, institutions, and funding agencies, the National Microbiome Data Collaborative (NMDC) hosted a workshop in October 2019. This report provides a summary of discussions that took place throughout the workshop, as well as outcomes of the working groups initiated at the workshop.

Immunological and proteomic analysis of preparative isoelectric focusing separated culture filtrate antigens of Mycobacterium tuberculosis.

Isolation of the secreted proteins and studying the immune response they induce is an essential prerequisite for understanding the pathogenesis of M. tuberculosis. In this study, preparative liquid-phase isoelectric focusing was used for the separation of culture filtrate protein (CFP) of M. tuberculosis. This procedure resolved culture filtrate proteins into 20 fractions with a pI range of 2.59 to 12.9. These 20 fractions were subjected to immunological analysis in healthy laboratory volunteers from our endemic area. Eleven fractions (Fractions 5, 6, 7, 8, 9, 10, 11, 13, 15, 16, and 19) showed increased interferon gamma (IFN-gamma) secretion and 5 fractions induced increased proliferative response, when compared to unfractionated CFP. In the 11 fractions which showed increased IFN-gamma secretion, mass spectrometric analysis identified 19 different proteins. Apart from the already reported immunodominant antigens like FbpB, CFP-10 and ESAT-6, two new T cell antigens (AcpM and PpiA) were also identified in the immunologically active fractions. Immunoinformatic analysis showed that PpiA was predicted to bind more number of class I and class II HLA alleles compared with the immunodominant ESAT-6 and CFP-10. Population coverage calculations also showed that PpiA protein (85%) had a higher population coverage compared with ESAT-6 (79%) and CFP-10 (77%). This result shows that the PpiA protein has a potential to be a novel T cell antigen.

HIV-1 transmitted drug resistance surveillance: shifting trends in study design and prevalence estimates.

INTRODUCTION: HIV-1 transmitted drug resistance (TDR) prevalence increased during the initial years of the antiretroviral therapy (ART) global scale-up. Few studies have examined recent trends in TDR prevalence using published genetic sequences and described the characteristics of ART-naïve persons from whom these published sequences have been obtained. METHODS: We identified 125 studies published between 2014 and 2019 for which HIV-1 reverse transcriptase (RT) with or without protease from ≥50 ART-naïve adult persons were submitted to the GenBank sequence database. The population characteristics and TDR prevalence were compared to those in 122 studies published in the preceding five years between 2009 and 2013. TDR prevalence was analysed using median study-level and person-level data. RESULTS AND DISCUSSION: The 2009 to 2013 and 2014 to 2019 studies reported sequence data from 32,866 and 41,724 ART-naïve persons respectively. Studies from the low- and middle-income country (LMIC) regions in sub-Saharan Africa, South/Southeast Asia and Latin America/Caribbean accounted for approximately two-thirds of the studies during each period. Between the two periods, the proportion of studies from sub-Saharan Africa and from South/Southeast Asia countries other than China decreased from 43% to 32% and the proportion of studies performed at sentinel sites for recent HIV-1 infection decreased from 33% to 22%. Between 2014 and 2019, median study-level TDR prevalence was 4.1% in South/Southeast Asia, 6.0% in sub-Saharan Africa, 9.1% in Latin America/Caribbean, 8.5% in Europe and 14.2% in North America. In the person-level analysis, there was an increase in overall, NNRTI and two-class NRTI/NNRTI resistance in sub-Saharan Africa; an increase in NNRTI resistance in Latin America/Caribbean, and an increase in overall, NNRTI and PI resistance in North America. CONCLUSIONS: Overall, NNRTI and dual NRTI/NNRTI-associated TDR prevalence was significantly higher in sub-Saharan Africa studies published between 2014 and 2019 compared with those published between 2009 and 2013. The decreasing proportion of studies from the hardest hit LMIC regions and the shift away from sentinel sites for recent infection suggests that global TDR surveillance efforts and publication of findings require renewed emphasis.

Mapping of Neutralizing Antibody Epitopes on the Envelope of Viruses Obtained from Plasma Samples Exhibiting Broad Cross-Clade Neutralization Potential Against HIV-1.

Several broadly neutralizing antibodies (bNAbs) that can target HIV strains with large degrees of variability have recently been identified. However, efforts to induce synthesis of such bNAbs that can protect against HIV infection have not met with much success. Identification of specific epitopes encoded in the HIV-1 envelope (Env) that can direct the host to synthesize bNAbs remains a challenge. In a previous study, we identified 12 antiretroviral therapy-naive HIV-1-infected individuals whose plasma exhibited broad cross-clade neutralization property against different clades of HIV-1. In this study, we sequenced the full-length HIV-1 gp160 from 11 of these individuals and analyzed the sequences to identify bNAb epitopes. We identified critical residues in the viral envelopes that contribute to the formation of conformational epitopes and possibly induce the production of bNAbs, using in silico methods. We found that many of the sequences had conserved glycans at positions N160 (10/11) and N332 (9/11), which are known to be critical for the binding of PG9/PG16-like and PGT128-like bNAbs, respectively. We also observed conservation of critical glycans at positions N234 and N276 critical for the interaction with CD4 binding site bNAbs in 8/11 and 11/11 sequences, respectively. We modeled the three-dimensional structure of the 11 HIV-1 envelopes and found that though each had structural differences, the critical residues were mostly present on the surface of the Env structures. The identified critical residues are proposed as candidates for further evaluation as bNAb epitopes.

Evidence for Highly Variable, Region-Specific Patterns of T-Cell Epitope Mutations Accumulating in Mycobacterium tuberculosis Strains.

Vaccines that confer protection through induction of adaptive T-cell immunity rely on understanding T-cell epitope (TCE) evolution induced by immune escape. This is poorly understood in tuberculosis (TB), an ancient, chronic disease, where CD4 T-cell immunity is of recognized importance. We probed 905 functionally validated, curated human CD4 T cell epitopes in 79 Mycobacterium tuberculosis (Mtb) whole genomes from India. This screen resulted in identifying 64 mutated epitopes in these strains initially using a computational pipeline and subsequently verified by single nucleotide polymorphism (SNP) analysis. SNP based phylogeny revealed the 79 Mtb strains to cluster to East African Indian (EAI), Central Asian Strain (CAS), and Beijing (BEI) lineages. Eighty-nine percent of the mutated T-cell epitopes (mTCEs) identified in the 79 Mtb strains from India has not previously been reported. These mTCEs were encoded by genes with high nucleotide diversity scores including seven mTCEs encoded by six antigens in the top 10% of rapidly divergent Mtb genes encoded by these strains. Using a T cell functional assay readout, we demonstrate 62% of mTCEs tested to significantly alter CD4 T-cell IFNγ and/or IL2 secretion with associated changes in predicted HLA-DR binding affinity: the gain of function mutations displayed higher predicted HLA-DR binding affinity and conversely mutations resulting in loss of function displayed lower predicted HLA-DR binding affinity. Most mutated antigens belonged to the cell wall/cell processes, and, intermediary metabolism and respiration families though all known Mtb proteins encoded mutations. Analysis of the mTCEs in an SNP database of 5,310 global Mtb strains identified 82% mTCEs to be significantly more prevalent in Mtb strains isolated from India, including 36 mTCEs identified exclusively in strains from India. These epitopes had a significantly higher predicted binding affinity to HLA-DR alleles that were highly prevalent in India compared to HLA-DR alleles rare in India, highlighting HLA-DR maybe an important driver of these mutations. This first evidence of region-specific TCE mutations potentially employed by Mtb to escape host immunity has important implications for TB vaccine design.

HLA based selection of epitopes offers a potential window of opportunity for vaccine design against HIV.

The pace of progression to AIDS after HIV infection varies from individual to individual. While some individuals develop AIDS quickly, others are protected from the onset of disease for more than a decade (elite controllers and long term non-progressors). The mechanisms of protection are not yet clearly understood, though various factors including host genetics, immune components and virus attenuation have been elucidated partly. The influence of HLA alleles on HIV-1 infection and disease outcome has been studied extensively. Several HLA alleles are known to be associated with resistance to infection or delayed progression to AIDS after infection. Similarly, certain HLA alleles are reported to be associated with rapid progression to disease. Since HLA alleles influence the outcome of HIV infection differentially, selection of epitopes specifically recognized by protective alleles could serve asa rational means for HIV vaccine design. In this review article, we discuss existing knowledge on HLA alleles and their association with resistance/susceptibility to HIV and its relevance to vaccine design.

Bioinformatics approach to prioritize known drugs towards repurposing for tuberculosis.

New drugs are urgently needed to cure tuberculosis (TB) in a short period of time without causing any adverse effects since currently used drugs for the treatment of multi drug-resistant TB cause several adverse effects with poor success rate. Therefore, we aimed to prioritize known drugs towards repurposing for TB by employing bioinformatics approach in the present study. A total of 1554 FDA approved drugs were obtained from DrugBank. Serine/threonine-protein kinase, pknB (Rv0014c) of Mycobacterium tuberculosis (Mtb) was selected as the drug target since it involves in several vital functions of the Mtb. All of the 1554 drugs were subjected to molecular docking with pknB. Glide and AutoDock Vina were employed using rigid docking followed by induced fit docking protocol for prioritization of drugs. Out of 14 drugs prioritized, six are suggested as high-confident drugs towards repurposing for TB as they were consistently found within top 10 ranks of both methods, and strongly binding in the active site of the pknB. We also found atorvastatin as one of the high-confident drugs, which has already been demonstrated to be active against Mtb under in vitro conditions by other researchers. Therefore, we propose that the prioritized six high-confident drugs as potential candidates for repurposing for TB and suggest for further experimental studies. We also suggest that the bioinformatics procedure we have employed in this study could be effectively applied for prioritization of drugs for other diseases.

Broad and potent cross clade neutralizing antibodies with multiple specificities in the plasma of HIV-1 subtype C infected individuals.

Broadly Cross clade Neutralizing (BCN) antibodies are recognized as potential therapeutic tools and leads for the design of a vaccine that can protect human beings against various clades of Human Immunodeficiency Virus (HIV). In the present study, we screened plasma of 88 HIV-1 infected ART naïve individuals for their neutralization potential using a standard panel of 18 pseudoviruses belonging to different subtypes and different levels of neutralization. We identified 12 samples with good breadth of neutralization (neutralized >90% of the viruses). Four of these samples neutralized even the difficult-to-neutralize tier-3 pseudoviruses with great potency (GMT > 600). Analysis of neutralization specificities indicated that four samples had antibodies with multiple epitope binding specificities, viz. CD4-binding site (CD4BS), glycans in the V1/V2 and V3 regions and membrane proximal external region (MPER). Our findings indicate the strong possibility of identifying highly potent bNAbs with known or novel specificities from HIV-1 subtype C infected individuals from India that can be exploited as therapeutic tools or lead molecules for the identification of potential epitopes for design of a protective HIV-1 vaccine.

Recent developments in genomics, bioinformatics and drug discovery to combat emerging drug-resistant tuberculosis.

Emergence of drug-resistant tuberculosis (DR-TB) is a big challenge in TB control. The delay in diagnosis of DR-TB leads to its increased transmission, and therefore prevalence. Recent developments in genomics have enabled whole genome sequencing (WGS) of Mycobacterium tuberculosis (M. tuberculosis) from 3-day-old liquid culture and directly from uncultured sputa, while new bioinformatics tools facilitate to determine DR mutations rapidly from the resulting sequences. The present drug discovery and development pipeline is filled with candidate drugs which have shown efficacy against DR-TB. Furthermore, some of the FDA-approved drugs are being evaluated for repurposing, and this approach appears promising as several drugs are reported to enhance efficacy of the standard TB drugs, reduce drug tolerance, or modulate the host immune response to control the growth of intracellular M. tuberculosis. Recent developments in genomics and bioinformatics along with new drug discovery collectively have the potential to result in synergistic impact leading to the development of a rapid protocol to determine the drug resistance profile of the infecting strain so as to provide personalized medicine. Hence, in this review, we discuss recent developments in WGS, bioinformatics and drug discovery to perceive how they would transform the management of tuberculosis in a timely manner.