Years of sequence feature curation by UniProtKB/Swiss-Prot, PIR-PSD, NCBI-CDD, RefSeq and other database biocurators has led to a rich repository of information on functional sites of genes and proteins. This information along with variation-related annotation can be used to scan human short sequence reads from next-generation sequencing (NGS) pipelines for presence of non-synonymous single-nucleotide variations (nsSNVs) that affect functional sites. This and similar workflows are becoming more important because thousands of NGS data sets are being made available through projects such as The Cancer Genome Atlas (TCGA), and researchers want to evaluate their biomarkers in genomic data. BioMuta, an integrated sequence feature database, provides a framework for automated and manual curation and integration of cancer-related sequence features so that they can be used in NGS analysis pipelines. Sequence feature information in BioMuta is collected from the Catalogue of Somatic Mutations in Cancer (COSMIC), ClinVar, UniProtKB and through biocuration of information available from publications. Additionally, nsSNVs identified through automated analysis of NGS data from TCGA are also included in the database. Because of the petabytes of data and information present in NGS primary repositories, a platform HIVE (High-performance Integrated Virtual Environment) for storing, analyzing, computing and curating NGS data and associated metadata has been developed. Using HIVE, 31 979 nsSNVs were identified in TCGA-derived NGS data from breast cancer patients. All variations identified through this process are stored in a Curated Short Read archive, and the nsSNVs from the tumor samples are included in BioMuta. Currently, BioMuta has 26 cancer types with 13 896 small-scale and 308 986 large-scale study-derived variations. Integration of variation data allows identifications of novel or common nsSNVs that can be prioritized in validation studies.Database URL: BioMuta: http://hive.biochemistry.gwu.edu/tools/biomuta/index.php; CSR: http://hive.biochemistry.gwu.edu/dna.cgi?cmd=csr; HIVE: http://hive.biochemistry.gwu.edu

ObjectiveAltered bacterial composition is associated with disease progression in cirrhosis but the role of virome, especially phages, is unclear.DesignCross-sectional and pre/post rifaximin cohorts were enrolled. Cross-sectional: controls and cirrhotic outpatients (compensated, on lactulose (Cirr-L), on rifaximin (Cirr-LR)) were included and followed for 90-day hospitalisations. Pre/post: compensated cirrhotics underwent stool collection pre/post 8 weeks of rifaximin. Stool metagenomics for bacteria and phages and their correlation networks were analysed in controls versus cirrhosis, within cirrhotics, hospitalised/not and pre/post rifaximin.ResultsCross-sectional: 40 controls and 163 cirrhotics (63 compensated, 43 Cirr-L, 57 Cirr-LR) were enrolled. Cirr-L/LR groups were similar on model for end-stage liver disease (MELD) score but Cirr-L developed greater hospitalisations versus Cirr-LR (56% vs 30%, p=0.008). Bacterial alpha/beta diversity worsened from controls through Cirr-LR. While phage alpha diversity was similar, beta diversity was different between groups. Autochthonous bacteria linked negatively, pathobionts linked positively with MELD but only modest phage-MELD correlations were seen. Phage–bacterial correlation network complexity was highest in controls, lowest in Cirr-L and increased in Cirr-LR. Microviridae and Faecalibacterium phages were linked with autochthonous bacteria in Cirr-LR, but not Cirr-L hospitalised patients had greater pathobionts, lower commensal bacteria and phages focused on Streptococcus, Lactococcus and Myoviridae. Pre/post: No changes in alpha/beta diversity of phages or bacteria were seen postrifaximin. Phage–bacterial linkages centred around urease-producing Streptococcus species collapsed postrifaximin.ConclusionUnlike bacteria, faecal phages are sparsely linked with cirrhosis characteristics and 90-day outcomes. Phage and bacterial linkages centred on urease-producing, ammonia-generating Streptococcus species were affected by disease progression and rifaximin therapy and were altered in patients who experienced 90-day hospitalisations.

BackgroundPolycystic ovarian syndrome (PCOS) is one of the most common reproductive disorders with strong association with both insulin resistance and non-alcoholic fatty liver disease (NAFLD). To untangle the complex relationship between PCOS and NAFLD, we analyzed serum biomarkers of apoptosis, some adipokines and mRNA profiles in the visceral adipose tissue of obese patients with NAFLD who were also diagnosed with PCOS and compared to a group with NAFLD only.MethodsWe included patients with biopsy-proven NAFLD and PCOS (N = 12) and BMI-matched biopsy-proven NAFLD patients without PCOS (N = 12). Expression levels of individual mRNAs and soluble serum biomarkers were compared by non-parametric Mann–Whitney test. The analysis also included Spearman rank correlation tests and multiple regression analysis. For co-correlated genes, the factor analysis was performed.ResultsThe total serum levels of apoptotic biomarker M30 were significantly elevated in PCOS patients with liver steatosis as compared to non-PCOS NAFLD controls (P < 0.02), pointing that androgen-dependent proapoptotic PCOS environment that may directly contribute to NAFLD progression in these patients. Similarly, hyperandrogenism may explain the observed PCOS-specific decrease (P < 0.04) in adipose LDLR mRNA expression that may be connected to the proneness of PCOS patients to NAFLD. The levels of mRNA encoding angiogenesis-associated GSK-3B interacting protein ninein were also significantly increased in the adipose tissue of NAFLD patients with PCOS (P < 0.007). Furthermore, the levels of resistin positively correlated with expression levels of LDLR and prothrombin time (PT).ConclusionAn androgen-dependent proapoptotic PCOS environment may directly contribute to NAFLD progression in these patients. Hyperandrogenism may explain an observed decrease in adipose LDLR mRNA expression. An inflammation-associated increase in the release of resistin into circulation might contribute to the prothrombotic state observed under conditions associated with insulin resistance, including PCOS. The studies of larger cohorts of NAFLD with and without PCOS patients are needed to further assess these potential interactions.

Antibiotic resistance leads to poor outcomes in cirrhosis. Fecal microbiota transplant (FMT) is associated with reduction in antibiotic resistance gene (ARG) burden in patients without cirrhosis; however, the impact in cirrhosis is unclear. We aimed to study the effect of capsule and enema FMT on ARG abundance in fecal samples, which were collected during two published FMT trials in patients with cirrhosis on rifaximin, lactulose, and proton pump inhibitors. ARGs were identified using metagenomics and mapped against the Comprehensive Antibiotic Resistance Database. Changes in ARG abundance were studied within/between groups. The capsule FMT trial involved a one-time FMT or placebo capsule administration with stool collection at baseline and week 4 postintervention. Antibiotics+enema FMT included preprocedure antibiotics followed by FMT enema versus standard-of-care (SOC). Stool was collected at baseline, postantibiotics, and day 7/15 postintervention. Both trials included 20 patients each. There was no safety/ infection signal linked to FMT. In the capsule trial, beta-lactamase (OXY/LEN) expression decreased post-FMT versus baseline. Compared to placebo, patients who were post-FMT had lower abundance of vancomycin (VanH), betalactamase (ACT), and rifamycin ARGs; the latter was associated with cognitive improvement. No changes were seen within patients treated with placebo. In the antibiotics+enema trial for postantibiotics at day 7 versus baseline, there was an increase in vancomycin and beta-lactamase ARGs, which decreased at day 15. However, quinolone resistance increased at day 15 versus baseline. Between SOC and FMT, day 7 had largely lower ARG (CfxA beta-lactamase, VanW, and VanX) that continued at day 15 (cepA beta-lactamase, VanW). No changes were seen within the SOC group. Conclusion: Despite differences in routes of administration and preintervention antibiotics, we found that ARG abundance is largely reduced after FMT compared to pre-FMT baseline and non-FMT groups in decompensated cirrhosis. (Hepatology Communications 2020;0:1-14).

BACKGROUND AND AIMS: Cirrhosis is associated with changes in intestinal microbiota that can lead to hepatic encephalopathy (HE) and infections, especially with antibioticresistant organisms. However, the impact of gut microbial antibiotic resistance gene (ARG) burden on clinical outcomes is unclear. The aims of the study were to determine the impact of ARGs in cirrhosis-related gut metagenome on outcomes and disease progression, study the effect of rifaximin on ARG burden, and compare ARGs in cirrhosis with chronic kidney disease (CKD) and diabetes. METHODS: In outpatients with cirrhosis who underwent metagenomics, we evaluated change in ARG abundances with progression and their multivariable impact on 90-day hospitalizations and deaths over 1 year. We also studied ARGs pre-and 8 weeks post-rifaximin in patients with compensated cirrhosis in an open-label trial. Finally, ARGs from CKD and diabetes studies were compared with cirrhosis on machine learning. RESULTS: A total of 163 patients with cirrhosis (43 compensated, 20 ascites-only, 30 HE-only, 70 both) and 40 controls were included. ARG abundances were higher in cirrhosis versus controls and worsened with advancing cirrhosis severity; 44 patients were hospitalized and 14 died. ARG abundances were associated with hospitalizations and mortality while controlling for cirrhosis complications, medications, and demographics. Rifaximin trial: ARG abundance patterns were minimally affected in 19 patients post-rifaximin. CKD/diabetes comparison: ARG abundance patterns in cirrhosis are distinguishable on machine learning and include more gram-positive ARGs. CONCLUSIONS: Cirrhosis is associated with high gut microbial ARG gene burden compared with controls, which worsens with disease progression and may be different from CKD and diabetes. ARGs are not affected by rifaximin and are associated with hospitalizations and death.

Human adenovirus type 14 (HAdV-B14p) was originally identified as an acute respiratory disease (ARD) pathogen in The Netherlands in 1955. For approximately fifty years, few sporadic infections were observed. In 2005, HAdV-B14p1, a genomic variant, re-emerged and was associated with several large ARD outbreaks across the U.S. and, subsequently, in Canada, the U.K., Ireland, and China. This strain was associated with an unusually higher fatality rate than previously reported for both this prototype and other HAdV types in general. In China, HAdV-B14 was first observed in 2010, when two unrelated HAdV-B14-associated ARD cases were reported in Southern China (GZ01) and Northern China (BJ430), followed by three subsequent outbreaks. While comparative genomic analysis, including indel analysis, shows that the three China isolates, with whole genome data available, are similar to the de Wit prototype, all are divergent from the U.S. strain (303600; 2007). Although the genomes of strains GZ01 and BJ430 are nearly identical, as per their genome type characterization and percent identities, they are subtly divergent in their genome mutation patterns. These genomes indicate possibly two lineages of HAdV-B14 and independent introductions into China from abroad, or subsequent divergence from one; CHN2012 likely represents a separate sub-lineage. Observations of these simultaneously reported emergent strains in China add to the understanding of the circulation, epidemiology, and evolution of these HAdV pathogens, as well as provide a foundation for developing effective vaccines and public health strategies, including nationwide surveillance in anticipation of larger outbreaks with potentially higher fatality rates associated with HAdV-B14p1.

Background A number of serious human adenovirus (HAdV) outbreaks have been recently reported: HAdV-B7 (Israel, Singapore, and USA), HAdV-B7d (USA and China), HAdV-D8, -D54, and -C2 (Japan), HAdV-B14p1 (USA, Europe, and China), and HAdV-B55 (China, Singapore, and France). Methods To understand the epidemiology of HAdV infections in Singapore, we studied 533 HAdV-positive clinical samples collected from 396 pediatric and 137 adult patients in Singapore from 2012 to 2018. Genome sequencing and phylogenetic analyses were performed to identify HAdV genotypes, clonal clusters, and recombinant or novel HAdVs. Results The most prevalent genotypes identified were HAdV-B3 (35.6%), HAdV-B7 (15.4%), and HAdV-E4 (15.2%). We detected 4 new HAdV-C strains and detected incursions with HAdV-B7 (odds ratio [OR], 14.6; 95% confidence interval [CI], 4.1–52.0) and HAdV-E4 (OR, 13.6; 95% CI, 3.9–46.7) among pediatric patients over time. In addition, immunocompromised patients (adjusted OR [aOR], 11.4; 95% CI, 3.8–34.8) and patients infected with HAdV-C2 (aOR, 8.5; 95% CI, 1.5–48.0), HAdV-B7 (aOR, 3.7; 95% CI, 1.2–10.9), or HAdV-E4 (aOR, 3.2; 95% CI, 1.1–8.9) were at increased risk for severe disease. Conclusions Singapore would benefit from more frequent studies of clinical HAdV genotypes to identify patients at risk for severe disease and help guide the use of new antiviral therapies, such as brincidofovir, and potential administration of HAdV 4 and 7 vaccine.

Identification of non-synonymous single nucleotide variations (nsSNVs) has exponentially increased due to advances in Next-Generation Sequencing technologies. The functional impacts of these variations have been difficult to ascertain because the corresponding knowledge about sequence functional sites is quite fragmented. It is clear that mapping of variations to sequence functional features can help us better understand the pathophysiological role of variations. In this study, we investigated the effect of nsSNVs on more than 17 common types of post-translational modification (PTM) sites, active sites and binding sites. Out of 1 705 285 distinct nsSNVs on 259 216 functional sites we identified 38 549 variations that significantly affect 10 major functional sites. Furthermore, we found distinct patterns of site disruptions due to germline and somatic nsSNVs. Pan-cancer analysis across 12 different cancer types led to the identification of 51 genes with 106 nsSNV affected functional sites found in 3 or more cancer types. 13 of the 51 genes overlap with previously identified Significantly Mutated Genes (Nature. 2013 Oct 17;502(7471)). 62 mutations in these 13 genes affecting functional sites such as DNA, ATP binding and various PTM sites occur across several cancers and can be prioritized for additional validation and investigations.