Integrating biogeographic and genetic approaches to investigate the history of bioluminescent colour alleles in the Jamaican click beetle, Pyrophorus plagiophthalamus

Bioluminescent colour in the Jamaican click beetle, Pyrophorus plagiophthalamus, is an ideal system for studies moving from gene to landscape to gain a holistic understanding of the molecular, ecological, and historical bases for adaptation. Previous studies have established the genetics of bioluminescent colour variation in the beetle to the level of the nucleotide base pair in the target gene luciferase. Three different luciferase colour alleles affecting ventral light organ colour [yellow-green (vYG), yellow (vYE), and orange (vOR)] were found segregating in P. plagiophthalamus populations. These alleles differ from each other in a number of replacement mutations (14 total), the majority of which (11) have a measurable effect on colour. Phylogenetic analysis revealed a long-term adaptive trend on Jamaica towards longer wavelength bioluminescence, culminating in the most recently derived vOR allele. Here, we further investigate the historical and geographic context of adaptive colour evolution by testing a vicariance model for the origins of the extant ventral light organ polymorphism: that the vOR allele arose and differentiated in an isolated deme on the east side of Jamaica before spreading westward. Comparisons of colour phenotypes, luciferase coding sequences, the third intron of the gene, mtDNA, and microsatellite data provided evidence for past population subdivision on Jamaica and ongoing gene flow, as has been found for other island endemics. However, the pattern of differentiation supported the allopatric divergence of vYG and vYE alleles. The vOR gene appears to have arisen relatively recently from a vYE precursor and postdates the period of major biogeographic isolation. We discuss the implications of the results for discerning ecological causation in the adaptive sequence from nucleotide to landscape to population change for bioluminescent colour.     

The Spectral Distribution of Firefly Light. II

The in vivo peak emission wavelengths of bioluminescence are reported for 15 species of American fireflies. A spectrophotometric study of the dorsal light organs of 155 specimens of the Jamaican firefly Pyrophorus plagiophthalamus showed three distinct color distributions with peak emission wavelengths at 550.1 ± 1.5 mµ, 556.8 ± 1.4 mµ, and 562.4 ± 1.0 mµ. Similar spectral measurements of 35 ventral light organs of the same insects gave peak emission wavelengths ranging from 547 through 594 mµ. This is a wider distribution than the total range of all 34 species of firefly studied to date. There was no obvious correlation between the colors of the ventral and dorsal light organs. It appears that P. plagiophthalamus is a special case in which the luciferase enzyme is not only different among members of the same species, but it may be different for the dorsal and ventral light organs in a single individual. A minimum of six different luciferase molecules for P. plagiophthalamus ventral light organs is proposed. The statistical precision in making these spectrophotometric measurements is discussed.     

Isolation and characterization of novel dinucleotide repeat microsatellites in the Jamaican click beetle Pyrophorus plagiophthalamus (Coleoptera: Elateridae)

We report 10 polymorphic microsatellite loci for the Jamaican click beetle, Pyrophorus plagiophthalamus. A survey of 36 individuals from three populations in Jamaica showed that these are highly variable, with three to 17 alleles per locus. Observed heterozygosity ranged from 0.250 to 0.917, and mean heterozygosity from 0.601 to 0.747. Most loci were in Hardy–Weinberg equilibrium, although excess of homozygotes was observed in four tests (out of 20), suggesting the possibility of null alleles. Significant linkage disequilibrium was observed for only one pair. These newly developed markers will be useful in understanding the population structure of click beetles in Jamaica, and in identifying possible selective factors responsible for bioluminescent colour variation.     

Functional divergence of gene duplicates through ectopic recombination

Gene duplication stimulates evolutionary innovation as the resulting paralogs acquire mutations that lead to sub‐ or neofunctionalization. A comprehensive in silico analysis of paralogs in Saccharomyces cerevisiae reveals that duplicates of cell‐surface and subtelomeric genes also undergo ectopic recombination, which leads to new chimaeric alleles. Mimicking such intergenic recombination events in the FLO (flocculation) family of cell‐surface genes shows that chimaeric FLO alleles confer different adhesion phenotypes than the parental genes. Our results indicate that intergenic recombination between paralogs can generate a large set of new alleles, thereby providing the raw material for evolutionary adaptation and innovation.     

Use of Controlled Luciferase Expression To Monitor Chemicals Affecting Protein Synthesis

In this article, we present a new bioluminescent test system for the screening of chemical compounds with an inhibitory effect on protein synthesis. The test is based on the measurement of real-time in vivo light production by Escherichia coli strains expressing different luciferase genes. The eukaryotic lucGR gene from Pyrophorus plagiophthalamus was found to be the best of three types of luciferase genes tested. Chemicals with known inhibitory effects on protein synthesis were used as test chemicals together with some general toxicants. The incubation of a test chemical with cells was performed either prior to or after the induction of protein synthesis, and the difference in the results of the two methods distinguishes the possible influence on protein synthesis from direct metabolic inhibition. Using lyophilized bacteria, the test is performed in less than an hour without any bacterial cultivation, which makes the test suitable for rapid and sensitive screening of chemicals or environmental samples. Compared with the standardized 50% inhibitory concentration calculation method of the bioluminescent cytotoxicity test, the more direct approach of calculation developed in this study proved to be more convenient than and as reliable as the standard method.     

Insertions, substitutions, and the origin of microsatellites

This paper uses data from the Human Gene Mutation Database to contrast two hypotheses for the origin of short DNA repeats: substitutions and insertions that duplicate adjacent sequences. Because substitutions are much more common than insertions, they are the dominant source of new 2-repeat loci. Insertions are rarer, but over 70% of the 2-4 base insertion mutations are duplications of adjacent sequences, and over half of these generate new repeat regions. Insertions contribute fewer new repeat loci than substitutions, but their relative importance increases rapidly with repeat number so that all new 4-5-repeat mutations come from insertions, as do all 3-repeat mutations of tetranucleotide repeats. This suggests that the process of repeat duplication that dominates microsatellite evolution at high repeat numbers is also important very early in microsatellite evolution. This result sheds light on the puzzle of the origin of short tandem repeats. It also suggests that most short insertion mutations derive from a slippage-like process during replication.     

Members of a dinoflagellate luciferase gene family differ in synonymous substitution rates.

Regulation and evolution of dinoflagellate luciferases are of particular interest since the enzyme is structurally unique and bioluminescence is under circadian control. In this study, three new members of the dinoflagellate luciferase gene family were identified and characterized from Pyrocystis lunula. These genes, lcfA, lcfB, and lcfC, also exhibit the unusual structure and organization previously reported for the luciferase gene of a related dinoflagellate, Lingulodinium polyedrum: three repeated domains, each encoding an active catalytic site, multiple gene copies, and tandem organization. The histidine residues involved in the pH regulation of L. polyedrum luciferase activity, and implicated in the regulation of flashing, are also fully conserved in P. lunula. The interspecific conservation between the individual luciferase domains of P. lunula and L. polyedrum is higher than among domains intramolecularly, indicating that this unique gene structure arose through duplication events that occurred prior to the divergence of these dinoflagellates. However, P. lunula luciferase genes differ from L. polyedrum in several respects, notably, the occurrence of an intron in one gene (lcfC), a 2.25-kb intergenic region connecting lcfA and lcfB, and, of particular interest, an invariant rate of synonymous (silent) substitutions along the repeat domains, in contrast to L. polyedrum luciferase, where the occurrence of synonymous substitutions is practically absent in the central region of the domains.     

Orthologous gene of beetle luciferase in non-luminous click beetle, Agrypnus binodulus (Elateridae), encodes a fatty acyl-CoA synthetase

A homologous gene of beetle luciferase, AbLL (Agrypnus binodulusluciferase-like gene) was isolated from a Japanese non-luminous click beetle, A. binodulus, and its gene product was characterized. The identity of amino acid sequence deduced from AbLL with the click beetle luciferase from the Jamaican luminous click beetle, Pyrophorus plagiophthalmus, is 55%, which is higher than that between click beetle luciferase and firefly luciferase (approximately 48%). Phylogenetic analysis indicated that AbLL places in a clade of beetle luciferases, suggesting that AbLL is an orthologous gene of beetle luciferase. The gene product of AbLL (AbLL) has medium- and long-chain fatty acyl-CoA synthetase activity, but not luciferase activity. The fatty acyl-CoA synthetic activity was slightly inhibited in the presence of beetle luciferin, suggesting that AbLL has poor affinity for beetle luciferin. By comparing the amino acid residues of the catalytic domains in beetle luciferases with AbLL, the key substitutions for the luminescence activity in beetle luciferase will be proposed.     

Luciferase from Fulgeochlizus bruchi (Coleoptera:Elateridae), a Brazilian click-beetle with a single abdominal lantern: molecular evolution, biological function and comparison with other click-beetle luciferases

Bioluminescent click-beetles emit a wide range of bioluminescence colors (λ(Max) = 534-594 nm) from thoracic and abdominal lanterns, which are used for courtship. Only the luciferases from Pyrophorus and Pyrearinus species were cloned and sequenced. The Brazilian Fulgeochlizus bruchi click-beetle, which inhabits the Central-west Cerrado (Savannas), is noteworthy because, differently from other click-beetles, the adult stage displays only a functional abdominal lantern, which produces a bright green bioluminescence for sexual attraction purposes, and lacks functional thoracic lanterns. We cloned the cDNA for the abdominal lantern luciferase of this species. Notably, the primary sequence of this luciferase showed slightly higher identity with the green emitting dorsal lantern luciferases of the Pyrophorus genus instead of the abdominal lanterns luciferases. This luciferase displays a blue-shifted spectrum (λ(Max) = 540 nm), which is pH-insensitive from pH 7.5 to 9.5 and undergoes a slight red shift and broadening above this pH; the lowest K(M) for luciferin among studied click-beetle luciferases, and the highest optimum pH (9.0) ever reported for a beetle luciferase. At pH 9.0, the K(M) for luciferin increases, showing a decrease of affinity for this substrate, despite the higher activity. The slow luminescence decay rate of F. bruchi luciferase in vitro reaction could be an adaptation of this luciferase for the long and sustained in vivo luminescence display of the click-beetle during the courtship, and could be useful for in vivo intracellular imaging.     

An alternative mechanism of bioluminescence color determination in firefly luciferase

Beetle luciferases (including those of the firefly) use the same luciferin substrate to naturally display light ranging in color from green (lambda(max) approximately 530 nm) to red (lambda(max) approximately 635 nm). In a recent communication, we reported (Branchini, B. R., Murtiashaw, M. H., Magyar, R. A., Portier, N. C., Ruggiero, M. C., and Stroh, J. G. (2002) J. Am. Chem. Soc. 124, 2112-2113) that the synthetic adenylate of firefly luciferin analogue D-5,5-dimethylluciferin was transformed into the emitter 5,5-dimethyloxyluciferin in bioluminescence reactions catalyzed by luciferases from Photinus pyralis and the click beetle Pyrophorus plagiophthalamus. 5,5-Dimethyloxyluciferin is constrained to exist in the keto form and fluoresces mainly in the red. However, bioluminescence spectra revealed that green light emission was produced by the firefly enzyme, and red light was observed with the click beetle protein. These results, augmented with steady-state kinetic studies, were taken as experimental support for mechanisms of firefly bioluminescence color that require only a single keto form of oxyluciferin. We report here the results of mutagenesis studies designed to determine the basis of the observed differences in bioluminescence color with the analogue adenylate. Mutants of P. pyralis luciferase putative active site residues Gly246 and Phe250, as well as corresponding click beetle residues Ala243 and Ser247 were constructed and characterized using bioluminescence emission spectroscopy and steady state kinetics with adenylate substrates. Based on an analysis of these and recently reported (Branchini, B. R., Southworth, T. L., Murtiashaw, M. H., Boije, H., and Fleet, S. E. (2003) Biochemistry 42, 10429-10436) data, we have developed an alternative mechanism of bioluminescence color. The basis of the mechanism is that luciferase modulates emission color by controlling the resonance-based charge delocalization of the anionic keto form of the oxyluciferin excited state.     

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.     

The evolution of the adenylate-forming protein family in beetles: multiple luciferase gene paralogues in fireflies and glow-worms

Bioluminescence in beetles is dependent upon the enzyme luciferase. It has been hypothesised luciferase evolved from a fatty acyl-CoA synthetase gene deriving a novel bioluminescent function (neofunctionalization) after a gene duplication event. We evaluated this hypothesis within a phylogenetic framework using independent evidence obtained from the genome of Tribolium castaneum, published luciferase genes and novel luciferase and luciferase-like sequences. This phylogenetic study provides evidence for a large gene family of luciferase and luciferase-like paralogues in bioluminescent and non-bioluminescent beetles. All luciferase sequences formed a clade supporting a protoluciferase existing prior to the divergence of the Lampyridae, Elateridae and Phengodidae (Elateroidea). Multiple luciferase genes were identified from members of the Photurinae and the Luciolinae indicating complex gene duplication events within lampyrid genomes. The majority of luciferase residues were identified to be under purifying selection as opposed to positive selection. We conclude that beetle luciferase may have arisen from a process of subfunctionalization as opposed to neofunctionalization early on in the evolution of the Elateroidea.     

Gene duplication at an isocitrate dehydrogenase locus in Scaphiopus

Spadefoot toads of the subgenus Scaphiopus have two isocitrate dehydrogenase loci, with no intergenic interaction between them. Toads of the subgenus Spea have three Idh loci, with intergenic enzymes formed between two of them, providing strong evidence for their homology and the origin of one through a duplication process. The Idh phenotype of interspecific hybrids is consistent with the theory of a gene duplication.     

Spectral correspondence between visual spectral sensitivity and bioluminescence emission spectra in the click beetle Pyrophorus punctatissimus (Coleoptera: Elateridae)

The presence of two spectral mechanisms, near-ultraviolet and green (lambda(max)=545nm), is strongly suggested by electroretinographic visual spectral sensitivity curves obtained under dark and red chromatic adaptation conditions in the compound eyes of the click beetle Pyrophorus punctatissimus. The bioluminescence emission of the dorsal prothoracic lanterns is deep green (lambda(max)=543nm) and that of the ventral abdominal lantern is lime green (lambda(max)=556nm) in colour in P. punctatissimus. A broad green visual receptor would detect both deep green and lime green bioluminescent optical signals.     

Genome-level evolution of resistance genes in Arabidopsis thaliana

Pathogen resistance genes represent some of the most abundant and diverse gene families found within plant genomes. However, evolutionary mechanisms generating resistance gene diversity at the genome level are not well understood. We used the complete Arabidopsis thaliana genome sequence to show that most duplication of individual NBS-LRR sequences occurs at close physical proximity to the parent sequence and generates clusters of closely related NBS-LRR sequences. Deploying the statistical strength of phylogeographic approaches and using chromosomal location as a proxy for spatial location, we show that apparent duplication of NBS-LRR genes to ectopic chromosomal locations is largely the consequence of segmental chromosome duplication and rearrangement, rather than the independent duplication of individual sequences. Although accounting for a smaller fraction of NBS-LRR gene duplications, segmental chromosome duplication and rearrangement events have a large impact on the evolution of this multigene family. Intergenic exchange is dramatically lower between NBS-LRR sequences located in different chromosome regions as compared to exchange between sequences within the same chromosome region. Consequently, once translocated to new chromosome locations, NBS-LRR gene copies have a greater likelihood of escaping intergenic exchange and adopting new functions than do gene copies located within the same chromosomal region. We propose an evolutionary model that relates processes of genome evolution to mechanisms of evolution for the large, diverse, NBS-LRR gene family.     

Gene duplication as a major force in evolution

Gene duplication is an important mechanism for acquiring new genes and creating genetic novelty in organisms. Many new gene functions have evolved through gene duplication and it has contributed tremendously to the evolution of developmental programmes in various organisms. Gene duplication can result from unequal crossing over, retroposition or chromosomal (or genome) duplication. Understanding the mechanisms that generate duplicate gene copies and the subsequent dynamics among gene duplicates is vital because these investigations shed light on localized and genomewide aspects of evolutionary forces shaping intra-specific and inter-specific genome contents, evolutionary relationships, and interactions. Based on whole-genome analysis of Arabidopsis thaliana, there is compelling evidence that angiosperms underwent two whole-genome duplication events early during their evolutionary history. Recent studies have shown that these events were crucial for creation of many important developmental and regulatory genes found in extant angiosperm genomes. Recent studies also provide strong indications that even yeast (Saccharomyces cerevisiae), with its compact genome, is in fact an ancient tetraploid. Gene duplication can provide new genetic material for mutation, drift and selection to act upon, the result of which is specialized or new gene functions. Without gene duplication the plasticity of a genome or species in adapting to changing environments would be severely limited. Whether a duplicate is retained depends upon its function, its mode of duplication, (i.e. whether it was duplicated during a whole-genome duplication event), the species in which it occurs, and its expression rate. The exaptation of preexisting secondary functions is an important feature in gene evolution, just as it is in morphological evolution.     

Evolutionary dynamics of duplicated genes in plants

Gene duplication, arising from region-specific duplication or genome-wide polyploidization, is a prominent feature in plant genome evolution. Understanding the mechanisms generating duplicate gene copies and the subsequent dynamics among gene duplicates is vital because these investigations shed light on regional and genome-wide aspects of evolutionary forces shaping intra- and interspecific genome contents, evolutionary relationships, and interactions. This review discusses recent gene duplication analyses in plants, focusing on the molecular and evolutionary dynamics occurring at three different timescales following duplication: (1). initial establishment and persistence of cytotypes, (2). interactions among duplicate gene copies, and (3). longer term differentiation between duplicated genes. These relative time points are presented in terms of their potential adaptive significance and impact on plant evolutionary genomics research.     

Synorth: exploring the evolution of synteny and long-range regulatory interactions in vertebrate genomes

Genomic regulatory blocks are chromosomal regions spanned by long clusters of highly conserved noncoding elements devoted to long-range regulation of developmental genes, often immobilizing other, unrelated genes into long-lasting syntenic arrangements. Synorth is a web resource for exploring and categorizing the syntenic relationships in genomic regulatory blocks across multiple genomes, tracing their evolutionary fate after teleost whole genome duplication at the level of genomic regulatory block loci, individual genes, and their phylogenetic context.     

Asymmetric Functional Divergence of Young, Dispersed Gene Duplicates in Arabidopsis thaliana

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K _sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno’s prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene duplicates.