Rendong Yang

Therapeutic targeting of BET bromodomain proteins in castration-resistant prostate cancer

Men who develop metastatic castration-resistant prostate cancer (CRPC) invariably succumb to the disease. Progression to CRPC after androgen ablation therapy is predominantly driven by deregulated androgen receptor (AR) signalling. Despite the success of recently approved therapies targeting AR signalling, such as abiraterone and second-generation anti-androgens including MDV3100 (also known as enzalutamide), durable responses are limited, presumably owing to acquired resistance. Recently, JQ1 and I-BET762 two selective small-molecule inhibitors that target the amino-terminal bromodomains of BRD4, have been shown to exhibit anti-proliferative effects in a range of malignancies. Here we show that AR-signalling-competent human CRPC cell lines are preferentially sensitive to bromodomain and extraterminal (BET) inhibition. BRD4 physically interacts with the N-terminal domain of AR and can be disrupted by JQ1 (refs 11, 13). Like the direct AR antagonist MDV3100, JQ1 disrupted AR recruitment to target gene loci. By contrast with MDV3100, JQ1 functions downstream of AR, and more potently abrogated BRD4 localization to AR target loci and AR-mediated gene transcription, including induction of the TMPRSS2-ERG gene fusion and its oncogenic activity. In vivo, BET bromodomain inhibition was more efficacious than direct AR antagonism in CRPC xenograft mouse models. Taken together, these studies provide a novel epigenetic approach for the concerted blockade of oncogenic drivers in advanced prostate cancer.

Analyzing circadian expression data by harmonic regression based on autoregressive spectral estimation

Motivation: Circadian rhythms are prevalent in most organisms. Identification of circadian-regulated genes is a crucial step in discovering underlying pathways and processes that are clock-controlled. Such genes are largely detected by searching periodic patterns in microarray data. However, temporal gene expression profiles usually have a short time-series with low sampling frequency and high levels of noise. This makes circadian rhythmic analysis of temporal microarray data very challenging. Results: We propose an algorithm named ARSER, which combines time domain and frequency domain analysis for extracting and characterizing rhythmic expression profiles from temporal microarray data. ARSER employs autoregressive spectral estimation to predict an expression profiles periodicity from the frequency spectrum and then models the rhythmic patterns by using a harmonic regression model to fit the time-series. ARSER describes the rhythmic patterns by four parameters: period, phase, amplitude and mean level, and measures the multiple testing significance by false discovery rate q-value. When tested on well defined periodic and non-periodic short time-series data, ARSER was superior to two existing and widely-used methods, COSOPT and Fishers G-test, during identification of sinusoidal and non-sinusoidal periodic patterns in short, noisy and non-stationary time-series. Finally, analysis of Arabidopsis microarray data using ARSER led to identification of a novel set of previously undetected non-sinusoidal periodic transcripts, which may lead to new insights into molecular mechanisms of circadian rhythms. Availability: ARSER is implemented by Python and R. All source codes are available from http://bioinformatics.cau.edu.cn/ARSER Contact: zhensu@cau.edu.cn

The central role of EED in the orchestration of polycomb group complexes

Polycomb Repressive Complexes 1 and 2 (PRC1 and 2) play a critical role in the epigenetic regulation of transcription during cellular differentiation, stem cell pluripotency, and neoplastic progression. Here we show that the Polycomb Group protein EED, a core component of PRC2, physically interacts with and functions as part of PRC1. Components of PRC1 and PRC2 compete for EED binding. EED functions to recruit PRC1 to H3K27me3 loci and enhances PRC1 mediated H2A ubiquitin E3 ligase activity. Taken together, we suggest an integral role for EED as an epigenetic exchange factor coordinating the activities of PRC1 and 2.

Truncation and constitutive activation of the androgen receptor by diverse genomic rearrangements in prostate cancer

Molecularly targeted therapies for advanced prostate cancer include castration modalities that suppress ligand-dependent transcriptional activity of the androgen receptor (AR). However, persistent AR signalling undermines therapeutic efficacy and promotes progression to lethal castration-resistant prostate cancer (CRPC), even when patients are treated with potent second-generation AR-targeted therapies abiraterone and enzalutamide. Here we define diverse AR genomic structural rearrangements (AR-GSRs) as a class of molecular alterations occurring in one third of CRPC-stage tumours. AR-GSRs occur in the context of copy-neutral and amplified AR and display heterogeneity in breakpoint location, rearrangement class and sub-clonal enrichment in tumours within and between patients. Despite this heterogeneity, one common outcome in tumours with high sub-clonal enrichment of AR-GSRs is outlier expression of diverse AR variant species lacking the ligand-binding domain and possessing ligand-independent transcriptional activity. Collectively, these findings reveal AR-GSRs as important drivers of persistent AR signalling in CRPC.

Impaired replication elongation in Tetrahymena mutants deficient in histone H3 Lys 27 monomethylation

Replication of nuclear DNA occurs in the context of chromatin and is influenced by histone modifications. In the ciliate Tetrahymena thermophila, we identified TXR1, encoding a histone methyltransferase. TXR1 deletion resulted in severe DNA replication stress, manifested by the accumulation of ssDNA, production of aberrant replication intermediates, and activation of robust DNA damage responses. Paired-end Illumina sequencing of ssDNA revealed intergenic regions, including replication origins, as hot spots for replication stress in ΔTXR1 cells. ΔTXR1 cells showed a deficiency in histone H3 Lys 27 monomethylation (H3K27me1), while ΔEZL2 cells, deleting a Drosophila E(z) homolog, were deficient in H3K27 di- and trimethylation, with no detectable replication stress. A point mutation in histone H3 at Lys 27 (H3 K27Q) mirrored the phenotype of ΔTXR1, corroborating H3K27me1 as a key player in DNA replication. Additionally, we demonstrated interactions between TXR1 and proliferating cell nuclear antigen (PCNA). These findings support a conserved pathway through which H3K27me1 facilitates replication elongation.

Cerebellar Transcriptome Profiles of ATXN1 Transgenic Mice Reveal SCA1 Disease Progression and Protection Pathways

SCA1, a fatal neurodegenerative disorder, is caused by a CAG expansion encoding a polyglutamine stretch in the protein ATXN1. We used RNA sequencing to profile cerebellar gene expression in Pcp2-ATXN1[82Q] mice with ataxia and progressive pathology and Pcp2-ATXN1[30Q]D776 animals having ataxia in absence of Purkinje cell progressive pathology. Weighted Gene Coexpression Network Analysis of the cerebellar expression data revealed two gene networks that significantly correlated with disease and have an expression profile correlating with disease progression in ATXN1[82Q] Purkinje cells. The Magenta Module provides a signature of suppressed transcriptional programs reflecting disease progression in Purkinje cells, while the Lt Yellow Module reflects transcriptional programs activated in response to disease in Purkinje cells as well as other cerebellar cell types. Furthermore, we found that upregulation of cholecystokinin (Cck) and subsequent interaction with the Cck1 receptor likely underlies the lack of progressive Purkinje cell pathology in Pcp2-ATXN1[30Q]D776 mice.

Lineage-Specific Effector Signatures of Invariant NKT Cells Are Shared amongst γδ T, Innate Lymphoid, and Th Cells

Invariant NKT cells differentiate into three predominant effector lineages in the steady state. To understand these lineages, we sorted undifferentiated invariant NK T progenitor cells and each effector population and analyzed their transcriptional profiles by RNAseq. Bioinformatic comparisons were made to effector subsets among other lymphocytes, specifically Th cells, innate lymphoid cells (ILC), and γδ T cells. Myc-associated signature genes were enriched in NKT progenitors, like in other hematopoietic progenitors. Only NKT1 cells, but not NKT2 and NKT17 cells, had transcriptome similarity to NK cells and were also similar to other IFN-γ–producing lineages such as Th1, ILC1, and intraepithelial γδ T cells. NKT2 and NKT17 cells were similar to their analogous subsets of γδ T cells and ILCs, but surprisingly, not to Th2 and Th17 cells. We identified a set of genes common to each effector lineage regardless of Ag receptor specificity, suggesting the use of conserved regulatory cores for effector function.

Androgen Receptor Variant AR-V9 Is Coexpressed with AR-V7 in Prostate Cancer Metastases and Predicts Abiraterone Resistance

Purpose: Androgen receptor (AR) variant AR-V7 is a ligand-independent transcription factor that promotes prostate cancer resistance to AR-targeted therapies. Accordingly, efforts are under way to develop strategies for monitoring and inhibiting AR-V7 in castration-resistant prostate cancer (CRPC). The purpose of this study was to understand whether other AR variants may be coexpressed with AR-V7 and promote resistance to AR-targeted therapies. Experimental Design: We utilized complementary short- and long-read sequencing of intact AR mRNA isoforms to characterize AR expression in CRPC models. Coexpression of AR-V7 and AR-V9 mRNA in CRPC metastases and circulating tumor cells was assessed by RNA-seq and RT-PCR, respectively. Expression of AR-V9 protein in CRPC models was evaluated with polyclonal antisera. Multivariate analysis was performed to test whether AR variant mRNA expression in metastatic tissues was associated with a 12-week progression-free survival endpoint in a prospective clinical trial of 78 CRPC-stage patients initiating therapy with the androgen synthesis inhibitor, abiraterone acetate. Results: AR-V9 was frequently coexpressed with AR-V7. Both AR variant species were found to share a common 3′ terminal cryptic exon, which rendered AR-V9 susceptible to experimental manipulations that were previously thought to target AR-V7 uniquely. AR-V9 promoted ligand-independent growth of prostate cancer cells. High AR-V9 mRNA expression in CRPC metastases was predictive of primary resistance to abiraterone acetate (HR = 4.0; 95% confidence interval, 1.31–12.2; P = 0.02). Conclusions: AR-V9 may be an important component of therapeutic resistance in CRPC. Clin Cancer Res; 23(16); 4704–15. ©2017 AACR.

ScanIndel: a hybrid framework for indel detection via gapped alignment, split reads and de novo assembly

Comprehensive identification of insertions/deletions (indels) across the full size spectrum from second generation sequencing is challenging due to the relatively short read length inherent in the technology. Different indel calling methods exist but are limited in detection to specific sizes with varying accuracy and resolution. We present ScanIndel, an integrated framework for detecting indels with multiple heuristics including gapped alignment, split reads and de novo assembly. Using simulation data, we demonstrate ScanIndel’s superior sensitivity and specificity relative to several state-of-the-art indel callers across various coverage levels and indel sizes. ScanIndel yields higher predictive accuracy with lower computational cost compared with existing tools for both targeted resequencing data from tumor specimens and high coverage whole-genome sequencing data from the human NIST standard NA12878. Thus, we anticipate ScanIndel will improve indel analysis in both clinical and research settings. ScanIndel is implemented in Python, and is freely available for academic use at https://github.com/cauyrd/ScanIndel.

Transcriptome Phase Distribution Analysis Reveals Diurnal Regulated Biological Processes and Key Pathways in Rice Flag Leaves and Seedling Leaves

Plant diurnal oscillation is a 24-hour period based variation. The correlation between diurnal genes and biological pathways was widely revealed by microarray analysis in different species. Rice (Oryza sativa) is the major food staple for about half of the worlds population. The rice flag leaf is essential in providing photosynthates to the grain filling. However, there is still no comprehensive view about the diurnal transcriptome for rice leaves. In this study, we applied rice microarray to monitor the rhythmically expressed genes in rice seedling and flag leaves. We developed a new computational analysis approach and identified 6,266 (10.96%) diurnal probe sets in seedling leaves, 13,773 (24.08%) diurnal probe sets in flag leaves. About 65% of overall transcription factors were identified as flag leaf preferred. In seedling leaves, the peak of phase distribution was from 2:00am to 4:00am, whereas in flag leaves, the peak was from 8:00pm to 2:00am. The diurnal phase distribution analysis of gene ontology (GO) and cis-element enrichment indicated that, some important processes were waken by the light, such as photosynthesis and abiotic stimulus, while some genes related to the nuclear and ribosome involved processes were active mostly during the switch time of light to dark. The starch and sucrose metabolism pathway genes also showed diurnal phase. We conducted comparison analysis between Arabidopsis and rice leaf transcriptome throughout the diurnal cycle. In summary, our analysis approach is feasible for relatively unbiased identification of diurnal transcripts, efficiently detecting some special periodic patterns with non-sinusoidal periodic patterns. Compared to the rice flag leaves, the gene transcription levels of seedling leaves were relatively limited to the diurnal rhythm. Our comprehensive microarray analysis of seedling and flag leaves of rice provided an overview of the rice diurnal transcriptome and indicated some diurnal regulated biological processes and key functional pathways in rice.