Evolution of bioluminescence in marine planktonic copepods

Copepods are the dominant taxa in zooplankton communities of the ocean worldwide. Although bioluminescence of certain copepods has been known for more than a 100 years, there is very limited information about the structure and evolutionary history of copepod luciferase genes. Here, we report the cDNA sequences of 11 copepod luciferases isolated from the superfamily Augaptiloidea in the order Calanoida. Highly conserved amino acid residues in two similar repeat sequences were confirmed by the multiple alignment of all known copepod luciferases. Copepod luciferases were classified into two groups of Metridinidae and Heterorhabdidae/Lucicutiidae families based on phylogenetic analyses, with confirmation of the interrelationships within the Calanoida using 18S ribosomal DNA sequences. The large diversity in the specific activity of planktonic homogenates and copepod luciferases that we were able to express in mammalian cultured cells illustrates the importance of bioluminescence as a protective function against predators. We also discuss the relationship between the evolution of copepod bioluminescence and the aspects of their ecological characteristics, such as swimming activity and vertical habitat.     

Firefly luciferase genes from the subfamilies Psilocladinae and Ototretinae (Lampyridae, Coleoptera)

Firefly luciferase genes have been isolated from approximately 20 species of Lampyrinae, Luciolinae, and Photurinae. These are mostly nocturnal luminescent species that use light signals for sexual communication. In this study, we isolated three cDNAs for firefly luciferase from Psilocladinae (Cyphonocerus ruficollis) and Ototretinae (Drilaster axillaris and Stenocladius azumai), which are diurnal non-luminescent or weakly luminescent species that may use pheromones for communication. The amino acid sequences deduced from the three cDNAs showed 81-89% identities to each other and 60-81% identities with known firefly luciferases. The three purified recombinant proteins showed luminescence and fatty acyl-CoA synthetic activities, as observed in other firefly luciferases. The emission maxima by the three firefly luciferases (λmax, 545-546 nm) were shorter than those by known luciferases from the nocturnal fireflies (λmax, 550-568 nm). These results suggest that the primary structures and enzymatic properties of luciferases are conserved in Lampyridae, but the luminescence colors were red-shifted in nocturnal species compared to diurnal species.     

Identification of two catalytic domains in a luciferase secreted by the copepod Gaussia princeps

Gaussia luciferase secreted by the copepod Gaussia princeps catalyzes the oxidation of coelenterazine to produce blue light. The primary structure of Gaussia luciferase deduced from the cDNA sequence shows two repeat sequences of 71 amino acid residues, suggesting the luciferase consists of two structural domains. Two domains in Gaussia luciferase were expressed independently in Escherichia coli cells, purified and characterized. We found that both domains have luminescence activity with coelenterazine, and the catalytic properties including luminescence spectrum, optimal pH, substrate specificity and luminescence stimulation by halogen ions (Cl⁻, Br⁻ and I⁻) are identical to intact Gaussia luciferase. Thus, Gaussia luciferase has two catalytic domains for the luminescence reaction.     

Functional conversion of fatty acyl-CoA synthetase to firefly luciferase by site-directed mutagenesis: A key substitution responsible for luminescence activity

We demonstrated that firefly luciferase has a catalytic function of fatty acyl-CoA synthesis [Oba, Y., Ojika, M. and Inouye, S. (2003) Firefly luciferase is a bifunctional enzyme: ATP-dependent monooxygenase and a long chain fatty acyl-CoA synthetase. FEBS Lett. 540, 251-254] and proposed that the evolutionary origin of beetle luciferase is a fatty acyl-CoA synthetase (FACS) in insect. In this study, we performed the functional conversion of FACS to luciferase by replacing a single amino acid to serine. This serine residue is conserved in luciferases and possibly interacts with luciferin. The mutants of FACSs in non-luminous click beetle Agrypnus binodulus (AbLL) and Drosophilamelanogaster (CG6178) gave luminescence enhancement, suggesting that the serine residue is a key substitution responsible for luminescence activity.     

Interconversion of ketoprofen recognition in firefly luciferase-catalyzed enantioselective thioesterification reaction using from Pylocoeria miyako (PmL) and Hotaria parvura (HpL) just by mutating two amino acid residues

We identified the critical amino acid residues for substrate recognition using two firefly luciferases from Pylocoeria miyako (PmL) and Hotaria parvura (HpL), as these two luciferase enzymes exhibit different activities toward ketoprofen. Specifically, PmL can catalyze the apparent enantioselective thioesterification reaction, while HpL cannot. By comparing the amino acid sequences around the active site, we identified two residues (I350 and M397 in PmL and F351 and S398 in HpL) that were different between the two enzymes, and the replacement of these amino acids resulted in changing the ketoprofen recognition pattern. The inactive HpL was converted to the active enzyme toward ketoprofen and vice versa for PmL. These residues also affected the enantioselectivity toward ketoprofen; however, the bioluminescent color was not affected. In addition, using molecular dynamics calculations, the replacement of these two amino acids induced changes in the state of hydrogen bonding between ketoprofen and the S349 side chain through the active site water. As S349 is not considered to influence color tuning, these changes specifically caused the differences in ketoprofen recognition in the enzyme.     

Positive selection in AvrP4 avirulence gene homologues across the genus Melampsora

Pathogen genes involved in interactions with their plant hosts are expected to evolve under positive Darwinian selection or balancing selection. In this study a single copy avirulence gene, AvrP4, in the plant pathogen Melampsora lini, was used to investigate the evolution of such a gene across species. Partial translation elongation factor 1-alpha sequences were obtained to establish phylogenetic relationships among the Melampsora species. We amplified AvrP4 homologues from species pathogenic on hosts from different plant families and orders, across the inferred phylogeny. Translations of the AvrP4 sequences revealed a predicted signal peptide and towards the C-terminus of the protein, six identically spaced cysteines were identified in all sequences. Maximum likelihood analysis of synonymous versus non-synonymous substitution rates indicated that positive selection played a role in the evolution of the gene during the diversification of the genus. Fourteen codons under significant positive selection reside in the C-terminal 28 amino acid region, suggesting that this region interacts with host molecules in most sequenced accessions. Selection pressures on the gene may be either due to the pathogenicity or avirulence function of the gene or both.     

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.     

Zonal distribution and seasonal vertical migration of copepod assemblages in the Scotia Sea

Large, biomass-dominant Southern Ocean copepod species have been much studied, but small and mesopelagic species also play major rôles in these ecosystems. However, little is known of some basic aspects of their ecology. To address this, the abundances of 23 copepod species and genera were analysed from 72 stations sampled during the Discovery Expeditions in the 1920s to 1950s. Stratified net samples, usually to a depth of 1000?m, provided year-round coverage in the Scotia Sea from the Subantarctic Front to the Weddell-Scotia Confluence. Small copepods (Microcalanus pygmaeus, Ctenocalanus spp., Oncaea spp. and Oithona spp.) formed ～75% of total copepod abundance in the top 1000?m across all major zones. Oithona spp. composed ～40% of copepod numbers in the Polar Front area and to its south: further north their importance declined. All mesopelagic taxa except for the warmer-water species Metridia lucens and Pleuromamma robusta, extended throughout the entire study area, with smaller regional differences than for the shallower-living species. The species showed a continuum of temperature ranges, and there was no evidence that the Polar Front was a major biogeographic boundary to their distribution. Indeed, several important species, including Oithona spp. (mainly Oithona similis), Ctenocalanus spp., Metridia lucens and Rhincalanus gigas reached maximum numbers in this area. Total copepod abundance was thus higher in the vicinity of the Polar Front than in any other region. Only two copepod families made pronounced seasonal vertical migrations: Eucalaniidae (Eucalanus longiceps and R. gigas) and Calaniidae (Neocalanus tonsus, Calanoides acutus, Calanus simillimus and Calanus propinquus). Some evidence for a winter descent was found for Ctenocalanus spp. and some deeper-living groups: Euchaeta spp. and the Metridiidae, although their migrations were not so great as for the eucalanids and calanids.     

Speciation slowing down in widespread and long-living tree taxa: insights from the tropical timber tree genus Milicia (Moraceae)

The long generation time and large effective size of widespread forest tree species can result in slow evolutionary rate and incomplete lineage sorting, complicating species delimitation. We addressed this issue with the African timber tree genus Milicia that comprises two morphologically similar and often confounded species: M. excelsa, widespread from West to East Africa, and M. regia, endemic to West Africa. We combined information from nuclear microsatellites (nSSRs), nuclear and plastid DNA sequences, and morphological systematics to identify significant evolutionary units and infer their evolutionary and biogeographical history. We detected five geographically coherent genetic clusters using nSSRs and three levels of genetic differentiation. First, one West African cluster matched perfectly with the morphospecies M. regia that formed a monophyletic clade at both DNA sequences. Second, a West African M. excelsa cluster formed a monophyletic group at plastid DNA and was more related to M. regia than to Central African M. excelsa, but shared many haplotypes with the latter at nuclear DNA. Third, three Central African clusters appeared little differentiated and shared most of their haplotypes. Although gene tree paraphyly could suggest a single species in Milicia following the phylogenetic species concept, the existence of mutual haplotypic exclusivity and nonadmixed genetic clusters in the contact area of the two taxa indicate strong reproductive isolation and, thus, two species following the biological species concept. Molecular dating of the first divergence events showed that speciation in Milicia is ancient (Tertiary), indicating that long-living tree taxa exhibiting genetic speciation may remain similar morphologically.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Adaptive evolution of tight junction protein claudin-14 in echolocating whales

Toothed whales and bats have independently evolved specialized ultrasonic hearing for echolocation. Recent findings have suggested that several genes including Prestin, Tmc1, Pjvk and KCNQ4 appear to have undergone molecular adaptations associated with the evolution of this ultrasonic hearing in mammals. Here we studied the hearing gene Cldn14, which encodes the claudin-14 protein and is a member of tight junction proteins that functions in the organ of Corti in the inner ear to maintain a cationic gradient between endolymph and perilymph. Particular mutations in human claudin-14 give rise to non-syndromic deafness, suggesting an essential role in hearing. Our results uncovered two bursts of positive selection, one in the ancestral branch of all toothed whales and a second in the branch leading to the delphinid, phocoenid and ziphiid whales. These two branches are the same as those previously reported to show positive selection in the Prestin gene. Furthermore, as with Prestin, the estimated hearing frequencies of whales significantly correlate with numbers of branch-wise non-synonymous substitutions in Cldn14, but not with synonymous changes. However, in contrast to Prestin, we found no evidence of positive selection in bats. Our findings from Cldn14, and comparisons with Prestin, strongly implicate multiple loci in the acquisition of echolocation in cetaceans, but also highlight possible differences in the evolutionary route to echolocation taken by whales and bats.     

Phylogenetic Relationships Among Selected Heteroderoidea Based on 18S and ITS Ribosomal DNA

In a study of relationships among selected cyst-forming and noncyst-forming species of Heteroderoidea, combined sequences comprised of DNA from part of the conserved 18S ribosomal RNA gene (rDNA) plus the complete ITS rDNA segment were more similar to analyses based on the ITS data alone than to analyses based on the 18S data alone. One of the two noncyst-forming species, Ekphymatodera thomasoni, grouped with cyst-forming species of Heteroderoidea. Bilobodera flexa, also a noncyst-forming species, was separated from all the other taxa by a long branch. Afenestrata koreana, with a weakly sclerotized cyst, grouped closely with H. bifenestra. These observations suggest that phylogenetic analyses using molecular data may aid in our understanding of the evolution of cyst formation in nematodes, including the possibility of secondary loss. The usefulness of molecular phylogenetic analyses in nematodes may depend more on the particular selection of taxa than on mere addition of data from additional genes.     

Description of Caurinus tlagu,new species,from Prince of Wales Island,Alaska (Mecoptera,Boreidae, Caurininae)

A new species of the cryptic, minute, wingless, and enigmatic taxon Caurinus, and the second for the subfamily Caurininae,is described from Prince of Wales Island in the Alexander Archipelago, Alaska. It is distinguished from its only congener, Caurinus dectes Russell, 1979b, which occurs 1,059 km southeast in Oregon and Washington, based on external morphology and sequences of the mitochondrial gene cytochrome oxidase II. These two species are probably evolutionary relicts – the only known members of a clade dating to the Late Jurassic or older.     

Identification of MicroRNAs in Helicoverpa armigera and Spodoptera litura Based on Deep Sequencing and Homology Analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.     

The making of a pest: the evolution of a fruit-penetrating ovipositor in Drosophila suzukii and related species

Evolutionary innovation can allow a species access to a new ecological niche, potentially reducing competition with closely related species. While the vast majority of Drosophila flies feed on rotting fruit and other decaying matter, and are harmless to human activity, Drosophila suzukii, which has a morphologically modified ovipositor, is capable of colonizing live fruit that is still in the process of ripening, causing massive agricultural damage. Here, we conducted the first comparative analysis of this species and its close relatives, analysing both ovipositor structure and fruit susceptibility. We found that the ovipositor of the species most closely related to D. suzukii, Drosophila subpulchrella, has a similar number of enlarged, evolutionarily derived bristles, but a notably different overall shape. Like D. suzukii, D. subpulchrella flies are capable of puncturing the skin of raspberries and cherries, but we found no evidence that they could penetrate the thicker skin of two varieties of grapes. More distantly related species, one of which has previously been mistaken for D. suzukii, have blunt ovipositors with small bristles. While they did not penetrate fruit skin in any of the assays, they readily colonized fruit interiors where the skin was broken. Our results suggest that considering evolutionary context may be beneficial to the management of invasive species.     

The mitochondrial genomes of Amphiascoides atopus and Schizopera knabeni (Harpacticoida: Miraciidae) reveal similarities between the copepod orders Harpacticoida and Poecilostomatoida

Members of subclass Copepoda are abundant, diverse, and—as a result of their variety of ecological roles in marine and freshwater environments—important, but their phylogenetic interrelationships are unclear. Recent studies of arthropods have used gene arrangements in the mitochondrial (mt) genome to infer phylogenies, but for copepods, only seven complete mt genomes have been published. These data revealed several within-order and few among-order similarities. To increase the data available for comparisons, we sequenced the complete mt genome (13,831 base pairs) of Amphiascoides atopus and 10,649 base pairs of the mt genome of Schizopera knabeni (both in the family Miraciidae of the order Harpacticoida). Comparison of our data to those for Tigriopus japonicus (family Harpacticidae, order Harpacticoida) revealed similarities in gene arrangement among these three species that were consistent with those found within and among families of other copepod orders. Comparison of the mt genomes of our species with those known from other copepod orders revealed the arrangement of mt genes of our Harpacticoida species to be more similar to that of Sinergasilus polycolpus (order Poecilostomatoida) than to that of T. japonicus. The similarities between S. polycolpus and our species are the first to be noted across the boundaries of copepod orders and support the possibility that mt-gene arrangement might be used to infer copepod phylogenies. We also found that our two species had extremely truncated transfer RNAs and that gene overlaps occurred much more frequently than has been reported for other copepod mt genomes.