Programmed loss of millions of base pairs from a vertebrate genome

Functional analyses of the polycomb-group genes in sea lamprey embryos undergoing programmed DNA loss

During early development, sea lamprey embryos undergo programmatic elimination of DNA from somatic progenitor cells in a process termed programmed genome rearrangement (PGR). Eliminated DNA eventually becomes condensed into micronuclei, which are then physically degraded and permanently lost from the cell. Previous studies indicated that many of the genes eliminated during PGR have mammalian homologs that are bound by polycomb repressive complex (PRC) in embryonic stem cells. To test whether PRC components play a role in the faithful elimination of germline-specific sequences, we used a combination of CRISPR/Cas9 and lightsheet microscopy to investigate the impact of gene knockouts on early development and the progression through stages of DNA elimination. Analysis of knockout embryos for the core PRC2 subunits EZH, SUZ12, and EED show that disruption of all three genes results in an increase in micronucleus number, altered distribution of micronuclei within embryos, and an increase in micronucleus volume in mutant embryos. While the upstream events of DNA elimination are not strongly impacted by loss of PRC2 components, this study suggests that PRC2 plays a role in the later stages of elimination related to micronucleus condensation and degradation. These findings also suggest that other genes/epigenetic pathways may work in parallel during DNA elimination to mediate chromatin structure, accessibility, and the ultimate loss of germline-specific DNA.

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.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

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.

Bioinformatic and fine-scale chromosomal mapping reveal the nature and evolution of eliminated chromosomes in the Japanese hagfish, Eptatretus burgeri, through analysis of repetitive DNA families

In the Japanese hagfish, Eptatretus burgeri, approximately 21% of the genomic DNA in germ cells (2n = 52) consists of 16 chromosomes (eliminated [E]-chromosomes) that are eliminated from presumptive somatic cells (2n = 36). To uncover the eliminated genome (E-genome), we have identified 16 eliminated repetitive DNA families from eight hagfish species, with 11 of these repeats being selectively amplified in the germline genome of E. burgeri. Furthermore, we have demonstrated that six of these sequences, namely EEEb1-6, are exclusively localized on all 16 E-chromosomes. This has led to the hypothesis that the eight pairs of E-chromosomes are derived from one pair of ancestral chromosomes via multiple duplication events over a prolonged evolutionary period. NGS analysis has recently facilitated the re-assembly of two distinct draft genomes of E. burgeri, derived from the testis and liver. This advancement allows for the prediction of not only nonrepetitive eliminated sequences but also over 100 repetitive and eliminated sequences, accomplished through K-mer-based analysis. In this study, we report four novel eliminated repetitive DNA sequences (designated as EEEb7-10) and confirm the relative chromosomal localization of all eliminated repeats (EEEb1-10) by fluorescence in situ hybridization (FISH). With the exception of EEEb10, all sequences were exclusively detected on EEEb1-positive chromosomes. Surprisingly, EEEb10 was detected as an intense signal on EEEb1-positive chromosomes and as a scattered signal on other chromosomes in germ cells. The study further divided the eight pairs of E-chromosomes into six groups based on the signal distribution of each DNA family, and fiber-FISH experiments showed that the EEEb2-10 family was dispersed in the EEEb1-positive extended chromatin fiber. These findings provide new insights into the mechanisms underlying chromosome elimination and the evolution of E-chromosomes, supporting our previous hypothesis.

Reduced steroid activation of elephant shark GR and MR after inserting four amino acids from the DNA-binding domain of lamprey corticoid receptor-1

Atlantic sea lamprey contains two corticoid receptors (CRs), CR1 and CR2, that have identical amino acid sequences, except for a four amino acid insert (Thr-Arg-Gln-Gly) in the CR1 DNA-binding domain (DBD). Steroids are stronger transcriptional activators of CR2 than of CR1 suggesting that the insert reduces the transcriptional response of lamprey CR1 to steroids. The DBD in elephant shark mineralocorticoid receptor (MR) and glucocorticoid receptor (GR), which are descended from a CR, lack these four amino acids, suggesting that a CR2 is their common ancestor. To determine if, similar to lamprey CR1, the presence of this insert in elephant shark MR and GR decreases transcriptional activation by corticosteroids, we inserted these four CR1-specific residues into the DBD of elephant shark MR and GR. Compared to steroid activation of wild-type elephant shark MR and GR, cortisol, corticosterone, aldosterone, 11-deoxycorticosterone and 11-deoxycortisol had lower transcriptional activation of these mutant MR and GR receptors, indicating that the absence of this four-residue segment in the DBD in wild-type elephant shark MR and GR increases transcriptional activation by corticosteroids.

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.

N-terminal domain influences steroid activation of the Atlantic sea lamprey corticoid receptor

Lampreys are jawless fish that evolved about 550 million years ago at the base of the vertebrate line. Modern lampreys contain a corticoid receptor (CR), the common ancestor of the glucocorticoid receptor (GR) and mineralocorticoid receptor (MR), which first appear in cartilaginous fish, such as sharks. Until recently, 344 amino acids at the amino terminus of adult lamprey CR were not present in the lamprey CR sequence in GenBank. A search of the recently sequenced lamprey germline genome identified two CR sequences, CR1 and CR2, containing the 344 previously un-identified amino acids. CR1 also contains a novel four amino acid insertion in the DNA-binding domain (DBD). We studied corticosteroid and progesterone activation of CR1 and CR2 and found their strongest response was to 11-deoxycorticosterone and 11-deoxycortisol, the two circulating corticosteroids in lamprey. Based on steroid specificity, both CRs are close to elephant shark MR and distant from elephant shark GR. HEK293 cells that were transfected with full-length CR1 or CR2 and the MMTV promoter have about 3-fold higher steroid-mediated activation compared to HEK293 cells transfected with these CRs and the TAT3 promoter. Deletion of the amino-terminal domain (NTD) of lamprey CR1 and CR2 to form truncated CRs decreased transcriptional activation by about 70% in HEK293 cells that were transfected with MMTV, but increased transcription by about 6-fold in cells transfected with TAT3. This indicated that the promoter has an important effect on NTD regulation of transcriptional activation of the CR by steroids. Our results also indicate that the entire lamprey CR sequence is needed for an accurate determination of steroid-mediated transcription.

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.

Novel selectively amplified DNA sequences in the germline genome of the Japanese hagfish, Eptatretus burgeri

In the Japanese hagfish Eptatretus burgeri, 16 chromosomes (eliminated [E]-chromosomes) have been lost in somatic cells (2n = 36), which is equivalent to approx. 21% of the genomic DNA in germ cells (2n = 52). At least seven of the 12 eliminated repetitive DNA families isolated in eight hagfish species were selectively amplified in the germline genome of this species. One of them, EEEb1 (eliminated element of E. burgeri 1) is exclusively localized on all E-chromosomes. Herein, we identified four novel eliminated repetitive DNA families (named EEEb3-6) through PCR amplification and suppressive subtractive hybridization (SSH) combined with Southern-blot hybridization. EEEb3 was mosaic for 5S rDNA and SINE elements. EEEb4 was GC-rich repeats and has one pair of direct and inverted repeats, whereas EEEb5 and EEEb6 were AT-rich repeats with one pair and two pairs of sub-repeats, respectively. Interestingly, all repeat classes except EEEb3 were transcribed in the testes, although no open reading frames (ORF) were identified. We conducted fluorescence in situ hybridization (FISH) to examine the chromosomal localizations of EEEb3-6 and EEEb2, which was previously isolated from the germline genome of E. burgeri. All sequences were only found on all EEEb1-positive E-chromosomes. Copy number estimation of the repeated elements by slot-blot hybridization revealed that (i) the EEEb1-6 family members occupied 39.9% of the total eliminated DNA, and (ii) a small number of repeats were retained in somatic cells, suggesting that there is incomplete elimination of the repeated elements. These results provide new insights into the mechanisms involved in the chromosome elimination and the evolution of E-chromosomes.

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.

Population genomics of New Zealand pouched lamprey (kanakana; piharau; Geotria australis)

Pouched lamprey (Geotria australis) or kanakana/piharau is a culturally and ecologically significant jawless fish that is distributed throughout Aotearoa New Zealand. Despite its importance, much remains unknown about historical relationships and gene flow between populations of this enigmatic species within New Zealand. To help inform management, we assembled a draft G. australis genome and completed the first comprehensive population genomics analysis of pouched lamprey within New Zealand using targeted gene sequencing (Cyt-b and COI) and restriction site-associated DNA sequencing (RADSeq) methods. Employing 16 000 genome-wide single nucleotide polymorphisms (SNPs) derived from RADSeq (n = 186) and sequence data from Cyt-b (766 bp, n = 94) and COI (589 bp, n = 20), we reveal low levels of structure across 10 sampling locations spanning the species range within New Zealand. F-statistics, outlier analyses, and STRUCTURE suggest a single panmictic population, and Mantel and EEMS tests reveal no significant isolation by distance. This implies either ongoing gene flow among populations or recent shared ancestry among New Zealand pouched lamprey. We can now use the information gained from these genetic tools to assist managers with monitoring effective population size, managing potential diseases, and conservation measures such as artificial propagation programs. We further demonstrate the general utility of these genetic tools for acquiring information about elusive species.

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.

Evolutionary Perspectives on Germline-Restricted Chromosomes in Flies (Diptera)

In some eukaryotes, germline soma differentiation involves elimination of parts of the genome from somatic cells. The portions of the genome restricted to the germline often contain genes that play a role in development and function of the germline. Lineages with germline-restricted DNA are taxonomically diverse, and the size of the germline-restricted genome varies substantially. Unfortunately, few of these lineages have been studied in detail. As a result, we understand little about the general evolutionary forces that drive the origin and maintenance of germline-restricted DNA. One of the taxonomic groups where germline-restricted DNA has been poorly studied are the flies (Diptera). In three Dipteran families, Chironomidae, Cecidomyiidae, and Sciaridae, entire chromosomes are eliminated from somatic cells early in embryonic development. Germline-restricted chromosomes are thought to have evolved independently in the Dipteran families and their size, number, and transmission patterns vary between families. Although there is a wealth of cytological studies on these chromosomes in flies, almost no genomic studies have been undertaken. As a result, very little is known about how and why they evolved and what genes they encode. This review summarizes the literature on germline-restricted chromosomes in Diptera, discusses hypotheses for their origin and function, and compares germline-restricted DNA in Diptera to other eukaryotes. Finally, we discuss why Dipteran lineages represent a promising system for the study of germline-restricted chromosomes and propose future avenues of research on this topic.

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.

The Cranial Neural Crest in a Multiomics Era

Neural crest ontogeny plays a prominent role in craniofacial development. In this Perspective article, we discuss recent advances to the understanding of mechanisms underlying the cranial neural crest gene regulatory network (cNC-GRN) stemming from omics-based studies. We briefly summarize how parallel considerations of transcriptome, interactome, and epigenome data significantly elaborated the roles of key players derived from pre-omics era studies. Furthermore, the growing cohort of cNC multiomics data revealed contribution of the non-coding genomic landscape. As technological improvements are constantly being developed, we reflect on key questions we are poised to address by taking advantage of the unique perspective a multiomics approach has to offer.

Activated Erk Is an Early Retrograde Signal After Spinal Cord Injury in the Lamprey

We previously reported that spinal cord transection (TX) in the lamprey causes mRNA to accumulate in the injured tips of large reticulospinal (RS) axons. We sought to determine whether this mRNA accumulation results from phosphorylation and transport of retrograde signals, similar to what has been reported in mammalian peripheral nerve. Extracellular signal-regulated protein kinase (Erk), mediates the neurite outgrowth-promoting effects of many neurotrophic factors. To assess the role of Erk in retrograde signaling of RS axon injury, we used immunoblot and immunohistochemistry to determine the changes in phosphorylated Erk (p-Erk) in the spinal cord after spinal cord TX. Immunostaining for p-Erk increased within axons and local cell bodies, most heavily within the 1-2 mm closest to the TX site, at between 3 and 6 h post-TX. In axons, p-Erk was concentrated in 3-5 μm granules that became less numerous with distance from the TX. The retrograde molecular motor dynein colocalized with p-Erk, but vimentin, which in peripheral nerve was reported to participate with p-Erk as part of a retrograde signal complex, did not colocalize with p-Erk, even though vimentin levels were elevated post-TX. The results suggest that p-Erk, but not vimentin, may function as a retrograde axotomy signal in lamprey central nervous system neurons, and that this signal may induce transcription of mRNA, which is then transported down the axon to its injured tip.

Evolutionary history of histamine receptors: Early vertebrate origin and expansion of the H3-H4 subtypes

Histamine receptors belonging to the superfamily of G protein-coupled receptors (GPCRs) mediate the diverse biological effects of biogenic histamine. They are classified into four phylogenetically distinct subtypes H₁-H₄, each with a different binding affinity for histamine and divergent downstream signaling pathways. Here we present the evolutionary history of the histamine receptors using a phylogenetic approach complemented with comparative genomics analyses of the sequences, gene structures, and synteny of gene neighborhoods. The data indicate the earliest emergence of histamine-mediated GPCR signaling by a H₂ in a prebilaterian ancestor. The analyses support a revised classification of the vertebrate H₃-H₄ receptor subtypes. We demonstrate the presence of the H₄ across vertebrates, contradicting the currently held notion that H₄ is restricted to mammals. These non-mammalian vertebrate H₄ orthologs have been mistaken for H₃. We also identify the presence of a new H₃ subtype (H₃B), distinct from the canonical H₃ (H₃A), and propose that the H₃A, H₃B, and H₄ likely emerged from a H₃ progenitor through the 1R/2R whole genome duplications in an ancestor of the vertebrates. It is apparent that the ability of the H₁, H₂, and H_3-4 to bind histamine was acquired convergently. We identified genomic signatures suggesting that the H₁ and H₃-H₄ shared a last common ancestor with the muscarinic receptor in a bilaterian predecessor whereas, the H₂ and the α-adrenoreceptor shared a progenitor in a prebilaterian ancestor. Furthermore, site-specific analysis of the vertebrate subtypes revealed potential residues that may account for the functional divergence between them.

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.

Micronuclei in germ cells of hybrid frogs from Pelophylax esculentus complex contain gradually eliminated chromosomes

In most organisms, cells typically maintain genome integrity, as radical genome reorganization leads to dramatic consequences. However, certain organisms, ranging from unicellular ciliates to vertebrates, are able to selectively eliminate specific parts of their genome during certain stages of development. Moreover, partial or complete elimination of one of the parental genomes occurs in interspecies hybrids reproducing asexually. Although several examples of this phenomenon are known, the molecular and cellular processes involved in selective elimination of genetic material remain largely undescribed for the majority of such organisms. Here, we elucidate the process of selective genome elimination in water frog hybrids from the Pelophylax esculentus complex reproducing through hybridogenesis. Specifically, in the gonads of diploid and triploid hybrids, but not those of the parental species, we revealed micronuclei in the cytoplasm of germ cells. In each micronucleus, only one centromere was detected with antibodies against kinetochore proteins, suggesting that each micronucleus comprises a single chromosome. Using 3D-FISH with species-specific centromeric probe, we determined the role of micronuclei in selective genome elimination. We found that in triploid LLR hybrids, micronuclei preferentially contain P. ridibundus chromosomes, while in diploid hybrids, micronuclei preferentially contain P. lessonae chromosomes. The number of centromere signals in the nuclei suggested that germ cells were aneuploid until they eliminate the whole chromosomal set of one of the parental species. Furthermore, in diploid hybrids, misaligned P. lessonae chromosomes were observed during the metaphase stage of germ cells division, suggesting their possible elimination due to the inability to attach to the spindle and segregate properly. Additionally, we described gonocytes with an increased number of P. ridibundus centromeres, indicating duplication of the genetic material. We conclude that selective genome elimination from germ cells of diploid and triploid hybrids occurs via the gradual elimination of individual chromosomes of one of the parental genomes, which are enclosed within micronuclei.

Intron and gene size expansion during nervous system evolution

Background

The evolutionary radiation of animals was accompanied by extensive expansion of gene and genome sizes, increased isoform diversity, and complexity of regulation.

Results

Here we show that the longest genes are enriched for expression in neuronal tissues of diverse vertebrates and of invertebrates. Additionally, we show that neuronal gene size expansion occurred predominantly through net gains in intron size, with a positional bias toward the 5′ end of each gene.

Conclusions

We find that intron and gene size expansion is a feature of many genes whose expression is enriched in nervous systems. We speculate that unique attributes of neurons may subject neuronal genes to evolutionary forces favoring net size expansion. This process could be associated with tissue-specific constraints on gene function and/or the evolution of increasingly complex gene regulation in nervous systems.

Regulation of Genomic Output and (Pluri)potency in Regeneration

Regeneration is a remarkable phenomenon that has been the subject of awe and bafflement for hundreds of years. Although regeneration competence is found in highly divergent organisms throughout the animal kingdom, recent advances in tools used for molecular and genomic characterization have uncovered common genes, molecular mechanisms, and genomic features in regenerating animals. In this review we focus on what is known about how genome regulation modulates cellular potency during regeneration. We discuss this regulation in the context of complex tissue regeneration in animals, from Hydra to humans, with reference to ex vivo-cultured cell models of pluripotency when appropriate. We emphasize the importance of a detailed molecular understanding of both the mechanisms that regulate genomic output and the functional assays that assess the biological relevance of such molecular characterizations.

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.

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.

Genomic Organization of the Repetitive Sequences in Centropomus undecimalis (Perciformes, Centropomidae): Implications for Hybridization and Aquaculture Programs

The family Centropomidae includes a number of fish species of high commercial value. One of these species, Centropomus undecimalis, is a target of artisanal, industrial, and sports fisheries and has also considerable potential for captive breeding, which has led to its inclusion in several aquaculture programs. While the biology and ecology of C. undecimalis are relatively well documented, few karyological data are available on this species, and they are still scarce for other centropomids. The few chromosomal data available on this family indicate a conserved karyotype 2n = 48, but it is unclear whether the chromosome microstructure is also conserved. In this study, new cytogenetic data are presented on C. undecimalis from the Amazon coastal zone, including C-banding, Ag-NOR, in situ hybridization with repetitive DNA probes (5S and 18S ribosomal genes), and telomeric (TTAGGG)n sequences. The diploid number of the species was 2n = 48, with heterochromatic blocks in the centromeric and pericentromeric regions, as well as distal signals; the nucleolus organizer regions (NORs) were associated with the heterochromatic region. The 18S and 5S recombinant DNA (rDNA) clusters were located in the distal region of chromosome pairs 1 and 11, respectively. The similarities of the karyotype macrostructure found among the centropomid species reinforce their exceptional chromosomal stability. However, the presence of heterochromatic blocks and location of NORs suggest the occurrence of structural rearrangements, which indicates that evolutionary dynamics at the microstructural level in this group may be relatively complex and should be evaluated carefully in any study that targets the production of hybrids for aquaculture.

The protective function of noncoding DNA in genome defense of eukaryotic male germ cells

Peripheral and abundant noncoding DNA has been hypothesized to protect the genome and the central protein-coding sequences against DNA damage in somatic genome. In the cytosol, invading exogenous nucleic acids may first be deactivated by small RNAs encoded by noncoding DNA via mechanisms similar to the prokaryotic CRISPR-Cas system. In the nucleus, the radicals generated by radiation in the cytosol, radiation energy and invading exogenous nucleic acids are absorbed, blocked and/or reduced by peripheral heterochromatin, and damaged DNA in heterochromatin is removed and excluded from the nucleus to the cytoplasm through nuclear pore complexes. To further strengthen the hypothesis, this review summarizes the experimental evidence supporting the protective function of noncoding DNA in the genome of male germ cells. Based on these data, this review provides evidence supporting the protective role of noncoding DNA in the genome defense of sperm genome through similar mechanisms to those of the somatic genome.

First Report of a Mitochondrial Pseudogene in Agnathan Vertebrates (Cyclostomata: Petromyzontidae)

We report herein the characterization of a nuclear paralog of a fragment of the mitochondrial genome (a numt) in two closely related species of lampreys (Ichthyomyzon spp.). Although numts have been characterized in several vertebrate taxa, numts have yet to be reported for fishes in general. Given the phylogenetic position of lampreys relative to other vertebrates, the presence of numts within the lamprey genome is either evidence of an ancestral trait lost in other fishes but uniquely retained in agnathans and amniotes, or (more intriguingly) a product of the genome rearrangements these animals undergo during development.

Expansions, diversification, and interindividual copy number variations of AID/APOBEC family cytidine deaminase genes in lampreys

Cytidine deaminases of the AID/APOBEC family mutate the genetic material of pathogens or contribute to the generation and diversification of antibody repertoires in jawed vertebrates. In the extant jawless vertebrate, the lamprey, two members of the AID/APOBEC family are implicated in the somatic diversification of variable lymphocyte receptor (VLR) repertoires. We discovered an unexpected diversity of cytidine deaminase genes within and among lamprey species. The cytidine deaminases with features comparable to jawed vertebrate AID are always present, suggesting that they are involved in essential processes, such as VLR assembly. In contrast, other genes show a remarkable copy number variation, like the APOBEC3 genes in mammals. This suggests an unexpected similarity in functional deployment of AID/APOBEC cytidine deaminases across all vertebrates.

Cytidine deaminases of the AID/APOBEC family catalyze C-to-U nucleotide transitions in mRNA or DNA. Members of the APOBEC3 branch are involved in antiviral defense, whereas AID contributes to diversification of antibody repertoires in jawed vertebrates via somatic hypermutation, gene conversion, and class switch recombination. In the extant jawless vertebrate, the lamprey, two members of the AID/APOBEC family are implicated in the generation of somatic diversity of the variable lymphocyte receptors (VLRs). Expression studies linked CDA1 and CDA2 genes to the assembly of VLRA/C genes in T-like cells and the VLRB genes in B-like cells, respectively. Here, we identify and characterize several CDA1-like genes in the larvae of different lamprey species and demonstrate that these encode active cytidine deaminases. Structural comparisons of the CDA1 variants highlighted substantial differences in surface charge; this observation is supported by our finding that the enzymes require different conditions and substrates for optimal activity in vitro. Strikingly, we also found that the number of CDA-like genes present in individuals of the same species is variable. Nevertheless, irrespective of the number of different CDA1-like genes present, all lamprey larvae have at least one functional CDA1-related gene encoding an enzyme with predicted structural and chemical features generally comparable to jawed vertebrate AID. Our findings suggest that, similar to APOBEC3 branch expansion in jawed vertebrates, the AID/APOBEC family has undergone substantial diversification in lamprey, possibly indicative of multiple distinct biological roles.

The Future of Multiplexed Eukaryotic Genome Engineering

Multiplex genome editing is the simultaneous introduction of multiple distinct modifications to a given genome. Though in its infancy, maturation of this field will facilitate powerful new biomedical research approaches and will enable a host of far-reaching biological engineering applications, including new therapeutic modalities and industrial applications, as well as "genome writing" and de-extinction efforts. In this Perspective, we focus on multiplex editing of large eukaryotic genomes. We describe the current state of multiplexed genome editing, the current limits of our ability to multiplex edits, and provide perspective on the many applications that fully realized multiplex editing technologies would enable in higher eukaryotic genomes. We offer a broad look at future directions, covering emergent CRISPR-based technologies, advances in intracellular delivery, and new DNA assembly approaches that may enable future genome editing on a massively multiplexed scale.

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.

Exploration of the Germline Genome of the Ciliate Chilodonella uncinata through Single-Cell Omics (Transcriptomics and Genomics)

Separate germline and somatic genomes are found in numerous lineages across the eukaryotic tree of life, often separated into distinct tissues (e.g., in plants, animals, and fungi) or distinct nuclei sharing a common cytoplasm (e.g., in ciliates and some foraminifera). In ciliates, germline-limited (i.e., micronuclear-specific) DNA is eliminated during the development of a new somatic (i.e., macronuclear) genome in a process that is tightly linked to large-scale genome rearrangements, such as deletions and reordering of protein-coding sequences. Most studies of germline genome architecture in ciliates have focused on the model ciliates Oxytricha trifallax, Paramecium tetraurelia, and Tetrahymena thermophila, for which the complete germline genome sequences are known. Outside of these model taxa, only a few dozen germline loci have been characterized from a limited number of cultivable species, which is likely due to difficulties in obtaining sufficient quantities of “purified” germline DNA in these taxa. Combining single-cell transcriptomics and genomics, we have overcome these limitations and provide the first insights into the structure of the germline genome of the ciliate Chilodonella uncinata, a member of the understudied class Phyllopharyngea. Our analyses reveal the following: (i) large gene families contain a disproportionate number of genes from scrambled germline loci; (ii) germline-soma boundaries in the germline genome are demarcated by substantial shifts in GC content; (iii) single-cell omics techniques provide large-scale quality germline genome data with limited effort, at least for ciliates with extensively fragmented somatic genomes. Our approach provides an efficient means to understand better the evolution of genome rearrangements between germline and soma in ciliates.

Our understanding of the distinctions between germline and somatic genomes in ciliates has largely relied on studies of a few model genera (e.g., Oxytricha, Paramecium, Tetrahymena). We have used single-cell omics to explore germline-soma distinctions in the ciliate Chilodonella uncinata, which likely diverged from the better-studied ciliates ~700 million years ago. The analyses presented here indicate that developmentally regulated genome rearrangements between germline and soma are demarcated by rapid transitions in local GC composition and lead to diversification of protein families. The approaches used here provide the basis for future work aimed at discerning the evolutionary impacts of germline-soma distinctions among diverse ciliates.

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.

Evolution of the vertebrate insulin receptor substrate (Irs) gene family

BACKGROUND

Insulin receptor substrate (Irs) proteins are essential for insulin signaling as they allow downstream effectors to dock with, and be activated by, the insulin receptor. A family of four Irs proteins have been identified in mice, however the gene for one of these, IRS3, has been pseudogenized in humans. While it is known that the Irs gene family originated in vertebrates, it is not known when it originated and which members are most closely related to each other. A better understanding of the evolution of Irs genes and proteins should provide insight into the regulation of metabolism by insulin.

RESULTS

Multiple genes for Irs proteins were identified in a wide variety of vertebrate species. Phylogenetic and genomic neighborhood analyses indicate that this gene family originated very early in vertebrae evolution. Most Irs genes were duplicated and retained in fish after the fish-specific genome duplication. Irs genes have been lost of various lineages, including Irs3 in primates and birds and Irs1 in most fish. Irs3 and Irs4 experienced an episode of more rapid protein sequence evolution on the ancestral mammalian lineage. Comparisons of the conservation of the proteins sequences among Irs paralogs show that domains involved in binding to the plasma membrane and insulin receptors are most strongly conserved, while divergence has occurred in sequences involved in interacting with downstream effector proteins.

CONCLUSIONS

The Irs gene family originated very early in vertebrate evolution, likely through genome duplications, and in parallel with duplications of other components of the insulin signaling pathway, including insulin and the insulin receptor. While the N-terminal sequences of these proteins are conserved among the paralogs, changes in the C-terminal sequences likely allowed changes in biological function.

Early vertebrate origin and diversification of small transmembrane regulators of cellular ion transport

KEY POINTS

Small transmembrane proteins such as FXYDs, which interact with Na⁺ ,K⁺ -ATPase, and the micropeptides that interact with sarco/endoplasmic reticulum Ca²⁺ -ATPase play fundamental roles in regulation of ion transport in vertebrates. Uncertain evolutionary origins and phylogenetic relationships among these regulators of ion transport have led to inconsistencies in their classification across vertebrate species, thus hampering comparative studies of their functions. We discovered the first FXYD homologue in sea lamprey, a basal jawless vertebrate, which suggests small transmembrane regulators of ion transport emerged early in the vertebrate lineage. We also identified 13 gene subfamilies of FXYDs and propose a revised, phylogeny-based FXYD classification that is consistent across vertebrate species. These findings provide an improved framework for investigating physiological and pathophysiological functions of small transmembrane regulators of ion transport.

ABSTRACT

Small transmembrane proteins are important for regulation of cellular ion transport. The most prominent among these are members of the FXYD family (FXYD1-12), which regulate Na⁺ ,K⁺ -ATPase, and phospholamban, sarcolipin, myoregulin and DWORF, which regulate the sarco/endoplasmic reticulum Ca²⁺ -ATPase (SERCA). FXYDs and regulators of SERCA are present in fishes, as well as terrestrial vertebrates; however, their evolutionary origins and phylogenetic relationships are obscure, thus hampering comparative physiological studies. Here we discovered that sea lamprey (Petromyzon marinus), a representative of extant jawless vertebrates (Cyclostomata), expresses an FXYD homologue, which strongly suggests that FXYDs predate the emergence of fishes and other jawed vertebrates (Gnathostomata). Using a combination of sequence-based phylogenetic analysis and conservation of local chromosome context, we determined that FXYDs markedly diversified in the lineages leading to cartilaginous fishes (Chondrichthyes) and bony vertebrates (Euteleostomi). Diversification of SERCA regulators was much less extensive, indicating they operate under different evolutionary constraints. Finally, we found that FXYDs in extant vertebrates can be classified into 13 gene subfamilies, which do not always correspond to the established FXYD classification. We therefore propose a revised classification that is based on evolutionary history of FXYDs and that is consistent across vertebrate species. Collectively, our findings provide an improved framework for investigating the function of ion transport in health and disease.

Molecular cytogenetic characterization of chromosome site-specific repetitive sequences in the Arctic lamprey (Lethenteron camtschaticum, Petromyzontidae)

All extant lamprey karyotypes are characterized by almost all dot-shaped microchromosomes. To understand the molecular basis of chromosome structure in lampreys, we performed chromosome C-banding and silver staining and chromosome mapping of the 18S–28S and 5S ribosomal RNA (rRNA) genes and telomeric TTAGGG repeats in the Arctic lamprey (Lethenteron camtschaticum). In addition, we cloned chromosome site-specific repetitive DNA sequences and characterized them by nucleotide sequencing, chromosome in situ hybridization, and filter hybridization. Three types of repetitive sequences were detected; a 200-bp AT-rich repetitive sequence, LCA-EcoRIa that co-localized with the 18S–28S rRNA gene clusters of 3 chromosomal pairs; a 364-bp AT-rich LCA-EcoRIb sequence that showed homology to the EcoRI sequence family from the sea lamprey (Petromyzon marinus), which contains short repeats as centromeric motifs; and a GC-rich 702-bp LCA-ApaI sequence that was distributed on nearly all chromosomes and showed significant homology with the integrase-coding region of a Ty3/Gypsy family long terminal repeat (LTR) retrotransposon. All three repetitive sequences are highly conserved within the Petromyzontidae or within Petromyzontidae and Mordaciidae. Molecular cytogenetic characterization of these site-specific repeats showed that they may be correlated with programed genome rearrangement (LCA-EcoRIa), centromere structure and function (LCA-EcoRIb), and site-specific amplification of LTR retroelements through homogenization between non-homologous chromosomes (LCA-ApaI).

DNA Remodeling by Strict Partial Endoreplication in Orchids, an Original Process in the Plant Kingdom

DNA remodeling during endoreplication appears to be a strong developmental characteristic in orchids. In this study, we analyzed DNA content and nuclei in 41 species of orchids to further map the genome evolution in this plant family. We demonstrate that the DNA remodeling observed in 36 out of 41 orchids studied corresponds to strict partial endoreplication. Such process is developmentally regulated in each wild species studied. Cytometry data analyses allowed us to propose a model where nuclear states 2C, 4E, 8E, etc. form a series comprising a fixed proportion, the euploid genome 2C, plus 2-32 additional copies of a complementary part of the genome. The fixed proportion ranged from 89% of the genome in Vanilla mexicana down to 19% in V. pompona, the lowest value for all 148 orchids reported. Insterspecific hybridization did not suppress this phenomenon. Interestingly, this process was not observed in mass-produced epiphytes. Nucleolar volumes grow with the number of endocopies present, coherent with high transcription activity in endoreplicated nuclei. Our analyses suggest species-specific chromatin rearrangement. Towards understanding endoreplication, V. planifolia constitutes a tractable system for isolating the genomic sequences that confer an advantage via endoreplication from those that apparently suffice at diploid level.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Identification of three somatostatin genes in lampreys

Somatostatins (SSs) are a structurally diverse family of neuropeptides that play important roles in the regulation of growth, development and metabolism in vertebrates. It has been recently proposed that the common ancestor of gnathostomes possessed three SS genes, namely SS1, SS2 and SS5. SS1 and SS2 are still present in most extant gnathostome species investigated so far while SS5 primarily occurs in chondrichthyes, actinopterygians and actinistia but not in tetrapods. Very little is known about the repertoire of SSs in cyclostomes, which are extant jawless vertebrates. In the present study, we report the cloning of the cDNAs encoding three distinct lamprey SS variants that we call SSa, SSb and SSc. SSa and SSb correspond to the two SS variants previously characterized in lamprey, while SSc appears to be a totally novel one. SSa exhibits the same sequence as gnathostome SS1. SSb differs from SSa by only one substitution (Thr¹²→Ser). SSc exhibits a totally unique structure (ANCRMFYWKTMAAC) that shares only 50% identity with SSa and SSb. SSa, SSb and SSc precursors do not exhibit any appreciable sequence similarity outside the C-terminal region containing the SS sequence. Phylogenetic analyses failed to clearly assign orthology relationships between lamprey and gnathostome SS genes. Synteny analysis suggests that the SSc gene arose before the split of the three gnathostome genes SS1, SS2 and SS5.

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.

Evolution and expression of the phosphodiesterase 6 genes unveils vertebrate novelty to control photosensitivity

Background

Phosphodiesterase 6 (PDE6) is a protein complex that hydrolyses cGMP and acts as the effector of the vertebrate phototransduction cascade. The PDE6 holoenzyme consists of catalytic and inhibitory subunits belonging to two unrelated gene families. Rods and cones express distinct genes from both families: PDE6A and PDE6B code for the catalytic and PDE6G the inhibitory subunits in rods while PDE6C codes for the catalytic and PDE6H the inhibitory subunits in cones. We performed phylogenetic and comparative synteny analyses for both gene families in genomes from a broad range of animals. Furthermore, gene expression was investigated in zebrafish.

Results

We found that both gene families expanded from one to three members in the two rounds of genome doubling (2R) that occurred at the base of vertebrate evolution. The PDE6 inhibitory subunit gene family appears to be unique to vertebrates and expanded further after the teleost-specific genome doubling (3R). We also describe a new family member that originated in 2R and has been lost in amniotes, which we have named pde6i. Zebrafish has retained two additional copies of the PDE6 inhibitory subunit genes after 3R that are highly conserved, have high amino acid sequence identity, are coexpressed in the same photoreceptor type as their amniote orthologs and, interestingly, show strikingly different daily oscillation in gene expression levels.

Conclusions

Together, these data suggest specialisation related to the adaptation to different light intensities during the day-night cycle, most likely maintaining the regulatory function of the PDE inhibitory subunits in the phototransduction cascade.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0695-z) contains supplementary material, which is available to authorized users.

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.