Comparative genomics

Involvement of tdTomato-Tagged RPE cells in a mouse PVR model with enzymatically compromised retina

Ocular trauma and surgery are considered the most common cause for proliferative vitreoretinopathy (PVR). Many retinal cell types are thought to be the cellular source for PVR although most risk factors for PVR are associated with intravitreal dispersion of the retinal pigment epithelium (RPE) cells. Major PVR animal models are rabbit and swine with an artificial implantation of exogenous cells into the vitreous to form epiretinal membrane (ERM) which does not recapitulate a real PVR pathology. To clarify and validate the participation of RPE cells, to mimic ocular trauma in situ, and to reveal the related macromolecule changes in PVR pathology, we utilized a dispase treatment to damage the retina in establishment of a reliable RPE-tagged PVR mouse model with ERM-like tissues formed within and on both surface of the retina. The immunostaining of patient epiretinal membranes with lineage markers confirms RPE is involved in PVR development. Quantitative PCR analysis indicates the dedifferentiation of RPE cells switches RPE from epithelial to mesenchymal phenotype to re-enter a proliferative and mobile state underlying PVR. Gene expression results of the mouse PVR model retinas are consistent with the microarray gene expression profile of human PVR retinas, validating that our mouse PVR model resembles human PVR and is thereby suitable for molecular mechanism and pharmaceutical studies.

Localization and Molecular Cloning of the ASMT Gene for Melatonin Synthesis in Pigs

Melatonin is synthesized in multiple tissues and organs of pigs, and existing studies have shown the presence of the melatonin-synthesizing enzyme ASMT protein. However, the genomic information for the ASMT gene has been lacking. The aim of this study was to locate the genomic information of the ASMT gene in pigs using comparative genomics analysis and then obtain the coding region information through molecular cloning. First, using the NCBI Genome Data Viewer, we found that in most animals, the AKAP17A gene is often located next to the ASMT gene, with both genes arranged in the same direction. Similarly, the P2RY8 gene is commonly adjacent to the ASMTL gene, also in the same orientation. We also discovered that the ASMTL gene is frequently adjacent to the ASMT gene and arranged in the opposite direction. Using the "three-point localization" principle, we inferred the position of the ASMT gene based on the coordinates of AKAP17A and ASMTL in pigs. Our results revealed that on the pig X chromosome, a gene called LOC110258194 is located next to the AKAP17A and ASMTL genes, and its arrangement aligns with the ASMT gene in other species. Additionally, Ensembl contains a gene, ENSSSCG00000032659, at the same position, with completely overlapping exons, though it is not annotated as ASMT. Further analysis using the TreeFam tool from EMBL-EBI and the CDD tool from NCBI revealed that LOC110258194 and ENSSSCG00000032659 do not contain the typical Maf domain of ASMTL and, thus, should not be annotated as ASMTL, but rather as the ASMT gene. Using a slow-down PCR method for high-GC content genes, we successfully cloned the full CDS region of the pig ASMT gene and identified a new transcript missing Exon 6 and Exon 7. This transcript was submitted to NCBI and assigned the GenBank accession number MW847601. Our results represent the first successful localization of the ASMT gene in pigs, the first cloning of the ASMT gene's coding region, and the first discovery of a new transcript of the pig ASMT gene.

Phylogenomics provides insights into the evolution of cactophily and host plant shifts in Drosophila

Cactophilic species of the Drosophila buzzatii cluster (repleta group) comprise an excellent model group to investigate genomic changes underlying adaptation to extreme climate conditions and host plants. In particular, these species form a tractable system to study the transition from chemically simpler breeding sites (like prickly pears of the genus Opuntia) to chemically more complex hosts (columnar cacti). Here, we report four highly contiguous genome assemblies of three species of the buzzatii cluster. Based on this genomic data and inferred phylogenetic relationships, we identified candidate taxonomically restricted genes (TRGs) likely involved in the evolution of cactophily and cactus host specialization. Functional enrichment analyses of TRGs within the buzzatii cluster identified genes involved in detoxification, water preservation, immune system response, anatomical structure development, and morphogenesis. In contrast, processes that regulate responses to stress, as well as the metabolism of nitrogen compounds, transport, and secretion were found in the set of species that are columnar cacti dwellers. These findings are in line with the hypothesis that those genomic changes brought about key mechanisms underlying the adaptation of the buzzatii cluster species to arid regions in South America.

Advanced Paternal Age and Future Generations

Paternal age at conception has been increasing. In this review, we first present the results from the major mammalian animal models used to establish that increasing paternal age does affect progeny outcome. These models provide several major advantages including the possibility to assess multi- transgenerational effects of paternal age on progeny in a relatively short time window. We then present the clinical observations relating advanced paternal age to fertility and effects on offspring with respect to perinatal health, cancer risk, genetic diseases, and neurodevelopmental effects. An overview of the potential mechanism operating in altering germ cells in advanced age is presented. This is followed by an analysis of the current state of management of reproductive risks associated with advanced paternal age. The numerous challenges associated with developing effective, practical strategies to mitigate the impact of advanced paternal age are outlined along with an approach on how to move forward with this important clinical quandary.

Selection on Accessible Chromatin Regions in Capsella grandiflora

Accurate estimates of genome-wide rates and fitness effects of new mutations are essential for an improved understanding of molecular evolutionary processes. Although eukaryotic genomes generally contain a large noncoding fraction, functional noncoding regions and fitness effects of mutations in such regions are still incompletely characterized. A promising approach to characterize functional noncoding regions relies on identifying accessible chromatin regions (ACRs) tightly associated with regulatory DNA. Here, we applied this approach to identify and estimate selection on ACRs in Capsella grandiflora, a crucifer species ideal for population genomic quantification of selection due to its favorable population demography. We describe a population-wide ACR distribution based on ATAC-seq data for leaf samples of 16 individuals from a natural population. We use population genomic methods to estimate fitness effects and proportions of positively selected fixations (α) in ACRs and find that intergenic ACRs harbor a considerable fraction of weakly deleterious new mutations, as well as a significantly higher proportion of strongly deleterious mutations than comparable inaccessible intergenic regions. ACRs are enriched for expression quantitative trait loci (eQTL) and depleted of transposable element insertions, as expected if intergenic ACRs are under selection because they harbor regulatory regions. By integrating empirical identification of intergenic ACRs with analyses of eQTL and population genomic analyses of selection, we demonstrate that intergenic regulatory regions are an important source of nearly neutral mutations. These results improve our understanding of selection on noncoding regions and the role of nearly neutral mutations for evolutionary processes in outcrossing Brassicaceae species.

Chromatography-Free Purification Strategies for Large Biological Macromolecular Complexes Involving Fractionated PEG Precipitation and Density Gradients

A complex interplay between several biological macromolecules maintains cellular homeostasis. Generally, the demanding chemical reactions which sustain life are not performed by individual macromolecules, but rather by several proteins that together form a macromolecular complex. Understanding the functional interactions amongst subunits of these macromolecular machines is fundamental to elucidate mechanisms by which they maintain homeostasis. As the faithful function of macromolecular complexes is essential for cell survival, their mis-function leads to the development of human diseases. Furthermore, detailed mechanistic interrogation of the function of macromolecular machines can be exploited to develop and optimize biotechnological processes. The purification of intact macromolecular complexes is an essential prerequisite for this; however, chromatographic purification schemes can induce the dissociation of subunits or the disintegration of the whole complex. Here, we discuss the development and application of chromatography-free purification strategies based on fractionated PEG precipitation and orthogonal density gradient centrifugation that overcomes existing limitations of established chromatographic purification protocols. The presented case studies illustrate the capabilities of these procedures for the purification of macromolecular complexes.

Pan-Genome Analysis Reveals Host-Specific Functional Divergences in Burkholderia gladioli

Burkholderia gladioli has high versatility and adaptability to various ecological niches. Here, we constructed a pan-genome using 14 genome sequences of B. gladioli, which originate from different niches, including gladiolus, rice, humans, and nature. Functional roles of core and niche-associated genomes were investigated by pathway enrichment analyses. Consequently, we inferred the uniquely important role of niche-associated genomes in (1) selenium availability during competition with gladiolus host; (2) aromatic compound degradation in seed-borne and crude oil-accumulated environments, and (3) stress-induced DNA repair system/recombination in the cystic fibrosis-niche. We also identified the conservation of the rhizomide biosynthetic gene cluster in all the B. gladioli strains and the concentrated distribution of this cluster in human isolates. It was confirmed the absence of complete CRISPR/Cas system in both plant and human pathogenic B. gladioli and the presence of the system in B. gladioli living in nature, possibly reflecting the inverse relationship between CRISPR/Cas system and virulence.

Haplotype-resolved genome assembly enables gene discovery in the red palm weevil Rhynchophorus ferrugineus

The red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is an economically-important invasive species that attacks multiple species of palm trees around the world. A better understanding of gene content and function in R. ferrugineus has the potential to inform pest control strategies and thereby mitigate economic and biodiversity losses caused by this species. Using 10x Genomics linked-read sequencing, we produced a haplotype-resolved diploid genome assembly for R. ferrugineus from a single heterozygous individual with modest sequencing coverage (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼ 62x). Benchmarking against conserved single-copy Arthropod orthologs suggests both pseudo-haplotypes in our R. ferrugineus genome assembly are highly complete with respect to gene content, and do not suffer from haplotype-induced duplication artifacts present in a recently published hybrid assembly for this species. Annotation of the larger pseudo-haplotype in our assembly provides evidence for 23,413 protein-coding loci in R. ferrugineus, including over 13,000 predicted proteins annotated with Gene Ontology terms and over 6000 loci independently supported by high-quality Iso-Seq transcriptomic data. Our assembly also includes 95% of R. ferrugineus chemosensory, detoxification and neuropeptide-related transcripts identified previously using RNA-seq transcriptomic data, and provides a platform for the molecular analysis of these and other functionally-relevant genes that can help guide management of this widespread insect pest.

Comparative Genomics, from the Annotated Genome to Valuable Biological Information: A Case Study

High availability of fast, cheap, and high-throughput next generation sequencing techniques resulted in acquisition of numerous de novo sequenced and assembled bacterial genomes. It rapidly became clear that digging out useful biological information from such a huge amount of data presents a considerable challenge. In this chapter we share our experience with utilization of several handy open source comparative genomic tools. All of them were applied in the studies focused on revealing inter- and intraspecies variation in pectinolytic plant pathogenic bacteria classified to Dickeya solani and Pectobacterium parmentieri. As the described software performed well on the species within the Pectobacteriaceae family, it presumably may be readily utilized on some closely related taxa from the Enterobacteriaceae family. First of all, implementation of various annotation software is discussed and compared. Then, tools computing whole genome comparisons including generation of circular juxtapositions of multiple sequences, revealing the order of synteny blocks or calculation of ANI or Tetra values are presented. Besides, web servers intended either for functional annotation of the genes of interest or for detection of genomic islands, plasmids, prophages, CRISPR/Cas are described. Last but not least, utilization of the software designed for pangenome studies and the further downstream analyses is explained. The presented work not only summarizes broad possibilities assured by the comparative genomic approach but also provides a user-friendly guide that might be easily followed by nonbioinformaticians interested in undertaking similar studies.

Bioinformatics as a Tool for the Structural and Evolutionary Analysis of Proteins

This chapter deals with the topic of bioinformatics, computational, mathematics, and statistics tools applied to biology, essential for the analysis and characterization of biological molecules, in particular proteins, which play an important role in all cellular and evolutionary processes of the organisms. In recent decades, with the next generation sequencing technologies and bioinformatics, it has facilitated the collection and analysis of a large amount of genomic, transcriptomic, proteomic, and metabolomic data from different organisms that have allowed predictions on the regulation of expression, transcription, translation, structure, and mechanisms of action of proteins as well as homology, mutations, and evolutionary processes that generate structural and functional changes over time. Although the information in the databases is greater every day, all bioinformatics tools continue to be constantly modified to improve performance that leads to more accurate predictions regarding protein functionality, which is why bioinformatics research remains a great challenge.

Unique genomic traits for cold adaptation in Naganishia vishniacii, a polyextremophile yeast isolated from Antarctica

Cold environments impose challenges to organisms. Polyextremophile microorganisms can survive in these conditions thanks to an array of counteracting mechanisms. Naganishia vishniacii, a yeast species hitherto only isolated from McMurdo Dry Valleys, Antarctica, is an example of a polyextremophile. Here we present the first draft genomic sequence of N. vishniacii. Using comparative genomics, we unraveled unique characteristics of cold associated adaptations. 336 putative genes (total: 6183) encoding solute transfers and chaperones, among others, were absent in sister species. Among genes shared by N. vishniacii and its closest related species we found orthologs encompassing possible evidence of positive selection (dN/dS > 1). Genes associated with photoprotection were found in agreement with high solar irradiation exposure. Also genes coding for desaturases and genomic features associated with cold tolerance (i.e. trehalose synthesis and lipid metabolism) were explored. Finally, biases in amino acid usage (namely an enrichment of glutamine and a trend in proline reduction) were observed, possibly conferring increased protein flexibility. To the best of our knowledge, such a combination of mechanisms for cold tolerance has not been previously reported in fungi, making N. vishniacii a unique model for the study of the genetic basis and evolution of cold adaptation strategies.

Endogenous or Exogenous Retinal Pigment Epithelial Cells: A Comparison of Two Experimental Animal Models of Proliferative Vitreoretinopathy

Purpose

Proliferative vitreoretinopathy (PVR) is a blinding condition that can occur following ocular penetrating injury and retinal detachment. To develop effective therapeutics for PVR, it is imperative to establish an animal model that is reproducible, closest in anatomy to the human eye, and most representative of the human disease. We compared two in vivo models of PVR in minipig eyes to assess reproducibility and consistency.

Methods

Six minipigs underwent PVR induction with procedure A and six underwent procedure B. In both procedures, PVR was induced with vitrectomy, bleb retinal detachment, retinotomy, and injection of platelet-rich plasma. In procedure A, retinal pigment epithelial (RPE) cells were harvested from cadaveric pig eyes and injected at the end of surgery. In procedure B, native RPE cells were released into the vitreous cavity by creating a RPE detachment and scraping the RPE layer. PVR severity was graded on fundoscopic examination with a modified Silicone Study Classification System for PVR. Severe PVR was defined as stages 2 to 5.

Results

Three eyes (50%) and five eyes (83.3%) developed re-detachment of the retina from severe PVR in procedures A and B, respectively (P = 0.55). Median PVR stage was higher in eyes that underwent procedure B compared to eyes that underwent procedure A, although the difference was not statistically significant (2.5 vs. 1.5, P = 0.26).

Conclusions

This new model utilizing native RPE cells achieved a high consistency in inducing severe PVR in the minipig.

Translational Relevance

Our model closely follows pathogenic events in human PVR, making it ideal for preclinical testing of novel therapeutics for PVR.

Bioinformatics of Metalloproteins and Metalloproteomes

Trace metals are inorganic elements that are required for all organisms in very low quantities. They serve as cofactors and activators of metalloproteins involved in a variety of key cellular processes. While substantial effort has been made in experimental characterization of metalloproteins and their functions, the application of bioinformatics in the research of metalloproteins and metalloproteomes is still limited. In the last few years, computational prediction and comparative genomics of metalloprotein genes have arisen, which provide significant insights into their distribution, function, and evolution in nature. This review aims to offer an overview of recent advances in bioinformatic analysis of metalloproteins, mainly focusing on metalloprotein prediction and the use of different metals across the tree of life. We describe current computational approaches for the identification of metalloprotein genes and metal-binding sites/patterns in proteins, and then introduce a set of related databases. Furthermore, we discuss the latest research progress in comparative genomics of several important metals in both prokaryotes and eukaryotes, which demonstrates divergent and dynamic evolutionary patterns of different metalloprotein families and metalloproteomes. Overall, bioinformatic studies of metalloproteins provide a foundation for systematic understanding of trace metal utilization in all three domains of life.

Population genetic models of GERP scores suggest pervasive turnover of constrained sites across mammalian evolution

Comparative genomic approaches have been used to identify sites where mutations are under purifying selection and of functional consequence by searching for sequences that are conserved across distantly related species. However, the performance of these approaches has not been rigorously evaluated under population genetic models. Further, short-lived functional elements may not leave a footprint of sequence conservation across many species. We use simulations to study how one measure of conservation, the Genomic Evolutionary Rate Profiling (GERP) score, relates to the strength of selection (N_es). We show that the GERP score is related to the strength of purifying selection. However, changes in selection coefficients or functional elements over time (i.e. functional turnover) can strongly affect the GERP distribution, leading to unexpected relationships between GERP and N_es. Further, we show that for functional elements that have a high turnover rate, adding more species to the analysis does not necessarily increase statistical power. Finally, we use the distribution of GERP scores across the human genome to compare models with and without turnover of sites where mutations are under purifying selection. We show that mutations in 4.51% of the noncoding human genome are under purifying selection and that most of this sequence has likely experienced changes in selection coefficients throughout mammalian evolution. Our work reveals limitations to using comparative genomic approaches to identify deleterious mutations. Commonly used GERP score thresholds miss over half of the noncoding sites in the human genome where mutations are under purifying selection.

One of the most significant and challenging tasks in modern genomics is to assess the functional consequences of a particular nucleotide change in a genome. A common approach to address this challenge prioritizes sequences that share similar nucleotides across distantly related species, with the rationale that mutations at such positions were deleterious and removed from the population by purifying natural selection. Our manuscript shows that one popular measure of sequence conservation, the GERP score, performs well at identifying selected mutations if mutations at a site were under selection across all of mammalian evolution. Changes in selection at a given site dramatically reduces the power of GERP to detect selected mutations in humans. We also combine population genetic models with the distribution of GERP scores at noncoding sites across the human genome to show that the degree of selection at individual sites has changed throughout mammalian evolution. Importantly, we demonstrate that at least 80 Mb of noncoding sequence under purifying selection in humans will not have extreme GERP scores and will likely be missed by modern comparative genomic approaches. Our work argues that new approaches, potentially based on genetic variation within species, will be required to identify deleterious mutations.

A genome alignment of 120 mammals highlights ultraconserved element variability and placenta-associated enhancers

BACKGROUND

Multiple alignments of mammalian genomes have been the basis of many comparative genomic studies aiming at annotating genes, detecting regions under evolutionary constraint, and studying genome evolution. A key factor that affects the power of comparative analyses is the number of species included in a genome alignment.

RESULTS

To utilize the increased number of sequenced genomes and to provide an accessible resource for genomic studies, we generated a mammalian genome alignment comprising 120 species. We used this alignment and the CESAR method to provide protein-coding gene annotations for 119 non-human mammals. Furthermore, we illustrate the utility of this alignment by 2 exemplary analyses. First, we quantified how variable ultraconserved elements (UCEs) are among placental mammals. Leveraging the high taxonomic coverage in our alignment, we estimate that UCEs contain on average 4.7%-15.6% variable alignment columns. Furthermore, we show that the center regions of UCEs are generally most constrained. Second, we identified enhancer sequences that are only conserved in placental mammals. We found that these enhancers are significantly associated with placenta-related genes, suggesting that some of these enhancers may be involved in the evolution of placental mammal-specific aspects of the placenta.

CONCLUSION

The 120-mammal alignment and all other data are available for analysis and visualization in a genome browser at https://genome-public.pks.mpg.de/and for download at https://bds.mpi-cbg.de/hillerlab/120MammalAlignment/.

Primate optogenetics: Progress and prognosis

Monkeys are a premier model organism for neuroscience research. Activity in the central nervous systems of monkeys can be recorded and manipulated while they perform complex perceptual, motor, or cognitive tasks. Conventional techniques for manipulating neural activity in monkeys are too coarse to address many of the outstanding questions in primate neuroscience, but optogenetics holds the promise to overcome this hurdle. In this article, we review the progress that has been made in primate optogenetics over the past 5 years. We emphasize the use of gene regulatory sequences in viral vectors to target specific neuronal types, and we present data on vectors that we engineered to target parvalbumin-expressing neurons. We conclude with a discussion of the utility of optogenetics for treating sensorimotor hearing loss and Parkinson's disease, areas of translational neuroscience in which monkeys provide unique leverage for basic science and medicine.

Mapping of Human FOXP2 Enhancers Reveals Complex Regulation

Mutations of the FOXP2 gene cause a severe speech and language disorder, providing a molecular window into the neurobiology of language. Individuals with FOXP2 mutations have structural and functional alterations affecting brain circuits that overlap with sites of FOXP2 expression, including regions of the cortex, striatum, and cerebellum. FOXP2 displays complex patterns of expression in the brain, as well as in non-neuronal tissues, suggesting that sophisticated regulatory mechanisms control its spatio-temporal expression. However, to date, little is known about the regulation of FOXP2 or the genomic elements that control its expression. Using chromatin conformation capture (3C), we mapped the human FOXP2 locus to identify putative enhancer regions that engage in long-range interactions with the promoter of this gene. We demonstrate the ability of the identified enhancer regions to drive gene expression. We also show regulation of the FOXP2 promoter and enhancer regions by candidate regulators - FOXP family and TBR1 transcription factors. These data point to regulatory elements that may contribute to the temporal- or tissue-specific expression patterns of human FOXP2. Understanding the upstream regulatory pathways controlling FOXP2 expression will bring new insight into the molecular networks contributing to human language and related disorders.

Comparative and evolutionary insights into CD4 gene across mammalian and avian taxa

The present day genetic architecture of a species bears much significance to its closely related species which is due to species-specific differences, shaped by different evolutionary forces across time scale. With the availability of whole genome sequence of several closely related species, it is now possible to infer evolutionary patterns of genes and genomes in specific lineages. To this respect, CD4 gene, primarily responsible for defensive mechanism in human, is conserved across a few taxa, and thus, comparative genomic studies could be useful for better understanding of host-pathogen biology. Comparative and evolutionary analyses were performed in eleven taxa (10 mammalian and avian) with different statistical algorithms. Phylogenetic inferences revealed recent divergence of human and chimpanzee, and pig was found to be diverged from rest of the taxa significantly. Additionally, gene length, microsatellites, and secondary structures were observed across taxa. The genetic architecture of CD4 gene and its evolutionary history in different mammalian taxa provide crucial evidence in support of the fact that this gene might have been evolving at a similar rate to other human immune system genes. Future population-based study and structural modeling would unravel the differential ability to interact with HIV virus and influence immune system in humans.

Parallel computing in genomic research: advances and applications

Today's genomic experiments have to process the so-called "biological big data" that is now reaching the size of Terabytes and Petabytes. To process this huge amount of data, scientists may require weeks or months if they use their own workstations. Parallelism techniques and high-performance computing (HPC) environments can be applied for reducing the total processing time and to ease the management, treatment, and analyses of this data. However, running bioinformatics experiments in HPC environments such as clouds, grids, clusters, and graphics processing unit requires the expertise from scientists to integrate computational, biological, and mathematical techniques and technologies. Several solutions have already been proposed to allow scientists for processing their genomic experiments using HPC capabilities and parallelism techniques. This article brings a systematic review of literature that surveys the most recently published research involving genomics and parallel computing. Our objective is to gather the main characteristics, benefits, and challenges that can be considered by scientists when running their genomic experiments to benefit from parallelism techniques and HPC capabilities.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: application to zebrafish

Many important model organisms for biomedical and evolutionary research have sequenced genomes, but occupy a phylogenetically isolated position, evolutionarily distant from other sequenced genomes. This phylogenetic isolation is exemplified for zebrafish, a vertebrate model for cis-regulation, development and human disease, whose evolutionary distance to all other currently sequenced fish exceeds the distance between human and chicken. Such large distances make it difficult to align genomes and use them for comparative analysis beyond gene-focused questions. In particular, detecting conserved non-genic elements (CNEs) as promising cis-regulatory elements with biological importance is challenging. Here, we develop a general comparative genomics framework to align isolated genomes and to comprehensively detect CNEs. Our approach integrates highly sensitive and quality-controlled local alignments and uses alignment transitivity and ancestral reconstruction to bridge large evolutionary distances. We apply our framework to zebrafish and demonstrate substantially improved CNE detection and quality compared with previous sets. Our zebrafish CNE set comprises 54 533 CNEs, of which 11 792 (22%) are conserved to human or mouse. Our zebrafish CNEs (http://zebrafish.stanford.edu) are highly enriched in known enhancers and extend existing experimental (ChIP-Seq) sets. The same framework can now be applied to the isolated genomes of frog, amphioxus, Caenorhabditis elegans and many others.

Conservation genomics of threatened animal species

The genomics era has opened up exciting possibilities in the field of conservation biology by enabling genomic analyses of threatened species that previously were limited to model organisms. Next-generation sequencing (NGS) and the collection of genome-wide data allow for more robust studies of the demographic history of populations and adaptive variation associated with fitness and local adaptation. Genomic analyses can also advance management efforts for threatened wild and captive populations by identifying loci contributing to inbreeding depression and disease susceptibility, and predicting fitness consequences of introgression. However, the development of genomic tools in wild species still carries multiple challenges, particularly those associated with computational and sampling constraints. This review provides an overview of the most significant applications of NGS and the implications and limitations of genomic studies in conservation.

Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo)

Background

The turkey (Meleagris gallopavo) is an important agricultural species and the second largest contributor to the world’s poultry meat production. Genetic improvement is attributed largely to selective breeding programs that rely on highly heritable phenotypic traits, such as body size and breast muscle development. Commercial breeding with small effective population sizes and epistasis can result in loss of genetic diversity, which in turn can lead to reduced individual fitness and reduced response to selection. The presence of genomic diversity in domestic livestock species therefore, is of great importance and a prerequisite for rapid and accurate genetic improvement of selected breeds in various environments, as well as to facilitate rapid adaptation to potential changes in breeding goals. Genomic selection requires a large number of genetic markers such as e.g. single nucleotide polymorphisms (SNPs) the most abundant source of genetic variation within the genome.

Results

Alignment of next generation sequencing data of 32 individual turkeys from different populations was used for the discovery of 5.49 million SNPs, which subsequently were used for the analysis of genetic diversity among the different populations. All of the commercial lines branched from a single node relative to the heritage varieties and the South Mexican turkey population. Heterozygosity of all individuals from the different turkey populations ranged from 0.17-2.73 SNPs/Kb, while heterozygosity of populations ranged from 0.73-1.64 SNPs/Kb. The average frequency of heterozygous SNPs in individual turkeys was 1.07 SNPs/Kb. Five genomic regions with very low nucleotide variation were identified in domestic turkeys that showed state of fixation towards alleles different than wild alleles.

Conclusion

The turkey genome is much less diverse with a relatively low frequency of heterozygous SNPs as compared to other livestock species like chicken and pig. The whole genome SNP discovery study in turkey resulted in the detection of 5.49 million putative SNPs compared to the reference genome. All commercial lines appear to share a common origin. Presence of different alleles/haplotypes in the SM population highlights that specific haplotypes have been selected in the modern domesticated turkey.

An encyclopedia of mouse DNA elements (Mouse ENCODE)

To complement the human Encyclopedia of DNA Elements (ENCODE) project and to enable a broad range of mouse genomics efforts, the Mouse ENCODE Consortium is applying the same experimental pipelines developed for human ENCODE to annotate the mouse genome.

Aluminum tolerance association mapping in triticale

Background

Crop production practices and industrialization processes result in increasing acidification of arable soils. At lower pH levels (below 5.0), aluminum (Al) remains in a cationic form that is toxic to plants, reducing growth and yield. The effect of aluminum on agronomic performance is particularly important in cereals like wheat, which has promoted the development of programs directed towards selection of tolerant forms. Even in intermediately tolerant cereals (i.e., triticale), the decrease in yield may be significant. In triticale, Al tolerance seems to be influenced by both wheat and rye genomes. However, little is known about the precise chromosomal location of tolerance-related genes, and whether wheat or rye genomes are crucial for the expression of that trait in the hybrid.

Results

A mapping population consisting of 232 advanced breeding triticale forms was developed and phenotyped for Al tolerance using physiological tests. AFLP, SSR and DArT marker platforms were applied to obtain a sufficiently large set of molecular markers (over 3000). Associations between the markers and the trait were tested using General (GLM) and Multiple (MLM) Linear Models, as well as the Statistical Machine Learning (SML) approach. The chromosomal locations of candidate markers were verified based on known assignments of SSRs and DArTs or by using genetic maps of rye and triticale.

Two candidate markers on chromosome 3R and 9, 15 and 11 on chromosomes 4R, 6R and 7R, respectively, were identified. The r² values were between 0.066 and 0.220 in most cases, indicating a good fit of the data, with better results obtained with the GML than the MLM approach. Several QTLs on rye chromosomes appeared to be involved in the phenotypic expression of the trait, suggesting that rye genome factors are predominantly responsible for Al tolerance in triticale.

Conclusions

The Diversity Arrays Technology was applied successfully to association mapping studies performed on triticale breeding forms. Statistical approaches allowed the identification of numerous markers associated with Al tolerance. Available rye and triticale genetic maps suggested the putative location of the markers and demonstrated that they formed several linked groups assigned to distinct chromosomes (3R, 4R, 6R and 7R). Markers associated with genomic regions under positive selection were identified and indirectly mapped in the vicinity of the Al-tolerant markers. The present findings were in agreement with prior reports.

Linking fold, function and phylogeny: a comparative genomics view on protein (domain) evolution

Domains are the building blocks of all globular proteins and present one of the most useful levels at which protein function can be understood. Through recombination and duplication of a limited set of domains, proteomes evolved and the collection of protein superfamilies in an organism formed. As such, the presence of a shared domain can be regarded as an indicator of similar function and evolutionary history, but it does not necessarily imply it since convergent evolution may give rise to similar gene functions as well as architectures.Through the wealth of sequences and annotation data brought about by genomics, evolutionary links can be sought for via homology relationships and comparative genomics, structural modeling and phylogenetics. The goal hereby is not only to predict the function of newly discovered proteins, but also to spell out their pathway of evolution and, possibly, identify their most likely origin. This can ultimately help to understand protein function and functional relationships of protein families. Additionally, through comparison with transcriptional data, evolutionary data can be linked to gene (and genome) activity and thus allow for the identification of common principles behind fast evolving proteins and relatively stable ones.In this review, we describe the basic principles of studying protein (domain) evolution and illustrate recent developments in molecular evolution and give valuable new insights in the field of comparative genomics. As an example, we include here molecular models of the multiple PDZ domain protein MUPP-1 and present a simple comparative genomic view on its structural course of evolution.

Evolutionary plasticity determination by orthologous groups distribution

Background

Genetic plasticity may be understood as the ability of a functional gene network to tolerate alterations in its components or structure. Usually, the studies involving gene modifications in the course of the evolution are concerned to nucleotide sequence alterations in closely related species. However, the analysis of large scale data about the distribution of gene families in non-exclusively closely related species can provide insights on how plastic or how conserved a given gene family is. Here, we analyze the abundance and diversity of all Eukaryotic Clusters of Orthologous Groups (KOG) present in STRING database, resulting in a total of 4,850 KOGs. This dataset comprises 481,421 proteins distributed among 55 eukaryotes.

Results

We propose an index to evaluate the evolutionary plasticity and conservation of an orthologous group based on its abundance and diversity across eukaryotes. To further KOG plasticity analysis, we estimate the evolutionary distance average among all proteins which take part in the same orthologous group. As a result, we found a strong correlation between the evolutionary distance average and the proposed evolutionary plasticity index. Additionally, we found low evolutionary plasticity in Saccharomyces cerevisiae genes associated with inviability and Mus musculus genes associated with early lethality. At last, we plot the evolutionary plasticity value in different gene networks from yeast and humans. As a result, it was possible to discriminate among higher and lower plastic areas of the gene networks analyzed.

Conclusions

The distribution of gene families brings valuable information on evolutionary plasticity which might be related with genetic plasticity. Accordingly, it is possible to discriminate among conserved and plastic orthologous groups by evaluating their abundance and diversity across eukaryotes.

Reviewers

This article was reviewed by Prof Manyuan Long, Hiroyuki Toh, and Sebastien Halary.

Reconstructing the history of large-scale genomic changes: biological questions and computational challenges

In addition to point mutations, larger-scale structural changes (including rearrangements, duplications, insertions, and deletions) are also prevalent between different mammalian genomes. Capturing these large-scale changes is critical to unraveling the history of mammalian evolution in order to better understand the human genome. It also has profound biomedical significance, because many human diseases are associated with structural genomic aberrations. The increasing number of mammalian genomes being sequenced as well as the recent advancement in DNA sequencing technologies are allowing us to identify these structural genomic changes with vastly greater accuracy. However, there are a considerable number of computational challenges related to these problems. In this article, we introduce the ancestral genome reconstruction problem, which enables us to explain the large-scale genomic changes between species in an evolutionary context. The application of these methods to within-species structural variation and disease genome analysis is also discussed. The target audience of this article is advanced undergraduate students in biology.

hzAnalyzer: detection, quantification, and visualization of contiguous homozygosity in high-density genotyping datasets

The analysis of contiguous homozygosity (runs of homozygous loci) in human genotyping datasets is critical in the search for causal disease variants in monogenic disorders, studies of population history and the identification of targets of natural selection. Here, we report methods for extracting homozygous segments from high-density genotyping datasets, quantifying their local genomic structure, identifying outstanding regions within the genome and visualizing results for comparative analysis between population samples.

Localization of pre- and postsynaptic cholinergic markers in rodent forebrain: a brief history and comparison of rat and mouse

Rat and mouse models are widely used for studies in cognition and pathophysiology, among others. Here, we sought to determine to what extent these two model species differ for cholinergic and cholinoceptive features. For this purpose, we focused on cholinergic innervation patterns based on choline acetyltransferase (ChAT) immunostaining, and the expression of muscarinic acetylcholine receptors (mAChRs) detected immunocytochemically. In this brief review we first place cholinergic and cholinoceptive markers in a historic perspective, and then provide an overview of recent publications on cholinergic studies and techniques to provide a literature survey of current research. Next, we compare mouse (C57Bl/J6) and rat (Wistar) cholinergic and cholinoceptive systems simultaneously stained, respectively, for ChAT (analyzed qualitatively) and mAChRs (analyzed qualitatively and quantitatively). In general, the topographic cholinergic innervation patterns of both rodent species are highly comparable, with only considerable (but region specific) differences in number of detectable cholinergic interneurons, which are more numerous in rat. In contrast, immunolabeling for mAChRs, detected by the monoclonal antibody M35, differs markedly in the forebrain between the two species. In mouse brain, basal levels of activated and/or internalized mAChRs (as a consequence of cholinergic neurotransmission) are significantly higher. This suggests a higher cholinergic tone in mouse than rat, and hence the animal model of choice may have consequences for cholinergic drug testing experiments.

Evolution of protein phosphorylation for distinct functional modules in vertebrate genomes

Recent publications have revealed that the evolution of phosphosites is influenced by the local protein structures and whether the phosphosites have characterized functions or not. With knowledge of the wide functional range of phosphorylation, we attempted to clarify whether the evolutionary conservation of phosphosites is different among distinct functional modules. We grouped the phosphosites in the human genome into the modules according to the functional categories of KEGG (Kyoto Encyclopedia of Genes and Genomes) and investigated their evolutionary conservation in vertebrate genomes from mouse to zebrafish. We have found that the phosphosites in the vertebrate-specific functional modules (VFMs), such as cellular signaling processes and responses to stimuli, are evolutionarily more conserved than those in the basic functional modules (BFMs), such as metabolic and genetic processes. The phosphosites in the VFMs are also significantly more conserved than their flanking regions, whereas those in the BFMs are not. These results hold for both serine/threonine and tyrosine residues, although the fraction of phosphorylated tyrosine residues is increased in the VFMs. Moreover, the difference in the evolutionary conservation of the phosphosites between the VFMs and BFMs could not be explained by the difference in the local protein structures. There is also a higher fraction of phosphosites with known functions in the VFMs than BFMs. Based on these findings, we have concluded that protein phosphorylation may play more dominant roles for the VFMs than BFMs during the vertebrate evolution. As phosphorylation is a quite rapid biological reaction, the VFMs that quickly respond to outer stimuli and inner signals might heavily depend on this regulatory mechanism. Our results imply that phosphorylation may have an essential role in the evolution of vertebrates.

Annotating non-coding regions of the genome

Most of the human genome consists of non-protein-coding DNA. Recently, progress has been made in annotating these non-coding regions through the interpretation of functional genomics experiments and comparative sequence analysis. One can conceptualize functional genomics analysis as involving a sequence of steps: turning the output of an experiment into a 'signal' at each base pair of the genome; smoothing this signal and segmenting it into small blocks of initial annotation; and then clustering these small blocks into larger derived annotations and networks. Finally, one can relate functional genomics annotations to conserved units and measures of conservation derived from comparative sequence analysis.

Functional analysis and identification of cis-regulatory elements of human chromosome 21 gene promoters

Given the inherent limitations of in silico studies relying solely on DNA sequence analysis, the functional characterization of mammalian promoters and associated cis-regulatory elements requires experimental support, which demands cloning and analysis of putative promoter regions. Focusing on human chromosome 21, we cloned 182 gene promoters of 2500 bp in length and conducted reporter gene assays on transfected-cell arrays. We found 56 promoters that were active in HEK293 cells, while another 49 promoters could be activated by treatment of cells with Trichostatin A or depletion of serum. We observed high correlations between promoter activities and endogenous transcript levels, RNA polymerase II occupancy, CpG islands and core promoter elements. Truncation of a subset of 62 promoters to ∼500 bp revealed that truncation rarely resulted in loss of activity, but rather in loss of responses to external stimuli, suggesting the presence of cis-regulatory response elements within distal promoter regions. In these regions, we found a strong enrichment of transcription factor binding sites that could potentially activate gene expression in the presence of stimuli. This study illustrates the modular functional architecture of chromosome 21 promoters and helps to reveal the complex mechanisms governing transcriptional regulation.

A time for atlases and atlases for time

Advances in neuroanatomy and computational power are leading to the construction of new digital brain atlases. Atlases are rising as indispensable tools for comparing anatomical data as well as being stimulators of new hypotheses and experimental designs. Brain atlases describe nervous systems which are inherently plastic and variable. Thus, the levels of brain plasticity and stereotypy would be important to evaluate as limiting factors in the context of static brain atlases. In this review, we discuss the extent of structural changes which neurons undergo over time, and how these changes would impact the static nature of atlases. We describe the anatomical stereotypy between neurons of the same type, highlighting the differences between invertebrates and vertebrates. We review some recent experimental advances in our understanding of anatomical dynamics in adult neural circuits, and how these are modulated by the organism's experience. In this respect, we discuss some analogies between brain atlases and the sequenced genome and the emerging epigenome. We argue that variability and plasticity of neurons are substantially high, and should thus be considered as integral features of high-resolution digital brain atlases.

[Reviews in comparative genomic research based on orthologs]

The orthologs have similar or even identical functions in different species, share the anological regulatory pathways, and play the close or even same role among species. Furthermore, the vast majority of the biological core functions were assumed to a considerable number of orthologous genes in organisms. Orthologs was the most reliable choices for functional annotation and analysis of genomic sequences, whose unique biological characteristics demonstrated that comparative genomics research based on orthologs will certainly provide clues for detecting the origin, expression, and loss of important functional genes during the biological evolution in different organisms. In this review, the fundamental characteristics of orthologous genes and the relationship between the orthologs and comparison genomics were recounted. The corresponding approaches and the current status in comparative genomic research based on the orthologs were summarized.

A Study of the Distribution of Phylogenetically Conserved Blocks within Clusters of Mammalian Homeobox Genes

Genome sequencing efforts of the last decade have produced a large amount of data, which has enabled whole-genome comparative analyses in order to locate potentially functional elements and study the overall patterns of phylogenetic conservation. In this paper we present a statistically based method for the characterization of these patterns in mammalian DNA sequences. We have applied this approach to the study of exceptionally well conserved homeobox gene clusters (Hox), based on the alignment of six species, and we have constructed a map of Hox cataloguing the conserved fragments, along with their locations in relation to the genes and other landmarks, sometimes showing unexpected layouts.

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.

Syntenator: multiple gene order alignments with a gene-specific scoring function

BACKGROUND

Identification of homologous regions or conserved syntenies across genomes is one crucial step in comparative genomics. This task is usually performed by genome alignment softwares like WABA or blastz. In case of conserved syntenies, such regions are defined as conserved gene orders. On the gene order level, homologous regions can even be found between distantly related genomes, which do not align on the nucleotide sequence level.

RESULTS

We present a novel approach to identify regions of conserved synteny across multiple genomes. Syntenator represents genomes and alignments thereof as partial order graphs (POGs). These POGs are aligned by a dynamic programming approach employing a gene-specific scoring function. The scoring function reflects the level of protein sequence similarity for each possible gene pair. Our method consistently defines larger homologous regions in pairwise gene order alignments than nucleotide-level comparisons. Our method is superior to methods that work on predefined homology gene sets (as implemented in Blockfinder). Syntenator successfully reproduces 80% of the EnsEMBL man-mouse conserved syntenic blocks. The full potential of our method becomes visible by comparing remotely related genomes and multiple genomes. Gene order alignments potentially resolve up to 75% of the EnsEMBL 1:many orthology relations and 27% of the many:many orthology relations.

CONCLUSION

We propose Syntenator as a software solution to reliably infer conserved syntenies among distantly related genomes. The software is available from http://www2.tuebingen.mpg.de/abt4/plone.

Transcriptional enhancement by GATA1-occupied DNA segments is strongly associated with evolutionary constraint on the binding site motif

Tissue development and function are exquisitely dependent on proper regulation of gene expression, but it remains controversial whether the genomic signals controlling this process are subject to strong selective constraint. While some studies show that highly constrained noncoding regions act to enhance transcription, other studies show that DNA segments with biochemical signatures of regulatory regions, such as occupancy by a transcription factor, are seemingly unconstrained across mammalian evolution. To test the possible correlation of selective constraint with enhancer activity, we used chromatin immunoprecipitation as an approach unbiased by either evolutionary constraint or prior knowledge of regulatory activity to identify DNA segments within a 66-Mb region of mouse chromosome 7 that are occupied by the erythroid transcription factor GATA1. DNA segments bound by GATA1 were identified by hybridization to high-density tiling arrays, validated by quantitative PCR, and tested for gene regulatory activity in erythroid cells. Whereas almost all of the occupied segments contain canonical WGATAR binding site motifs for GATA1, in only 45% of the cases is the motif deeply preserved (found at the orthologous position in placental mammals or more distant species). However, GATA1-bound segments with high enhancer activity tend to be the ones with an evolutionarily preserved WGATAR motif, and this relationship was confirmed by a loss-of-function assay. Thus, GATA1 binding sites that regulate gene expression during erythroid maturation are under strong selective constraint, while nonconstrained binding may have only a limited or indirect role in regulation.

Transcriptional regulatory elements in the human genome

The faithful execution of biological processes requires a precise and carefully orchestrated set of steps that depend on the proper spatial and temporal expression of genes. Here we review the various classes of transcriptional regulatory elements (core promoters, proximal promoters, distal enhancers, silencers, insulators/boundary elements, and locus control regions) and the molecular machinery (general transcription factors, activators, and coactivators) that interacts with the regulatory elements to mediate precisely controlled patterns of gene expression. The biological importance of transcriptional regulation is highlighted by examples of how alterations in these transcriptional components can lead to disease. Finally, we discuss the methods currently used to identify transcriptional regulatory elements, and the ability of these methods to be scaled up for the purpose of annotating the entire human genome.

Approaches to comparative sequence analysis: towards a functional view of vertebrate genomes

The comparison of genomic sequences is now a common approach to identifying and characterizing functional regions in vertebrate genomes. However, for theoretical reasons and because of practical issues, the generation of these data sets is non-trivial and can have many pitfalls. We are currently seeing an explosion of comparative sequence data, the benefits and limitations of which need to be disseminated to the scientific community. This Review provides a critical overview of the different types of sequence data that are available for analysis and of contemporary comparative sequence analysis methods, highlighting both their strengths and limitations. Approaches to determining the biological significance of constrained sequence are also explored.

MicroRNomics: a newly emerging approach for disease biology

Genomic evidence reveals that gene expression in humans is precisely controlled in cellular, tissue-type, temporal, and condition-specific manners. Completely understanding the regulatory mechanisms of gene expression is therefore one of the most important issues in genomic medicine. Surprisingly, recent analyses of the human and animal genomes have demonstrated that the majority of RNA transcripts are relatively small, noncoding RNAs (sncRNAs), rather than large, protein coding message RNAs (mRNAs). Moreover, these sncRNAs may represent a novel important layer of regulation for gene expression. The most important breakthrough in this new area is the discovery of microRNAs (miRNAs). miRNAs comprise a novel class of endogenous, small, noncoding RNAs that negatively regulate gene expression via degradation or translational inhibition of their target mRNAs. As a group, miRNAs may directly regulate approximately 30% of the genes in the human genome. In keeping with the nomenclature of RNomics, which is to study sncRNAs on the genomic scale, "microRNomics" is coined here to describe a novel subdiscipline of genomics that studies the identification, expression, biogenesis, structure, regulation of expression, targets, and biological functions of miRNAs on the genomic scale. A growing body of exciting evidence suggests that miRNAs are important regulators of cell differentiation, proliferation/growth, mobility, and apoptosis. These miRNAs therefore play important roles in development and physiology. Consequently, dysregulation of miRNA function may lead to human diseases such as cancer, cardiovascular disease, liver disease, immune dysfunction, and metabolic disorders. microRNomics may be a newly emerging approach for human disease biology.

Repetitive sequences in complex genomes: structure and evolution

Eukaryotic genomes contain vast amounts of repetitive DNA derived from transposable elements (TEs). Large-scale sequencing of these genomes has produced an unprecedented wealth of information about the origin, diversity, and genomic impact of what was once thought to be "junk DNA." This has also led to the identification of two new classes of DNA transposons, Helitrons and Polintons, as well as several new superfamilies and thousands of new families. TEs are evolutionary precursors of many genes, including RAG1, which plays a role in the vertebrate immune system. They are also the driving force in the evolution of epigenetic regulation and have a long-term impact on genomic stability and evolution. Remnants of TEs appear to be overrepresented in transcription regulatory modules and other regions conserved among distantly related species, which may have implications for our understanding of their impact on speciation.