Genome size, GC percentage and 5mC level in the Indonesian coelacanth Latimeria menadoensis

The living fossil Latimeria menadoensis is important to understand sarcopterygian evolution. To gain further insights into this fish species we studied its genome size, GC% and 5mC level. The genome size and the GC% of the Indonesian coelacanth seem to be very similar to those of the African coelacanth. Moreover the GC%, the CpG frequency and the 5mC level of L. menadoensis are more similar to those of fish and amphibians than to those of mammals, birds and reptiles and this is in line with the hypothesis that two different DNA methylation and CpG shortage equilibria arose during vertebrate evolution. Our results suggest that the genome of L. menadoensis has remained unchanged for several million years, maybe since the origin of the lineage which from lobe-finned fish led to tetrapods. These data fit a conservative evolutionary landscape and suggest that the genome of the extant crossopterygians may be a sort of evolutionarily frozen genome.     

An Endogenous Foamy-like Viral Element in the Coelacanth Genome

Little is known about the origin and long-term evolutionary mode of retroviruses. Retroviruses can integrate into their hosts'' genomes, providing a molecular fossil record for studying their deep history. Here we report the discovery of an endogenous foamy virus-like element, which we designate ‘coelacanth endogenous foamy-like virus’ (CoeEFV), within the genome of the coelacanth (Latimeria chalumnae). Phylogenetic analyses place CoeEFV basal to all known foamy viruses, strongly suggesting an ancient ocean origin of this major retroviral lineage, which had previously been known to infect only land mammals. The discovery of CoeEFV reveals the presence of foamy-like viruses in species outside the Mammalia. We show that foamy-like viruses have likely codiverged with their vertebrate hosts for more than 407 million years and underwent an evolutionary transition from water to land with their vertebrate hosts. These findings suggest an ancient marine origin of retroviruses and have important implications in understanding foamy virus biology.     

Conservation of globin genes in the “living fossil” Latimeria chalumnae and reconstruction of the evolution of the vertebrate globin family

The (hemo-)globins are among the best-investigated proteins in biomedical sciences. These small heme-proteins play an important role in oxygen supply, but may also have other functions. In addition to well known hemoglobin and myoglobin, six other vertebrate globin types have been identified in recent years: neuroglobin, cytoglobin, globin E, globin X, globin Y, and androglobin. Analyses of the genome of the “living fossil” Latimeria chalumnae show that the coelacanth is the only known vertebrate that includes all eight globin types. Thus, Latimeria can also be considered as a “globin fossil”. Analyses of gene synteny and phylogenetic reconstructions allow us to trace the evolution and the functional changes of the vertebrate globin family. Neuroglobin and globin X diverged from the other globin types before the separation of Protostomia and Deuterostomia. The cytoglobins, which are unlikely to be involved in O₂ supply, form the earliest globin branch within the jawed vertebrates (Gnathostomata), but do not group with the agnathan hemoglobins, as it has been proposed before. There is strong evidence from phylogenetic reconstructions and gene synteny that the eye-specific globin E and muscle-specific myoglobin constitute a common clade, suggesting a similar role in intracellular O₂ supply. Latimeria possesses two α- and two β-hemoglobin chains, of which one α-chain emerged prior to the divergence of Actinopterygii and Sarcopterygii, but has been retained only in the coelacanth. Notably, the embryonic hemoglobin α-chains of Gnathostomata derive from a common ancestor, while the embryonic β-chains – with the exception of a more complex pattern in the coelacanth and amphibians – display a clade-specific evolution. Globin Y is associated with the hemoglobin gene cluster, but its phylogenetic position is not resolved. Our data show an early divergence of distinct globin types in the vertebrate evolution before the emergence of tetrapods. The subsequent loss of globins in certain taxa may be associated with changes in the oxygen-dependent metabolism. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Predicting suitable environments and potential occurrences for coelacanths (Latimeria spp.)

Extant coelacanths (Latimeria chalumnae) were first discovered in the western Indian Ocean in 1938; in 1998, a second species of coelacanth, Latimeria menadoensis, was discovered off the north coast of Sulawesi, Indonesia, expanding the known distribution of the genus across the Indian Ocean Basin. This study uses ecological niche modeling techniques to estimate dimensions of realized niches of coelacanths and generate hypotheses for additional sites where they might be found. Coelacanth occurrence information was integrated with environmental and oceanographic data using the Genetic Algorithm for Rule-set Production (GARP) and a maximum entropy algorithm (Maxent). Resulting models were visualized as maps of relative suitability of sites for coelacanths throughout the Indian Ocean, as well as scatterplots of ecological variables. Our findings suggest that the range of coelacanths could extend beyond their presently known distribution and suggests alternative mechanisms for currently observed distributions. Further investigation into these hypotheses could aid in forming a more complete picture of the distributions and populations of members of genus Latimeria, which in turn could aid in developing conservation strategies, particularly in the case of L. menadoensis.     

Molecules,fossils, and the origin of tetrapods

Summary Since the discovery of the coelacanth, Latimeria chalumnae, more than 50 years ago, paleontologists and comparative morphologists have debated whether coelacanths or lungfishes, two groups of lobe-finned fishes, are the closest living relatives of land vertebrates (Tetrapoda). Previously, Meyer and Wilson (1990) determined partial DNA sequences from two conservative mitochondrial genes and found support for a close relationship of lungfishes to tetrapods. We present additional DNA sequences from the 12S rRNA mitochondria gene for three species of the two lineages of lungfishes that were not represented in the first study: Protopterus annectens and Protopterus aethiopicus from Africa and Neoceratodus forsteri (kindly provided by B. Hedges and L. Maxson) from Australia. This extended data set tends to group the two lepidosirenid lungfish lineages (Lepidosiren and Protopterus) with Neoceratodus as their sister group. All lungfishes seem to be more closely related to tetrapods than the coelacanth is. This result appears to rule out the possibility that the coelacanth lineage gave rise to land vertebrates. The common ancestor of lungfishes and tetrapods might have possessed multiple morphological traits that are shared by lungfishes and tetrapods [Meyer and Wilson (1990) listed 14 such traits]. Those traits that seem to link Latimeria and tetrapods are arguably due to convergent evolution or reversals and not to common descent. In this way, the molecular tree facilitates an evolutionary interpretation of the morphological differences among the living forms. We recommended that the extinct groups of lobe-finned fishes be placed onto the molecular tree that has lungfishes and not the coelacanth more closely related to tetrapods. The placement of fossils would help to further interpret the sequence of morphological events and innovations associated with the origin of tetrapods but appears to be problematic because the quality of fossils is not always high enough, and differences among paleontologists in the interpretation of the fossils have stood in the way of a consensus opinion for the branching order among lobefinned fishes. Marshall and Schultze (1992) criticized the morphological analysis presented by Meyer and Wilson (1990) and suggest that 13 of the 14 morphological traits that support the sister group relationship of lungfishes and tetrapods are not shared derived characters. Here we present further alternative viewpoints to the ones of Marshall and Schultze (1992) from the paleontological literature. We argue that all available information (paleontological, neontological, and molecular data) and rigorous cladistic methodology should be used when relating fossils and extant taxa in a phylogenetic framework. Offprint requests to: Axel Meyer     

Draft genomes and genomic divergence of two Lepidurus tadpole shrimp species (Crustacea,Branchiopoda, Notostraca)

Crustaceans of the order Notostraca (Branchiopoda) are distributed worldwide and are known for the remarkable morphological stasis between their extant and Permian fossil species. Moreover, these crustaceans show relevant ecological traits and a wide range of reproductive strategies. However, genomic studies on notostracans are fairly limited. Here, we present the genome sequences of two notostracan taxa, Lepidurus arcticus and Lepidurus apus lubbocki. Taking advantage of the small genome sizes (~0.11 pg) of these taxa, genomes were sequenced for one individual per species with one run on the Illumina HiSeq X platform. We finally assembled 73.2 Mbp (L. arcticus) and 90.3 Mbp (L. apus lubbocki) long genomes. Assemblies cover up to 84% of the estimated genome size, with a gene completeness >97% for both genomes. In total, 13%–16% of the assembled genomes consist of repeats, and based on read mapping, L. apus lubbocki shows a significantly lower transposable element content than L. arcticus. The analysis of 2,376 orthologous genes indicates an ~7% divergence between the two Lepidurus taxa, with a nucleotide substitution rate significantly lower than that of Daphnia taxa. K_a/K_s analysis suggests purifying selection in both branchiopod lineages, raising the question of whether the low substitution rate of Lepidurus is correlated with morphological conservation or is linked to specific biological traits. Our analysis demonstrates that, in these organisms, it is possible to obtain high‐quality draft genomes from single individuals with a relatively low sequencing effort. This result makes Lepidurus and Notostraca interesting models for genomic studies at taxonomic, ecological and evolutionary levels.     

Genome evolution and biodiversity in teleost fish

Teleost fish, which roughly make up half of the extant vertebrate species, exhibit an amazing level of biodiversity affecting their morphology, ecology and behaviour as well as many other aspects of their biology. This huge variability makes fish extremely attractive for the study of many biological questions, particularly of those related to evolution. New insights gained from different teleost species and sequencing projects have recently revealed several peculiar features of fish genomes that might have played a role in fish evolution and speciation. There is now substantial evidence that a round of tetraploidization/rediploidization has taken place during the early evolution of the ray-finned fish lineage, and that hundreds of duplicate pairs generated by this event have been maintained over hundreds of millions of years of evolution. Differential loss or subfunction partitioning of such gene duplicates might have been involved in the generation of fish variability. In contrast to mammalian genomes, teleost genomes also contain multiple families of active transposable elements, which might have played a role in speciation by affecting hybrid sterility and viability. Finally, the amazing diversity of sex determination systems and the plasticity of sex chromosomes observed in teleost might have been involved in both pre- and postmating reproductive isolation. Comparison of data generated by current and future genome projects as well as complementary studies in other species will allow one to approach the molecular and evolutionary mechanisms underlying genome diversity in fish, and will certainly significantly contribute to our understanding of gene evolution and function in humans and other vertebrates.     

Why coelacanths are not ‘living fossils’

A series of recent studies on extant coelacanths has emphasised the slow rate of molecular and morphological evolution in these species. These studies were based on the assumption that a coelacanth is a ‘living fossil’ that has shown little morphological change since the Devonian, and they proposed a causal link between low molecular evolutionary rate and morphological stasis. Here, we have examined the available molecular and morphological data and show that: (i) low intra‐specific molecular diversity does not imply low mutation rate, (ii) studies not showing low substitution rates in coelacanth are often neglected, (iii) the morphological stability of coelacanths is not supported by paleontological evidence. We recall that intra‐species levels of molecular diversity, inter‐species genome divergence rates and morphological divergence rates are under different constraints and they are not necessarily correlated. Finally, we emphasise that concepts such as ‘living fossil’, ‘basal lineage’, or ‘primitive extant species’ do not make sense from a tree‐thinking perspective. Editor's suggested further reading in BioEssays Tree thinking for all biology: the problem with reading phylogenies as ladders of progress Abstract     

Gene Structure and Amino Acid Sequence of Latimeria chalumnae (Coelacanth) Myelin DM20: Phylogenetic Relation of the Fish

The structure of Latimeria chalumnae (coelacanth) proteolipid protein/DM20 gene excluding exon 1 was determined, and the amino acid sequence of Latimeria DM20 corresponding to exons 2–7 was deduced. The nucleotide sequence of exon 3 suggests that only DM20 isoform is expressed in Latimeria. The structure of proteolipid protein/DM20 gene is well preserved among human, dog, mouse, and Latimeria. Southern blot analysis indicates that Latimeria DM20 gene is a single-copy gene. When the amino acid sequences of DM20 were compared among various species, Latimeria was more similar to tetrapods than other fishes including lungfish, confirming the previous finding by immunoreactivity (Waehneldt and Malotka 1989 J. Neurochem. 52:1941–1943). However, when phylogenetic trees were constructed from the DM20 sequences, lungfish was clearly the closest to tetrapods. Latimeria was situated outside of lungfish by the maximum likelihood method. The apparent similarity of Latimeria DM20 to tetrapod proteolipid protein/DM20 is explained by the slow amino acid substitution rate of Latimeria DM20.     

Analysis of collinear regions of Oryza AA and CC genomes

Comparative analyses of genome structure and sequence of closely related species have yielded insights into the evolution and function of plant genomes. A total of 103,844 BAC end sequences delegated -73.8 Mb of O. officinalis that belongs to the CC genome type of the rice genus Oryza were obtained and compared with the genome sequences office cultivar, O. sativa ssp.japonica cv. Nipponbare. We found that more than 45% of O. officinalis genome consists of repeat sequences, which is higher than that of Nipponbare cultivar. To further investigate the evolutionary divergence of AA and CC genomes, two BAC-contigs of O. officinalis were compared with the collinear genomic regions of Nipponbare. Of 57 genes predicted in the AA genome orthologous regions, 39 had orthologs in the regions of the CC genome. Alignment of the orthologous regions indicated that the CC genome has undergone expansion in both genic and intergenic regions through primarily retroelement insertion. Particularly, the density of RNA transposable elements was 17.95% and 1.78% in O. officinalis and O. sativa, respectively. This explains why the orthologous region is about 100 kb longer in the CC genome in comparison to the AA genome.     

The Complete DNA Sequence of the Mitochondrial Genome of a ``living Fossil,'''' the Coelacanth (Latimeria Chalumnae)

The complete nucleotide sequence of the 16,407-bp mitochondrial genome of the coelacanth (Latimeria chalumnae) was determined. The coelacanth mitochondrial genome order is identical to the consensus vertebrate gene order which is also found in all ray-finned fishes, the lungfish, and most tetrapods. Base composition and codon usage also conform to typical vertebrate patterns. The entire mitochondrial genome was PCR-amplified with 24 sets of primers that are expected to amplify homologous regions in other related vertebrate species. Analyses of the control region of the coelacanth mitochondrial genome revealed the existence of four 22-bp tandem repeats close to its 3' end. The phylogenetic analyses of a large data set combining genes coding for rRNAs, tRNAs, and proteins (16,140 characters) confirmed the phylogenetic position of the coelacanth as a lobe-finned fish; it is more closely related to tetrapods than to ray-finned fishes. However, different phylogenetic methods applied to this largest available molecular data set were unable to resolve unambiguously the relationship of the coelacanth to the two other groups of extant lobe-finned fishes, the lungfishes and the tetrapods. Maximum parsimony favored a lungfish/coelacanth or a lungfish/tetrapod sistergroup relationship depending on which transversion:transition weighting is assumed. Neighbor-joining and maximum likelihood supported a lungfish/tetrapod sistergroup relationship.     

Evolution of Genome size in Hawaiian Endemic Genus Schiedea (Caryophyllaceae)

Genome sizes may vary by orders of magnitude among relatively closely related species. Gene and genome duplications and movements of transposable elements (TEs) can quickly inflate genome sizes. Whole genome duplications (polyploidization) may be an important source of evolutionary innovation and many fast evolving island genera are known or assumed to be polyploids. Our main aim is to shed light on the question of how genome size evolved within a rapidly diversifying island lineage. We report the estimates of DNA content for 27 species of the Hawaiian endemic plant genus Schiedea and its widespread sister genus Honckenya (Caryophyllaceae: Alsinoideae). Unexpectedly, genomes of Schiedea species appeared to be relatively compact (1.41 to 3.74 pg/cell), compared to Honckenya (8.57 to 10.66 pg/cell). Interestingly, Schiedea species from younger islands tended to have larger genomes than species from older islands, which may be explained by activation of TE transpositions in small populations after colonization events that resulted in the formation of new species on younger islands. To test whether the Schiedea genome has undergone recent polyploidization events we measured divergence between 62 pairs of paralogous genes in S. globosa. The distribution of divergence values was unimodal with a mode of 7%, supporting a single polyploidization event. Dating this event using Schiedea/Honckenya divergence (2%) and Schiedea/Silene divergence (11%) we estimate that it might have occurred in the ancestor of the genera Schiedea and Honckenya, but after the split between the subfamilies Alsinoideae and Silenoideae.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Ciliated Sensory Cells in the Rostral Organ of the Coelacanth Latimeria chalumnae (Smith 1939)

The rostral organ of the coelacanth Latimeria chalumnae contains a sensory epithelium with luminated ciliated sensory cells, which very much resemble electroreceptors from other non-teleost fish. This, as well as the innervation and a dorsal nucleus, indicate that this interesting species is electroreceptive.     

Evolutionary sequence analysis of complete eukaryote genomes

Background  

Gene duplication and gene loss during the evolution of eukaryotes have hindered attempts to estimate phylogenies and divergence times of species. Although current methods that identify clusters of orthologous genes in complete genomes have helped to investigate gene function and gene content, they have not been optimized for evolutionary sequence analyses requiring strict orthology and complete gene matrices. Here we adopt a relatively simple and fast genome comparison approach designed to assemble orthologs for evolutionary analysis. Our approach identifies single-copy genes representing only species divergences (panorthologs) in order to minimize potential errors caused by gene duplication. We apply this approach to complete sets of proteins from published eukaryote genomes specifically for phylogeny and time estimation.     

Characterization of a Tc1-like transposable element in zebrafish (Danio rerio)

We have characterized Tdr1, a family of Tc1-like transposable elements found in the genome of zebrafish (Danio rerio). The copy number and distribution of the sequence in the zebrafish genome have been determined, and by these criteria Tdr1 can be classified as a moderately repetitive, interspersed element. Examination of the sequences and structures of several copies of Tdr1 revealed that a particular deletion derivative, 1250 by long, of the transposon has been amplified to become the dominant form of Tdr1. The deletion in these elements encompasses sequences encoding the N-terminal portion of the putative Tdr1 transposase. Sequences corresponding to the deleted region were also detected, and thus allowed prediction of the nucleotide sequence of a hypothetical full-length element. Well conserved segments of Tc1-like transposons were found in the flanking regions of known fish genes, suggesting that these elements have a long evolutionary history in piscine genomes. Tdr1 elements have long, 208 by inverted repeats, with a short DNA motif repeated four times at the termini of the inverted repeats. Although different from that of the prototype C. elegans transposon Tc1, this inverted repeat structure is shared by transposable elements from salmonid fish species and two Drosophila species. We propose that these transposons form a subgroup within the Tc1-like family. Comparison of Tc1-like transposons supports the hypothesis that the transposase genes and their flanking sequences have been shaped by independent evolutionary constraints. Although Tc1-like sequences are present in the genomes of several strains of zebrafish and in salmonid fishes, these sequences are not conserved in the genus Danio, thus raising the possibility that these elements can be exploited for gene tagging and genome mapping.     

Current views of the origin and diversification of tetrapods

In several recent decades, a wealth of evidence was obtained (including new fossils and new nucleotide sequences) that allows a revision of the current evolutionary theories. The primary structure was determined for orthologous regions from a variety of genomes. Complete genome sequences were established for some organisms. The review considers the modern concepts of the origin and diversification of tetrapods and of the formation of modern vertebrate classes.     

The nuclear DNA content and chromatin ultrastructure of the coelacanth Latimeria chalumnae

By measurement of 731 erythrocytes by Feulgen cytophotometry, the nuclear DNA content of the coelacanth Latimeria chalumnae Smith is determined to be 7.22 picograms (pg). This value is high among fishes but is closely comparable to that of man and most other mammals. The average mass of erythrocyte nuclear chromatin, measured by quantitative electron microscopy, is 15.2 pg. This chromatin is in the form of fibers having a mean diameter of 202 Å. The average weight of the chromatin fiber is 6.75 × 10^?16 g/μm. Thus, the nucleus contains 22 500 μm of chromatin fiber. Dividing the nuclear DNA content of Latimeria by the known mass of the DNA double helix (3.26 × 10^?18 g/μm) gives a total length of 2 215 000 μm of DNA double helix. In comparing these two measurements of structural length, it is found that 98.4 lengths of double helix are packed into one length of chromatin fiber. This packing ratio is over three times greater than that of human G1 lymphocytes. The difference may be attributable to the difference between the two tissues and thus reflect a functional distinction, or it may be due to the difference between the two species and reflect an evolutionary distinction.     

Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

Recent sequencing of the Brassica rapa and Brassica oleracea genomes revealed extremely contrasting genomic features such as the abundance and distribution of transposable elements between the two genomes. However, whether and how these structural differentiations may have influenced the evolutionary rates of the two genomes since their split from a common ancestor are unknown. Here, we investigated and compared the rates of nucleotide substitution between two long terminal repeats (LTRs) of individual orthologous LTR‐retrotransposons, the rates of synonymous and non‐synonymous substitution among triplicated genes retained in both genomes from a shared whole genome triplication event, and the rates of genetic recombination estimated/deduced by the comparison of physical and genetic distances along chromosomes and ratios of solo LTRs to intact elements. Overall, LTR sequences and genic sequences showed more rapid nucleotide substitution in B. rapa than in B. oleracea. Synonymous substitution of triplicated genes retained from a shared whole genome triplication was detected at higher rates in B. rapa than in B. oleracea. Interestingly, non‐synonymous substitution was observed at lower rates in the former than in the latter, indicating shifted densities of purifying selection between the two genomes. In addition to evolutionary asymmetry, orthologous genes differentially regulated and/or disrupted by transposable elements between the two genomes were also characterized. Our analyses suggest that local genomic and epigenomic features, such as recombination rates and chromatin dynamics reshaped by independent proliferation of transposable elements and elimination between the two genomes, are perhaps partially the causes and partially the outcomes of the observed inter‐specific asymmetric evolution.