Phylogenetic analysis of the fast-evolving SuUR gene in insects

Different genomic regions replicate at a different fixed time during the S phase. Late-replicating sequences are often underreplicated in the Drosophila salivary-gland polytene chromosomes. The SuUR gene, whose mutation changes the replication time of late-replicating regions in salivary-gland cells, has been identified in Drosophila melanogaster. The SUUR protein lacks homologs by a BLAST search, and only moderate similarity is observed between its N-terminal part and chromatin-remodeling proteins of the SWI2/SNF2 family. The gene and the protein were analyzed in insects. Orthologs of the SuUR gene were found in all annotated Drosophila species. The number of amino acid substitutions in the SUUR protein proved to be extremely high, corresponding to that of fast-evolving genes. Orthologs with low homology were found in mosquitoes Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus. No orthologs of the SuUR gene were detected beyond Diptera.     

QTL analysis of a rapidly evolving speciation phenotype in the Hawaiian cricket Laupala

In mate recognition systems, the functional necessity to coordinate traits involved in sexual communication should result in reduced pairing potential for new variants outside the distribution of common reproductive signals. Yet, many closely related, sexual species differ in mate recognition traits, suggesting that directional selection influences the divergence of mate recognition systems. Species of the endemic Hawaiian cricket genus Laupala are morphologically and ecologically cryptic, although both male calling song and female acoustic preference have diverged rapidly between closely related species. These mate recognition phenotypes are therefore often likely to be speciation phenotypes, i.e. traits whose divergence contributes, directly or indirectly, to a reduction of gene flow during speciation, given their frequent participation in early lineage divergence. We conducted a replicated, quantitative trait loci (QTL) mapping study of the genetic basis of differences in male calling song between two closely related species, Laupala paranigra and Laupala kohalensis, allowing us to examine the genetic basis of traits involved in rapid speciation. We found statistical support for eight QTL in one replicate, with at least four of these QTL mapping to the same regions in a second replicate. QTL effects ranged between 3.0% and 10.7% of the difference in pulse rate between L. paranigra and L. kohalensis, and are thus of moderate to small effect. All QTL identified show directional effects consistent with the hypothesis of directional selection. Thus, we conclude that rapid speciation can occur under the influence of many genes of moderate to small effect. This study implicates the role of directional selection in the divergence of mate recognition traits and speciation the Hawaiian cricket genus Laupala.     

Selection on rapidly evolving proteins in the Arabidopsis genome

Genes that have undergone positive or diversifying selection are likely to be associated with adaptive divergence between species. One indicator of adaptive selection at the molecular level is an excess of amino acid replacement fixed differences per replacement site relative to the number of synonymous fixed differences per synonymous site (omega = K(a)/K(s)). We used an evolutionary expressed sequence tag (EST) approach to estimate the distribution of omega among 304 orthologous loci between Arabidopsis thaliana and A. lyrata to identify genes potentially involved in the adaptive divergence between these two Brassicaceae species. We find that 14 of 304 genes (approximately 5%) have an estimated omega > 1 and are candidates for genes with increased selection intensities. Molecular population genetic analyses of 6 of these rapidly evolving protein loci indicate that, despite their high levels of between-species nonsynonymous divergence, these genes do not have elevated levels of intraspecific replacement polymorphisms compared to previously studied genes. A hierarchical Bayesian analysis of protein-coding region evolution within and between species also indicates that the selection intensities of these genes are elevated compared to previously studied A. thaliana nuclear loci.     

Phylogenetic analysis of the isopenicillin-N-synthetase horizontal gene transfer

A phylogenetic study of the isopenicillin-N-synthetase (IPNS) gene sequence from prokaryotic and lower eukaryotic producers of β-lactam antibiotics by means of a maximum-likelihood approach has been carried out. After performing an extensive search, rather than invoking a global molecular clock, the results obtained are best explained by a model with three rates of evolution. Grouped in decreasing order, these correspond toA. nidulans and then to the rest of the eukaryotes and prokaryotes, respectively. The estimated branching date between prokaryotic and fungal IPNS sequences (852 ±106 MY) strongly supports the hypothesis that the IPNS gene was horizontally transferred from bacterial β-lactam producers to filamentous fungi. Correspondence to: A. Moya     

The rapidly evolving field of plant centromeres

Meiotic and mitotic chromosome segregation are highly conserved in eukaryotic organisms, yet centromeres--the chromosomal sites that mediate segregation--evolve extremely rapidly. Plant centromeres have DNA elements that are shared across species, yet they diverge rapidly through large- and small-scale changes. Over evolutionary time-scales, centromeres migrate to non-centromeric regions and, in plants, heterochromatic knobs can acquire centromere activity. Discerning the functional significance of these changes will require comparative analyses of closely related species. Combined with functional assays, continued efforts in plant genomics will uncover key DNA elements that allow centromeres to retain their role in chromosome segregation while allowing rapid evolution.     

Vsx1, a rapidly evolving paired-like homeobox gene expressed in cone bipolar cells.

The paired-like homeodomain (HD) protein Chx10 is distinguished by the presence of the CVC domain, a conserved 56 amino acid sequence C-terminal to the HD. In mammals, Chx10 is essential both for the proliferation of retinal progenitor cells and for the formation or survival of retinal bipolar interneurons. We describe the cloning and characterization of a mouse Chx10 homologue, Vsx1; phylogenetic analysis suggests that Vsx1 and its putative vertebrate orthologues have evolved rapidly. Vsx1 expression in the adult is predominantly retinal. Whereas Chx10 is expressed both in retinal progenitors in the developing eye and apparently in all bipolar cells of the mature retina, Vsx1 expression is first detected in the eye at postnatal day 5, where it is restricted to cone bipolar cells.     

Influence of the SuUR gene on intercalary heterochromatin in Drosophila melanogaster polytene chromosomes

Salivary gland polytene chromosomes of Drosophila melanogaster have a reproducible set of intercalary heterochromatin (IH) sites, characterized by late DNA replication, underreplicated DNA, breaks and frequent ectopic contacts. The SuUR mutation has been shown to suppress underreplication, and wild-type SuUR protein is found at late-replicating IH sites and in pericentric heterochromatin. Here we show that the SuUR gene influences all four IH features. The SuUR mutation leads to earlier completion of DNA replication. Using transgenic strains with two, four or six additional SuUR(+) doses (4-8xSuUR(+)) we show that wild-type SuUR is an enhancer of DNA underreplication, causing many late-replicating sites to become underreplicated. We map the underreplication sites and show that their number increases from 58 in normal strains (2xSuUR(+)) to 161 in 4-8xSuUR(+) strains. In one of these new sites (1AB) DNA polytenization decreases from 100% in the wild type to 51%-85% in the 4xSuUR (+) strain. In the 4xSuUR(+) strain, 60% of the weak points coincide with the localization of Polycomb group (PcG) proteins. At the IH region 89E1-4 (the Bithorax complex), a typical underreplication site, the degree of underreplication increases with four doses of SuUR(+) but the extent of the underreplicated region is the same as in wild type and corresponds to the region containing PcG binding sites. We conclude that the polytene chromosome regions known as IH are binding sites for SuUR protein and in many cases PcG silencing proteins. We propose that these stable silenced regions are late replicated and, in the presence of SuUR protein, become underreplicated.     

Phylogenetic discordance at the species boundary: comparative gene genealogies among rapidly radiating Heliconius butterflies

Recent adaptive radiations provide excellent model systems for understanding speciation, but rapid diversification can cause problems for phylogenetic inference. Here we use gene genealogies to investigate the phylogeny of recent speciation in the heliconiine butterflies. We sequenced three gene regions, intron 3 ( approximately 550 bp) of sex-linked triose-phosphate isomerase (Tpi), intron 3 ( approximately 450 bp) of autosomal mannose-phosphate isomerase (Mpi), and 1,603 bp of mitochondrial cytochrome oxidase subunits I and II (COI and COII), for 37 individuals from 25 species of Heliconius and related genera. The nuclear intron sequences evolved at rates similar to those of mitochondrial coding sequences, but the phylogenetic utility of introns was restricted to closely related geographic populations and species due to high levels of indel variation. For two sister species pairs, Heliconius erato-Heliconius himera and Heliconius melpomene-Heliconius cydno, there was highly significant discordance between the three genes. At mtDNA and Tpi, the hypotheses of reciprocal monophyly and paraphyly of at least one species with respect to its sister could not be distinguished. In contrast alleles sampled from the third locus, Mpi, showed polyphyletic relationships between both species pairs. In all cases, recent coalescence of mtDNA lineages within species suggests that polyphyly of nuclear genes is not unexpected. In addition, very similar alleles were shared between melpomene and cydno, implying recent gene flow. Our finding of discordant genealogies between genes is consistent with models of adaptive speciation with ongoing gene flow and highlights the need for multiple locus comparisons to resolve phylogeny among closely related species.     

Phylogenetic and evolutionary analysis of the PLUNC gene family

The PLUNC family of human proteins are candidate host defense proteins expressed in the upper airways. The family subdivides into short (SPLUNC) and long (LPLUNC) proteins, which contain domains predicted to be structurally similar to one or both of the domains of bactericidal/permeability-increasing protein (BPI), respectively. In this article we use analysis of the human, mouse, and rat genomes and other sequence data to examine the relationships between the PLUNC family proteins from humans and other species, and between these proteins and members of the BPI family. We show that PLUNC family clusters exist in the mouse and rat, with the most significant diversification in the locus occurring for the short PLUNC family proteins. Clear orthologous relationships are established for the majority of the proteins, and ambiguities are identified. Completion of the prediction of the LPLUNC4 proteins reveals that these proteins contain approximately a 150-residue insertion encoded by an additional exon. This insertion, which is predicted to be largely unstructured, replaces the structure homologous to the 40s hairpin of BPI. We show that the exon encoding this region is anomalously variable in size across the LPLUNC proteins, suggesting that this region is key to functional specificity. We further show that the mouse and human PLUNC family orthologs are evolving rapidly, which supports the hypothesis that these proteins are involved in host defense. Intriguingly, this rapid evolution between the human and mouse sequences is replaced by intense purifying selection in a large portion of the N-terminal domain of LPLUNC4. Our data provide a basis for future functional studies of this novel protein family.     

Phylogenetic analysis of the 5-aminolevulinate synthase gene.

The evolution of 5-aminolevulinate synthase (ALS) was studied by acquiring sequence data and generating phylogenetic trees. Gene sequences were already available for a variety of vertebrates (which have both a housekeeping and an erythroid form of the gene), fungi, alpha-proteobacteria, and one protist and one protostome. In order to generate representative trees, ALS sequence data were acquired from various deuterostomes and protostomes. The species and tissues selected for study were beluga whale liver, hagfish blood, sea urchin gonadal tissue, cuttlefish hepatopancreas, horseshoe crab hepatopancreas, and bloodworm blood. The new sequences and those previously published were examined for the presence of heme-regulatory motifs (HRMs) and iron-responsive elements (IREs). The HRMs are present in almost all eukaryotic species, which suggests their fundamental role in the regulation of ALS. The IREs are present in all vertebrate erythroid forms of ALS, which indicates that in those animals, expression of the erythroid form of the enzyme and, hence, hemoglobin production can be influenced by the intracellular content of iron. The new sequences were aligned with previously reported ALS sequences, and phylogenetic analyses were performed. The resulting trees provided evidence regarding the timing of the gene duplication event that led to the two forms of the ALS gene in vertebrates. It appears that the housekeeping and erythroid forms of ALS probably arose before the divergence of hagfish from the deuterostome line leading to the vertebrates. The data also add to the evidence indicating that alpha-proteobacteria are the nearest contemporary relatives of mitochondria.     

Phylogenetic distribution of TTAGG telomeric repeats in insects.

We examined the presence of TTAGG telomeric repeats in 22 species from 20 insect orders with no or inconclusive information on the telomere composition by single-primer polymerase chain reaction with (TTAGG)6 primers, Southern hybridization of genomic DNAs, and fluorescence in situ hybridization of chromosomes with (TTAGG)n probes. The (TTAGG)n sequence was present in 15 species and absent in 7 species. In a compilation of new and published data, we combined the distribution of (TTAGG)n telomere motif with the insect phylogenetic tree. The pattern of phylogenetic distribution of the TTAGG repeats clearly supported a hypothesis that the sequence was an ancestral motif of insect telomeres but was lost repeatedly during insect evolution. The motif was conserved in the "primitive" apterous insect orders, the Archaeognatha and Zygentoma, in the "lower" Neoptera (Plecoptera, Phasmida, Orthoptera, Blattaria, Mantodea, and Isoptera) with the exception of Dermaptera, and in Paraneoptera (Psocoptera, Thysanoptera, Auchenorrhyncha, and Sternorrhyncha) with the exception of Heteroptera. Surprisingly, the (TTAGG)n motif was not found in the "primitive" pterygotes, the Palaeoptera (Ephemeroptera and Odonata). The Endopterygota were heterogeneous for the occurrence of TTAGG repeats. The motif was conserved in Hymenoptera, Lepidoptera, and Trichoptera but was lost in one clade formed by Diptera, Siphonaptera, and Mecoptera. It was also lost in Raphidioptera, whereas it was present in Megaloptera. In contrast with previous authors, we did not find the motif in Neuroptera. Finally, both TTAGG-positive and TTAGG-negative species were reported in Coleoptera. The repeated losses of TTAGG in different branches of the insect phylogenetic tree and, in particular, in the most successful lineage of insect evolution, the Endopterygota, suggest a backup mechanism in the genome of insects that enabled them frequent evolutionary changes in telomere composition.     

Comparing the qualitatively different effects rapidly evolving and rapidly induced defences have on predator-prey interactions

Interspecific interactions depend not only on the population densities of the interacting species, but on their phenotypes as well. Phenotypic variation can be plastic or heritable and both mechanisms can drive phenotypic change at rates comparable to or faster than those of ecological dynamics (e.g. changes in population abundances or spatial distributions). In this study, we compare the effects rapidly induced and rapidly evolving defences have on community dynamics by considering the fast phenotypic change limit using fast-slow systems theory. Our approach allows us to study phenotypically plastic and evolving systems with one overarching theory, thus capturing the effects rapidly induced defences have on ecological dynamics and how those effects differ from the effects of evolving defences. Our results show that rapidly induced defences tend to stabilize community dynamics and that some behaviours observed in rapidly evolving systems cannot be produced by phenotypic plasticity.     

Evolution in the fast lane: rapidly evolving sex-related genes in Drosophila

A large portion of the annotated genes in Drosophila melanogaster show sex-biased expression, indicating that sex and reproduction-related genes (SRR genes) represent an appreciable component of the genome. Previous studies, in which subsets of genes were compared among few Drosophila species, have found that SRR genes exhibit unusual evolutionary patterns. Here, we have used the newly released genome sequences from 12 Drosophila species, coupled to a larger set of SRR genes, to comprehensively test the generality of these patterns. Among 2505 SRR genes examined, including ESTs with biased expression in reproductive tissues and genes characterized as involved in gametogenesis, we find that a relatively high proportion of SRR genes have experienced accelerated divergence throughout the genus Drosophila. Several testis-specific genes, male seminal fluid proteins (SFPs), and spermatogenesis genes show lineage-specific bursts of accelerated evolution and positive selection. SFP genes also show evidence of lineage-specific gene loss and/or gain. These results bring us closer to understanding the details of the evolutionary dynamics of SRR genes with respect to species divergence.     

Phylogenetic indices for measuring the diet breadths of phytophagous insects

Prevailing methods of measuring diet breadth of phytophagous insects are not consistent between studies and generally rely on counts of a variety of higher plant taxa (e.g. genera, families, orders). Results derived from them can be inconsistent if different taxonomic levels are used between studies. In any case, such indices do not include information from the whole branching structure of the host plant phylogeny, and do not address the fact that higher taxa are not necessarily phylogenetically equivalent. Here we present novel phylogeny-based methods which address these shortcomings. Although a previously proposed index (the Phylogenetic Diversity index) may be employed, it cannot be used to measure diets of strictly monophagous insects (i.e. those which utilise a single host species). We therefore introduce a modification of this index (the Root Phylogenetic Diversity index) which may be applied to all diets. In addition, we propose a Clade Dispersion index as a branch-length-independent measure of the degree to which hosts are scattered across the host phylogeny. We describe how these indices could be employed in studies of insect diet breadth and discuss potential problems which may be encountered in their use. Received: 16 November 1998 / Accepted: 10 February 1999     

Phylogenetic analysis of cubilin (CUBN) gene

Cubilin, (CUBN; also known as intrinsic factor-cobalamin receptor [Homo sapiens Entrez Pubmed ref {"type":"entrez-nucleotide","attrs":{"text":"NM_001081.3","term_id":"126091151","term_text":"NM_001081.3"}}NM_001081.3; {"type":"entrez-nucleotide","attrs":{"text":"NG_008967.1","term_id":"212549718","term_text":"NG_008967.1"}}NG_008967.1; GI: 119606627]), located in the epithelium of intestine and kidney acts as a receptor for intrinsic factor – vitamin B12 complexes. Mutations in CUBN may play a role in autosomal recessive megaloblastic anemia. The current study investigated the possible role of CUBN in evolution using phylogenetic testing. A total of 588 BLAST hits were found for the cubilin query sequence and these hits showed putative conserved domain, CUB superfamily (as on 27^th Nov 2012). A first-pass phylogenetic tree was constructed to identify the taxa which most often contained the CUBN sequences. Following this, we narrowed down the search by manually deleting sequences which were not CUBN. A repeat phylogenetic analysis of 25 taxa was performed using PhyML, RAxML and TreeDyn softwares to confirm that CUBN is a conserved protein emphasizing its importance as an extracellular domain and being present in proteins mostly known to be involved in development in many chordate taxa but not found in prokaryotes, plants and yeast.. No horizontal gene transfers have been found between different taxa.     

Phylogenetic analysis of chemokine-binding protein gene from orthopoxviruses

Poxviruses, belonging to the genus of Orthopoxviruses and Leporipoxviruses, express a protein that binds a wide spectrum of CC-chenokines. The above viral protein (vCCI) has no homology by the amino acid sequence with the known cellular receptors of chemokines or with any other proteins represented in databases. Sixty-nine nucleotide sequences of the vCCI protein gene of 4 Orthopoxviruses, which are pathogenic to man, and of 2 species, to whom man is immune, as well as 16 of such sequences borrowed from GenBank were subjected to the phylogenetic analysis. Strains of Orthopoxviruses, except for cowpox viruses (CPXV), were found to cluster together in groups according to their specific attribute. CPXV isolates can be conditionally divided into 2 subgroups, which are much remote from one another evolutionally. Four subgroups were defined for smallpox virus with their strains originating from different geographic regions.     

Delineating slowly and rapidly evolving fractions of the Drosophila genome.

Evolutionary conservation is an important indicator of function and a major component of bioinformatic methods to identify non-protein-coding genes. We present a new Bayesian method for segmenting pairwise alignments of eukaryotic genomes while simultaneously classifying segments into slowly and rapidly evolving fractions. We also describe an information criterion similar to the Akaike Information Criterion (AIC) for determining the number of classes. Working with pairwise alignments enables detection of differences in conservation patterns among closely related species. We analyzed three whole-genome and three partial-genome pairwise alignments among eight Drosophila species. Three distinct classes of conservation level were detected. Sequences comprising the most slowly evolving component were consistent across a range of species pairs, and constituted approximately 62-66% of the D. melanogaster genome. Almost all (>90%) of the aligned protein-coding sequence is in this fraction, suggesting much of it (comprising the majority of the Drosophila genome, including approximately 56% of non-protein-coding sequences) is functional. The size and content of the most rapidly evolving component was species dependent, and varied from 1.6% to 4.8%. This fraction is also enriched for protein-coding sequence (while containing significant amounts of non-protein-coding sequence), suggesting it is under positive selection. We also classified segments according to conservation and GC content simultaneously. This analysis identified numerous sub-classes of those identified on the basis of conservation alone, but was nevertheless consistent with that classification. Software, data, and results available at www.maths.qut.edu.au/-keithj/. Genomic segments comprising the conservation classes available in BED format.