Search

The cichlids of Lake Victoria are a textbook example of adaptive radiation, as >500 endemic species arose in just 14,600 years. The degree of genetic differentiation among species is very low due to the short period of time after the radiation, which allows us to ascertain highly differentiated genes that are strong candidates for driving speciation and adaptation. Previous studies have revealed the critical contribution of vision to speciation by showing the existence of highly differentiated alleles in the visual opsin gene among species with different habitat depths. In contrast, the processes of species-specific adaptation to different ecological backgrounds remain to be investigated. Here, we used genome-wide comparative analyses of three species of Lake Victoria cichlids that inhabit different environments-Haplochromis chilotes, H. sauvagei, and Lithochromis rufus-to elucidate the processes of adaptation by estimating population history and by searching for candidate genes that contribute to adaptation. The patterns of changes in population size were quite distinct among the species according to their habitats. We identified many novel adaptive candidate genes, some of which had surprisingly long divergent haplotypes between species, thus showing the footprint of selective sweep events. Molecular phylogenetic analyses revealed that a large fraction of the allelic diversity among Lake Victoria cichlids was derived from standing genetic variation that originated before the adaptive radiation. Our analyses uncovered the processes of species-specific adaptation of Lake Victoria cichlids and the complexity of the genomic substrate that facilitated this adaptation., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Reproductive isolation that prevents interspecific hybridization between closely related coexisting species maintains sympatric species diversity. One of the reproductive isolations is mate choice based on color signals (breeding color perceived by color vision). This is well known in several animal taxa, yet little is known about its genetic and molecular mechanism. Lake Victoria cichlid fishes are thought to be an example of sympatric species diversity. In the species inhabiting different light environments in rocky shore, speciation by sensory drive through color signals has been proposed by analyses of the long wavelength-sensitive (LWS) opsin gene and the male nuptial coloration. However, the genetic and molecular mechanism of how diversity of sympatric species occurring in the same habitat is maintained remains unknown. To address this issue, we determined nucleotide sequences of eight opsins of six sympatric species collected from a sandy-muddy shore--an ideal model system for studying sympatric species. Among eight opsins, the LWS and RH1 alleles were diversified and one particular allele is dominant or fixed in each species, and we propose that this is due to natural selection. The functions of their LWS alleles were also diversified as shown by absorption measurements of reconstituted visual pigments. To analyze the relationship between nuptial coloration and the absorption of LWS pigments, we systematically evaluated and defined nuptial coloration. We showed that the coloration was species specific with respect to hue and significantly differentiated by the index values of hue (dominant wavelength: λ(d)). The λ(d) value of the male nuptial coloration correlated with the absorption of LWS pigments from all the species, suggesting that reproductive isolation through mate choice using color signals may prevent sympatric interspecific hybridization, thereby maintaining the species diversity in sympatric species in Lake Victoria.

Reverse evolution is a widespread phenomenon in biology, but the genetic mechanism for the reversal of a genetic change for adaptation to the ancestral state is not known. Here, we report the first case of complete reverse evolution of two amino acids, serine and alanine, at a single position in RH1 opsin pigment for adaptation to water depth. We determined RH1 sequences of cichlid fishes from four tribes of Lake Tanganyika with different habitat depths. Most of the species were divided into two types: RH1 with 292A for species in shallow water or 292S for species in deep water. Both types were adapted to their ambient light environments as indicated by the absorption spectra of the RH1 pigments. Based on the RH1 locus tree and ecological data, we inferred the ancestral amino acids at position 292 and the distribution of the depth ranges (shallow or deep) of ancestral species of each tribe. According to these estimates, we identified two distinct parallel adaptive evolutions: The replacement A292S occurred at least four times for adaptation from shallow to deep water, and the opposite replacement S292A occurred three times for adaptation from deep to shallow water. The latter parallelism represents the complete reverse evolution from the derived to the ancestral state, following back adaptive mutation with reversal of the RH1 pigment function accompanied by reversal of the species habitat shift.

Rhodopsin is a photoreceptive protein present in vertebrate rod photoreceptor cells, which are responsible for scotopic vision. Recent molecular studies have shown that several aquatic vertebrate species have independently acquired rhodopsin containing Asp83Asn, Glu122Gln, and Ala292Ser substitutions, causing a blue shift in the rhodopsin absorption spectra for adaptation to the blue-green photic environment in deep water. Here, we provide new evidence for the evolutionary and functional relevance of the Asp83Asn substitution. Spectroscopic and kinetic analyses of rhodopsins in six cichlid fishes from the East African Great Lakes using charge-coupled device spectrophotometer revealed that the Asp83Asn substitution accelerated the formation of meta-II, a rhodopsin intermediate crucial for activation of the G-protein transducin. Because rapid formation of meta-II likely results in effective transduction of photic signals, it is reasonable to assume that deep-water cichlid species have acquired rhodopsin containing Asn83 to adapt to dim lighting. Remarkably, rhodopsin containing Asn83 has been identified in terrestrial vertebrates such as bats, and these rhodopsin variants also exhibit accelerated meta-II formation. Our results indicated that the Asp83Asn substitution observed in a variety of animal species was acquired independently in many different lineages during vertebrate evolution for adaptation to dimly lit environments.

Human long interspersed nuclear element-1 (L1) occupies one-sixth of our genome and has contributed to genome evolution in various ways. Approximately 10% of human L1 copies are composed of two L1 segments; the 5' segment and 3' segment are in head-to-head (i.e., 5'-inverted) orientation. Besides mediating their own retrotransposition, L1 has the ability to mobilize mRNA "in trans," and the number of retrotransposed mRNA sequences (retrocopies) is estimated to be >6,000. In this study, we identified 48 human-specific retrocopies and 95 chimpanzee-specific retrocopies by comparing the human and chimpanzee genomes. Among these retrocopies, 12 were 5'-inverted. The characteristics of these 5'-inverted retrocopies were similar to those of 5'-inverted L1 copies, indicating that the 5' inversion is generated by the same mechanism. With these findings, we examined the possibility that 5' inversion of the retrocopy generates a new gene that codes for a peptide with a different N terminus. We identified several potential 5'-inverted retrogenes, including those of thymopoietin beta (TMPO) and eukaryotic translation initiation factor 3 subunit 5 (EIF3F). The most interesting candidate was the 5'-inverted retrocopy of small nuclear ribonucleoprotein polypeptide N (SNRPN). This retrocopy was transcribed in the reverse orientation in several organs, had multiple transcript variants, and encoded a protein containing a peptide fragment derived from the N-terminal portion of SNRPN. Our results suggest that mRNA retrotransposition coupled with 5' inversion may be a mechanism to generate new genes distinct from parental genes.

Autonomous non-long terminal repeat retrotransposons (NLRs) are ubiquitous mobile genetic elements that insert their DNA copies at new locations by retrotransposition. In vertebrates, there are 4 NLR clades, L1, L2, CR1, and RTE, which diverged in the Precambrian era. It has been demonstrated that retrotransposition of L1 and L2 members proceeds via coordinated reactions of targeted DNA cleavage and reverse transcription catalyzed by the NLR-encoded proteins, which are followed by the joining of the 5' (upstream) junction. However, the study on the mobility pathways for vertebrate NLRs is so far limited to L1 and L2. In this report, using target analysis of nested transposons for genomic copies, we studied retrotransposition pathways for a variety of vertebrate NLRs, including those of the L1, L2, CR1, and RTE clades in the human, cow, opossum, chicken, and zebrafish genomes. Thus, this study constitutes the first comprehensive analysis of NLR retrotransposition products in vertebrates. Our data revealed that these elements share similar mechanisms for the cleavages of the 2 target DNA strands and for the initiation of reverse transcription. Possible endonuclease-independent insertions were also identified. Overall, our results suggest the existence of multiple retrotransposition pathways that are conserved among the diverse NLR clades in various vertebrate hosts.

L1 is the most proliferative autonomous retroelement that comprises about 20% of mammalian genomes. Why L1s have proliferated so extensively in mammalian genomes is an important yet unsolved question. L1 copies are amplified via retrotransposition, in which the DNA cleavage specificity by the L1-encoded endonuclease (EN) primarily dictates sites of insertion. Whereas mammalian L1s show target preference for 5'-TTAAAA-3', other L1-like elements exhibit various degrees of target specificity. To gain insights on diversification of the EN specificity during L1 evolution, ENs of zebrafish L1 elements were analyzed here. We revealed that they form 3 discrete clades, M, F, and Tx1, which is in stark contrast to a single L1 clade in mammalian species. Interestingly, zebrafish clade M elements cluster as a sister group of mammalian L1s and show target-site preference for 5'-TTAAAA-3'. In contrast, elements of the clade F, the immediate outgroup of the clade M, show little specificity. We identified certain clade-specific amino acid residues in EN, many of which are located in the cleft that recognizes the substrate, suggesting that these amino acid alterations have generated 2 types of ENs with different substrate specificities. The distribution pattern of the 3 clades suggests a possibility that the acquisition of target specificity by the L1 ENs improved the L1 fitness under the circumstances in mammalian hosts.

Short interspersed elements (SINEs) constitute a group of retroposons propagating in the genome via a mechanism of reverse transcription, in which they depend on the enzymatic machinery of long retroposons (LINEs). Over 70 SINE families have been described to date from the genomes of various eukaryotes. Here, we characterize two novel SINEs from salmons (Actinopterygii: Salmonoidei). The first family, termed SlmI, was shown to be widespread among all genera of the suborder. These SINEs have a tRNA(Leu)-related promoter region at their 5'-end, a unique central conserved domain with a subfamily-specific region, and an end with RSg-1-LINE-derived 3'-terminus preceding the A/T-rich tail. The same LINE-related segment is also shared by two other salmonid SINEs: HpaI and OS-SINE1. The structural peculiarities and overall sequence identity of the SlmI 3'-terminus suggest that it has been acquired from HpaI SINEs but not directly from the partner LINE. This region plays a crucial role in the process of retrotransposition of short interspersed elements, and the case of its SINE-to-SINE transmission is the first recorded to date. Possible scenarios and potential evolutionary implications of the observed interaction between short retroposons are discussed. Apart from the above, we found a copy of the SlmI SINE in the GenBank entry for the blood fluke, Schistosoma japonicum (Trematoda: Strigeiformes) -- a trematode causing one of the most important human helminth infections, with its genome known to host other groups of salmonoid retroposons. In the present article, we suggest our views with regard to possible ways in which such an intensive horizontal transfer of salmonoid retroposons to the schistosomal genome occurs. The second novel SINE family, termed SlmII, originates from one of the SlmI subfamilies, with which it shares the same tRNA-related region, central domain, and a part of RSg-1-derived segment, but has a different 3'-tail of unidentified origin. Its distribution among salmonids validates Parahucho (Japanese huchen) as a distinct monotypic genus.

Baleen whales (suborder Mysticeti) comprise 11 extant species that are classified into four families. Although several phylogenetic hypotheses about these taxa have been proposed, their phylogenetic relationships remain confused. We addressed this problem using short interspersed repetitive element (SINE) insertion data, which now are regarded as almost ideal shared, derived characters at the molecular level. We reconstructed the phylogenetic relationships of baleen whales by characterizing 36 informative SINE loci. One of the intriguing conclusions is that balaenopterids and eschrichtiids radiated very rapidly during a very short evolutionary period. During this period, speciation occurred in balaenopterids and eschrichtiids while newly inserted SINE loci remains polymorphic. Later on, these SINEs were sorted incompletely into each lineage. Thus, there are now inconsistencies among species regarding the presence or absence of a given SINE. This is in sharp contrast to the phylogeny of toothed whales, for which no SINE inconsistencies have been found. Furthermore, we found monophyletic groupings between humpback and fin whales as well as between (sei+Bryde's) whales and blue whales, both of which have not previously been recognized. The comprehensive SINE insertion data, together with the mitochondrial DNA phylogeny that was recently completed (Sasaki, T., M. Nikaido, H. Healy et al. 2005. Mitochondrial phylogenetics and evolution of mysticete whales. Syst. Biol. 56:77-90; Rychel, A. L., T. W. Reeder, and A. Berta. 2004. Phylogeny of mysticete whales based on mitochondrial and nuclear data. Mol. Phylogenet. Evol. 32:892-901), provide a nearly complete picture of the evolutionary history of baleen whales.

Recent comprehensive studies of DNA sequences support the monophyly of Afrotheria, comprising elephants, sirenians (dugongs and manatees), hyraxes, tenrecs, golden moles, aardvarks, and elephant shrews, as well as that of Paenungulata, comprising elephants, sirenians, and hyraxes. However, phylogenetic relationships among paenungulates, as well as among nonpaenungulates, have remained ambiguous. Here we applied an extensive retroposon analysis to these problems to support the monophyly of aardvarks, tenrecs, and golden moles, with elephant shrews as their sister group. Regarding phylogenetic relationships in Paenungulata, we could characterize only one informative locus, although we could isolate many insertions specific to each of three lineages, namely, Proboscidea, Sirenia, and Hyracoidea. These data prompted us to reexamine phylogenetic relationships among Paenungulata using 19 nuclear gene sequences resulting in three different analyses, namely, short interspersed element (SINE) insertions, nuclear sequence analyses, and morphological cladistics, supporting different respective phylogenies. We concluded that these three lineages diverged very rapidly in a very short evolutionary period, with the consequence that ancestral polymorphism present in the last common ancestor of Paenungulata results in such incongruence. Our results suggest the rapid fixation of many large-scale morphological synapomorphies for Tethytheria; implications of this in relation to the morphological evolution in Paenungulata are discussed.

Lake Victoria cichlid fishes are excellent examples of explosive adaptive radiation. Although Lake Victoria cichlids are believed to have arisen during a short period (approximately 14,000 years), they have various species-specific phenotypes. One important phenotype that distinguishes each species is the shape of the jaw, which has diverged to adapt to the wide variety of trophic habitats present in the lake. Here we demonstrate a new approach to investigate the diversification of cichlid jaw morphology at the genetic level by examining differentially expressed genes. We used a DNA chip to compare gene expression levels between closely related cichlid fishes. This analysis indicated that the expression of some genes differed in the larvae of two cichlid species. One such clone encodes a new astacin family metalloproteinase. The expression level of the isolated gene, named cimp1, was analyzed in more detail by real-time quantitative reverse transcription-polymerase chain reaction. A significant difference in cimp1 expression was observed between two Haplochromis cichlid species during development. Using in situ hybridization, we found that this gene is expressed only in head and gill epithelia. Biochemical analysis showed that cichlid metalloproteinase 1 (CiMP1) has proteolytic activity, a common attribute of all astacin family proteins. Because some astacin family proteins contribute to morphogenesis in animals, CiMP1 is expected to participate in species-specific head morphogenesis in cichlids. This is the first study to demonstrate that differentially expressed genes among cichlids can be identified using a DNA chip.

Long interspersed elements (LINEs) and short interspersed elements (SINEs) are retrotransposons. These elements can mobilize by the "copy-and-paste" mechanism, in which their own RNA is reverse-transcribed into complementary DNA (cDNA). LINEs and SINEs not only are components of eukaryotic genomes but also drivers of genomic evolution. Thus, studies of the amplification mechanism of LINEs and SINEs are important for understanding eukaryotic genome evolution. Here we report the characterization of one LINE family (UnaL2) and two SINE families (UnaSINE1 and UnaSINE2) from the eel (Anguilla japonica) genome. UnaL2 is approximately 3.6 kilobases (kb) and encodes only one open reading frame (ORF). UnaL2 belongs to the stringent type--thought to be a major group of LINEs--and can mobilize in HeLa cells. We also show that UnaL2 and the two UnaSINEs have similar 3' tails, and that both UnaSINE1 and UnaSINE2 can be mobilized by UnaL2 in HeLa cells. These elements are thus useful for delineating the amplification mechanism of stringent type LINEs as well as that of SINEs.

Although turtles (order Testudines) constitute one of the major reptile groups, their phylogenetic relationships remain largely unresolved. Hence, we attempted to elucidate their phylogeny using the SINE (short interspersed repetitive element) method, in which the sharing of a SINE at orthologous loci is indicative of synapomorphy. First, a detailed characterization of the tortoise polIII/SINE was conducted using 10 species from eight families of hidden-necked turtles (suborder Cryptodira). Our analysis of 382 SINE sequences newly isolated in the present study revealed two subgroups, namely Cry I and Cry II, which were distinguishable according to diagnostic nucleotides in the 3' region. Furthermore, four (IA-ID) and five (IIA-IIE) different SINE types were identified within Cry I and Cry II subgroups, respectively, based on features of insertions/deletions located in the middle of the SINE sequences. The relative frequency of occurrence of the subgroups and the types of SINEs in this family were highly variable among different lineages of turtles, suggesting active differential retroposition in each lineage. Further application of the SINE method using the most retrotranspositionally active types, namely IB and IC, challenged the established phylogenetic relationships of Bataguridae and its related families. The data for 11 orthologous loci demonstrated a close relationship between Bataguridae and Testudinidae, as well as the presence of the three clades within Bataguridae. Although the SINE method has been used to infer the phylogenies of a number of vertebrate groups, it has never been applied to reptiles. The present study represents the first application of this method to a phylogenetic analysis of this class of vertebrates, and it provides detailed information on the SINE subgroups and types. This information may be applied to the phylogenetic resolution of relevant turtle lineages.

Short interspersed nuclear elements (SINEs) are dispersed repetitive DNA sequences that are major components of all mammalian genomes. They have been described in almost all lineages of Euarchontoglires (rodents, rabbits, primates, flying lemurs, and tree shrews), except in flying lemurs. Most SINE family members are composed of three distinct regions: a 5' tRNA-related region, a tRNA-unrelated region, and a short tandem repeat at the 3' end that is AT-rich. The newly discovered SINE family in Cynocephalus deviates from this common structure. All 30 SINE loci analyzed in this family lack a tRNA-unrelated region and are composed exclusively of tRNA-related elements. Therefore, this novel SINE structure, described for the first time in mammalian genomes, was designated as t-SINE. The t-SINE family exhibits a high copy number and is specific to flying lemurs. Three major t-SINE subfamilies could be distinguished on the basis of characteristic nucleotides, deletions, insertions, and duplications. These sequence-specific characteristics within subfamilies and sub-subfamilies reveal that they are derived copies of distinct progenitors. We present evolutionary relationships between subfamilies and compare relationships between the subfamilies and the isoleucine tRNA gene. t-SINE amplification occurred through multiple sources and is supposedly mobilized via the L1-encoded reverse transcriptase-dependent retrotranspositional mechanism in trans.

Cichlid fishes of the east African Great Lakes represent a paradigm of adaptive radiation. We conducted a phylogenetic analysis of cichlids including pan-African and west African species by using insertion patterns of short interspersed elements (SINEs) at orthologous loci. The monophyly of the east African cichlids was consistently supported by seven independent insertions of SINE sequences that are uniquely shared by these species. In addition, data from four other loci indicated that the genera Tilapia (pan-African) and Steatocranus (west African) are the closest relatives to east African cichlids. However, relationships among Tilapia, Steatocranus, and the east African clade were ambiguous because of incongruencies among topologies suggested by insertion patterns of SINEs at six other loci. One plausible explanation for this phenomenon is incomplete lineage sorting of alleles containing or missing a SINE insertion at these loci during ancestral speciation. Such incomplete sorting may have taken place earlier than 14 MYA, followed by random and stochastic fixation of the alleles in subsequent lineages. These observations prompted us to consider the possibility that cichlid speciation occurred at an accelerated rate during this period when the African Great Lakes did not exist. The SINE method could be useful for detecting ancient exclusive speciation events that tend to remain hidden during conventional sequence analyses because of accumulated point mutations.

Afrotheria is a newly recognized taxon comprising elephants, hyraxes, sea cows, aardvarks, golden moles, tenrecs, and elephant shrews, each of which originated in Africa. Although some members of this taxon were once classified into distantly related groups, recent molecular studies have demonstrated their close relationships. It was suggested that this group emerged as a result of physical isolation of the African continent during the successive breakup events of Gondowanaland. In this study, a novel family of SINEs, designated AfroSINEs, was isolated and characterized from the genomes of afrotherians. This SINE family is distributed exclusively among the afrotherian species, confirming their monophyletic relationships. Furthermore, a distinct subfamily, which shares a deletion in the middle region of the SINE, was identified. The distribution of this subfamily is apparently restricted to the genomes of hyraxes, elephants, and sea cows, suggesting monophyly of these three groups, which was previously proposed as Paenungulata. We characterized the structures of the AfroSINEs from all afrotherian representatives by PCR, and we discuss how they were generated as well as the phylogenetic relationships of their host species.

Until recently, phylogenetic analyses have been routinely based on homologous sequences of a single gene. Given the vast number of gene sequences now available, phylogenetic studies are now based on the analysis of multiple genes. Thus, it has become necessary to devise statistical methods to combine multiple molecular data sets. Here, we compare several models for combining different genes for the purpose of evaluating the likelihood of tree topologies. Three methods of branch length estimation were studied: assuming all genes have the same branch lengths (concatenate model), assuming that branch lengths are proportional among genes (proportional model), or assuming that each gene has a separate set of branch lengths (separate model). We also compared three models of among-site rate variation: the homogenous model, a model that assumes one gamma parameter for all genes, and a model that assumes one gamma parameter for each gene. On the basis of two nuclear and one mitochondrial amino acid data sets, our results suggest that, depending on the data set chosen, either the separate model or the proportional model represents the most appropriate method for branch length analysis. For all the data sets examined, one gamma parameter for each gene represents the best model for among-site rate variation. Using these models we analyzed alternative mammalian tree topologies, and we describe the effect of the assumed model on the maximum likelihood tree. We show that the choice of the model has an impact on the best phylogeny obtained.

We characterized two novel 7SL RNA-derived short interspersed nuclear element (SINE) families (Tu types I and II) and a novel tRNA-derived SINE family (Tu type III) from the tree shrew (Tupaia belangeri). Tu type I contains a monomer unit of a 7SL RNA-derived Alu-like sequence and a tRNA-derived region that includes internal RNA polymerase III promoters. Tu type II has a similar hybrid structure, although the monomer unit of the 7SL RNA-derived sequence is replaced by a dimer. Along with the primate Alu, the galago Alu type II, and the rodent B1, these two families represent the fourth and fifth 7SL RNA-derived SINE families to be identified. Furthermore, comparison of the Alu domains of Tu types I and II with those of other 7SL RNA-derived SINEs reveals that the nucleotides responsible for stabilization of the Alu domain have been conserved during evolution, providing the possibility that these conserved nucleotides play an indispensable role in retropositional activity. Evolutionary relationships among these 7SL RNA-derived SINE families, as well as phylogenetic relationships of their host species, are discussed.

The African cichlid (AFC) family of short interspersed elements (SINEs) is found in the genomes of cichlid fish. The alignment of the sequences of 70 members of this family, isolated from such fish in Africa, revealed the presence of correlated changes in specific nucleotides (diagnostic nucleotides) that allowed us to categorize the various members into six subfamilies, which were designated Af1 through Af6. Dividing the SINE consensus sequence into a 5'-head and 3'-tail region, these subfamilies were defined by various combinations of four types of head region (A-D) and three types of tail region [X, Y, and (YX)], with each region of each type including unique diagnostic nucleotides. The observed structures of the subfamilies Af1 through Af6 were AX, AY, CY, A(YX), BY, and DX, respectively. The formation of such structures might have involved the shuffling of head or tail regions among preexisting and existing (or both) subfamilies of the AFC family (and, probably, even another SINE family or a pseudogene for a tRNA in the case of the Af6 subfamily) by recombination at the so-called core region during the course of evolution. By plotting the timing of the retroposition of individual members of each subfamily on a phylogenetic tree of AFCs, we found that the Af3 and Af6 subfamilies became active only recently in the evolutionary history of these fish. The integrity of the 3'-tails of SINEs, which are, apparently, recognized by reverse transcriptase, has been reported to be indispensable for retention of retropositional activity. Therefore, we postulate that recombination might have been involved in the apparent recent activation of the retroposition of the Af3 and Af6 subfamilies via introduction of active tails (types Y and X, respectively) into potential ancestral sequences that might have had inactive tails. If this hypothesis is correct, shuffling of tail regions among subfamilies by recombination at the core region might have played a role in the recycling of dead copies of AFC SINEs.

Lake Tanganyika harbors numerous endemic species of extremely diverse cichlid fish that have been classified into 12 major taxonomic groups known as tribes. Analysis of short interspersed element (SINE) insertion data has been acknowledged to be a powerful tool for the elucidation of phylogenetic relationships, and we applied this method in an attempt to clarify such relationships among these cichlids. We studied insertion patterns of 38 SINEs in total, 24 of which supported the monophyly of three clades. The other 14 loci revealed extensive incongruence in terms of the patterns of SINE insertions. These incongruencies most likely stem from a period of adaptive radiation. One possible explanation for this phenomenon is the extensive incomplete lineage sorting of alleles for the presence or absence of a SINE during successive speciation events which took place about 5-10 MYA. The present study is the first to report the successful application of the SINE method in demonstrating the existence of such possible "ancient" incomplete lineage sorting. We discuss the possibility that it might potentially be very difficult to resolve the species phylogeny of a group that radiated explosively, even by resolving the genealogies of more than 10 nuclear loci, as a consequence of incomplete lineage sorting during speciation.

Short interspersed nuclear elements (SINEs) have been used to generate unambiguous phylogenetic topologies relating eukaryotic taxa. The irreversible nature of SINE retroposition is supported by a large body of comparative genome data and is a fundamental assumption inherent in the value of this qualitative method of inference. Here, we assess the key assumption of unidirectional SINE insertion by comparing the SINE insertion-derived topology and the phylogenetic tree based on seven independent loci of five taxa in the order Cetartiodactyla (Cetacea + Artiodactyla). The data sets and analyses were largely independent, but the loci were, by definition, linked, and thus their consistency supported an irreversible pattern of SINE retroposition. Moreover, our analyses of the flanking sequences provided estimates of divergence times among cetartiodactyl lineages unavailable from SINE insertion analysis alone. Unexpected rate heterogeneity among sites of SINE-flanking sequences and other noncoding DNA sequences were observed. Sequence simulations suggest that this rate heterogeneity may be an artifact resulting from the inaccuracies of the substitution model used.

Some previously unidentified short interspersed repetitive elements (SINEs) and long interspersed repetitive element (LINEs) were isolated from various higher elasmobranchs (sharks, skates, and rays) and characterized. These SINEs, members of the HE1 SINE family, were tRNA-derived and were widespread in higher elasmobranches. The 3'-tail region of this SINE family was strongly conserved among elasmobranchs. The LINEs, members of the HER1 LINE family, encoded an amino acid sequence similar to that encoded by the chicken CR1 LINE family, and they contained a strongly conserved 3'-tail region in the 3' untranslated region. This tail region of the HER1 LINE family was almost identical to that of the HE1 SINE family. Thus, the HE1 SINE family and the HER1 LINE family provide a clear example of a pair of SINEs and LINEs that share the same tail region. Conservation of the secondary structures of the tail regions, as well as of the nucleotide sequences, between the HE1 SINE family and HER1 LINE family during evolution suggests that SINEs utilize the enzymatic machinery for retroposition of LINEs through the recognition of higher-order structures of the conserved 3'-tail region. A discussion is presented of the parasitism of SINEs on LINEs during the evolution of these retroposons.

Several novel (sub)families of SINEs were isolated from the genomes of cetaceans and artiodactyls, and their sequences were determined. From comparisons of diagnostic nucleotides among the short interspersed repetitive elements (SINEs) in these (sub)families, we were able to draw the following conclusions. (1) After the divergence of the suborder Tylopoda (camels), the CHRS family of SINEs was newly created from tRNA(Glu) in a common ancestor of the lineages of the Suina (pigs and peccaries), Ruminantia (cows and deer), and Cetacea (whales and dolphins). (2) After divergence of the Suina lineage, the CHR-1 SINE and the CHR-2 SINE were generated successively in a common ancestor of ruminants, hippopotamuses, and cetaceans. (3) In the Ruminantia lineage, the Bov-tA SINE was generated by recombination between the CHR-2 SINE and Bov-A. (4) In the Suina lineage, the CHRS-S SINE was generated from the CHRS SINE. (5) In this latter lineage, the PRE-1 family of SINEs was created by insertion of part of the gene for tRNA(Arg) into the 5' region of the CHRS-S family. The distribution of a particular family of SINEs among species of artiodactyls and cetaceans confirmed the most recent conclusion for paraphyly of the order Artiodactyla. The present study also revealed that a newly created tRNA(Glu)-derived family of SINEs was subjected both to recombination with different units and to duplication of an internal sequence within a SINE unit to generate, during evolution, a huge superfamily of tRNA(Glu)-related families of SINEs that are now found in the genomes of artiodactyls and cetaceans.

A remarkable example of a misleading mitochondrial protein tree is presented, involving ray-finned fishes, coelacanths, lungfishes, and tetrapods, with sea lampreys as an outgroup. In previous molecular phylogenetic studies on the origin of tetrapods, ray-finned fishes have been assumed as an outgroup to the tetrapod/lungfish/coelacanth clade, an assumption supported by morphological evidence. Standard methods of molecular phylogenetics applied to the protein-encoding genes of mitochondria, however, give a bizarre tree in which lamprey groups with lungfish and, therefore, ray-finned fishes are not the outgroup to a tetrapod/lungfish/coelacanth clade. All of the dozens of published phylogenetic methods, including every possible modification to maximum likelihood known to us (such as inclusion of site heterogeneity and exclusion of potentially misleading hydrophobic amino acids), fail to place the ray-finned fishes in a biologically acceptable position. A likely cause of this failure may be the use of an inappropriate outgroup. Accordingly, we have determined the complete mitochondrial DNA sequence from the shark, Mustelus manazo, which we have used as an alternative and more proximal outgroup than the lamprey. Using sharks as the outgroup, lungfish appear to be the closest living relative of tetrapods, although the possibility of a lungfish/coelacanth clade being the sister group of tetrapods cannot be excluded.

AFC short interspersed elements (SINEs) were isolated from cichlids from Madagascar, the New World, and Africa and characterized. A new family of long interspersed elements (LINEs), designated the CiLINE2 family, was also isolated from African cichlids, and its consensus sequence was deduced. Upon aligning all of the consensus sequences, we found that the 3'-tail regions of the AFC SINEs and the CiLINE2 family were very similar, providing another example in which a reverse transcriptase responsible for retroposition of SINEs might be contributed in trans by a LINE. Sequence comparisons showed that CiLINE2 in cichlids was closely related to LINE2 in mammals. Furthermore, we found that the 3'-tail sequence shared by the AFC SINEs and CiLINE2 in cichlids was very similar to the 3'-tail sequence shared by the MIR SINEs and LINE2 in mammals, even though the remaining parts of the AFC SINEs and the MIR SINEs were totally different from each other. Thus, the present report not only describes a new pair of SINEs and LINEs with the same 3' tail in cichlids, but also provides a new example of the phenomenon whereby the 3' ends of LINEs with the same genealogical origin can be incorporated into the 3'-end tails of different families of SINEs that have been generated independently in two different lineages during evolution.

Short interspersed repetitive elements (SINEs) have been shown to be excellent markers of molecular phylogeny, since the integration of a SINE at a particular position in a genome can be considered an unambiguous derived homologous character. In the present study, we isolated a new family of SINEs from cichlids in Lake Tanganyika, whose speciation and diversification have been regarded as prime examples of explosive adaptive radiation. Members of this new SINE family, which we named the AFC family, are about 320 bp in length, and each has a tRNA-related region in its 5' region, as do most of the members of SINE families reported to date. A dot blot hybridization experiment showed that this family is distributed extensively in the genomes of cichlids in Africa, with estimated copy numbers of 2 x 10(3)-2 x 10(4) per haploid genome. Our investigations of the patterns of insertion of members of this family at six orthologous loci demonstrated clearly that four previously identified tribes, namely, the Lamprologini, Ectodini, Tropheini, and Perissodini, each form a monophyletic group. These results provide a basis for the elucidation of the phylogenetic framework of the cichlid fishes in Lake Tanganyika.

CR1 elements are a family of retroposons. They are classified as long interspersed elements (LINEs) or non-long-terminal-repeat (non-LTR) retrotransposons, and they have been found in the genomes of many vertebrates. However, they have been only partially characterized, and only a 2-kb region of the 3' end of chicken CR1 has been sequenced. In the present study, we determined the entire consensus sequence of CR1 elements in the turtle genome, designated PsCR1. The first open reading frame (ORF1) of PsCR1 has two unusual arrangements of Cys residues. One of them includes a zinc finger motif, CX2CX14CX2C. The putative zinc finger has cysteine residues with identical spacing and a similar amino acid composition to those found in the species-specific transcription initiation factors SL1 and TIF-IB. The 5' untranslated region (5' UTR) of PsCR1 contains a sequence similar to part of the human L1 promoter, L1 site A, and several cis elements of the type found in eukaryotic genes. Within a region of about 500 bp, there are nine "E boxes," cis elements that are recognized by the basic helix-loop-helix (bHLH) family of proteins. This observation raises the possibility that cellular transcription factors that bind to these sequences might act in concert to regulate the expression of PsCR1. The extent of the sequence divergence of the 3' UTR of CR1 between species was found to be lower than the rate of nonsynonymous substitutions per site in ORF2, suggesting that a strict functional constraint must exist for this region. This result strongly suggests that the conserved 3'-end sequence of CR1 is the recognition site for the reverse transcriptase of CR1. A discussion is presented of a possible mechanism for the integration of CR1 elements and also of the intriguing possible recruitment of the reverse transcriptase for the retroposition of SINEs.