Pediatric brain tumor cell lines exhibit miRNA-depleted, Y RNA-enriched extracellular vesicles.

BACKGROUND: Medulloblastoma (MB) and diffuse infiltrative pontine glioma (DIPG) are malignant pediatric tumors. Extracellular vesicles (EVs) and their bioactive cargoes have been implicated in tumorigenesis. Most studies have focused on adult tumors, therefore the role of EVs and the noncoding RNA (ncRNA) landscape in pediatric brain tumors is not fully characterized. The overall aim of this pilot study was to isolate EVs from MB and DIPG patient-derived cell lines and to explore the small ncRNA transcriptome. METHODS: EVs from 3 DIPG and 4 MB patient-derived cell lines were analyzed. High-throughput next generation sequencing interrogated the short non-coding RNA (ncRNA) transcriptome. Known and novel miRNAs were quantified. Differential expression analysis, in silico target prediction, and functional gene enrichment were performed. RESULTS: EV secretomes from MB and DIPG patient-derived cell lines demonstrated discrete ncRNA biotypes. Notably, miRNAs were depleted and Y RNAs were enriched in EV samples. Hierarchical cluster analysis revealed high discrimination in miRNA expression between DIPG and MB cell lines and RNA-Seq identified novel miRNAs not previously implicated in MB or DIPG pathogenesis. Known and putative target genes of dysregulated miRNAs were identified. Functional annotation analysis of the target genes for differentially expressed EV-and parental-derived miRNAs revealed significant cancer-related pathway involvement. CONCLUSIONS: This hypothesis-generating study demonstrated that pediatric brain tumor-derived cell lines secrete EVs comprised of various ncRNA cargoes. Validation of these findings in patient samples may provide new insights into the pediatric brain tumor microenvironment and identification of novel therapeutic candidates.

Natural underlying mtDNA heteroplasmy as a potential source of intra-person hiPSC variability.

Functional variability among human clones of induced pluripotent stem cells (hiPSCs) remains a limitation in assembling high-quality biorepositories. Beyond inter-person variability, the root cause of intra-person variability remains unknown. Mitochondria guide the required transition from oxidative to glycolytic metabolism in nuclear reprogramming. Moreover, mitochondria have their own genome (mitochondrial DNA [mtDNA]). Herein, we performed mtDNA next-generation sequencing (NGS) on 84 hiPSC clones derived from a cohort of 19 individuals, including mitochondrial and non-mitochondrial patients. The analysis of mtDNA variants showed that low levels of potentially pathogenic mutations in the original fibroblasts are revealed through nuclear reprogramming, generating mutant hiPSCs with a detrimental effect in their differentiated progeny. Specifically, hiPSC-derived cardiomyocytes with expanded mtDNA mutations non-related with any described human disease, showed impaired mitochondrial respiration, being a potential cause of intra-person hiPSC variability. We propose mtDNA NGS as a new selection criterion to ensure hiPSC quality for drug discovery and regenerative medicine.

Short non-coding RNA sequencing of glioblastoma extracellular vesicles.

INTRODUCTION: Like all nucleated cells, glioblastoma (GBM) cells shed small membrane-encapsulated particles called extracellular vesicles (EVs). EVs can transfer oncogenic components and promote tumor growth by transferring short non-coding RNAs, altering target cell gene expression. Furthermore, GBM-derived EVs can be detected in blood and have potential to serve as liquid biopsies. METHODS: EVs were harvested from culture supernatants from human GBM cell lines, purified via sequential centrifugation, and quantified by nanoparticle tracking. RNA was isolated and short non-coding RNA was sequenced. Data was analyzed via the OASIS-2.0 platform using HG38. MirTarBase and MirDB interrogated validated/predicted miRNA-gene interactions respectively. RESULTS: Many short non-coding RNA's were identified within GBM EV's. In keeping with earlier reports utilizing GBM EV micro-RNA (miRNA) arrays, these included abundant micro-RNA's including miR-21. However, RNA sequencing revealed a total of 712 non-coding RNA sequences most of which have not been associated with GBM EV's previously. These included many RNA species (piRNA, snoRNA, snRNA, rRNA and yRNAs) in addition to miRNA's. miR-21-5p, let-7b-5p, miR-3182, miR-4448, let-7i-5p constituted highest overall expression. Top genes targeted by non-coding RNA's were highly conserved and specific for cell cycle, PI3K/Akt signaling, p53 and Glioma curated KEGG pathways. CONCLUSIONS: Next generation short non-coding RNA sequencing on GBM EV's validates findings from earlier studies using miRNA arrays but also demonstrates expression of many additional non-coding RNA sequences and classes previously unassociated with GBM. This may yield important insights into pathophysiology, point to new therapeutic targets, and help develop new biomarkers for disease burden and treatment response.

Long-lived particles at the energy frontier: the MATHUSLA physics case.

We examine the theoretical motivations for long-lived particle (LLP) signals at the LHC in a comprehensive survey of standard model (SM) extensions. LLPs are a common prediction of a wide range of theories that address unsolved fundamental mysteries such as naturalness, dark matter, baryogenesis and neutrino masses, and represent a natural and generic possibility for physics beyond the SM (BSM). In most cases the LLP lifetime can be treated as a free parameter from the [Formula: see text]m scale up to the Big Bang Nucleosynthesis limit of [Formula: see text] m. Neutral LLPs with lifetimes above [Formula: see text]100 m are particularly difficult to probe, as the sensitivity of the LHC main detectors is limited by challenging backgrounds, triggers, and small acceptances. MATHUSLA is a proposal for a minimally instrumented, large-volume surface detector near ATLAS or CMS. It would search for neutral LLPs produced in HL-LHC collisions by reconstructing displaced vertices (DVs) in a low-background environment, extending the sensitivity of the main detectors by orders of magnitude in the long-lifetime regime. We study the LLP physics opportunities afforded by a MATHUSLA-like detector at the HL-LHC, assuming backgrounds can be rejected as expected. We develop a model-independent approach to describe the sensitivity of MATHUSLA to BSM LLP signals, and compare it to DV and missing energy searches at ATLAS or CMS. We then explore the BSM motivations for LLPs in considerable detail, presenting a large number of new sensitivity studies. While our discussion is especially oriented towards the long-lifetime regime at MATHUSLA, this survey underlines the importance of a varied LLP search program at the LHC in general. By synthesizing these results into a general discussion of the top-down and bottom-up motivations for LLP searches, it is our aim to demonstrate the exceptional strength and breadth of the physics case for the construction of the MATHUSLA detector.

Recessive TAF1A mutations reveal ribosomopathy in siblings with end-stage pediatric dilated cardiomyopathy.

Non-ischemic dilated cardiomyopathy (DCM) has been recognized as a heritable disorder for over 25 years, yet clinical genetic testing is non-diagnostic in >50% of patients, underscoring the ongoing need for DCM gene discovery. Here, whole exome sequencing uncovered a novel molecular basis for idiopathic end-stage heart failure in two sisters who underwent cardiac transplantation at three years of age. Compound heterozygous recessive mutations in TAF1A, encoding an RNA polymerase I complex protein, were associated with marked fibrosis of explanted hearts and gene-specific nucleolar segregation defects in cardiomyocytes, indicative of impaired ribosomal RNA synthesis. Knockout of the homologous gene in zebrafish recapitulated a heart failure phenotype with pericardial edema, decreased ventricular systolic function, and embryonic mortality. These findings expand the clinical spectrum of ribosomopathies to include pediatric DCM.

Noncardiac genetic predisposition in sudden infant death syndrome.

PURPOSE: Sudden infant death syndrome (SIDS) is the commonest cause of sudden death of an infant; however, the genetic basis remains poorly understood. We aimed to identify noncardiac genes underpinning SIDS and determine their prevalence compared with ethnically matched controls. METHODS: Using exome sequencing we assessed the yield of ultrarare nonsynonymous variants (minor allele frequency [MAF] ≤0.00005, dominant model; MAF ≤0.01, recessive model) in 278 European SIDS cases (62% male; average age =2.7 ± 2 months) versus 973 European controls across 61 noncardiac SIDS-susceptibility genes. The variants were classified according to American College of Medical Genetics and Genomics criteria. Case-control, gene-collapsing analysis was performed in eight candidate biological pathways previously implicated in SIDS pathogenesis. RESULTS: Overall 43/278 SIDS cases harbored an ultrarare single-nucleotide variant compared with 114/973 controls (15.5 vs. 11.7%, p=0.10). Only 2/61 noncardiac genes were significantly overrepresented in cases compared with controls (ECE1, 3/278 [1%] vs. 1/973 [0.1%] p=0.036; SLC6A4, 2/278 [0.7%] vs. 1/973 [0.1%] p=0.049). There was no difference in yield of pathogenic or likely pathogenic variants between cases and controls (1/278 [0.36%] vs. 4/973 [0.41%]; p=1.0). Gene-collapsing analysis did not identify any specific biological pathways to be significantly associated with SIDS. CONCLUSIONS: A monogenic basis for SIDS amongst the previously implicated noncardiac genes and their encoded biological pathways is negligible.

A Comprehensive Approach to Sequence-oriented IsomiR annotation (CASMIR): demonstration with IsomiR profiling in colorectal neoplasia.

BACKGROUND: MicroRNA (miRNA) profiling is an important step in studying biological associations and identifying marker candidates. miRNA exists in isoforms, called isomiRs, which may exhibit distinct properties. With conventional profiling methods, limitations in assay and analysis platforms may compromise isomiR interrogation. RESULTS: We introduce a comprehensive approach to sequence-oriented isomiR annotation (CASMIR) to allow unbiased identification of global isomiRs from small RNA sequencing data. In this approach, small RNA reads are maintained as independent sequences instead of being summarized under miRNA names. IsomiR features are identified through step-wise local alignment against canonical forms and precursor sequences. Through customizing the reference database, CASMIR is applicable to isomiR annotation across species. To demonstrate its application, we investigated isomiR profiles in normal and neoplastic human colorectal epithelia. We also ran miRDeep2, a popular miRNA analysis algorithm to validate isomiRs annotated by CASMIR. With CASMIR, specific and biologically relevant isomiR patterns could be identified. We note that specific isomiRs are often more abundant than their canonical forms. We identify isomiRs that are commonly up-regulated in both colorectal cancer and advanced adenoma, and illustrate advantages in targeting isomiRs as potential biomarkers over canonical forms. CONCLUSIONS: Studying miRNAs at the isomiR level could reveal new insight into miRNA biology and inform assay design for specific isomiRs. CASMIR facilitates comprehensive annotation of isomiR features in small RNA sequencing data for isomiR profiling and differential expression analysis.

Exome-Wide Rare Variant Analyses in Sudden Infant Death Syndrome.

OBJECTIVE: To determine whether a monogenic basis explains sudden infant death syndrome (SIDS) using an exome-wide focus. STUDY DESIGN: A cohort of 427 unrelated cases of SIDS (257 male; average age = 2.7 ± 1.9 months) underwent whole-exome sequencing. Exome-wide rare variant analyses were carried out with 278 SIDS cases of European ancestry (173 male; average age = 2.7 ± 1.98 months) and 973 ethnic-matched controls based on 6 genetic models. Ingenuity Pathway Analysis also was performed. The cohort was collected in collaboration with coroners, medical examiners, and pathologists by St George's University of London, United Kingdom, and Mayo Clinic, Rochester, Minnesota. Whole-exome sequencing was performed at the Genomic Laboratory, Kings College London, United Kingdom, or Mayo Clinic's Medical Genome Facility, Rochester, Minnesota. RESULTS: Although no exome-wide significant (P < 2.5 × 10-6) difference in burden of ultra-rare variants was detected for any gene, 405 genes had a greater prevalence (P < .05) of ultra-rare nonsynonymous variants among cases with 17 genes at P < .005. Some of these potentially overrepresented genes may represent biologically plausible novel candidate genes for a monogenic basis for a portion of patients with SIDS. The top canonical pathway identified was glucocorticoid biosynthesis (P = .01). CONCLUSIONS: The lack of exome-wide significant genetic associations indicates an extreme heterogeneity of etiologies underlying SIDS. Our approach to understanding the genetic mechanisms of SIDS has far reaching implications for the SIDS research community as a whole and may catalyze new evidence-based SIDS research across multiple disciplines. Perturbations in glucocorticoid biosynthesis may represent a novel SIDS-associated biological pathway for future SIDS investigative research.

Enhancer of Zeste Homolog 2 Inhibition Stimulates Bone Formation and Mitigates Bone Loss Caused by Ovariectomy in Skeletally Mature Mice.

Perturbations in skeletal development and bone degeneration may result in reduced bone mass and quality, leading to greater fracture risk. Bone loss is mitigated by bone protective therapies, but there is a clinical need for new bone-anabolic agents. Previous work has demonstrated that Ezh2 (enhancer of zeste homolog 2), a histone 3 lysine 27 (H3K27) methyltransferase, suppressed differentiation of osteogenic progenitors. Here, we investigated whether inhibition of Ezh2 can be leveraged for bone stimulatory applications. Pharmacologic inhibition and siRNA knockdown of Ezh2 enhanced osteogenic commitment of MC3T3 preosteoblasts. Next generation RNA sequencing of mRNAs and real time quantitative PCR profiling established that Ezh2 inactivation promotes expression of bone-related gene regulators and extracellular matrix proteins. Mechanistically, enhanced gene expression was linked to decreased H3K27 trimethylation (H3K27me3) near transcriptional start sites in genome-wide sequencing of chromatin immunoprecipitations assays. Administration of an Ezh2 inhibitor modestly increases bone density parameters of adult mice. Furthermore, Ezh2 inhibition also alleviated bone loss in an estrogen-deficient mammalian model for osteoporosis. Ezh2 inhibition enhanced expression of Wnt10b and Pth1r and increased the BMP-dependent phosphorylation of Smad1/5. Thus, these data suggest that inhibition of Ezh2 promotes paracrine signaling in osteoblasts and has bone-anabolic and osteoprotective potential in adults.

Homozygous/Compound Heterozygous Triadin Mutations Associated With Autosomal-Recessive Long-QT Syndrome and Pediatric Sudden Cardiac Arrest: Elucidation of the Triadin Knockout Syndrome.

BACKGROUND: Long-QT syndrome (LQTS) may result in syncope, seizures, or sudden cardiac arrest. Although 16 LQTS-susceptibility genes have been discovered, 20% to 25% of LQTS remains genetically elusive. METHODS AND RESULTS: We performed whole-exome sequencing child-parent trio analysis followed by recessive and sporadic inheritance modeling and disease-network candidate analysis gene ranking to identify a novel underlying genetic mechanism for LQTS. Subsequent mutational analysis of the candidate gene was performed with polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing on a cohort of 33 additional unrelated patients with genetically elusive LQTS. After whole-exome sequencing and variant filtration, a homozygous p.D18fs*13 TRDN-encoded triadin frameshift mutation was discovered in a 10-year-old female patient with LQTS with a QTc of 500 milliseconds who experienced recurrent exertion-induced syncope/cardiac arrest beginning at 1 year of age. Subsequent mutational analysis of TRDN revealed either homozygous or compound heterozygous frameshift mutations in 4 of 33 unrelated cases of LQTS (12%). All 5 TRDN-null patients displayed extensive T-wave inversions in precordial leads V1 through V4, with either persistent or transient QT prolongation and severe disease expression of exercise-induced cardiac arrest in early childhood (≤3 years of age) and required aggressive therapy. The overall yield of TRDN mutations was significantly greater in patients ≤10 years of age (5 of 10, 50%) compared with older patients (0 of 24, 0%; P=0.0009). CONCLUSIONS: We identified TRDN as a novel underlying genetic basis for recessively inherited LQTS. All TRDN-null patients had strikingly similar phenotypes. Given the recurrent nature of potential lethal arrhythmias, patients fitting this phenotypic profile should undergo cardiac TRDN genetic testing.

TNNI3K mutation in familial syndrome of conduction system disease, atrial tachyarrhythmia and dilated cardiomyopathy.

Locus mapping has uncovered diverse etiologies for familial atrial fibrillation (AF), dilated cardiomyopathy (DCM), and mixed cardiac phenotype syndromes, yet the molecular basis for these disorders remains idiopathic in most cases. Whole-exome sequencing (WES) provides a powerful new tool for familial disease gene discovery. Here, synergistic application of these genomic strategies identified the pathogenic mutation in a familial syndrome of atrial tachyarrhythmia, conduction system disease (CSD), and DCM vulnerability. Seven members of a three-generation family exhibited the variably expressed phenotype, three of whom manifested CSD and clinically significant arrhythmia in childhood. Genome-wide linkage analysis mapped two equally plausible loci to chromosomes 1p3 and 13q12. Variants from WES of two affected cousins were filtered for rare, predicted-deleterious, positional variants, revealing an unreported heterozygous missense mutation disrupting the highly conserved kinase domain in TNNI3K. The G526D substitution in troponin I interacting kinase, with the most deleterious SIFT and Polyphen2 scores possible, resulted in abnormal peptide aggregation in vitro and in silico docking models predicted altered yet energetically favorable wild-type mutant dimerization. Ventricular tissue from a mutation carrier displayed histopathological hallmarks of DCM and reduced TNNI3K protein staining with unique amorphous nuclear and sarcoplasmic inclusions. In conclusion, mutation of TNNI3K, encoding a heart-specific kinase previously shown to modulate cardiac conduction and myocardial function in mice, underlies a familial syndrome of electrical and myopathic heart disease. The identified substitution causes a TNNI3K aggregation defect and protein deficiency, implicating a dominant-negative loss of function disease mechanism.

De novo RRAGC mutation activates mTORC1 signaling in syndromic fetal dilated cardiomyopathy.

Idiopathic dilated cardiomyopathy (DCM) is a heritable, genetically heterogeneous disorder with variable age-dependent penetrance. We sought to identify the genetic underpinnings of syndromic, sporadic DCM in a newborn female diagnosed in utero. Postnatal evaluation revealed ventricular dilation and systolic dysfunction, bilateral cataracts, and mild facial dysmorphisms. Comprehensive metabolic and genetic testing, including chromosomal microarray, mitochondrial DNA and targeted RASopathy gene sequencing, and clinical whole exome sequencing for known cardiomyopathy genes was non-diagnostic. Following exclusion of asymptomatic DCM in the parents, trio-based whole exome sequencing was carried out on a research basis, filtering for rare, predicted deleterious de novo and recessive variants. An unreported de novo S75Y mutation was discovered in RRAGC, encoding Ras-related GTP binding C, an essential GTPase in nutrient-activated mechanistic target of rapamycin complex 1 (mTORC1) signaling. In silico protein modeling and molecular dynamics simulation predicted the mutation to disrupt ligand interactions and increase the GDP-bound state. Overexpression of RagC(S75Y) rendered AD293 cells partially insensitive to amino acid deprivation, resulting in increased mTORC1 signaling compared to wild-type RagC. These findings implicate mTORC1 dysregulation through a gain-of-function mutation in RagC as a novel molecular basis for syndromic forms of pediatric heart failure, and expand genotype-phenotype correlation in RASopathy-related syndromes.

Compound heterozygous NOTCH1 mutations underlie impaired cardiogenesis in a patient with hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) that necessitates staged, single ventricle surgical palliation. An increased frequency of bicuspid aortic valve (BAV) has been observed among relatives. We postulated number of mutant alleles as a molecular basis for variable CHD expression in an extended family comprised of an HLHS proband and four family members who underwent echocardiography and whole-genome sequencing (WGS). Dermal fibroblast-derived induced pluripotent stem cells (iPSC) were procured from the proband-parent trio and bioengineered into cardiomyocytes. Cardiac phenotyping revealed aortic valve atresia and a slit-like left ventricular cavity in the HLHS proband, isolated bicuspid pulmonary valve in his mother, BAV in a maternal 4° relative, and no CHD in his father or sister. Filtering of WGS for rare, functional variants that segregated with CHD and were compound heterozygous in the HLHS proband identified NOTCH1 as the sole candidate gene. An unreported missense mutation (P1964L) in the cytoplasmic domain, segregating with semilunar valve malformation, was maternally inherited and a rare missense mutation (P1256L) in the extracellular domain, clinically silent in the heterozygous state, was paternally inherited. Patient-specific iPSCs exhibited diminished transcript levels of NOTCH1 signaling pathway components, impaired myocardiogenesis, and a higher prevalence of heterogeneous myofilament organization. Extended, phenotypically characterized families enable WGS-derived variant filtering for plausible Mendelian modes of inheritance, a powerful strategy to discover molecular underpinnings of CHD. Identification of compound heterozygous NOTCH1 mutations and iPSC-based functional modeling implicate mutant allele burden and impaired myogenic potential as mechanisms for HLHS.

PatternCNV: a versatile tool for detecting copy number changes from exome sequencing data.

MOTIVATION: Exome sequencing (exome-seq) data, which are typically used for calling exonic mutations, have also been utilized in detecting DNA copy number variations (CNVs). Despite the existence of several CNV detection tools, there is still a great need for a sensitive and an accurate CNV-calling algorithm with built-in QC steps, and does not require a paired reference for each sample. RESULTS: We developed a novel method named PatternCNV, which (i) accounts for the read coverage variations between exons while leveraging the consistencies of this variability across different samples; (ii) reduces alignment BAM files to WIG format and therefore greatly accelerates computation; (iii) incorporates multiple QC measures designed to identify outlier samples and batch effects; and (iv) provides a variety of visualization options including chromosome, gene and exon-level views of CNVs, along with a tabular summarization of the exon-level CNVs. Compared with other CNV-calling algorithms using data from a lymphoma exome-seq study, PatternCNV has higher sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: The software for PatternCNV is implemented using Perl and R, and can be used in Mac or Linux environments. Software and user manual are available at http://bioinformaticstools.mayo.edu/research/patterncnv/, and R package at https://github.com/topsoil/patternCNV/.

Exome sequencing establishes diagnosis of Alström syndrome in an infant presenting with non-syndromic dilated cardiomyopathy.

Idiopathic dilated cardiomyopathy is a heritable, genetically heterogeneous disorder characterized by progressive heart failure. Dilated cardiomyopathy typically exhibits autosomal dominant inheritance, yet frequently remains clinically silent until adulthood. We sought to discover the molecular basis of idiopathic, non-syndromic dilated cardiomyopathy in a one-month-old male presenting with severe heart failure. Previous comprehensive testing of blood, urine, and skin biopsy specimen was negative for metabolic, mitochondrial, storage, and infectious etiologies. Ophthalmologic examination was normal. Chromosomal microarray and commercial dilated cardiomyopathy gene panel testing failed to identify a causative mutation. Parental screening echocardiograms revealed no evidence of clinically silent dilated cardiomyopathy. Whole exome sequencing was carried out on the family trio on a research basis, filtering for rare, deleterious, recessive and de novo genetic variants. Pathogenic compound heterozygous truncating mutations were identified in ALMS1, diagnostic of Alström syndrome and prompting disclosure of genetic findings. Alström syndrome is a known cause for dilated cardiomyopathy in children yet delayed and mis-diagnosis are common owing to its rarity and age-dependent emergence of multisystem clinical manifestations. At six months of age the patient ultimately developed bilateral nystagmus and hyperopia, features characteristic of the syndrome. Early diagnosis is guiding clinical monitoring of other organ systems and allowing for presymptomatic intervention. Furthermore, recognition of recessive inheritance as the mechanism for sporadic disease has informed family planning. This case highlights a limitation of standard gene testing panels for pediatric dilated cardiomyopathy and exemplifies the potential for whole exome sequencing to solve a diagnostic dilemma and enable personalized care.

Association between magnitude of the virus-specific plasmablast response and disease severity in dengue patients.

Dengue is a globally expanding disease caused by infection with dengue virus (DENV) that ranges from febrile illness to acute disease with serious complications. Secondary infection predisposes individuals to more severe disease, and B lymphocytes may play a role in this phenomenon through production of Ab that enhance infection. To better define the acute B cell response during dengue, we analyzed peripheral B cells from an adult Brazilian hospital cohort with primary and secondary DENV infections of varying clinical severity. Circulating B cells in dengue patients were proliferating, activated, and apoptotic relative to individuals with other febrile illnesses. Severe secondary DENV infection was associated with extraordinary peak plasmablast frequencies between 4 and 7 d of illness, averaging 46% and reaching 87% of B cells, significantly greater than those seen in mild illness or primary infections. On average >70% of IgG-secreting cells in individuals with severe secondary DENV infection were DENV specific. Plasmablasts produced Ab that cross-reacted with heterotypic DENV serotypes, but with a 3-fold greater reactivity to DENV-3, the infecting serotype. Plasmablast frequency did not correlate with acute serum-neutralizing Ab titers to any DENV serotype regardless of severity of disease. These findings indicate that massive expansion of DENV-specific and serotype cross-reactive plasmablasts occurs in acute secondary DENV infection of adults in Brazil, which is associated with increasing disease severity.

Integrated analysis of genome-wide DNA methylation and gene expression profiles in molecular subtypes of breast cancer.

Aberrant DNA methylation of CpG islands, CpG island shores and first exons is known to play a key role in the altered gene expression patterns in all human cancers. To date, a systematic study on the effect of DNA methylation on gene expression using high resolution data has not been reported. In this study, we conducted an integrated analysis of MethylCap-sequencing data and Affymetrix gene expression microarray data for 30 breast cancer cell lines representing different breast tumor phenotypes. As well-developed methods for the integrated analysis do not currently exist, we created a series of four different analysis methods. On the computational side, our goal is to develop methylome data analysis protocols for the integrated analysis of DNA methylation and gene expression data on the genome scale. On the cancer biology side, we present comprehensive genome-wide methylome analysis results for differentially methylated regions and their potential effect on gene expression in 30 breast cancer cell lines representing three molecular phenotypes, luminal, basal A and basal B. Our integrated analysis demonstrates that methylation status of different genomic regions may play a key role in establishing transcriptional patterns in molecular subtypes of human breast cancer.

HiChIP: a high-throughput pipeline for integrative analysis of ChIP-Seq data.

BACKGROUND: Chromatin immunoprecipitation (ChIP) followed by next-generation sequencing (ChIP-Seq) has been widely used to identify genomic loci of transcription factor (TF) binding and histone modifications. ChIP-Seq data analysis involves multiple steps from read mapping and peak calling to data integration and interpretation. It remains challenging and time-consuming to process large amounts of ChIP-Seq data derived from different antibodies or experimental designs using the same approach. To address this challenge, there is a need for a comprehensive analysis pipeline with flexible settings to accelerate the utilization of this powerful technology in epigenetics research. RESULTS: We have developed a highly integrative pipeline, termed HiChIP for systematic analysis of ChIP-Seq data. HiChIP incorporates several open source software packages selected based on internal assessments and published comparisons. It also includes a set of tools developed in-house. This workflow enables the analysis of both paired-end and single-end ChIP-Seq reads, with or without replicates for the characterization and annotation of both punctate and diffuse binding sites. The main functionality of HiChIP includes: (a) read quality checking; (b) read mapping and filtering; (c) peak calling and peak consistency analysis; and (d) result visualization. In addition, this pipeline contains modules for generating binding profiles over selected genomic features, de novo motif finding from transcription factor (TF) binding sites and functional annotation of peak associated genes. CONCLUSIONS: HiChIP is a comprehensive analysis pipeline that can be configured to analyze ChIP-Seq data derived from varying antibodies and experiment designs. Using public ChIP-Seq data we demonstrate that HiChIP is a fast and reliable pipeline for processing large amounts of ChIP-Seq data.

Histone deacetylase inhibition promotes osteoblast maturation by altering the histone H4 epigenome and reduces Akt phosphorylation.

Bone has remarkable regenerative capacity, but this ability diminishes during aging. Histone deacetylase inhibitors (HDIs) promote terminal osteoblast differentiation and extracellular matrix production in culture. The epigenetic events altered by HDIs in osteoblasts may hold clues for the development of new anabolic treatments for osteoporosis and other conditions of low bone mass. To assess how HDIs affect the epigenome of committed osteoblasts, MC3T3 cells were treated with suberoylanilide hydroxamic acid (SAHA) and subjected to microarray gene expression profiling and high-throughput ChIP-Seq analysis. As expected, SAHA induced differentiation and matrix calcification of osteoblasts in vitro. ChIP-Seq analysis revealed that SAHA increased histone H4 acetylation genome-wide and in differentially regulated genes, except for the 500 bp upstream of transcriptional start sites. Pathway analysis indicated that SAHA increased the expression of insulin signaling modulators, including Slc9a3r1. SAHA decreased phosphorylation of insulin receptor ß, Akt, and the Akt substrate FoxO1, resulting in FoxO1 stabilization. Thus, SAHA induces genome-wide H4 acetylation and modulates the insulin/Akt/FoxO1 signaling axis, whereas it promotes terminal osteoblast differentiation in vitro.

CAP-miRSeq: a comprehensive analysis pipeline for microRNA sequencing data.

BACKGROUND: miRNAs play a key role in normal physiology and various diseases. miRNA profiling through next generation sequencing (miRNA-seq) has become the main platform for biological research and biomarker discovery. However, analyzing miRNA sequencing data is challenging as it needs significant amount of computational resources and bioinformatics expertise. Several web based analytical tools have been developed but they are limited to processing one or a pair of samples at time and are not suitable for a large scale study. Lack of flexibility and reliability of these web applications are also common issues. RESULTS: We developed a Comprehensive Analysis Pipeline for microRNA Sequencing data (CAP-miRSeq) that integrates read pre-processing, alignment, mature/precursor/novel miRNA detection and quantification, data visualization, variant detection in miRNA coding region, and more flexible differential expression analysis between experimental conditions. According to computational infrastructure, users can install the package locally or deploy it in Amazon Cloud to run samples sequentially or in parallel for a large number of samples for speedy analyses. In either case, summary and expression reports for all samples are generated for easier quality assessment and downstream analyses. Using well characterized data, we demonstrated the pipeline's superior performances, flexibility, and practical use in research and biomarker discovery. CONCLUSIONS: CAP-miRSeq is a powerful and flexible tool for users to process and analyze miRNA-seq data scalable from a few to hundreds of samples. The results are presented in the convenient way for investigators or analysts to conduct further investigation and discovery.