Megabase deletions of gene deserts result in viable mice

FOXO1-FGFR1 fusion and amplification in a solid variant of alveolar rhabdomyosarcoma

Rhabdomyosarcoma is the most common pediatric soft tissue malignancy. Two major subtypes, alveolar rhabdomyosarcoma and embryonal rhabdomyosarcoma, constitute 20 and 60% of all cases, respectively. Approximately 80% of alveolar rhabdomyosarcoma carry two signature chromosomal translocations, t(2;13)(q35;q14) resulting in PAX3-FOXO1 fusion, and t(1;13)(p36;q14) resulting in PAX7-FOXO1 fusion. Whether the remaining cases are truly negative for gene fusion has been questioned. We are reporting the case of a 9-month-old girl with a metastatic neck mass diagnosed histologically as solid variant alveolar rhabdomyosarcoma. Chromosome analysis showed a t(8;13;9)(p11.2;q14;9q32) three-way translocation as the sole clonal aberration. Fluorescent in situ hybridization (FISH) demonstrated a rearrangement at the FOXO1 locus and an amplification of its centromeric region. Single-nucleotide polymorphism-based microarray analysis illustrated a co-amplification of the FOXO1 gene at 13q14 and the FGFR1 gene at 8p12p11.2, suggesting formation and amplification of a chimerical FOXO1-FGFR1 gene. This is the first report to identify a novel fusion partner FGFR1 for the known anchor gene FOXO1 in alveolar rhabdomyosarcoma.

What fraction of the human genome is functional?

Many evolutionary studies over the past decade have estimated α(sel), the proportion of all nucleotides in the human genome that are subject to purifying selection because of their biological function. Most of these studies have estimated the nucleotide substitution rates from genome sequence alignments across many diverse mammals. Some α(sel) estimates will be affected by the heterogeneity of substitution rates in neutral sequence across the genome. Most will also be inaccurate if change in the functional sequence repertoire occurs rapidly relative to the separation of lineages that are being compared. Evidence gathered from both evolutionary and experimental analyses now indicate that rates of "turnover" of functional, predominantly noncoding, sequence are, indeed, high. They are sufficiently high that an estimated 50% of mouse constrained noncoding sequence is predicted not to be shared with rat, a closely related rodent. The rapidity of turnover results in, at least, a twofold underestimate of α(sel) by analyses that measure constraint across the eutherian phylogeny. Approaches that take account of turnover estimate that the steady-state value of α(sel) lies between 10% and 15%. Experimental studies corroborate the predicted rates of loss and gain of noncoding functional sites. These studies show the limitations inherent in the use of deep sequence conservation for identifying functional sequence. Experimental investigations focusing on lineage-specific, noncoding, and functional sequence are now essential if we are to appreciate the complete functional repertoire of the human genome.

[Cis-ruptions of highly conserved non-coding genomic elements distant from the SOX9 gene in the Pierre Robin sequence]

Major developmental genes, exhibiting complex expression patterns, are often embedded within a genic desert particularly rich in regions, which though non-coding are highly conserved. The developmental expression of these genes in many areas requires coordinated regulation in time and space, which is orchestrated by some of these conserved non-coding regions, acting as transcriptional regulators. SOX9 is an essential gene for many developmental processes, such as chondrogenesis, migration and differentiation of neural crest cells and testis development. In agreement with these major expression areas, SOX9 haploinsufficiency, linked to alterations in coding sequence, leads to a polymorphic malformation syndrome - campomelic dysplasia - whose major symptoms are a bone anomaly, a Pierre Robin sequence, and a sexual differentiation anomaly (Disorder of Sex Development, DSD). SOX9 is located in a  ~2.5 Mb gene desert extremely rich in conserved sequences. We have used the SOX9 locus and campomelic dysplasia as a model to show that one or several endophenotypes within a complex syndrome may arise from a tissue-specific deregulation of a major developmental gene transcription. Our work has focused on one of these endophenotypes, SPR, characterized by the triad micro- and/or retrognathy, glossoptosis and cleft palate. Here we report in detail how we identified alterations (translocations, deletions, point mutations) in non-coding regions, located far away (more than 1.2 Mb) upstream and downstream of SOX9, in clustered or sporadic SPR cases.

Molecular genetic analyses of polytene chromosome region 72A-D in Drosophila melanogaster reveal a gene desert in 72D

We have investigated a region of ∼310 kb of genomic DNA within polytene chromosome subdivisions 72A to 72D of Drosophila melanogaster. This region includes 57 predicted protein-coding genes. Seventeen of these genes are in six clusters that appear to have arisen by tandem duplication. Within this region we found 23 complementation groups that are essential for zygotic viability, and we have identified the transcription units for 18 of the 23. We also found a 55 kb region in 72D that is nonessential. Flies deficient for this region are viable and fertile. Within this nonessential region are 48 DNA sequences of 12 to 33 base pairs that are completely conserved among 12 distantly related Drosophila species. These sequences do not have the evolutionary signature of conserved protein-coding DNA sequences, nor do they appear to encode microRNAs, however, the strong selection suggests functions in wild populations that are not apparent in laboratory cultures. This region resembles dispensable gene deserts previously characterized in the mouse genome.

Controlled somatic and germline copy number variation in the mouse model

Changes in the number of chromosomes, but also variations in the copy number of chromosomal regions have been described in various pathological conditions, such as cancer and aneuploidy, but also in normal physiological condition. Our classical view of DNA replication and mitotic preservation of the chromosomal integrity is now challenged as new technologies allow us to observe such mosaic somatic changes in copy number affecting regions of chromosomes with various sizes. In order to go further in the understanding of copy number influence in normal condition we could take advantage of the novel strategy called Targeted Asymmetric Sister Chromatin Event of Recombination (TASCER) to induce recombination during the G2 phase so that we can generate deletions and duplications of regions of interest prior to mitosis. Using this approach in the mouse we could address the effects of copy number variation and segmental aneuploidy in daughter cells and allow us to explore somatic mosaics for large region of interest in the mouse.

Functional characterization of liver enhancers that regulate drug-associated transporters

Little is known about how genetic variations in enhancers influence drug response. In this study, we investigated whether nucleotide variations in enhancers that regulate drug transporters can alter their expression levels. Using comparative genomics and liver-specific transcription factor binding site (TFBS) analyses, we identified evolutionary conserved regions (ECRs) surrounding nine liver membrane transporters that interact with commonly used pharmaceuticals. The top 50 ECRs were screened for enhancer activity in vivo, of which five--located around ABCB11, SLC10A1, SLCO1B1, SLCO1A2, and SLC47A1--exhibited significant enhancer activity. Common variants identified in a large ethnically diverse cohort (n = 272) were assayed for differential enhancer activity, and three variants were found to have significant effects on reporter activity as compared with the reference allele. In addition, one variant was associated with reduced SLCO1A2 mRNA expression levels in human liver tissues, and another was associated with increased methotrexate (MTX) clearance in patients. This work provides a general model for the rapid characterization of liver enhancers and identifies associations between enhancer variants and drug response.

Potential novel mechanism for Axenfeld-Rieger syndrome: deletion of a distant region containing regulatory elements of PITX2

PURPOSE

Mutations in PITX2 are associated with Axenfeld-Rieger syndrome (ARS), which involves ocular, dental, and umbilical abnormalities. Identification of cis-regulatory elements of PITX2 is important to better understand the mechanisms of disease.

METHODS

Conserved noncoding elements surrounding PITX2/pitx2 were identified and examined through transgenic analysis in zebrafish; expression pattern was studied by in situ hybridization. Patient samples were screened for deletion/duplication of the PITX2 upstream region using arrays and probes.

RESULTS

Zebrafish pitx2 demonstrates conserved expression during ocular and craniofacial development. Thirteen conserved noncoding sequences positioned within a gene desert as far as 1.1 Mb upstream of the human PITX2 gene were identified; 11 have enhancer activities consistent with pitx2 expression. Ten elements mediated expression in the developing brain, four regions were active during eye formation, and two sequences were associated with craniofacial expression. One region, CE4, located approximately 111 kb upstream of PITX2, directed a complex pattern including expression in the developing eye and craniofacial region, the classic sites affected in ARS. Screening of ARS patients identified an approximately 7600-kb deletion that began 106 to 108 kb upstream of the PITX2 gene, leaving PITX2 intact while removing regulatory elements CE4 to CE13.

CONCLUSIONS

These data suggest the presence of a complex distant regulatory matrix within the gene desert located upstream of PITX2 with an essential role in its activity and provides a possible mechanism for the previous reports of ARS in patients with balanced translocations involving the 4q25 region upstream of PITX2 and the current patient with an upstream deletion.

When needles look like hay: how to find tissue-specific enhancers in model organism genomes

A major prerequisite for the investigation of tissue-specific processes is the identification of cis-regulatory elements. No generally applicable technique is available to distinguish them from any other type of genomic non-coding sequence. Therefore, researchers often have to identify these elements by elaborate in vivo screens, testing individual regions until the right one is found. Here, based on many examples from the literature, we summarize how functional enhancers have been isolated from other elements in the genome and how they have been characterized in transgenic animals. Covering computational and experimental studies, we provide an overview of the global properties of cis-regulatory elements, like their specific interactions with promoters and target gene distances. We describe conserved non-coding elements (CNEs) and their internal structure, nucleotide composition, binding site clustering and overlap, with a special focus on developmental enhancers. Conflicting data and unresolved questions on the nature of these elements are highlighted. Our comprehensive overview of the experimental shortcuts that have been found in the different model organism communities and the new field of high-throughput assays should help during the preparation phase of a screen for enhancers. The review is accompanied by a list of general guidelines for such a project.

Enrichment of mRNA-like noncoding RNAs in the divergence of Drosophila males

With the advent of transcriptome data, it has become clear that mRNA-like noncoding RNAs (mlncRNAs) are widespread in eukaryotes. Although their functions are poorly understood, these transcripts may play an important role in development and could thus be involved in determining developmental complexity and phenotypic diversification. However, few studies have assessed their potential roles in the divergence of closely related species. Here, we identify and study patterns of sequence and expression divergence in ten novel candidate mlncRNAs from Drosophila pseudoobscura and its close relative D. persimilis. The candidate mlncRNAs were identified by randomly sequencing a group of 734 cDNA clones from a microarray that showed either no difference in expression (187 clones) or differential expression (547 clones) in comparisons between D. pseudoobscura and D. persimilis and between these two species and their F(1) hybrids. Candidate mlncRNAs are overrepresented among differentially expressed transcripts between males of D. pseudoobscura and D. persimilis, and although they have high sequence conservation between these two species, seven of them have no putative homologs in any of the other ten Drosophila species whose genomes have been sequenced. Expression of eight of the ten candidate mlncRNAs was detected either in whole bodies (adults) or testes using a custom-designed oligonucleotide microarray. Three of the ten candidate mlncRNAs are highly expressed (in the top 4% of the male transcriptome), differentially expressed between species, and show extreme levels of sex-bias, with one transcript having the highest level of male bias in the whole transcriptome. Proteomic data from testes show no traces of any predicted peptides from the candidate mlncRNAs. Our results suggest that these mlncRNAs may be important in male-specific processes related to sexual dimorphism and species divergence in this species group.

Genomics of long-range regulatory elements

Transcriptional regulation of gene expression plays a significant role in establishing the diversity of human cell types and biological functions from a common set of genes. The components of regulatory control in the human genome include cis-acting elements that act across immense genomic distances to influence the spatial and temporal distribution of gene expression. Here we review the established categories of distant-acting regulatory elements, discussing the classical and contemporary evidence of their regulatory potential and clinical importance. Current efforts to identify regulatory sequences throughout the genome and elucidate their biological significance depend heavily on advances in sequence conservation-based analyses and on increasingly large-scale efforts applying transgenic technologies in model organisms. We discuss the advantages and limitations of sequence conservation as a predictor of regulatory function and present complementary emerging technologies now being applied to annotate regulatory elements in vertebrate genomes.

Alu and b1 repeats have been selectively retained in the upstream and intronic regions of genes of specific functional classes

Alu and B1 repeats are mobile elements that originated in an initial duplication of the 7SL RNA gene prior to the primate-rodent split about 80 million years ago and currently account for a substantial fraction of the human and mouse genome, respectively. Following the primate-rodent split, Alu and B1 elements spread independently in each of the two genomes in a seemingly random manner, and, according to the prevailing hypothesis, negative selection shaped their final distribution in each genome by forcing the selective loss of certain Alu and B1 copies. In this paper, contrary to the prevailing hypothesis, we present evidence that Alu and B1 elements have been selectively retained in the upstream and intronic regions of genes belonging to specific functional classes. At the same time, we found no evidence for selective loss of these elements in any functional class. A subset of the functional links we discovered corresponds to functions where Alu involvement has actually been experimentally validated, whereas the majority of the functional links we report are novel. Finally, the unexpected finding that Alu and B1 elements show similar biases in their distribution across functional classes, despite having spread independently in their respective genomes, further supports our claim that the extant instances of Alu and B1 elements are the result of positive selection.

Despite their fundamental role in cell regulation, genes account for less than 1% of the human genome. Recent studies have shown that non-genic regions of our DNA may also play an important functional role in human cells. In this paper, we study Alu and B elements, a specific class of such non-genic elements that account for ∼10% of the human genome and ∼7% of the mouse genome respectively. We show that, contrary to the prevailing hypothesis, Alu and B elements have been preferentially retained in the proximity of genes that perform specific functions in the cell. In contrast, we found no evidence for selective loss of these elements in any functional class. Several of the functional classes that we have linked to Alu and B elements are central to the proper working of the cell, and their disruption has previously been shown to lead to the onset of disease. Interestingly, the DNA sequences of Alu and B elements differ substantially between human and mouse, thus hinting at the existence of a potentially large number of non-conserved regulatory elements.

Establishing legitimacy and function in the new transcriptome

The last decade has seen an explosion of interest in new classes of non-coding RNA. While some are now firmly established as new categories of legitimate functional RNAs, the purpose and even existence of others remain to be solidified. Here, we discuss the challenges associated with discovery and characterization of non-traditional categories of non-coding RNA.

Genome desertification in eutherians: can gene deserts explain the uneven distribution of genes in placental mammalian genomes?

The evolution of genome size as well as structure and organization of genomes belongs among the key questions of genome biology. Here we show, based on a comparative analysis of 30 genomes, that there is generally a tight correlation between the number of genes per chromosome and the length of the respective chromosome in eukaryotic genomes. The surprising exceptions to this pattern are placental mammalian genomes. We identify the number and, more importantly, the uneven distribution of gene deserts among chromosomes, i.e., long (>500 kb) stretches of DNA that do not encode for genes, as the main contributing factor for the observed anomaly of eutherian genomes. Gene-rich placental mammalian chromosomes have smaller proportions of gene deserts and vice versa. We show that the uneven distribution of gene deserts is a derived character state of eutherians. The functional and evolutionary significance of this particular feature of eutherian genomes remains to be explained.

Comparative genomics and evolution of conserved noncoding elements (CNE) in rainbow trout

Background

Recent advances in the accumulation of genetic mapping and DNA sequence information from several salmonid species support the long standing view of an autopolyploid origin of these fishes (i.e., 4R). However, the paralogy relationships of the chromosomal segments descendent from earlier polyploidization events (i.e., 2R/3R) largely remain unknown, mainly due to an unbalanced pseudogenization of paralogous genes that were once resident on the ancient duplicated segments. Inter-specific conserved noncoding elements (CNE) might hold the key in identifying these regions, if they are associated with arrays of genes that have been highly conserved in syntenic blocks through evolution. To test this hypothesis, we investigated the chromosomal positions of subset of CNE in the rainbow trout genome using a comparative genomic framework.

Results

Through a genome wide analysis, we selected 41 pairs of adjacent CNE located on various chromosomes in zebrafish and obtained their intervening, less conserved, sequence information from rainbow trout. We identified 56 distinct fragments corresponding to about 150 Kbp of sequence data that were localized to 67 different chromosomal regions in the rainbow trout genome. The genomic positions of many duplicated CNE provided additional support for some previously suggested homeologies in this species. Additionally, we now propose 40 new potential paralogous affinities by analyzing the variation in the segregation patterns of some multi-copy CNE along with the synteny association comparison using several model vertebrates. Some of these regions appear to carry signatures of the 1R, 2R or 3R duplications. A subset of these CNE markers also demonstrated high utility in identifying homologous chromosomal segments in the genomes of Atlantic salmon and Arctic charr.

Conclusion

CNE seem to be more efficacious than coding sequences in providing insights into the ancient paralogous affinities within the vertebrate genomes. Such a feature makes these elements extremely attractive for comparative genomics studies, as they can be treated as 'anchor' markers to investigate the association of distally located candidate genes on the homologous genomic segments of closely or distantly related organisms.

Identifying regulatory elements in eukaryotic genomes

Proper development and functioning of an organism depends on precise spatial and temporal expression of all its genes. These coordinated expression-patterns are maintained primarily through the process of transcriptional regulation. Transcriptional regulation is mediated by proteins binding to regulatory elements on the DNA in a combinatorial manner, where particular combinations of transcription factor binding sites establish specific regulatory codes. In this review, we survey experimental and computational approaches geared towards the identification of proximal and distal gene regulatory elements in the genomes of complex eukaryotes. Available approaches that decipher the genetic structure and function of regulatory elements by exploiting various sources of information like gene expression data, chromatin structure, DNA-binding specificities of transcription factors, cooperativity of transcription factors, etc. are highlighted. We also discuss the relevance of regulatory elements in the context of human health through examples of mutations in some of these regions having serious implications in misregulation of genes and being strongly associated with human disorders.

The problem of the eukaryotic genome size

The current state of knowledge concerning the unsolved problem of the huge interspecific eukaryotic genome size variations not correlating with the species phenotypic complexity (C-value enigma also known as C-value paradox) is reviewed. Characteristic features of eukaryotic genome structure and molecular mechanisms that are the basis of genome size changes are examined in connection with the C-value enigma. It is emphasized that endogenous mutagens, including reactive oxygen species, create a constant nuclear environment where any genome evolves. An original quantitative model and general conception are proposed to explain the C-value enigma. In accordance with the theory, the noncoding sequences of the eukaryotic genome provide genes with global and differential protection against chemical mutagens and (in addition to the anti-mutagenesis and DNA repair systems) form a new, third system that protects eukaryotic genetic information. The joint action of these systems controls the spontaneous mutation rate in coding sequences of the eukaryotic genome. It is hypothesized that the genome size is inversely proportional to functional efficiency of the anti-mutagenesis and/or DNA repair systems in a particular biological species. In this connection, a model of eukaryotic genome evolution is proposed.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells

The fraction of experimentally active conserved non-coding sequences within any given cell type is low, so classical assays are unlikely to expose their potential.

Background

Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian.

Results

We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types.

Conclusions

The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Identification and characterization of new long conserved noncoding sequences in vertebrates

Comparative sequence analyses have identified highly conserved genomic DNA sequences, including noncoding sequences, between humans and other species. By performing whole-genome comparisons of human and mouse, we have identified 611 conserved noncoding sequences longer than 500 bp, with more than 95% identity between the species. These long conserved noncoding sequences (LCNS) include 473 new sequences that do not overlap with previously reported ultraconserved elements (UCE), which are defined as aligned sequences longer than 200 bp with 100% identity in human, mouse, and rat. The LCNS were distributed throughout the genome except for the Y chromosome and often occurred in clusters within regions with a low density of coding genes. Many of the LCNS were also highly conserved in other mammals, chickens, frogs, and fish; however, we were unable to find orthologous sequences in the genomes of invertebrate species. In order to examine whether these conserved sequences are functionally important or merely mutational cold spots, we directly measured the frequencies of ENU-induced germline mutations in the LCNS of the mouse. By screening about 40.7 Mb, we found 35 mutations, including mutations at nucleotides that were conserved between human and fish. The mutation frequencies were equivalent to those found in other genomic regions, including coding sequences and introns, suggesting that the LCNS are not mutational cold spots at all. Taken together, these results suggest that mutations occur with equal frequency in LCNS but are eliminated by natural selection during the course of evolution.

Noncoding RNA in development

Non-protein-coding sequences increasingly dominate the genomes of multicellular organisms as their complexity increases, in contrast to protein-coding genes, which remain relatively static. Most of the mammalian genome and indeed that of all eukaryotes is expressed in a cell- and tissue-specific manner, and there is mounting evidence that much of this transcription is involved in the regulation of differentiation and development. Different classes of small and large noncoding RNAs (ncRNAs) have been shown to regulate almost every level of gene expression, including the activation and repression of homeotic genes and the targeting of chromatin-remodeling complexes. ncRNAs are involved in developmental processes in both simple and complex eukaryotes, and we illustrate this in the latter by focusing on the animal germline, brain, and eye. While most have yet to be systematically studied, the emerging evidence suggests that there is a vast hidden layer of regulatory ncRNAs that constitutes the majority of the genomic programming of multicellular organisms and plays a major role in controlling the epigenetic trajectories that underlie their ontogeny.

Dispensability of mammalian DNA

In the lab, the cis-regulatory network seems to exhibit great functional redundancy. Many experiments testing enhancer activity of neighboring cis-regulatory elements show largely overlapping expression domains. Of recent interest, mice in which cis-regulatory ultraconserved elements were knocked out showed no obvious phenotype, further suggesting functional redundancy. Here, we present a global evolutionary analysis of mammalian conserved nonexonic elements (CNEs), and find strong evidence to the contrary. Given a set of CNEs conserved between several mammals, we characterize functional dispensability as the propensity for the ancestral element to be lost in mammalian species internal to the spanned species tree. We show that ultraconserved-like elements are over 300-fold less likely than neutral DNA to have been lost during rodent evolution. In fact, many thousands of noncoding loci under purifying selection display near uniform indispensability during mammalian evolution, largely irrespective of nucleotide conservation level. These findings suggest that many genomic noncoding elements possess functions that contribute noticeably to organism fitness in naturally evolving populations.

Characterizing the mouse ES cell transcriptome with Illumina sequencing

Large datasets generated by Illumina sequencing are ideally suited to transcriptome characterization. We generated 3,052,501 27-mer reads from F1 mouse embryonic stem (ES) cell cDNA. Using the ELAND alignment tool, 74.5% of reads matched sequenced mouse resources, <1% were contaminants, and 3.7% failed quality control. Of the reads, 21.6% did not match mouse sequences using ELAND, but most of them were successfully aligned with mouse mRNAs using MegaBLAST. We conclude that most of the reads in the dataset are derived from mouse transcripts. A total of 14,434 mouse RefSeq genes were represented by at least 1 read. A Pearson correlation coefficient of 0.7 between Illumina sequencing and Illumina array expression data suggested similar results for both technologies. A weak 3' bias of reads was found. Reads from genes with low expression had lower GC content than the corresponding RefSeq genes, indicating a GC bias. Biases were confirmed with further Illumina read datasets generated with cDNA from mouse brain and from mutagen-treated F1 ES cells. We calculated relative expression values, because transcript length and read number were correlated. In the absence of signal saturation or background noise, we believe that short-read sequencing technologies will have a major impact on gene expression studies in the near future.

How can we realize the promise of personalized antidepressant medicines?

Personalized medication that is based on pharmacogenetic data has long been expected to improve the efficacy of treatments for neurological and psychiatric disorders, including depression. However, the complexity of the regulation of gene transcription and its interactions with environmental factors means that straightforward translation of individual genetic information into tailored treatment is unlikely. Nevertheless, when data from genomics, proteomics, metabolomics, neuroimaging and neuroendocrinology are used in combination, they could lead to the development of effective personalized antidepressant treatment that is based on both genotypes and biomarkers. This process will require many further steps and collaboration between basic and clinical neuroscience.

Analysis of the effects of exposure to acute hypoxia on oxidative lesions and tumour progression in a transgenic mouse breast cancer model

Background

Tumour hypoxia is known to be a poor prognostic indicator, predictive of increased risk of metastatic disease and reduced survival. Genomic instability has been proposed as one of the potential mechanisms for hypoxic tumour progression. Both of these features are commonly found in many cancer types, but their relationship and association with tumour progression has not been examined in the same model.

Methods

To address this issue, we determined the effects of 6 week in vivo acute hypoxic exposure on the levels of mutagenic lipid peroxidation product, malondialdehyde, and 8-oxo-7,8-dihydro-2'-deoxyguanosine DNA (8-oxo-dG) lesions in the transgenic polyomavirus middle T (PyMT) breast cancer mouse model.

Results

We observed significantly increased plasma lipid peroxidation and 8-oxo-dG lesion levels in the hypoxia-exposed mice. Consumption of malondialdehyde also induced a significant increase in the PyMT tumour DNA lesion levels, however, these increases did not translate into enhanced tumour progression. We further showed that the in vivo exposure to acute hypoxia induced accumulation of F4/80 positive tumour-associated macrophages (TAMs), demonstrating a relationship between hypoxia and macrophages in an experimental model.

Conclusion

These data suggest that although exposure to acute hypoxia causes an increase in 8-oxo-dG lesions and TAMs in the PyMT tumours, these increases do not translate into significant changes in tumour progression at the primary or metastatic levels in this strong viral oncogene-driven breast cancer model.

New developments in ancient genomics

Ancient DNA research is on the crest of a 'third wave' of progress due to the introduction of a new generation of DNA sequencing technologies. Here we review the advantages and disadvantages of the four new DNA sequencers that are becoming available to researchers. These machines now allow the recovery of orders of magnitude more DNA sequence data, albeit as short sequence reads. Hence, the potential reassembly of complete ancient genomes seems imminent, and when used to screen libraries of ancient sequences, these methods are cost effective. This new wealth of data is also likely to herald investigations into the functional properties of extinct genes and gene complexes and will improve our understanding of the biological basis of extinct phenotypes.

Human and mouse introns are linked to the same processes and functions through each genome's most frequent non-conserved motifs

We identified the most frequent, variable-length DNA sequence motifs in the human and mouse genomes and sub-selected those with multiple recurrences in the intergenic and intronic regions and at least one additional exonic instance in the corresponding genome. We discovered that these motifs have virtually no overlap with intronic sequences that are conserved between human and mouse, and thus are genome-specific. Moreover, we found that these motifs span a substantial fraction of previously uncharacterized human and mouse intronic space. Surprisingly, we found that these genome-specific motifs are over-represented in the introns of genes belonging to the same biological processes and molecular functions in both the human and mouse genomes even though the underlying sequences are not conserved between the two genomes. In fact, the processes and functions that are linked to these genome-specific sequence-motifs are distinct from the processes and functions which are associated with intronic regions that are conserved between human and mouse. The findings show that intronic regions from different genomes are linked to the same processes and functions in the absence of underlying sequence conservation. We highlight the ramifications of this observation with a concrete example that involves the microsatellite instability gene MLH1.

Eukaryotic transcription factor binding sites--modeling and integrative search methods

A comprehensive knowledge of transcription factor binding sites (TFBS) is important for a mechanistic understanding of transcriptional regulation as well as for inferring gene regulatory networks. Because the DNA motif recognized by a transcription factor is typically short and degenerate, computational approaches for identifying binding sites based only on the sequence motif inevitably suffer from high error rates. Current state-of-the-art techniques for improving computational identification of binding sites can be broadly categorized into two classes: (1) approaches that aim to improve binding motif models by extracting maximal sequence information from experimentally determined binding sites and (2) approaches that supplement binding motif models with additional genomic or other attributes (such as evolutionary conservation). In this review we will discuss recent attempts to improve computational identification of TFBS through these two types of approaches and conclude with thoughts on future development.

Knockout of cytochrome P450 3A yields new mouse models for understanding xenobiotic metabolism

Cytochrome P450 3A (CYP3A) enzymes constitute an important detoxification system that contributes to primary metabolism of more than half of all prescribed medications. To investigate the physiological and pharmacological roles of CYP3A, we generated Cyp3a-knockout (Cyp3a-/-) mice lacking all functional Cyp3a genes. Cyp3a-/- mice were viable, fertile, and without marked physiological abnormalities. However, these mice exhibited severely impaired detoxification capacity when exposed to the chemotherapeutic agent docetaxel, displaying higher exposure levels in response to both oral and intravenous administration. These mice also demonstrated increased sensitivity to docetaxel toxicity, suggesting a primary role for Cyp3a in xenobiotic detoxification. To determine the relative importance of intestinal versus hepatic Cyp3a in first-pass metabolism, we generated transgenic Cyp3a-/- mice expressing human CYP3A4 in either the intestine or the liver. Expression of CYP3A4 in the intestine dramatically decreased absorption of docetaxel into the bloodstream, while hepatic expression aided systemic docetaxel clearance. These results suggest that CYP3A expression determines impairment of drug absorption and efficient systemic clearance in a tissue-specific manner. The genetic models used in this study provide powerful tools to further study CYP3A-mediated xenobiotic metabolism, as well as interactions between CYP3A and other detoxification systems.

Organization of conserved elements near key developmental regulators in vertebrate genomes

Sequence conservation has traditionally been used as a means to target functional regions of complex genomes. In addition to its use in identifying coding regions of genes, the recent availability of whole genome data for a number of vertebrates has permitted high-resolution analyses of the noncoding "dark matter" of the genome. This has resulted in the identification of a large number of highly conserved sequence elements that appear to be preserved in all bony vertebrates. Further positional analysis of these conserved noncoding elements (CNEs) in the genome demonstrates that they cluster around genes involved in developmental regulation. This chapter describes the identification and characterization of these elements, with particular reference to their composition and organization.

Long-range gene control and genetic disease

The past two decades have seen great progress in the elucidation of the genetic basis of human genetic disease. Many clinical phenotypes have been linked with mutations or deletions in specific causative genes. However, it is often less recognized that in addition to the integrity of the protein-coding sequences, human health critically also depends on the spatially, temporally, and quantitatively correct expression of those genes. Genetic disease can therefore equally be caused by disruption of the regulatory mechanisms that ensure proper gene expression. The term "position effect" is used in those situations where the expression level of a gene is deleteriously affected by an alteration in its chromosomal environment, while maintaining an intact transcription unit. Here, we review recent advances in our understanding of the possible mechanisms of a number of "position effect" disease cases and discuss the findings with respect to current models for genome organization and long-range control of gene expression.

Scaffold/matrix attachment regions (S/MARs): relevance for disease and therapy

There is increasing awareness that processes, such as development, aging and cancer, are governed, to a considerable extent, by epigenetic processes, such as DNA and histone modifications. The sites of these modifications in turn reflect their position and role in the nuclear architecture. Since epigenetic changes are easier to reverse than mutations, drugs that remove or add the chemical tags are at the forefront of research for the treatment of cancerous and inflammatory diseases. This review will use selected examples to develop a unified view that might assist the systematic development of novel therapeutic regimens.

Annotating noncoding RNA genes

Noncoding RNA genes produce a functional RNA product rather than a translated protein. More than 1500 homologs of known "classical" RNA genes can be annotated in the human genome sequence, and automatic homology-based methods predict up to 5000 related sequences. Methods to predict novel RNA genes on a whole-genome scale are immature at present, but their use hints at tens of thousands of such genes in the human genome. Messenger RNA-like transcripts with no protein-coding potential are routinely discovered by high-throughput transcriptome analyses. Meanwhile, various experimental studies have suggested that the vast majority of the human genome is transcribed, although the proportion of the detected RNAs that is functional remains unknown.

TREMOR--a tool for retrieving transcriptional modules by incorporating motif covariance

A transcriptional module (TM) is a collection of transcription factors (TF) that as a group, co-regulate multiple, functionally related genes. The task of identifying TMs poses an important biological challenge. Since TFs belong to evolutionarily and structurally related families, TF family members often bind to similar DNA motifs and can confound sequence-based approaches to TM identification. A previous approach to TM detection addresses this issue by pre-selecting a single representative from each TF family. One problem with this approach is that closely related transcription factors can still target sufficiently distinct genes in a biologically meaningful way, and thus, pre-selecting a single family representative may in principle miss certain TMs. Here we report a method—TREMOR (Transcriptional Regulatory Module Retriever). This method uses the Mahalanobis distance to assess the validity of a TM and automatically incorporates the inter-TF binding similarity without resorting to pre-selecting family representatives. The application of TREMOR on human muscle-specific, liver-specific and cell-cycle-related genes reveals TFs and TMs that were validated from literature and also reveals additional related genes.

Structural variation in the human genome: the impact of copy number variants on clinical diagnosis

Over the past few years, the application of whole-genome scanning array technologies has catalyzed the appreciation of a new form of submicroscopic genomic imbalances, referred to as copy number variants. Copy number variants contribute substantially to genetic diversity and result from gains and losses of genomic regions that are 1000 base pairs in size or larger, sometimes encompassing millions of bases of contiguous DNA sequences. As genome-wide scanning techniques become more widely used in diagnostic laboratories, a major challenge is how to accurately interpret which submicroscopic genomic imbalances are pathogenic in nature and which are benign. Herein, we review the literature from the past 3 years on this new source of genomic variability and comment on factors that should be considered when trying to differentiate between a pathogenic and a benign copy number variant.

Mutational and selective effects on copy-number variants in the human genome

Comprehensive descriptions of large insertion/deletion or segmental duplication polymorphisms (SDs) in the human genome have recently been generated. These annotations, known collectively as structural or copy-number variants (CNVs), include thousands of discrete genomic regions and span hundreds of millions of nucleotides. Here we review the genomic distribution of CNVs, which is strongly correlated with gene, repeat and segmental duplication content. We explore the evolutionary mechanisms giving rise to this nonrandom distribution, considering the available data on both human polymorphisms and the fixed changes that differentiate humans from other species. It is likely that mutational biases, selective effects and interactions between these forces all contribute substantially to the spectrum of human copy-number variation. Although defining these variants with nucleotide-level precision remains a largely unmet but critical challenge, our understanding of their potential medical impact and evolutionary importance is rapidly emerging.

Ultraconserved regions encoding ncRNAs are altered in human leukemias and carcinomas

Noncoding RNA (ncRNA) transcripts are thought to be involved in human tumorigenesis. We report that a large fraction of genomic ultraconserved regions (UCRs) encode a particular set of ncRNAs whose expression is altered in human cancers. Genome-wide profiling revealed that UCRs have distinct signatures in human leukemias and carcinomas. UCRs are frequently located at fragile sites and genomic regions involved in cancers. We identified certain UCRs whose expression may be regulated by microRNAs abnormally expressed in human chronic lymphocytic leukemia, and we proved that the inhibition of an overexpressed UCR induces apoptosis in colon cancer cells. Our findings argue that ncRNAs and interaction between noncoding genes are involved in tumorigenesis to a greater extent than previously thought.

Wide-scale analysis of human functional transcription factor binding reveals a strong bias towards the transcription start site

Background

Transcription factors (TF) regulate expression by binding to specific DNA sequences. A binding event is functional when it affects gene expression. Functionality of a binding site is reflected in conservation of the binding sequence during evolution and in over represented binding in gene groups with coherent biological functions. Functionality is governed by several parameters such as the TF-DNA binding strength, distance of the binding site from the transcription start site (TSS), DNA packing, and more. Understanding how these parameters control functionality of different TFs in different biological contexts is a must for identifying functional TF binding sites and for understanding regulation of transcription.

Methodology/Principal Findings

We introduce a novel method to screen the promoters of a set of genes with shared biological function (obtained from the functional Gene Ontology (GO) classification) against a precompiled library of motifs, and find those motifs which are statistically over-represented in the gene set. More than 8000 human (and 23,000 mouse) genes, were assigned to one of 134 GO sets. Their promoters were searched (from 200 bp downstream to 1000 bp upstream the TSS) for 414 known DNA motifs. We optimized the sequence similarity score threshold, independently for every location window, taking into account nucleotide heterogeneity along the promoters of the target genes. The method, combined with binding sequence and location conservation between human and mouse, identifies with high probability functional binding sites for groups of functionally-related genes. We found many location-sensitive functional binding events and showed that they clustered close to the TSS. Our method and findings were tested experimentally.

Conclusions/Significance

We identified reliably functional TF binding sites. This is an essential step towards constructing regulatory networks. The promoter region proximal to the TSS is of central importance for regulation of transcription in human and mouse, just as it is in bacteria and yeast.

Finding cis-regulatory elements using comparative genomics: some lessons from ENCODE data

Identification of functional genomic regions using interspecies comparison will be most effective when the full span of relationships between genomic function and evolutionary constraint are utilized. We find that sets of putative transcriptional regulatory sequences, defined by ENCODE experimental data, have a wide span of evolutionary histories, ranging from stringent constraint shown by deep phylogenetic comparisons to recent selection on lineage-specific elements. This diversity of evolutionary histories can be captured, at least in part, by the suite of available comparative genomics tools, especially after correction for regional differences in the neutral substitution rate. Putative transcriptional regulatory regions show alignability in different clades, and the genes associated with them are enriched for distinct functions. Some of the putative regulatory regions show evidence for recent selection, including a primate-specific, distal promoter that may play a novel role in regulation.

Widely distributed noncoding purifying selection in the human genome

It is widely assumed that human noncoding sequences comprise a substantial reservoir for functional variants impacting gene regulation and other chromosomal processes. Evolutionarily conserved noncoding sequences (CNSs) in the human genome have attracted considerable attention for their potential to simplify the search for functional elements and phenotypically important human alleles. A major outstanding question is whether functionally significant human noncoding variation is concentrated in CNSs or distributed more broadly across the genome. Here, we combine whole genome sequence data from four nonhuman species (chimp, dog, mouse, and rat) with recently available comprehensive human polymorphism data to analyze selection at single-nucleotide resolution. We show that a substantial fraction of active purifying selection in human noncoding sequences occurs outside of CNSs and is diffusely distributed across the genome. This finding suggests the existence of a large complement of human noncoding variants that may impact gene expression and phenotypic traits, the majority of which will escape detection with current approaches to genome analysis.

Diploid genome reconstruction of Ciona intestinalis and comparative analysis with Ciona savignyi

One of the main goals in genome sequencing projects is to determine a haploid consensus sequence even when clone libraries are constructed from homologous chromosomes. However, it has been noticed that haplotypes can be inferred from genome assemblies by investigating phase conservation in sequenced reads. In this study, we seek to infer haplotypes, a diploid consensus sequence, from the genome assembly of an organism, Ciona intestinalis. The Ciona intestinalis genome is an ideal resource from which haplotypes can be inferred because of the high polymorphism rate (1.2%). The haplotype estimation scheme consists of polymorphism detection and phase estimation. The core step of our method is a Gibbs sampling procedure. The mate-pair information from two-end sequenced clone inserts is exploited to provide long-range continuity. We estimate the polymorphism rate of Ciona intestinalis to be 1.2% and 1.5%, according to two different polymorphism counting schemes. The distribution of heterozygosity number is well fit by a compound Poisson distribution. The N50 length of haplotype segments is 37.9 kb in our assembly, while the N50 scaffold length of the Ciona intestinalis assembly is 190 kb. We also infer diploid gene sequences from haplotype segments. According to our reconstruction, 85.4% of predicted gene sequences are continuously covered by single haplotype segments. Our results indicate 97% accuracy in haplotype estimation, based on a simulated data set. We conduct a comparative analysis with Ciona savignyi, and discover interesting patterns of conserved DNA elements in chordates.

Shuffling of cis-regulatory elements is a pervasive feature of the vertebrate lineage

Alignment of orthologous vertebrate loci reveals that a significant proportion of conserved cis-regulatory elements have undergone shuffling during evolution.

Background

All vertebrates share a remarkable degree of similarity in their development as well as in the basic functions of their cells. Despite this, attempts at unearthing genome-wide regulatory elements conserved throughout the vertebrate lineage using BLAST-like approaches have thus far detected noncoding conservation in only a few hundred genes, mostly associated with regulation of transcription and development.

Results

We used a unique combination of tools to obtain regional global-local alignments of orthologous loci. This approach takes into account shuffling of regulatory regions that are likely to occur over evolutionary distances greater than those separating mammalian genomes. This approach revealed one order of magnitude more vertebrate conserved elements than was previously reported in over 2,000 genes, including a high number of genes found in the membrane and extracellular regions. Our analysis revealed that 72% of the elements identified have undergone shuffling. We tested the ability of the elements identified to enhance transcription in zebrafish embryos and compared their activity with a set of control fragments. We found that more than 80% of the elements tested were able to enhance transcription significantly, prevalently in a tissue-restricted manner corresponding to the expression domain of the neighboring gene.

Conclusion

Our work elucidates the importance of shuffling in the detection of cis-regulatory elements. It also elucidates how similarities across the vertebrate lineage, which go well beyond development, can be explained not only within the realm of coding genes but also in that of the sequences that ultimately govern their expression.

Intrinsic promoter activities of primary DNA sequences in the human genome

In order to understand an overview of promoter activities intrinsic to primary DNA sequences in the human genome within a particular cell type, we carried out systematic quantitative luciferase assays of DNA fragments corresponding to putative promoters for 472 human genes which are expressed in HEK (human embryonic kidney epithelial) 293 cells. We observed the promoter activities of them were distributed in a bimodal manner; putative promoters belonging to the first group (with strong promoter activities) were designated as P1 and the latter (with weak promoter activities) as P2. The frequencies of the TATA-boxes, the CpG islands, and the overall G + C-contents were significantly different between these two populations, indicating there are two separate groups of promoters. Interestingly, similar analysis using 251 randomly isolated genomic DNA fragments showed that P2-type promoter occasionally occurs within the human genome. Furthermore, 35 DNA fragments corresponding to putative promoters of non-protein-coding transcripts (ncRNAs) shared similar features with the P2 in both promoter activities and sequence compositions. At least, a part of ncRNAs, which have been massively identified by full-length cDNA projects with no functional relevance inferred, may have originated from those sporadic promoter activities of primary DNA sequences inherent to the human genome.

Position and distance specificity are important determinants of cis-regulatory motifs in addition to evolutionary conservation

Computational discovery of cis-regulatory elements remains challenging. To cope with the high false positives, evolutionary conservation is routinely used. However, conservation is only one of the attributes of cis-regulatory elements and is neither necessary nor sufficient. Here, we assess two additional attributes—positional and inter-motif distance specificity—that are critical for interactions between transcription factors. We first show that for a greater than expected fraction of known motifs, the genes that contain the motifs in their promoters in a position-specific or distance-specific manner are related, both in function and/or in expression pattern. We then use the position and distance specificity to discover novel motifs. Our work highlights the importance of distance and position specificity, in addition to the evolutionary conservation, in discovering cis-regulatory motifs.

Enhancer identification through comparative genomics

With the availability of genomic sequence from numerous vertebrates, a paradigm shift has occurred in the identification of distant-acting gene regulatory elements. In contrast to traditional gene-centric studies in which investigators randomly scanned genomic fragments that flank genes of interest in functional assays, the modern approach begins electronically with publicly available comparative sequence datasets that provide investigators with prioritized lists of putative functional sequences based on their evolutionary conservation. However, although a large number of tools and resources are now available, application of comparative genomic approaches remains far from trivial. In particular, it requires users to dynamically consider the species and methods for comparison depending on the specific biological question under investigation. While there is currently no single general rule to this end, it is clear that when applied appropriately, comparative genomic approaches exponentially increase our power in generating biological hypotheses for subsequent experimental testing. It is anticipated that cardiac-related genes and the identification of their distant-acting transcriptional enhancers are particularly poised to benefit from these modern capabilities.

Cre/loxP-mediated chromosome engineering of the mouse genome

Together with numerous other genome modifications, chromosome engineering offers a very powerful tool to accelerate the functional analysis of the mammalian genome. The technology, based on the Cre/loxP system, is used more and more in the scientific community in order to generate new chromosomes carrying deletions, duplications, inversions and translocations in targeted regions of interest. In this review, we will present the basic principle of the technique either in vivo or in vitro and we will briefly describe some applications to provide highly valuable genetic tools, to decipher the mammalian genome organisation and to analyze human diseases in the mouse.

How repetitive are genomes?

Background

Genome sequences vary strongly in their repetitiveness and the causes for this are still debated. Here we propose a novel measure of genome repetitiveness, the index of repetitiveness, I_r, which can be computed in time proportional to the length of the sequences analyzed. We apply it to 336 genomes from all three domains of life.

Results

The expected value of I_r is zero for random sequences of any G/C content and greater than zero for sequences with excess repeats. We find that the I_r of archaea is significantly smaller than that of eubacteria, which in turn is smaller than that of eukaryotes. Mouse chromosomes have a significantly higher I_r than human chromosomes and within each genome the Y chromosome is most repetitive. A sliding window analysis reveals that the human HOXA cluster and two surrounding genes are characterized by local minima in I_r. A program for calculating the I_r is freely available at .

Conclusion

The general measure of DNA repetitiveness proposed in this paper can be efficiently computed on a genomic scale. This reveals a broad spectrum of repetitiveness among diverse genomes which agrees qualitatively with previous studies of repeat content. A sliding window analysis helps to analyze the intragenomic distribution of repeats.

New techniques to understand chromosome dosage: mouse models of aneuploidy

Aberrations in human chromosome copy number and structure are common and extremely deleterious. Their downstream effects on phenotype are caused by aberrant dosage of sequences in the affected regions. However, we know little about why the abnormal gene copy number causes disease or why specific features result from deficits in specific chromosomes. Mice are the organism of choice to help us try to tease apart the complex relationships between genotype and phenotype in aneuploidy and segmental aneusomy syndromes. As new technologies such as chromosome engineering and the creation of transchromosomic mice become routine, these will help us identify individual dosage-sensitive genes that are causative in specific syndromes and will enable us to produce mouse models to accurately recapitulate human chromosomal disorders.

Parathyroid hormone controls receptor activator of NF-kappaB ligand gene expression via a distant transcriptional enhancer

RANKL, a protein essential for osteoclast development and survival, is stimulated by parathyroid hormone (PTH) via a PTH receptor 1/cyclic AMP (cAMP)/protein kinase A (PKA)/CREB cascade, exclusively in osteoblastic cells. We report that a bacterial artificial chromosome-based transcriptional reporter construct containing 120 kb of RANKL 5'-flanking region was stimulated by dibutyryl-cAMP in stromal/osteoblastic cells, but not other cell types. Full cAMP responsiveness was dependent upon a conserved 715-bp region located 76 kb upstream from the transcription start site, which we identified by sequential deletion analysis and by comparison of human and mouse genomic sequences in silico. This region contained conserved consensus sequences which bound CREB and the osteoblast-specific transcription factor Runx2, and when mutated blunted cAMP responsiveness. Overexpression of Runx2 potentiated cAMP responsiveness of the endogenous RANKL gene in a cell-type-specific manner. Lastly, PTH responsiveness of the endogenous RANKL gene was abrogated in mice from which we deleted this conserved upstream region. Thus, PTH responsiveness of the RANKL gene is determined by a distant regulatory region that responds to cAMP in a cell-type-specific manner and Runx2 may contribute to such cell-type specificity.