Search

Summary Productive splicing of human precursor messenger RNAs (pre-mRNAs) requires the correct selection of authentic splice sites (SS) from the large pool of potential SS. Although SS consensus sequence and splicing regulatory proteins are known to influence SS usage, the mechanisms ensuring the effective suppression of cryptic SS are insufficiently explored. Here, we find that many aberrant exonic SS are efficiently silenced by the exon junction complex (EJC), a multi-protein complex that is deposited on spliced mRNA near the exon-exon junction. Upon depletion of EJC proteins, cryptic SS are de-repressed, leading to the mis-splicing of a broad set of mRNAs. Mechanistically, the EJC-mediated recruitment of the splicing regulator RNPS1 inhibits cryptic 5′SS usage, while the deposition of the EJC core directly masks reconstituted 3′SS, thereby precluding transcript disintegration. Thus, the EJC protects the transcriptome of mammalian cells from inadvertent loss of exonic sequences and safeguards the expression of intact, full-length mRNAs.

The phylum Nematoda is biologically diverse; it includes parasites of plants and animals in addition to free-living taxa. To date, the genomes of six nematodes have been sequenced. Comparative analyses of these ecologically diverse nematodes are beginning to reveal the mechanisms by which parasites arise and how they evolve. Here, we discuss some emerging principles for the mechanisms and evolution of parasitism. First, horizontal gene transfer represents a common theme in nematode parasites. Second, the human parasite Brugia malayi lost otherwise essential genes most probably owing to the mutualistic relationship with a bacterial endosymbiont. Finally, some parasitic features evolved under free-living conditions. A recent study revealed a conserved endocrine mechanism controlling the formation of dauer and infective larvae in nematodes.

Runx2 is an essential factor for skeletogenesis and heterozygous loss causes cleidocranial dysplasia in humans and a corresponding phenotype in the mouse. Homozygous Runx2-deficient mice lack hypertrophic cartilage and bone. We compared the expression profiles of E14.5 wildtype and Runx2−/− murine embryonal humeri to identify new transcripts potentially involved in cartilage and bone development. Seventy-one differentially expressed genes were identified by two independent oligonucleotide-microarray hybridizations and quantitative RT-PCR experiments. Gene Ontology analysis demonstrated an enrichment of the differentially regulated genes in annotations to terms such as extracellular, skeletal development, and ossification. In situ hybridization on E15.5 limb sections was performed for all 71 differentially regulated genes. For 54 genes conclusive in situ hybridization results were obtained and all of them showed skeletal expression. Co-expression with Runx2 was demonstrated for 44 genes. While 41 of the 71 differentially expressed genes have a known role in bone and cartilage, we identified 21 known genes that have not yet been implicated in skeletal development and 9 entirely new transcripts. Expression in the developing skeleton was demonstrated for 21 of these genes.

Germline removal provokes longevity in several species and shifts resources towards survival and repair. Several Caenorhabditis elegans transcription factors regulate longevity arising from germline removal; yet, how they work together is unknown. Here we identify a Myc-like HLH transcription factor network comprised of Mondo/Max-like complex (MML-1/MXL-2) to be required for longevity induced by germline removal, as well as by reduced TOR, insulin/IGF signalling and mitochondrial function. Germline removal increases MML-1 nuclear accumulation and activity. Surprisingly, MML-1 regulates nuclear localization and activity of HLH-30/TFEB, a convergent regulator of autophagy, lysosome biogenesis and longevity, by downregulating TOR signalling via LARS-1/leucyl-transfer RNA synthase. HLH-30 also upregulates MML-1 upon germline removal. Mammalian MondoA/B and TFEB show similar mutual regulation. MML-1/MXL-2 and HLH-30 transcriptomes show both shared and preferential outputs including MDL-1/MAD-like HLH factor required for longevity. These studies reveal how an extensive interdependent HLH transcription factor network distributes responsibility and mutually enforces states geared towards reproduction or survival., Removal of the C. elegans germline substantially increases organismal lifespan. Here, Nakamura et al. show that the transcription factors MML-1 and MXL-2 coordinate this process in that they reduce TOR signalling and increase autophagy by regulating activity of HLH-30.

Background The development of multicellular organisms is accompanied by gene expression changes in differentiating cells. Profiling stage-specific expression during development may reveal important insights into gene sets that contributed to the morphological diversity across the animal kingdom. Results We sequenced RNA-seq libraries throughout a developmental timecourse of the nematode Pristionchus pacificus. The transcriptomes reflect early larval stages, adult worms including late larvae, and growth-arrested dauer larvae and allowed the identification of developmentally regulated gene clusters. Our data reveals similar trends as previous transcriptome profiling of dauer worms and represents the first expression data for early larvae in P. pacificus. Gene expression clusters characterizing early larval stages show most significant enrichments of chaperones, while collagens are most significantly enriched in transcriptomes of late larvae and adult worms. By combining expression data with phylogenetic analysis, we found that developmentally regulated genes are found in paralogous clusters that have arisen through lineage-specific duplications after the split from the Caenorhabditis elegans branch. Conclusions We propose that gene duplications of developmentally regulated genes represent a plausible evolutionary mechanism to increase the dosage of stage-specific expression. Consequently, this may contribute to the substantial divergence in expression profiles that has been observed across larger evolutionary time scales. Electronic supplementary material The online version of this article (doi:10.1186/s12862-015-0466-2) contains supplementary material, which is available to authorized users.

Growing evidence suggests a key role for RNA binding proteins (RBPs) in genome stability programs. Additionally, recent developments in RNA sequencing technologies, as well as mass-spectrometry techniques, have greatly expanded our knowledge on protein-RNA interactions. We here use full transcriptome sequencing and label-free LC/MS/MS to identify global changes in protein-RNA interactions in response to etoposide-induced genotoxic stress. We show that RBPs have distinct binding patterns in response to genotoxic stress and that inactivation of the RBP regulator module, p38/MK2, can affect the entire spectrum of protein-RNA interactions that take place in response to stress. In addition to validating the role of known RBPs like Srsf1, Srsf2, Elavl1 in the genotoxic stress response, we add a new collection of RBPs to the DNA damage response. We identify Khsrp as a highly regulated RBP in response to genotoxic stress and further validate its role as a driver of the G(1/)S transition through the suppression of Cdkn1a(P21) transcripts. Finally, we identify KHSRP as an indicator of overall survival, as well as disease free survival in glioblastoma multiforme.

Healthy placental development is essential for reproductive success; failure of the feto-maternal interface results in pre-eclampsia and intrauterine growth retardation. We found that grainyhead-like 2 (GRHL2), a CP2-type transcription factor, is highly expressed in chorionic trophoblast cells, including basal chorionic trophoblast (BCT) cells located at the chorioallantoic interface in murine placentas. Placentas from Grhl2-deficient mouse embryos displayed defects in BCT cell polarity and basement membrane integrity at the chorioallantoic interface, as well as a severe disruption of labyrinth branching morphogenesis. Selective Grhl2 inactivation only in epiblast-derived cells rescued all placental defects but phenocopied intraembryonic defects observed in global Grhl2 deficiency, implying the importance of Grhl2 activity in trophectoderm-derived cells. ChIP-seq identified 5282 GRHL2 binding sites in placental tissue. By integrating these data with placental gene expression profiles, we identified direct and indirect Grhl2 targets and found a marked enrichment of GRHL2 binding adjacent to genes downregulated in Grhl2−/− placentas, which encoded known regulators of placental development and epithelial morphogenesis. These genes included that encoding the serine protease inhibitor Kunitz type 1 (Spint1), which regulates BCT cell integrity and labyrinth formation. In human placenta, we found that human orthologs of murine GRHL2 and its targets displayed co-regulation and were expressed in trophoblast cells in a similar domain as in mouse placenta. Our data indicate that a conserved Grhl2-coordinated gene network controls trophoblast branching morphogenesis, thereby facilitating development of the site of feto-maternal exchange. This might have implications for syndromes related to placental dysfunction.

The E3 ubiquitin ligase Ube3a is an important regulator of activity-dependent synapse development and plasticity. Ube3a mutations cause Angelman syndrome and have been associated with autism spectrum disorders (ASD). However, the biological significance of alternative Ube3a transcripts generated in mammalian neurons remains unknown. We report here that Ube3a1 RNA, a transcript that encodes a truncated Ube3a protein lacking catalytic activity, prevents exuberant dendrite growth and promotes spine maturation in rat hippocampal neurons. Surprisingly, Ube3a1 RNA function was independent of its coding sequence but instead required a unique 3' untranslated region and an intact microRNA pathway. Ube3a1 RNA knockdown increased activity of the plasticity-regulating miR-134, suggesting that Ube3a1 RNA acts as a dendritic competing endogenous RNA. Accordingly, the dendrite-growth-promoting effect of Ube3a1 RNA knockdown in vivo is abolished in mice lacking miR-134. Taken together, our results define a noncoding function of an alternative Ube3a transcript in dendritic protein synthesis, with potential implications for Angelman syndrome and ASD.

Stem cells reside in specialized niches that are critical for their function. Quiescent hair follicle stem cells (HFSCs) are confined within the bulge niche, but how the molecular composition of the niche regulates stem cell behaviour is poorly understood. Here we show that integrin-linked kinase (ILK) is a key regulator of the bulge extracellular matrix microenvironment, thereby governing the activation and maintenance of HFSCs. ILK mediates deposition of inverse laminin (LN)-332 and LN-511 gradients within the basement membrane (BM) wrapping the hair follicles. The precise BM composition tunes activities of Wnt and transforming growth factor-β pathways and subsequently regulates HFSC activation. Notably, reconstituting an optimal LN microenvironment restores the altered signalling in ILK-deficient cells. Aberrant stem cell activation in ILK-deficient epidermis leads to increased replicative stress, predisposing the tissue to carcinogenesis. Overall, our findings uncover a critical role for the BM niche in regulating stem cell activation and thereby skin homeostasis., Integrin-linked kinase (ILK) is known to modulate the extracellular matrix and hair follicle morphogenesis. Here, Morgner et al. show that lack of ILK causes an aberrant ratio of basement membrane laminins, activating stem cells and predisposing skin to carcinogenesis.

Genetic analysis using experimentally induced mutations has been a most valuable tool in the analysis of various organisms. However, genetic analysis of endoparasitic organisms tends to be difficult because of the limited accessibility of the sexually reproducing adults, which are normally located within the host. Nematodes of the genera Strogyloides and Parastrongyloides represent an exception to this because they can form facultative free-living sexually reproducing generations in between parasitic generations. Here we present a protocol for the chemical mutagenesis of Strongyloides ratti. Further we evaluate the feasibility of identifying the induced mutations by whole genome re-sequencing.

RhoA signaling regulates the activity of the transcription factor SRF (serum response factor) during muscle differentiation. How RhoA signaling is integrated at SRF target promoters to achieve muscle-lineage-specific expression is largely unknown. Using large-scale expression profiling combined with bioinformatic and biochemical approaches, we identified several SRF target genes, including Fhl2, encoding a transcriptional cofactor that is highly expressed in the heart. SRF binds the Fhl2 promoter in vivo and regulates Fhl2 expression in response to RhoA activation. FHL2 protein and SRF interact physically, and FHL2 binds the promoters of SRF-responsive smooth muscle (SM) genes, but not the promoters of immediate-early genes (IEGs), in response to RhoA. FHL2 antagonizes induction of SM genes, but not IEGs or cardiac genes, by competing with the coactivator MAL/MRTF-A for SRF binding. Our findings identify an autoregulatory mechanism to selectively regulate subsets of RhoA-activated SRF target genes.

The transcription factor ZEB1 has gained attention in tumor biology of epithelial cancers because of its function in epithelial-mesenchymal transition, DNA repair, stem cell biology and tumor-induced immunosuppression, but its role in gliomas with respect to invasion and prognostic value is controversial. We characterized ZEB1 expression at single cell level in 266 primary brain tumors and present a comprehensive dataset of high grade gliomas with Ki67, p53, IDH1, and EGFR immunohistochemistry, as well as EGFR FISH. ZEB1 protein expression in glioma stem cell lines was compared to their parental tumors with respect to gene expression subtypes based on RNA-seq transcriptomic profiles. ZEB1 is widely expressed in glial tumors, but in a highly variable fraction of cells. In glioblastoma, ZEB1 labeling index is higher in tumors with EGFR amplification or IDH1 mutation. Co-labeling studies showed that tumor cells and reactive astroglia, but not immune cells contribute to the ZEB1 positive population. In contrast, glioma cell lines constitutively express ZEB1 irrespective of gene expression subtype. In conclusion, our data indicate that immune infiltration likely contributes to differential labelling of ZEB1 and confounds interpretation of bulk ZEB1 expression data. publishedVersion

Discovery of trans-splicing in multiple metazoan lineages led to the identification of operon-like gene organization in diverse organisms, including trypanosomes, tunicates, and nematodes, but the functional significance of such operons is not completely understood. To see whether the content or organization of operons serves similar roles across species, we experimentally defined operons in the nematode model Pristionchus pacificus. We performed affinity capture experiments on mRNA pools to specifically enrich for transcripts that are trans-spliced to either the SL1- or SL2-spliced leader, using spliced leader–specific probes. We obtained distinct trans-splicing patterns from the analysis of three mRNA pools (total mRNA, SL1 and SL2 fraction) by RNA-seq. This information was combined with a genome-wide analysis of gene orientation and spacing. We could confirm 2219 operons by RNA-seq data out of 6709 candidate operons, which were predicted by sequence information alone. Our gene order comparison of the Caenorhabditis elegans and P. pacificus genomes shows major changes in operon organization in the two species. Notably, only 128 out of 1288 operons in C. elegans are conserved in P. pacificus. However, analysis of gene-expression profiles identified conserved functions such as an enrichment of germline-expressed genes and higher expression levels of operonic genes during recovery from dauer arrest in both species. These results provide support for the model that a necessity for increased transcriptional efficiency in the context of certain developmental processes could be a selective constraint for operon evolution in metazoans. Our method is generally applicable to other metazoans to see if similar functional constraints regulate gene organization into operons.

Background The phylum Apicomplexa comprises important unicellular human parasites such as Toxoplasma and Plasmodium. Eimeria is the largest and most diverse genus of apicomplexan parasites and some species of the genus are the causative agent of coccidiosis, a disease economically devastating in poultry. We report a complete genome sequence of the mouse parasite Eimeria falciformis. We assembled and annotated the genome sequence to study host-parasite interactions in this understudied genus in a model organism host. Results The genome of E. falciformis is 44 Mb in size and contains 5,879 predicted protein coding genes. Comparative analysis of E. falciformis with Toxoplasma gondii shows an emergence and diversification of gene families associated with motility and invasion mainly at the level of the Coccidia. Many rhoptry kinases, among them important virulence factors in T. gondii, are absent from the E. falciformis genome. Surface antigens are divergent between Eimeria species. Comparisons with T. gondii showed differences between genes involved in metabolism, N-glycan and GPI-anchor synthesis. E. falciformis possesses a reduced set of transmembrane transporters and we suggest an altered mode of iron uptake in the genus Eimeria. Conclusions Reduced diversity of genes required for host-parasite interaction and transmembrane transport allow hypotheses on host adaptation and specialization of a single host parasite. The E. falciformis genome sequence sheds light on the evolution of the Coccidia and helps to identify determinants of host-parasite interaction critical for drug and vaccine development. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-696) contains supplementary material, which is available to authorized users.

BACKGROUND: RNA-binding proteins (RBPs) mediate mRNA biogenesis, translation and decay. We recently developed an approach to profile transcriptome-wide RBP contacts on polyadenylated transcripts by next-generation sequencing. A comparison of such profiles from different biological conditions has the power to unravel dynamic changes in protein-contacted cis-regulatory mRNAs regions without a priori knowledge of the regulatory protein component. RESULTS: We compared protein occupancy profiles of polyadenylated transcripts in MCF7 and HEK293 cells. Briefly, we developed a bioinformatics workflow to identify differential crosslinking sites in cDNA reads of 4-thiouridine crosslinked polyadenylated RNA samples. We identified 30,000 differential crosslinking sites between MCF7 and HEK293 cells at an estimated false-discovery rate of 10%. 73% of all reported differential protein-RNA contact sites cannot be explained by local changes in exon usage as indicated by complementary RNA-seq data. The majority of differentially crosslinked positions are located in 3[prime] UTRs, show distinct secondary-structure characteristics and overlap with binding sites of known RBPs, such as ELAVL1. Importantly, mRNA transcripts with the most significant occupancy changes show elongated mRNA half-lives in MCF7 cells. CONCLUSIONS: We present a global comparison of protein occupancy profiles from different cell types, and provide evidence for altered mRNA metabolism as a result of differential protein-RNA contacts. Additionally, we introduce POPPI, a bioinformatics workflow for the analysis of protein occupancy profiling experiments. Our work demonstrates the value of protein occupancy profiling for assessing cis-regulatory RNA sequence space and its dynamics in growth, development and disease.

Studies of X chromosome evolution in various organisms have indicated that sex-biased genes are nonrandomly distributed between the X and autosomes. Here, to extend these studies to nematodes, we annotated and analyzed X chromosome gene content in four Caenorhabditis species and in Pristionchus pacificus. Our gene expression analyses comparing young adult male and female mRNA-seq data indicate that, in general, nematode X chromosomes are enriched for genes with high female-biased expression and depleted of genes with high male-biased expression. Genes with low sex-biased expression do not show the same trend of X chromosome enrichment and depletion. Combined with the observation that highly sex-biased genes are primarily expressed in the gonad, differential distribution of sex-biased genes reflects differences in evolutionary pressures linked to tissue-specific regulation of X chromosome transcription. Our data also indicate that X dosage imbalance between males (XO) and females (XX) is influential in shaping both expression and gene content of the X chromosome. Predicted upregulation of the single male X to match autosomal transcription (Ohno’s hypothesis) is supported by our observation that overall transcript levels from the X and autosomes are similar for highly expressed genes. However, comparison of differentially located one-to-one orthologs between C. elegans and P. pacificus indicates lower expression of X-linked orthologs, arguing against X upregulation. These contradicting observations may be reconciled if X upregulation is not a global mechanism but instead acts locally on a subset of tissues and X-linked genes that are dosage sensitive.

Animal development is complex yet surprisingly robust. Animals may develop alternative phenotypes conditional on environmental changes. Under unfavorable conditions, Caenorhabditis elegans larvae enter the dauer stage, a developmentally arrested, long-lived, and stress-resistant state. Dauer larvae of free-living nematodes and infective larvae of parasitic nematodes share many traits including a conserved endocrine signaling module (DA/DAF-12), which is essential for the formation of dauer and infective larvae. We speculated that conserved post-transcriptional regulatory mechanism might also be involved in executing the dauer and infective larvae fate. We used an unbiased sequencing strategy to characterize the microRNA (miRNA) gene complement in C. elegans, Pristionchus pacificus, and Strongyloides ratti. Our study raised the number of described miRNA genes to 257 for C. elegans, tripled the known gene set for P. pacificus to 362 miRNAs, and is the first to describe miRNAs in a Strongyloides parasite. Moreover, we found a limited core set of 24 conserved miRNA families in all three species. Interestingly, our estimated expression fold changes between dauer versus nondauer stages and infective larvae versus free-living stages reveal that despite the speed of miRNA gene set evolution in nematodes, homologous gene families with conserved "dauer-infective" expression signatures are present. These findings suggest that common post-transcriptional regulatory mechanisms are at work and that the same miRNA families play important roles in developmental arrest and long-term survival in free-living and parasitic nematodes.

Background An organism can respond to changing environmental conditions by adjusting gene regulation and by forming alternative phenotypes. In nematodes, these mechanisms are coupled because many species will form dauer larvae, a stress-resistant and non-aging developmental stage, when exposed to unfavorable environmental conditions, and execute gene expression programs that have been selected for the survival of the animal in the wild. These dauer larvae represent an environmentally induced, homologous developmental stage across many nematode species, sharing conserved morphological and physiological properties. Hence it can be expected that some core components of the associated transcriptional program would be conserved across species, while others might diverge over the course of evolution. However, transcriptional and metabolic analysis of dauer development has been largely restricted to Caenorhabditis elegans. Here, we use a transcriptomic approach to compare the dauer stage in the evolutionary model system Pristionchus pacificus with the dauer stage in C. elegans. Results We have employed Agilent microarrays, which represent 20,446 P. pacificus and 20,143 C. elegans genes to show an unexpected divergence in the expression profiles of these two nematodes in dauer and dauer exit samples. P. pacificus and C. elegans differ in the dynamics and function of genes that are differentially expressed. We find that only a small number of orthologous gene pairs show similar expression pattern in the dauers of the two species, while the non-orthologous fraction of genes is a major contributor to the active transcriptome in dauers. Interestingly, many of the genes acquired by horizontal gene transfer and orphan genes in P. pacificus, are differentially expressed suggesting that these genes are of evolutionary and functional importance. Conclusion Our data set provides a catalog for future functional investigations and indicates novel insight into evolutionary mechanisms. We discuss the limited conservation of core developmental and transcriptional programs as a common aspect of animal evolution.

The identification of cis-acting elements on DNA is crucial for the understanding of the complex regulatory networks that govern many cell mechanisms. However, this task is very complex since it is estimated that there are 1500 different transcription factors (TFs) in the human genome, each of which can bind to multiple loci directly or indirectly. The standard computational approach is the use of a position weight matrix (PWM) to represent the binding preference of a transcription factor and the use of statistical procedures to detect genomic regions with high binding scores. Given the small and degenerate signals of most PWMs, such approach suffers from a very high number of false positive hits. Current research has proven that genome wide assays reflecting open chromatin, such as DNase digestion or histone modifications, can improve sequence based detection of the binding location of transcription factors that are active in a particular cell type. We propose here a Multivariate Hidden Markov Model that is able to improve the prediction of transcription factor binding locations by integrating DNase digestion and histone modification data. Our methodology improves sensitivity, in comparison to existing methods, with little or no effect at specificity rates. This study shows that it is possible to improve predictability power of cis-acting elements by correctly integrating DNase and histone modification data, allowing for more sophisticated studies using a larger set of epigenetic signals.

In animals, RNA binding proteins (RBPs) and microRNAs (miRNAs) post-transcriptionally regulate the expression of virtually all genes by binding to RNA. Recent advances in experimental and computational methods facilitate transcriptome-wide mapping of these interactions. It is thought that the combinatorial action of RBPs and miRNAs on target mRNAs form a post-transcriptional regulatory code. We provide a database that supports the quest for deciphering this regulatory code. Within doRiNA, we are systematically curating, storing and integrating binding site data for RBPs and miRNAs. Users are free to take a target (mRNA) or regulator (RBP and/or miRNA) centric view on the data. We have implemented a database framework with short query response times for complex searches (e.g. asking for all targets of a particular combination of regulators). All search results can be browsed, inspected and analyzed in conjunction with a huge selection of other genome-wide data, because our database is directly linked to a local copy of the UCSC genome browser. At the time of writing, doRiNA encompasses RBP data for the human, mouse and worm genomes. For computational miRNA target site predictions, we provide an update of PicTar predictions.

Freshwater planaria are a very attractive model system for stem cell biology, tissue homeostasis, and regeneration. The genome of the planarian Schmidtea mediterranea has recently been sequenced and is estimated to contain >20,000 protein-encoding genes. However, the characterization of its transcriptome is far from complete. Furthermore, not a single proteome of the entire phylum has been assayed on a genome-wide level. We devised an efficient sequencing strategy that allowed us to de novo assemble a major fraction of the S. mediterranea transcriptome. We then used independent assays and massive shotgun proteomics to validate the authenticity of transcripts. In total, our de novo assembly yielded 18,619 candidate transcripts with a mean length of 1118 nt after filtering. A total of 17,564 candidate transcripts could be mapped to 15,284 distinct loci on the current genome reference sequence. RACE confirmed complete or almost complete 5′ and 3′ ends for 22/24 transcripts. The frequencies of frame shifts, fusion, and fission events in the assembled transcripts were computationally estimated to be 4.2%–13%, 0%–3.7%, and 2.6%, respectively. Our shotgun proteomics produced 16,135 distinct peptides that validated 4200 transcripts (FDR ≤1%). The catalog of transcripts assembled in this study, together with the identified peptides, dramatically expands and refines planarian gene annotation, demonstrated by validation of several previously unknown transcripts with stem cell-dependent expression patterns. In addition, our robust transcriptome characterization pipeline could be applied to other organisms without genome assembly. All of our data, including homology annotation, are freely available at SmedGD, the S. mediterranea genome database.

Pristionchus pacificus is a nematode model organism whose genome has recently been sequenced. To refine the genome annotation we performed transcriptome and proteome analysis and gathered comprehensive experimental information on gene expression. Transcriptome analysis on a 454 Life Sciences (Roche) FLX platform generated >700,000 expressed sequence tags (ESTs) from two normalized EST libraries, whereas proteome analysis on an LTQ-Orbitrap mass spectrometer detected >27,000 nonredundant peptide sequences from more than 4000 proteins at sub-parts-per-million (ppm) mass accuracy and a false discovery rate of P. pacificus proteome contains a high proportion of small proteins with no known homologs in other species (“pioneer” proteins). Some of these proteins appear to be products of highly homologous genes, pointing to their common origin. We show that >50% of all pioneer genes are transcribed under standard culture conditions and that pioneer proteins significantly contribute to a unimodal distribution of predicted protein sizes in P. pacificus, which has an unusually low median size of 240 amino acids (26.8 kDa). In contrast, the predicted proteome of Caenorhabditis elegans follows a distinct bimodal protein size distribution, with significant functional differences between small and large protein populations. Combined, these results provide the first catalog of the expressed genome of P. pacificus, refinement of its genome annotation, and the first comparison of related nematode models at the proteome level.

Whole genome gene order evolution in higher eukaryotes was initially considered as a random process. Gene order conservation or conserved synteny was seen as a feature of common descent and did not imply the existence of functional constraints. This view had to be revised in the light of results from sequencing dozens of vertebrate genomes. It became apparent that other factors exist that constrain gene order in some genomic regions over long evolutionary time periods. Outside of these regions, genomes diverge more rapidly in terms of gene content and order. We have developed CYNTENATOR, a progressive gene order alignment software, to identify genomic regions of conserved synteny over a large set of diverging species. CYNTENATOR does not depend on nucleotide-level alignments and a priori homology assignment. Our software implements an improved scoring function that utilizes the underlying phylogeny. In this manuscript, we report on our progressive gene order alignment approach, a and give a comparison to previous software and an analysis of 17 vertebrate genomes for conservation in gene order. CYNTENATOR has a runtime complexity of and a space complexity of with being the gene number in a genome. CYNTENATOR performs as good as state-of-the-art software on simulated pairwise gene order comparisons, but is the only algorithm that works in practice for aligning dozens of vertebrate-sized gene orders. Lineage-specific characterization of gene order across 17 vertebrate genomes revealed mechanisms for maintaining conserved synteny such as enhancers and coregulation by bidirectional promoters. Genes outside conserved synteny blocks show enrichments for genes involved in responses to external stimuli, stimuli such as immunity and olfactory response in primate genome comparisons. We even see significant gene ontology term enrichments for breakpoint regions of ancestral nodes close to the root of the phylogeny. Additionally, our analysis of transposable elements has revealed a significant accumulation of LINE-1 elements in mammalian breakpoint regions. In summary, CYNTENATOR is a flexible and scalable tool for the identification of conserved gene orders across multiple species over long evolutionary distances.

We present a highly accurate gene-prediction system for eukaryotic genomes, called mGene. It combines in an unprecedented manner the flexibility of generalized hidden Markov models (gHMMs) with the predictive power of modern machine learning methods, such as Support Vector Machines (SVMs). Its excellent performance was proved in an objective competition based on the genome of the nematode Caenorhabditis elegans. Considering the average of sensitivity and specificity, the developmental version of mGene exhibited the best prediction performance on nucleotide, exon, and transcript level for ab initio and multiple-genome gene-prediction tasks. The fully developed version shows superior performance in 10 out of 12 evaluation criteria compared with the other participating gene finders, including Fgenesh++ and Augustus. An in-depth analysis of mGene's genome-wide predictions revealed that ≈2200 predicted genes were not contained in the current genome annotation. Testing a subset of 57 of these genes by RT-PCR and sequencing, we confirmed expression for 24 (42%) of them. mGene missed 300 annotated genes, out of which 205 were unconfirmed. RT-PCR testing of 24 of these genes resulted in a success rate of merely 8%. These findings suggest that even the gene catalog of a well-studied organism such as C. elegans can be substantially improved by mGene's predictions. We also provide gene predictions for the four nematodes C. briggsae, C. brenneri, C. japonica, and C. remanei. Comparing the resulting proteomes among these organisms and to the known protein universe, we identified many species-specific gene inventions. In a quality assessment of several available annotations for these genomes, we find that mGene's predictions are most accurate.

SUMMARYNematodes and bacteria are major components of the soil ecosystem. Many nematodes use bacteria for food, whereas others evolved specialized bacterial interactions ranging from mutualism to parasitism. Little is known about the biological mechanisms by which nematode–bacterial interactions are achieved, largely because in the laboratory nematodes are often cultured under artificial conditions. We investigated the bacterial interactions of nematodes from the genus Pristionchus that have a strong association with scarab beetles. Pristionchus has a different feeding strategy than Caenorhabditis and meta-genomic 16S sequence analysis of Pristionchus individuals showed a diversity of living bacteria within the nematode gut and on the nematode cuticle. Twenty-three different bacterial strains were isolated from three Pristionchus–beetle associations and were used to study nematode–bacterial interactions under controlled laboratory conditions. We show a continuum of bacterial interactions from dissemination, to reduction in brood size and nematode mortality caused by bacteria derived from insect hosts. Olfactory discrimination experiments show distinct chemoattraction and fitness profiles of Pristionchus nematodes when exposed to different bacteria. For example, Pristionchus pacificus avoids Serratia marcescenspossibly because of pathogenicity. Also, P. pacificus avoids Bacillus thuringiensis and insect pathogenic bacteria but is resistant to the human pathogens Staphylococcus aureus and Pseudomonas aeruginosa, unlike Caenorhabditis elegans. Pristionchus specifically recognize and respond to bacteria that cause ill health. Bringing the nematode–bacterial interaction into the laboratory allows detailed functional studies, including the genetic manipulation of the interaction in both nematodes and bacteria.

Here we present a draft genome sequence of the nematode Pristionchus pacificus, a species that is associated with beetles and is used as a model system in evolutionary biology. With 169 Mb and 23,500 predicted protein-coding genes, the P. pacificus genome is larger than those of Caenorhabditis elegans and the human parasite Brugia malayi. Compared to C. elegans, the P. pacificus genome has more genes encoding cytochrome P450 enzymes, glucosyltransferases, sulfotransferases and ABC transporters, many of which were experimentally validated. The P. pacificus genome contains genes encoding cellulase and diapausin, and cellulase activity is found in P. pacificus secretions, indicating that cellulases can be found in nematodes beyond plant parasites. The relatively higher number of detoxification and degradation enzymes in P. pacificus is consistent with its necromenic lifestyle and might represent a preadaptation for parasitism. Thus, comparative genomics analysis of three ecologically distinct nematodes offers a unique opportunity to investigate the association between genome structure and lifestyle.

Background Controlling gene expression is fundamental to biological complexity. The nematode Caenorhabditis elegans is an important model for studying principles of gene regulation in multi-cellular organisms. A comprehensive parts list of putative regulatory motifs was yet missing for this model system. In this study, we compile a set of putative regulatory motifs by combining evidence from conservation and expression data. Description We present an unbiased comparative approach to a regulatory motif compendium for Caenorhabditis species. This involves the assembly of a new nematode genome, whole genome alignments and assessment of conserved k-mers counts. Candidate motifs are selected from a set of 9,500 randomly picked genes by three different motif discovery strategies. Motif candidates have to pass a conservation enrichment filter. Motif degeneracy and length are optimized. Retained motif descriptions are evaluated by expression data using a non-parametric test, which assesses expression changes due to the presence/absence of individual motifs. Finally, we also provide condition-specific motif ensembles by conditional tree analysis. Conclusion The nematode genomes align surprisingly well despite high neutral substitution rates. Our pipeline delivers motif sets by three alternative strategies. Each set contains less than 400 motifs, which are significantly conserved and correlated with 214 out of 270 tested gene expression conditions. This motif compendium is an entry point to comprehensive studies on nematode gene regulation. The website: http://corg.eb.tuebingen.mpg.de/CMC has extensive query capabilities, supplements this article and supports the experimental list.

Comparative studies have been of invaluable importance to the understanding of evolutionary biology. The evolution of developmental programs can be studied in nematodes at a single cell resolution given their fixed cell lineage. We have established Pristionchus pacificus as a major satellite organism for evolutionary developmental biology relative to Caenorhabditis elegans, the model nematode. Online genomic information to support studies in this satellite system can be accessed at http://www.pristionchus.org. Our web resource offers diverse content covering genome browsing, genetic and physical maps, similarity searches, a community platform and assembly details. Content will be continuously improved as we annotate the P.pacificus genome, and will be an indispensable resource for P.pacificus genomics.

OBJECTIVE: The transcription factor ZEB1 has gained attention in tumor biology of epithelial cancers because of its suggested roles in epithelial-mesenchymal transition, DNA repair, mediation of stem cell properties and tumor-induced immunosuppression, but its role in CNS tumors, especially glioblastoma multiforme is not well understood. METHODS: We comprehensively characterized ZEB1 expression at single cell level by automated classification of 1.6 million nuclei in biopsy samples from 273 human astroglial tumors and correlated its expression with EGFR status, IDH1 mutation and Ki67, CD44, GFAP and CD68 immunohistochemistry. Interaction of ZEB1 and EGFR signalling was interrogated in vitro by realtime quantitative PCR and immunocytochemistry in glioma stem cell lines from subtype-stratified parental tumors based on RNAseq transcriptomic gene expression profiles. RESULTS: ZEB1 is expressed in nearly all examined tumors but in highly variable fraction of cells. Co-staining of reactive CNS cell types in non-tumor tissue show that besides bona fide tumor cells, only reactive astrocytes but not leucocytes, microglia and macrophages contribute to the ZEB1 positive population. In GBM, ZEB1 labelling index correlates with EGFR amplification and IDH1 mutation, but EGFR modulation by cetuximab treatment did not alter ZEB1 expression in vitro. CONCLUSION: ZEB1 is expressed ubiquitously and across molecularily defined subtypes in glioblastoma multiforme as well as in non-GBM astroglial tumors. ZEB1 labelling varies across GBM subtypes, but a causal link remains questionable. We speculate that ZEB1 is a pan-astroglial tumor marker that reflects tumor purity.

The CORG resource (Comparative Regulatory Genomics, http://corg.eb.tuebingen.mpg.de) provides extensive cross-species comparisons of promoter regions in particular and whole gene loci in general. Pairwise as well as multiple alignments of 10 vertebrate species form the key component of CORG. We implemented a rapid alignment approach based on weight matrix motif anchors to ensure efficient computation and biologically informative alignments. All CORG workbench components have been enhanced towards more flexibility and interactivity. Reference sequence based data presentation and analysis was put into the well-known and modular Generic Genome Browser framework. Herein, various plugins facilitate online data analysis and integration with static conservation data. Main emphasis was put on the design of a new JAVA WebStart application for comparative data display. Flexible data import and export options for standard formats complete the provided services.

CORG is a versatile web-based workbench for comparative promoter analysis in vertebrate model organisms. Two kinds of information are explicitly considered in the automated annotation process. First, local conservation patterns in upstream regions of homologous genes: These phylogenetic footprints are likely to stem from sequence elements that are under selective pressure. The CORG pipeline detects and exploits patterns of local similarity to annotate promoter regions. Second, experimental data on transcription start sites: exon positions and DNA binding site descriptions complete the promoter annotation. These data are made available via an interactive web portal. Individual promoter studies are supported by a JAVA applet that supplies all data down to the nucleotide level.

CORG is a versatile web-based workbench for comparative promoter analysis in vertebrate model organisms. Two kinds of information are explicitly considered in the automated annotation process. First, local conservation patterns in upstream regions of homologous genes: These phylogenetic footprints are likely to stem from sequence elements that are under selective pressure. The CORG pipeline detects and exploits patterns of local similarity to annotate promoter regions. Second, experimental data on transcription start sites: exon positions and DNA binding site descriptions complete the promoter annotation. These data are made available via an interactive web portal. Individual promoter studies are supported by a JAVA applet that supplies all data down to the nucleotide level.

This paper presents a method for predicting biologically meaningful modules of transcription factors. For this purpose, we employ the CORG database of conserved transcription factor binding sites. We aim at enhancing the power of in-silico binding site predictions by employing three crucial constraints. First, we rely on conserved promoter regions of orthologous genes in human and mouse, second we look for synergistic transcription factor modules which bind upstream regions preferentially together, and finally we restrict our results to those modules, whose genes have a significant functional overlap. Many of our predicted binding sites coincide with known biological facts as is evidenced by a direct comparison with a single large-scale experiment for E2F binding. We also identified known combinations of transcription factors with a functional enrichment in the set of their shared target genes. Several new modules are suggested for experimental investigation. Finally we study the transcription factor network and suggest a classification of transcription factors according to their regulatory power and control.

Background Promoters are key players in gene regulation. They receive signals from various sources (e.g. cell surface receptors) and control the level of transcription initiation, which largely determines gene expression. In vertebrates, transcription start sites and surrounding regulatory elements are often poorly defined. To support promoter analysis, we present CORG http://corg.molgen.mpg.de, a framework for studying upstream regions including untranslated exons (5' UTR). Description The automated annotation of promoter regions integrates information of two kinds. First, statistically significant cross-species conservation within upstream regions of orthologous genes is detected. Pairwise as well as multiple sequence comparisons are computed. Second, binding site descriptions (position-weight matrices) are employed to predict conserved regulatory elements with a novel approach. Assembled EST sequences and verified transcription start sites are incorporated to distinguish exonic from other sequences. As of now, we have included 5 species in our analysis pipeline (man, mouse, rat, fugu and zebrafish). We characterized promoter regions of 16,127 groups of orthologous genes. All data are presented in an intuitive way via our web site. Users are free to export data for single genes or access larger data sets via our DAS server http://tomcat.molgen.mpg.de:8080/das. The benefits of our framework are exemplarily shown in the context of phylogenetic profiling of transcription factor binding sites and detection of microRNAs close to transcription start sites of our gene set. Conclusion The CORG platform is a versatile tool to support analyses of gene regulation in vertebrate promoter regions. Applications for CORG cover a broad range from studying evolution of DNA binding sites and promoter constitution to the discovery of new regulatory sequence elements (e.g. microRNAs and binding sites).

Motivation: Comparative sequence analysis is the essence of many approaches to genome annotation. Heuristic alignment algorithms utilize similar seed pairs to anchor an alignment. Some applications of local alignment algorithms (e.g. phylogenetic footprinting) would benefit from including prior knowledge (e.g. binding site motifs) in the alignment building process. Results: We introduce predefined sequence patterns as anchor points into a heuristic local alignment strategy. We extended the BLASTZ program for this purpose. A set of seed patterns is either given as consensus sequences in IUPAC code or position-weight-matrices. Phylogenetic footprinting of promoter regions is one of many potential applications for the SITEBLAST software. Availability: The source code is freely available to the academic community from http://corg.molgen.mpg.de/software Contact: christoph.dieterich@molgen.mpg.de

Sequence alignment algorithms have a long standing tradition in bioinformatics. In this paper, we formulate an extension to existing local alignment algorithms: local alignments across multiple scoring functions. For this purpose, we use the Waterman-Eggert algorithm for suboptimal local alignments as template and introduce two new features therein: 1) an alignment of two strings over a set of score functions and 2) a switch cost function 5 for penalizing jumps into a different scoring scheme within an alignment. Phylogenetic footprinting, as one potential application of this algorithm, was studied in greater detail. In this context, the right evolutionary distance and thus the scoring scheme is often not known a priori. We measured sensitivity and specificity on a test set of 21 human-rodent promoter pairs. Ultimately, we could attain a 4.5-fold enrichment of verified binding sites in our alignments.

Many cellular signaling pathways induce gene expression by activating specific transcription factor complexes. Conventional approaches to the prediction of transcription factor binding sites lead to a notoriously high number of false discoveries. To alleviate this problem, we consider only binding sites that are conserved in man-mouse genomic sequence comparisons. We employ two alternative methods for predicting binding sites: exact matches to validated binding site sequences and weight matrix scans. We then ask the question whether there is a characteristic association between a transcription factor or set thereof to a particular group of genes. Our approach is tested on genes, which are induced in dendritic cells in response to the cells’ exposure to LPS. We chose this example because the underlying signaling pathways are well understood. We demonstrate the benefit of conserved predicted binding sites in interpreting the LPS experiment. Additionally, we find that both methods for the prediction of conserved binding sites complement one another. Finally, our results suggest a distinct role for SRF in the context of LPS-induced gene expression. Contact: christoph.dieterich@molgen.mpg.de *To whom correspondence should be addressed.

Sequence conservation in non-coding, upstream regions of orthologous genes from man and mouse is likely to reflect common regulatory DNA sites. Motivated by this assumption we have delineated a catalogue of conserved non-coding sequence blocks and provide the CORG-'COmparative Regulatory Genomics'-database. The data were computed based on statistically significant local suboptimal alignments of 15 kb regions upstream of the translation start sites of, currently, 10 793 pairs of orthologous genes. The resulting conserved non-coding blocks were annotated with EST matches for easier detection of non-coding mRNA and with hits to known transcription factor binding sites. CORG data are accessible from the ENSEMBL web site via a DAS service as well as a specially developed web service (http://corg.molgen.mpg.de) for query and interactive visualization of the conserved blocks and their annotation.

Non-coding DNA segments that are conserved between the human and mouse genomic sequence are good indicators of possible regulatory sequences. Here we report on a systematic approach to delineate such conserved elements from upstream regions of orthologous gene pairs from man and mouse. We focus on orthologous genes in order to maximize our chances to find functionally similar regulatory elements. The identification of conserved elements is effected using the Waterman–Eggert local suboptimal alignment algorithm. We have modified an implementation of this algorithm such that it integrates the determination of statistical significance for the local suboptimal alignments. This has the effect of outputting a dynamically determined number of suboptimal alignments that are deemed statistically significant. Comparison with experimentally determined annotation shows a striking enrichement of regulatory sites among the conserved regions. Furthermore, the conserved regions tend to cover the promotor region described in the EPD database. Contact: martin.vingron@molgen.mpg.de

The individual and synergistic contributions of two transcription factors,EFG1andCPH1, have been characterized with regard to adhesion to, and invasion of, human epithelia byCandida albicans. For this purpose twoin vitroreconstructed tissue models were developed. A multi-layered model of human epidermis was used to simulate superficial infections of the skin, whereas a reconstructed human intestinal model was used to mimic the first steps of systemic infections. It was shown thatC. albicansdeleted for both transcription factorsCPH1andEFG1, in contrast to the congenic clinical isolate Sc5314, was neither able to adhere to, nor to penetrate, either of the model systems. A strain deleted forEFG1alone showed significant reduction in adhesion and was not able to penetrate through the stratum corneum. However, strains deleted forCPH1showed phenotypes paralleling the phenotypes of the clinical isolate Sc5314. Using different types of multi-layered human tissues reconstructedin vitrothe individual contributions of Efg1p and Cph1p to two important virulence factors ofC. albicans, namely adhesion and invasion, could be defined.

Summary: Next generation sequencing technologies facilitate genome-wide analysis of several biological processes. We are interested in whole-genome genotyping. To our knowledge, none of the existing single nucleotide polymorphism (SNP) callers consider the quality of the reference genome, which is not necessary for high-quality assemblies of well-studied model organisms. However, most genome projects will remain in draft status with little to no genome assembly improvement due to time and financial constraints. Here, we present a simple yet elegant solution (‘ACCUSA’) that considers both the read qualities as well as the reference genome's quality using a Bayesian framework. We demonstrate that ACCUSA is as good as the current SNP calling software in detecting true SNPs. More importantly, ACCUSA does not call spurious SNPs, which originate from a poor reference sequence. Availability: ACCUSA is available free of charge to academic users and may be obtained from ftp://bbc.mdc-berlin.de/software. ACCUSA is programmed in JAVA 6 and runs on any platform with JAVA support. Contact: christoph.dieterich@mdc-berlin.de Supplementary information: Supplementary data are available at Bioinformatics online.

Introduction Premature senescence is a permanent proliferative arrest that occurs in response to oncogenic signaling or DNA-damaging chemotherapy. Although tumor cell senescence has been recognized as a prognostically relevant contribution to long-term outcome post-therapy in hematological tumor models, therapeutic utilization of senescence is still hampered by an incomplete understanding of biological properties and long-term fate of senescent tumor cells in patients. Interestingly, senescence-regulating factors have recently been shown to limit reprogramming of somatic cells to pluripotency, and to protect stem cell compartments from premature exhaustion. Hence, we explore here whether cellular senescence and stemness may functionally overlap, thereby potentially equipping arrested cells with latent self-renewing potential. Methods Stem cell-related features (stem cell gene signatures, Sca-1 expression, ALDH and ABC transporter activity) were analyzed in primary apoptosis-blocked Eµ-myc transgenic B-cell lymphomas, which enter treatment-induced senescence (TIS) in response to standard antineoplastic agents. Several (e.g. Suv39h1- or p53-based) genetic models were established, in which TIS occurred in a conditional and reversible fashion. Clonogenicity, proliferative and repopulating assays were performed in vitro and in vivo, comparing individual lymphomas that grew out of senescence (“Previously Senescent”, PS) with matched lymphoma cells that equally received chemotherapy, but were incapable of entering TIS (“Never Senescent”, NS). In a mouse model of T-cell acute lymphoblastic leukemia (T-ALL), stem cell markers and tumor initiating potential were assessed in a flow-sorted non-self-renewing leukemia cell population after senescence induction by chemotherapy. Human hematological cancer cell lines and tumor samples obtained from B-cell lymphoma and acute myeloid leukemia (AML) patients were analyzed for stem cell features after exposure to senescence-inducing chemotherapy in vitro. Results Senescent mouse lymphomas were strongly skewed towards an increased expression of an adult tissue stem cell signature, distinct stem cell markers and functional stemness properties. Upon release from conditional senescence, PS cells resumed proliferation and rapidly exceeded the proliferative, clonogenic and tumor-initiating capacity of NS cells. Interrogation of self-renewal-relevant cascades revealed activation of and dependence on canonical Wnt signaling in senescence, as blocking of this pathway reduced the growth of PS, but not NS cells. Moreover, TIS-related stemness occurred independent of secretable factors. Strikingly, in a murine T-ALL model, temporary senescence enforcement re-programmed non-stem leukemia cells into leukemia stem cells, allowing PS bulk leukemia cells to de novo initiate leukemias in recipient mice. These results were supported by consistent findings in human hematological cancer cell lines as well as primary human B-cell lymphoma and acute myeloid leukemia samples. Conclusions Our findings uncover senescence-associated stemness as a detrimental capability which is latently enriched for by chemotherapy in lymphoma and leukemia. The aggressive growth potential might become evident when senescent cells occasionally acquire alterations that allow them to re-enter the cell-cycle, thereby unleashing the tumor-promoting potential of a biological program so far considered to operate as a tumor-suppressive mechanism. However, targeted intervention at stemness-related signaling cascades in senescence may open novel therapeutic options for apoptosis-resistant lymphoma and leukemia. Disclosures: No relevant conflicts of interest to declare.

The skeleton is one of the most important features for the reconstruction of vertebrate phylogeny but few data are available to understand its molecular origin. In mammals the Runt genes are central regulators of skeletogenesis. Runx2 was shown to be essential for osteoblast differentiation, tooth development, and bone formation. Both Runx2 and Runx3 are essential for chondrocyte maturation. Furthermore, Runx2 directly regulates Indian hedgehog expression, a master coordinator of skeletal development. To clarify the correlation of Runt gene evolution and the emergence of cartilage and bone in vertebrates, we cloned the Runt genes from hagfish as representative of jawless fish (MgRunxA, MgRunxB) and from dogfish as representative of jawed cartilaginous fish (ScRunx1–3). According to our phylogenetic reconstruction the stem species of chordates harboured a single Runt gene and thereafter Runt locus duplications occurred during early vertebrate evolution. All newly isolated Runt genes were expressed in cartilage according to quantitative PCR. In situ hybridisation confirmed high MgRunxA expression in hard cartilage of hagfish. In dogfish ScRunx2 and ScRunx3 were expressed in embryonal cartilage whereas all three Runt genes were detected in teeth and placoid scales. In cephalochordates (lancelets) Runt, Hedgehog and SoxE were strongly expressed in the gill bars and expression of Runt and Hedgehog was found in endo- as well as ectodermal cells. Furthermore we demonstrate that the lancelet Runt protein binds to Runt binding sites in the lancelet Hedgehog promoter and regulates its activity. Together, these results suggest that Runt and Hedgehog were part of a core gene network for cartilage formation, which was already active in the gill bars of the common ancestor of cephalochordates and vertebrates and diversified after Runt duplications had occurred during vertebrate evolution. The similarities in expression patterns of Runt genes support the view that teeth and placoid scales evolved from a homologous developmental module., Author Summary Important molecular mechanisms underlying mammalian skeletogenesis have been described but knowledge about the evolutionary origin of these gene networks is limited. The Runt gene family (Runx1–3) is of extraordinary importance for skeletogenesis. Runx2 deficient mice completely lack bone. Runx2 and Runx3 are essential for cartilage development and Runx2 regulates the key factor Indian hedgehog, which coordinates skeletogenesis. Here, we reconstructed Runt gene evolution in correlation to skeletal evolution. By analyzing lancelets, one of the closest living relatives of vertebrates, we revealed that the single Runt and Hedgehog family founder genes were co-expressed in primitive skeletal elements of the chordate stem species. Interestingly, at this stage the Runt and Hedgehog pathways were already directly linked to one another. Furthermore we isolated two Runt genes from a representative of jawless cartilaginous fish (hagfish) and three Runt genes from jawed cartilaginous fish (dogfish) which were all expressed in cartilage. The dogfish Runt genes were also found in teeth and placoid scales. This study suggests that Runt genes were involved in all ancient processes of chordate skeletogenesis. Furthermore the analysis supports the theory that most likely the gut was the tissue that originally secreted an acellular gill gut skeleton in the chordate ancestor.

Understanding transcriptional regulation of gene expression is one of the greatest challenges of modern molecular biology. A central role in this mechanism is played by transcription factors, which typically bind to specific, short DNA sequence motifs usually located in the upstream region of the regulated genes. We discuss here a simple and powerful approach for the ab initio identification of these cis-regulatory motifs. The method we present integrates several elements: human-mouse comparison, statistical analysis of genomic sequences and the concept of coregulation. We apply it to a complete scan of the human genome. By using the catalogue of conserved upstream sequences collected in the CORG database we construct sets of genes sharing the same overrepresented motif (short DNA sequence) in their upstream regions both in human and in mouse. We perform this construction for all possible motifs from 5 to 8 nucleotides in length and then filter the resulting sets looking for two types of evidence of coregulation: first, we analyze the Gene Ontology annotation of the genes in the set, searching for statistically significant common annotations; second, we analyze the expression profiles of the genes in the set as measured by microarray experiments, searching for evidence of coexpression. The sets which pass one or both filters are conjectured to contain a significant fraction of coregulated genes, and the upstream motifs characterizing the sets are thus good candidates to be the binding sites of the TF's involved in such regulation. In this way we find various known motifs and also some new candidate binding sites. We have discussed a new integrated algorithm for the "ab initio" identification of transcription factor binding sites in the human genome. The method is based on three ingredients: comparative genomics, overrepresentation, different types of coregulation. The method is applied to a full-scan of the human genome, giving satisfactory results.

Human embryonic stem cells (hESCs) exhibit high levels of proteasome activity, an intrinsic characteristic required for their self-renewal, pluripotency and differentiation. However, the mechanisms by which enhanced proteasome activity maintains hESC identity are only partially understood. Besides its essential role for the ability of hESCs to suppress misfolded protein aggregation, we hypothesize that enhanced proteasome activity could also be important to degrade endogenous regulatory factors. Since E3 ubiquitin ligases are responsible for substrate selection, we first define which E3 enzymes are increased in hESCs compared with their differentiated counterparts. Among them, we find HECT-domain E3 ligases such as HERC2 and UBE3A as well as several RING-domain E3s, including UBR7 and RNF181. Systematic characterization of their interactome suggests a link with hESC identity. Moreover, loss of distinct up-regulated E3s triggers significant changes at the transcriptome and proteome level of hESCs. However, these alterations do not dysregulate pluripotency markers and differentiation ability. On the contrary, global proteasome inhibition impairs diverse processes required for hESC identity, including protein synthesis, rRNA maturation, telomere maintenance and glycolytic metabolism. Thus, our data indicate that high proteasome activity is coupled with other determinant biological processes of hESC identity.

Dilated cardiomyopathy (DCM) is the second most common cause for heart failure with no cure except a high-risk heart transplantation. Approximately 30% of DCM patients harbor heritable mutations which are amenable to CRISPR-based gene therapy1. However, challenges related to delivery of the editing complex and off-target concerns hamper the broad applicability of CRISPR agents in the heart2. We employed a combination of the viral gene transfer vector AAVMYO with superior targeting specificity of heart muscle tissue3and CRISPR base editors to repair patient mutations in the cardiac splice factorRbm20, which cause aggressive and arrhythmogenic DCM4. Using optimized conditions, we could improve splice defects in human iPSC-derived cardiomyocytes (iPSC-CMs) and repair >70% of cardiomyocytes in twoRbm20knock-in mouse models that we generated to serve as anin vivoplatform of our editing strategy. Treatment of juvenile mice restored the localization defect of RBM20 in 75% of cells and splicing of RBM20 targets including TTN. Three months after injection, cardiac dilation and ejection fraction reached wildtype levels. Single-nuclei RNA sequencing (snRNA-seq) uncovered restoration of the transcriptional profile across all major cardiac cell types and whole-genome sequencing (WGS) revealed no evidence for aberrant off-target editing. Our study highlights the potential of base editors combined with AAVMYO to achieve gene repair for treatment of hereditary cardiac diseases.