Reconstruction of ancestral sequences by the inferential method,a tool for protein engineering studies

This paper describes the inferential method, an approach for reconstructing protein and nucleotide sequences of ancestral species, starting from known, homologous, contemporary sequences. The method requires knowledge of the topology of the phylogenetic tree, whose nodes are the species to whom the reconstructed sequences belong.The method has been tested by computer simulation of speciation and nucleotide substitutions, starting from a single ancestral sequence, and by subsequent reconstruction of nodal sequences. Results have shown that reconstructions obtained by the inferential method are affected by limited error frequencies, which (1) are proportional to the squares of nucleotide substitution rates and of internodal distances, and (2) are little influenced by non-uniformity of transformation rates of nucleotides.Furthermore, good agreement of the results has been obtained by comparing protein-sequence reconstructions carried out with the inferential method with those obtained using the maximum parsimony method in two different cases: e.g., a reconstruction of simulated sequences and a reconstruction of mammalian ribonuclease sequences.Abbreviations used MP maximum parsimony method - ML maximum likelihood method - IM inferential method - MY millions of years - N-tree natural-like phylogenetic tree - E-tree equibranched phylogenetic tree - EA percentage number of erroneous amino acids in a reconstructed sequence - EC percentage number of erroneous codons in a reconstructed sequence - t _n time interval between a P- and its - F-sequence nucleotides and amino acids are indicated by their I.U.B. codes (N.C.-I.U.B., 1985) Correspondence to: A. Di Donato     

Direct maximum parsimony phylogeny reconstruction from genotype data

Background  

Maximum parsimony phylogenetic tree reconstruction from genetic variation data is a fundamental problem in computational genetics with many practical applications in population genetics, whole genome analysis, and the search for genetic predictors of disease. Efficient methods are available for reconstruction of maximum parsimony trees from haplotype data, but such data are difficult to determine directly for autosomal DNA. Data more commonly is available in the form of genotypes, which consist of conflated combinations of pairs of haplotypes from homologous chromosomes. Currently, there are no general algorithms for the direct reconstruction of maximum parsimony phylogenies from genotype data. Hence phylogenetic applications for autosomal data must therefore rely on other methods for first computationally inferring haplotypes from genotypes.     

Strong signature of selection in seeder populations but not in resprouters of the fynbos heath Erica coccinea (Ericaceae)

A higher frequency of natural selection is expected in populations of organisms with shorter generation times. In fire‐prone ecosystems, populations of seeder plants behave as functionally semelparous populations, with short generation times compared to populations of resprouter plants, which are truly iteroparous. Therefore, a stronger signature of natural selection should be detected in seeder populations, favoured by their shorter generation times and higher rates of population turnover. Here we test this idea in Erica coccinea from the Cape Floristic Region, which is dimorphic for post‐fire regeneration mode. We measured three floral traits supposedly subject to natural selection in seeder and resprouter populations. We then compared phenotypic trait variation with neutral genetic variation in each group of populations using P_ST–F_ST comparisons to detect signatures of natural selection in seeders and resprouters. We found a strong signature of selection in seeder populations, but not in resprouters. Furthermore, anthers of seeders were more exserted (and larger) than those of resprouters. These differences were maintained at sites where seeders and resprouters co‐occurred, suggesting that phenotypic plasticity or adaptation to different growth environments are unlikely explanations for trait variation. These results provide empirical support for the hypothesis that the genetic signature of natural selection is certainly more intense in seeder than in resprouter populations, favoured by their comparatively faster generation turnovers. Increased frequency of natural selection would increase differentiation among populations, thus promoting speciation in pyrophytic seeder lineages of the Cape flora.     

Development and testing of a genetic marker‐based pedigree reconstruction system ‘PR‐genie’ incorporating size‐class data

For wildlife populations, it is often difficult to determine biological parameters that indicate breeding patterns and population mixing, but knowledge of these parameters is essential for effective management. A pedigree encodes the relationship between individuals and can provide insight into the dynamics of a population over its recent history. Here, we present a method for the reconstruction of pedigrees for wild populations of animals that live long enough to breed multiple times over their lifetime and that have complex or unknown generational structures. Reconstruction was based on microsatellite genotype data along with ancillary biological information: sex and observed body size class as an indicator of relative age of individuals within the population. Using body size‐class data to infer relative age has not been considered previously in wildlife genealogy and provides a marked improvement in accuracy of pedigree reconstruction. Body size‐class data are particularly useful for wild populations because it is much easier to collect noninvasively than absolute age data. This new pedigree reconstruction system, PR‐genie, performs reconstruction using maximum likelihood with optimization driven by the cross‐entropy method. We demonstrated pedigree reconstruction performance on simulated populations (comparing reconstructed pedigrees to known true pedigrees) over a wide range of population parameters and under assortative and intergenerational mating schema. Reconstruction accuracy increased with the presence of size‐class data and as the amount and quality of genetic data increased. We provide recommendations as to the amount and quality of data necessary to provide insight into detailed familial relationships in a wildlife population using this pedigree reconstruction technique.     

Establishing the foundation for an applied molecular taxonomy of otters in Southeast Asia

Four species of otters (Mustelidae, Lutrinae) occur in Southeast Asia and are considered to be of conservation concern: Aonyx cinerea (Asian small-clawed otter), Lutra lutra (Eurasian otter), Lutra sumatrana (Hairy-nosed otter), and Lutrogale perspicillata (Smooth-coated otter). Among these, L. sumatrana is endemic to the region, yet little is known about its biology, and the precise distribution of all four species in Southeast Asia is not well known. Furthermore, the taxonomy and systematics of L. sumatrana and L. perspicillata have been the subject of controversy, which has implications for the legal protection and for conservation programs of these taxa. To resolve these controversies, we used a multigene data set comprised of segments from 13 nuclear and 5 mitochondrial loci (11,180 nucleotides) to evaluate the phylogenetic relationships of Asian Old World otters. Phylogenies were also estimated using two mitochondrial loci (1,832 nucleotides) obtained from two or more individuals of the four Southeast Asian species. The results from maximum parsimony, maximum likelihood and Bayesian inference showed that L. sumatrana and L. lutra are sister taxa, whereas L. perspicillata is sister to A. cinerea. Furthermore, the results from the two-mitochondrial gene analyses indicate that L. sumatrana is reciprocally monophyletic with respect to L. lutra, supporting the specific validity of the former taxon. Signs such as tracks and feces are often used in field surveys to provide information on the distribution and abundance of otters, but the accuracy of these methods may be compromised when several closely related species occur sympatrically. Therefore, the two-gene data set was used to develop a provisional set of diagnostic nucleotides that can be potentially used to identify the four species of Southeast Asian otters from noninvasively collected biological samples, such as feces. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Miami heat: Urban heat islands influence the thermal suitability of habitats for ectotherms

The urban heat island effect, where urban areas exhibit higher temperatures than less‐developed suburban and natural habitats, occurs in cities across the globe and is well understood from a physical perspective and at broad spatial scales. However, very little is known about how thermal variation caused by urbanization influences the ability of organisms to live in cities. Ectotherms are sensitive to environmental changes that affect thermal conditions, and therefore, increased urban temperatures may pose significant challenges to thermoregulation and alter temperature‐dependent activity. To evaluate whether these changes to the thermal environment affect the persistence and dispersal of ectothermic species in urban areas, we studied two species of Anolis lizards (Anolis cristatellus and Anolis sagrei) introduced to Miami‐Dade County, FL, USA, where they occur in both urban and natural habitats. We calculated canopy openness and measured operative temperature (T_e), which estimates the distribution of body temperatures in a non‐thermoregulating population, in four urban and four natural sites. We also captured lizards throughout the day and recorded their internal body temperature (T_b). We found that urban areas had more open canopies and higher T_e compared to natural habitats. Laboratory trials showed that A. cristatellus preferred lower temperatures than A. sagrei. Urban sites currently occupied by each species appear to lower thermoregulatory costs for both species, but only A. sagreihad field T_b that were more often within their preferred temperature range in urban habitats compared to natural areas. Furthermore, based on available T_e within each species' preferred temperature range, urban sites with only A. sagrei appear less suitable for A. cristatellus, whereas natural sites with only A. cristatellus are less suitable for A. sagrei. These results highlight how the thermal properties of urban areas contribute to patterns of persistence and dispersal, particularly relevant for studying species invasions worldwide.     

Globin evolution was apparently very rapid in early vertebrates: A reasonable case against the rate-constancy hypothesis

Kimura mistook ambiguous maximum parsimony codons for wrong codons. The maximum parsimony method performed well as judged by the two classes of serine codons (which can not be connected by silent mutations) on comparing the parsimony codons for serines in human, rabbit, and mouse hemoglobin chains to actual codons determined by nucleotide sequencing. In genealogical reconstructions involving 247 eucaryotic globins, the maximum parsimony distances separating the contemporary sequences show that Kimura's Poisson and Dayhoff's PAM estimates of rate of globin evolution miss most of the superimposed replacements and are therefore seriously in error. Nor is Kimura's constant rate assumption and his belief in a single origin of myoglobin supported. Lamprey myoglobin appears to be most like lamprey hemoglobin, while gnathostome myoglobin seems closest to gnathostome hemoglobin. It was found that the three types of gnathostome globins (Mb, Hb, Hb) evolved between the shark-boney vertebrate and bird-mammal ancestors at a much faster rate than from the latter ancestor to the present. The data indicate that rates were exceedingly fast during the origin of these globin chains because a high proportion of substitutions were adaptive. It was concluded that wherever strong stabilizing selection acts on a protein, somewhere in the past positive Darwinian selection must have spread the amino acid substitutions now being preserved.     

Rapid chromosomal evolution in the mesic four‐striped grass rat Rhabdomys dilectus (Rodentia,Muridae) revealed by mtDNA phylogeographic analysis

The mesic four‐striped grass rat Rhabdomys dilectus De Winton, 1897 is distributed in mesic regions of southern and eastern Africa. We carried out a molecular and chromosomal study of the northernmost populations of the species to provide insight into the subspecific boundaries identified within the species and to describe its genetic structure in Eastern Africa. Maximum likelihood, maximum parsimony and neighbour‐joining methods were used to construct phylogenetic relationships among all the haplotypes belonging to the large part of the species range. Times of divergences were estimated assuming a relaxed molecular clock with two calibration points. We identified three well‐supported clades within R. dilectus. One basal clade corresponding to Rhabdomys d. chakae (2n = 48) is found in South Africa. Two additional sister clades corresponding to R. d. dilectus (2n = 48 and 2n = 46) are allopatrically distributed in southern and northern parts of the species range. Genetic divergence among the three clades is relatively high (ranges 4.2–5.7%). A very divergent new karyotype 2n = 38, FNa = 60 was found in two high‐altitude populations on Mt. Meru and Mt. Kilimanjaro. The karyotype differences consist in three Robertsonian fusions and one whole‐arm reciprocal translocation. Interestingly, the mtDNA phylogeny does not match with the diploid numbers. In fact, the 2n = 38 specimens form a monophyletic group within a clade that includes specimens with the 2n = 46 karyotype that appears as paraphyletic. We estimated the new karyotype originated in peripatric condition during the last phases of the Pleistocene. This study confirms the importance of chromosomal analysis in detecting taxonomic units and cryptic diversity in rodents.     

Allelic diversity at the Mhc-DQA locus of woodmouse populations (Apodemus sylvaticus) present in the islands and mainland of the northern Mediterranean

Aim Polymorphism at neutral markers and at MHC loci in rodent populations living on islands is generally low. The main genetic factors that may contribute to a reduced level of genetic variability are genetic drift, reduced gene flow and founder events. We investigated the pattern of polymorphism at the second exon of the Mhc‐DQA gene in island populations of Apodemus sylvaticus and in their mainland counterparts to investigate the pattern of MHC polymorphism in a phylogeographical context and to assess the impact of insularity on diversity at this locus. Location Eight north Mediterranean populations of Apodemus sylvaticus were studied, including five island populations (Majorca, Minorca, Porquerolles, Port‐Cros and Sicily) and three mainland populations. Methods cDNA sequencing and nucleotide sequences analyses. Synonymous and non‐synonymous substitutions were examined at the PBR and non‐PBR sites. The DQA allelic distribution in populations was compared with the woodmouse phylogeography. Results This study presents novel DQA alleles. High polymorphism of the DQA locus is recorded in natural populations of A. sylvaticus with 13 alleles being widely distributed irrespective of the geographical origin and palaeoclimatic history of populations. The DQA locus does not show the expected pattern for non‐synonymous substitutions at the PBR sites. However, island populations show a weak loss of polymorphism in comparison with their mainland counterparts. Main conclusions The DQA locus in the woodmouse seems to be subject to weak selection and does not allow resolution of phylogeographical relationships among European woodmouse populations. The presence of at least three alleles in island populations and the maintenance of five alleles between the two European lineages over 1.5 Myr suggest that balancing selection may act within populations, and more precisely within island populations, to maintain genetic variability. This study shows that phylogeographical studies are a prerequisite for any genetic investigation of selected genes in natural populations.     

Analysis on the reconstruction accuracy of the Fitch method for inferring ancestral states

Background  

As one of the most widely used parsimony methods for ancestral reconstruction, the Fitch method minimizes the total number of hypothetical substitutions along all branches of a tree to explain the evolution of a character. Due to the extensive usage of this method, it has become a scientific endeavor in recent years to study the reconstruction accuracies of the Fitch method. However, most studies are restricted to 2-state evolutionary models and a study for higher-state models is needed since DNA sequences take the format of 4-state series and protein sequences even have 20 states.     

Taxonomic implications for diaptomid copepods based on contrasting patterns of mitochondrial DNA sequence divergences in four morphospecies

Morphological identification methods do not provide reliable and meaningful species identifications for taxa where morphological differences among distinct species are either absent or overlooked (i.e., cryptic species). For example, due to the minute nature of the morphological characters used to delineate diaptomid copepod species and the apparent potential for copepod speciation to occur with little or no morphological change (i.e., morphological stasis), morphological identifications of diaptomid species may not adequately capture their true species diversity. Here, we present results from a geographic survey of mtDNA sequences from populations across the geographic ranges of four North American diaptomid species—Leptodiaptomus minutus, Skistodiaptomus pallidus, Skistodiaptomus reighardi, and Onychodiaptomus sanguineus. Shallow mitochondrial DNA sequence divergences (maximum of 1.1%) among haplotypes of L. minutus from across its geographic range suggest that current morphological identification techniques reliably identify this species. In contrast, we found large mitochondrial DNA sequence divergences (14–22%) among populations within the currently recognized morphospecies of S. pallidus, S. reighardi, and O. sanguineus. However, pairwise sequence divergences within four distinct S. pallidus clades and within populations of S. reighardi and O. sanguineus were similarly low (maximum of 1.5%) as found within L. minutus as a whole. Thus, the S. pallidus, S. reighardi, and O. sanguineus morphospecies may be considered best as cryptic species complexes. Our study therefore indicates that morphological identifications, while sufficient for some species, likely underestimate the true species diversity of diaptomid copepods. As such, we stress the need for extensive taxonomic revision that integrates genetic, morphological, reproductive, and ecological analyses of this diverse and important group of freshwater zooplankton. Furthermore, we believe an extensive taxonomic revision will shed important insight into major questions regarding the roles of geography, phylogeny, and habitat on the frequency of cryptic species on earth. Handling editor: S. I. Dodson     

Genetic structure of the genera Psorophora (Diptera: Culicidae) in Columbian and north American populations using isoenzymes and ITS2 sequences

Two a priori Psorophora columbiae and one a priori Ps. toltecum populations from Colombia were studied by means of eleven isoenzyme loci. The levels of genetic diversity, Hardy-Weinberg equilibrium, gene flow estimates, isolation by distance and genetic relationships among these three populations were studied and the results were as follows: (1) The gene diversity levels were high as well as very similar among the three populations indicating no differences between them. (2) Within each population and from a hierarchical standpoint several loci were not in Hardy-Weinberg equilibrium simultaneously by excess and defect of heterozygous. This could reflect that several natural selection forces are acting upon these mosquito populations. (3) The gene flow estimate for these populations showed the existence of this event between them. This agrees quite well with the fact that the three populations really belong to a single species. In addition, maximum parsimony analyses with 16 isoenzyme for several individuals from four Colombian and two United States Psorophora populations showed that effectively only one species was present in the Colombian area surveyed, and that this species is Ps. columbiae, which hardly contrast with that previously suggested by other authors. The text was submitted by the authors in English.     

Comparative impacts of aboveground and belowground enemies on an invasive thistle

Most research examining how herbivores and pathogens affect performance of invasive plants focuses on aboveground interactions. Although important, the role of belowground communities remains poorly understood, and the relative impact of aboveground and belowground interactions is still debated. As well, most studies of belowground interactions have been carried out in controlled environments, so little is known about the role of these interactions under natural conditions or how these relationships may change across a plant's range. Using the invasive plant Cirsium arvense, we performed a reciprocal transplant experiment to test the relative impacts of above‐ and belowground interactions at three sites across a 509‐km latitudinal gradient in its invaded range in Ontario, Canada. At each site, C. arvense seedlings were protected with above‐ and/or belowground exclosures in a factorial design. Plant performance (biomass, height, stem thickness, number of leaves, length of longest leaf, maximum rhizome length) was greatest when both above‐ and belowground exclosures were applied and lowest when no exclosures were applied. When only one type of exclosure was applied, biomass generally improved more with belowground exclosures than with aboveground exclosures. Despite site‐to‐site differences in foliar damage, root damage, and mesofaunal populations, belowground interactions generally had a greater negative impact on performance than aboveground herbivory alone. These results stress the importance of including both aboveground enemy interactions and plant–soil interactions in studies of plant community dynamics and invader performance.     

A phylogenetic analysis of Pseudopaludicola (Anura) providing evidence of progressive chromosome reduction

Here, we present a molecular phylogenetic analysis of the Neotropical genus Pseudopaludicola focusing on species relationships including 11 of the 17 known species of Pseudopaludicola; several samples of Pseudopaludicola are not assigned to any species; and 34 terminal species as an outgroup. The study was based on the analysis of approximately 2.3 kb of the sequence of the mitochondrial 12S rRNA, tRNA^val and 16S rRNA genes through maximum parsimony and Bayesian phylogenetic reconstruction approaches. Our results showed that Pseudopaludicola is a well‐supported monophyletic group organized into four major clades and confirmed that the assemblage of species that lack T‐shaped terminal phalanges is paraphyletic with respect to the P. pusilla Group. Chromosomal data mapped on the cladogram showed a direct correlation among the four clades and observed chromosome numbers (2n = 22, 20, 18 and 16) with a progressive reduction in the chromosome number. Overall, our findings suggest that some taxonomic changes are necessary and reinforce the need for a revision of the genus Pseudopaludicola.     

Efficient inference of paternity and sibship inference given known maternity via hierarchical clustering

Pedigree and sibship reconstruction are important methods in quantifying relationships and fitness of individuals in natural populations. Current methods employ a Markov chain‐based algorithm to explore plausible possible pedigrees iteratively. This provides accurate results, but is time‐consuming. Here, we develop a method to infer sibship and paternity relationships from half‐sibling arrays of known maternity using hierarchical clustering. Given 50 or more unlinked SNP markers and empirically derived error rates, the method performs as well as the widely used package Colony, but is faster by two orders of magnitude. Using simulations, we show that the method performs well across contrasting mating scenarios, even when samples are large. We then apply the method to open‐pollinated arrays of the snapdragon Antirrhinum majus and find evidence for a high degree of multiple mating. Although we focus on diploid SNP data, the method does not depend on marker type and as such has broad applications in nonmodel systems.     

Widespread positive selection in the photosynthetic Rubisco enzyme

Background  

Rubisco enzyme catalyzes the first step in net photosynthetic CO₂ assimilation and photorespiratory carbon oxidation and is responsible for almost all carbon fixation on Earth. The large subunit of Rubisco is encoded by the chloroplast rbcL gene, which is widely used for reconstruction of plant phylogenies due to its conservative nature. Plant systematicists have mainly used rbcL paying little attention to its function, and the question whether it evolves under Darwinian selection has received little attention. The purpose of our study was to evaluate how common is positive selection in Rubisco among the phototrophs and where in the Rubisco structure does positive selection occur.     

Surprisingly diverged populations of Saccharomyces cerevisiae in natural environments remote from human activity

The budding yeast, Saccharomyces cerevisiae, is a leading system in genetics, genomics and molecular biology and is becoming a powerful tool to illuminate ecological and evolutionary principles. However, little is known of the ecology and population structure of this species in nature. Here, we present a field survey of this yeast at an unprecedented scale and have performed population genetics analysis of Chinese wild isolates with different ecological and geographical origins. We also included a set of worldwide isolates that represent the maximum genetic variation of S. cerevisiae documented so far. We clearly show that S. cerevisiae is a ubiquitous species in nature, occurring in highly diversified substrates from human‐associated environments as well as habitats remote from human activity. Chinese isolates of S. cerevisiae exhibited strong population structure with nearly double the combined genetic variation of isolates from the rest of the world. We identified eight new distinct wild lineages (CHN I–VIII) from a set of 99 characterized Chinese isolates. Isolates from primeval forests occur in ancient and significantly diverged basal lineages, while those from human‐associated environments generally cluster in less differentiated domestic or mosaic groups. Basal lineages from primeval forests are usually inbred, exhibit lineage‐specific karyotypes and are partially reproductively isolated. Our results suggest that greatly diverged populations of wild S. cerevisiae exist independently of and predate domesticated isolates. We find that China harbours a reservoir of natural genetic variation of S. cerevisiae and perhaps gives an indication of the origin of the species.     

Repeated evolution of local adaptation in swimming performance: population‐level trade‐offs between burst and endurance swimming in Brachyrhaphis freshwater fish

Specialization is fundamentally important in biology because specialized traits allow species to expand into new environments, in turn promoting population differentiation and speciation. Specialization often results in trade‐offs between traits that maximize fitness in one environment but not others. Despite the ubiquity of trade‐offs, we know relatively little about how consistently trade‐offs evolve between populations when multiple sets of populations experience similarly divergent selective regimes. In the present study, we report a case study on Brachyrhaphis fishes from different predation environments. We evaluate apparent within/between population trade‐offs in burst‐speed and endurance at two levels of evolutionary diversification: high‐ and low‐predation populations of Brachyrhaphis rhabdophora, and sister species Brachyrhaphis roseni and Brachyrhaphis terrabensis, which occur in high‐ and low‐predation environments, respectively. Populations of Brachyrhaphis experiencing different predation regimes consistently evolved swimming specializations indicative of a trade‐off between two swimming forms that are likely highly adaptive in the environment in which they occur. We show that populations have become similarly locally adapted at both levels of diversification, suggesting that swimming specialization has evolved rather rapidly and persisted post‐speciation. Our findings provide valuable insight into how local adaptation evolves at different stages of evolutionary divergence.