New heuristic methods for joint species delimitation and species tree inference.
O'Meara Brian C
Systematic biology
Species delimitation and species tree inference are difficult problems in cases of recent divergence, especially when different loci have different histories. This paper quantifies the difficulty of jointly finding the division of samples to species and estimating a species tree without constraining the possible assignments a priori. It introduces a parametric and a nonparametric method, including new heuristic search strategies, to do this delimitation and tree inference using individual gene trees as input. The new methods were evaluated using thousands of simulations and 4 empirical data sets. These analyses suggest that the new methods, especially the nonparametric one, may provide useful insights for systematists working at the species level with molecular data. However, they still often return incorrect results.
10.1093/sysbio/syp077
Locust retinoid X receptors: 9-Cis-retinoic acid in embryos from a primitive insect.
Nowickyj Shaun M,Chithalen James V,Cameron Don,Tyshenko Michael G,Petkovich Martin,Wyatt Gerard R,Jones Glenville,Walker Virginia K
Proceedings of the National Academy of Sciences of the United States of America
The retinoid X receptor (RXR) is activated by its often elusive cognate ligand, 9-cis-retinoic acid (9-cis-RA). In flies and moths, molting is mediated by a heterodimer ecdysone receptor consisting of the ecdysone monomer (EcR) and an RXR homolog, ultraspiracle (USP); the latter is believed to have diverged from its RXR origin. In the more primitive insect, Locusta migratoria (Lm), RXR is more similar to human RXRs than to USPs. LmRXR was detected in early embryos when EcR transcripts were absent, suggesting another role apart from ecdysone signaling. Recombinant LmRXRs bound 9-cis-RA and all-trans-RA with high affinity (IC(50) = 61.2-107.7 nM; K(d) = 3 nM), similar to human RXR. To determine whether specific binding had functional significance, the presence of endogenous retinoids was assessed. Embryos were extracted by using modified Bligh and Dyer and solid-phase protocols to avoid the oily precipitate that makes this material unsuitable for assay. These extracts contained retinoids (5.4 nM) as assessed by RA-inducible Cyp26A1-promoter luciferase reporter cell lines. Furthermore, the use of HPLC and MS confirmed the presence of retinoids and identified in any embryo, 9-cis-RA, in addition to all-trans-RA. We estimate that whole embryos contain 3 nM RA, including 9-cis-RA at a concentration of 1.6 nM. These findings strongly argue for a functional role for retinoids in primitive insects and favor a model where signaling through the binding of 9-cis-RA to its RXR is established relatively early in evolution and embryonic development.
10.1073/pnas.0712132105
Coupling genetic and ecological-niche models to examine how past population distributions contribute to divergence.
Knowles L Lacey,Carstens Bryan C,Keat Marcia L
Current biology : CB
Understanding the impact of climate-induced distributional shifts on species divergence, like those accompanying the Pleistocene glacial cycles [1, 2], requires tools that explicitly incorporate the geographic configuration of past distributions into analyses of genetic differentiation. Depending on the historical distribution of species, genetic differences may accumulate among ancestral source populations, but there is long-standing debate whether displacements into glacial refugia promoted divergence. Here we integrate coalescent-based genetic models [3, 4] with ecological-niche modeling [5, 6] to generate expectations for patterns of genetic variation based on an inferred past distribution of a species. Reconstruction of the distribution of a montane grasshopper species during the last glacial maximum suggests that Melanoplus marshalli populations from the sky islands of Colorado and Utah were likely colonized from multiple ancestral source populations. The genetic analyses provide compelling evidence that the historical distribution of M. marshalli-namely, spatial separation of multiple refugia-was conducive to genetic differentiation. The coupling of genetic and ecological-niche modeling provides a new and flexible tool for integrating paleoenvironmental details into species-specific predictions of population structure that can increase our understanding of why the glacial cycles promoted speciation in some taxa and yet inhibited diversification in others [7, 8].
10.1016/j.cub.2007.04.033
Horizontal gene transfer allowed the emergence of broad host range entomopathogens.
Zhang Qiangqiang,Chen Xiaoxuan,Xu Chuan,Zhao Hong,Zhang Xing,Zeng Guohong,Qian Ying,Liu Ran,Guo Na,Mi Wubin,Meng Yamin,Leger Raymond J St,Fang Weiguo
Proceedings of the National Academy of Sciences of the United States of America
The emergence of new pathogenic fungi has profoundly impacted global biota, but the underlying mechanisms behind host shifts remain largely unknown. The endophytic insect pathogen evolved from fungi that were plant associates, and entomopathogenicity is a more recently acquired adaptation. Here we report that the broad host-range entomopathogen has 18 genes that are derived via horizontal gene transfer (HGT). The necessity of degrading insect cuticle served as a major selective pressure to retain these genes, as 12 are up-regulated during penetration; 6 were confirmed to have a role in penetration, and their collective actions are indispensable for infection. Two lipid-carrier genes are involved in utilizing epicuticular lipids, and a third (MrNPC2a) facilitates hemocoel colonization. Three proteases degraded the procuticular protein matrix, which facilitated up-regulation of other cuticle-degrading enzymes. The three lipid carriers and one of the proteases are present in all analyzed species and are essential for entomopathogenicity. Acquisition of another protease (MAA_01413) in an ancestor of broad host-range lineages contributed to their host-range expansion, as heterologous expression in the locust specialist enabled it to kill caterpillars. Our work reveals that HGT was a key mechanism in the emergence of entomopathogenicity in from a plant-associated ancestor and in subsequent host-range expansion by some lineages.
10.1073/pnas.1816430116
Ecological and evolutionary drivers of the elevational gradient of diversity.
Laiolo Paola,Pato Joaquina,Obeso José Ramón
Ecology letters
Ecological, evolutionary, spatial and neutral theories make distinct predictions and provide distinct explanations for the mechanisms that control the relationship between diversity and the environment. Here, we test predictions of the elevational diversity gradient focusing on Iberian bumblebees, grasshoppers and birds. Processes mediated by local abundance and regional diversity concur in explaining local diversity patterns along elevation. Effects expressed through variation in abundance were similar among taxa and point to the overriding role of a physical factor, temperature. This determines how energy is distributed among individuals and ultimately how the resulting pattern of abundance affects species incidence. Effects expressed through variation in regional species pools depended instead on taxon-specific evolutionary history, and lead to diverging responses under similar environmental pressures. Local filters and regional variation also explain functional diversity gradients, in line with results from species richness that indicate an (local) ecological and (regional) historical unfolding of diversity-elevation relationships.
10.1111/ele.12967
Estimating species phylogeny from gene-tree probabilities despite incomplete lineage sorting: an example from Melanoplus grasshoppers.
Carstens Bryan C,Knowles L Lacey
Systematic biology
Estimating phylogenetic relationships among closely related species can be extremely difficult when there is incongruence among gene trees and between the gene trees and the species tree. Here we show that incorporating a model of the stochastic loss of gene lineages by genetic drift into the phylogenetic estimation procedure can provide a robust estimate of species relationships, despite widespread incomplete sorting of ancestral polymorphism. This approach is applied to a group of montane Melanoplus grasshoppers for which genealogical discordance among loci and incomplete lineage sorting obscures any obvious phylogenetic relationships among species. Unlike traditional treatments where gene trees estimated using standard phylogenetic methods are implicitly equated with the species tree, with the coalescent-based approach the species tree is modeled probabilistically from the estimated gene trees. The estimated species phylogeny (the ESP) is calculated for the grasshoppers from multiple gene trees reconstructed for nuclear loci and a mitochondrial gene. This empirical application is coupled with a simulation study to explore the performance of the coalescent-based approach. Specifically, we test the accuracy of the ESP given the data based on analyses of simulated data matching the multilocus data collected in Melanoplus (i.e., data were simulated for each locus with the same number of base pairs and locus-specific mutational models). The results of the study show that ESPs can be computed using the coalescent-based approach long before reciprocal monophyly has been achieved, and that these statistical estimates are accurate. This contrasts with analyses of the empirical data collected in Melanoplus and simulated data based on concatenation of multiple loci, for which the incomplete lineage sorting of recently diverged species posed significant problems. The strengths and potential challenges associated with incorporating an explicit model of gene-lineage coalescence into the phylogenetic procedure to obtain an ESP, as illustrated by application to Melanoplus, versus concatenation and consensus approaches are discussed. This study represents a fundamental shift in how species relationships are estimated - the relationship between the gene trees and the species phylogeny is modeled probabilistically rather than equating gene trees with a species tree.
10.1080/10635150701405560
Sequence evolution in mitochondrial ribosomal and ND-1 genes in lepidoptera: implications for phylogenetic analyses.
Pashley D P,Ke L D
Molecular biology and evolution
A 2,256-bp sequence of the mitochondrial genome of a lepidopteran (Spodoptera frugiperda) contains tRNAs for valine and leucine, the 16S rRNA, and three-quarters of the ND-1 presumptive protein-coding gene. A 64-bp stretch of unknown function was located between the rRNA and leucine tRNA. Sequence divergence in the 16S rRNA obtained from alignment with published insect sequences is consistent with phylogenetic hypotheses, in that Diptera and Lepidoptera are more closely related to each other (24% sequence divergence) than either is to Hymenoptera (31%). Within the ND-1 gene, sequences for four additional Lepidoptera were generated for a 314-bp region and contrasted with published sequences for the locust and Drosophila. Sequence divergence in this region was consistent with accepted phylogenetic relationships, but results of parsimony analyses were not. Cladograms consistently recovered accepted higher level relationships (monophyly of Lepidoptera), despite high homoplasy, but were unable to resolve superfamily and family relationships within Lepidoptera, regardless of the outgroup or character subset analyzed. Character analysis indicated that homoplasy was decreased at higher levels when first- and second-codon sites were used exclusively. At the lowest level (families), resolution was enhanced by inclusion of third-codon sites. Inability of molecular data to recover a well-established phylogeny may be rectified by additional characters or taxa, but it is clear that homoplasy is sufficiently high to caution against the acceptance of relationships generated with this molecular region that are not extremely robust.
10.1093/oxfordjournals.molbev.a040778
Changing role of even-skipped during the evolution of insect pattern formation.
Patel N H,Ball E E,Goodman C S
Nature
The development of Drosophila is typical of the so-called long germband mode of insect development, in which the pattern of segments is established by the end of the blastoderm stage. Short germband insects, such as the grasshopper Schistocerca americana, by contrast, generate all or most of their metameric pattern after the blastoderm stage by the sequential addition of segments during caudal elongation. This difference is discernible at the molecular level in the pattern of initiation of the segment polarity gene engrailed, and the homeotic gene abdominal-A (ref. 5). For example, in both types of insects, engrailed is expressed by the highly conserved germband stage in a pattern of regularly spaced stripes, one stripe per segment. In Drosophila, the complete pattern is visible by the end of the blastoderm stage, although engrailed appears initially in alternate segments in a pair-rule pattern that reflects its known control by pair-rule genes such as even-skipped. In contrast, in the grasshopper, the engrailed stripes appear one at a time after the blastoderm stage as the embryo elongates. To address the molecular basis for this difference, we have cloned the grasshopper homologue of the Drosophila pair-rule gene even-skipped and show that it does not serve a pair-rule function in early development, although it does have a similar function in both insects during neurogenesis later in development.
10.1038/357339a0
Molecular characterization and phylogenetic relationships of a protein with potential oxygen-binding capabilities in the grasshopper embryo. A hemocyanin in insects?
Sánchez D,Ganfornina M D,Gutiérrez G,Bastiani M J
Molecular biology and evolution
Arthropodan hemocyanins, prophenoloxidases (PPOs), and insect hexamerins form a superfamily of hemolymph proteins that we propose to call the AHPH superfamily. The evolutionary and functional relationships of these proteins are illuminated by a new embryonic hemolymph protein (EHP) that is expressed during early stages of development in the grasshopper embryo. EHP is a 78-kDa soluble protein present initially in the yolk sac content, and later in the embryonic hemolymph. Protein purification and peptide sequencing were used to identify an embryonic cDNA clone coding for EHP. In situ hybridization identifies hemocytes as EHP-expressing cells. As deduced from the cDNA clone, EHP is a secreted protein with two potential glycosylation sites. Sequence analysis defines EHP as a member of the AHPH superfamily. Phylogenetic analyses with all the currently available AHPH proteins, including EHP, were performed to ascertain the evolutionary history of this protein superfamily. We used both the entire protein sequence and each of the three domains present in the AHPH proteins. The phylogenies inferred for each of the domains suggest a mosaic evolution of these protein modules. Phylogenetic and multivariate analyses consistently group EHP with crustacean hemocyanins and, less closely, with insect hexamerins, relative to cheliceratan hemocyanins and PPOs. The grasshopper protein rigorously preserves the residues involved in oxygen binding, oligomerization, and allosteric regulation of the oxygen transport proteins. Although insects were thought not to have hemocyanins, we propose that EHP functions as an oxygen transport or storage protein during embryonic development.
10.1093/oxfordjournals.molbev.a025938
Increased chromosomal mutation rate after hybridization between two subspecies of grasshoppers.
Shaw D D,Wilkinson P,Coates D J
Science (New York, N.Y.)
Hybridization between two chromosomally distinct subspecies of the grasshopper Caledia captiva results in a high incidence of novel chromosomal rearrangements among the backcross progeny. Rearrangements are restricted to those chromosomes derived from the F1 hybrid parent. Chromosomal involvement is nonrandom with the same rearrangement occurring repeatedly in different backcrosses. A single individual can also generate an array of different rearrangements among its offspring. Several of the rearrangements have also been found in natural populations. The nonrandom and recurrent nature of these chromosomal mutations at high frequencies provides a plausible explanation for the establishment and fixation of chromosomal rearrangements in natural populations.
10.1126/science.6407107
The heterochromatin of grasshoppers from the Caledia captiva species complex. I. Sequence evolution and conservation in a highly repeated DNA family.
Arnold M L,Appels R,Shaw D D
Molecular biology and evolution
The restriction enzyme TaqI digests 0.2% of the genomic DNA from the grasshopper Caledia captiva to a family of sequences 168 bp in length (length of consensus sequence). The sequence variation of this "Taq family" of repeat units was examined among four races from C. captiva to assay the pattern of evolution within this highly repeated DNA. The Taq-family repeats are located in C-banded heterochromatin on at least one member of each homologous pair of chromosomes; the locations range from centromeric to telomeric. Thirty-nine cloned repeats isolated from two population 1A individuals along with 11 clones from seven populations taken from three of the races demonstrated sequence variation at 72 positions. Pairwise comparisons of the cloned repeats, both within an individual and between different races, indicate that levels of intraspecific divergence, as measured by reproductive incompatibility, do not correlate with sequence divergence among the 168-bp repeats. A number of subsequences within the repeat remain unchanged among all 50 clones; the longest of these is 18 bp. That the same 18-bp subsequence is present in all clones examined is a finding that departs significantly (P less than 0.01) from what would be expected to occur at random. Two other cloned repeats, from a reproductively isolated race of C. captiva, have sequences that show 56% identity with this 18-bp conserved region. An analysis showed that the frequency of occurrence of an RsaI recognition site within the 168-bp repeat in the entire Taq family agreed with that found in the cloned sequences. These data, along with a partial sequence for the entire Taq family obtained by sequencing uncloned repeats, suggest that the consensus sequence from the cloned copies is representative of this highly repeated family and is not a biased sample resulting from the cloning procedure. The 18-bp conserved sequence is part of a 42-bp sequence that possesses dyad symmetry typical of protein-binding sites. We speculate that this may be significant in the evolution of the Taq family of sequences.
10.1093/oxfordjournals.molbev.a040380
Genomic gigantism: DNA loss is slow in mountain grasshoppers.
Bensasson D,Petrov D A,Zhang D X,Hartl D L,Hewitt G M
Molecular biology and evolution
Several studies have shown DNA loss to be inversely correlated with genome size in animals. These studies include a comparison between Drosophila and the cricket, Laupala, but there has been no assessment of DNA loss in insects with very large genomes. Podisma pedestris, the brown mountain grasshopper, has a genome over 100 times as large as that of Drosophila and 10 times as large as that of Laupala. We used 58 paralogous nuclear pseudogenes of mitochondrial origin to study the characteristics of insertion, deletion, and point substitution in P. pedestris and Italopodisma. In animals, these pseudogenes are "dead on arrival"; they are abundant in many different eukaryotes, and their mitochondrial origin simplifies the identification of point substitutions accumulated in nuclear pseudogene lineages. There appears to be a mononucleotide repeat within the 643-bp pseudogene sequence studied that acts as a strong hot spot for insertions or deletions (indels). Because the data for other insect species did not contain such an unusual region, hot spots were excluded from species comparisons. The rate of DNA loss relative to point substitution appears to be considerably and significantly lower in the grasshoppers studied than in Drosophila or Laupala. This suggests that the inverse correlation between genome size and the rate of DNA loss can be extended to comparisons between insects with large or gigantic genomes (i.e., Laupala and Podisma). The low rate of DNA loss implies that in grasshoppers, the accumulation of point mutations is a more potent force for obscuring ancient pseudogenes than their loss through indel accumulation, whereas the reverse is true for Drosophila. The main factor contributing to the difference in the rates of DNA loss estimated for grasshoppers, crickets, and Drosophila appears to be deletion size. Large deletions are relatively rare in Podisma and Italopodisma.
10.1093/oxfordjournals.molbev.a003798
Phylogenetic dispersion and diversity in regional assemblages of seed plants in China.
Qian Hong,Deng Tao,Jin Yi,Mao Lingfeng,Zhao Dan,Ricklefs Robert E
Proceedings of the National Academy of Sciences of the United States of America
Species assemble into communities through ecological and evolutionary processes. Phylogenetic niche conservatism-the tendency of species to retain ancestral ecological distributions-is thought to influence which species from a regional species pool can persist in a particular environment. We analyzed data for seed plants in China to test hypotheses about the distribution of species within regional floras. Of 16 environmental variables, actual evapotranspiration, minimum temperature of the coldest month, and annual precipitation most strongly influenced regional species richness, phylogenetic dispersion, and phylogenetic diversity for both gymnosperms (cone-bearing plants) and angiosperms (flowering plants). For most evolutionary clades at, and above, the family level, the relationships between metrics of phylogenetic dispersion (i.e., average phylogenetic distance among species), or phylogenetic diversity, and the 3 environmental variables were consistent with the tropical niche conservatism hypothesis, which predicts closer phylogenetic relatedness and reduced phylogenetic diversity with increasing environmental stress. The slopes of the relationships between phylogenetic relatedness and the 3 environmental drivers identified in this analysis were steeper for primarily tropical clades, implying greater niche conservatism, than for primarily temperate clades. These observations suggest that the distributions of seed plants across large-scale environmental gradients in China are constrained by conserved adaptations to the physical environment, i.e., phylogenetic niche conservatism.
10.1073/pnas.1822153116