Ewart de Bruijn

In situ detection of miRNAs in animal embryos using LNA-modified oligonucleotide probes

MicroRNAs (miRNAs) are 20–23 nucleotide (nt) RNA molecules that regulate gene expression post-transcriptionally. A key step toward understanding the function of the hundreds of miRNAs identified in animals is to determine their expression during development. Here we performed a detailed analysis of conditions for in situ detection of miRNAs in the zebrafish embryo using locked nucleic acid (LNA)-modified DNA probes and report expression patterns for 15 miRNAs in the mouse embryo.

MicroRNAs show a wide diversity of expression profiles in the developing and mature central nervous system

BackgroundMicroRNA (miRNA) encoding genes are abundant in vertebrate genomes but very few have been studied in any detail. Bioinformatic tools allow prediction of miRNA targets and this information coupled with knowledge of miRNA expression profiles facilitates formulation of hypotheses of miRNA function. Although the central nervous system (CNS) is a prominent site of miRNA expression, virtually nothing is known about the spatial and temporal expression profiles of miRNAs in the brain. To provide an overview of the breadth of miRNA expression in the CNS, we performed a comprehensive analysis of the neuroanatomical expression profiles of 38 abundant conserved miRNAs in developing and adult zebrafish brain.ResultsOur results show miRNAs have a wide variety of different expression profiles in neural cells, including: expression in neuronal precursors and stem cells (for example, miR-92b); expression associated with transition from proliferation to differentiation (for example, miR-124); constitutive expression in mature neurons (miR-124 again); expression in both proliferative cells and their differentiated progeny (for example, miR-9); regionally restricted expression (for example, miR-222 in telencephalon); and cell-type specific expression (for example, miR-218a in motor neurons).ConclusionThe data we present facilitate prediction of likely modes of miRNA function in the CNS and many miRNA expression profiles are consistent with the mutual exclusion mode of function in which there is spatial or temporal exclusion of miRNAs and their targets. However, some miRNAs, such as those with cell-type specific expression, are more likely to be co-expressed with their targets. Our data provide an important resource for future functional studies of miRNAs in the CNS.

A systematic genome-wide analysis of zebrafish protein-coding gene function

Since the publication of the human reference genome, the identities of specific genes associated with human diseases are being discovered at a rapid rate. A central problem is that the biological activity of these genes is often unclear. Detailed investigations in model vertebrate organisms, typically mice, have been essential for understanding the activities of many orthologues of these disease-associated genes. Although gene-targeting approaches and phenotype analysis have led to a detailed understanding of nearly 6,000 protein-coding genes, this number falls considerably short of the more than 22,000 mouse protein-coding genes. Similarly, in zebrafish genetics, one-by-one gene studies using positional cloning, insertional mutagenesis, antisense morpholino oligonucleotides, targeted re-sequencing, and zinc finger and TAL endonucleases have made substantial contributions to our understanding of the biological activity of vertebrate genes, but again the number of genes studied falls well short of the more than 26,000 zebrafish protein-coding genes. Importantly, for both mice and zebrafish, none of these strategies are particularly suited to the rapid generation of knockouts in thousands of genes and the assessment of their biological activity. Here we describe an active project that aims to identify and phenotype the disruptive mutations in every zebrafish protein-coding gene, using a well-annotated zebrafish reference genome sequence, high-throughput sequencing and efficient chemical mutagenesis. So far we have identified potentially disruptive mutations in more than 38% of all known zebrafish protein-coding genes. We have developed a multi-allelic phenotyping scheme to efficiently assess the effects of each allele during embryogenesis and have analysed the phenotypic consequences of over 1,000 alleles. All mutant alleles and data are available to the community and our phenotyping scheme is adaptable to phenotypic analysis beyond embryogenesis.

Limitations and possibilities of small RNA digital gene expression profiling

1. Li, H., Ruan, J. & Durbin, R. Genome Res. 18, 1851–1858 (2008) 2. Hormozdiari, F. et al. Genome Res. (in the press). 3. Korbel, J.O. et al. Genome Biol. 10, R23 (2009). 4. Lee, S., Cheran, E. & Brudno, M. Bioinformatics 24, i59–i67 (2008). 5. Mills, R.E. et al. Genome Res. 16, 1182–1190 (2006). 6. Bentley, D.R. et al. Nature 456, 53–59 (2008). 7. Kidd, J.M. et al. Nature 453, 56–64 (2008). Gaussian distribution with mean μ and standard deviation σ where

Chromothripsis as a mechanism driving complex de novo structural rearrangements in the germline

A variety of mutational mechanisms shape the dynamic architecture of human genomes and occasionally result in congenital defects and disease. Here, we used genome-wide long mate-pair sequencing to systematically screen for inherited and de novo structural variation in a trio including a child with severe congenital abnormalities. We identified 4321 inherited structural variants and 17 de novo rearrangements. We characterized the de novo structural changes to the base-pair level revealing a complex series of balanced inter- and intra-chromosomal rearrangements consisting of 12 breakpoints involving chromosomes 1, 4 and 10. Detailed inspection of breakpoint regions indicated that a series of simultaneous double-stranded DNA breaks caused local shattering of chromosomes. Fusion of the resulting chromosomal fragments involved non-homologous end joining, since junction points displayed limited or no homology and small insertions and deletions. The pattern of random joining of chromosomal fragments that we observe here strongly resembles the somatic rearrangement patterns--termed chromothripsis--that have recently been described in deranged cancer cells. We conclude that a similar mechanism may also drive the formation of de novo structural variation in the germline.

Extensive localization of long noncoding RNAs to the cytosol and mono- and polyribosomal complexes.

BackgroundLong noncoding RNAs (lncRNAs) form an abundant class of transcripts, but the function of the majority of them remains elusive. While it has been shown that some lncRNAs are bound by ribosomes, it has also been convincingly demonstrated that these transcripts do not code for proteins. To obtain a comprehensive understanding of the extent to which lncRNAs bind ribosomes, we performed systematic RNA sequencing on ribosome-associated RNA pools obtained through ribosomal fractionation and compared the RNA content with nuclear and (non-ribosome bound) cytosolic RNA pools.ResultsThe RNA composition of the subcellular fractions differs significantly from each other, but lncRNAs are found in all locations. A subset of specific lncRNAs is enriched in the nucleus but surprisingly the majority is enriched in the cytosol and in ribosomal fractions. The ribosomal enriched lncRNAs include H19 and TUG1.ConclusionsMost studies on lncRNAs have focused on the regulatory function of these transcripts in the nucleus. We demonstrate that only a minority of all lncRNAs are nuclear enriched. Our findings suggest that many lncRNAs may have a function in cytoplasmic processes, and in particular in ribosome complexes.

Primary Colorectal Cancers and Their Subsequent Hepatic Metastases Are Genetically Different: Implications for Selection of Patients for Targeted Treatment

Purpose: In the era of DNA-guided personalized cancer treatment, it is essential to conduct predictive analysis on the tissue that matters. Here, we analyzed genetic differences between primary colorectal adenocarcinomas (CRC) and their respective hepatic metastasis. Experimental Design: The primary CRC and the subsequent hepatic metastasis of 21 patients with CRC were analyzed using targeted deep-sequencing of DNA isolated from formalin-fixed, paraffin-embedded archived material. Results: We have interrogated the genetic constitution of a designed “Cancer Mini-Genome” consisting of all exons of 1,264 genes associated with pathways relevant to cancer. In total, 6,696 known and 1,305 novel variations were identified in 1,174 and 667 genes, respectively, including 817 variants that potentially altered protein function. On average, 83 (SD = 69) potentially function-impairing variations were gained in the metastasis and 70 (SD = 48) variations were lost, showing that the primary tumor and hepatic metastasis are genetically significantly different. Besides novel and known variations in genes such as KRAS, BRAF, KDR, FLT1, PTEN, and PI3KCA, aberrations in the up/downstream genes of EGFR/PI3K/VEGF-pathways and other pathways (mTOR, TGFβ, etc.) were also detected, potentially influencing therapeutic responsiveness. Chemotherapy between removal of the primary tumor and the metastasis (N = 11) did not further increase the amount of genetic variation. Conclusion: Our study indicates that the genetic characteristics of the hepatic metastases are different from those of the primary CRC tumor. As a consequence, the choice of treatment in studies investigating targeted therapies should ideally be based on the genetic properties of the metastasis rather than on those of the primary tumor. Clin Cancer Res; 18(3); 688–99. ©2011 AACR.

Accurate SNP and mutation detection by targeted custom microarray-based genomic enrichment of short-fragment sequencing libraries

Microarray-based enrichment of selected genomic loci is a powerful method for genome complexity reduction for next-generation sequencing. Since the vast majority of exons in vertebrate genomes are smaller than 150 nt, we explored the use of short fragment libraries (85–110 bp) to achieve higher enrichment specificity by reducing carryover and adverse effects of flanking intronic sequences. High enrichment specificity (60–75%) was obtained with a relative even base coverage. Up to 98% of the target-sequence was covered more than 20× at an average coverage depth of about 200×. To verify the accuracy of SNP/mutation detection, we evaluated 384 known non-reference SNPs in the targeted regions. At ∼200× average sequence coverage, we were able to survey 96.4% of 1.69 Mb of genomic sequence with only 4.2% false negative calls, mostly due to low coverage. Using the same settings, a total of 1197 novel candidate variants were detected. Verification experiments revealed only eight false positive calls, indicating an overall false positive rate of less than 1 per ∼200 000 bp. Taken together, short fragment libraries provide highly efficient and flexible enrichment of exonic targets and yield relatively even base coverage, which facilitates accurate SNP and mutation detection. Raw sequencing data, alignment files and called SNPs have been submitted into GEO database http://www.ncbi.nlm.nih.gov/geo/ with accession number GSE18542.

Small RNA expression and strain specificity in the rat

BackgroundDigital gene expression (DGE) profiling has become an established tool to study RNA expression. Here, we provide an in-depth analysis of small RNA DGE profiles from two different rat strains (BN-Lx and SHR) from six different rat tissues (spleen, liver, brain, testis, heart, kidney). We describe the expression patterns of known and novel micro (mi)RNAs and piwi-interacting (pi)RNAs.ResultsWe confirmed the expression of 588 known miRNAs (54 in antisense orientation) and identified 56 miRNAs homologous to known human or mouse miRNAs, as well as 45 new rat miRNAs. Furthermore, we confirmed specific A to I editing in brain for mir-376a/b/c and identified mir-377 as a novel editing target. In accordance with earlier findings, we observed a highly tissue-specific expression pattern for all tissues analyzed. The brain was found to express the highest number of tissue-specific miRNAs, followed by testis. Notably, our experiments also revealed robust strain-specific differential miRNA expression in the liver that is caused by genetic variation between the strains. Finally, we identified two types of germline-specific piRNAs in testis, mapping either to transposons or in strand-specific clusters.ConclusionsTaken together, the small RNA compendium described here advances the annotation of small RNAs in the rat genome. Strain and tissue-specific expression patterns furthermore provide a strong basis for studying the role of small RNAs in regulatory networks as well as biological process like physiology and neurobiology that are extensively studied in this model system.

Efficient double fragmentation ChIP-seq provides nucleotide resolution protein-DNA binding profiles

Immunoprecipitated crosslinked protein-DNA fragments typically range in size from several hundred to several thousand base pairs, with a significant part of chromatin being much longer than the optimal length for next-generation sequencing (NGS) procedures. Because these larger fragments may be non-random and represent relevant biology that may otherwise be missed, but also because they represent a significant fraction of the immunoprecipitated material, we designed a double-fragmentation ChIP-seq procedure. After conventional crosslinking and immunoprecipitation, chromatin is de-crosslinked and sheared a second time to concentrate fragments in the optimal size range for NGS. Besides the benefits of increased chromatin yields, the procedure also eliminates a laborious size-selection step. We show that the double-fragmentation ChIP-seq approach allows for the generation of biologically relevant genome-wide protein-DNA binding profiles from sub-nanogram amounts of TCF7L2/TCF4, TBP and H3K4me3 immunoprecipitated material. Although optimized for the AB/SOLiD platform, the same approach may be applied to other platforms.