Deep ancestry of programmed genome rearrangement in lampreys

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.

The First Analysis of Synaptonemal Complexes in Jawless Vertebrates: Chromosome Synapsis and Transcription Reactivation at Meiotic Prophase I in the Lamprey Lampetra fluviatilis (Petromyzontiformes, Cyclostomata)

Transcription is known to be substage-specific in meiotic prophase I. If transcription is reactivated in the mid pachytene stage in mammals when synapsis is completed, then this process is observed in the zygotene stage in insects. The process of transcriptional reactivation has been studied in a small number of different taxa of invertebrates and vertebrates. Here, for the first time, we investigate synapsis and transcription in prophase I in the European river lamprey Lampetra fluviatilis (Petromyzontiformes, Cyclostomata), which is representative of jawless vertebrates that diverged from the main branch of vertebrates between 535 and 462 million years ago. We found that not all chromosomes complete synapsis in telomeric regions. Rounded structures were detected in chromatin and in some synaptonemal complexes, but their nature could not be determined conclusively. An analysis of RNA polymerase II distribution led to the conclusion that transcriptional reactivation in lamprey prophase I is not associated with the completion of chromosome synapsis. Monomethylated histone H3K4 is localized in meiotic chromatin throughout prophase I, and this pattern has not been previously detected in animals. Thus, the findings made it possible to identify synaptic and epigenetic patterns specific to this group and to expand knowledge about chromatin epigenetics in prophase I.

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.

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.

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.

Transcriptomic analysis of gonadal development in parasitic and non-parasitic lampreys (Ichthyomyzon spp.), with a comparison of genomic resources in these non-model species

Lampreys are jawless fishes that diverged ∼550 million years ago from other vertebrates. Sequencing of the somatic and the germline genomes of the sea lamprey (Petromyzon marinus) in 2013 and 2018, respectively, has helped to improve our understanding of the genes and gene networks that control many aspects of lamprey development. However, little is known about the genetic basis of gonadal differentiation in lampreys, partly due to the prolonged period during which their gonads remain sexually indeterminate. We performed RNA-sequencing on gonadal samples from four chestnut lamprey (Ichthyomyzon castaneus) and six northern brook lamprey (I. fossor) to identify differentially expressed genes (DEG’s) and pathways associated with transcriptomic differences in: (1) larvae during early gonadal differentiation versus definitive females (i.e., with oocytes in the slow cytoplasmic growth phase); and (2) females versus definitive males undergoing spermatogonial proliferation. We compared the mapping percentages of these transcriptomes to the two available sea lamprey reference genomes and three annotation files (Ensembl and UCSC for the somatic genome and SIMRbase for the germline genome). We found that mapping the RNA-seq reads to the germline genome gave superior results and, using Trinotate, we provided new putative annotations for 8161 genes in the somatic assembly and 880 genes for the germline assembly. We identified >2000 DEG’s between stages and sexes, as well as biological pathways associated with each. Interestingly, some of the upregulated genes (e.g., DEG’s associated with spermiation) suggest that changes in gene expression can precede morphological changes by several months. In contrast, only 81 DEG’s were evident between the chestnut lamprey (that remains sexually immature during an extended post-metamorphic parasitic feeding phase) and the nonparasitic northern brook lamprey (that undergoes sexual maturation near the end of metamorphosis), but few replicates were available for comparable stages and sexes. This work lays the foundation for identifying and confirming the orthology and the function of genes involved in gonadal development in these and other lamprey species across more developmental stages.

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.

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.

Miniscule differences between sex chromosomes in the giant genome of a salamander

In the Mexican axolotl (Ambystoma mexicanum), sex is determined by a single Mendelian factor, yet its sex chromosomes do not exhibit morphological differentiation typical of many vertebrate taxa that possess a single sex-determining locus. As sex chromosomes are theorized to differentiate rapidly, species with undifferentiated sex chromosomes provide the opportunity to reconstruct early events in sex chromosome evolution. Whole genome sequencing of 48 salamanders, targeted chromosome sequencing and in situ hybridization were used to identify the homomorphic sex chromosome that carries an A. mexicanum sex-determining factor and sequences that are present only on the W chromosome. Altogether, these sequences cover ~300 kb of validated female-specific (W chromosome) sequence, representing ~1/100,000^th of the 32 Gb genome. Notably, a recent duplication of ATRX, a gene associated with mammalian sex-determining pathways, is one of few functional (non-repetitive) genes identified among these W-specific sequences. This duplicated gene (ATRW) was used to develop highly predictive markers for diagnosing sex and represents a strong candidate for a recently-acquired sex determining locus (or sexually antagonistic gene) in A. mexicanum.

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 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.