Octopus-toolkit: a workflow to automate mining of public epigenomic and transcriptomic next-generation sequencing data.

Octopus-toolkit is a stand-alone application for retrieving and processing large sets of next-generation sequencing (NGS) data with a single step. Octopus-toolkit is an automated set-up-and-analysis pipeline utilizing the Aspera, SRA Toolkit, FastQC, Trimmomatic, HISAT2, STAR, Samtools, and HOMER applications. All the applications are installed on the user's computer when the program starts. Upon the installation, it can automatically retrieve original files of various epigenomic and transcriptomic data sets, including ChIP-seq, ATAC-seq, DNase-seq, MeDIP-seq, MNase-seq and RNA-seq, from the gene expression omnibus data repository. The downloaded files can then be sequentially processed to generate BAM and BigWig files, which are used for advanced analyses and visualization. Currently, it can process NGS data from popular model genomes such as, human (Homo sapiens), mouse (Mus musculus), dog (Canis lupus familiaris), plant (Arabidopsis thaliana), zebrafish (Danio rerio), fruit fly (Drosophila melanogaster), worm (Caenorhabditis elegans), and budding yeast (Saccharomyces cerevisiae) genomes. With the processed files from Octopus-toolkit, the meta-analysis of various data sets, motif searches for DNA-binding proteins, and the identification of differentially expressed genes and/or protein-binding sites can be easily conducted with few commands by users. Overall, Octopus-toolkit facilitates the systematic and integrative analysis of available epigenomic and transcriptomic NGS big data.

Transposable element-mediated structural variation analysis in dog breeds using whole-genome sequencing.

Naturally occurring diseases in dogs provide an important animal model for studying human disease including cancer, heart disease, and autoimmune disorders. Transposable elements (TEs) make up ~ 31% of the dog (Canis lupus familiaris) genome and are one of main drivers to cause genomic variations and alter gene expression patterns of the host genes, which could result in genetic diseases. To detect structural variations (SVs), we conducted whole-genome sequencing of three different breeds, including Maltese, Poodle, and Yorkshire Terrier. Genomic SVs were detected and visualized using BreakDancer program. We identified a total of 2328 deletion SV events in the three breeds compared with the dog reference genome of Boxer. The majority of the genetic variants were found to be TE insertion polymorphism (1229) and the others were TE-mediated deletion (489), non-TE-mediated deletion (542), simple repeat-mediated deletion (32), and other indel (36). Among the TE insertion polymorphism, 286 elements were full-length LINE-1s (L1s). In addition, the 49 SV candidates located in the genic regions were experimentally verified and their polymorphic rates within each breed were examined using PCR assay. Polymorphism analysis of the genomic variants revealed that some of the variants exist polymorphic in the three dog breeds, suggesting that their SV events recently occurred in the dog genome. The findings suggest that TEs have contributed to the genomic variations among the three dog breeds of Maltese, Poodle, and Yorkshire Terrier. In addition, the polymorphic events between the dog breeds indicate that TEs were recently retrotransposed in the dog genome.

Barrier to autointegration factor 1, procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3, and splicing factor 3b subunit 4 as early-stage cancer decision markers and drivers of hepatocellular carcinoma.

An accurate tool enabling early diagnosis of hepatocellular carcinoma (HCC) is clinically important, given that early detection of HCC markedly improves survival. We aimed to investigate the molecular markers underlying early progression of HCC that can be detected in precancerous lesions. We designed a gene selection strategy to identify potential driver genes by integrative analysis of transcriptome and clinicopathological data of human multistage HCC tissues, including precancerous lesions, low- and high-grade dysplastic nodules. The gene selection process was guided by detecting the selected molecules in both HCC and precancerous lesion. Using various computational approaches, we selected 10 gene elements as a candidate and, through immunohistochemical staining, showed that barrier to autointegration factor 1 (BANF1), procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3 (PLOD3), and splicing factor 3b subunit 4 (SF3B4) are HCC decision markers with superior capability to diagnose early-stage HCC in a large cohort of HCC patients, as compared to the currently popular trio of HCC diagnostic markers: glypican 3, glutamine synthetase, and heat-shock protein 70. Targeted inactivation of BANF1, PLOD3, and SF3B4 inhibits in vitro and in vivo liver tumorigenesis by selectively modulating epithelial-mesenchymal transition and cell-cycle proteins. Treatment of nanoparticles containing small-interfering RNAs of the three genes suppressed liver tumor incidence as well as tumor growth rates in a spontaneous mouse HCC model. We also demonstrated that SF3B4 overexpression triggers SF3b complex to splice tumor suppressor KLF4 transcript to nonfunctional skipped exon transcripts. This contributes to malignant transformation and growth of hepatocyte through transcriptional inactivation of p27Kip1 and simultaneously activation of Slug genes. CONCLUSION: The findings suggest molecular markers of BANF1, PLOD3, and SF3B4 indicating early-stage HCC in precancerous lesion, and also suggest drivers for understanding the development of hepatocarcinogenesis. (Hepatology 2018;67:1360-1377).

The T47D cell line is an ideal experimental model to elucidate the progesterone-specific effects of a luminal A subtype of breast cancer.

Cell lines are often used as in vitro tools to mimic certain types of in vivo system; several cell lines, including MCF-7 and T47D, have been widely used in breast cancer studies without investigating the cell lines' characteristics. In this study, we compared the genome-wide binding profiles of ERα, PR, and P300, and the gene expression changes between MCF-7 and T47D cell lines that represent the luminal A subtype of breast cancer. Surprisingly, several thousand genes were differentially expressed under estrogenic condition. In addition, ERα, PR, and P300 binding to regulatory elements showed distinct genomic landscapes between MCF-7 and T47D cell lines in the same hormonal states. In particular, the T47D cell line was markedly susceptible to progesterone, whereas the MCF-7 cell line did not respond to progesterone in the presence of estrogen. Consistently, changes in the expression level of the PR-target gene, STAT5A, were only observed in the T47D cell line, not the MCF-7 cell line, when treated with progesterone. Overall, the results highlight the importance of selecting appropriate cell lines for breast cancer studies and suggest that T47D cell lines can be an ideal experimental model to elucidate the progesterone-specific effects of a luminal A subtype of breast cancer.

Comparative characterization of 3D chromatin organization in triple-negative breast cancers.

Triple-negative breast cancer (TNBC) is a malignant cancer subtype with a high risk of recurrence and an aggressive phenotype compared to other breast cancer subtypes. Although many breast cancer studies conducted to date have investigated genetic variations and differential target gene expression, how 3D chromatin architectures are reorganized in TNBC has been poorly elucidated. Here, using in situ Hi-C technology, we characterized the 3D chromatin organization in cells representing five distinct subtypes of breast cancer (including TNBC) compared to that in normal cells. We found that the global and local 3D architectures were severely disrupted in breast cancer. TNBC cell lines (especially BT549 cells) showed the most dramatic changes relative to normal cells. Importantly, we detected CTCF-dependent TNBC-susceptible losses/gains of 3D chromatin organization and found that these changes were strongly associated with perturbed chromatin accessibility and transcriptional dysregulation. In TNBC tissue, 3D chromatin disorganization was also observed relative to the 3D chromatin organization in normal tissues. We observed that the perturbed local 3D architectures found in TNBC cells were partially conserved in TNBC tissues. Finally, we discovered distinct tissue-specific chromatin loops by comparing normal and TNBC tissues. In this study, we elucidated the characteristics of the 3D chromatin organization in breast cancer relative to normal cells/tissues at multiple scales and identified associations between disrupted structures and various epigenetic features and transcriptomes. Collectively, our findings reveal important 3D chromatin structural features for future diagnostic and therapeutic studies of TNBC.

The chromatin remodeler Ino80 mediates RNAPII pausing site determination.

BACKGROUND: Promoter-proximal pausing of RNA polymerase II (RNAPII) is a critical step for the precise regulation of gene expression. Despite the apparent close relationship between promoter-proximal pausing and nucleosome, the role of chromatin remodeler governing this step has mainly remained elusive. RESULTS: Here, we report highly confined RNAPII enrichments downstream of the transcriptional start site in Saccharomyces cerevisiae using PRO-seq experiments. This non-uniform distribution of RNAPII exhibits both similar and different characteristics with promoter-proximal pausing in Schizosaccharomyces pombe and metazoans. Interestingly, we find that Ino80p knockdown causes a significant upstream transition of promoter-proximal RNAPII for a subset of genes, relocating RNAPII from the main pausing site to the alternative pausing site. The proper positioning of RNAPII is largely dependent on nucleosome context. We reveal that the alternative pausing site is closely associated with the + 1 nucleosome, and nucleosome architecture around the main pausing site of these genes is highly phased. In addition, Ino80p knockdown results in an increase in fuzziness and a decrease in stability of the + 1 nucleosome. Furthermore, the loss of INO80 also leads to the shift of promoter-proximal RNAPII toward the alternative pausing site in mouse embryonic stem cells. CONCLUSIONS: Based on our collective results, we hypothesize that the highly conserved chromatin remodeler Ino80p is essential in establishing intact RNAPII pausing during early transcription elongation in various organisms, from budding yeast to mouse.

A compendium of chromatin contact maps reflecting regulation by chromatin remodelers in budding yeast.

We herein employ in situ Hi-C with an auxin-inducible degron (AID) system to examine the effect of chromatin remodeling on 3D genome organization in yeast. Eight selected ATP-dependent chromatin remodelers representing various subfamilies contribute to 3D genome organization differently. Among the studied remodelers, the temporary depletions of Chd1p, Swr1p, and Sth1p (a catalytic subunit of the Remodeling the Structure of Chromatin [RSC] complex) cause the most significant defects in intra-chromosomal contacts, and the regulatory roles of these three remodelers in 3D genome organization differ depending on the chromosomal context and cell cycle stage. Furthermore, even though Chd1p and Isw1p are known to share functional similarities/redundancies, their depletions lead to distinct effects on 3D structures. The RSC and cohesin complexes also differentially modulate 3D genome organization within chromosome arm regions, whereas RSC appears to support the function of cohesin in centromeric clustering at G2 phase. Our work suggests that the ATP-dependent chromatin remodelers control the 3D genome organization of yeast through their chromatin-remodeling activities.

ETS1 Suppresses Tumorigenesis of Human Breast Cancer via Trans-Activation of Canonical Tumor Suppressor Genes.

ETS1 has shown dichotomous roles as an oncogene and a tumor suppressor gene in diverse cancers, but its functionality in breast cancer tumorigenesis still remains unclear. We utilized the Cancer Genome Atlas (TCGA) database to analyze comprehensive functions of ETS1 in human breast cancer (BRCA) patients by investigating its expression patterns and methylation status in relation to clinical prognosis. ETS1 expression was significantly diminished by hyper-methylation of the ETS1 promoter region in specimens from BRCA patients compared to a healthy control group. Moreover, ETS1 high BRCA patients showed better prognosis and longer survival compared to ETS1 low BRCA patients. Consistent with clinical evidence, comparative transcriptome analysis combined with CRISPR/Cas9 or shRNA based perturbation of ETS1 expression revealed direct as well as indirect mechanisms of ETS1 that hinder tumorigenesis of BRCA cells. Taken together, our study enlightens a novel function of ETS1 as a tumor suppressor in breast cancer cells.

Ets1 suppresses atopic dermatitis by suppressing pathogenic T cell responses.

Atopic dermatitis (AD) is a complex inflammatory skin disease mediated by immune cells of both adaptive and innate types. Among them, CD4+ Th cells are one of major players of AD pathogenesis. Although the pathogenic role of Th2 cells has been well characterized, Th17/Th22 cells are also implicated in the pathogenesis of AD. However, the molecular mechanisms underlying pathogenic immune responses in AD remain unclear. We sought to investigate how the defect in the AD susceptibility gene, Ets1, is involved in AD pathogenesis in human and mice and its clinical relevance in disease severity by identifying Ets1 target genes and binding partners. Consistent with the decrease in ETS1 levels in severe AD patients and the experimental AD-like skin inflammation model, T cell-specific Ets1-deficient mice (Ets1ΔdLck) developed severe AD-like symptoms with increased pathogenic Th cell responses. A T cell-intrinsic increase of gp130 expression upon Ets1 deficiency promotes the gp130-mediated IL-6 signaling pathway, thereby leading to the development of severe AD-like symptoms. Functional blocking of gp130 by selective inhibitor SC144 ameliorated the disease pathogenesis by reducing pathogenic Th cell responses. Our results reveal a protective role of Ets1 in restricting pathogenic Th cell responses and suggest a potential therapeutic target for AD treatment.

Upregulation of Ets1 expression by NFATc2 and NFKB1/RELA promotes breast cancer cell invasiveness.

Breast cancer is highly aggressive and is the leading cause of cancer-related mortality in women in developed countries. The ETS proto-oncogene 1 (Ets1) has versatile roles during the cellular processes of cancer development. It is often highly expressed in breast cancers and mediates migration and invasion of human breast cancer cells. However, underlying mechanisms of Ets1 gene expression is still ambiguous. Here, we identified a core-regulatory element (CRE) located in the Ets1 promoter region (-540/-80 bp from TSS) that contains elements responsible for associating with NFATs and NF-κBs. Compared with the less metastatic breast cancer cells, metastatic breast cancer cells (MDA-MB-231) show open chromatin configurations in the CRE, which facilitates direct binding of NFATc2 and/or NFKB1/RELA complex to trans-activate Ets1 transcription. Moreover, enhanced level of Nfatc2 and Nfkb1 positively correlated with Ets1 expression in the human breast cancer specimens. Deletion of the CRE region by CRISPR/Cas9 system resulted in significant reduction in Ets1 expression, which led to alterations of Ets1-mediated transcription programs including tumor invasiveness-related genes. Proper regulation of Ets1 gene expression by targeting the NFATc2 and NFKB1/RELA interaction could be a potential therapeutic target for Ets1-mediated metastatic breast cancer.

Janus Kinase 1 Is Essential for Inflammatory Cytokine Signaling and Mammary Gland Remodeling.

Despite a wealth of knowledge about the significance of individual signal transducers and activators of transcription (STATs), essential functions of their upstream Janus kinases (JAKs) during postnatal development are less well defined. Using a novel mammary gland-specific JAK1 knockout model, we demonstrate here that this tyrosine kinase is essential for the activation of STAT1, STAT3, and STAT6 in the mammary epithelium. The loss of JAK1 uncouples interleukin-6-class ligands from their downstream effector, STAT3, which leads to the decreased expression of STAT3 target genes that are associated with the acute-phase response, inflammation, and wound healing. Consequently, JAK1-deficient mice exhibit impaired apoptosis and a significant delay in mammary gland remodeling. Using RNA sequencing, we identified several new JAK1 target genes that are upregulated during involution. These include Bmf and Bim, which are known regulators of programmed cell death. Using a BMF/BIM-double-knockout epithelial transplant model, we further validated that the synergistic action of these proapoptotic JAK1 targets is obligatory for the remodeling of the mammary epithelium. The collective results of this study suggest that JAK1 has nonredundant roles in the activation of particular STAT proteins and this tyrosine kinase is essential for coupling inflammatory cytokine signals to the cell death machinery in the differentiated mammary epithelium.