Comparison of sequence-dependent tiling array normalization approaches

<p>Abstract</p> <p>Background</p> <p>The detection of enriched DNA or RNA fragments by tiling microarrays has become more and more popular. These microarrays contain a high number of small probes covering genomic loci. However, to achieve high coverage the probe sequences cannot be selected for their hybridization properties. The affinity of the probes towards their targets varies in a sequence-dependent manner. In order to remove this bias a number of approaches have been developed and shown to increase the detection of enriched DNA or RNA fragments. However, these approaches also employ a peak detection algorithm that is different from the one used previously. Thus, it seems possible that the enhancement of detection is due to the peak detection algorithm rather than the sequence-dependent normalization.</p> <p>Results</p> <p>We compared three different sequence-dependent probe level normalization procedures to a naïve sequence-independent normalization technique. In order to achieve maximal comparability, we used the normalized intensity values as input to a single peak detection algorithm. A so-called "spike-in" data set served as benchmark for the performance. We will show that the sequence-dependent normalization procedures do not perform better than the naïve approach, suggesting that the benefit of using these normalization approaches is limited. Furthermore, we will show that the naïve approach does well, because it effectively removes the sequence-dependent component of the measured intensities with the help of the control hybridization experiment.</p> <p>Conclusion</p> <p>Sequence-dependent normalization of microarray data hardly improves the detection of enriched DNA or RNA fragments. The "success" of the sequence-independent naïve approach is only possible due to the control experiment and requires proper scaling of the measured intensities.</p

Total Empiricism: Learning from Data

Statistical analysis is an important tool to distinguish systematic from chance findings. Current statistical analyses rely on distributional assumptions reflecting the structure of some underlying model, which if not met lead to problems in the analysis and interpretation of the results. Instead of trying to fix the model or "correct" the data, we here describe a totally empirical statistical approach that does not rely on ad hoc distributional assumptions in order to overcome many problems in contemporary statistics. Starting from elementary combinatorics, we motivate an information-guided formalism to quantify knowledge extracted from the given data. Subsequently, we derive model-agnostic methods to identify patterns that are solely evidenced by the data based on our prior knowledge. The data-centric character of empiricism allows for its universal applicability, particularly as sample size grows larger. In this comprehensive framework, we re-interpret and extend model distributions, scores and statistical tests used in different schools of statistics.Comment: Keywords: effective description, large-N, operator formalism, statistical testing, inference, information divergenc

The Zinc Finger Associated Domain of Drosophila melanogaster, its Evolution and Phylogenetic Restriction

The Drosophila melanogaster genome contains 359 C2H2 zinc finger protein (ZFP) coding genes. During an in-depth analysis of these proteins I identified a sequence motif that can be found at the N-terminus of 94 ZFP coding genes (26.1%). It has been named zinc finger associated domain (ZAD), since it is almost exclusively asociated with C2H2 zinc finger motifs, with the exception of four genes that code for proteins with isolated ZADs. The ZAD represents an independently folding domain that is stabilised by zinc coordination. It mediates homodimer formation and probably association of closely related ZAD family members. The ZAD coding genes fall into two large subsets: subset 1 where the ZAD is encoded by a single exon and subset 2 where the ZAD is coded by two exons. I propose that subset 1 ZADs have lost their intron at multiple time points by distinct retroposition events. The majority of ZAD coding genes in Drosophila melanogaster were generated by local gene duplication events, which is consistent with the observation that the ZAD coding genes participates in lineage-specific expansions in three fly species. The finding that many ZAD coding genes are specific for a given species prompted the idea that these genes may be involved in mechanisms that maintain genetic divergence during speciation. My results show that ZAD coding genes are overrepresented in a class of genes that represent putative "speciation genes".Das Drosophila melanogaster Genom enthält 359 C2H2 Zink Finger Protein (ZFP) kodierende Gene. Während der Analyse dieser Proteine konnte ich ein Sequenzmotiv identifizieren, welches am N-Terminus von 94 ZFP kodierenden Genen (26,1%) auftritt. Es wurde Zink Finger assoziierte Domäne genannt, da es fast auschließlich mit C2H2 Zink Finger Motiven assoziiert ist und nur bei vier Genen auftritt, die Proteine mit isolierten ZADs kodieren. Die ZAD ist eine unabhängig faltende Domäne, die durch Zink Koordination stabilisiert wird. Sie ermöglicht die Bildung von Homodimeren und wahrscheinlich die Bindung zwischen nah verwandten ZADs. Die ZADs kann man in zwei Klassen einteilen: Klasse 1, bei welchen die ZAD durch ein Exon kodiert wird, und Klasse 2, welche durch zwei Exons kodiert werden. Die ZADs der Klasse 1 haben ihr Intron wahrscheinlich mehrfach und unabhängig durch Retropositionsereignisse verloren. Die meisten ZAD kodierenden Gene von Drosophila melanogaster wurden durch lokale Genduplikationen erzeugt, was mit der Beobachtung übereinstimmt, dass die ZAD an "lineage-specific expansions" in drei Fliegenarten teilnimmt. Die Beobachtung, dass viele ZAD kodierende Gene spezifisch für eine gegebene Art sind, führte zu der Idee, dass diese an Mechanismen beteiligt sind, die die genetische Divergenz während des Artbildungsprozesses aufrecht erhalten. Meine Ergebnisse zeigen, dass ZAD kodierende Gene überproportional häufig in einer Klasse von möglichen "Speziationsgenen" vorkommen

The biological significance of non-enzymatic reaction of menadione with plasma thiols: enhancement of menadione-induced cytotoxicity to platelets by the presence of blood plasma

AbstractTo test the hypothesis that the non-enzymatic reaction of quinones with thiols in plasma can generate reactive oxygens (ROS), thereby leading to potentiated cellular toxicity, we have studied the effect of a representative quinone compound, menadione, on plasma isolated from rats. The experimental results are as follows: (1) menadione generated ROS via non-enzymatic reaction with protein thiols in plasma; (2) the presence of plasma increased menadione-induced cytotoxicity to platelets; (3) pretreatment of plasma with a thiol-depleting agent significantly suppressed menadione-induced ROS and cytotoxicity. These results suggest that the non-enzymatic reaction of menadione with plasma thiols could be an important process in quinone-induced cellular toxicity

reChIP-seq reveals widespread bivalency of H3K4me3 and H3K27me3 in CD4+ memory T cells

The combinatorial action of co-localizing chromatin modifications and regulators determines chromatin structure and function. However, identifying co-localizing chromatin features in a high-throughput manner remains a technical challenge. Here we describe a novel reChIP-seq approach and tailored bioinformatic analysis tool, normR that allows for the sequential enrichment and detection of co-localizing DNA-associated proteins in an unbiased and genome-wide manner. We illustrate the utility of the reChIP-seq method and normR by identifying H3K4me3 or H3K27me3 bivalently modified nucleosomes in primary human CD4+ memory T cells. We unravel widespread bivalency at hypomethylated CpG-islands coinciding with inactive promoters of developmental regulators. reChIP-seq additionally uncovered heterogeneous bivalency in the population, which was undetectable by intersecting H3K4me3 and H3K27me3 ChIP- seq tracks. Finally, we provide evidence that bivalency is established and stabilized by an interplay between the genome and epigenome. Our reChIP-seq approach augments conventional ChIP-seq and is broadly applicable to unravel combinatorial modes of action

Prevention of mitochondrial impairment by inhibition of protein phosphatase 1 activity in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease caused by progressive loss of motor neurons (MNs) and subsequent muscle weakness. These pathological features are associated with numerous cellular changes, including alteration in mitochondrial morphology and function. However, the molecular mechanisms associating mitochondrial structure with ALS pathology are poorly understood. In this study, we found that Dynamin-related protein 1 (Drp1) was dephosphorylated in several ALS models, including those with SOD1 and TDP-43 mutations, and the dephosphorylation was mediated by the pathological induction of protein phosphatase 1 (PP1) activity in these models. Suppression of the PP1-Drp1 cascade effectively prevented ALS-related symptoms, including mitochondrial fragmentation, mitochondrial complex I impairment, axonal degeneration, and cell death, in primary neuronal culture models, iPSC-derived human MNs, and zebrafish models in vivo. These results suggest that modulation of PP1-Drp1 activity may be a therapeutic target for multiple pathological features of ALS

Genetic Study in a Case of Birt-Hogg-Dubé Syndrome

Birt-Hogg-Dubé syndrome (BHDS) is an autosomal dominantly inherited disorder characterized by multiple trichodiscomas, fibrofolliculomas, and acrocordons. There is also an increased risk of developing renal neoplasms and lung cysts/spontaneous pneumothorax. We present a 43-year-old man with multiple, 2~4 mm sized, dome-shaped, and skin-colored papules on his cheek and neck. On the basis of clinical finding and histopathologic examination on the cheek lesion, it was diagnosed as multiple trichodiscomas. Subsequently, molecular analysis revealed a mutation in the folliculin gene. We report a rare case of BHDS with a proved gene mutation

Genome-wide H4 K16 acetylation by SAS-I is deposited independently of transcription and histone exchange

The MYST HAT Sas2 is part of the SAS-I complex that acetylates histone H4 lysine 16 (H4 K16Ac) and blocks the propagation of heterochromatin at the telomeres of Saccharomyces cerevisiae. In this study, we investigated Sas2-mediated H4 K16Ac on a genome-wide scale. Interestingly, H4 K16Ac loss in sas2Δ cells outside of the telomeric regions showed a distinctive pattern in that there was a pronounced decrease of H4 K16Ac within the majority of open reading frames (ORFs), but little change in intergenic regions. Furthermore, regions of low histone H3 exchange and low H3 K56 acetylation showed the most pronounced loss of H4 K16Ac in sas2Δ, indicating that Sas2 deposited this modification on chromatin independently of histone exchange. In agreement with the effect of Sas2 within ORFs, sas2Δ caused resistance to 6-azauracil, indicating a positive effect on transcription elongation in the absence of H4 K16Ac. In summary, our data suggest that Sas2-dependent H4 K16Ac is deposited into chromatin independently of transcription and histone exchange, and that it has an inhibitory effect on the ability of PolII to travel through the body of the gene

Gain-of-function screen for genes that affect Drosophila muscle pattern formation.

This article reports the production of an EP-element insertion library with more than 3,700 unique target sites within the Drosophila melanogaster genome and its use to systematically identify genes that affect embryonic muscle pattern formation. We designed a UAS/GAL4 system to drive GAL4-responsive expression of the EP-targeted genes in developing apodeme cells to which migrating myotubes finally attach and in an intrasegmental pattern of cells that serve myotubes as a migration substrate on their way towards the apodemes. The results suggest that misexpression of more than 1.5% of the Drosophila genes can interfere with proper myotube guidance and/or muscle attachment. In addition to factors already known to participate in these processes, we identified a number of enzymes that participate in the synthesis or modification of protein carbohydrate side chains and in Ubiquitin modifications and/or the Ubiquitin-dependent degradation of proteins, suggesting that these processes are relevant for muscle pattern formation