BISCUIT: an efficient, standards-compliant tool suite for simultaneous genetic and epigenetic inference in bulk and single-cell studies.

Data from both bulk and single-cell whole-genome DNA methylation experiments are under-utilized in many ways. This is attributable to inefficient mapping of methylation sequencing reads, routinely discarded genetic information, and neglected read-level epigenetic and genetic linkage information. We introduce the BISulfite-seq Command line User Interface Toolkit (BISCUIT) and its companion R/Bioconductor package, biscuiteer, for simultaneous extraction of genetic and epigenetic information from bulk and single-cell DNA methylation sequencing. BISCUIT's performance, flexibility and standards-compliant output allow large, complex experimental designs to be characterized on clinical timescales. BISCUIT is particularly suited for processing data from single-cell DNA methylation assays, with its excellent scalability, efficiency, and ability to greatly enhance mappability, a key challenge for single-cell studies. We also introduce the epiBED format for single-molecule analysis of coupled epigenetic and genetic information, facilitating the study of cellular and tissue heterogeneity from DNA methylation sequencing.

A molecular switch for neuroprotective astrocyte reactivity.

The intrinsic mechanisms that regulate neurotoxic versus neuroprotective astrocyte phenotypes and their effects on central nervous system degeneration and repair remain poorly understood. Here we show that injured white matter astrocytes differentiate into two distinct C3-positive and C3-negative reactive populations, previously simplified as neurotoxic (A1) and neuroprotective (A2)1,2, which can be further subdivided into unique subpopulations defined by proliferation and differential gene expression signatures. We find the balance of neurotoxic versus neuroprotective astrocytes is regulated by discrete pools of compartmented cyclic adenosine monophosphate derived from soluble adenylyl cyclase and show that proliferating neuroprotective astrocytes inhibit microglial activation and downstream neurotoxic astrocyte differentiation to promote retinal ganglion cell survival. Finally, we report a new, therapeutically tractable viral vector to specifically target optic nerve head astrocytes and show that raising nuclear or depleting cytoplasmic cyclic AMP in reactive astrocytes inhibits deleterious microglial or macrophage cell activation and promotes retinal ganglion cell survival after optic nerve injury. Thus, soluble adenylyl cyclase and compartmented, nuclear- and cytoplasmic-localized cyclic adenosine monophosphate in reactive astrocytes act as a molecular switch for neuroprotective astrocyte reactivity that can be targeted to inhibit microglial activation and neurotoxic astrocyte differentiation to therapeutic effect. These data expand on and define new reactive astrocyte subtypes and represent a step towards the development of gliotherapeutics for the treatment of glaucoma and other optic neuropathies.

Effects of biomechanical and biochemical stimuli on angio- and vasculogenesis in a complex microvasculature-on-chip.

The endothelium of blood vessels is a vital organ that reacts differently to subtle changes in stiffness and mechanical forces exerted on its environment (extracellular matrix (ECM)). Upon alteration of these biomechanical cues, endothelial cells initiate signaling pathways that govern vascular remodeling. The emerging organs-on-chip technologies allow the mimicking of complex microvasculature networks, identifying the combined or singular effects of these biomechanical or biochemical stimuli. Here, we present a microvasculature-on-chip model to investigate the singular effect of ECM stiffness and mechanical cyclic stretch on vascular development. Following two different approaches for vascular growth, the effect of ECM stiffness on sprouting angiogenesis and the effect of cyclic stretch on endothelial vasculogenesis are studied. Our results indicate that ECM hydrogel stiffness controls the size of the patterned vasculature and the density of sprouting angiogenesis. RNA sequencing shows that the cellular response to stretching is characterized by the upregulation of certain genes such as ANGPTL4+5, PDE1A, and PLEC.

Mithramycin induces promoter reprogramming and differentiation of rhabdoid tumor.

Rhabdoid tumor (RT) is a pediatric cancer characterized by the inactivation of SMARCB1, a subunit of the SWI/SNF chromatin remodeling complex. Although this deletion is the known oncogenic driver, there are limited effective therapeutic options for these patients. Here we use unbiased screening of cell line panels to identify a heightened sensitivity of rhabdoid tumor to mithramycin and the second-generation analogue EC8042. The sensitivity of MMA and EC8042 was superior to traditional DNA damaging agents and linked to the causative mutation of the tumor, SMARCB1 deletion. Mithramycin blocks SMARCB1-deficient SWI/SNF activity and displaces the complex from chromatin to cause an increase in H3K27me3. This triggers chromatin remodeling and enrichment of H3K27ac at chromHMM-defined promoters to restore cellular differentiation. These effects occurred at concentrations not associated with DNA damage and were not due to global chromatin remodeling or widespread gene expression changes. Importantly, a single 3-day infusion of EC8042 caused dramatic regressions of RT xenografts, recapitulated the increase in H3K27me3, and cellular differentiation described in vitro to completely cure three out of eight mice.

Practical Calling Approach for Exome Array-Based Genome-Wide Association Studies in Korean Population.

Exome-based genotyping arrays are cost-effective and have recently been used as alternative platforms to whole-exome sequencing. However, the automated clustering algorithm in an exome array has a genotype calling problem in accuracy for identifying rare and low-frequency variants. To address these shortcomings, we present a practical approach for accurate genotype calling using the Illumina Infinium HumanExome BeadChip. We present comparison results and a statistical summary of our genotype data sets. Our data set comprises 14,647 Korean samples. To solve the limitation of automated clustering, we performed manual genotype clustering for the targeted identification of 46,076 variants that were identified using GenomeStudio software. To evaluate the effects of applying custom cluster files, we tested cluster files using 804 independent Korean samples and the same platform. Our study firstly suggests practical guidelines for exome chip quality control in Asian populations and provides valuable insight into an association study using exome chip.

Combinatorial approach to estimate copy number genotype using whole-exome sequencing data.

Copy number variations (CNVs) are known risk factors in complex diseases. Array-based approaches have been widely used to detect CNVs, but limitations of array-based CNV detection methods, such as noisy signal and low resolution, have hindered detection of small CNVs. Recently, the development of next-generation sequencing techniques has increased rapidly owing to declines in cost. Particularly, whole-exome sequencing has proved useful for finding causal genes and variants in complex diseases. Because gene copy number may affect expression, CNV genotyping can be very valuable in disease association studies. However, almost all current CNV detection tools consider only two types of CNV genotypes. In this study, we propose a CNV genotype estimation approach using a combination of existing methods. Our approach was comprehensively compared with the customized Agilent array-comparative genomic hybridization. We found that our genotyping approach proved to be accurate, and reproducible, suggesting that it can complement existing CNV genotyping methods.

Recapitulation of previous genome-wide association studies with two distinct pathophysiological entities of gastric cancer in the Korean population.

Gastric cancer (GC) is the most common malignancy. The incidence rates remain remarkably high in East Asians. Although genome-wide association studies in the Han Chinese and Japanese populations have so far yielded susceptibility loci for GC, these findings need to be validated in an independent ethnic group. To identify the potential heterogeneity by histological classified subtypes (intestinal and diffuse), we examined the previously reported associations in the Korean population. PRKAA1 at 5p13.1 was found to be more strongly associated with intestinal type (odds ratio, OR=1.39, 95% CI (confidence interval) =1.22-1.58, P=3.77 × 10(-7)) than diffuse type. In addition, PSCA at 8q23.3 was significantly replicated in diffuse type (OR=1.49, 95% CI=1.32-1.67, P=2.43 × 10(-11)) but far less significant in intestinal type. In conclusion, these findings could bring additional insights into the etiologic heterogeneity in gastric carcinogenesis mechanisms.

A genome-wide association study of a coronary artery disease risk variant.

Although over 30 common genetic susceptibility loci have been identified to be independently associated with coronary artery disease (CAD) risk through genome-wide association studies (GWAS), genetic risk variants reported to date explain only a small fraction of heritability. To identify novel susceptibility variants for CAD and confirm those previously identified in European population, GWAS and a replication study were performed in the Koreans and Japanese. In the discovery stage, we genotyped 2123 cases and 3591 controls with 521 786 SNPs using the Affymetrix SNP Array 6.0 chips in Korean. In the replication, direct genotyping was performed using 3052 cases and 4976 controls from the KItaNagoya Genome study of Japan with 14 selected SNPs. To maximize the coverage of the genome, imputation was performed based on 1000 Genome JPT+CHB and 5.1 million SNPs were retained. CAD association was replicated for three GWAS-identified loci (1p13.3/SORT1 (rs599839), 9p21.3/CDKN2A/2B (rs4977574), and 11q22.3/ PDGFD (rs974819)) in Koreans. From GWAS and a replication, SNP rs3782889 showed a strong association (combined P=3.95 × 10(-14)), although the association of SNP rs3782889 doesn't remain statistically significant after adjusting for SNP rs11066015 (proxy SNP with BRAP (r(2)=1)). But new possible CAD-associated variant was observed for rs9508025 (FLT1), even though its statistical significance did marginally reach at the genome-wide a significance level (combined P=6.07 × 10(-7)). This study shows that three CAD susceptibility loci, which were previously identified in European can be directly replicated in Koreans and also provides additional evidences implicating suggestive loci as risk variants for CAD in East Asian.