Alzheimer's disease rewires gene coexpression networks coupling different brain regions.

Connectome studies have shown how Alzheimer's disease (AD) disrupts functional and structural connectivity among brain regions. But the molecular basis of such disruptions is less studied, with most genomic/transcriptomic studies performing within-brain-region analyses. To inspect how AD rewires the correlation structure among genes in different brain regions, we performed an Inter-brain-region Differential Correlation (Inter-DC) analysis of RNA-seq data from Mount Sinai Brain Bank on four brain regions (frontal pole, superior temporal gyrus, parahippocampal gyrus and inferior frontal gyrus, comprising 264 AD and 372 control human post-mortem samples). An Inter-DC network was assembled from all pairs of genes across two brain regions that gained (or lost) correlation strength in the AD group relative to controls at FDR 1%. The differentially correlated (DC) genes in this network complemented known differentially expressed genes in AD, and likely reflects cell-intrinsic changes since we adjusted for cell compositional effects. Each brain region used a distinctive set of DC genes when coupling with other regions, with parahippocampal gyrus showing the most rewiring, consistent with its known vulnerability to AD. The Inter-DC network revealed master dysregulation hubs in AD (at genes ZKSCAN1, SLC5A3, RCC1, IL17RB, PLK4, etc.), inter-region gene modules enriched for known AD pathways (synaptic signaling, endocytosis, etc.), and candidate signaling molecules that could mediate region-region communication. The Inter-DC network generated in this study is a valuable resource of gene pairs, pathways and signaling molecules whose inter-brain-region functional coupling is disrupted in AD, thereby offering a new perspective of AD etiology.

The genome sequence of Molossusnigricans (Chiroptera, Molossidae; Miller, 1902).

We present a genome assembly from an individual male Molossus nigricans (Chordata; Mammalia; Chiroptera; Molossidae). The genome sequence is 2.41 gigabases in span. The majority of the assembly is scaffolded into 24 chromosomal pseudomolecules, with the X sex chromosome assembled.

iCOMIC: a graphical interface-driven bioinformatics pipeline for analyzing cancer omics data.

Despite the tremendous increase in omics data generated by modern sequencing technologies, their analysis can be tricky and often requires substantial expertise in bioinformatics. To address this concern, we have developed a user-friendly pipeline to analyze (cancer) genomic data that takes in raw sequencing data (FASTQ format) as input and outputs insightful statistics. Our iCOMIC toolkit pipeline featuring many independent workflows is embedded in the popular Snakemake workflow management system. It can analyze whole-genome and transcriptome data and is characterized by a user-friendly GUI that offers several advantages, including minimal execution steps and eliminating the need for complex command-line arguments. Notably, we have integrated algorithms developed in-house to predict pathogenicity among cancer-causing mutations and differentiate between tumor suppressor genes and oncogenes from somatic mutation data. We benchmarked our tool against Genome In A Bottle benchmark dataset (NA12878) and got the highest F1 score of 0.971 and 0.988 for indels and SNPs, respectively, using the BWA MEM-GATK HC DNA-Seq pipeline. Similarly, we achieved a correlation coefficient of r = 0.85 using the HISAT2-StringTie-ballgown and STAR-StringTie-ballgown RNA-Seq pipelines on the human monocyte dataset (SRP082682). Overall, our tool enables easy analyses of omics datasets, significantly ameliorating complex data analysis pipelines.

CBP Regulates Recruitment and Release of Promoter-Proximal RNA Polymerase II.

Transcription activation involves RNA polymerase II (Pol II) recruitment and release from the promoter into productive elongation, but how specific chromatin regulators control these steps is unclear. Here, we identify a novel activity of the histone acetyltransferase p300/CREB-binding protein (CBP) in regulating promoter-proximal paused Pol II. We find that Drosophila CBP inhibition results in "dribbling" of Pol II from the pause site to positions further downstream but impedes transcription through the +1 nucleosome genome-wide. Promoters strongly occupied by CBP and GAGA factor have high levels of paused Pol II, a unique chromatin signature, and are highly expressed regardless of cell type. Interestingly, CBP activity is rate limiting for Pol II recruitment to these highly paused promoters through an interaction with TFIIB but for transit into elongation by histone acetylation at other genes. Thus, CBP directly stimulates both Pol II recruitment and the ability to traverse the first nucleosome, thereby promoting transcription of most genes.

Erratum to: CBP binding outside of promoters and enhancers in Drosophila melanogaster.

[This corrects the article DOI: 10.1186/s13072-015-0042-4.].

CBP binding outside of promoters and enhancers in Drosophila melanogaster.

BACKGROUND: CREB-binding protein (CBP, also known as nejire) is a transcriptional co-activator that is conserved in metazoans. CBP plays an important role in embryonic development and cell differentiation and mutations in CBP can lead to various diseases in humans. In addition, CBP and the related p300 protein have successfully been used to predict enhancers in both humans and flies when they occur with monomethylation of histone H3 on lysine 4 (H3K4me1). RESULTS: Here, we compare CBP chromatin immunoprecipitation sequencing data from Drosophila S2 cells with modENCODE data and show that CBP is bound at genomic sites with a wide range of functions. As expected, we find that CBP is bound at active promoters and enhancers. In addition, we find that the strongest CBP sites in the genome are found at Polycomb response elements embedded in histone H3 lysine 27 trimethylated (H3K27me3) chromatin, where they correlate with binding of the Pho repressive complex. Interestingly, we find that CBP also binds to most insulators in the genome. At a subset of these, CBP may regulate insulating activity, measured as the ability to prevent repressive H3K27 methylation from spreading into adjacent chromatin. CONCLUSIONS: We conclude that CBP could be involved in a much wider range of functions than has previously been appreciated, including Polycomb repression and insulator activity. In addition, we discuss the possibility that a common role for CBP at all functional elements may be to regulate interactions between distant chromosomal regions and speculate that CBP is controlling higher order chromatin organization.

Non-coding roX RNAs prevent the binding of the MSL-complex to heterochromatic regions.

Long non-coding RNAs contribute to dosage compensation in both mammals and Drosophila by inducing changes in the chromatin structure of the X-chromosome. In Drosophila melanogaster, roX1 and roX2 are long non-coding RNAs that together with proteins form the male-specific lethal (MSL) complex, which coats the entire male X-chromosome and mediates dosage compensation by increasing its transcriptional output. Studies on polytene chromosomes have demonstrated that when both roX1 and roX2 are absent, the MSL-complex becomes less abundant on the male X-chromosome and is relocated to the chromocenter and the 4th chromosome. Here we address the role of roX RNAs in MSL-complex targeting and the evolution of dosage compensation in Drosophila. We performed ChIP-seq experiments which showed that MSL-complex recruitment to high affinity sites (HAS) on the X-chromosome is independent of roX and that the HAS sequence motif is conserved in D. simulans. Additionally, a complete and enzymatically active MSL-complex is recruited to six specific genes on the 4th chromosome. Interestingly, our sequence analysis showed that in the absence of roX RNAs, the MSL-complex has an affinity for regions enriched in Hoppel transposable elements and repeats in general. We hypothesize that roX mutants reveal the ancient targeting of the MSL-complex and propose that the role of roX RNAs is to prevent the binding of the MSL-complex to heterochromatin.

Male X-linked genes in Drosophila melanogaster are compensated independently of the Male-Specific Lethal complex.

BACKGROUND: In organisms where the two sexes have unequal numbers of X-chromosomes, the expression of X-linked genes needs to be balanced not only between the two sexes, but also between X and the autosomes. In Drosophila melanogaster, the Male-Specific Lethal (MSL) complex is believed to produce a 2-fold increase in expression of genes on the male X, thus restoring this balance. RESULTS: Here we show that almost all the genes on the male X are effectively compensated. However, many genes are compensated without any significant recruitment of the MSL-complex. These genes are very weakly, if at all, affected by mutations or RNAi against MSL-complex components. In addition, even the genes that are strongly bound by MSL rely on mechanisms other than the MSL-complex for proper compensation. We find that long, non-ubiquitously expressed genes tend to rely less on the MSL-complex for their compensation and genes that in addition are far from High Affinity Sites tend to not bind the complex at all or very weakly. CONCLUSIONS: We conclude that most of the compensation of X-linked genes is produced by an MSL-independent mechanism. Similar to the case of the MSL-mediated compensation we do not yet know the mechanism behind the MSL-independent compensation that appears to act preferentially on long genes. Even if we observe similarities, it remains to be seen if the mechanism is related to the buffering that is observed in autosomal aneuploidies.

HP1a recruitment to promoters is independent of H3K9 methylation in Drosophila melanogaster.

Heterochromatin protein 1 (HP1) proteins, recognized readers of the heterochromatin mark methylation of histone H3 lysine 9 (H3K9me), are important regulators of heterochromatin-mediated gene silencing and chromosome structure. In Drosophila melanogaster three histone lysine methyl transferases (HKMTs) are associated with the methylation of H3K9: Su(var)3-9, Setdb1, and G9a. To probe the dependence of HP1a binding on H3K9me, its dependence on these three HKMTs, and the division of labor between the HKMTs, we have examined correlations between HP1a binding and H3K9me patterns in wild type and null mutants of these HKMTs. We show here that Su(var)3-9 controls H3K9me-dependent binding of HP1a in pericentromeric regions, while Setdb1 controls it in cytological region 2L:31 and (together with POF) in chromosome 4. HP1a binds to the promoters and within bodies of active genes in these three regions. More importantly, however, HP1a binding at promoters of active genes is independent of H3K9me and POF. Rather, it is associated with heterochromatin protein 2 (HP2) and open chromatin. Our results support a hypothesis in which HP1a nucleates with high affinity independently of H3K9me in promoters of active genes and then spreads via H3K9 methylation and transient looping contacts with those H3K9me target sites.

Preferential genome targeting of the CBP co-activator by Rel and Smad proteins in early Drosophila melanogaster embryos.

CBP and the related p300 protein are widely used transcriptional co-activators in metazoans that interact with multiple transcription factors. Whether CBP/p300 occupies the genome equally with all factors or preferentially binds together with some factors is not known. We therefore compared Drosophila melanogaster CBP (nejire) ChIP-seq peaks with regions bound by 40 different transcription factors in early embryos, and we found high co-occupancy with the Rel-family protein Dorsal. Dorsal is required for CBP occupancy in the embryo, but only at regions where few other factors are present. CBP peaks in mutant embryos lacking nuclear Dorsal are best correlated with TGF-ß/Dpp-signaling and Smad-protein binding. Differences in CBP occupancy in mutant embryos reflect gene expression changes genome-wide, but CBP also occupies some non-expressed genes. The presence of CBP at silent genes does not result in histone acetylation. We find that Polycomb-repressed H3K27me3 chromatin does not preclude CBP binding, but restricts histone acetylation at CBP-bound genomic sites. We conclude that CBP occupancy in Drosophila embryos preferentially overlaps factors controlling dorso-ventral patterning and that CBP binds silent genes without causing histone hyperacetylation.

Sequence signatures involved in targeting the Male-Specific Lethal complex to X-chromosomal genes in Drosophila melanogaster.

BACKGROUND: In Drosophila melanogaster, the dosage-compensation system that equalizes X-linked gene expression between males and females, thereby assuring that an appropriate balance is maintained between the expression of genes on the X chromosome(s) and the autosomes, is at least partially mediated by the Male-Specific Lethal (MSL) complex. This complex binds to genes with a preference for exons on the male X chromosome with a 3' bias, and it targets most expressed genes on the X chromosome. However, a number of genes are expressed but not targeted by the complex. High affinity sites seem to be responsible for initial recruitment of the complex to the X chromosome, but the targeting to and within individual genes is poorly understood. RESULTS: We have extensively examined X chromosome sequence variation within five types of gene features (promoters, 5' UTRs, coding sequences, introns, 3' UTRs) and intergenic sequences, and assessed its potential involvement in dosage compensation. Presented results show that: the X chromosome has a distinct sequence composition within its gene features; some of the detected variation correlates with genes targeted by the MSL-complex; the insulator protein BEAF-32 preferentially binds upstream of MSL-bound genes; BEAF-32 and MOF co-localizes in promoters; and that bound genes have a distinct sequence composition that shows a 3' bias within coding sequence. CONCLUSIONS: Although, many strongly bound genes are close to a high affinity site neither our promoter motif nor our coding sequence signatures show any correlation to HAS. Based on the results presented here, we believe that there are sequences in the promoters and coding sequences of targeted genes that have the potential to direct the secondary spreading of the MSL-complex to nearby genes.

Normalization of high dimensional genomics data where the distribution of the altered variables is skewed.

Genome-wide analysis of gene expression or protein binding patterns using different array or sequencing based technologies is now routinely performed to compare different populations, such as treatment and reference groups. It is often necessary to normalize the data obtained to remove technical variation introduced in the course of conducting experimental work, but standard normalization techniques are not capable of eliminating technical bias in cases where the distribution of the truly altered variables is skewed, i.e. when a large fraction of the variables are either positively or negatively affected by the treatment. However, several experiments are likely to generate such skewed distributions, including ChIP-chip experiments for the study of chromatin, gene expression experiments for the study of apoptosis, and SNP-studies of copy number variation in normal and tumour tissues. A preliminary study using spike-in array data established that the capacity of an experiment to identify altered variables and generate unbiased estimates of the fold change decreases as the fraction of altered variables and the skewness increases. We propose the following work-flow for analyzing high-dimensional experiments with regions of altered variables: (1) Pre-process raw data using one of the standard normalization techniques. (2) Investigate if the distribution of the altered variables is skewed. (3) If the distribution is not believed to be skewed, no additional normalization is needed. Otherwise, re-normalize the data using a novel HMM-assisted normalization procedure. (4) Perform downstream analysis. Here, ChIP-chip data and simulated data were used to evaluate the performance of the work-flow. It was found that skewed distributions can be detected by using the novel DSE-test (Detection of Skewed Experiments). Furthermore, applying the HMM-assisted normalization to experiments where the distribution of the truly altered variables is skewed results in considerably higher sensitivity and lower bias than can be attained using standard and invariant normalization methods.

RoBuST: an integrated genomics resource for the root and bulb crop families Apiaceae and Alliaceae.

BACKGROUND: Root and bulb vegetables (RBV) include carrots, celeriac (root celery), parsnips (Apiaceae), onions, garlic, and leek (Alliaceae)--food crops grown globally and consumed worldwide. Few data analysis platforms are currently available where data collection, annotation and integration initiatives are focused on RBV plant groups. Scientists working on RBV include breeders, geneticists, taxonomists, plant pathologists, and plant physiologists who use genomic data for a wide range of activities including the development of molecular genetic maps, delineation of taxonomic relationships, and investigation of molecular aspects of gene expression in biochemical pathways and disease responses. With genomic data coming from such diverse areas of plant science, availability of a community resource focused on these RBV data types would be of great interest to this scientific community. DESCRIPTION: The RoBuST database has been developed to initiate a platform for collecting and organizing genomic information useful for RBV researchers. The current release of RoBuST contains genomics data for 294 Alliaceae and 816 Apiaceae plant species and has the following features: (1) comprehensive sequence annotations of 3663 genes 5959 RNAs, 22,723 ESTs and 11,438 regulatory sequence elements from Apiaceae and Alliaceae plant families; (2) graphical tools for visualization and analysis of sequence data; (3) access to traits, biosynthetic pathways, genetic linkage maps and molecular taxonomy data associated with Alliaceae and Apiaceae plants; and (4) comprehensive plant splice signal repository of 659,369 splice signals collected from 6015 plant species for comparative analysis of plant splicing patterns. CONCLUSIONS: RoBuST, available at http://robust.genome.com, provides an integrated platform for researchers to effortlessly explore and analyze genomic data associated with root and bulb vegetables.

AspAlt: A tool for inter-database, inter-genomic and user-specific comparative analysis of alternative transcription and alternative splicing in 46 eukaryotes.

We have developed AspAlt-a web-based comparative analytical platform for exploring the variations in alternative transcription (AT) events and alternative splicing (AS) events in eukaryotes. AspAlt provides integrated access to 2.1 million AT-AS annotations from 1,58,876 multi-isoform genes and has the following user-friendly analytical features: (1) advanced graphical display to visualize and analyze AT-AS events in 46 eukaryotic genomes; (2) compare and identify the differences in AT-AS patterns among a group of genes specified by the user or among homologous gene groups; (3) inter-database comparative viewer to analyze the differences in the AT-AS annotations for the same gene among Ensembl, RefSeq and AceView databases; (4) dynamically classify and generate graphical plots of AT-AS events from mRNA annotations submitted by the user; and (5) download genomic AT-AS annotations of 46 eukaryotes in XML and tab-delimited formats. The AspAlt resource is available at http://66.170.16.154/AspAlt.

ExDom: an integrated database for comparative analysis of the exon-intron structures of protein domains in eukaryotes.

We have developed ExDom, a unique database for the comparative analysis of the exon-intron structures of 96 680 protein domains from seven eukaryotic organisms (Homo sapiens, Mus musculus, Bos taurus, Rattus norvegicus, Danio rerio, Gallus gallus and Arabidopsis thaliana). ExDom provides integrated access to exon-domain data through a sophisticated web interface which has the following analytical capabilities: (i) intergenomic and intragenomic comparative analysis of exon-intron structure of domains; (ii) color-coded graphical display of the domain architecture of proteins correlated with their corresponding exon-intron structures; (iii) graphical analysis of multiple sequence alignments of amino acid and coding nucleotide sequences of homologous protein domains from seven organisms; (iv) comparative graphical display of exon distributions within the tertiary structures of protein domains; and (v) visualization of exon-intron structures of alternative transcripts of a gene correlated to variations in the domain architecture of corresponding protein isoforms. These novel analytical features are highly suited for detailed investigations on the exon-intron structure of domains and make ExDom a powerful tool for exploring several key questions concerning the function, origin and evolution of genes and proteins. ExDom database is freely accessible at: http://66.170.16.154/ExDom/.