Search

Molecular interactions are key drivers of biological function. Providing interaction resources to the research community is important since they allow functional interpretation and network-based analysis of molecular data. ConsensusPathDB (http://consensuspathdb.org) is a meta-database combining interactions of diverse types from 31 public resources for humans, 16 for mice and 14 for yeasts. Using ConsensusPathDB, researchers commonly evaluate lists of genes, proteins and metabolites against sets of molecular interactions defined by pathways, Gene Ontology and network neighborhoods and retrieve complex molecular neighborhoods formed by heterogeneous interaction types. Furthermore, the integrated protein–protein interaction network is used as a basis for propagation methods. Here, we present the 2022 update of ConsensusPathDB, highlighting content growth, additional functionality and improved database stability. For example, the number of human molecular interactions increased to 859 848 connecting 200 499 unique physical entities such as genes/proteins, metabolites and drugs. Furthermore, we integrated regulatory datasets in the form of transcription factor–, microRNA– and enhancer–gene target interactions, thus providing novel functionality in the context of overrepresentation and enrichment analyses. We specifically emphasize the use of the integrated protein–protein interaction network as a scaffold for network inferences, present topological characteristics of the network and discuss strengths and shortcomings of such approaches.

Genetic predisposition affects the penetrance of tumor-initiating mutations, such as APC mutations that stabilize β-catenin and cause intestinal tumors in mice and humans. However, the mechanisms involved in genetically predisposed penetrance are not well understood. Here, we analyzed tumor multiplicity and gene expression in tumor-prone ApcMin/+ mice on highly variant C57BL/6J (B6) and PWD/Ph (PWD) genetic backgrounds. (B6 × PWD) F1 APCMin offspring mice were largely free of intestinal adenoma, and several chromosome substitution (consomic) strains carrying single PWD chromosomes on the B6 genetic background displayed reduced adenoma numbers. Multiple dosage-dependent modifier loci on PWD chromosome 5 each contributed to tumor suppression. Activation of β-catenin–driven and stem cell–specific gene expression in the presence of ApcMin or following APC loss remained moderate in intestines carrying PWD chromosome 5, suggesting that PWD variants restrict adenoma initiation by controlling stem cell homeostasis. Gene expression of modifier candidates and DNA methylation on chromosome 5 were predominantly cis controlled and largely reflected parental patterns, providing a genetic basis for inheritance of tumor susceptibility. Human SNP variants of several modifier candidates were depleted in colorectal cancer genomes, suggesting that similar mechanisms may also affect the penetrance of cancer driver mutations in humans. Overall, our analysis highlights the strong impact that multiple genetic variants acting in networks can exert on tumor development. Significance: These findings in mice show that, in addition to accidental mutations, cancer risk is determined by networks of individual gene variants.

The generation of peroxynitrite (ONOO−) is associated with several diseases, including atherosclerosis, hypertension, neurodegeneration, cancer, inflammation, and sepsis. Alpha-ketoglutarate (αKG) is a known potential highly antioxidative agent for radical oxidative species such as peroxides. The question arises as to whether αKG is also a potential scavenger of ONOO− and a potential protector against ONOO−-mediated nitration of proteins. NMR studies of 1 mM αKG in 100 mM phosphate-buffered saline at pH 7.4 and pH 6.0 were carried out in the presence or absence of a final concentration of 2 mM ONOO−. An ONOO−–luminol-induced chemiluminescence reaction was used to measure the scavenging function of several concentrations of αKG, quantification of αKG was performed via spectrophotometric enzymatic assay of αKG in the absence or presence of 0, 1, or 2 mM ONOO−. The nitration of tyrosine residues on proteins was measured on ONOO−-treated bovine serum albumin (BSA) in the presence or absence of 0–24 mM αKG by an ELISA technique using a specific anti-IgG against nitro-tyrosine. The addition of ONOO− to αKG led to the formation of succinic acid and nitrite at pH 7.0, but not at pH 6.0, as αKG was stable against ONOO−. The absorbance of enzymatically estimated αKG at the time point of 30 min was significantly lower in favour of ONOO− (1 mM: 0.21 ± 0.03, 2 mM: 0.12 ± 0.05 vs. 0 mM: 0.32 ± 0.02, p &lt, 0.001). The luminol technique showed an inverse logarithmic correlation of the ONOO− and αKG concentrations (y = −2 × 105 ln(x) + 1 × 106, r2 = 0.99). The usage of 4 mM αKG showed a significant reduction by nearly half in the chemiluminescence signal (284,456 ± 29,293 cps, p &lt, 0.001) compared to the control (474,401 ± 18,259), for 20 and 200 mM αKG, there were further reductions to 163,546 ± 26,196 cps (p &lt, 0.001) and 12,658 ± 1928 cps (p &lt, 0.001). Nitrated tyrosine residues were estimated using the ELISA technique. A negative linear correlation was obtained by estimating nitrated tyrosine residues in the presence of αKG (r2 = 0.94): a reduction by half of nitrated tyrosine was estimated using 12 mM αKG compared to the control (326.1 ± 39.6 nmol vs. 844.5 ± 128.4 nmol, 0.001).

Mutations in splicing factor genes have a severe impact on the survival of cancer patients. Splicing factor 3b subunit 1 (SF3B1) is one of the most frequently mutated genes in chronic lymphocytic leukemia (CLL), patients carrying these mutations have a poor prognosis. Since the splicing machinery and the epigenome are closely interconnected, we investigated whether these alterations may affect the epigenomes of CLL patients. While an overall hypomethylation during CLL carcinogenesis has been observed, the interplay between the epigenetic stage of the originating B cells and SF3B1 mutations, and the subsequent effect of the mutations on methylation alterations in CLL, have not been investigated. We profiled the genome-wide DNA methylation patterns of 27 CLL patients with and without SF3B1 mutations and identified local decreases in methylation levels in SF3B1mut CLL patients at 67 genomic regions, mostly in proximity to telomeric regions. These differentially methylated regions (DMRs) were enriched in gene bodies of cancer-related signaling genes, e.g., NOTCH1, HTRA3, and BCL9L. In our study, SF3B1 mutations exclusively emerged in two out of three epigenetic stages of the originating B cells. However, not all the DMRs could be associated with the methylation programming of B cells during development, suggesting that mutations in SF3B1 cause additional epigenetic aberrations during carcinogenesis.

SARS-CoV-2 has evolved rapidly towards higher infectivity and partial immune escape over the course of the pandemic. This evolution is driven by the enormous virus population, that has infected close to 200 million people by now. Therefore, cost effective and scalable methods are needed to monitor viral evolution globally. Mutation-specific PCR approaches have become inadequate to distinguish the variety of circulating SARS-CoV-2 variants and are unable to detect novel ones. Conversely, whole genome sequencing protocols remain too labor- and cost-intensive to monitor SARS-CoV-2 at the required density. By adapting SARSeq we present a simple, fast, and scalable S-gene tiling pipeline for focused sequencing of the S-gene encoding for the spike protein. This method reports on all sequence positions with known importance for infectivity and immunity, yet scales to >20K samples per run. S-gene tiling is used for nationwide surveillance of SARS-CoV-2 at a density of 10% to 50% of all cases of infection in Austria. SARSeq S-tiling uncovered several infection clusters with variants of concern such as the biggest known cluster of Beta/B.1.351 outside Africa and successfully informed public health measures in a timely manner, allowing their successful implementation. Our close monitoring of mutations further highlighted evolutionary constraints and freedom of the spike protein ectodomain and sheds light on foreseeable evolutionary trajectories.

Current molecular tumor diagnostics encompass panel sequencing to detect mutations, copy number alterations, and rearrangements. However, tumor suppressor genes can also be inactivated by methylation within their promoter region. These epigenetic alterations are so far rarely assessed in the clinical setting. Therefore, we established the AllCap protocol facilitating the combined detection of mutations and DNA methylation at the coding and promoter regions of 342 DNA repair genes in one experiment. We demonstrate the use of the protocol by applying it to ovarian cancer cell lines with different responsiveness to poly(ADP-ribose) polymerase inhibition. BRCA1, ATM, ATR, and EP300 mutations and methylation of the BRCA1 promoter were detected as potential predictors for therapy response. The required amount of input DNA was optimized, and the application to formalin-fixed, paraffin-embedded tissue samples was verified to improve the clinical applicability. Thus, by adding DNA methylation values to panel resequencings, the AllCap assay will add another important level of information to clinical tests and will improve stratification of patients for systemic therapies.

To fully understand human genetic variation and its functional consequences, the specific distribution of variants between the two chromosomal homologues of genes must be known. The ‘phase’ of variants can significantly impact gene function and phenotype. To assess patterns of phase at large scale, we have analyzed 18 121 autosomal genes in 1092 statistically phased genomes from the 1000 Genomes Project and 184 experimentally phased genomes from the Personal Genome Project. Here we show that genes with cis-configurations of coding variants are more frequent than genes with trans-configurations in a genome, with global cis/trans ratios of ∼60:40. Significant cis-abundance was observed in virtually all genomes in all populations. Moreover, we identified a large group of genes exhibiting cis-configurations of protein-changing variants in excess, so-called ‘cis-abundant genes’, and a smaller group of ‘trans-abundant genes’. These two gene categories were functionally distinguishable, and exhibited strikingly different distributional patterns of protein-changing variants. Underlying these phenomena was a shared set of phase-sensitive genes of importance for adaptation and evolution. This work establishes common patterns of phase as key characteristics of diploid human exomes and provides evidence for their functional significance, highlighting the importance of phase for the interpretation of protein-coding genetic variation and gene function.

Uncovering cellular responses from heterogeneous genomic data is crucial for molecular medicine in particular for drug safety. This can be realized by integrating the molecular activities in networks of interacting proteins. As proof-of-concept we challenge network modeling with time-resolved proteome, transcriptome and methylome measurements in iPSC-derived human 3D cardiac microtissues to elucidate adverse mechanisms of anthracycline cardiotoxicity measured with four different drugs (doxorubicin, epirubicin, idarubicin and daunorubicin). Dynamic molecular analysis at in vivo drug exposure levels reveal a network of 175 disease-associated proteins and identify common modules of anthracycline cardiotoxicity in vitro, related to mitochondrial and sarcomere function as well as remodeling of extracellular matrix. These in vitro-identified modules are transferable and are evaluated with biopsies of cardiomyopathy patients. This to our knowledge most comprehensive study on anthracycline cardiotoxicity demonstrates a reproducible workflow for molecular medicine and serves as a template for detecting adverse drug responses from complex omics data., Using a network propagation approach with integrated multi-omic data, Selevsek et al. develop a reproducible workflow for identifying drug toxicity effects in cellular systems. This is demonstrated with the analysis of anthracycline cardiotoxicity in cardiac microtissues under the effect of multiple drugs.

Objective Physical exercise training is associated with increased glucose uptake in skeletal muscle and improved glycemic control. HDAC5, a class IIa histone deacetylase, has been shown to regulate transcription of the insulin-responsive glucose transporter GLUT4 in cultured muscle cells. In this study, we analyzed the contribution of HDAC5 to the transcriptional network in muscle and the beneficial effect of muscle contraction and regular exercise on glucose metabolism. Methods HDAC5 knockout mice (KO) and wild-type (WT) littermates were trained for 8 weeks on treadmills, metabolically phenotyped, and compared to sedentary controls. Hdac5-deficient skeletal muscle and cultured Hdac5-knockdown (KD) C2C12 myotubes were utilized for studies of gene expression and glucose metabolism. Chromatin immunoprecipitation (ChIP) studies were conducted to analyze Il6 promoter activity using H3K9ac and HDAC5 antibodies. Results Global transcriptome analysis of Hdac5 KO gastrocnemius muscle demonstrated activation of the IL-6 signaling pathway. Accordingly, knockdown of Hdac5 in C2C12 myotubes led to higher expression and secretion of IL-6 with enhanced insulin-stimulated activation of AKT that was reversed by Il6 knockdown. Moreover, Hdac5-deficient myotubes exhibited enhanced glucose uptake, glycogen synthesis, and elevated expression levels of the glucose transporter GLUT4. Transcription of Il6 was further enhanced by electrical pulse stimulation in Hdac5-deficient C2C12 myotubes. ChIP identified a ∼1 kb fragment of the Il6 promoter that interacts with HDAC5 and demonstrated increased activation-associated histone marker AcH3K9 in Hdac5-deficient muscle cells. Exercise intervention of HDAC5 KO mice resulted in improved systemic glucose tolerance as compared to WT controls. Conclusions We identified HDAC5 as a negative epigenetic regulator of IL-6 synthesis and release in skeletal muscle. HDAC5 may exert beneficial effects through two different mechanisms, transcriptional control of genes required for glucose disposal and utilization, and HDAC5-dependent IL-6 signaling cross-talk to improve glucose uptake in muscle in response to exercise.

To identify novel disease genes for type 2 diabetes (T2D) we generated two backcross populations of obese and diabetes-susceptible New Zealand Obese (NZO/HI) mice with the two lean mouse strains 129P2/OlaHsd and C3HeB/FeJ. Subsequent whole-genome linkage scans revealed 30 novel quantitative trait loci (QTL) for T2D-associated traits. The strongest association with blood glucose [12 cM, logarithm of the odds (LOD) 13.3] and plasma insulin (17 cM, LOD 4.8) was detected on proximal chromosome 7 (designated Nbg7p, NZO blood glucose on proximal chromosome 7) exclusively in the NZOxC3H crossbreeding, suggesting that the causal gene is contributed by the C3H genome. Introgression of the critical C3H fragment into the genetic NZO background by generating recombinant congenic strains and metabolic phenotyping validated the phenotype. For the detection of candidate genes in the critical region (30–46 Mb), we used a combined approach of haplotype and gene expression analysis to search for C3H-specific gene variants in the pancreatic islets, which appeared to be the most likely target tissue for the QTL. Two genes, Atp4a and Pop4, fulfilled the criteria from our candidate gene approaches. The knockdown of both genes in MIN6 cells led to decreased glucose-stimulated insulin secretion, indicating a regulatory role of both genes in insulin secretion, thereby possibly contributing to the phenotype linked to Nbg7p. In conclusion, our combined- and comparative-cross analysis approach has successfully led to the identification of two novel diabetes susceptibility candidate genes, and thus has been proven to be a valuable tool for the discovery of novel disease genes.

Though clinical trials for medical applications of dimethyl sulfoxide (DMSO) reported toxicity in the 1960s, later, the FDA classified DMSO in the safest solvent category. DMSO became widely used in many biomedical fields and biological effects were overlooked. Meanwhile, biomedical science has evolved towards sensitive high-throughput techniques and new research areas, including epigenomics and microRNAs. Considering its wide use, especially for cryopreservation and in vitro assays, we evaluated biological effect of DMSO using these technological innovations. We exposed 3D cardiac and hepatic microtissues to medium with or without 0.1% DMSO and analyzed the transcriptome, proteome and DNA methylation profiles. In both tissue types, transcriptome analysis detected >2000 differentially expressed genes affecting similar biological processes, thereby indicating consistent cross-organ actions of DMSO. Furthermore, microRNA analysis revealed large-scale deregulations of cardiac microRNAs and smaller, though still massive, effects in hepatic microtissues. Genome-wide methylation patterns also revealed tissue-specificity. While hepatic microtissues demonstrated non-significant changes, findings from cardiac microtissues suggested disruption of DNA methylation mechanisms leading to genome-wide changes. The extreme changes in microRNAs and alterations in the epigenetic landscape indicate that DMSO is not inert. Its use should be reconsidered, especially for cryopreservation of embryos and oocytes, since it may impact embryonic development.

Autosomal recessive (AR) gene defects are the leading genetic cause of intellectual disability (ID) in countries with frequent parental consanguinity, which account for about 1/7th of the world population. Yet, compared to autosomal dominant de novo mutations, which are the predominant cause of ID in Western countries, the identification of AR-ID genes has lagged behind. Here, we report on whole exome and whole genome sequencing in 404 consanguineous predominantly Iranian families with two or more affected offspring. In 219 of these, we found likely causative variants, involving 77 known and 77 novel AR-ID (candidate) genes, 21 X-linked genes, as well as 9 genes previously implicated in diseases other than ID. This study, the largest of its kind published to date, illustrates that high-throughput DNA sequencing in consanguineous families is a superior strategy for elucidating the thousands of hitherto unknown gene defects underlying AR-ID, and it sheds light on their prevalence.

The assessment of the carcinogenic potential of chemicals with alternative, human-based in vitro systems has become a major goal of toxicogenomics. The central read-out of these assays is the transcriptome, and while many studies exist that explored the gene expression responses of such systems, reports on robustness and reproducibility, when testing them independently in different laboratories, are still uncommon. Furthermore, there is limited knowledge about variability induced by the data analysis protocols. We have conducted an inter-laboratory study for testing chemical carcinogenicity evaluating two human in vitro assays: hepatoma-derived cells and hTERT-immortalized renal proximal tubule epithelial cells, representing liver and kidney as major target organs. Cellular systems were initially challenged with thirty compounds, genome-wide gene expression was measured with microarrays, and hazard classifiers were built from this training set. Subsequently, each system was independently established in three different laboratories, and gene expression measurements were conducted using anonymized compounds. Data analysis was performed independently by two separate groups applying different protocols for the assessment of inter-laboratory reproducibility and for the prediction of carcinogenic hazard. As a result, both workflows came to very similar conclusions with respect to (1) identification of experimental outliers, (2) overall assessment of robustness and inter-laboratory reproducibility and (3) re-classification of the unknown compounds to the respective toxicity classes. In summary, the developed bioinformatics workflows deliver accurate measures for inter-laboratory comparison studies, and the study can be used as guidance for validation of future carcinogenicity assays in order to implement testing of human in vitro alternatives to animal testing.

Deletions at the C-terminus of the proto-oncogene protein c-Src kinase are found in the viral oncogene protein v-Src as well as in some advanced human colon cancers. They are associated with increased kinase activity and cellular invasiveness. Here, we analyzed the mRNA expression signature of a constitutively active C-terminal mutant of c-Src, c-Src(mt), in comparison with its wild-type protein, c-Src(wt), in the human non-transformed breast epithelial cell line MCF-10A. We demonstrated previously that the mutant altered migratory and metastatic properties. Genome-wide transcriptome analysis revealed that c-Src(mt) de-regulated the expression levels of approximately 430 mRNAs whose gene products are mainly involved in the cellular processes of migration and adhesion, apoptosis and protein synthesis. 82.9% of these genes have previously been linked to cellular migration, while the others play roles in RNA transport and splicing processes, for instance. Consistent with the transcriptome data, cells expressing c-Src(mt), but not those expressing c-Src(wt), showed the capacity to metastasize into the lungs of mice in vivo. The mRNA expression profile of c-Src(mt)-expressing cells shows significant overlap with that of various primary human tumor samples, possibly reflecting elevated Src activity in some cancerous cells. Expression of c-Src(mt) led to elevated migratory potential. We used this model system to analyze the transcriptional changes associated with an invasive cellular phenotype. These genes and pathways de-regulated by c-Src(mt) may provide suitable biomarkers or targets of therapeutic approaches for metastatic cells. DATABASE: This project was submitted to the National Center for Biotechnology Information BioProject under ID PRJNA288540. The Illumina RNA-Seq reads are available in the National Center for Biotechnology Information Sequence Read Archive under study ID SRP060008 with accession numbers SRS977414 for MCF-10A cells, SRS977717 for mock cells, SRS978053 for c-Src(wt) cells and SRS978046 for c-Src(mt) cells.

We present NetCore, a novel network propagation approach based on node coreness, for phenotype–genotype associations and module identification. NetCore addresses the node degree bias in PPI networks by using node coreness in the random walk with restart procedure, and achieves improved re-ranking of genes after propagation. Furthermore, NetCore implements a semi-supervised approach to identify phenotype-associated network modules, which anchors the identification of novel candidate genes at known genes associated with the phenotype. We evaluated NetCore on gene sets from 11 different GWAS traits and showed improved performance compared to the standard degree-based network propagation using cross-validation. Furthermore, we applied NetCore to identify disease genes and modules for Schizophrenia GWAS data and pan-cancer mutation data. We compared the novel approach to existing network propagation approaches and showed the benefits of using NetCore in comparison to those. We provide an easy-to-use implementation, together with a high confidence PPI network extracted from ConsensusPathDB, which can be applied to various types of genomics data in order to obtain a re-ranking of genes and functionally relevant network modules.

To fully understand human genetic variation, one must assess the specific distribution of variants between the two chromosomal homologues of genes, and any functional units of interest, as the phase of variants can significantly impact gene function and phenotype. To this end, we have systematically analyzed 18,121 autosomal protein-coding genes in 1,092 statistically phased genomes from the 1000 Genomes Project, and an unprecedented number of 184 experimentally phased genomes from the Personal Genome Project. Here we show that mutations predicted to functionally alter the protein, and coding variants as a whole, are not randomly distributed between the two homologues of a gene, but do occur significantly more frequently in cis-than trans-configurations, with cis/trans ratios of ∼60:40. Significant cis-abundance was observed in virtually all individual genomes in all populations. Nearly all variable genes exhibited either cis, or trans configurations of protein-altering mutations in significant excess, allowing distinction of cis- and trans-abundant genes. These common patterns of phase were largely constituted by a shared, global set of phase-sensitive genes. We show significant enrichment of this global set with gene sets indicating its involvement in adaptation and evolution. Moreover, cis- and trans-abundant genes were found functionally distinguishable, and exhibited strikingly different distributional patterns of protein-altering mutations. This work establishes common patterns of phase as key characteristics of diploid human exomes and provides evidence for their potential functional significance. Thus, it highlights the importance of phase for the interpretation of protein-coding genetic variation, challenging the current conceptual and functional interpretation of autosomal genes.

Current testing models for predicting drug-induced liver injury are inadequate, as they basically under-report human health risks. We present here an approach towards developing pathways based on hepatotoxicity-associated gene groups derived from two types of publicly accessible hepatotoxicity databases, in order to develop drug-induced liver injury biomarker profiles. One human liver ‘omics-based and four text-mining-based databases were explored for hepatotoxicity-associated gene lists. Over-representation analysis of these gene lists with a hepatotoxicant-exposed primary human hepatocytes data set showed that human liver ‘omics gene lists performed better than text-mining gene lists and the results of the latter differed strongly between databases. However, both types of databases contained gene lists demonstrating biomarker potential. Visualizing those in pathway format may aid in interpreting the biomolecular background. We conclude that exploiting existing and openly accessible databases in a dedicated manner seems promising in providing venues for translational research in toxicology and biomarker development.

Actinobacillus pleuropneumoniae is among the most important pathogens worldwide in pig production. The agent can cause severe economic losses due to decreased performance, acute or chronic pleuropneumonia and an increased incidence of death. Therapeutics cannot be used in a sustainable manner, and vaccination is not always available, but discovering more about host defence and disease mechanisms might lead to new methods of prophylaxis. The aim of the present study was to detect quantitative trait loci (QTL) associated with resistance/susceptibility to A. pleuropneumoniae. Under controlled conditions, 170 F2 animals of a Hampshire/Landrace family, with known differences in founder populations regarding A. pleuropneumoniae resistance, were challenged with an A. pleuropneumoniae serotype 7 aerosol followed by a detailed clinical, radiographic, ultrasonographic, pathological and bacteriological examination. F2 pigs were genotyped with 159 microsatellite markers. Significant QTL were identified on Sus scrofa chromosomes (SSC) 2, 6, 12, 13, 16, 17 and 18. They explained 6–22 % of phenotypic variance. One QTL on SSC2 reached significance on a genome-wide level for five associated phenotypic traits. A multiple regression analysis revealed a combinatory effect of markers SWR345 (SSC2) and S0143 (SSC12) on Respiratory Health Score, Clinical Score and the occurrence of death. The results indicate the genetic background of A. pleuropneumoniae resistance in swine and provide new insights into the genetic architecture of resistance/susceptibility to porcine pleuropneumonia. The results will be helpful in identifying the underlying genes and mechanisms.

The yeast two-hybrid (Y2H) system is the most widely applied methodology for systematic protein-protein interaction (PPI) screening and the generation of comprehensive interaction networks. We developed a novel Y2H interaction screening procedure using DNA microarrays for high-throughput quantitative PPI detection. Applying a global pooling and selection scheme to a large collection of human open reading frames, proof-of-principle Y2H interaction screens were performed for the human neurodegenerative disease proteins huntingtin and ataxin-1. Using systematic controls for unspecific Y2H results and quantitative benchmarking, we identified and scored a large number of known and novel partner proteins for both huntingtin and ataxin-1. Moreover, we show that this parallelized screening procedure and the global inspection of Y2H interaction data are uniquely suited to define specific PPI patterns and their alteration by disease-causing mutations in huntingtin and ataxin-1. This approach takes advantage of the specificity and flexibility of DNA microarrays and of the existence of solid-related statistical methods for the analysis of DNA microarray data, and allows a quantitative approach toward interaction screens in human and in model organisms.

The use of hybridisation of synthetic oligonucleotides to cDNAs under high stringency to characterise gene sequences has been demonstrated by a number of groups. We have used two cDNA libraries of 9 and 12 day mouse embryos (24 133 and 34 783 clones respectively) in a pilot study to characterise expressed genes by hybridisation with 110 hybridisation probes. We have identified 33 369 clusters of cDNA clones, that ranged in representation from 1 to 487 copies (0.7%). 737 were assigned to known rodent genes, and a further 13 845 showed significant homologies. A total of 404 clusters were identified as significantly differentially represented (P < 0.01) between the two cDNA libraries. This study demonstrates the utility of the fingerprinting approach for the generation of comparative gene expression profiles through the analysis of cDNAs derived from different biological materials.

When evaluating compound similarity, addressing multiple sources of information to reach conclusions about common pharmaceutical and/or toxicological mechanisms of action is a crucial strategy. In this chapter, we describe a systems biology approach that incorporates analyses of hepatotoxicant data for 33 compounds from three different sources: a chemical structure similarity analysis based on the 3D Tanimoto coefficient, a chemical structure-based protein target prediction analysis, and a cross-study/cross-platform meta-analysis of in vitro and in vivo human and rat transcriptomics data derived from public resources (i.e., the diXa data warehouse). Hierarchical clustering of the outcome scores of the separate analyses did not result in a satisfactory grouping of compounds considering their known toxic mechanism as described in literature. However, a combined analysis of multiple data types may hypothetically compensate for missing or unreliable information in any of the single data types. We therefore performed an integrated clustering analysis of all three data sets using the R-based tool iClusterPlus. This indeed improved the grouping results. The compound clusters that were formed by means of iClusterPlus represent groups that show similar gene expression while simultaneously integrating a similarity in structure and protein targets, which corresponds much better with the known mechanism of action of these toxicants. Using an integrative systems biology approach may thus overcome the limitations of the separate analyses when grouping liver toxicants sharing a similar mechanism of toxicity.

Hepatocyte-like cells derived from the differentiation of human embryonic stem cells (hES-Hep) have potential to provide a human relevant in vitro test system in which to evaluate the carcinogenic hazard of chemicals. In this study, we have investigated this potential using a panel of 15 chemicals classified as noncarcinogens, genotoxic carcinogens, and nongenotoxic carcinogens and measured whole-genome transcriptome responses with gene expression microarrays. We applied an ANOVA model that identified 592 genes highly discriminative for the panel of chemicals. Supervised classification with these genes achieved a cross-validation accuracy of > 95%. Moreover, the expression of the response genes in hES-Hep was strongly correlated with that in human primary hepatocytes cultured in vitro. In order to infer mechanistic information on the consequences of chemical exposure in hES-Hep, we developed a computational method that measures the responses of biochemical pathways to the panel of treatments and showed that these responses were discriminative for the three toxicity classes and linked to carcinogenesis through p53, mitogen-activated protein kinases, and apoptosis pathway modules. It could further be shown that the discrimination of toxicity classes was improved when analyzing the microarray data at the pathway level. In summary, our results demonstrate, for the first time, the potential of human embryonic stem cell--derived hepatic cells as an in vitro model for hazard assessment of chemical carcinogenesis, although it should be noted that more compounds are needed to test the robustness of the assay.

Aims/hypothesis Numerous new genes have recently been identified in genome-wide association studies for type 2 diabetes. Most are highly expressed in beta cells and presumably play important roles in their function. However, these genes account for only a small proportion of total risk and there are likely to be additional candidate genes not detected by current methodology. We therefore investigated islets from the polygenic New Zealand mouse (NZL) model of diet-induced beta cell dysfunction to identify novel genes and pathways that may play a role in the pathogenesis of diabetes. Methods NZL mice were fed a diabetogenic high-fat diet (HF) or a diabetes-protective carbohydrate-free HF diet (CHF). Pancreatic islets were isolated by laser capture microdissection (LCM) and subjected to genome-wide transcriptome analyses. Results In the prediabetic state, 2,109 islet transcripts were differentially regulated (>1.5-fold) between HF and CHF diets. Of the genes identified, 39 (e.g. Cacna1d, Chd2, Clip2, Igf2bp2, Dach1, Tspan8) correlated with data from the Diabetes Genetics Initiative and Wellcome Trust Case Control Consortium genome-wide scans for type 2 diabetes, thus validating our approach. HF diet induced early changes in gene expression associated with increased cell-cycle progression, proliferation and differentiation of islet cells, and oxidative stress (e.g. Cdkn1b, Tmem27, Pax6, Cat, Prdx4 and Txnip). In addition, pathway analysis identified oxidative phosphorylation as the predominant gene-set that was significantly upregulated in response to the diabetogenic HF diet. Conclusions/interpretation We demonstrated that LCM of pancreatic islet cells in combination with transcriptional profiling can be successfully used to identify novel candidate genes for diabetes. Our data strongly implicate glucose-induced oxidative stress in disease progression. Electronic supplementary material The online version of this article (doi:10.1007/s00125-009-1576-4) contains supplementary material, which is available to authorised users.

ConsensusPathDB is a database system for the integration of human functional interactions. Current knowledge of these interactions is dispersed in more than 200 databases, each having a specific focus and data format. ConsensusPathDB currently integrates the content of 12 different interaction databases with heterogeneous foci comprising a total of 26,133 distinct physical entities and 74,289 distinct functional interactions (protein-protein interactions, biochemical reactions, gene regulatory interactions), and covering 1738 pathways. We describe the database schema and the methods used for data integration. Furthermore, we describe the functionality of the ConsensusPathDB web interface, where users can search and visualize interaction networks, upload, modify and expand networks in BioPAX, SBML or PSI-MI format, or carry out over-representation analysis with uploaded identifier lists with respect to substructures derived from the integrated interaction network. The ConsensusPathDB database is available at: http://cpdb.molgen.mpg.de.

The present study investigated mRNA expression profiles of bovine oocytes and blastocysts by using a cross-species hybridization approach employing an array consisting of 15,529 human cDNAs as probe, thus enabling the identification of conserved genes during human and bovine preimplantation development. Our analysis revealed 419 genes that were expressed in both oocytes and blastocysts. The expression of 1,324 genes was detected exclusively in the blastocyst, in contrast to 164 in the oocyte including a significant number of novel genes. Genes indicative for transcriptional and translational control ( ELAVL4, TACC3) were overexpressed in the oocyte, whereas cellular trafficking ( SLC2A14, SLC1A3), proteasome ( PSMA1, PSMB3), cell cycle ( BUB3, CCNE1, GSPT1), and protein modification and turnover ( TNK1, UBE3A) genes were found to be overexpressed in blastocysts. Transcripts implicated in chromatin remodeling were found in both oocytes ( NASP, SMARCA2) and blastocysts ( H2AFY, HDAC7A). The trophectodermal markers PSG2 and KRT18 were enriched 5- and 50-fold in the blastocyst. Pathway analysis revealed differential expression of genes involved in 107 distinct signaling and metabolic pathways. For example, phosphatidylinositol signaling and gluconeogenesis were prominent pathways identified in the blastocyst. Expression patterns in bovine and human blastocysts were to a large extent identical. This analysis compared the transcriptomes of bovine oocytes and blastocysts and provides a solid foundation for future studies on the first major differentiation events in blastocysts and identification of a set of markers indicative for regular mammalian development.

Modelling and simulation techniques are valuable tools for the understanding of complex biological systems. The design of a computer model necessarily has many diverse inputs, such as information on the model topology, reaction kinetics and experimental data, derived either from the literature, databases or direct experimental investigation. In this review, we describe different data resources, standards and modelling and simulation tools that are relevant to integrative systems biology.

Nodal signals induce mesodermal and endodermal progenitors during vertebrate development. To determine the role of Nodal signaling at a genomic level, we isolated Nodal-regulated genes by expression profiling using macroarrays and gene expression databases. Putative Nodal-regulated genes were validated by in situ hybridization screening in wild type and Nodal signaling mutants. 46 genes were identified, raising the currently known number of Nodal-regulated genes to 72. Based on their expression patterns along the dorsoventral axis, most of these genes can be classified into two groups. One group is expressed in the dorsal margin, whereas the other group is expressed throughout the margin. In addition to transcription factors and signaling components, the screens identified several new functional classes of Nodal-regulated genes, including cytoskeletal components and molecules involved in protein secretion or endoplasmic reticulum stress. We found that x-box binding protein-1 (xbp1) is a direct target of Nodal signaling and required for the terminal differentiation of the hatching gland, a specialized secretory organ whose specification is also dependent on Nodal signaling. These results indicate that Nodal signaling regulates not only specification genes but also differentiation genes.

Background Leishmania (L.) are intracellular protozoan parasites able to survive and replicate in the hostile phagolysosomal environment of infected macrophages. They cause leishmaniasis, a heterogeneous group of worldwide-distributed affections, representing a paradigm of neglected diseases that are mainly embedded in impoverished populations. To establish successful infection and ensure their own survival, Leishmania have developed sophisticated strategies to subvert the host macrophage responses. Despite a wealth of gained crucial information, these strategies still remain poorly understood. MicroRNAs (miRNAs), an evolutionarily conserved class of endogenous 22-nucleotide non-coding RNAs, are described to participate in the regulation of almost every cellular process investigated so far. They regulate the expression of target genes both at the levels of mRNA stability and translation; changes in their expression have a profound effect on their target transcripts. Methodology/Principal Findings We report in this study a comprehensive analysis of miRNA expression profiles in L. major-infected human primary macrophages of three healthy donors assessed at different time-points post-infection (three to 24 h). We show that expression of 64 out of 365 analyzed miRNAs was consistently deregulated upon infection with the same trends in all donors. Among these, several are known to be induced by TLR-dependent responses. GO enrichment analysis of experimentally validated miRNA-targeted genes revealed that several pathways and molecular functions were disturbed upon parasite infection. Finally, following parasite infection, miR-210 abundance was enhanced in HIF-1α-dependent manner, though it did not contribute to inhibiting anti-apoptotic pathways through pro-apoptotic caspase-3 regulation. Conclusions/Significance Our data suggest that alteration in miRNA levels likely plays an important role in regulating macrophage functions following L. major infection. These results could contribute to better understanding of the dynamics of gene expression in host cells during leishmaniasis., Author Summary Leishmania parasites belong to different species, each one characterized by specific vectors and reservoirs, and causes cutaneous or visceral disease(s) of variable clinical presentation and severity. In its mammalian host, the parasite is an obligate intracellular pathogen infecting the monocyte/macrophage lineage. Leishmania have developed ambiguous relationships with macrophages. Indeed, these cells are the shelter of invading parasites, where they will grow and eventually will reside in a silent state for life. But macrophages are also the cells that participate, through the induction of several pro-inflammatory mediators and antigen presentation, to shape the host immune response and ultimately kill the invader. To subvert these anti-parasite responses, Leishmania manipulate the host machinery for their own differentiation and survival. We aimed to evaluate the impact of L. major (the causative agent of zoonotic cutaneous leishmaniasis) infection on deregulation of non-coding miRNAs, a class of important regulators of gene expression. Our results revealed the implication of several miRNAs on macrophage fate upon parasite infection through regulation of different pathways, including cell death. Our findings provided a new insight for understanding mechanisms governing this miRNA deregulation by parasite infection and will help to provide clues for the development of control strategies for this disease.

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies., Wellcome Trust (London, England) (Core Award 090532/Z/09/Z), Wellcome Trust (London, England) (Senior Investigator Award 095552/Z/11/Z ), Wellcome Trust (London, England) (WT095908), Wellcome Trust (London, England) (WT109497), Wellcome Trust (London, England) (WT098051), Wellcome Trust (London, England) (WT086084/Z/08/Z), Wellcome Trust (London, England) (WT100956/Z/13/Z ), Wellcome Trust (London, England) (WT097307), Wellcome Trust (London, England) (WT0855322/Z/08/Z ), Wellcome Trust (London, England) (WT090770/Z/09/Z ), Wellcome Trust (London, England) (Major Overseas program in Vietnam grant 089276/Z.09/Z), Medical Research Council (Great Britain) (grant G0801823), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I02593X/1), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I021213/1), Zhongguo ke xue ji shu qing bao yan jiu suo. Office of 863 Programme of China (2012AA02A201), National Basic Research Program of China (2011CB809201), National Basic Research Program of China (2011CB809202), National Basic Research Program of China (2011CB809203), National Natural Science Foundation of China (31161130357), Shenzhen Municipal Government of China (grant ZYC201105170397A), Canadian Institutes of Health Research (grant 136855), Quebec Ministry of Economic Development, Innovation, and Exports (PSR-SIIRI-195), Germany. Bundesministerium für Bildung und Forschung (0315428A), Germany. Bundesministerium für Bildung und Forschung (01GS08201), Germany. Bundesministerium für Bildung und Forschung (BMBF-EPITREAT grant 0316190A), Deutsche Forschungsgemeinschaft (Emmy Noether Grant KO4037/1-1), Beatriu de Pinos Program (2006 BP-A 10144), Beatriu de Pinos Program (2009 BP-B 00274), Spanish National Institute for Health (grant PRB2 IPT13/0001-ISCIII-SGEFI/FEDER), Japan Society for the Promotion of Science (fellowship number PE13075), Marie Curie Actions Career Integration (grant 303772), Fonds National Suisse del la Recherche, SNSF, Scientifique (31003A_130342), National Center for Biotechnology Information (U.S.) (U54HG3067), National Center for Biotechnology Information (U.S.) (U54HG3273), National Center for Biotechnology Information (U.S.) (U01HG5211), National Center for Biotechnology Information (U.S.) (U54HG3079), National Center for Biotechnology Information (U.S.) (R01HG2898), National Center for Biotechnology Information (U.S.) (R01HG2385), National Center for Biotechnology Information (U.S.) (RC2HG5552), National Center for Biotechnology Information (U.S.) (U01HG6513), National Center for Biotechnology Information (U.S.) (U01HG5214), National Center for Biotechnology Information (U.S.) (U01HG5715), National Center for Biotechnology Information (U.S.) (U01HG5718), National Center for Biotechnology Information (U.S.) (U01HG5728), National Center for Biotechnology Information (U.S.) (U41HG7635), National Center for Biotechnology Information (U.S.) (U41HG7497), National Center for Biotechnology Information (U.S.) (R01HG4960), National Center for Biotechnology Information (U.S.) (R01HG5701), National Center for Biotechnology Information (U.S.) (R01HG5214), National Center for Biotechnology Information (U.S.) (R01HG6855), National Center for Biotechnology Information (U.S.) (R01HG7068), National Center for Biotechnology Information (U.S.) (R01HG7644), National Center for Biotechnology Information (U.S.) (DP2OD6514), National Center for Biotechnology Information (U.S.) (DP5OD9154), National Center for Biotechnology Information (U.S.) (R01CA166661), National Center for Biotechnology Information (U.S.) (R01CA172652), National Center for Biotechnology Information (U.S.) (P01GM99568), National Center for Biotechnology Information (U.S.) (R01GM59290), National Center for Biotechnology Information (U.S.) (R01GM104390), National Center for Biotechnology Information (U.S.) (T32GM7790), National Center for Biotechnology Information (U.S.) (R01HL87699), National Center for Biotechnology Information (U.S.) (R01HL104608), National Center for Biotechnology Information (U.S.) (T32HL94284), National Center for Biotechnology Information (U.S.) (HHSN268201100040C), National Center for Biotechnology Information (U.S.) (HHSN272201000025C), Lundbeck Foundation (grant R170-2014-1039, Simons Foundation (SFARI award SF51), National Science Foundation (U.S.) (Research Fellowship DGE-1147470)

The importance of hormones in endogenous aging has been displayed by recent studies performed on animal models and humans. To decipher the molecular mechanisms involved in aging we maintained human sebocytes at defined hormone-substituted conditions that corresponded to average serum levels of females from 20 (f20) to 60 (f60) years of age. The corresponding hormone receptor expression was demonstrated by reverse transcription-polymerase chain reaction (RT-PCR), Western blotting and immunocytochemistry. Cells at f60 produced significantly lower lipids than at f20. Increased mRNA and protein levels of c-Myc and increased protein levels of FN1, which have been associated with aging, were detected in SZ95 sebocytes at f60 compared to those detected at f20 after 5 days of treatment. Expression profiling employing a cDNA microarray composed of 15 529 cDNAs identified 899 genes with altered expression levels at f20 vs. f60. Confirmation of gene regulation was performed by real-time RT-PCR. The functional annotation of these genes according to the Gene Ontology identified pathways related to mitochondrial function, oxidative stress, ubiquitin-mediated proteolysis, cell cycle, immune responses, steroid biosynthesis and phospholipid degradation - all hallmarks of aging. Twenty-five genes in common with those identified in aging kidneys and several genes involved in neurodegenerative diseases were also detected. This is the first report describing the transcriptome of human sebocytes and its modification by a cocktail of hormones administered in age-specific levels and provides an in vitro model system, which approximates some of the hormone-dependent changes in gene transcription that occur during aging in humans.

Acne vulgaris is a skin disorder of the sebaceous follicles, involving hyperkeratinization and perifollicular inflammation. Matrix metalloproteinases (MMP) have a predominant role in inflammatory matrix remodeling and hyperproliferative skin disorders. We investigated the expression of MMP and tissue inhibitors of MMP (TIMP) in facial sebum specimens from acne patients, before and after treatment with isotretinoin. Gelatin zymography and Western-blot analysis revealed that sebum contains proMMP-9, which was decreased following per os or topical treatment with isotretinoin and in parallel to the clinical improvement of acne. Sebum also contains MMP-1, MMP-13, TIMP-1, and TIMP-2, as assessed by ELISA and western blot, but only MMP-13 was decreased following treatment with isotretinoin. The origin of MMP and TIMP in sebum is attributed to keratinocytes and sebocytes, since we found that HaCaT keratinocytes in culture secrete proMMP-2, proMMP-9, MMP-1, MMP-13, TIMP-1, and TIMP-2. SZ95 sebocytes in culture secreted proMMP-2 and proMMP-9, which was also confirmed by microarray analysis. Isotretinoin inhibited the arachidonic acid-induced secretion and mRNA expression of proMMP-2 and -9 in both cell types and of MMP-13 in HaCaT keratinocytes. These data indicate that MMP and TIMP of epithelial origin may be involved in acne pathogenesis, and that isotretinoin-induced reduction in MMP-9 and -13 may contribute to the therapeutic effects of the agent in acne.

Patients: In order to determine if there are areas of major and minor perfusion in a single testicle and if the quality of sperm is correlated with quantity of perfusion we collected testicle tissue for TESE in accordance to the local testicle tissue perfusion. Methods: A patient undergoing TESE underwent testicular perfusion mapping using contrast enhanced ultrasound. The exposed tissue was scanned with a Laser Doppler scanner and perfusion rates were determined measuring tissue perfusion units (TPUs). Tissue was biopsied and sperm were selected and prepared for assisted reproduction. Results: The total amount of isolated sperm correlated highly with the intensity of tissue perfusion showing high number of sperm in areas with high TPUs. Conclusions: This is the first demonstration that sperm quality and quantity is depending on tissue perfusion within the testicle. To further improve infertility treatment we propose that random biopsies could be replaced by perfusion-dependent collection of testicular tissue.

Differential gene transcription is a fundamental regulatory mechanism of biological systems during development, body homeostasis, and disease. Comparative genomics is believed to be a rapid means for the identification of regulatory sequences in genomes. We tested this approach to identify regulatory sequences that control expression in the midline of the zebrafish embryo. We first isolated a set of genes that are coexpressed in the midline of the zebrafish embryo during somitogenesis stages by gene array analysis and subsequent rescreens by in situ hybridization. We subjected 45 of these genes to Compare and DotPlot analysis to detect conserved sequences in noncoding regions of orthologous loci in the zebrafish and Takifugu genomes. The regions of homology that were scored in nonconserved regions were inserted into expression vectors and tested for their regulatory activity by transient transgenesis in the zebrafish embryo. We identified one conserved region from the connective tissue growth factor gene (ctgf), which was able to drive expression in the midline of the embryo. This region shares sequence similarity with other floor plate/notochord-specific regulatory regions. Our results demonstrate that an unbiased comparative approach is a relevant method for the identification of tissue-specific cis-regulatory sequences in the zebrafish embryo.

Together with the hemichordates, sea urchins represent basal groups of nonchordate invertebrate deuterostomes that occupy a key position in bilaterian evolution. Because sea urchin embryos are also amenable to functional studies, the sea urchin system has emerged as one of the leading models for the analysis of the function of genomic regulatory networks that control development. We have analyzed a total of 107,283 cDNA clones of libraries that span the development of the sea urchin Strongylocentrotus purpuratus. Normalization by oligonucleotide fingerprinting, EST sequencing and sequence clustering resulted in an EST catalog comprised of 20,000 unique genes or gene fragments. Around 7000 of the unique EST consensus sequences were associated with molecular and developmental functions. Phylogenetic comparison of the identified genes to the genome of the urochordate Ciona intestinalis indicate that at least one quarter of the genes thought to be chordate specific were already present at the base of deuterostome evolution. Comparison of the number of gene copies in sea urchins to those in chordates and vertebrates indicates that the sea urchin genome has not undergone extensive gene or complete genome duplications. The established unique gene set represents an essential tool for the annotation and assembly of the forthcoming sea urchin genome sequence. All cDNA clones and filters of all analyzed libraries are available from the resource center of the German genome project at http://www.rzpd.de.

The 2R hypothesis predicting two genome duplications at the origin of vertebrates is highly controversial. Studies published so far include limited sequence data from organisms close to the hypothesized genome duplications. Through the comparison of a gene catalog from amphioxus, the closest living invertebrate relative of vertebrates, to 3453 single-copy genes orthologous between Caenorhabditis elegans (C), Drosophila melanogaster (D), and Saccharomyces cerevisiae (Y), and to Ciona intestinalis ESTs, mouse, and human genes, we show with a large number of genes that the gene duplication activity is significantly higher after the separation of amphioxus and the vertebrate lineages, which we estimate at 650 million years (Myr). The majority of human orthologs of 195 CDY groups that could be dated by the molecular clock appear to be duplicated between 300 and 680 Myr with a mean at 488 million years ago (Mya). We detected 485 duplicated chromosomal segments in the human genome containing CDY orthologs, 331 of which are found duplicated in the mouse genome and within regions syntenic between human and mouse, indicating that these were generated earlier than the human–mouse split. Model based calculations of the codon substitution rate of the human genes included in these segments agree with the molecular clock duplication time-scale prediction. Our results favor at least one large duplication event at the origin of vertebrates, followed by smaller scale duplication closer to the bird–mammalian split.

Access to the complete gene inventory of an organism is crucial to understanding physiological processes like development, differentiation, pathogenesis, or adaptation to the environment. Transcripts from many active genes are present at low copy numbers. Therefore, procedures that rely on random EST sequencing or on normalisation and subtraction methods have to produce massively redundant data to get access to low-abundance genes. Here, we present an improved oligonucleotide fingerprinting (ofp) approach to the genome of sugar beet (Beta vulgaris), a plant for which practically no molecular information has been available. To identify distinct genes and to provide a representative 'unigene' cDNA set for sugar beet, 159 936 cDNA clones were processed utilizing large-scale, high-throughput data generation and analysis methods. Data analysis yielded 30 444 ofp clusters reflecting the number of different genes in the original cDNA sample. A sample of 10 961 cDNA clones, each representing a different cluster, were selected for sequencing. Standard sequence analysis confirmed that 89% of these EST sequences did represent different genes. These results indicate that the full set of 30 444 ofp clusters represent up to 25 000 genes. We conclude that the ofp analysis pipeline is an accurate and effective way to construct large representative 'unigene' sets for any plant of interest with no requirement for prior molecular sequence data.

We developed a novel efficient scheme, DEFOG (for "deciphering families of genes"), for determining sequences of numerous genes from a family of interest. The scheme provides a powerful means to obtain a gene family composition in species for which high-throughput genomic sequencing data are not available. DEFOG uses two key procedures. The first is a novel algorithm for designing highly degenerate primers based on a set of known genes from the family of interest. These primers are used in PCR reactions to amplify the members of the gene family. The second combines oligofingerprinting of the cloned PCR products with clustering of the clones based on their fingerprints. By selecting members from each cluster, a low-redundancy clone subset is chosen for sequencing. We applied the scheme to the human olfactory receptor (OR) genes. OR genes constitute the largest gene superfamily in the human genome, as well as in the genomes of other vertebrate species. DEFOG almost tripled the size of the initial repertoire of human ORs in a single experiment, and only 7% of the PCR clones had to be sequenced. Extremely high degeneracies, reaching over a billion combinations of distinct PCR primer pairs, proved to be very effective and yielded only 0.4% nonspecific products.

PURPOSE: The assessment of plasma cytokine levels adds a useful tool to the diagnostic measures in severe inflammatory diseases. Proinflammatory cytokine levels in abdominal fluid after abdominal surgery have been shown to far exceed plasma cytokine levels. Thus, we investigated the local release of interleukin 1β, interleukin 6, and tumor necrosis factor-α in patients after colorectal surgery during the early postoperative period to evaluate whether it may serve as an indicator of evolving peritonitis. METHOD: In a prospective, observational pilot study, the first 12 consecutive patients who did not develop any postoperative complications (Group I), and the first 12 patients with secondary peritonitis caused by an anastomotic leakage (Group II), were included in the study. Interleukin 6, interleukin 1β, and tumor necrosis factor-α levels were determined in the abdominal exudate and compared between the groups within the first four days after colorectal surgery. RESULTS: Abdominal fluid interleukin 6 levels in Group II patients were higher (162,500 ± 105,800 pg/ml) as early as the first postoperative day compared with Group I (27,940 ± 13,860 pg/ml; P < 0.0001); this lasted for the whole observation period. The same applies to tumor necrosis factor-α levels (461.4 ± 167.8 pg/ml vs. 175.8 ± 178.6 pg/ml on day 1; P = 0.0007). The difference in interleukin 1β cytokine levels became statistically significant on the third postoperative day. Moreover, abdominal fluid cytokine levels rose in Group II, whereas they remained virtually unchanged or even tended to decrease over time in Group I. CONCLUSION: We suggest that the estimation of the peritoneal cytokine levels might be an additional diagnostic tool that can support the early recognition of peritonitic complications in colorectal surgery.

Objectives The therapeutic scheme for non-small cell lung cancer (NSCLC) patients can be improved if adapted to the individual response. For example, 60–70% of adenocarcinoma patients show response to EGFR-tyrosine kinase inhibitors in the presence of mutated EGFR. We searched for additional target molecules involved in the action of the EGFR-tyrosine kinase inhibitor erlotinib in the absence of EGFR mutations, which might be suitable for combinatorial therapy approaches. Materials and Methods Erlotinib-response associated proteins were investigated in patient-derived NSCLC mouse xenografts by reverse-phase protein array technology (RPPA) and Western blotting. A combinatorial treatment approach was carried out in NSCLC cell lines and H1299 mouse xenografts, and subsequently analyzed for consequences in cell growth and signal transduction. Results AMP-activated protein kinase (AMPK) expression was increased in erlotinib responders before and after treatment. In a combinatorial approach, activation of AMPK by A-769662 and erlotinib treatment showed a synergistic effect in cell growth reduction and apoptosis activation in H1299 cells compared to the single drugs. AMPK pathway analyses revealed an effective inhibition of mTOR signaling by drug combination. In H1299 xenografts, the tumor size was significantly decreased after combinatorial treatment. Conclusion Our results suggest that AMPK activation status affects response to erlotinib in distinct lung tumor models.

The computational prediction of alternative splicing from high-throughput sequencing data is inherently difficult and necessitates robust statistical measures because the differential splicing signal is overlaid by influencing factors such as gene expression differences and simultaneous expression of multiple isoforms amongst others. In this work we describe ARH-seq, a discovery tool for differential splicing in case–control studies that is based on the information-theoretic concept of entropy. ARH-seq works on high-throughput sequencing data and is an extension of the ARH method that was originally developed for exon microarrays. We show that the method has inherent features, such as independence of transcript exon number and independence of differential expression, what makes it particularly suited for detecting alternative splicing events from sequencing data. In order to test and validate our workflow we challenged it with publicly available sequencing data derived from human tissues and conducted a comparison with eight alternative computational methods. In order to judge the performance of the different methods we constructed a benchmark data set of true positive splicing events across different tissues agglomerated from public databases and show that ARH-seq is an accurate, computationally fast and high-performing method for detecting differential splicing events.

Actinobacillus (A.) pleuropneumoniae is among the most important pathogens in pig. The agent causes severe economic losses due to decreased performance, the occurrence of acute or chronic pleuropneumonia, and an increase in death incidence. Since therapeutics cannot be used in a sustainable manner, and vaccination is not always available, new prophylactic measures are urgently needed. Recent research has provided evidence for a genetic predisposition in susceptibility to A. pleuropneumoniae in a Hampshire × German Landrace F2 family with 170 animals. The aim of the present study is to characterize the expression response in this family in order to unravel resistance and susceptibility mechanisms and to prioritize candidate genes for future fine mapping approaches. F2 pigs differed distinctly in clinical, pathological, and microbiological parameters after challenge with A. pleuropneumoniae. We monitored genome-wide gene expression from the 50 most and 50 least susceptible F2 pigs and identified 171 genes differentially expressed between these extreme phenotypes. We combined expression QTL analyses with network analyses and functional characterization using gene set enrichment analysis and identified a functional hotspot on SSC13, including 55 eQTL. The integration of the different results provides a resource for candidate prioritization for fine mapping strategies, such as TF, TFRC, RUNX1, TCN1, HP, CD14, among others.

Motivation: DNA enrichment followed by sequencing is a versatile tool in molecular biology, with a wide variety of applications including genome-wide analysis of epigenetic marks and mechanisms. A common requirement of these diverse applications is a comparison of read coverage between experimental conditions. The amount of samples generated for such comparisons ranges from few replicates to hundreds of samples per condition for epigenome-wide association studies. Consequently, there is an urgent need for software that allows for fast and simple processing and comparison of sequencing data derived from enriched DNA. Results: Here, we present a major update of the R/Bioconductor package MEDIPS, which allows for an arbitrary number of replicates per group and integrates sophisticated statistical methods for the detection of differential coverage between experimental conditions. Our approach can be applied to a diversity of quantitative sequencing data. In addition, our update adds novel functionality to MEDIPS, including correlation analysis between samples, and takes advantage of Bioconductor’s annotation databases to facilitate annotation of specific genomic regions. Availability and implementation: The latest version of MEDIPS is available as version 1.12.0 and part of Bioconductor 2.13. The package comes with a manual containing detailed description of its functionality and is available at http://www.bioconductor.org. Contact: lienhard@molgen.mpg.de Supplementary information: Supplementary data are available at Bioinformatics online.

Nodal signaling controls germ layer formation, left-right asymmetry, and patterning of the brain in the vertebrate embryo. Cellular responses to Nodal signals are complex and include changes in gene expression, cell morphology, and migratory behavior. Only little is known about the genes regulated by Nodal signaling. We designed a subtractive screening strategy by using a constitutively active Nodal receptor to identify putative target genes of Nodal signals in the early gastrula of zebrafish embryos. By quantitative analysis of macro-array hybridizations, 132 genes corresponding to 1.4% of genes on the entire macro-array were identified, which were enriched in the Nodal-induced probe pool. These genes encode components of signal transduction pathways, transcription regulators, proteins involved in protein metabolism but also cytoskeletal components and metabolic enzymes, suggesting dramatic changes of cell physiology in gastrula cells in response to Nodal signals. © 2001 Wiley-Liss, Inc.

The zebrafish is a powerful system for understanding the vertebrate genome, allowing the combination of genetic, molecular, and embryological analysis. Expressed sequence tags (ESTs) provide a rapid means of identifying an organism's genes for further analysis, but any EST project is limited by the availability of suitable libraries. Such cDNA libraries must be of high quality and provide a high rate of gene discovery. However, commonly used normalization and subtraction procedures tend to select for shorter, truncated, and internally primed inserts, seriously affecting library quality. An alternative procedure is to use oligonucleotide fingerprinting (OFP) to precluster clones before EST sequencing, thereby reducing the re-sequencing of common transcripts. Here, we describe the use of OFP to normalize and subtract 75,000 clones from two cDNA libraries, to a minimal set of 25,102 clones. We generated 25,788 ESTs (11,380 3′ and 14,408 5′) from over 16,000 of these clones. Clustering of 10,654 high-quality 3′ ESTs from this set identified 7232 clusters (likely genes), corresponding to a 68% gene diversity rate, comparable to what has been reported for the best normalized human cDNA libraries, and indicating that the complete set of 25,102 clones contains as many as 17,000 genes. Yet, the library quality remains high. The complete set of 25,102 clones is available for researchers as glycerol stocks, filters sets, and as individual EST clones. These resources have been used for radiation hybrid, genetic, and physical mapping of the zebrafish genome, as well as positional cloning and candidate gene identification, molecular marker, and microarray development.[The sequence data described in this paper have been submitted to the dbEST/GenBank data library under accession nos. AA497144–AA497369, AA542435–AA542678,AA545709–AA545724, AI384176–AI384205, AI384761–AI384796,AI396646–AI396663, AI396733–AI396777, AI396895–AI396938,AI397015–AI397130, AI397219–AI397252, AI397388–AI397484,AI415743–AI416403, AI436854–AI437493, AI444118–AI444540,AI461280–AI461395, AI476823–AI478024, AI496677–AI497576,AI522337–AI522810, AI544445–AI546083, AI558267–AI558995,AI584192–AI585023, AI585025–AI585238, AI588088–AI588836,AI601277–AI601868, AI626134–AI626875, AI629052–AI629398,AI641018–AI641780, AI657549–AI658347, AI666929–AI667197,AI667264–AI667414, AI667488–AI667567, AI721460–AI721747,AI721839–AI721978, and AI722283–AI722483.]

Clustering is one of the main mathematical challenges in large-scale gene expression analysis. We describe a clustering procedure based on a sequential k-means algorithm with additional refinements that is able to handle high-throughput data in the order of hundreds of thousands of data items measured on hundreds of variables. The practical motivation for our algorithm is oligonucleotide fingerprinting—a method for simultaneous determination of expression level for every active gene of a specific tissue—although the algorithm can be applied as well to other large-scale projects like EST clustering and qualitative clustering of DNA-chip data. As a pairwise similarity measure between two p-dimensional data points,x and y, we introduce mutual information that can be interpreted as the amount of information about x iny, and vice versa. We show that for our purposes this measure is superior to commonly used metric distances, for example, Euclidean distance. We also introduce a modified version of mutual information as a novel method for validating clustering results when the true clustering is known. The performance of our algorithm with respect to experimental noise is shown by extensive simulation studies. The algorithm is tested on a subset of 2029 cDNA clones coming from 15 different genes from a cDNA library derived from human dendritic cells. Furthermore, the clustering of these 2029 cDNA clones is demonstrated when the entire set of 76,032 cDNA clones is processed.

We describe the use of oligonucleotide fingerprinting for the generation of a normalized cDNA library from unfertilized sea urchin eggs and report the preliminary analysis of this library, which resulted in the establishment of a partial gene catalogue of the sea urchin egg. In an analysis of 21,925 cDNA clones by hybridization with 217 oligonucleotide probes, we were able to identify 6291 clusters corresponding to different transcripts, ranging in size from 1 to 265 clones. This corresponds to an average 3.5-fold normalization of the starting library. The normalized library represents about one-third of all genes expressed in the sea urchin egg. To generate sequence information for the transcripts represented by the clusters, representative clones selected from 711 clusters were sequenced. The construction and preliminary analysis of the normalized library are the first steps in the assembly of an increasingly complete collection of maternal genes expressed in the sea urchin egg, which will provide a number of insights into the early development of this well-characterized model organism.

Colorectal carcinoma is a common malignant disease with a high mortality rate. It arises most frequently in adenomas of the colorectum with different grades of dysplasia. Endoscopy and biopsy are among the most reliable diagnostic tools currently available. Diagnosis of malignancy at an early stage is sometimes difficult. This study reports on a new method, "immunoscopy", that combines endoscopy and immunofluorescent diagnostic procedures; it is the first reported use of locally applied fluorescein-labeled antibodies for detection of colorectal carcinomas.A monoclonal antibody against carcinoembryonic antigen was fluorescein labeled. In phase I, formalin-fixed tissue samples, and in phase II, postoperative fresh tissue samples from tumorous and nontumorous areas of resected colorectal carcinomas were studied. After being incubated for 10 minutes, specific fluorescence was visualized with a conventional endoscope whose range was increased by means of two narrow band filters.Because of high levels of autofluorescence, evaluation of immunoscopic investigations using formalin fixed tissue (phase I) was not carried out. Immunoscopic investigation with postoperative fresh tissue samples could differentiate between tumorous and nontumorous areas (p0.001). Immunoscopic results were compared with data from immunohistochemical investigations with anti-carcinoembryonic antigen on the same tissue samples.Immunoscopy can differentiate between malignant and benign mucosal areas in fresh tissue samples. The high sensitivity of immunoscopy could potentially make it a useful diagnostic complement to routine endoscopy.

The group of dinucleoside polyphosphates encompasses a large number of molecules consisting of two nucleosides which are connected by a phosphate chain of variable length. While the receptors activated by dinucleoside polyphosphates as well as their degradation have been studied in detail, its biosynthesis has not been elucidated so far. Since endothelial cells released the dinucleoside polyphosphate uridine adenosine tetraphosphate (Up4A), we tested cytosolic proteins of human endothelial cells obtained from dermal vessels elicited for enzymatic activity. When incubated with ADP and UDP, these cells showed increasing concentrations of Up4A. The underlying enzyme was isolated by chromatography and the mass spectrometric analysis revealed that the enzymatic activity was caused by the vascular endothelial growth factor receptor 2 (VEGFR2). Since VEGFR2 but neither VEGFR1 nor VEGFR3 were capable to synthesise dinucleoside polyphosphates, Tyr-1175 of VEGFR2 is most likely essential for the enzymatic activity of interest. Further, VEGFR2-containing cells like HepG2, THP-1 and RAW264.7 were capable of synthesising dinucleoside polyphosphates. VEGFR2-transfected HEK 293T/17 but not native HEK 293T/17 cells synthesised dinucleoside polyphosphates in vivo too. The simultaneous biosynthesis of dinucleoside polyphosphates could amplify the response to VEGF, since dinucleoside polyphosphates induce cellular growth via P2Y purinergic receptors. Thus the biosynthesis of dinucleoside polyphosphates by VEGFR2 may enhance the proliferative response to VEGF. Given that VEGFR2 is primarily expressed in endothelial cells, the biosynthesis of dinucleoside polyphosphates is mainly located in the vascular system. Since the vasculature is also the main site of action of dinucleoside polyphosphates, activating vascular purinoceptors, blood vessels appear as an autocrine system with respect to dinucleoside polyphosphates. We conclude that VEGFR2 receptor is capable of synthesising dinucleoside polyphosphates. These mediators may modulate the effects of VEGFR2 due to their proliferative effects.