Localized Expression of Olfactory Receptor Genes in the Olfactory Organ of Common Minke Whales

Baleen whales (Mysticeti) possess the necessary anatomical structures and genetic elements for olfaction. Nevertheless, the olfactory receptor gene (OR) repertoire has undergone substantial degeneration in the cetacean lineage following the divergence of the Artiodactyla and Cetacea. The functionality of highly degenerated mysticete ORs within their olfactory epithelium remains unknown. In this study, we extracted total RNA from the nasal mucosae of common minke whales (Balaenoptera acutorostrata) to investigate ORs' localized expression. All three sections of the mucosae examined in the nasal chamber displayed comparable histological structure. However, the posterior portion of the frontoturbinal region exhibited notably high OR expression. Neither the olfactory bulb nor the external skin exhibited the expression of these genes. Although this species possesses four intact non-class-2 ORs, all the ORs expressed in the nasal mucosae belong to class-2, implying the loss of aversion to specific odorants. These anatomical and genomic analyses suggest that ORs are still responsible for olfaction within the nasal region of baleen whales, enabling them to detect desirable scents such as prey and potential mating partners.

In situ hybridization analysis of odorant receptor expression in the olfactory organ of the pig-nosed turtle Carettochelys insculpta

The turtle olfactory organ consists of upper (UCE) and lower (LCE) chamber epithelium, which send axons to the ventral and dorsal portions of the olfactory bulbs, respectively. Generally, the UCE is associated with glands and contains ciliated olfactory receptor neurons (ORNs), while the LCE is devoid of glands and contains microvillous ORNs. However, the olfactory organ of the pig-nosed turtle Carettochelys insculpta appears to be a single olfactory system morphologically: there are no associated glands; ciliated ORNs are distributed throughout the olfactory organ; and the olfactory bulb is not divided into ventral and dorsal portions. In this study, we analyzed the expression of odorant receptors (ORs), the major olfactory receptors in turtles, in the pig-nosed turtle olfactory organ, via in situ hybridization. Of 690 ORs, 375 were classified as class I and 315 as class II. Some class II ORs were expressed predominantly in the posterior dorsomedial walls of the nasal cavity, while other class II ORs and all class I ORs examined were expressed in the remaining region. These results suggest that the pig-nosed turtle olfactory organ can be divided into two regions according to the expression of ORs.

Unveiling the expansion of keratin genes in lungfishes: a possible link to terrestrial adaptation

Keratins are intermediate filament proteins that are important for epidermal strength and protection from desiccation. Keratin genes are highly duplicated and have diversified by forming two major clusters in the genomes of terrestrial vertebrates. The keratin genes of lungfishes, the closest fish to tetrapods, have not been studied at the genomic level, despite the importance of lungfishes in terrestrial adaptation. Here, we identified keratin genes in the genomes of two lungfish species and performed syntenic and phylogenetic analyses. Additionally, we identified keratin genes from two gobies and two mudskippers, inhabiting underwater and terrestrial environments. We found that in lungfishes, keratin genes were duplicated and diversified within two major clusters, similar to but independent of terrestrial vertebrates. By contrast, keratin genes were not notably duplicated in mudskippers. The results indicate that keratin gene duplication occurred repeatedly in lineages close to tetrapods, but not in teleost fish, even in species adapted to terrestrial environments.

Expression of type 1 vomeronasal receptors in the olfactory organ of the African lungfish, Protopterus dolloi

The vomeronasal organ is an olfactory organ found in amphibians and higher vertebrates. Type 1 vomeronasal receptors, one of the major olfactory receptors in vertebrates, are expressed in the vomeronasal organ in mammals. In amphibians and fish, they are expressed in the olfactory epithelium. The lungfish, which is the species of fish most closely related to amphibians, has a primitive vomeronasal organ: the recess epithelium. Expression of type 1 vomeronasal receptors has been reported in both the olfactory epithelium and the recess epithelium in three species of African lungfish and one species of South American lungfish. However, a previous study suggested that in the African lungfish Protopterus dolloi these receptors are expressed only in the olfactory epithelium. In this study, we identified 21 type 1 vomeronasal receptor genes in P. dolloi and examined the expression sites in the olfactory organ. In P. dolloi, most cells expressing the type 1 vomeronasal receptor were distributed in the olfactory epithelium, but a few were also found in the recess epithelium. This implies that the functions of the olfactory epithelium and the primitive vomeronasal organ are incompletely separated, and that all extant African and South American lungfish share this trait.

Ancient standing genetic variation facilitated the adaptive radiation of Lake Victoria cichlids

Cichlid fishes are textbook examples of explosive speciation and adaptive radiation, providing a great opportunity to understand how the genomic substrate yields extraordinary species diversity. Recently, we performed comparative genomic analyses of three Lake Victoria cichlids to reveal the genomic substrates underlying their rapid speciation and adaptation. We found that long divergent haplotypes derived from large-scale standing genetic variation, which originated before the adaptive radiation of Lake Victoria cichlids, may have contributed to their rapid diversification. In addition, the present study on genomic data from other East African cichlids suggested the reuse of alleles that may have originated in the ancestral lineages of Lake Tanganyika cichlids during cichlid evolution. Therefore, our results highlight that the primary factor that could drive repeated adaptive radiation across East African cichlids was allelic reuse from standing genetic variation to adapt to their own specific environment. In this report, we summarize the main results and discuss the evolutionary mechanisms of cichlids, based on our latest findings.

Plastic loss of motile cilia in the gills of Polypterus in response to high CO2 or terrestrial environments

The evolutionary transition of vertebrates from water to land during the Devonian period was accompanied by major changes in animal respiratory systems in terms of physiology and morphology. Indeed, the fossil record of the early tetrapods has revealed the existence of internal gills, which are vestigial fish-like traits used underwater. However, the fossil record provides only limited data on the process of the evolutionary transition of gills from fish to early tetrapods. This study investigated the gills of Polypterus senegalus, a basal ray-finned/amphibious fish which shows many ancestral features of stem Osteichthyes. Based on scanning electron microscopy observations and transcriptome analysis, the existence of motile cilia in the gills was revealed which may create a flow on the gill surface leading to efficient ventilation or remove particles from the surface. Interestingly, these cilia were observed to disappear after rearing in terrestrial or high CO₂ environments, which mimics the environmental changes in the Devonian period. The cilia re-appeared after being returned to the original aquatic environment. The ability of plastic changes of gills in Polypterus revealed in this study may allow them to survive in fluctuating environments, such as shallow swamps. The ancestor of Osteichthyes is expected to have possessed such plasticity in the gills, which may be one of the driving forces behind the transition of vertebrates from water to land.

Type 1 vomeronasal receptor expression in juvenile and adult lungfish olfactory organ

Lungfish are the most closely related fish to tetrapods. The olfactory organ of lungfish contains lamellae and abundant recesses at the base of lamellae. Based on the ultrastructural and histochemical characteristics, the lamellar olfactory epithelium (OE), covering the surface of lamellae, and the recess epithelium, contained in the recesses, are thought to correspond to the OE of teleosts and the vomeronasal organ (VNO) of tetrapods. With increasing body size, the recesses increase in number and distribution range in the olfactory organ. In tetrapods, the expression of olfactory receptors is different between the OE and VNO; for instance, the type 1 vomeronasal receptor (V1R) is expressed only in the OE in amphibians and mainly in the VNO in mammals. We recently reported that V1R-expressing cells are contained mainly in the lamellar OE but also rarely in the recess epithelium in the olfactory organ of lungfish of approximately 30 cm body length. However, it is unclear whether the distribution of V1R-expressing cells in the olfactory organ varies during development. In this study, we compared the expression of V1Rs in the olfactory organs between juveniles and adults of the African lungfish Protopterus aethiopicus and South American lungfish, Lepidosiren paradoxa. The density of V1R-expressing cells was higher in the lamellae than in the recesses in all specimens evaluated, and this pattern was more pronounced in juveniles than adults. In addition, the juveniles showed a higher density of V1R-expressing cells in the lamellae compared with the adults. Our results imply that differences in lifestyle between juveniles and adults are related to differences in the density of V1R-expressing cells in the lamellae of lungfish.

c-fos expression in the olfactory epithelium of the East African cichlid (Haplochromis chilotes) in response to odorant exposure

Fishes use olfaction to gain varied information vital for survival and communication. To understand biodiversity in fishes, it is important to identify what receptors individual fish use to detect specific chemical compounds. However, studies of fish olfactory receptors and their ligands are still limited to a few model organisms represented primarily by zebrafish. Here, we tested the c-fos expression of olfactory sensory neurons (OSNs) in an East African cichlid, the most diversified teleost lineage, by in situ hybridization with a c-fos riboprobe. We confirmed that microvillous neurons contributed the most to the detection of amino acids, as in other fishes. Conversely, we found that ciliated neurons contributed the most to detection of conjugated steroids, known as pheromone candidates. We also found that V2Rs, the major receptor type in microvillous neurons, exhibited differential responsiveness to amino acids, and further suggest that the cichlid-specific duplication of V2R led to ligand differentiation by demonstrating a differential response to arginine. Finally, we established a nonlethal method to collect cichlid urine and showed how various OSNs, including V1R+ neurons, respond to male urine. This study provides an experimental basis for understanding how cichlids encode natural odors, which ultimately lead to provide insight into how olfaction has contributed to the diversification of cichlids.

Type 1 vomeronasal receptors expressed in the olfactory organs of two African lungfish, Protopterus annectens and Protopterus amphibius

Lungfish are the fish related most closely to tetrapods. The olfactory organ of lungfish contains two distinct sensory epithelia: the lamellar olfactory epithelium (OE) and the recess epithelium (RecE). Based on their ultrastructural and histological characteristics, the lamellar OE and the RecE are considered to correspond respectively to the teleost OE and a primitive vomeronasal organ (VNO). In tetrapods, the OE and VNO have been shown to express different families of olfactory receptors; for example, in mammals, the OE expresses odorant receptors and trace amine-associated receptors, while the VNO expresses type 1 (V1Rs) and type 2 (V2Rs) vomeronasal receptors. In the present study, we examined the expression of V1Rs in the olfactory organs of two African lungfish, Protopterus annectens and Protopterus amphibius. RNA sequencing and phylogenetic analyses identified 29 V1R genes in P. annectens and 50 V1R genes in P. amphibius. Most V1Rs identified in these lungfish were classified as the tetrapod-type V1Rs initially found in tetrapods and distinct from fish-type V1Rs. In teleost, which all lack a VNO, all olfactory receptors are expressed in the OE, while in Xenopus V1Rs are expressed exclusively in the OE, and not in the VNO. In situ hybridization analysis indicated that lungfish V1Rs were expressed mainly in the lamellar OE and rarely in the RecE. These results imply that V1R expression in lungfish represents an intermediate step toward the complete segregation of V1R expression between the OE and VNO, reflecting the phylogenetic position of lungfish between teleosts and amphibians.

Alternative splicing plays key roles in response to stress across different stages of fighting in the fish Betta splendens

Background

Aggression is an evolutionarily conserved behavior critical for animal survival. In the fish Betta splendens, across different stages of fighting interactions, fighting opponents suffer from various stressors, especially from the great demand for oxygen. Using RNA sequencing, we profiled differential alternative splicing (DAS) events in the brains of fish collected before fighting, during fighting, and after fighting to study the involvement of alternative splicing (AS) in the response to stress during the fight.

Results

We found that fighting interactions induced the greatest increase in AS in the ‘during-fighting’ fish, followed by that of the ‘after-fighting’ fish. Intron retention (IR) was the most enriched type among all the basic AS events. DAS genes were mainly associated with synapse assembly, ion transport, and regulation of protein secretion. We further observed that IR events significantly differentiated between winners and losers for 19 genes, which were associated with messenger RNA biogenesis, DNA repair, and transcription machinery. These genes share many common features, including shorter intron length and higher GC content.

Conclusions

This study is the first comprehensive view of AS induced by fighting interactions in a fish species across different stages of those interactions, especially with respect to IR events in winners and losers. Together, these findings facilitate future investigations into transcriptome complexity and AS regulation in response to stress under the context of aggression in vertebrates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08609-2.

New Sex Chromosomes in Lake Victoria Cichlid Fishes (Cichlidae: Haplochromini)

African cichlid fishes harbor an extraordinary diversity of sex-chromosome systems. Within just one lineage, the tribe Haplochromini, at least 6 unique sex-chromosome systems have been identified. Here we focus on characterizing sex chromosomes in cichlids from the Lake Victoria basin. In Haplochromis chilotes, we identified a new ZW system associated with the white blotch color pattern, which shows substantial sequence differentiation over most of LG16, and is likely to be present in related species. In Haplochromis sauvagei, we found a coding polymorphism in amh that may be responsible for an XY system on LG23. In Pundamilia nyererei, we identified a feminizing effect of B chromosomes together with XY- and ZW-patterned differentiation on LG23. In Haplochromis latifasciatus, we identified a duplication of amh that may be present in other species of the Lake Victoria superflock. We further characterized the LG5-14 XY system in Astatotilapia burtoni and identified the oldest stratum on LG14. This species also showed ZW differentiation on LG2. Finally, we characterized an XY system on LG7 in Astatoreochromis alluaudi. This report brings the number of distinct sex-chromosome systems in haplochromine cichlids to at least 13, and highlights the dynamic evolution of sex determination and sex chromosomes in this young lineage.

A transcriptomic reevaluation of the accessory olfactory organ in Bichir (Polypterus senegalus)

Fish possess one olfactory organ called the olfactory epithelium (OE), by which various chemical substances are detected. On the other hand, tetrapods possess two independent olfactory organs called the main olfactory epithelium (MOE) and vomeronasal organ (VNO), each of which mainly detects general odorants and pheromones, respectively. Traditionally, the VNO, so-called concentrations of vomeronasal neurons, was believed to have originated in tetrapods. However, recent studies have identified a primordial VNO in lungfish, implying that the origin of the VNO was earlier than traditionally expected. In this study, we examined the presence/absence of the VNO in the olfactory organ of bichir (Polypterus senegalus), which is the most ancestral group of extant bony vertebrates. In particular, we conducted a transcriptomic evaluation of the accessory olfactory organ (AOO), which is anatomically separated from the main olfactory organ (MOO) in bichir. As a result, several landmark genes specific to the VNO and MOE in tetrapods were both expressed in the MOO and AOO, suggesting that these organs were not functionally distinct in terms of pheromone and odorant detection. Instead, differentially expressed gene (DEG) analysis showed that DEGs in AOO were enriched in genes for cilia movement, implying its additional and specific function in efficient water uptake into the nasal cavity other than chemosensing. This transcriptomic study provides novel insight into the long-standing question of AOO function in bichir and suggests that VNO originated in the lineage of lobe-finned fish during vertebrate evolution.

Data of RNA-seq transcriptomes in the brain associated with aggression in males of the fish Betta splendens

Siamese fighting fish Betta splendens are notorious for their aggressiveness and males of this fish have been widely used to study aggression. However, an understanding of brain transcriptome signature associated with aggression in the context of male-male interaction in this fish remains to be understood. Herein, RNA-Seq transcriptome data from 37 brains samples collected at different fighting stages are described. These brain samples were collected before fighting (B), during fighting (D20 and D60), and after fighting (A0 and A30). The raw data were analyzed for differential gene expression using edgeR package in R. A criterion of FDR cut-off ≤ 0.05 and an absolute fold change (FC) of 0 or greater were used to identify top upregulated and downregulated genes in fighting groups (D20, D60, A0, and A30) relative to non-fighting group (B). The data presented hereafter enable fundamental studies on genes and molecular events mediating aggressive behavior in this fish and will lay a valuable foundation for future research on the aggression of vertebrates.

Inactivation of ancV1R as a Predictive Signature for the Loss of Vomeronasal System in Mammals

The vomeronasal organ (VNO) plays a key role in sensing pheromonal cues, which elicits social and reproductive behaviors. Although the VNO is highly conserved across mammals, it has been lost in some species that have evolved alternate sensing systems during diversification. In this study, we investigate a newly identified VNO-specific gene, ancV1R, in the extant 261 species of mammals to examine the correlation between genotype (ancV1R) and phenotype (VNO). As a result, we found signatures for the relaxation of purifying selection (inactivating mutations and the elevation of dN/dS) on ancV1Rs in VNO-lacking mammals, such as catarrhine primates, cetaceans, the manatees, and several bat lineages, showing the distinct correlation between genotype and phenotype. Interestingly, we further revealed signatures for the relaxation of purifying selection on ancV1R in true seals, otters, the fossa, the owl monkey, and alcelaphine antelopes in which the existence of a functional VNO is still under debate. Our additional analyses on TRPC2, another predictive marker gene for the functional VNO, showed a relaxation of purifying selection, supporting the possibility of VNO loss in these species. The results of our present study invite more in-depth neuroanatomical investigation in mammals for which VNO function remains equivocal.

Platypus and echidna genomes reveal mammalian biology and evolution

Egg-laying mammals (monotremes) are the only extant mammalian outgroup to therians (marsupial and eutherian animals) and provide key insights into mammalian evolution1,2. Here we generate and analyse reference genomes of the platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus), which represent the only two extant monotreme lineages. The nearly complete platypus genome assembly has anchored almost the entire genome onto chromosomes, markedly improving the genome continuity and gene annotation. Together with our echidna sequence, the genomes of the two species allow us to detect the ancestral and lineage-specific genomic changes that shape both monotreme and mammalian evolution. We provide evidence that the monotreme sex chromosome complex originated from an ancestral chromosome ring configuration. The formation of such a unique chromosome complex may have been facilitated by the unusually extensive interactions between the multi-X and multi-Y chromosomes that are shared by the autosomal homologues in humans. Further comparative genomic analyses unravel marked differences between monotremes and therians in haptoglobin genes, lactation genes and chemosensory receptor genes for smell and taste that underlie the ecological adaptation of monotremes.

A unique neurogenomic state emerges after aggressive confrontations in males of the fish Betta splendens

Territorial defense involves frequent aggressive confrontations with competitors, but little is known about how brain-transcriptomic profiles change between individuals competing for territory establishment. Our previous study elucidated that when two fish Betta splendens males interact, transcriptomes across their brains synchronize in a way that reflects a mutual assessment process between them at the gene expression level. Here we aim to evaluate how the brain-transcriptomic profiles of opponents change immediately after shifting their social status (i.e., the winner/loser has emerged) and 30 min after this shift. We showed that changes in the expression of certain genes are unique to different fighting stages and the expression patterns of certain genes are transiently or persistently changed across all fighting stages. These brain transcriptomic responses are in accordance with behavioral changes across the fight. Strikingly, the specificity of the brain-transcriptomic synchronization of a pair during fighting was gradually lost after fighting ceased, leading to the emergence of a basal neurogenomic state in which the changes in gene expression were reduced to minimum and consistent across all individuals. This state shares common characteristics with the hibernation state that animals adopt to minimize their metabolic rates to save energy. Interestingly, expression changes for genes related to metabolism, autism spectrum disorder, and long-term memory still differentiated losers from winners. Together, the fighting system using male B. splendens provides a promising platform for investigating neurogenomic states of aggression in vertebrates.

Genomic Signatures for Species-Specific Adaptation in Lake Victoria Cichlids Derived from Large-Scale Standing Genetic Variation

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.

Conserved keratin gene clusters in ancient fish: An evolutionary seed for terrestrial adaptation

Type I and type II keratins are subgroups of intermediate filament proteins that provide toughness to the epidermis and protect it from water loss. In terrestrial vertebrates, the keratin genes form two major clusters, clusters 1 and 2, each of which is dominated by type I and II keratin genes. By contrast, such clusters are not observed in teleost fish. Although the diversification of keratins is believed to have made a substantial contribution to terrestrial adaptation, its evolutionary process has not been clarified. Here, we performed a comprehensive genomic survey of the keratin genes of a broad range of vertebrates. As a result, we found that ancient fish lineages such as elephant shark, reedfish, spotted gar, and coelacanth share both keratin gene clusters. We also discovered an expansion of keratin genes that form a novel subcluster in reedfish. Syntenic and phylogenetic analyses revealed that two pairs of krt18/krt8 keratin genes were shared among all vertebrates, thus implying that they encode ancestral type I and II keratin protein sets. We further revealed that distinct keratin gene subclusters, which show specific expressions in the epidermis of adult amphibians, stemmed from canonical keratin genes in non-terrestrial ancestors. Molecular evolutionary analyses suggested that the selective constraints were relaxed in the adult epidermal subclusters of amphibians as well as the novel subcluster of reedfish. The results of the present study represent the process of diversification of keratins through a series of gene duplications that could have facilitated the terrestrial adaptation of vertebrates.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister–taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault’s rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Evolution of V1R pheromone receptor genes in vertebrates: diversity and commonality

The vomeronasal organ (VNO) plays a key role in sensing pheromonal cues, which elicit innate responses and induce social and sexual behaviors. The vomeronasal receptor 1 genes, V1Rs, encode members of a pheromone receptor family that are mainly expressed in the VNO. Previous studies have revealed that the V1R family shows extraordinary variety among mammalian species owing to successive gene gains and losses. Because species-specific pheromonal interaction may facilitate species-specific reproductive behaviors, understanding the evolution of V1Rs in terms of their origin, repertoire and phylogeny should provide insight into the mechanisms of animal diversification. Here I summarize recent studies about the V1R family from its initial discovery in the rat genome to extensive comparative analyses among vertebrates. I further introduce our recent findings for V1Rs in a broad range of vertebrates, which reveal unexpected diversity as well as shared features common among lineages.

A Single Pheromone Receptor Gene Conserved across 400 My of Vertebrate Evolution

Pheromones are crucial for eliciting social and sexual behaviors in diverse animal species. The vomeronasal receptor type-1 (V1R) genes, encoding members of a pheromone receptor family, are highly variable in number and repertoire among mammals due to extensive gene gain and loss. Here, we report a novel pheromone receptor gene belonging to the V1R family, named ancient V1R (ancV1R), which is shared among most Osteichthyes (bony vertebrates) from the basal lineage of ray-finned fishes to mammals. Phylogenetic and syntenic analyses of ancV1R using 115 vertebrate genomes revealed that it represents an orthologous gene conserved for >400 My of vertebrate evolution. Interestingly, the loss of ancV1R in some tetrapods is coincident with the degeneration of the vomeronasal organ in higher primates, cetaceans, and some reptiles including birds and crocodilians. In addition, ancV1R is expressed in most mature vomeronasal sensory neurons in contrast with canonical V1Rs, which are sparsely expressed in a manner that is consistent with the "one neuron-one receptor" rule. Our results imply that a previously undescribed V1R gene inherited from an ancient Silurian ancestor may have played an important functional role in the evolution of vertebrate vomeronasal organ.

Patterns of genomic differentiation between two Lake Victoria cichlid species, Haplochromis pyrrhocephalus and H. sp. 'macula'

BACKGROUND

The molecular basis of the incipient stage of speciation is still poorly understood. Cichlid fish species in Lake Victoria are a prime example of recent speciation events and a suitable system to study the adaptation and reproductive isolation of species.

RESULTS

Here, we report the pattern of genomic differentiation between two Lake Victoria cichlid species collected in sympatry, Haplochromis pyrrhocephalus and H. sp. 'macula,' based on the pooled genome sequences of 20 individuals of each species. Despite their ecological differences, population genomics analyses demonstrate that the two species are very close to a single panmictic population due to extensive gene flow. However, we identified 21 highly differentiated short genomic regions with fixed nucleotide differences. At least 15 of these regions contained genes with predicted roles in adaptation and reproductive isolation, such as visual adaptation, circadian clock, developmental processes, adaptation to hypoxia, and sexual selection. The nonsynonymous fixed differences in one of these genes, LWS, were reported as substitutions causing shift in absorption spectra of LWS pigments. Fixed differences were found in the promoter regions of four other differentially expressed genes, indicating that these substitutions may alter gene expression levels.

CONCLUSIONS

These diverged short genomic regions may have contributed to the differentiation of two ecologically different species. Moreover, the origins of adaptive variants within the differentiated regions predate the geological formation of Lake Victoria; thus Lake Victoria cichlid species diversified via selection on standing genetic variation.

Evolution of vomeronasal receptor 1 (V1R) genes in the common marmoset (Callithrix jacchus)

Pheromones are crucial for eliciting innate responses and inducing social and sexual behaviors in mammals. The vomeronasal receptor 1 genes, V1Rs, encode members of a pheromone receptor family that are mainly expressed in the vomeronasal organ (VNO). The V1R family shows extraordinary variety in gene number among vertebrates owing to successive gene gains and losses during evolution. Such diversity is believed to reflect a degree of dependence on the VNO. We investigated V1R evolution in primate lineages closely related to humans because these VNOs show a trend toward degeneration. We performed extensive phylogenetic analyses for V1Rs from a broad range of primate species. Although the decline of intact genes was evident in anthropoids (hominoids, Old World monkeys and New World monkeys), we found that a certain number of intact genes persist in New World monkeys. In one New World monkey species, the common marmoset (Callithrix jacchus), we examined seven putatively functional V1Rs using in situ hybridization and reverse transcription-PCR. Based on their mRNA expression patterns in the VNO and other organs, two types of V1Rs emerged: the canonical class with VNO-specific expression, and a second group having more ubiquitous expression in various organs as well as VNO. Furthermore, phylogenetic analysis revealed that the class with the more widespread expression had been retained longer in evolution than the VNO-specific type. We propose that the acquisition of a novel non-VNO-related function(s) may have led to the survival of a small but persistent number of V1Rs in anthropoid primates.

Stim1 Regulates Enamel Mineralization and Ameloblast Modulation

Loss-of-function mutations in the Ca²⁺ release-activated Ca²⁺ channel genes ORAI1 and STIM1 abolish store-operated Ca²⁺ entry (SOCE) and result in ectodermal dysplasia with amelogenesis imperfecta. However, because of the limited availability of patient tissue, analyses of enamel mineralization or possible changes in ameloblast function or morphology have not been possible. Here, we generated mice with ectodermal tissue-specific deletion of Stim1 ( Stim1 cKO [conditional knockout]), Stim2 ( Stim2 cKO), and Stim1 and Stim2 ( Stim1/2 cKO) and analyzed their enamel phenotypes as compared with those of control ( Stim1/2^fl/fl) animals. Ablation of Stim1 and Stim1/2 but not Stim2 expression resulted in chalky enamel and severe attrition at the incisor tips and molar cusps. Stim1 and Stim1/2 cKO, but not Stim2 cKO, demonstrated inferior enamel mineralization with impaired structural integrity, whereas the shape of the teeth and enamel thickness appeared to be normal in all animals. The gene expression levels of the enamel matrix proteins Amelx and Ambn and the enamel matrix proteases Mmp20 and Klk4 were not altered by the abrogation of SOCE in Stim1/2 cKO mice. The morphology of ameloblasts during the secretory and maturation stages was not significantly altered in either the incisors or molars of the cKO animals. However, in Stim1 and Stim1/2 cKO incisors, the alternating modulation of maturation-stage ameloblasts between the smooth- and ruffle-ended cell types continued beyond the regular cycle and extended to the areas corresponding to the zone of postmodulation ameloblasts in the teeth of control animals. These results indicate that SOCE is essential for proper enamel mineralization, in which Stim1 plays a critical role during the maturation process.

Environmental DNA enables detection of terrestrial mammals from forest pond water

Terrestrial animals must have frequent contact with water to survive, implying that environmental DNA (eDNA) originating from those animals should be detectable from places containing water in terrestrial ecosystems. Aiming to detect the presence of terrestrial mammals using forest water samples, we applied a set of universal PCR primers (MiMammal, a modified version of fish universal primers) for metabarcoding mammalian eDNA. The versatility of MiMammal primers was tested in silico and by amplifying DNAs extracted from tissues. The results suggested that MiMammal primers are capable of amplifying and distinguishing a diverse group of mammalian species. In addition, analyses of water samples from zoo cages of mammals with known species composition suggested that MiMammal primers could successfully detect mammalian species from water samples in the field. Then, we performed an experiment to detect mammals from natural ecosystems by collecting five 500-ml water samples from ponds in two cool-temperate forests in Hokkaido, northern Japan. MiMammal amplicon libraries were constructed using eDNA extracted from water samples, and sequences generated by Illumina MiSeq were subjected to data processing and taxonomic assignment. We thereby detected multiple species of mammals common to the sampling areas, including deer (Cervus nippon), mouse (Mus musculus), vole (Myodes rufocanus), raccoon (Procyon lotor), rat (Rattus norvegicus) and shrew (Sorex unguiculatus). Many previous applications of the eDNA metabarcoding approach have been limited to aquatic/semiaquatic systems, but the results presented here show that the approach is also promising even for forest mammal biodiversity surveys.

Distinct functions of two olfactory marker protein genes derived from teleost-specific whole genome duplication

BACKGROUND

Whole genome duplications (WGDs) have been proposed to have made a significant impact on vertebrate evolution. Two rounds of WGD (1R and 2R) occurred in the common ancestor of Gnathostomata and Cyclostomata, followed by the third-round WGD (3R) in a common ancestor of all modern teleosts. The 3R-derived paralogs are good models for understanding the evolution of genes after WGD, which have the potential to facilitate phenotypic diversification. However, the recent studies of 3R-derived paralogs tend to be based on in silico analyses. Here we analyzed the paralogs encoding teleost olfactory marker protein (OMP), which was shown to be specifically expressed in mature olfactory sensory neurons and is expected to be involved in olfactory transduction.

RESULTS

Our genome database search identified two OMPs (OMP1 and OMP2) in teleosts, whereas only one was present in other vertebrates. Phylogenetic and synteny analyses suggested that OMP1 and 2 were derived from 3R. Both OMPs showed distinct expression patterns in zebrafish; OMP1 was expressed in the deep layer of the olfactory epithelium (OE), which is consistent with previous studies of mice and zebrafish, whereas OMP2 was sporadically expressed in the superficial layer. Interestingly, OMP2 was expressed in a very restricted region of the retina as well as in the OE. In addition, the analysis of transcriptome data of spotted gar, a non-teleost fish, revealed that single OMP gene was expressed in the eyes.

CONCLUSION

We found distinct expression patterns of zebrafish OMP1 and 2 at the tissue and cellular level. These differences in expression patterns may be explained by subfunctionalization as the model of molecular evolution. Namely, single OMP gene was speculated to be originally expressed in the OE and the eyes in the common ancestor of all Osteichthyes (bony fish including tetrapods). Then, two OMP gene paralogs derived from 3R-WGD reduced and specialized the expression patterns. This study provides a good example for analyzing a functional subdivision of the teleost OE and eyes as revealed by 3R-derived paralogs of OMPs.

Multiple episodic evolution events in V1R receptor genes of East-African cichlids

Fish use olfaction to detect a variety of nonvolatile chemical signals, and thus, this sense is key to survival and communication. However, the contribution of the olfactory sense to social—especially reproductive—interactions in cichlids is still controversial. To obtain insights into this issue, we investigated the genes encoding V1Rs—possible candidates for reproductive pheromone receptors—among East-African cichlids. Interestingly, we found an excess of nonsynonymous over synonymous substitutions in four of six V1R genes in multiple cichlid lineages. First, we found that highly dimorphic V1R2 allele groups were shared among the cichlids inhabiting all East-African Great Lakes emerged through the episodic accumulation of the nonsynonymous substitutions prior to the radiation of the Lake Tanganyika species flock. We further detected such episodic events in V1R1 of the tribe Tropheini, and in V1R3 and V1R6 of the tribe Trematocarini. The excess of nonsynonymous substitutions in these examples were indicated as dN/dS > 1, which were all statistically significant by Fisher’s exact test. Furthermore, we speculate that the amino acid changes in these episodic events are likely functional switch because they occurred in the putative ligand-binding pocket. Our finding of the occurrence of multiple episodic events and the unexpected gene diversity in one unique gene family is suggestive of the contribution of the V1R to the species diversification and the social interaction in cichlids.

Lineage-specific expansion of vomeronasal type 2 receptor-like (OlfC) genes in cichlids may contribute to diversification of amino acid detection systems

Fish use olfaction to sense a variety of nonvolatile chemical signals in water. However, the evolutionary importance of olfaction in species-rich cichlids is controversial. Here, we determined an almost complete sequence of the vomeronasal type 2 receptor-like (OlfC: putative amino acids receptor in teleosts) gene cluster using the bacterial artificial chromosome library of the Lake Victoria cichlid, Haplochromis chilotes. In the cluster region, we found 61 intact OlfC genes, which is the largest number of OlfC genes identified among the seven teleost fish investigated to date. Data mining of the Oreochromis niloticus (Nile tilapia) draft genome sequence, and genomic Southern hybridization analysis revealed that the ancestor of all modern cichlids had already developed almost the same OlfC gene repertoire, which was accomplished by lineage-specific gene expansions. Furthermore, comparison of receptor sequences showed that recently duplicated paralogs are more variable than orthologs of different species at particular sites that were predicted to be involved in amino acid selectivity. Thus, the increase of paralogs through gene expansion may lead to functional diversification in detection of amino acids. This study implies that cichlids have developed a potent capacity to detect a variety of amino acids (and their derivatives) through OlfCs, which may have contributed to the extraordinary diversity of their feeding habitats.

Coelacanth genomes reveal signatures for evolutionary transition from water to land

Coelacanths are known as “living fossils,” as they show remarkable morphological resemblance to the fossil record and belong to the most primitive lineage of living Sarcopterygii (lobe-finned fishes and tetrapods). Coelacanths may be key to elucidating the tempo and mode of evolution from fish to tetrapods. Here, we report the genome sequences of five coelacanths, including four Latimeria chalumnae individuals (three specimens from Tanzania and one from Comoros) and one L. menadoensis individual from Indonesia. These sequences cover two African breeding populations and two known extant coelacanth species. The genome is ∼2.74 Gbp and contains a high proportion (∼60%) of repetitive elements. The genetic diversity among the individuals was extremely low, suggesting a small population size and/or a slow rate of evolution. We found a substantial number of genes that encode olfactory and pheromone receptors with features characteristic of tetrapod receptors for the detection of airborne ligands. We also found that limb enhancers of bmp7 and gli3, both of which are essential for limb formation, are conserved between coelacanth and tetrapods, but not ray-finned fishes. We expect that some tetrapod-like genes may have existed early in the evolution of primitive Sarcopterygii and were later co-opted to adapt to terrestrial environments. These coelacanth genomes will provide a cornerstone for studies to elucidate how ancestral aquatic vertebrates evolved into terrestrial animals.

The complete mitochondrial genomes of deep-sea squid (Bathyteuthis abyssicola), bob-tail squid (Semirossia patagonica) and four giant cuttlefish (Sepia apama, S. latimanus, S. lycidas and S. pharaonis), and their application to the phylogenetic analysis of Decapodiformes

We determined the complete mitochondrial (mt) genomes of the deep-sea squid (Bathyteuthis abyssicola; supperfamily Bathyteuthoidea), the bob-tail squid (Semirossia patagonica; order Sepiolida) and four giant cuttlefish (Sepia apama, S. latimanus, S. lycidas and S. pharaonis; order Sepiida). The unique structures of the mt genomes of Bathyteuthis and Semirossia provide new information about the evolution of decapodiform mt genomes. We show that the mt genome of B. abyssicola, like those of other oegopsids studied so far, has two long duplicated regions that include seven genes (COX1-3, ATP6 and ATP8, tRNA(Asn), and either ND2 or ND3) and that one of the duplicated COX3 genes has lost its function. The mt genome of S. patagonica is unlike any other decapodiforms and, like Nautilus, its ATP6 and ATP8 genes are not adjacent to each other. The four giant cuttlefish have identical mt gene order to other cuttlefish determined to date. Molecular phylogenetic analyses using maximum likelihood and Bayesian methods suggest that traditional order Sepioidea (Sepiolida+Sepiida) is paraphyletic and Sepia (cuttlefish) has the sister-relationship with all other decapodiforms. Taking both the phylogenetic analyses and the mt gene order analyses into account, it is likely that the octopus-type mt genome is an ancestral state and that it had maintained from at least the Cephalopoda ancestor to the common ancestor of Oegopsida, Myopsida and Sepiolida.

Characterization of V1R receptor (ora) genes in Lake Victoria cichlids

Although olfaction could play a crucial role in underwater habitats by allowing fish to sense a variety of nonvolatile chemical signals, the importance of olfaction in species-rich cichlids is still controversial. In particular, examining whether cichlids rely on olfaction for reproduction is of primary interest to understand the mechanisms of speciation. In the present study, we explored the V1R (also known as ora) genes, which are believed to encode reproductive pheromone receptors in fish, in the genomes of Lake Victoria cichlids. By screening a bacterial artificial chromosome library, we identified all six intact V1R genes (V1R1 to V1R6) that have been reported in other teleost fish. Furthermore, RT-PCR and in situ hybridization analyses showed that all of the V1R genes were expressed in the olfactory epithelium, indicating that these receptors are functional in cichlids. These observations indicate that cichlids use V1R-mediated olfaction in some ways for their social behaviors.

Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last approximately 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events.

Characterization of a novel SINE superfamily from invertebrates: "Ceph-SINEs" from the genomes of squids and cuttlefish

Five tRNA-derived short interspersed repetitive elements (SINEs), named SepiaSINE, Sepioth-SINE1, Sepioth-SINE2A, Sepioth-SINE2B and OegopSINE, were isolated from the genomes of three decabrachian species [Sepia officinalis (order Sepiida), Sepiotheuthis lessoniana (suborder Myopsida), and Mastigoteuthis cordiformes (suborder Oegopsida)], by random sequencing and genome screening. In addition, two tRNA-derived SINEs, named IdioSINE1 and IdioSINE2, were further detected from EST (expressed sequence tag) data of Idiosepius paradoxus (order Idiosepiida), using a GenBank FASTA search with a conserved sequence of the SepiaSINE as the query. All the isolated SINEs had a common and unique highly conserved 149-bp sequence in their central structures (Sepioth-SINE2B and IdioSINEs, however, had a continuous 73-bp deletion in the conserved region.), and are therefore grouped as the fourth SINE superfamily "Ceph-SINEs", following the CORE-SINE, V-SINE, and DeuSINE superfamilies. Our analysis suggested that the central conserved region called the "Ceph-domain" might have originated before the diversification of cephalopods (505 myr ago). A sequence alignment of Sepioth-SINE1, Sepioth-SINE2A, and Sepioth-SINE2B demonstrated that Sepioth-SINE2A has a chimeric structure shared with two other SINEs. The above relationship suggests possible template switching in the central conserved domain during reverse transcription for the birth of Sepioth-SINE2A, providing the possibility that the presence of the conserved domain contributed to yield a variety of SINEs during evolution. Furthermore, the distributions of the isolated SINEs showed that order Sepiida, suborders Oegopsida and Myopsida, and order Idiosepiida have their own independent SINE(s), and suggest that order Sepiida can be largely separated into two groups, with clarification of the phylogenetic relatedness between subfamily Sepioteuthinae and the other loliginid squids.

Vertebrate rhodopsin adaptation to dim light via rapid meta-II intermediate formation

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.

Conserved repertoire of orthologous vomeronasal type 1 receptor genes in ruminant species

Background

In mammals, pheromones play an important role in social and innate reproductive behavior within species. In rodents, vomeronasal receptor type 1 (V1R), which is specifically expressed in the vomeronasal organ, is thought to detect pheromones. The V1R gene repertoire differs dramatically between mammalian species, and the presence of species-specific V1R subfamilies in mouse and rat suggests that V1R plays a profound role in species-specific recognition of pheromones. In ruminants, however, the molecular mechanism(s) for pheromone perception is not well understood. Interestingly, goat male pheromone, which can induce out-of-season ovulation in anestrous females, causes the same pheromone response in sheep, and vice versa, suggesting that there may be mechanisms for detecting "inter-species" pheromones among ruminant species.

Results

We isolated 23 goat and 21 sheep intact V1R genes based on sequence similarity with 32 cow V1R genes in the cow genome database. We found that all of the goat and sheep V1R genes have orthologs in their cross-species counterparts among these three ruminant species and that the sequence identity of V1R orthologous pairs among these ruminants is much higher than that of mouse-rat V1R orthologous pairs. Furthermore, all goat V1Rs examined thus far are expressed not only in the vomeronasal organ but also in the main olfactory epithelium.

Conclusion

Our results suggest that, compared with rodents, the repertoire of orthologous V1R genes is remarkably conserved among the ruminants cow, sheep and goat. We predict that these orthologous V1Rs can detect the same or closely related chemical compound(s) within each orthologous set/pair. Furthermore, all identified goat V1Rs are expressed in the vomeronasal organ and the main olfactory epithelium, suggesting that V1R-mediated ligand information can be detected and processed by both the main and accessory olfactory systems. The fact that ruminant and rodent V1Rs have distinct features suggests that ruminant and rodent V1Rs have evolved distinct functions.

Chronology of the extant African elephant species and case study of the species identification of the small African elephant with the molecular phylogenetic method

Despite vigorous genetic studies of African elephants (Loxodonta africana and L. cyclotis) during the last decade, their evolutionary history is still obscure. Phylogenetic studies and coalescence time estimation using longer nucleotide sequence data from denser samplings are necessary to better understand the natural history of African elephants. Further, species identification among African elephants is sometimes very difficult using only the external morphological characteristics. This is a serious problem for making an adequate breeding plan in zoological gardens. In this paper, we investigated the continent-wide phylogeographical pattern of the African elephants and estimated the coalescence times among them. From these molecular data and geological evidence, we proposed an evolutionary scenario for the African elephants. We further demonstrated the effectiveness of molecular phylogenetic methods in species identification.

Newly discovered young CORE-SINEs in marsupial genomes

Although recent mammalian genome projects have uncovered a large part of genomic component of various groups, several repetitive sequences still remain to be characterized and classified for particular groups. The short interspersed repetitive elements (SINEs) distributed among marsupial genomes are one example. We have identified and characterized two new SINEs from marsupial genomes that belong to the CORE-SINE family, characterized by a highly conserved "CORE" domain. PCR and genomic dot blot analyses revealed that the distribution of each SINE shows distinct patterns among the marsupial genomes, implying different timing of their retroposition during the evolution of marsupials. The members of Mar3 (Marsupialia 3) SINE are distributed throughout the genomes of all marsupials, whereas the Mac1 (Macropodoidea 1) SINE is distributed specifically in the genomes of kangaroos. Sequence alignment of the Mar3 SINEs revealed that they can be further divided into four subgroups, each of which has diagnostic nucleotides. The insertion patterns of each SINE at particular genomic loci, together with the distribution patterns of each SINE, suggest that the Mar3 SINEs have intensively amplified after the radiation of diprotodontians, whereas the Mac1 SINE has amplified only slightly after the divergence of hypsiprimnodons from other macropods. By compiling the information of CORE-SINEs characterized to date, we propose a comprehensive picture of how SINE evolution occurred in the genomes of marsupials.

Molecular phylogenetic study on the origin and evolution of Mustelidae

The family Mustelidae, which consists of Mustelinae, Lutrinae, Melinae, and Taxidiinae, is the largest family among Carnivora and is a highly diverse group. Recent molecular phylogenetic studies have clarified the phylogenetic relations among Mustelidae, but there remain several unresolved problems, particularly concerning the deep branchings. Whereas many studies support the monophyly of Mustelidae+Procyonidae among Musteloidea, the relations between Mustelidae+Procyonidae, Ailuridae, and Miphitidae are still unclear. To address these problems, we inferred a tree on the basis of the sequences of mitochondrial genomes and of multiple nuclear genes using the maximum likelihood method. Our results strongly support the hypothesis that the Taxidiinae branched at first, followed by the branching of the Melinae. After that, Mustelinae diversified, and Lutrinae evolved within Mustelinae. With respect to the deep branchings in Musteloidea, the Ailuridae/Mephitidae monophyly tree and the Mephitidae-basal tree are indistinguishable in log-likelihood score, and this problem remains unresolved.

MyrSINEs: a novel SINE family in the anteater genomes

Recent rapid generation of genomic sequence data has allowed many researchers to perform comparative analyses in various mammalian species. However, characterization of transposable elements, such as short interspersed repetitive elements (SINEs), has not been reported for several mammalian groups. Because SINEs occupy a large portion of the mammalian genome, they are believed to have contributed to the constitution and diversification of the host genomes during evolution. In the present study, we characterized a novel SINE family in the anteater genomes and designated it the MyrSINE family. Typical SINEs consist of a tRNA-related, a tRNA-unrelated and an AT-rich (or poly-A) region. MyrSINEs have only tRNA-related and poly-A regions; they are included in a group called t-SINE. The tRNA-related regions of the MyrSINEs were found to be derived from tRNA(Gly). We demonstrate that the MyrSINE family can be classified into three subfamilies. Two of the MyrSINE subfamilies are distributed in the genomes of both giant anteater and tamandua, while the other is present only in the giant anteater. We discuss the evolutionary history of MyrSINEs and their relationship to the evolution of anteaters. We also speculate that the simple structure of t-SINEs may be a potential evolutionary source for the generation of the typical SINE structure.

Genomic expansion of the Bov-A2 retroposon relating to phylogeny and breed management

Bov-A2 is a retroposon that is widely distributed among the genomes of ruminants (e.g., cow, deer, giraffe, pronghorn, musk deer, and chevrotain). This retroposon is composed of two monomers, called Bov-A units, which are joined by a linker sequence. The structure and origin of Bov-A2 has been well characterized but a genome-level exploration of this retroposon has not been implemented. In this study we performed an extensive search for Bov-A2 using all available genome sequence data on Bos taurus. We found unique Bov-A2-derived sequences that were longer than Bov-A2 due to amplification of three to six Bov-A units arranged in tandem. Detailed analysis of these elongated Bov-A2-derived sequences revealed that they originated through unequal crossing-over of Bov-A2. We found a large number of these elongated Bov-A2-derived sequences in cattle genomes, indicating that unequal crossing-over of Bov-A2 occurred very frequently. We found that this type of elongation is not observed in wild bovine and is therefore specific to the domesticated cattle genome. Furthermore, at specific loci, the number of Bov-A units was also polymorphic between alleles, implying that the elongation of Bov-A units might have occurred very recently. For these reasons, we speculate that genomic instability in bovine genomes can lead to extensive unequal crossing-over of Bov-A2 and levels of polymorphism might be generated in part by repeated outbreeding.

Balaenoptera omurai is a newly discovered baleen whale that represents an ancient evolutionary lineage

Balaenoptera omurai, formerly classified as a small form of Bryde's whale, was recently reclassified as a new baleen whale species of the family Balaenopteridae. Although researchers have investigated the evolutionary history of Balaenopteridae and their relatives using molecular phylogenetic methods, the taxonomy of the ordinary Bryde's whale (Balaenoptera brydei) and small-form Bryde's whales (Balaenoptera edeni and B. omurai) remains unclear. We have used complete mtDNA sequences and short interspersed repetitive element (SINE) insertion patterns to construct the evolutionary history of both B. omurai and the taxonomically redefined species, B. edeni. The combined results demonstrate that B. omurai forms a monophyletic lineage with B. musculus, B. brydei, B. edeni and B. borealis and that B. omurai and B. musculus successively diverged from their common ancestor. In addition, we also showed that B. edeni constitutes a sister taxon to B. brydei. Our data suggest that B. omurai evolved as an ancient independent lineage that diverged much earlier than B. borealis, B. brydei and B. edeni, which were previously believed to be closely related to B. omurai.

Toothed whale monophyly reassessed by SINE insertion analysis: the absence of lineage sorting effects suggests a small population of a common ancestral species

Morphological data have indicated that toothed whales form a monophyletic group. However, research published in the last several years has made the issue of the monophyly or paraphyly of toothed whales a subject of debate. Our group previously characterized three independent loci in which SINE insertions were shared among dolphins and sperm whales, thus supporting the traditional, morphologically based hypothesis of toothed whale monophyly. Although in recent years a few additional molecular works proposed this topology, there is still skepticism over this monophyly from the view point of molecular systematics. When the phylogeny of rapidly radiated taxa is examined using the SINE method, it is important to consider the ascertainment bias that arises when choosing a particular taxon for SINE loci screening. To overcome this methodological problem specific to the SINE method, we examined all possible topologies among sperm whales, dolphins and baleen whales by extensively screening SINE loci from species of all three lineages. We characterized nine independent SINE loci from the genomes of sperm whales and dolphins, all of which cluster sperm whales and dolphins but exclude baleen whales. Furthermore, we characterized ten independent loci from baleen whales, all of which were amplified in a common ancestor of these whales. From these observations, we conclude that toothed whales form a monophyletic group and that no ancestral SINE polymorphisms hinder their phylogenetic assignment despite the short divergence times of the major lineages of extant whales during evolution. These results suggest that a small population of common ancestors of all toothed whales ultimately diverged into the lineages of sperm whales and dolphins.

New candidate species most closely related to penguins

The phylogenetic position of the order Spenisciformes in Aves remains unclear despite several independent analyses based on morphological and molecular data. To address this issue, we determined the complete mtDNA sequence of rockhopper penguins. The mitochondrial genome, excluding the region from the D-loop to 12SrRNA, was also sequenced for petrel, albatross, frigatebird, loon and grebe, which previous studies suggest are related to penguins. A maximum likelihood analysis of the phylogenetic placement of penguins with 23 birds, including 17 species whose mtDNA sequences were previously reported, suggested that storks are the closest extant relatives of penguins, with 78% and 56% bootstrap supports, depending on the choice of outgroup species. Thus, ciconiiform birds constitute new candidates as the closest extant relatives of penguins (previously proposed candidates were either gaviiform, podicipediform, or procellariiform birds). In addition to this new evidence, our analysis gave evidence to some of ambiguous relationships in the avian tree: our analysis supported a basal split between passerines and other neoavians within Neoaves, and rejected the monophyly of Falconiformes as well as that of loons and grebes.

Extensive mitochondrial gene arrangements in coleoid Cephalopoda and their phylogenetic implications

We determined the complete mitochondrial genomes of five cephalopods of the Subclass Coleoidea (Suborder Oegopsida: Watasenia scintillans, Todarodes pacificus, Suborder Myopsida: Sepioteuthis lessoniana, Order Sepiida: Sepia officinalis, and Order Octopoda: Octopus ocellatus) and used them to infer phylogenetic relationships. In our Maximum Likelihood (ML) tree, sepiids (cuttlefish) are at the most basal position of all decapodiformes, and oegopsids and myopsids form a monophyletic clade, thus supporting the traditional classification of the Order Teuthida. We detected extensive gene rearrangements in the mitochondrial genomes of broad cephalopod groups. It is likely that the arrangements of mitochondrial genes in Oegopsida and Sepiida were derived from those of Octopoda, which is thought to be the ancestral order, by entire gene duplication and random gene loss. Oegopsida in particular has undergone long-range gene duplications. We also found that the mitochondrial gene arrangement of Sepioteuthis lessoniana differs from that of Loligo bleekeri, although they belong to the same family. Analysis of both the phylogenetic tree and mitochondrial gene rearrangements of coleoid Cephalopoda suggests that each mitochondrial gene arrangement was acquired after the divergence of each lineage.

Metal levels in sugar cane (Saccharum spp.) samples from an area under the influence of a municipal landfill and a medical waste treatment system in Brazil

In July 2003, duplicated samples of roots, stems and leaves of sugar cane (Saccharum spp.) were collected in 25 points of an area under direct influence of the municipal landfill site (MLS) and medical waste treatment system (MWTS) of Ribeirao Preto, São Paulo, Brazil. Cadmium (Cd), chromium (Cr), copper (Cu), mercury (Hg), manganese (Mn), lead (Pb) and zinc (Zn) were determined by atomic absorption spectrophotometry. The following concentrations (mg/kg) were found in roots: Cd, 0.22+/-0.12; Cr, 64.3+/-48.7; Cu, 140.6+/-27.7; Hg, 0.04+/-0.02; Mn, 561.6+/-283.3; Pb, 7.9+/-2.1 and Zn, 177.4+/-64.9. For some metals, these levels are higher than the concentrations previously reported for different plants, reaching, in some cases, values that might be considered toxic for vegetables. Metal levels in stems were 80-90% of those found in roots, while the concentrations detected in leaves were significantly lower than those in roots. The present results suggest that MLS and MWTS activities might have been increasing metal concentrations in edible tissues of sugar cane grown in the area under their influence. Moreover, the traditional agricultural practices in the production of sugar cane could be also another determinant factor to reach the current metal levels. The results of this study indicate that sugar cane is a crop that is able to grow in areas where metals in soils are accumulated.

Phylogenetic Analysis of Diprotodontian Marsupials Based on Complete Mitochondrial Genomes

Australidelphia is the cohort, originally named by Szalay, of all Australian marsupials and the South American Dromiciops. A lot of mitochondria and nuclear genome studies support the hypothesis of a monophyly of Australidelphia, but some familial relationships in Australidelphia are still unclear. In particular, the familial relationships among the order Diprotodontia (koala, wombat, kangaroos and possums) are ambiguous. These Diprotodontian families are largely grouped into two suborders, Vombatiformes, which contains Phascolarctidae (koala) and Vombatidae (wombat), and Phalangerida, which contains Macropodidae, Potoroidae, Phalangeridae, Petauridae, Pseudocheiridae, Acrobatidae, Tarsipedidae and Burramyidae. Morphological evidence and some molecular analyses strongly support monophyly of the two families in Vombatiformes. The monophyly of Phalangerida as well as the phylogenetic relationships of families in Phalangerida remains uncertain, however, despite searches for morphological synapomorphy and mitochondrial DNA sequence analyses. Moreover, phylogenetic relationships among possum families (Phalangeridae, Petauridae, Pseudocheiridae, Acrobatidae, Tarsipedidae and Burramyidae) as well as a sister group of Macropodoidea (Macropodidae and Potoroidae) remain unclear. To evaluate familial relationships among Dromiciops and Australian marsupials as well as the familial relationships in Diprotodontia, we determined the complete mitochondrial sequence of six Diprotodontian species. We used Maximum Likelihood analyses with concatenated amino acid and codon sequences of 12 mitochondrial protein genomes. Our analysis of mitochondria amino acid sequence supports monophyly of Australian marsupials+Dromiciops and monophyly of Phalangerida. The close relatedness between Macropodidae and Phalangeridae is also weakly supported by our analysis.

The rise and fall of the CR1 subfamily in the lineage leading to penguins

The evolution of penguins has been investigated extensively, although inconclusively, by morphologists, biogeographers and molecular phylogeneticists. We investigated this issue using retroposon analysis of insertions of CR1, which is a member of the LINE (long interspersed element) family, in the genomes of penguins and penguin relatives. The retroposon method is a powerful tool for identifying monophyletic groups. Because retroposons often show different relative frequencies of retroposition during evolution, it is first necessary to identify a certain subgroup that was specifically active during the period when the species in question diverged. Hence, we systematically analyzed many CR1 members isolated from penguin and penguin-related genomes. These CR1s are divided into at least three distinct subgroups that share diagnostic nucleotide insertions and/or deletions, namely, penguin CR1 Sph I, Sph II type A and Sph II type B. The analysis of the inserted retroposons by PCR revealed that different CR1 subfamilies or types had amplified at different rates among different periods during penguin evolution. Namely, the penguin CR1 Sph I subfamily had higher rates of retroposition in a common ancestor of all orders examined in this study or at least in a common ancestor of all extant penguins, and the subfamily Sph II type A also had the same tendency. Therefore, these CR1 members can be used to elucidate the phylogenetic relationships of Sphenisciformes (penguins) among different avian orders. In contrast, the penguin CR1 Sph II type B subfamily had higher rates of retroposition just before and after the emergence of the extant genera in Spheniscidae, suggesting that they are useful for elucidating the intra-relationships among extant penguins. This is the first report for the characterization among the members of CR1 family in avian genomes excluding those of chickens. Hence, this work will be a cornerstone for elucidating the phylogenetic relationships in penguin evolution using the retroposon method.

Baleen Whale Phylogeny and a Past Extensive Radiation Event Revealed by SINE Insertion Analysis

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.