Stochastic segmentation models for array-based comparative genomic hybridization data analysis

Array-based comparative genomic hybridization (array-CGH) is a high throughput, high resolution technique for studying the genetics of cancer. Analysis of array-CGH data typically involves estimation of the underlying chromosome copy numbers from the log fluorescence ratios and segmenting the chromosome into regions with the same copy number at each location. We propose for the analysis of array-CGH data, a new stochastic segmentation model and an associated estimation procedure that has attractive statistical and computational properties. An important benefit of this Bayesian segmentation model is that it yields explicit formulas for posterior means, which can be used to estimate the signal directly without performing segmentation. Other quantities relating to the posterior distribution that are useful for providing confidence assessments of any given segmentation can also be estimated by using our method. We propose an approximation method whose computation time is linear in sequence length which makes our method practically applicable to the new higher density arrays. Simulation studies and applications to real array-CGH data illustrate the advantages of the proposed approach.     

A comparative genomic analysis of the cow,pig, and human CFTR genes identifies potential intronic regulatory elements

The identification of sequences within noncoding regions of genes that are conserved between several species may indicate potential regulatory elements. This is important for genes with complex control mechanisms such as the cystic fibrosis transmembrane conductance regulator (CFTR). CFTR demonstrates similar patterns of temporal and spatial expression in human and sheep, but these differ significantly in mouse cftr. The complete sheep CFTR sequence is unavailable so we annotated BAC clones encompassing the CFTR gene from two other artiodactyl species (cow and pig) for comparative sequence analysis. Regions of introns 2, 3, 10, 17a, 18, and 21 and 3' flanking sequence corresponding to human CFTR DNase I hypersensitive sites (DHS) showed high homology in the cow and pig. Cross-species sequence conservation also enabled finer mapping of other human DHS, including those in introns 1, 16, and 20. Additional potential regulatory elements not associated with human DHS were also identified.     

Functional impact of transposable elements using bioinformatic analysis and a comparative genomic approach

A dual coding event, which is the translation of different isoforms from a single gene, is one of the special patterns among the alternative splicing events. This is an important mechanism for the regulation of protein diversity in human and mouse genomes. Although the regulation for dual coding events has been characterized in a few genes, the individual mechanism remains unclear. Numerous studies have described the exonization of transposable elements, which is the splicing mediated insertion of transposable element sequence fragments into mature mRNAs. Therefore, in this study, we investigated the number of transposable element (TE)-derived dual coding genes in human, chimpanzee and mouse genomes. TE fusion exons appeared in the dual coding regions of 309 human genes. Functional protein domain alterations by TE-derived dual coding events were observed in 129 human genes. Comparative TE-derived dual coding events were also analyzed in chimpanzee and mouse orthologs. Seventy chimpanzee orthologs had TE-derived dual coding events, but mouse orthologs did not have any TE-derived dual coding events. Taken together, our analyses listed the number of TE-derived dual coding genes which could be investigated by experimental analysis and suggested that TE-derived dual coding events were major sources for the functional diversity of human genes, but not mouse genes.     

EMERGE: a flexible modelling framework to predict genomic regulatory elements from genomic signatures

Regulatory DNA elements, short genomic segments that regulate gene expression, have been implicated in developmental disorders and human disease. Despite this clinical urgency, only a small fraction of the regulatory DNA repertoire has been confirmed through reporter gene assays. The overall success rate of functional validation of candidate regulatory elements is low. Moreover, the number and diversity of datasets from which putative regulatory elements can be identified is large and rapidly increasing. We generated a flexible and user-friendly tool to integrate the information from different types of genomic datasets, e.g. ATAC-seq, ChIP-seq, conservation, aiming to increase the ease and success rate of functional prediction. To this end, we developed the EMERGE program that merges all datasets that the user considers informative and uses a logistic regression framework, based on validated functional elements, to set optimal weights to these datasets. ROC curve analysis shows that a combination of datasets leads to improved prediction of tissue-specific enhancers in human, mouse and Drosophila genomes. Functional assays based on this prediction can be expected to have substantially higher success rates. The resulting integrated signal for prediction of functional elements can be plotted in a build-in genome browser or exported for further analysis.     

Nomadic enhancers: tissue-specific cis-regulatory elements of yellow have divergent genomic positions among Drosophila species

cis-regulatory DNA sequences known as enhancers control gene expression in space and time. They are central to metazoan development and are often responsible for changes in gene regulation that contribute to phenotypic evolution. Here, we examine the sequence, function, and genomic location of enhancers controlling tissue- and cell-type specific expression of the yellow gene in six Drosophila species. yellow is required for the production of dark pigment, and its expression has evolved largely in concert with divergent pigment patterns. Using Drosophila melanogaster as a transgenic host, we examined the expression of reporter genes in which either 5' intergenic or intronic sequences of yellow from each species controlled the expression of Green Fluorescent Protein. Surprisingly, we found that sequences controlling expression in the wing veins, as well as sequences controlling expression in epidermal cells of the abdomen, thorax, and wing, were located in different genomic regions in different species. By contrast, sequences controlling expression in bristle-associated cells were located in the intron of all species. Differences in the precise pattern of spatial expression within the developing epidermis of D. melanogaster transformants usually correlated with adult pigmentation in the species from which the cis-regulatory sequences were derived, which is consistent with cis-regulatory evolution affecting yellow expression playing a central role in Drosophila pigmentation divergence. Sequence comparisons among species favored a model in which sequential nucleotide substitutions were responsible for the observed changes in cis-regulatory architecture. Taken together, these data demonstrate frequent changes in yellow cis-regulatory architecture among Drosophila species. Similar analyses of other genes, combining in vivo functional tests of enhancer activity with in silico comparative genomics, are needed to determine whether the pattern of regulatory evolution we observed for yellow is characteristic of genes with rapidly evolving expression patterns.     

Discovery of regulatory elements in vertebrates through comparative genomics

We have analyzed issues of reliability in studies in which comparative genomic approaches have been applied to the discovery of regulatory elements at a genome-wide level in vertebrates. We point out some potential problems with such studies, including difficulties in accurately identifying orthologous promoter regions. Many of these subtle analytical problems have become apparent only when studying the more complex vertebrate genomes. By determining motif reliability, we compared existing tools when applied to the discovery of vertebrate regulatory elements. We then used a statistical clustering method to produce a computational catalog of high quality putative regulatory elements from vertebrates, some of which are widely conserved among vertebrates and many of which are novel regulatory elements. The results provide a glimpse into the wealth of information that comparative genomics can yield and suggest the need for further improvement of genome-wide comparative computational techniques.     

Functional analysis of chicken Sox2 enhancers highlights an array of diverse regulatory elements that are conserved in mammals

Sox2 expression marks neural and sensory primordia at various stages of development. A 50 kb genomic region of chicken Sox2 was isolated and scanned for enhancer activity utilizing embryo electroporation, resulting in identification of a battery of enhancers. Although Sox2 expression in the early embryonic CNS appears uniform, it is actually pieced together by five separate enhancers with distinct spatio-temporal specificities, including the one activated by the neural induction signals emanating from Hensen's node. Enhancers for Sox2 expression in the lens and nasal/otic placodes and in the neural crest were also determined. These functionally identified Sox2 enhancers exactly correspond to the extragenic sequence blocks conspicuously conserved between chicken and mammals, which are not discernible by sequence comparison among mammals.