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

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

Emergence of a Homo sapiens-specific gene family and chromosome 16p11.2 CNV susceptibility.

Genetic differences that specify unique aspects of human evolution have typically been identified by comparative analyses between the genomes of humans and closely related primates, including more recently the genomes of archaic hominins. Not all regions of the genome, however, are equally amenable to such study. Recurrent copy number variation (CNV) at chromosome 16p11.2 accounts for approximately 1% of cases of autism and is mediated by a complex set of segmental duplications, many of which arose recently during human evolution. Here we reconstruct the evolutionary history of the locus and identify bolA family member 2 (BOLA2) as a gene duplicated exclusively in Homo sapiens. We estimate that a 95-kilobase-pair segment containing BOLA2 duplicated across the critical region approximately 282 thousand years ago (ka), one of the latest among a series of genomic changes that dramatically restructured the locus during hominid evolution. All humans examined carried one or more copies of the duplication, which nearly fixed early in the human lineage--a pattern unlikely to have arisen so rapidly in the absence of selection (P < 0.0097). We show that the duplication of BOLA2 led to a novel, human-specific in-frame fusion transcript and that BOLA2 copy number correlates with both RNA expression (r = 0.36) and protein level (r = 0.65), with the greatest expression difference between human and chimpanzee in experimentally derived stem cells. Analyses of 152 patients carrying a chromosome 16p11. rearrangement show that more than 96% of breakpoints occur within the H. sapiens-specific duplication. In summary, the duplicative transposition of BOLA2 at the root of the H. sapiens lineage about 282 ka simultaneously increased copy number of a gene associated with iron homeostasis and predisposed our species to recurrent rearrangements associated with disease.

Assessing Protein Sequence Database Suitability Using De Novo Sequencing.

The analysis of samples from unsequenced and/or understudied species as well as samples where the proteome is derived from multiple organisms poses two key questions. The first is whether the proteomic data obtained from an unusual sample type even contains peptide tandem mass spectra. The second question is whether an appropriate protein sequence database is available for proteomic searches. We describe the use of automated de novo sequencing for evaluating both the quality of a collection of tandem mass spectra and the suitability of a given protein sequence database for searching that data. Applications of this method include the proteome analysis of closely related species, metaproteomics, and proteomics of extinct organisms.

The African coelacanth genome provides insights into tetrapod evolution.

The discovery of a living coelacanth specimen in 1938 was remarkable, as this lineage of lobe-finned fish was thought to have become extinct 70 million years ago. The modern coelacanth looks remarkably similar to many of its ancient relatives, and its evolutionary proximity to our own fish ancestors provides a glimpse of the fish that first walked on land. Here we report the genome sequence of the African coelacanth, Latimeria chalumnae. Through a phylogenomic analysis, we conclude that the lungfish, and not the coelacanth, is the closest living relative of tetrapods. Coelacanth protein-coding genes are significantly more slowly evolving than those of tetrapods, unlike other genomic features. Analyses of changes in genes and regulatory elements during the vertebrate adaptation to land highlight genes involved in immunity, nitrogen excretion and the development of fins, tail, ear, eye, brain and olfaction. Functional assays of enhancers involved in the fin-to-limb transition and in the emergence of extra-embryonic tissues show the importance of the coelacanth genome as a blueprint for understanding tetrapod evolution.

Discovery of large genomic inversions using long range information.

BACKGROUND: Although many algorithms are now available that aim to characterize different classes of structural variation, discovery of balanced rearrangements such as inversions remains an open problem. This is mainly due to the fact that breakpoints of such events typically lie within segmental duplications or common repeats, which reduces the mappability of short reads. The algorithms developed within the 1000 Genomes Project to identify inversions are limited to relatively short inversions, and there are currently no available algorithms to discover large inversions using high throughput sequencing technologies. RESULTS: Here we propose a novel algorithm, VALOR, to discover large inversions using new sequencing methods that provide long range information such as 10X Genomics linked-read sequencing, pooled clone sequencing, or other similar technologies that we commonly refer to as long range sequencing. We demonstrate the utility of VALOR using both pooled clone sequencing and 10X Genomics linked-read sequencing generated from the genome of an individual from the HapMap project (NA12878). We also provide a comprehensive comparison of VALOR against several state-of-the-art structural variation discovery algorithms that use whole genome shotgun sequencing data. CONCLUSIONS: In this paper, we show that VALOR is able to accurately discover all previously identified and experimentally validated large inversions in the same genome with a low false discovery rate. Using VALOR, we also predicted a novel inversion, which we validated using fluorescent in situ hybridization. VALOR is available at https://github.com/BilkentCompGen/VALOR.

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.

Unusual Diversity of Myoglobin Genes in the Lungfish.

Myoglobin is a respiratory protein that serves as a model system in a variety of biological fields. Its main function is to deliver and store O2 in the heart and skeletal muscles, but myoglobin is also instrumental in homeostasis of nitric oxide (NO) and detoxification of reactive oxygen species (ROS). Almost every vertebrate harbors a single myoglobin gene; only some cyprinid fishes have two recently duplicated myoglobin genes. Here we show that the West African lungfish Protopterus annectens has at least seven distinct myoglobin genes (PanMb1-7), which diverged early in the evolution of lungfish and showed an enhanced evolutionary rate. These myoglobins are lungfish specific, and no other globin gene was found amplified. The myoglobins are differentially expressed in various lungfish tissues, and the brain is the main site of myoglobin expression. The typical myoglobin-containing tissues, the skeletal muscle and the heart, have much lower myoglobin mRNA levels. Muscle and heart express distinct myoglobins (PanMb1 and PanMb3, respectively). In cell culture, lungfish myoglobins improved cellular survival under hypoxia albeit with different efficiencies and reduced the production of reactive oxygen species. Only Mb2 and Mb6 enhanced the energy status of the cells. The unexpected diversity of myoglobin hints to a functional diversification of this gene: some myoglobins may have adapted to the O2 requirements of the specific tissue and help the lungfish to survive hypoxic periods; other myoglobins may have taken over the roles of neuroglobin and cytoglobin, which appear to be missing in the West African lungfish.

Morphological And Molecular Analyses of an Anatomical Novelty: The Pelvic Fin Filaments of the South American Lungfish.

The pelvic fins of male South American lungfish, Lepidosiren paradoxa, are adorned with a distinctive array of filaments, which grow and become highly vascularized during the breeding season. The resemblance between these pelvic fin filaments (PFFs) and external gills of other vertebrates suggested that this gill-like structure was used for physiological gas exchange. It has been proposed that the unique pelvic fin of male L. paradoxa is used for release of oxygen from its blood into the environment in order to aerate its nesting brood, or, conversely, as an auxiliary respiratory organ by absorbing oxygen from the environment into its bloodstream. Here, we employed histology, scanning electron microscopy (SEM) and quantitative PCR (qPCR) to assess whether the morphology and molecular profile of PFFs are compatible with a role in gas exchange. First, we closely examined its external morphology and showed that PFFs develop from short papillae during the rainy season, but remain covered by a thick nonvascularized epithelium. Histological examination confirmed that capillaries within the filaments are separated from the exterior by a basement membrane and a stratified epithelium composed of four to five cell layers. In addition, SEM analysis revealed significant differences between the fin filament epithelium and typical gill epithelium. Finally, our qPCR results showed that five genes commonly expressed in gills were downregulated in PFFs relative to their expression in regular pectoral fin epidermis. Collectively, our results do not directly support a role for PFFs, commonly referred to as "limb gills", in oxygen release or uptake.

Molecular cytogenetic differentiation of paralogs of Hox paralogs in duplicated and re-diploidized genome of the North American paddlefish (Polyodon spathula).

BACKGROUND: Acipenseriformes is a basal lineage of ray-finned fishes and comprise 27 extant species of sturgeons and paddlefishes. They are characterized by several specific genomic features as broad ploidy variation, high chromosome numbers, presence of numerous microchromosomes and propensity to interspecific hybridization. The presumed palaeotetraploidy of the American paddlefish was recently validated by molecular phylogeny and Hox genes analyses. A whole genome duplication in the paddlefish lineage was estimated at approximately 42 Mya and was found to be independent from several genome duplications evidenced in its sister lineage, i.e. sturgeons. We tested the ploidy status of available chromosomal markers after the expected rediploidization. Further we tested, whether paralogs of Hox gene clusters originated from this paddlefish specific genome duplication are cytogenetically distinguishable. RESULTS: We found that both paralogs HoxA alpha and beta were distinguishable without any overlapping of the hybridization signal - each on one pair of large metacentric chromosomes. Of the HoxD, only the beta paralog was unequivocally identified, whereas the alpha paralog did not work and yielded only an inconclusive diffuse signal. Chromosomal markers on three diverse ploidy levels reflecting different stages of rediploidization were identified: quadruplets retaining their ancestral tetraploid condition, semi-quadruplets still reflecting the ancestral tetraploidy with clear signs of advanced rediploidization, doublets were diploidized with ancestral tetraploidy already blurred. Also some of the available microsatellite data exhibited diploid allelic band patterns at their loci whereas another locus showed more than two alleles. CONCLUSIONS: Our exhaustive staining of paddlefish chromosomes combined with cytogenetic mapping of ribosomal genes and Hox paralogs and with microsatellite data, brings a closer look at results of the process of rediploidization in the course of paddlefish genome evolution. We show a partial rediploidization represented by a complex mosaic structure comparable with segmental paleotetraploidy revealed in sturgeons (Acipenseridae). Sturgeons and paddlefishes with their high propensity for whole genome duplication thus offer suitable animal model systems to further explore evolutionary processes that were shaping the early evolution of all vertebrates.

Establishing and maintaining primary cell cultures derived from the ctenophore Mnemiopsis leidyi.

We have developed an efficient method for the preparation and maintenance of primary cell cultures isolated from adult Mnemiopsis leidyi, a lobate ctenophore. Our primary cell cultures are derived from tissue explants or enzymatically dissociated cells, and maintained in a complex undefined ctenophore mesogleal serum. These methods can be used to isolate, maintain and visually monitor ctenophore cells to assess proliferation, cellular morphology and cell differentiation in future studies. Exemplar cell types that can be easily isolated from primary cultures include proliferative ectodermal and endodermal cells, motile amebocyte-like cells, and giant smooth muscle cells that exhibit inducible contractile properties. We have also derived 'tissue envelopes' containing sections of endodermal canal surrounded by mesoglea and ectoderm that can be used to monitor targeted cell types in an in vivo context. Access to efficient and reliably generated primary cell cultures will facilitate the analysis of ctenophore development, physiology and morphology from a cell biological perspective.

Genomic restructuring in the Tasmanian devil facial tumour: chromosome painting and gene mapping provide clues to evolution of a transmissible tumour.

Devil facial tumour disease (DFTD) is a fatal, transmissible malignancy that threatens the world's largest marsupial carnivore, the Tasmanian devil, with extinction. First recognised in 1996, DFTD has had a catastrophic effect on wild devil numbers, and intense research efforts to understand and contain the disease have since demonstrated that the tumour is a clonal cell line transmitted by allograft. We used chromosome painting and gene mapping to deconstruct the DFTD karyotype and determine the chromosome and gene rearrangements involved in carcinogenesis. Chromosome painting on three different DFTD tumour strains determined the origins of marker chromosomes and provided a general overview of the rearrangement in DFTD karyotypes. Mapping of 105 BAC clones by fluorescence in situ hybridisation provided a finer level of resolution of genome rearrangements in DFTD strains. Our findings demonstrate that only limited regions of the genome, mainly chromosomes 1 and X, are rearranged in DFTD. Regions rearranged in DFTD are also highly rearranged between different marsupials. Differences between strains are limited, reflecting the unusually stable nature of DFTD. Finally, our detailed maps of both the devil and tumour karyotypes provide a physical framework for future genomic investigations into DFTD.

SCPP genes in the coelacanth: tissue mineralization genes shared by sarcopterygians.

The coelacanth is the basal-most extant sarcopterygian that has teeth and tooth-like structures, comprising bone, dentin, and enamel or enameloid. Formation of these tissues involves many members of the secretory calcium-binding protein (SCPP) family. In tetrapods, acidic-residue-rich SCPPs are used in mineralization of bone and dentin, whereas Pro/Gln-rich SCPPs participate in enamel formation. Teleosts also employ many SCPPs for tissue mineralization. Nevertheless, the repertoire of SCPPs is largely different in teleosts and tetrapods; hence, filling this gap would be critical to elucidate early evolution of mineralized tissues in osteichthyans. In the present study, we searched for SCPP genes in the coelacanth genome and identified 11, of which two have clear orthologs in both tetrapods and teleosts, seven only in tetrapods, and two in neither of them. Given the divergence times of these vertebrate lineages, our discovery of this many SCPP genes shared between the coelacanth and tetrapods, but not with teleosts, suggests a complicated evolutionary scheme of SCPP genes in early osteichthyans. Our investigation also revealed both conserved and derived characteristics of SCPPs in the coelacanth and other vertebrates. Notably, acidic SCPPs independently evolved various acidic repeats in different lineages, while maintaining high acidity, presumably important for interactions with calcium. Furthermore, the three Pro/Gln-rich SCPP genes, required for mineralizing enamel matrix and confirmed only in tetrapods, were all identified in the coelacanth, strongly suggesting that enamel is equivalent in the coelacanth and tetrapods. This finding corroborates the previous proposition that true enamel evolved much earlier than the origin of tetrapods.

Atypical RNAs in the coelacanth transcriptome.

Circular and apparently trans-spliced RNAs have recently been reported as abundant types of transcripts in mammalian transcriptome data. Both types of non-colinear RNAs are also abundant in RNA-seq of different tissue from both the African and the Indonesian coelacanth. We observe more than 8,000 lincRNAs with normal gene structure and several thousands of circularized and trans-spliced products, showing that such atypical RNAs form a substantial contribution to the transcriptome. Surprisingly, the majority of the circularizing and trans-connecting splice junctions are unique to atypical forms, that is, are not used in normal isoforms.

Genome complexity in the coelacanth is reflected in its adaptive immune system.

We have analyzed the available genome and transcriptome resources from the coelacanth in order to characterize genes involved in adaptive immunity. Two highly distinctive IgW-encoding loci have been identified that exhibit a unique genomic organization, including a multiplicity of tandemly repeated constant region exons. The overall organization of the IgW loci precludes typical heavy chain class switching. A locus encoding IgM could not be identified either computationally or by using several different experimental strategies. Four distinct sets of genes encoding Ig light chains were identified. This includes a variant sigma-type Ig light chain previously identified only in cartilaginous fishes and which is now provisionally denoted sigma-2. Genes encoding α/ß and γ/δ T-cell receptors, and CD3, CD4, and CD8 co-receptors also were characterized. Ig heavy chain variable region genes and TCR components are interspersed within the TCR α/δ locus; this organization previously was reported only in tetrapods and raises questions regarding evolution and functional cooption of genes encoding variable regions. The composition, organization and syntenic conservation of the major histocompatibility complex locus have been characterized. We also identified large numbers of genes encoding cytokines and their receptors, and other genes associated with adaptive immunity. In terms of sequence identity and organization, the adaptive immune genes of the coelacanth more closely resemble orthologous genes in tetrapods than those in teleost fishes, consistent with current phylogenomic interpretations. Overall, the work reported described herein highlights the complexity inherent in the coelacanth genome and provides a rich catalog of immune genes for future investigations.

Evolution of adaptive immune recognition in jawless vertebrates.

All extant vertebrates possess an adaptive immune system wherein diverse immune receptors are created and deployed in specialized blood cell lineages. Recent advances in DNA sequencing and developmental resources for basal vertebrates have facilitated numerous comparative analyses that have shed new light on the molecular and cellular bases of immune defense and the mechanisms of immune receptor diversification in the "jawless" vertebrates. With data from these key species in hand, it is becoming possible to infer some general aspects of the early evolution of vertebrate adaptive immunity. All jawed vertebrates assemble their antigen-receptor genes through combinatorial recombination of different "diversity" segments into immunoglobulin or T-cell receptor genes. However, the jawless vertebrates employ an analogous, but independently derived set of immune receptors in order to recognize and bind antigens: the variable lymphocyte receptors (VLRs). The means by which this locus generates receptor diversity and achieves antigen specificity is of considerable interest because these mechanisms represent a completely independent strategy for building a large immune repertoire. Therefore, studies of the VLR system are providing insight into the fundamental principles and evolutionary potential of adaptive immune recognition systems. Here we review and synthesize the wealth of data that have been generated towards understanding the evolution of the adaptive immune system in the jawless vertebrates.

A living fossil in the genome of a living fossil: Harbinger transposons in the coelacanth genome.

Emerging data from the coelacanth genome are beginning to shed light on the origin and evolution of tetrapod genes and noncoding elements. Of particular relevance is the realization that coelacanth retains active copies of transposable elements that once served as raw material for the evolution of new functional sequences in the vertebrate lineage. Recognizing the evolutionary significance of coelacanth genome in this regard, we employed an ab initio search strategy to further classify its repetitive complement. This analysis uncovered a class of interspersed elements (Latimeria Harbinger 1-LatiHarb1) that is a major contributor to coelacanth genome structure and gene content (∼1% to 4% or the genome). Sequence analyses indicate that 1) each ∼8.7 kb LatiHarb1 element contains two coding regions, a transposase gene and a gene whose function is as yet unknown (MYB-like) and 2) copies of LatiHarb1 retain biological activity in the coelacanth genome. Functional analyses verify transcriptional and enhancer activities of LatiHarb1 in vivo and reveal transcriptional decoupling that could permit MYB-like genes to play functional roles not directly linked to transposition. Thus, LatiHarb1 represents the first known instance of a harbinger-superfamily transposon with contemporary activity in a vertebrate genome. Analyses of LatiHarb1 further corroborate the notion that exaptation of anciently active harbinger elements gave rise to at least two vertebrate genes (harbi1 and naif1) and indicate that the vertebrate gene tsnare1 also traces its ancestry to this transposon superfamily. Based on our analyses of LatiHarb1, we speculate that several functional features of harbinger elements may predispose the transposon superfamily toward recurrent exaptive evolution of cellular coding genes. In addition, these analyses further reinforce the broad utility of the coelacanth genome and other "outgroup" genomes in understanding the ancestry and evolution of vertebrate genes and genomes.

Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences.

We report a high-quality draft of the genome sequence of the grey, short-tailed opossum (Monodelphis domestica). As the first metatherian ('marsupial') species to be sequenced, the opossum provides a unique perspective on the organization and evolution of mammalian genomes. Distinctive features of the opossum chromosomes provide support for recent theories about genome evolution and function, including a strong influence of biased gene conversion on nucleotide sequence composition, and a relationship between chromosomal characteristics and X chromosome inactivation. Comparison of opossum and eutherian genomes also reveals a sharp difference in evolutionary innovation between protein-coding and non-coding functional elements. True innovation in protein-coding genes seems to be relatively rare, with lineage-specific differences being largely due to diversification and rapid turnover in gene families involved in environmental interactions. In contrast, about 20% of eutherian conserved non-coding elements (CNEs) are recent inventions that postdate the divergence of Eutheria and Metatheria. A substantial proportion of these eutherian-specific CNEs arose from sequence inserted by transposable elements, pointing to transposons as a major creative force in the evolution of mammalian gene regulation.

Mouse and zebrafish Hoxa3 orthologues have nonequivalent in vivo protein function.

Hox genes play evolutionarily conserved roles in specifying axial position during embryogenesis. A prevailing paradigm is that changes in Hox gene expression drive evolution of metazoan body plans. Conservation of Hox function across species, and among paralogous Hox genes within a species, supports a model of functional equivalence. In this report, we demonstrate that zebrafish hoxa3a (zfhoxa3a) expressed from the mouse Hoxa3 locus can substitute for mouse Hoxa3 in some tissues, but has distinct or null phenotypes in others. We further show, by using an allele encoding a chimeric protein, that this difference maps primarily to the zfhoxa3a C-terminal domain. Our data imply that the mouse and zebrafish proteins have diverged considerably since their last common ancestor, and that the major difference between them resides in the C-terminal domain. Our data further show that Hox protein function can evolve independently in different cell types or for specific functions. The inability of zfhoxa3a to perform all of the normal roles of mouse Hoxa3 illustrates that Hox orthologues are not always functionally interchangeable.