Single-copy nuclear genes resolve the phylogeny of the holometabolous insects

Background  

Evolutionary relationships among the 11 extant orders of insects that undergo complete metamorphosis, called Holometabola, remain either unresolved or contentious, but are extremely important as a context for accurate comparative biology of insect model organisms. The most phylogenetically enigmatic holometabolan insects are Strepsiptera or twisted wing parasites, whose evolutionary relationship to any other insect order is unconfirmed. They have been controversially proposed as the closest relatives of the flies, based on rDNA, and a possible homeotic transformation in the common ancestor of both groups that would make the reduced forewings of Strepsiptera homologous to the reduced hindwings of Diptera. Here we present evidence from nucleotide sequences of six single-copy nuclear protein coding genes used to reconstruct phylogenetic relationships and estimate evolutionary divergence times for all holometabolan orders.     

Single-copy nuclear DNA sequences obtained from noninvasively collected primate feces

Noninvasively collected primate feces have been shown to provide a useful source of mitochondrial DNA for sequencing and nuclear microsatellite DNA for size analysis. In this study, single-copy nuclear DNA sequences were obtained from noninvasively collected fecal samples of two species of wild tamarins, Saguinus fuscicollis and S. mystax, in the context of a project on the functional utility of color vision. Noninvasive genotyping of the X-linked opsin gene is important for future studies of selection and adaptation at this locus in a number of primate species. The wide range of techniques that can now be applied successfully to DNA extracted from feces introduces a broad spectrum of potential genetic studies that can be undertaken on primates, without the need for intrusive or invasive methods.     

Phylogeny of colletid bees (Hymenoptera: Colletidae) inferred from four nuclear genes

Colletidae comprise approximately 2500 species of bees primarily distributed in the southern continents (only two colletid genera are widely distributed: Colletes and Hylaeus). Previously published studies have failed to resolve phylogenetic relationships on a worldwide basis and this has been a major barrier to the progress of research regarding systematics and evolution of colletid bees. For this study, data from four nuclear gene loci: elongation factor-1alpha (F2 copy), opsin, wingless, and 28S rRNA were analyzed for 122 species of colletid bees, representing all subfamilies and tribes currently recognized; 22 species belonging to three other bee families were used as outgroups. Bayesian, maximum likelihood, and parsimony methods were employed to investigate the phylogenetic relationships within Colletidae and resulted in highly congruent and well-resolved trees. The phylogenetic results show that Colletidae are monophyletic and that all traditionally recognized subfamilies (except Paracolletinae) are also strongly supported as monophyletic. Our phylogenetic hypothesis provides a framework within which broad questions related to the taxonomy, biogeography, morphology, evolution, and ecology of colletid bees can be addressed.     

Nuclear genes resolve mesozoic-aged divergences in the insect order Lepidoptera

Compared to the number of genes available for study of both younger and older divergences, few genes have yet been identified that can strongly resolve phylogenetic splits of Mesozoic age ( approximately 65-250 mya). Thus, reconstruction of Mesozoic-age phylogenies, exemplified by basal divergences within the major orders of holometabolous insects, is likely to be especially dependent on combining multiple lines of evidence. This study tests the potential of the 18S ribosomal RNA gene for reconstructing Mesozoic-aged divergences within the insect order Lepidoptera and its ability when combined with a second, previously analyzed nuclear gene (phosphoenolpyruvate carboxykinase, PEPCK) to strongly resolve these relationships. 18S sequences were obtained for 21 taxa, representing major clades of Lepidoptera plus outgroups from the other "panorpoid orders. A well-corroborated morphology-based "test phylogeny was used to evaluate the effects of partitioning the 18S gene according to variable versus conserved domains, paired versus unpaired sites in the secondary structure, and transition versus transversion substitutions. Likelihood and unweighted parsimony analyses of the 18S data recover the "test phylogeny" almost completely, with no improvement of agreement or support provided by any form of weighting or partitioning. No conflict in signal between 18S and PEPCK was detected by the partition homogeneity test. Combined parsimony analysis yielded strong bootstrap support for nearly all relationships, much higher than for either gene alone, thereby also providing strong evidence on several hypotheses about the early evolution of lepidopteran-plant interactions. These genes in combination may be widely useful for resolving insect divergences of comparable age.     

Genetic divergences and phylogenetic relationships among the Fejervarya limnocharis complex in Thailand and neighboring countries revealed by mitochondrial and nuclear genes

To clarify the genetic divergence in the F. limnocharis complex from Thailand and neighboring countries and to elucidate the phylogenetic problems of this taxon, we analyzed partial sequences of the mitochondrial 12S and 16S rRNA genes and the nuclear CXCR4, NCX1, RAG-1, and tyrosinase genes. The F. limnocharis complex from Thailand had three distinct haplotypes for 12S and 16S rRNA genes. Nucleotide similarities and the phylogenetic relationships indicated that the haplotype 1 group corresponded to the real "F. limnocharis", the haplotype 2 group was F. orissaensis or closely related to it, and the haplotype 3 group was possibly an undescribed species. Mitochondrial gene data also showed two major clades of the genus Fejervarya, the Southeastern and South Asian groups. Although F. orissaensis is so far known only from Orissa in India, the haplotype 2 group was observed in Thailand. This distribution pattern and the phylogeny suggested that the origin of F. orissaensis and the haplotype 2 group might lie in Southeast Asia. There was also evidence suggesting that the haplotype 3 group originated in the South Asian area and has spread to northern Thailand. The nuclear gene data did not support the monophyly of the haplotypes recognized by mitochondrial genes. This incongruence between the mitochondrial and nuclear data seems to be caused by ancestral polymorphic sites contained in nuclear genes. Although neither the mitochondrial nor the nuclear data clarified intergeneric relationships, the nuclear data rejected the monophyly of the genus Fejervarya.     

Evolution and phylogenetic utility of alignment gaps within intron sequences of three nuclear genes in bumble bees (Bombus)

To test whether gaps resulting from sequence alignment contain phylogenetic signal concordant with those of base substitutions, we analyzed the occurrence of indel mutations upon a well-resolved, substitution-based tree for three nuclear genes in bumble bees (Bombus, Apidae: Bombini). The regions analyzed were exon and intron sequences of long-wavelength rhodopsin (LW Rh), arginine kinase (ArgK), and elongation factor-1alpha (EF-1alpha) F2 copy genes. LW Rh intron had only a few uninformative gaps, ArgK intron had relatively long gaps that were easily aligned, and EF-1alpha intron had many short gaps, resulting in multiple optimal alignments. The unambiguously aligned gaps within ArgK intron sequences showed no homoplasy upon the substitution-based tree, and phylogenetic signals within ambiguously aligned regions of EF-1alpha intron were highly congruent with those of base substitutions. We further analyzed the contribution of gap characters to phylogenetic reconstruction by incorporating them in parsimony analysis. Inclusion of gap characters consistently improved support for nodes recovered by substitutions, and inclusion of ambiguously aligned regions of EF-1alpha intron resolved several additional nodes, most of which were apical on the phylogeny. We conclude that gaps are an exceptionally reliable source of phylogenetic information that can be used to corroborate and refine phylogenies hypothesized by base substitutions, at least at lower taxonomic levels. At present, full use of gaps in phylogenetic reconstruction is best achieved in parsimony analysis, pending development of well-justified and generally applicable methods for incorporating indels in explicitly model-based methods.     

Hemolymph proteome changes during worker brood development match the biological divergences between western honey bees (Apis mellifera) and eastern honey bees (Apis cerana)

Background

Hemolymph plays key roles in honey bee molecule transport, immune defense, and in monitoring the physiological condition. There is a lack of knowledge regarding how the proteome achieves these biological missions for both the western and eastern honey bees (Apis mellifera and Apis cerana). A time-resolved proteome was compared using two-dimensional electrophoresis-based proteomics to reveal the mechanistic differences by analysis of hemolymph proteome changes between the worker bees of two bee species during the larval to pupal stages.

Results

The brood body weight of Apis mellifera was significantly heavier than that of Apis cerana at each developmental stage. Significantly, different protein expression patterns and metabolic pathways were observed in 74 proteins (166 spots) that were differentially abundant between the two bee species. The function of hemolymph in energy storage, odor communication, and antioxidation is of equal importance for the western and eastern bees, indicated by the enhanced expression of different protein species. However, stronger expression of protein folding, cytoskeletal and developmental proteins, and more highly activated energy producing pathways in western bees suggests that the different bee species have developed unique strategies to match their specific physiology using hemolymph to deliver nutrients and in immune defense.

Conclusions

Our disparate findings constitute a proof-of-concept of molecular details that the ecologically shaped different physiological conditions of different bee species match with the hemolymph proteome during the brood stage. This also provides a starting point for future research on the specific hemolymph proteins or pathways related to the differential phenotypes or physiology.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-563) contains supplementary material, which is available to authorized users.     

Basal divergences in birds and the phylogenetic utility of the nuclear RAG-1 gene.

The single-copy RAG-1 gene is found throughout higher vertebrates and consists of a single 3.1-kb exon without intervening introns. A 2.9-kb region of the RAG-1 locus was sequenced for 14 basal taxa of birds plus the crocodylian outgroups Alligator and Gavialis. Phylogenetic analysis of the sequences supported the hypothesis that the deepest evolutionary split in extant birds separates paleognaths from neognaths. A deep division among neognaths separates the chicken- and duck-like birds ("galloanserines") from a clade consisting of all other birds ("plethornithines"). The relationships of these three basal clades in Aves were supported by high bootstrap (98 to 100%) and large decay index values (above 14). Additionally, the plethornithine clade is characterized by a 15-bp (five-codon) synapomorphic deletion relative to all other birds. RAG-1 evolves slowly, with a number of properties favoring its phylogenetic utility, including rarity of indels, minimal saturation of transition changes at 3rd positions of codons, nearly constant base composition across taxa, and no asymmetry in directional patterns of reconstructed change. However, RAG-1 does not evolve in a clocklike manner, suggesting that this gene cannot easily be used for estimating ages of ancient lineages.     

Detecting correlation between sequence and expression divergences in a comparative analysis of human serpin genes

Physiological functions and characteristic structures of the serpin gene superfamily have been studied extensively, yet the evolution of the serpin genes remains unclear. Gene duplication in this superfamily may shed light on this issue. Two models are used to predict the preservation of duplicated genes: the classical model and the duplication-degeneration-complementation (DDC) model. In this study, we analyzed the phylogenetic relationships of 33 human serpin genes and the expression data of some members of the serpin superfamily from a DNA microarray of human leukemia U937 cells with stably inducible expression of the leukemia-related AML1-ETO gene. We then determined the utility of the DDC model by mapping serpin superfamily expression data to the phylogenetic tree. The correlation between sequence and expression divergences as measured by the Pearson correlation coefficient indicated that human serpin genes evolved under the DDC model. Our study provides a new strategy for comparative analysis of gene sequences and microarray data.     

Single-copy genes define a conserved order between rice and wheat for understanding differences caused by duplication, deletion, and transposition of genes

The high-quality rice genome sequence is serving as a reference for comparative genome analysis in crop plants, especially cereals. However, early comparisons with bread wheat showed complex patterns of conserved synteny (gene content) and colinearity (gene order). Here, we show the presence of ancient duplicated segments in the progenitor of wheat, which were first identified in the rice genome. We also show that single-copy (SC) rice genes, those representing unique matches with wheat expressed sequence tag (EST) unigene contigs in the whole rice genome, show more than twice the proportion of genes mapping to syntenic wheat chromosome as compared to the multicopy (MC) or duplicated rice genes. While 58.7% of the 1,244 mapped SC rice genes were located in single syntenic wheat chromosome groups, the remaining 41.3% were distributed randomly to the other six non-syntenic wheat groups. This could only be explained by a background dispersal of genes in the genome through transposition or other unknown mechanism. The breakdown of rice–wheat synteny due to such transpositions was much greater near the wheat centromeres. Furthermore, the SC rice genes revealed a conserved primordial gene order that gives clues to the origin of rice and wheat chromosomes from a common ancestor through polyploidy, aneuploidy, centromeric fusions, and translocations. Apart from the bin-mapped wheat EST contigs, we also compared 56,298 predicted rice genes with 39,813 wheat EST contigs assembled from 409,765 EST sequences and identified 7,241 SC rice gene homologs of wheat. Based on the conserved colinearity of 1,063 mapped SC rice genes across the bins of individual wheat chromosomes, we predicted the wheat bin location of 6,178 unmapped SC rice gene homologs and validated the location of 213 of these in the telomeric bins of 21 wheat chromosomes with 35.4% initial success. This opens up the possibility of directed mapping of a large number of conserved SC rice gene homologs in wheat. Overall, only 46.4% of these SC genes code for proteins with known functional domains; the remaining 53.6% have unknown function, and hence, represent an important, but yet, under explored category of genes. Electronic supplementary material Supplementary material is available for this article at and is accessible for authorized users.     

Detecting correlation between sequence and expression divergences in a comparative analysis of human serpin genes

Physiological functions and characteristic structures of the serpin gene superfamily have been studied extensively, yet the evolution of the serpin genes remains unclear. Gene duplication in this superfamily may shed light on this issue. Two models are used to predict the preservation of duplicated genes: the classical model and the duplication–degeneration–complementation (DDC) model. In this study, we analyzed the phylogenetic relationships of 33 human serpin genes and the expression data of some members of the serpin superfamily from a DNA microarray of human leukemia U937 cells with stably inducible expression of the leukemia-related AML1-ETO gene. We then determined the utility of the DDC model by mapping serpin superfamily expression data to the phylogenetic tree. The correlation between sequence and expression divergences as measured by the Pearson correlation coefficient indicated that human serpin genes evolved under the DDC model. Our study provides a new strategy for comparative analysis of gene sequences and microarray data.     

Single-copy and amplified CAD genes in Syrian hamster chromosomes localized by a highly sensitive method for in situ hybridization.

Syrian hamster cells resistant to N-(phosphonacetyl)-L-aspartate (PALA), a specific inhibitor of the aspartate transcarbamylase activity of the multifunctional protein CAD, overproduce this protein as a result of amplification of the CAD gene. We have used a sensitive in situ hybridization technique to localize CAD genomes in spreads of metaphase chromosomes from several independent PALA-resistant lines and from wild-type PALA-sensitive cells. The amplified genes were always found within chromosomes, usually in an expanded region of the short arm of chromosome B9. In wild-type cells, the CAD gene was also on the short arm of chromosome B9. In one mutant line, 90 to 100 CAD genes were found within an expanded B9 chromosome and 10 to 15 more were near the distal end of one arm of several different chromosomes. Another line contained most the genes in a telomeric chromosome or large chromosome fragment. The amplified genes were in chromosomal regions that were stained in a banded pattern by trypsin-Giemsa. A few double minute chromosomes were observed in a very small fraction of the total spreads examined. The it situ hybridizations were performed in the presence of 10% dextral sulfate 500, which increases the signal by as much as 100-fold. Using recombinant DNA plasmids nick-translated with [125I]dCTP to high specific radioactivity, 10 CAD genes in a single chromosomal region were revealed after 1 week of autoradiographic exposure, and the position of the unique gene could be seen after 1 month.     

Variation in synonymous substitution rates among mammalian genes and the correlation between synonymous and nonsynonymous divergences

Using mammalian gene sequences, the variances in the numbers of synonymous and nonsynonymous substitutions among genes were estimated together with the correlation coefficient between the two. The expected correlation coefficient can be obtained under the neutral theory using these estimated values of the variances. The expected coefficient is found to often be one-half to two-thirds of the observed value. Possible causes for the disagreement were discussed, such as correlated selective constraints on the two types of substitutions and excess doublet mutations. The variance of mutation rate and that of selective constraint were also estimated. The results show that the coefficient of variation of the former is 0.2–0.3, whereas that of the latter is 0.7–0.9. Correspondence to: T. Ohta     

Conservation of the drought-inducible DS2 genes and divergences from their ASR paralogues in solanaceous species.

The drought-inducible DS2 genes of potatoes are members of the ASR (abscisic acid, stress and ripening) gene family. Previously it was shown that expression of DS2 genes is highly dehydration-specific in potato leaves, is not inducible by cold, heat, salt, hypoxia or oxidative stresses, and is independent of abscisic acid (ABA). Now it is shown that StDS2 does not respond either to sucrose or any plant hormones. Conservation of DS2 genes with this unique mode of regulation was studied in the solanaceous species with different relationships to potatoes. DS2 orthologues were identified by DNA sequence alignment in the closely related Lycopersicon and Capsicum species but not in the more distantly related Nicotiana sp. DNA and RNA gel blot analysis revealed the presence of a gene highly homologous to the potato gene StDS2 in tomato (LeDS2) with the same desiccation-specific expression in leaves and organ-specific expression in flowers and green fruits. The LeDS2 promoter was isolated and found to be almost identical in sequence with the promoter of StDS2, except for a 45-bp insertion in tomato. In contrast, no gene highly similar to StDS2 was detected in Nicotiana species on DNA gel blots. Neither StDS2 nor LeDS2 promoter regions were able to confer expression for the beta-glucuronidase (GUS) reporter gene in transgenic tobacco plants indicating that the trans regulatory factors necessary for DS2 expression are not conserved either in Nicotiana tabacum. These data suggest a narrow species-specificity and late evolution of the DS2-type genes within the family Solanaceae.     

Analysis of family-level relationships in bees (Hymenoptera: Apiformes) using 28S and two previously unexplored nuclear genes: CAD and RNA polymerase II

We analyzed a combined data set of two protein-coding nuclear genes (CAD and RNA polymerase II) and a nuclear ribosomal gene (28S D2-D4 region) for 68 bee species and 11 wasp outgroups. Our taxon sampling included all seven extant bee families, 17 of 20 subfamilies, and diverse tribes. Wasp outgroups included the two families most closely related to bees: Crabronidae and Sphecidae. We analyzed the combined and single gene data sets using parsimony and Bayesian methods, which yielded largely congruent results. Our results provide reasonably strong support for family and subfamily-level relationships among bees. Our data set strongly supports the sister-group relationship of the Colletidae and Stenotritidae, and places Halictidae as sister to this clade combined. Our analyses place the Melittidae and the long-tongued (LT) bee clade (Apidae+Megachilidae) near the base of the tree with Colletidae (and Stenotritidae) in a fairly highly derived position. This topology ("Melittidae-LT basal") was obtained in previous morphological studies under certain methods of character coding. A more widely accepted tree topology that places Colletidae (and/or Stenotritidae) as sister to all other bees ("Colletidae basal") is not supported by our data. The "Melittidae-LT basal" hypothesis may better explain patterns in the bee fossil record as well as historical biogeography of certain bee groups. Our results provide new insights into higher-level bee phylogeny and indicate that CAD, RNA polymerase II, and 28S are useful data sets for resolving Cretaceous-age divergences in bees and other Hymenoptera.     

A molecular phylogeny of bivalve mollusks: ancient radiations and divergences as revealed by mitochondrial genes

Background

Bivalves are very ancient and successful conchiferan mollusks (both in terms of species number and geographical distribution). Despite their importance in marine biota, their deep phylogenetic relationships were scarcely investigated from a molecular perspective, whereas much valuable work has been done on taxonomy, as well as phylogeny, of lower taxa.

Methodology/Principal Findings

Here we present a class-level bivalve phylogeny with a broad sample of 122 ingroup taxa, using four mitochondrial markers (MT-RNR1, MT-RNR2, MT-CO1, MT-CYB). Rigorous techniques have been exploited to set up the dataset, analyze phylogenetic signal, and infer a single final tree. In this study, we show the basal position of Opponobranchia to all Autobranchia, as well as of Palaeoheterodonta to the remaining Autobranchia, which we here propose to call Amarsipobranchia. Anomalodesmata were retrieved as monophyletic and basal to (Heterodonta + Pteriomorphia).

Conclusions/Significance

Bivalve morphological characters were traced onto the phylogenetic trees obtained from the molecular analysis; our analysis suggests that eulamellibranch gills and heterodont hinge are ancestral characters for all Autobranchia. This conclusion would entail a re-evaluation of bivalve symplesiomorphies.     

Rates of synonymous substitution in plant nuclear genes

Summary The rate of synonymous nucleotide substitution in nuclear genes of higher plants has been estimated. The rate varies among genes by a factor of up to two, in a manner that is not immediately explicable in terms of base composition or codon usage bias. The average rate, in both monocots and dicots, is about four times higher than that in chloroplast genes. This leads to an estimated absolute silent substitution rate of 6 × 10^–9 substitutions per site per year that falls within the range of average rates (2–8 × 10^–9) seen in different mammalian nuclear genomes.