Genetic consequences of programmed genome rearrangement

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

C-value paradox: Genesis in misconception that natural selection follows anthropocentric parameters of 'economy' and 'optimum'

C-value paradox refers to the lack of correlation between biological complexity and the intuitively expected protein-coding genomic information or DNA content. Here I discuss five questions about this paradox: i) Do biologically complex organisms carry more protein-coding genes? ii) Does variable accumulation of selfish/ junk/ parasitic DNA underlie the c-value paradox? iii) Can nucleoskeletal or nucleotypic function of DNA explain the enigma of orders of magnitude high levels of DNA in some 'lower' taxa or in taxonomically related species? iv) Can the newly understood noncoding but functional DNA explain the c-value paradox? and, v) Does natural selection uniformly apply the anthropocentric parameters for 'optimum' and 'economy'? Answers to Q.1-5 are largely negative. Biology presents numerous 'anomalous' examples where the same end function/ phenotype is attained in different organisms through astoundingly diverse ways that appear 'illogical' in our perceptions. Such evolutionary oddities exist because natural selection, unlike a designer, exploits random and stochastic events to modulate the existing system. Consequently, persistence of the new-found 'solution/s' often appear bizarre, uneconomic, and therefore, paradoxical to human logic. The unexpectedly high c-values in diverse organisms are irreversible evolutionary accidents that persisted, and the additional DNA often got repurposed over the evolutionary time scale. Therefore, the c-value paradox is a redundant issue. Future integrative biological studies should address evolutionary mechanisms and processes underlying sporadic DNA expansions/ contractions, and how the newly acquired DNA content has been repurposed in diverse groups.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

An improved germline genome assembly for the sea lamprey Petromyzon marinus illuminates the evolution of germline-specific chromosomes

Programmed DNA loss is a gene silencing mechanism that is employed by several vertebrate and nonvertebrate lineages, including all living jawless vertebrates and songbirds. Reconstructing the evolution of somatically eliminated (germline-specific) sequences in these species has proven challenging due to a high content of repeats and gene duplications in eliminated sequences and a corresponding lack of highly accurate and contiguous assemblies for these regions. Here, we present an improved assembly of the sea lamprey (Petromyzon marinus) genome that was generated using recently standardized methods that increase the contiguity and accuracy of vertebrate genome assemblies. This assembly resolves highly contiguous, somatically retained chromosomes and at least one germline-specific chromosome, permitting new analyses that reconstruct the timing, mode, and repercussions of recruitment of genes to the germline-specific fraction. These analyses reveal major roles of interchromosomal segmental duplication, intrachromosomal duplication, and positive selection for germline functions in the long-term evolution of germline-specific chromosomes.

Comparative genomics reveals insight into the evolutionary origin of massively scrambled genomes

Ciliates are microbial eukaryotes that undergo extensive programmed genome rearrangement, a natural genome editing process that converts long germline chromosomes into smaller gene-rich somatic chromosomes. Three well-studied ciliates include Oxytricha trifallax, Tetrahymena thermophila, and Paramecium tetraurelia, but only the Oxytricha lineage has a massively scrambled genome, whose assembly during development requires hundreds of thousands of precisely programmed DNA joining events, representing the most complex genome dynamics of any known organism. Here we study the emergence of such complex genomes by examining the origin and evolution of discontinuous and scrambled genes in the Oxytricha lineage. This study compares six genomes from three species, the germline and somatic genomes for Euplotes woodruffi, Tetmemena sp., and the model ciliate O. trifallax. We sequenced, assembled, and annotated the germline and somatic genomes of E. woodruffi, which provides an outgroup, and the germline genome of Tetmemena sp. We find that the germline genome of Tetmemena is as massively scrambled and interrupted as Oxytricha's: 13.6% of its gene loci require programmed translocations and/or inversions, with some genes requiring hundreds of precise gene editing events during development. This study revealed that the earlier diverged spirotrich, E. woodruffi, also has a scrambled genome, but only roughly half as many loci (7.3%) are scrambled. Furthermore, its scrambled genes are less complex, together supporting the position of Euplotes as a possible evolutionary intermediate in this lineage, in the process of accumulating complex evolutionary genome rearrangements, all of which require extensive repair to assemble functional coding regions. Comparative analysis also reveals that scrambled loci are often associated with local duplications, supporting a gradual model for the origin of complex, scrambled genomes via many small events of DNA duplication and decay.

Pervasive male-biased expression throughout the germline-specific regions of the sea lamprey genome supports key roles in sex differentiation and spermatogenesis

Sea lamprey undergo programmed genome rearrangement (PGR) in which ∼20% of the genome is jettisoned from somatic cells during embryogenesis. Although the role of PGR in embryonic development has been studied, the role of the germline-specific region (GSR) in gonad development is unknown. We analysed RNA-sequence data from 28 sea lamprey gonads sampled across life-history stages, generated a genome-guided de novo superTranscriptome with annotations, and identified germline-specific genes (GSGs). Overall, we identified 638 GSGs that are enriched for reproductive processes and exhibit 36x greater odds of being expressed in testes than ovaries. Next, while 55% of the GSGs have putative somatic paralogs, the somatic paralogs are not differentially expressed between sexes. Further, putative orthologs of some the male-biased GSGs have known functions in sex determination or differentiation in other vertebrates. We conclude that the GSR of sea lamprey plays an important role in testicular differentiation and potentially sex determination.

RNA-sequencing of sea lamprey gonads at different life-history stage identifies germline-specific genes which are highly expressed in males during spermatogenesis. This suggests a link between male-biased germline expression and sex differentiation in the sea lamprey.

Mendelian nightmares: the germline-restricted chromosome of songbirds

Germline-restricted chromosomes (GRCs) are accessory chromosomes that occur only in germ cells. They are eliminated from somatic cells through programmed DNA elimination during embryo development. GRCs have been observed in several unrelated animal taxa and show peculiar modes of non-Mendelian inheritance and within-individual elimination. Recent cytogenetic and phylogenomic evidence suggests that a GRC is present across the species-rich songbirds, but absent in non-passerine birds, implying that over half of all 10,500 bird species have extensive germline/soma genome differences. Here, we review recent insights gained from genomic, transcriptomic, and cytogenetic approaches with regard to the genetic content, phylogenetic distribution, and inheritance of the songbird GRC. While many questions remain unsolved in terms of GRC inheritance, elimination, and function, we discuss plausible scenarios and future directions for understanding this widespread form of programmed DNA elimination.

Gene-rich germline-restricted chromosomes in black-winged fungus gnats evolved through hybridization

Germline-restricted DNA has evolved in diverse animal taxa and is found in several vertebrate clades, nematodes, and flies. In these lineages, either portions of chromosomes or entire chromosomes are eliminated from somatic cells early in development, restricting portions of the genome to the germline. Little is known about why germline-restricted DNA has evolved, especially in flies, in which 3 diverse families, Chironomidae, Cecidomyiidae, and Sciaridae, carry germline-restricted chromosomes (GRCs). We conducted a genomic analysis of GRCs in the fungus gnat Bradysia (Sciara) coprophila (Diptera: Sciaridae), which has 2 large germline-restricted "L" chromosomes. We sequenced and assembled the genome of B. coprophila and used differences in sequence coverage and k-mer frequency between somatic and germline tissues to identify GRC sequence and compare it to the other chromosomes in the genome. We found that the GRCs in B. coprophila are large, gene rich, and have many genes with divergent homologs on other chromosomes in the genome. We also found that 2 divergent GRCs exist in the population we sequenced. GRC genes are more similar in sequence to genes from another Dipteran family (Cecidomyiidae) than to homologous genes from Sciaridae. This unexpected finding suggests that these chromosomes likely arose in Sciaridae through hybridization with a related lineage. These results provide a foundation from which to answer many questions about the evolution of GRCs in Sciaridae, such as how this hybridization event resulted in GRCs and what features on these chromosomes cause them to be restricted to the germline.

Delete and survive: strategies of programmed genetic material elimination in eukaryotes

Genome stability is a crucial feature of eukaryotic organisms because its alteration drastically affects the normal development and survival of cells and the organism as a whole. Nevertheless, some organisms can selectively eliminate part of their genomes from certain cell types during specific stages of ontogenesis. This review aims to describe the phenomenon of programmed DNA elimination, which includes chromatin diminution (together with programmed genome rearrangement or DNA rearrangements), B and sex chromosome elimination, paternal genome elimination, parasitically induced genome elimination, and genome elimination in animal and plant hybrids. During programmed DNA elimination, individual chromosomal fragments, whole chromosomes, and even entire parental genomes can be selectively removed. Programmed DNA elimination occurs independently in different organisms, ranging from ciliate protozoa to mammals. Depending on the sequences destined for exclusion, programmed DNA elimination may serve as a radical mechanism of dosage compensation and inactivation of unnecessary or dangerous genetic entities. In hybrids, genome elimination results from competition between parental genomes. Despite the different consequences of DNA elimination, all genetic material destined for elimination must be first recognised, epigenetically marked, separated, and then removed and degraded.

Developmental timing of programmed DNA elimination in Paramecium tetraurelia recapitulates germline transposon evolutionary dynamics

With its nuclear dualism, the ciliate Paramecium constitutes a unique model to study how host genomes cope with transposable elements (TEs). P. tetraurelia harbors two germline micronuclei (MICs) and a polyploid somatic macronucleus (MAC) that develops from one MIC at each sexual cycle. Throughout evolution, the MIC genome has been continuously colonized by TEs and related sequences that are removed from the somatic genome during MAC development. Whereas TE elimination is generally imprecise, excision of approximately 45,000 TE-derived internal eliminated sequences (IESs) is precise, allowing for functional gene assembly. Programmed DNA elimination is concomitant with genome amplification. It is guided by noncoding RNAs and repressive chromatin marks. A subset of IESs is excised independently of this epigenetic control, raising the question of how IESs are targeted for elimination. To gain insight into the determinants of IES excision, we established the developmental timing of DNA elimination genome-wide by combining fluorescence-assisted nuclear sorting with high-throughput sequencing. Essentially all IESs are excised within only one endoreplication round (32C to 64C), whereas TEs are eliminated at a later stage. We show that DNA elimination proceeds independently of replication. We defined four IES classes according to excision timing. The earliest excised IESs tend to be independent of epigenetic factors, display strong sequence signals at their ends, and originate from the most ancient integration events. We conclude that old IESs have been optimized during evolution for early and accurate excision by acquiring stronger sequence determinants and escaping epigenetic control.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

New Perspectives on the Evolution of Within-Individual Genome Variation and Germline/Soma Distinction

Genomes can vary significantly even within the same individual. The underlying mechanisms are manifold, ranging from somatic mutation and recombination, development-associated ploidy changes and genetic bottlenecks, over to programmed DNA elimination during germline/soma differentiation. In this perspective piece, we briefly review recent developments in the study of within-individual genome variation in eukaryotes and prokaryotes. We highlight a Society for Molecular Biology and Evolution 2020 virtual symposium entitled "Within-individual genome variation and germline/soma distinction" and the present Special Section of the same name in Genome Biology and Evolution, together fostering cross-taxon synergies in the field to identify and tackle key open questions in the understanding of within-individual genome variation.

Transcribed germline-limited coding sequences in Oxytricha trifallax

The germline-soma divide is a fundamental distinction in developmental biology, and different genes are expressed in germline and somatic cells throughout metazoan life cycles. Ciliates, a group of microbial eukaryotes, exhibit germline-somatic nuclear dimorphism within a single cell with two different genomes. The ciliate Oxytricha trifallax undergoes massive RNA-guided DNA elimination and genome rearrangement to produce a new somatic macronucleus (MAC) from a copy of the germline micronucleus (MIC). This process eliminates noncoding DNA sequences that interrupt genes and also deletes hundreds of germline-limited open reading frames (ORFs) that are transcribed during genome rearrangement. Here, we update the set of transcribed germline-limited ORFs (TGLOs) in O. trifallax. We show that TGLOs tend to be expressed during nuclear development and then are absent from the somatic MAC. We also demonstrate that exposure to synthetic RNA can reprogram TGLO retention in the somatic MAC and that TGLO retention leads to transcription outside the normal developmental program. These data suggest that TGLOs represent a group of developmentally regulated protein-coding sequences whose gene expression is terminated by DNA elimination.

Chromosome-level genome assembly of Lethenteron reissneri provides insights into lamprey evolution

The reissner lamprey Lethenteron reissneri, belonging to the class Cyclostomata, serves as a bridge between invertebrates and jawed vertebrates, and is considered the sister group of jawed vertebrates. However, despite this evolutionary significance, the genetic mechanisms underlying the adaptive evolution of the lamprey lineage remain unclear. Here, we assembled a 1.06 Gb chromosome-level draft genome of L. reissneri, with 72 chromosomes (ranging in length from 4.5 Mb to 25.9 Mb) and a scaffold N50 length of 13.23 Mb. Genome quality comparisons revealed that the reissner lamprey genome has higher completeness and contiguity than the previously published sea lamprey and Japanese lamprey genomes. Moreover, reissner lamprey, sea lamprey, and Japanese lamprey species share similar transposable element profiles and Hox gene cluster compositions, suggesting that a burst of transposable element activity and whole genome duplication occurred before their divergence. Additionally, the Lip gene copy numbers, which have been studied for their functions in the host defence system, were found to be expanded uniquely in lamprey lineages, suggesting key roles for these genes in lamprey evolution and adaptation. We also identified two neural-related genes, Nrn1 and Unc13a, with copy number expansions in jawed vertebrates, which may be functionally relevant to the origin of lamprey brains. Hence, this study not only provides the first chromosome-level reference genome for Cyclostomata, but also highlights features of the unique biology and adaptive evolution of the lamprey lineage.

Programmed DNA Elimination in Vertebrates

Over the last few decades, an increasing number of vertebrate taxa have been identified that undergo programmed genome rearrangement, or programmed DNA loss, during development. In these organisms, the genome of germ cells is often reproducibly different from the genome of all other cells within the body. Although we clearly have not identified all vertebrate taxa that undergo programmed genome loss, the list of species known to undergo loss now represents ∼10% of vertebrate species, including several basally diverging lineages. Recent studies have shed new light on the targets and mechanisms of DNA loss and their association with canonical modes of DNA silencing. Ultimately, expansion of these studies into a larger collection of taxa will aid in reconstructing patterns of shared/independent ancestry of programmed DNA loss in the vertebrate lineage, as well as more recent evolutionary events that have shaped the structure and content of eliminated DNA.

Comprehensive Chromosome End Remodeling during Programmed DNA Elimination

Germline and somatic genomes are in general the same in a multicellular organism. However, programmed DNA elimination leads to a reduced somatic genome compared to germline cells. Previous work on the parasitic nematode Ascaris demonstrated that programmed DNA elimination encompasses high-fidelity chromosomal breaks and loss of specific genome sequences including a major tandem repeat of 120 bp and ~1,000 germline-expressed genes. However, the precise chromosomal locations of these repeats, breaks regions, and eliminated genes remained unknown. We used PacBio long-read sequencing and chromosome conformation capture (Hi-C) to obtain fully assembled chromosomes of Ascaris germline and somatic genomes, enabling a complete chromosomal view of DNA elimination. We found that all 24 germline chromosomes undergo comprehensive chromosome end remodeling with DNA breaks in their subtelomeric regions and loss of distal sequences including the telomeres at both chromosome ends. All new Ascaris somatic chromosome ends are recapped by de novo telomere healing. We provide an ultrastructural analysis of Ascaris DNA elimination and show that eliminated DNA is incorporated into double membrane-bound structures, similar to micronuclei, during telophase of a DNA elimination mitosis. These micronuclei undergo dynamic changes including loss of active histone marks and localize to the cytoplasm following daughter nuclei formation and cytokinesis where they form autophagosomes. Comparative analysis of nematode chromosomes suggests that chromosome fusions occurred, forming Ascaris sex chromosomes that become independent chromosomes following DNA elimination breaks in somatic cells. These studies provide the first chromosomal view and define novel features and functions of metazoan programmed DNA elimination.

Graph based analysis for gene segment organization In a scrambled genome

DNA recombinant processes can involve gene segments that overlap or interleave with gene segments of another gene. Such gene segment appearances relative to each other are called here gene segment organization. We use graphs to represent the gene segment organization in a chromosome locus. Vertices of the graph represent contigs resulting after the recombination and the edges represent the gene segment organization prior to rearrangement. To each graph we associate a vector whose entries correspond to graph properties, and consider this vector as a point in a higher dimensional Euclidean space such that cluster formations and analysis can be performed with a hierarchical clustering method. The analysis is applied to a recently sequenced model organism Oxytricha trifallax, a species of ciliate with highly scrambled genome that undergoes massive rearrangement process after conjugation. The analysis shows some emerging star-like graph structures indicating that segments of a single gene can interleave, or even contain all of the segments from fifteen or more other genes in between its segments. We also observe that as many as six genes can have their segments mutually interleaving or overlapping.

ADFinder: accurate detection of programmed DNA elimination using NGS high-throughput sequencing data

Motivation

Programmed DNA elimination (PDE) plays a crucial role in the transitions between germline and somatic genomes in diverse organisms ranging from unicellular ciliates to multicellular nematodes. However, software specific for the detection of DNA splicing events is scarce. In this paper, we describe Accurate Deletion Finder (ADFinder), an efficient detector of PDEs using high-throughput sequencing data. ADFinder can predict PDEs with relatively low sequencing coverage, detect multiple alternative splicing forms in the same genomic location and calculate the frequency for each splicing event. This software will facilitate research of PDEs and all down-stream analyses.

Results

By analyzing genome-wide DNA splicing events in two micronuclear genomes of Oxytricha trifallax and Tetrahymena thermophila, we prove that ADFinder is effective in predicting large scale PDEs.

Availability and implementation

The source codes and manual of ADFinder are available in our GitHub website: https://github.com/weibozheng/ADFinder.

Supplementary information

Supplementary data are available at Bioinformatics online.

Programmed DNA elimination of germline development genes in songbirds

In some eukaryotes, germline and somatic genomes differ dramatically in their composition. Here we characterise a major germline–soma dissimilarity caused by a germline-restricted chromosome (GRC) in songbirds. We show that the zebra finch GRC contains >115 genes paralogous to single-copy genes on 18 autosomes and the Z chromosome, and is enriched in genes involved in female gonad development. Many genes are likely functional, evidenced by expression in testes and ovaries at the RNA and protein level. Using comparative genomics, we show that genes have been added to the GRC over millions of years of evolution, with embryonic development genes bicc1 and trim71 dating to the ancestor of songbirds and dozens of other genes added very recently. The somatic elimination of this evolutionarily dynamic chromosome in songbirds implies a unique mechanism to minimise genetic conflict between germline and soma, relevant to antagonistic pleiotropy, an evolutionary process underlying ageing and sexual traits.

Songbirds have extensive germline–soma genome differences due to developmental elimination of a germline-specific chromosome (GRC). Here, the authors show that the GRC contains dozens of expressed developmental genes, some of which have been on the GRC since the ancestor of all songbirds.

Germline-Specific Repetitive Elements in Programmatically Eliminated Chromosomes of the Sea Lamprey (Petromyzon marinus)

The sea lamprey (Petromyzon marinus) is one of few vertebrate species known to reproducibly eliminate large fractions of its genome during normal embryonic development. This germline-specific DNA is lost in the form of large fragments, including entire chromosomes, and available evidence suggests that DNA elimination acts as a permanent silencing mechanism that prevents the somatic expression of a specific subset of "germline" genes. However, reconstruction of eliminated regions has proven to be challenging due to the complexity of the lamprey karyotype. We applied an integrative approach aimed at further characterization of the large-scale structure of eliminated segments, including: (1) in silico identification of germline-enriched repeats; (2) mapping the chromosomal location of specific repetitive sequences in germline metaphases; and (3) 3D DNA/DNA-hybridization to embryonic lagging anaphases, which permitted us to both verify the specificity of elements to physically eliminated chromosomes and characterize the subcellular organization of these elements during elimination. This approach resulted in the discovery of several repetitive elements that are found exclusively on the eliminated chromosomes, which subsequently permitted the identification of 12 individual chromosomes that are programmatically eliminated during early embryogenesis. The fidelity and specificity of these highly abundant sequences, their distinctive patterning in eliminated chromosomes, and subcellular localization in elimination anaphases suggest that these sequences might contribute to the specific targeting of chromosomes for elimination or possibly in molecular interactions that mediate their decelerated poleward movement in chromosome elimination anaphases, isolation into micronuclei and eventual degradation.

Programmed DNA Elimination: Keeping Germline Genes in Their Place

Each of our cells contains a full set of instructions needed to make an entire human: the genome. But a few special species buck this trend. A new study now identifies the first germline-specific gene in zebra finch, one of a small number of vertebrates that are known to undergo developmentally programmed DNA elimination.

Characterization and Evolution of Germ1, an Element that Undergoes Diminution in Lampreys (Cyclostomata: Petromyzontidae)

The sea lamprey (Petromyzon marinus) undergoes substantial genomic alterations during embryogenesis in which specific sequences are deleted from the genome of somatic cells yet retained in cells of the germ line. One element that undergoes diminution in P. marinus is Germ1, which consists of a somatically rare (SR) region and a fragment of 28S rDNA. Although the SR-region has been used as a marker for genomic alterations in lampreys, the evolutionary significance of its diminution is unknown. We examined the Germ1 element in five additional species of lamprey to better understand its evolutionary significance. Each representative species contained sequences similar enough to the Germ1 element of P. marinus to be detected via PCR and Southern hybridizations, although the SR-regions of Lampetra aepyptera and Lethenteron appendix are quite divergent from the homologous sequences of Petromyzon and three species of Ichthyomyzon. Lamprey Germ1 sequences have a number of features characteristic of the R2 retrotransposon, a mobile element that specifically targets 28S rDNA. Phylogenetic analyses of the SR-regions revealed patterns generally consistent with relationships among the species included in our study, although the 28S-fragments of each species/genus were most closely related to its own functional rDNA, suggesting that the two components of Germ1 were assembled independently in each lineage. Southern hybridizations showed evidence of genomic alterations involving Germ1 in each species. Our results suggest that Germ1 is a R2 retroelement that occurs in the genome of P. marinus and other petromyzontid lampreys, and that its diminution is incidental to the reduction in rDNA copies during embryogenesis.

Hagfish and lamprey Hox genes reveal conservation of temporal colinearity in vertebrates

Hox genes exert fundamental roles for proper regional specification along the main rostro-caudal axis of animal embryos. They are generally expressed in restricted spatial domains according to their position in the cluster (spatial colinearity)-a feature that is conserved across bilaterians. In jawed vertebrates (gnathostomes), the position in the cluster also determines the onset of expression of Hox genes (a feature known as whole-cluster temporal colinearity (WTC)), while in invertebrates this phenomenon is displayed as a subcluster-level temporal colinearity. However, little is known about the expression profile of Hox genes in jawless vertebrates (cyclostomes); therefore, the evolutionary origin of WTC, as seen in gnathostomes, remains a mystery. Here, we show that Hox genes in cyclostomes are expressed according to WTC during development. We investigated the Hox repertoire and Hox gene expression profiles in three different species-a hagfish, a lamprey and a shark-encompassing the two major groups of vertebrates, and found that these are expressed following a whole-cluster, temporally staggered pattern, indicating that WTC has been conserved during the past 500 million years despite drastically different genome evolution and morphological outputs between jawless and jawed vertebrates.

The sea lamprey germline genome provides insights into programmed genome rearrangement and vertebrate evolution

The sea lamprey (Petromyzon marinus) serves as a comparative model for reconstructing vertebrate evolution. To enable more informed analyses, we developed a new assembly of the lamprey germline genome that integrates several complementary datasets. Analysis of this highly contiguous (chromosome-scale) assembly reveals that both chromosomal and whole-genome duplications have played significant roles in the evolution of ancestral vertebrate and lamprey genomes, including chromosomes that carry the six lamprey HOX clusters. The assembly also contains several hundred genes that are reproducibly eliminated from somatic cells during early development in lamprey. Comparative analyses show that gnathostome (mouse) homologs of these genes are frequently marked by Polycomb Repressive Complexes (PRCs) in embryonic stem cells, suggesting overlaps in the regulatory logic of somatic DNA elimination and repressive/bivalent states that are regulated by early embryonic PRCs. This new assembly will enhance diverse studies that are informed by lampreys’ unique biology and evolutionary/comparative perspective.

Comparative genome analysis of programmed DNA elimination in nematodes

Programmed DNA elimination is a developmentally regulated process leading to the reproducible loss of specific genomic sequences. DNA elimination occurs in unicellular ciliates and a variety of metazoans, including invertebrates and vertebrates. In metazoa, DNA elimination typically occurs in somatic cells during early development, leaving the germline genome intact. Reference genomes for metazoa that undergo DNA elimination are not available. Here, we generated germline and somatic reference genome sequences of the DNA eliminating pig parasitic nematode Ascaris suum and the horse parasite Parascaris univalens. In addition, we carried out in-depth analyses of DNA elimination in the parasitic nematode of humans, Ascaris lumbricoides, and the parasitic nematode of dogs, Toxocara canis. Our analysis of nematode DNA elimination reveals that in all species, repetitive sequences (that differ among the genera) and germline-expressed genes (approximately 1000–2000 or 5%–10% of the genes) are eliminated. Thirty-five percent of these eliminated genes are conserved among these nematodes, defining a core set of eliminated genes that are preferentially expressed during spermatogenesis. Our analysis supports the view that DNA elimination in nematodes silences germline-expressed genes. Over half of the chromosome break sites are conserved between Ascaris and Parascaris, whereas only 10% are conserved in the more divergent T. canis. Analysis of the chromosomal breakage regions suggests a sequence-independent mechanism for DNA breakage followed by telomere healing, with the formation of more accessible chromatin in the break regions prior to DNA elimination. Our genome assemblies and annotations also provide comprehensive resources for analysis of DNA elimination, parasitology research, and comparative nematode genome and epigenome studies.

A germline-limited piggyBac transposase gene is required for precise excision in Tetrahymena genome rearrangement

Developmentally programmed genome rearrangement accompanies differentiation of the silent germline micronucleus into the transcriptionally active somatic macronucleus in the ciliated protozoan Tetrahymena thermophila. Internal eliminated sequences (IES) are excised, followed by rejoining of MAC-destined sequences, while fragmentation occurs at conserved chromosome breakage sequences, generating macronuclear chromosomes. Some macronuclear chromosomes, referred to as non-maintained chromosomes (NMC), are lost soon after differentiation. Large NMC contain genes implicated in development-specific roles. One such gene encodes the domesticated piggyBac transposase TPB6, required for heterochromatin-dependent precise excision of IES residing within exons of functionally important genes. These conserved exonic IES determine alternative transcription products in the developing macronucleus; some even contain free-standing genes. Examples of precise loss of some exonic IES in the micronucleus and retention of others in the macronucleus of related species suggest an evolutionary analogy to introns. Our results reveal that germline-limited sequences can encode genes with specific expression patterns and development-related functions, which may be a recurring theme in eukaryotic organisms experiencing programmed genome rearrangement during germline to soma differentiation.

Deep ancestry of programmed genome rearrangement in lampreys

In most multicellular organisms, the structure and content of the genome is rigorously maintained over the course of development. However some species have evolved genome biologies that permit, or require, developmentally regulated changes in the physical structure and content of the genome (programmed genome rearrangement: PGR). Relatively few vertebrates are known to undergo PGR, although all agnathans surveyed to date (several hagfish and one lamprey: Petromyzon marinus) show evidence of large scale PGR. To further resolve the ancestry of PGR within vertebrates, we developed probes that allow simultaneous tracking of nearly all sequences eliminated by PGR in P. marinus and a second lamprey species (Entosphenus tridentatus). These comparative analyses reveal conserved subcellular structures (lagging chromatin and micronuclei) associated with PGR and provide the first comparative embryological evidence in support of the idea that PGR represents an ancient and evolutionarily stable strategy for regulating inherent developmental/genetic conflicts between germline and soma.

Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

A second estrogen receptor from Japanese lamprey (Lethenteron japonicum) does not have activities for estrogen binding and transcription

Estrogens regulate many physiological responses in vertebrates by binding to the estrogen receptor (ER), a ligand-activated transcription factor. To understand the evolution of vertebrate ERs and to investigate how estrogen acts in a jawless vertebrate, we used degenerate primer sets and PCR to isolate DNA fragments encoding two distinct ER subtypes, Esr1a and Esr1b from the Japanese lamprey, Lethenteron japonicum. Phylogenetic analysis indicates that these two ERs are the result of lineage-specific gene duplication within the jawless fishes, different from the previous duplication event of Esr1 (ERα) and Esr2 (ERβ) within the jawed vertebrates. Reporter gene assays show that lamprey Esr1a displays both constitutive and estrogen-dependent activation of gene transcription. Domain swapping experiments indicate that constitutive activity resides in the A/B domain of lamprey Esr1a. Unexpectedly, lamprey Esr1b does not bind estradiol and is not stimulated by other estrogens, androgens or corticosteroids. A 3D model of lamprey Esr1b suggests that although estradiol fits into the steroid binding site, some stabilizing contacts between the ligand and side chains that are found in human Esr1 and Esr2 are missing in lamprey Esr1b.

The Nuclear DNA Content and Genetic Diversity of Lampetra morii

We investigated the nuclear DNA content and genetic diversity of a river lamprey, the Korean lamprey Lampetra morii, which is distributed in the northeast of China. L. morii spends its whole life cycle in fresh water, and its adult size is relatively small (~160 mm long) compared with that of other lampreys. The haploid nuclear DNA content of L. morii is 1.618 pg (approximately 1.582 Gb) in germline cells, and there is ~15% germline DNA loss in somatic cells. These values are significantly smaller than those of Petromyzon marinus, a lamprey with a published draft genome. The chromosomes of L. morii are small and acrocentric, with a diploid modal number of 2n = 132, lower than some other lampreys. Sequence and AFLP analyses suggest that the allelic polymorphism rate (~0.14% based on examined nuclear and mitochondrial DNA sequences) of L. morii is much lower than that (~2%) of P. marinus. Phylogenetic analysis based on a mitochondrial DNA fragment confirms that L. morii belongs to the genus Lampetra, which, together with the genus Lethenteron, forms a sister group to P. marinus. These genetic background data are valuable for subsequent genetic and genomic research on L. morii.

The Enigma of Progressively Partial Endoreplication: New Insights Provided by Flow Cytometry and Next-Generation Sequencing

In many plant species, somatic cell differentiation is accompanied by endoreduplication, a process during which cells undergo one or more rounds of DNA replication cycles in the absence of mitosis, resulting in nuclei with multiples of 2C DNA amounts (4C, 8C, 16C, etc.). In some orchids, a disproportionate increase in nuclear DNA contents has been observed, where successive endoreduplication cycles result in DNA amounts 2C + P, 2C + 3P, 2C + 7P, etc., where P is the DNA content of the replicated part of the 2C nuclear genome. This unique phenomenon was termed "progressively partial endoreplication" (PPE). We investigated processes behind the PPE in Ludisia discolor using flow cytometry (FCM) and Illumina sequencing. In particular, we wanted to determine whether chromatin elimination or incomplete genome duplication was involved, and to identify types of DNA sequences that were affected. Cell cycle analysis of root tip cell nuclei pulse-labeled with EdU revealed two cell cycles, one ending above the population of nuclei with 2C + P content, and the other with a typical "horseshoe" pattern of S-phase nuclei ranging from 2C to 4C DNA contents. The process leading to nuclei with 2C + P amounts therefore involves incomplete genome replication. Subsequent Illumina sequencing of flow-sorted 2C and 2C + P nuclei showed that all types of repetitive DNA sequences were affected during PPE; a complete elimination of any specific type of repetitive DNA was not observed. We hypothesize that PPE is part of a highly controlled transition mechanism from proliferation phase to differentiation phase of plant tissue development.

Cellular and Molecular Features of Developmentally Programmed Genome Rearrangement in a Vertebrate (Sea Lamprey: Petromyzon marinus)

The sea lamprey (Petromyzon marinus) represents one of the few vertebrate species known to undergo large-scale programmatic elimination of genomic DNA over the course of its normal development. Programmed genome rearrangements (PGRs) result in the reproducible loss of ~20% of the genome from somatic cell lineages during early embryogenesis. Studies of PGR hold the potential to provide novel insights related to the maintenance of genome stability during the cell cycle and coordination between mechanisms responsible for the accurate distribution of chromosomes into daughter cells, yet little is known regarding the mechanistic basis or cellular context of PGR in this or any other vertebrate lineage. Here we identify epigenetic silencing events that are associated with the programmed elimination of DNA and describe the spatiotemporal dynamics of PGR during lamprey embryogenesis. In situ analyses reveal that the earliest DNA methylation (and to some extent H3K9 trimethylation) events are limited to specific extranuclear structures (micronuclei) containing eliminated DNA. During early embryogenesis a majority of micronuclei (~60%) show strong enrichment for repressive chromatin modifications (H3K9me3 and 5meC). These analyses also led to the discovery that eliminated DNA is packaged into chromatin that does not migrate with somatically retained chromosomes during anaphase, a condition that is superficially similar to lagging chromosomes observed in some cancer subtypes. Closer examination of “lagging” chromatin revealed distributions of repetitive elements, cytoskeletal contacts and chromatin contacts that provide new insights into the cellular mechanisms underlying the programmed loss of these segments. Our analyses provide additional perspective on the cellular and molecular context of PGR, identify new structures associated with elimination of DNA and reveal that PGR is completed over the course of several successive cell divisions.

Lampreys possess a fascinating genome biology wherein large portions of the genome, including large numbers of genes, are programmatically deleted during development. The lamprey therefore represents a uniquely informative system with respect to several broad areas of biology, including genome stability/rearrangement, epigenetic silencing, and the establishment and maintenance of pluripotency. However, little is known regarding the cellular context or mechanism of deletion, partly due to the challenges of observing rearrangements in situ. Here we present analyses and new techniques that significantly advance our understanding of the subcellular context of programmed rearrangements and interactions between programmed deletion and canonical DNA silencing mechanisms. These analyses demonstrate that DNA elimination occurs earlier in embryogenesis than was previously recognized and reveal several new cellular and molecular aspects of programmed DNA loss. Specifically we show that eliminated DNA exhibits a unique migration pattern during cell division, is packaged into discreet subcellular structures later in the cell cycle, and undergoes epigenetic silencing through DNA and histone methylation. These observations provide new insight into the mechanisms underlying programmed DNA loss and suggest a functional link between programmed DNA loss and other, more conserved gene silencing pathways.

Characterization of Somatically-Eliminated Genes During Development of the Sea Lamprey (Petromyzon marinus)

The sea lamprey (Petromyzon marinus) is a basal vertebrate that undergoes developmentally programmed genome rearrangements (PGRs) during early development. These events facilitate the elimination of ∼20% of the genome from the somatic cell lineage, resulting in distinct somatic and germline genomes. Thus far only a handful of germline-specific genes have been definitively identified within the estimated 500 Mb of DNA that is deleted during PGR, although a few thousand germline-specific genes are thought to exist. To improve our understanding of the evolutionary/developmental logic of PGR, we generated computational predictions to identify candidate germline-specific genes within a new transcriptomic dataset derived from adult germline and the early embryonic stages during which PGR occurs. Follow-up validation studies identified 44 germline-specific genes and further characterized patterns of transcription and DNA loss during early embryogenesis. Expression analyses reveal that many of these genes are differentially expressed during early embryogenesis and presumably function in the early development of the germline. Ontology analyses indicate that many of these germline-specific genes play known roles in germline development, pluripotency, and oncogenesis (when misexpressed). These studies provide support for the theory that PGR serves to segregate molecular functions related to germline development/pluripotency in order to prevent their potential misexpression in somatic cells. This larger set of eliminated genes also allows us to extend the evolutionary/developmental breadth of this theory, as some deleted genes (or their gnathostome homologs) appear to be associated with the early development of somatic lineages, perhaps through the evolution of novel functions within gnathostome lineages.

Programmed Rearrangement in Ciliates: Paramecium

Programmed genome rearrangements in the ciliate Paramecium provide a nice illustration of the impact of transposons on genome evolution and plasticity. During the sexual cycle, development of the somatic macronucleus involves elimination of ∼30% of the germline genome, including repeated DNA (e.g., transposons) and ∼45,000 single-copy internal eliminated sequences (IES). IES excision is a precise cut-and-close process, in which double-stranded DNA cleavage at IES ends depends on PiggyMac, a domesticated piggyBac transposase. Genome-wide analysis has revealed that at least a fraction of IESs originate from Tc/mariner transposons unrelated to piggyBac. Moreover, genomic sequences with no transposon origin, such as gene promoters, can be excised reproducibly as IESs, indicating that genome rearrangements contribute to the control of gene expression. How the system has evolved to allow elimination of DNA sequences with no recognizable conserved motif has been the subject of extensive research during the past two decades. Increasing evidence has accumulated for the participation of noncoding RNAs in epigenetic control of elimination for a subset of IESs, and in trans-generational inheritance of alternative rearrangement patterns. This chapter summarizes our current knowledge of the structure of the germline and somatic genomes for the model species Paramecium tetraurelia, and describes the DNA cleavage and repair factors that constitute the IES excision machinery. We present an overview of the role of specialized RNA interference machineries and their associated noncoding RNAs in the control of DNA elimination. Finally, we discuss how RNA-dependent modification and/or remodeling of chromatin may guide PiggyMac to its cognate cleavage sites.

Programmed Genome Rearrangements in the Ciliate Oxytricha

The ciliate Oxytricha is a microbial eukaryote with two genomes, one of which experiences extensive genome remodeling during development. Each round of conjugation initiates a cascade of events that construct a transcriptionally active somatic genome from a scrambled germline genome, with considerable help from both long and small noncoding RNAs. This process of genome remodeling entails massive DNA deletion and reshuffling of remaining DNA segments to form functional genes from their interrupted and scrambled germline precursors. The use of Oxytricha as a model system provides an opportunity to study an exaggerated form of programmed genome rearrangement. Furthermore, studying the mechanisms that maintain nuclear dimorphism and mediate genome rearrangement has demonstrated a surprising plasticity and diversity of noncoding RNA pathways, with new roles that go beyond conventional gene silencing. Another aspect of ciliate genetics is their unorthodox patterns of RNA-mediated, epigenetic inheritance that rival Mendelian inheritance. This review takes the reader through the key experiments in a model eukaryote that led to fundamental discoveries in RNA biology and pushes the biological limits of DNA processing.

The p53-Mdm2 interaction and the E3 ligase activity of Mdm2/Mdm4 are conserved from lampreys to humans

Here, Coffill et al. characterize Tp53, Tp63, and Tp73 in a jawless vertebrate, the Japanese lamprey, as well as the Mdm2 and Mdm4 genes using genome analysis. Functional analysis reveals conservation of p63 and p73 compared with p53, which shows substantial variability within the C-terminal and N-terminal domains, and that lamprey Mdm2 degrades human p53 with great efficiency; however, this interaction is not inhibited by currently available small molecule inhibitors of the human HDM2 protein.

The extant jawless vertebrates, represented by lampreys and hagfish, are the oldest group of vertebrates and provide an interesting genomic evolutionary pivot point between invertebrates and jawed vertebrates. Through genome analysis of one of these jawless vertebrates, the Japanese lamprey (Lethenteron japonicum), we identified all three members of the important p53 transcription factor family—Tp53, Tp63, and Tp73—as well as the Mdm2 and Mdm4 genes. These genes and their products are significant cellular regulators in human cancer, and further examination of their roles in this most distant vertebrate relative sheds light on their origin and coevolution. Their important role in response to DNA damage has been highlighted by the discovery of multiple copies of the Tp53 gene in elephants. Expression of lamprey p53, Mdm2, and Mdm4 proteins in mammalian cells reveals that the p53–Mdm2 interaction and the Mdm2/Mdm4 E3 ligase activity existed in the common ancestor of vertebrates and have been conserved for >500 million years of vertebrate evolution. Lamprey Mdm2 degrades human p53 with great efficiency, but this interaction is not blocked by currently available small molecule inhibitors of the human HDM2 protein, suggesting utility of lamprey Mdm2 in the study of the human p53 signaling pathway.

The transcriptional repressor Blimp-1 acts downstream of BMP signaling to generate primordial germ cells in the cricket Gryllus bimaculatus

Segregation of the germ line from the soma is an essential event for transmission of genetic information across generations in all sexually reproducing animals. Although some well-studied systems such as Drosophila and Xenopus use maternally inherited germ determinants to specify germ cells, most animals, including mice, appear to utilize zygotic inductive cell signals to specify germ cells during later embryogenesis. Such inductive germ cell specification is thought to be an ancestral trait of Bilateria, but major questions remain as to the nature of an ancestral mechanism to induce germ cells, and how that mechanism evolved. We previously reported that BMP signaling-based germ cell induction is conserved in both the mouse Mus musculus and the cricket Gryllus bimaculatus, which is an emerging model organism for functional studies of induction-based germ cell formation. In order to gain further insight into the functional evolution of germ cell specification, here we examined the Gryllus ortholog of the transcription factor Blimp-1 (also known as Prdm1), which is a widely conserved bilaterian gene known to play a crucial role in the specification of germ cells in mice. Our functional analyses of the Gryllus Blimp-1 ortholog revealed that it is essential for Gryllus primordial germ cell development, and is regulated by upstream input from the BMP signaling pathway. This functional conservation of the epistatic relationship between BMP signaling and Blimp-1 in inductive germ cell specification between mouse and cricket supports the hypothesis that this molecular mechanism regulated primordial germ cell specification in a last common bilaterian ancestor.

Lampreys as Diverse Model Organisms in the Genomics Era

Lampreys, one of the two surviving groups of ancient vertebrates, have become important models for study in diverse fields of biology. Lampreys (of which there are approximately 40 species) are being studied, for example, (a) to control pest sea lamprey in the North American Great Lakes and to restore declining populations of native species elsewhere; (b) in biomedical research, focusing particularly on the regenerative capability of lampreys; and (c) by developmental biologists studying the evolution of key vertebrate characters. Although a lack of genetic resources has hindered research on the mechanisms regulating many aspects of lamprey life history and development, formerly intractable questions are now amenable to investigation following the recent publication of the sea lamprey genome. Here, we provide an overview of the ways in which genomic tools are currently being deployed to tackle diverse research questions and suggest several areas that may benefit from the availability of the sea lamprey genome.

Evolution of the CNS myelin gene regulatory program

Myelin is a specialized subcellular structure that evolved uniquely in vertebrates. A myelinated axon conducts action potentials many times faster than an unmyelinated axon of the same diameter; for the same conduction speed, the unmyelinated axon would need a much larger diameter and volume than its myelinated counterpart. Hence myelin speeds information transfer and saves space, allowing the evolution of a powerful yet portable brain. Myelination in the central nervous system (CNS) is controlled by a gene regulatory program that features a number of master transcriptional regulators including Olig1, Olig2 and Myrf. Olig family genes evolved from a single ancestral gene in non-chordates. Olig2, which executes multiple functions with regard to oligodendrocyte identity and development in vertebrates, might have evolved functional versatility through post-translational modification, especially phosphorylation, as illustrated by its evolutionarily conserved serine/threonine phospho-acceptor sites and its accumulation of serine residues during more recent stages of vertebrate evolution. Olig1, derived from a duplicated copy of Olig2 in early bony fish, is involved in oligodendrocyte development and is critical to remyelination in bony vertebrates, but is lost in birds. The origin of Myrf orthologs might be the result of DNA integration between an invading phage or bacterium and an early protist, producing a fusion protein capable of self-cleavage and DNA binding. Myrf seems to have adopted new functions in early vertebrates - initiation of the CNS myelination program as well as the maintenance of mature oligodendrocyte identity and myelin structure - by developing new ways to interact with DNA motifs specific to myelin genes. This article is part of a Special Issue entitled SI: Myelin Evolution.

The evolution of animal genomes

Genome sequences are now available for hundreds of species sampled across the animal phylogeny, bringing key features of animal genome evolution into sharper focus. The field of animal evolutionary genomics has focused on identifying and classifying the diversity genomic features, reconstructing the history of evolutionary changes in animal genomes, and testing hypotheses about the evolutionary relationships of animals. The grand challenges moving forward are to connect evolutionary changes in genomes with particular evolutionary changes in phenotypes, and to determine which changes are driven by selection. This will require far greater genome sampling both across and within species, extensive phenotype data, a well resolved animal phylogeny, and advances in comparative methods.

The sea lamprey meiotic map improves resolution of ancient vertebrate genome duplications

It is generally accepted that many genes present in vertebrate genomes owe their origin to two whole-genome duplications that occurred deep in the ancestry of the vertebrate lineage. However, details regarding the timing and outcome of these duplications are not well resolved. We present high-density meiotic and comparative genomic maps for the sea lamprey (Petromyzon marinus), a representative of an ancient lineage that diverged from all other vertebrates ∼550 million years ago. Linkage analyses yielded a total of 95 linkage groups, similar to the estimated number of germline chromosomes (1n ∼ 99), spanning a total of 5570.25 cM. Comparative mapping data yield strong support for the hypothesis that a single whole-genome duplication occurred in the basal vertebrate lineage, but do not strongly support a hypothetical second event. Rather, these comparative maps reveal several evolutionarily independent segmental duplications occurring over the last 600+ million years of chordate evolution. This refined history of vertebrate genome duplication should permit more precise investigations of vertebrate evolution.

KLF/SP Transcription Factor Family Evolution: Expansion, Diversification, and Innovation in Eukaryotes

The Krüppel-like factor and specificity protein (KLF/SP) genes play key roles in critical biological processes including stem cell maintenance, cell proliferation, embryonic development, tissue differentiation, and metabolism and their dysregulation has been implicated in a number of human diseases and cancers. Although many KLF/SP genes have been characterized in a handful of bilaterian lineages, little is known about the KLF/SP gene family in nonbilaterians and virtually nothing is known outside the metazoans. Here, we analyze and discuss the origins and evolutionary history of the KLF/SP transcription factor family and associated transactivation/repression domains. We have identified and characterized the complete KLF/SP gene complement from the genomes of 48 species spanning the Eukarya. We have also examined the phylogenetic distribution of transactivation/repression domains associated with this gene family. We report that the origin of the KLF/SP gene family predates the divergence of the Metazoa. Furthermore, the expansion of the KLF/SP gene family is paralleled by diversification of transactivation domains via both acquisitions of pre-existing ancient domains as well as by the appearance of novel domains exclusive to this gene family and is strongly associated with the expansion of cell type complexity.

Expression patterns of Eph genes in the "dual visual development" of the lamprey and their significance in the evolution of vision in vertebrates

Image-forming vision is crucial to animals for recognizing objects in their environment. In vertebrates, this type of vision is achieved with paired camera eyes and topographic projection of the optic nerve. Topographic projection is established by an orthogonal gradient of axon guidance molecules, such as Ephs. To explore the evolution of image-forming vision in vertebrates, lampreys, which belong to the basal lineage of vertebrates, are key animals because they show unique "dual visual development." In the embryonic and pre-ammocoete larval stage (the "primary" phase), photoreceptive "ocellus-like" eyes develop, but there is no retinotectal optic nerve projection. In the late ammocoete larval stage (the "secondary" phase), the eyes grow and form into camera eyes, and retinotectal projection is newly formed. After metamorphosis, this retinotectal projection in adult lampreys is topographic, similar to that of gnathostomes. In this study, we explored the involvement of Ephs in lamprey "dual visual development" and establishment of the image-form vision. We found that gnathostome-like orthogonal gradient expression was present in the retina during the "secondary" phase; i.e., EphB showed a gradient of expression along the dorsoventral axis, while EphC was expressed along the anteroposterior axis. However, no orthogonal gradient expression was observed during the "primary" phase. These observations suggest that Ephs are likely recruited de novo for the guidance of topographical "second" optic nerve projection. Transformations during lamprey "dual visual development" may represent "recapitulation" from a protochordate-like ancestor to a gnathostome-like vertebrate ancestor.

Evolution of retinoic acid receptors in chordates: insights from three lamprey species, Lampetra fluviatilis, Petromyzon marinus, and Lethenteron japonicum

Background

Retinoic acid (RA) signaling controls many developmental processes in chordates, from early axis specification to late organogenesis. The functions of RA are chiefly mediated by a subfamily of nuclear hormone receptors, the retinoic acid receptors (RARs), that act as ligand-activated transcription factors. While RARs have been extensively studied in jawed vertebrates (that is, gnathostomes) and invertebrate chordates, very little is known about the repertoire and developmental roles of RARs in cyclostomes, which are extant jawless vertebrates. Here, we present the first extensive study of cyclostome RARs focusing on three different lamprey species: the European freshwater lamprey, Lampetra fluviatilis, the sea lamprey, Petromyzon marinus, and the Japanese lamprey, Lethenteron japonicum.

Results

We identified four rar paralogs (rar1, rar2, rar3, and rar4) in each of the three lamprey species, and phylogenetic analyses indicate a complex evolutionary history of lamprey rar genes including the origin of rar1 and rar4 by lineage-specific duplication after the lamprey-hagfish split. We further assessed their expression patterns during embryonic development by in situ hybridization. The results show that lamprey rar genes are generally characterized by dynamic and highly specific expression domains in different embryonic tissues. In particular, lamprey rar genes exhibit combinatorial expression domains in the anterior central nervous system (CNS) and the pharyngeal region.

Conclusions

Our results indicate that the genome of lampreys encodes at least four rar genes and suggest that the lamprey rar complement arose from vertebrate-specific whole genome duplications followed by a lamprey-specific duplication event. Moreover, we describe a combinatorial code of lamprey rar expression in both anterior CNS and pharynx resulting from dynamic and highly specific expression patterns during embryonic development. This ‘RAR code’ might function in regionalization and patterning of these two tissues by differentially modulating the expression of downstream effector genes during development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0016-4) contains supplementary material, which is available to authorized users.

The architecture of a scrambled genome reveals massive levels of genomic rearrangement during development

Programmed DNA rearrangements in the single-celled eukaryote Oxytricha trifallax completely rewire its germline into a somatic nucleus during development. This elaborate, RNA-mediated pathway eliminates noncoding DNA sequences that interrupt gene loci and reorganizes the remaining fragments by inversions and permutations to produce functional genes. Here, we report the Oxytricha germline genome and compare it to the somatic genome to present a global view of its massive scale of genome rearrangements. The remarkably encrypted genome architecture contains >3,500 scrambled genes, as well as >800 predicted germline-limited genes expressed, and some posttranslationally modified, during genome rearrangements. Gene segments for different somatic loci often interweave with each other. Single gene segments can contribute to multiple, distinct somatic loci. Terminal precursor segments from neighboring somatic loci map extremely close to each other, often overlapping. This genome assembly provides a draft of a scrambled genome and a powerful model for studies of genome rearrangement.

The lineage-specific evolution of aquaporin gene clusters facilitated tetrapod terrestrial adaptation

A major physiological barrier for aquatic organisms adapting to terrestrial life is dessication in the aerial environment. This barrier was nevertheless overcome by the Devonian ancestors of extant Tetrapoda, but the origin of specific molecular mechanisms that solved this water problem remains largely unknown. Here we show that an ancient aquaporin gene cluster evolved specifically in the sarcopterygian lineage, and subsequently diverged into paralogous forms of AQP2, -5, or -6 to mediate water conservation in extant Tetrapoda. To determine the origin of these apomorphic genomic traits, we combined aquaporin sequencing from jawless and jawed vertebrates with broad taxon assembly of >2,000 transcripts amongst 131 deuterostome genomes and developed a model based upon Bayesian inference that traces their convergent roots to stem subfamilies in basal Metazoa and Prokaryota. This approach uncovered an unexpected diversity of aquaporins in every lineage investigated, and revealed that the vertebrate superfamily consists of 17 classes of aquaporins (Aqp0 - Aqp16). The oldest orthologs associated with water conservation in modern Tetrapoda are traced to a cluster of three aqp2-like genes in Actinistia that likely arose >500 Ma through duplication of an aqp0-like gene present in a jawless ancestor. In sea lamprey, we show that aqp0 first arose in a protocluster comprised of a novel aqp14 paralog and a fused aqp01 gene. To corroborate these findings, we conducted phylogenetic analyses of five syntenic nuclear receptor subfamilies, which, together with observations of extensive genome rearrangements, support the coincident loss of ancestral aqp2-like orthologs in Actinopterygii. We thus conclude that the divergence of sarcopterygian-specific aquaporin gene clusters was permissive for the evolution of water conservation mechanisms that facilitated tetrapod terrestrial adaptation.