Genomic Landscape Survey Identifies SRSF1 as a Key Oncodriver in Small Cell Lung Cancer.

Small cell lung cancer (SCLC) is an aggressive disease with poor survival. A few sequencing studies performed on limited number of samples have revealed potential disease-driving genes in SCLC, however, much still remains unknown, particularly in the Asian patient population. Here we conducted whole exome sequencing (WES) and transcriptomic sequencing of primary tumors from 99 Chinese SCLC patients. Dysregulation of tumor suppressor genes TP53 and RB1 was observed in 82% and 62% of SCLC patients, respectively, and more than half of the SCLC patients (62%) harbored TP53 and RB1 mutation and/or copy number loss. Additionally, Serine/Arginine Splicing Factor 1 (SRSF1) DNA copy number gain and mRNA over-expression was strongly associated with poor survival using both discovery and validation patient cohorts. Functional studies in vitro and in vivo demonstrate that SRSF1 is important for tumorigenicity of SCLC and may play a key role in DNA repair and chemo-sensitivity. These results strongly support SRSF1 as a prognostic biomarker in SCLC and provide a rationale for personalized therapy in SCLC.

BubbleTree: an intuitive visualization to elucidate tumoral aneuploidy and clonality using next generation sequencing data.

Tumors are characterized by properties of genetic instability, heterogeneity, and significant oligoclonality. Elucidating this intratumoral heterogeneity is challenging but important. In this study, we propose a framework, BubbleTree, to characterize the tumor clonality using next generation sequencing (NGS) data. BubbleTree simultaneously elucidates the complexity of a tumor biopsy, estimating cancerous cell purity, tumor ploidy, allele-specific copy number, and clonality and represents this in an intuitive graph. We further developed a three-step heuristic method to automate the interpretation of the BubbleTree graph, using a divide-and-conquer strategy. In this study, we demonstrated the performance of BubbleTree with comparisons to similar commonly used tools such as THetA2, ABSOLUTE, AbsCN-seq and ASCAT, using both simulated and patient-derived data. BubbleTree outperformed these tools, particularly in identifying tumor subclonal populations and polyploidy. We further demonstrated BubbleTree's utility in tracking clonality changes from patients' primary to metastatic tumor and dating somatic single nucleotide and copy number variants along the tumor clonal evolution. Overall, the BubbleTree graph and corresponding model is a powerful approach to provide a comprehensive spectrum of the heterogeneous tumor karyotype in human tumors. BubbleTree is R-based and freely available to the research community (https://www.bioconductor.org/packages/release/bioc/html/BubbleTree.html).

HERC5 is a prognostic biomarker for post-liver transplant recurrent human hepatocellular carcinoma.

BACKGROUND AND AIMS: Orthotopic liver transplantation (OLT) can be an effective treatment option for certain patients with early stage hepatocellular carcinoma (HCC) meeting Milan, UCSF, or Hangzhou criteria. However, HCC recurrence rates post-OLT range from 20 to 40 %, with limited follow-up options. Elucidating genetic drivers common to primary and post-OLT recurrent tumors may further our understanding and help identify predictive biomarkers of recurrence-both to ultimately help manage clinical decisions for patients undergoing OLT. METHODS: Whole exome and RNA sequencing in matched primary and recurrent tumors, normal adjacent tissues, and blood from four Chinese HCC patients was conducted. SiRNA knockdown and both qRT-PCR and Western assays were performed on PLCPRF5, SNU449 and HEPG2 cell lines; immunohistochemistry and RNA Sequencing were conducted on the primary tumors of Chinese HCC patients who experienced tumor recurrence post-OLT (n = 9) or did not experience tumor recurrence (n = 12). RESULTS: In three independent HCC studies of patients undergoing transplantation (n = 21) or surgical resection (n = 242, n = 44) of primary tumors (total n = 307), HERC5 mRNA under-expression correlated with shorter: time to tumor recurrence (p = 0.007 and 0.02) and overall survival (p = 0.0063 and 0.023), even after adjustment for relevant clinical variables. HERC5 loss drives CCL20 mRNA and protein over-expression and associates with regulatory T cell infiltration as measured by FOXP3 expression. Further, matched primary and recurrent tumors from the 4 HCC patients indicated clonal selection advantage of Wnt signaling activation and CDKN2A inactivation. CONCLUSIONS: HERC5 plays a crucial role in HCC immune evasion and has clinical relevance as a reproducible prognostic marker for risk of tumor recurrence and survival in patients.

New rapid scheme for distinguishing the subspecies of the Mycobacterium abscessus group and identifying Mycobacterium massiliense isolates with inducible clarithromycin resistance.

Mycobacterium abscessus (M. abscessus sensu lato, or the M. abscessus group) comprises three closely related taxa whose taxonomic statuses are under revision, i.e., M. abscessus sensu stricto, Mycobacterium bolletii, and Mycobacterium massiliense. We describe here a simple, robust, and cost-effective PCR-based method for distinguishing among M. abscessus, M. massiliense, and M. bolletii. Based on the M. abscessus ATCC 19977(T) genome, regions that discriminated between M. abscessus and M. massiliense were identified through array-based comparative genomic hybridization. A typing scheme using PCR primers designed for four of these locations was applied to 46 well-characterized clinical isolates comprising 29 M. abscessus, 15 M. massiliense, and 2 M. bolletii isolates previously identified by multitarget sequencing. Interestingly, 2 isolates unequivocally identified as M. massiliense were shown to have a full-length erm(41) gene instead of the expected gene deletion and showed inducible clarithromycin resistance after 14 days. We propose using this PCR-based typing scheme combined with erm(41) PCR for straightforward identification of M. abscessus, M. massiliense, and M. bolletii and the assessment of inducible clarithromycin resistance. This method can be easily integrated into a routine workflow to provide subspecies-level identification within 24 h after isolation of the M. abscessus group.

The Protein Naming Utility: a rules database for protein nomenclature.

Generation of syntactically correct and unambiguous names for proteins is a challenging, yet vital task for functional annotation processes. Proteins are often named based on homology to known proteins, many of which have problematic names. To address the need to generate high-quality protein names, and capture our significant experience correcting protein names manually, we have developed the Protein Naming Utility (PNU, http://www.jcvi.org/pn-utility). The PNU is a web-based database for storing and applying naming rules to identify and correct syntactically incorrect protein names, or to replace synonyms with their preferred name. The PNU allows users to generate and manage collections of naming rules, optionally building upon the growing body of rules generated at the J. Craig Venter Institute (JCVI). Since communities often enforce disparate conventions for naming proteins, the PNU supports grouping rules into user-managed collections. Users can check their protein names against a selected PNU rule collection, generating both statistics and corrected names. The PNU can also be used to correct GenBank table files prior to submission to GenBank. Currently, the database features 3080 manual rules that have been entered by JCVI Bioinformatics Analysts as well as 7458 automatically imported names.

Pathema: a clade-specific bioinformatics resource center for pathogen research.

Pathema (http://pathema.jcvi.org) is one of the eight Bioinformatics Resource Centers (BRCs) funded by the National Institute of Allergy and Infectious Disease (NIAID) designed to serve as a core resource for the bio-defense and infectious disease research community. Pathema strives to support basic research and accelerate scientific progress for understanding, detecting, diagnosing and treating an established set of six target NIAID Category A-C pathogens: Category A priority pathogens; Bacillus anthracis and Clostridium botulinum, and Category B priority pathogens; Burkholderia mallei, Burkholderia pseudomallei, Clostridium perfringens and Entamoeba histolytica. Each target pathogen is represented in one of four distinct clade-specific Pathema web resources and underlying databases developed to target the specific data and analysis needs of each scientific community. All publicly available complete genome projects of phylogenetically related organisms are also represented, providing a comprehensive collection of organisms for comparative analyses. Pathema facilitates the scientific exploration of genomic and related data through its integration with web-based analysis tools, customized to obtain, display, and compute results relevant to ongoing pathogen research. Pathema serves the bio-defense and infectious disease research community by disseminating data resulting from pathogen genome sequencing projects and providing access to the results of inter-genomic comparisons for these organisms.

Comparative systeomics to elucidate physiological differences between CHO and SP2/0 cell lines.

Omics-based tools were coupled with bioinformatics for a systeomics analysis of two biopharma cell types: Chinese hamster ovary (M-CHO and CHO-K1) and SP2/0. Exponential and stationary phase samples revealed more than 10,000 transcripts and 6000 proteins across these two manufacturing cell lines. A statistical comparison of transcriptomics and proteomics data identified downregulated genes involved in protein folding, protein synthesis and protein metabolism, including PPIA-cyclophilin A, HSPD1, and EIF3K, in M-CHO compared to SP2/0 while cell cycle and actin cytoskeleton genes were reduced in SP2/0. KEGG pathway comparisons revealed glycerolipids, glycosphingolipids, ABC transporters, calcium signaling, cell adhesion, and secretion pathways depleted in M-CHO while retinol metabolism was upregulated. KEGG and IPA also indicated apoptosis, RNA degradation, and proteosomes enriched in CHO stationary phase. Alternatively, gene ontology analysis revealed an underrepresentation in ion and potassium channel activities, membrane proteins, and secretory granules including Stxbpt2, Syt1, Syt9, and Cma1 proteins in M-CHO. Additional enrichment strategies involving ultracentrifugation, biotinylation, and hydrazide chemistry identified over 4000 potential CHO membrane and secretory proteins, yet many secretory and membrane proteins were still depleted. This systeomics pipeline has revealed bottlenecks and potential opportunities for cell line engineering in CHO and SP2/0 to improve their production capabilities.

High-throughput mediation analysis of human proteome and metabolome identifies mediators of post-bariatric surgical diabetes control.

To improve the power of mediation in high-throughput studies, here we introduce High-throughput mediation analysis (Hitman), which accounts for direction of mediation and applies empirical Bayesian linear modeling. We apply Hitman in a retrospective, exploratory analysis of the SLIMM-T2D clinical trial in which participants with type 2 diabetes were randomized to Roux-en-Y gastric bypass (RYGB) or nonsurgical diabetes/weight management, and fasting plasma proteome and metabolome were assayed up to 3 years. RYGB caused greater improvement in HbA1c, which was mediated by growth hormone receptor (GHR). GHR's mediation is more significant than clinical mediators, including BMI. GHR decreases at 3 months postoperatively alongside increased insulin-like growth factor binding proteins IGFBP1/BP2; plasma GH increased at 1 year. Experimental validation indicates (1) hepatic GHR expression decreases in post-bariatric rats; (2) GHR knockdown in primary hepatocytes decreases gluconeogenic gene expression and glucose production. Thus, RYGB may induce resistance to diabetogenic effects of GH signaling.Trial Registration: Clinicaltrials.gov NCT01073020.

Genome sequence of the deep-rooted Yersinia pestis strain Angola reveals new insights into the evolution and pangenome of the plague bacterium.

To gain insights into the origin and genome evolution of the plague bacterium Yersinia pestis, we have sequenced the deep-rooted strain Angola, a virulent Pestoides isolate. Its ancient nature makes this atypical isolate of particular importance in understanding the evolution of plague pathogenicity. Its chromosome features a unique genetic make-up intermediate between modern Y. pestis isolates and its evolutionary ancestor, Y. pseudotuberculosis. Our genotypic and phenotypic analyses led us to conclude that Angola belongs to one of the most ancient Y. pestis lineages thus far sequenced. The mobilome carries the first reported chimeric plasmid combining the two species-specific virulence plasmids. Genomic findings were validated in virulence assays demonstrating that its pathogenic potential is distinct from modern Y. pestis isolates. Human infection with this particular isolate would not be diagnosed by the standard clinical tests, as Angola lacks the plasmid-borne capsule, and a possible emergence of this genotype raises major public health concerns. To assess the genomic plasticity in Y. pestis, we investigated the global gene reservoir and estimated the pangenome at 4,844 unique protein-coding genes. As shown by the genomic analysis of this evolutionary key isolate, we found that the genomic plasticity within Y. pestis clearly was not as limited as previously thought, which is strengthened by the detection of the largest number of isolate-specific single-nucleotide polymorphisms (SNPs) currently reported in the species. This study identified numerous novel genetic signatures, some of which seem to be intimately associated with plague virulence. These markers are valuable in the development of a robust typing system critical for forensic, diagnostic, and epidemiological studies.

Draft genome sequences of Yersinia pestis isolates from natural foci of endemic plague in China.

To gain insights into the evolutionary origin, emergence, and pathogenicity of the etiologic agent of plague, we have sequenced the genomes of four Yersinia pestis strains isolated from the zoonotic rodent reservoir in foci of endemic plague in China. These resources enable in-depth studies of Y. pestis sequence variations and detailed whole-genome comparisons of very closely related genomes from the supposed site of the origin and the emergence of global pandemics of plague.

Baseline Plasma Cell Gene Signature Predicts Improvement in Systemic Sclerosis Skin Scores Following Treatment With Inebilizumab (MEDI-551) and Correlates With Disease Activity in Systemic Lupus Erythematosus and Chronic Obstructive Pulmonary Disease.

OBJECTIVE: B cells impact the progression of systemic sclerosis (SSc; scleroderma) through multiple pathogenic mechanisms. CD19 inhibition in mice reduced skin thickness, collagen production, and autoantibody levels, consistent with CD19 expression on plasma cells (PCs), the source of antibody production. PC depletion could effectively reduce collagen deposition and inflammation in SSc; therefore, we investigated the effects of PC depletion on SSc disease activity. METHODS: A PC gene signature was evaluated in SSc skin biopsy samples in 2 phase I clinical trials. We assessed microarray data from tissue from public studies of chronic obstructive pulmonary disease (COPD), idiopathic pulmonary fibrosis (IPF), dermatomyositis (DM), systemic lupus erythematosus (SLE), and atopic dermatitis, as well as blood from a phase IIb clinical trial in SLE. RESULTS: The PC signature was elevated in SSc skin specimens compared to healthy donor skin (P = 2.28 × 10-6 ) and correlated with the baseline modified Rodnan skin thickness score (MRSS) (r = 0.64, P = 0.0004). Patients with a high PC signature at baseline showed greater improvement in the MRSS (mean ± SD change 35 ± 16%; P = 6.30 × 10-4 ) following anti-CD19 treatment with inebilizumab (MEDI-551) than did patients with a low PC signature at baseline (mean ± SD change 8 ± 12%; P = 0.104). The PC signature was overexpressed in tissue from patients with SLE, DM, COPD, interstitial lung disease, and IPF relative to controls (all fold change >2; P < 0.001). The PC signature also differed significantly between SLE patients with mild-to-moderate disease and those with severe disease (SLE Disease Activity Index cutoff at 10) (fold change 1.44; P = 3.90 × 10-3 ) and correlated significantly with the degree of emphysema in COPD (r = 0.53, P = 7.55 × 10-8 ). CONCLUSION: Our results support the notion that PCs have a role in the pathogenesis of SSc and other autoimmune or pulmonary indications. An elevated pretreatment PC signature was associated with increased benefit from MEDI-551 in SSc.

High-throughput RNA sequencing reveals distinct gene signatures in active IgG4-related disease.

We aimed to characterize the molecular differences and effects from prednisone treatment among IgG4-related disease with salivary gland lesions (RD-SG), without SG lesions (RD-nonSG), and IgG4-related retroperitoneal fibrosis (RF). RNA sequencing was conducted on blood from 25 RD-SG, 11 RD-nonSG, 3 RF and 10 control subjects. Among these, 8 RD-nonSG and 12 RD-SG patients were subjected to treatment with prednisone and/or glucocorticoid-sparing agents. Six RD patients had a longitudinal time point. The mRNA levels of IgG4 and IgE, genes specific for Th2 cells, eosinophils, and neutrophils were over-expressed in RD-SG and RD-nonSG. A B-cell signature was suppressed in patients group versus controls, while Th1, Th2, Treg, and eosinophil gene signatures were increased in patients without treatment. Interestingly, Tfh genes and B cell signature were decreased at flare disease state. Prednisone treatment led to increased neutrophil, but decreased Treg signatures. Serum IgG4 levels correlated with the eosinophil and neutrophil gene signatures in RD-SG patients, and with a B cell signature in only RD-nonSG patients. IgG4, IgE, and cell-specific signatures are regulated in patients, suggesting the imbalance of immune and inflammatory cells in IgG4-related disease. Prednisone treatment selectively modulates Treg, eosinophil, and neutrophil signatures.

A high density of tertiary lymphoid structure B cells in lung tumors is associated with increased CD4+ T cell receptor repertoire clonality.

T and B cell receptor (TCR and BCR, respectively) Vß or immunoglobulin heavy chain complementarity-determining region 3 sequencing allows monitoring of repertoire changes through recognition, clonal expansion, affinity maturation, and T or B cell activation in response to antigen. TCR and BCR repertoire analysis can advance understanding of antitumor immune responses in the tumor microenvironment. TCR and BCR repertoires of sorted CD4+, CD8+ or CD19+ cells in tumor, non-tumoral distant tissue (NT), and peripheral compartments (blood/draining lymph node [P]) from 47 non-small cell lung cancer (NSCLC) patients (agemedian = 68 y) were sequenced. The clonotype spectra were assessed among different tissues and correlated with clinical and immunological parameters. In all tissues, CD4+ and CD8+ TCR repertoires had greater clonality relative to CD19+ BCR. CD4+ T cells exhibited greater clonality in NT compared to tumor (p = 0.002) and P (p < 0.001), concentrated among older patients (age > 68). Younger patients exhibited greater CD4+ T cell diversity in P compared to older patients (p = 0.05), and greater CD4+ T cell clonality in tumor relative to P (p < 0.001), with fewer shared clonotypes between tumor and P than older patients (p = 0.04). More interestingly, greater CD4+ and CD8+ T cell clonality in tumor and P, respectively (both p = 0.05), correlated with high density of tumor-associated tertiary lymphoid structure (TLS) B cells, a biomarker of higher overall survival in NSCLC. Results indicate distinct adaptive immune responses in NSCLC, where peripheral T cell diversity is modulated by age, and tumor T cell clonal expansion is favored by the presence of TLSs in the tumor microenvironment.