There are currently 0 users and 42 guests online.
BMC Evolutionary Biology
The latest research articles published by BMC Evolutionary Biology
Last update31 min 38 sec ago
December 21, 2014
Phylogeny of seven Bulinus species originating from endemic areas in three African countries, in relation to the human blood fluke Schistosoma haematobium
Background: Snails species belonging to the genus Bulinus (Planorbidae) serve as intermediate host for flukes belonging to the genus Schistosoma (Digenea, Platyhelminthes). Despite its importance in the transmission of these parasites, the evolutionary history of this genus is still obscure. In the present study, we used the partial mitochondrial cytochrome oxidase subunit I (cox1) gene, and the nuclear ribosomal ITS, 18S and 28S genes to investigate the haplotype diversity and phylogeny of seven Bulinus species originating from three endemic countries in Africa (Cameroon, Senegal and Egypt). Results: The cox1 region showed much more variation than the ribosomal markers within Bulinus sequences. High levels of genetic diversity were detected at all loci in the seven studied species, with clear segregation between individuals and appearance of different haplotypes, even within same species from the same locality. Sequences clustered into two lineages; (A) groups Bulinus truncatus, B. tropicus, B. globosus and B. umbilicatus; while (B) groups B. forskalii, B. senegalensis and B. camerunensis. Interesting patterns emerge regarding schistosome susceptibility: Bulinus species with lower genetic diversity are predicted to have higher infection prevalence than those with greater diversity in host susceptibility. Conclusion: The results reported in this study are very important since a detailed understanding of the population genetic structure of Bulinus is essential to understand the epidemiology of many schistosome parasites.
The GPCR repertoire in the demosponge Amphimedon queenslandica : insights into the GPCR system at the early divergence of animals
Background: G protein-coupled receptors (GPCRs) play a central role in eukaryotic signal transduction. However, the GPCR component of this signalling system, at the early origins of metazoans is not fully understood. Here we aim to identify and classify GPCRs in Amphimedon queenslandica (sponge), a member of an earliest diverging metazoan lineage (Porifera). Furthermore, phylogenetic comparisons of sponge GPCRs with eumetazoan and bilaterian GPCRs will be essential to our understanding of the GPCR system at the roots of metazoan evolution. Results: We present a curated list of 220 GPCRs in the sponge genome after excluding incomplete sequences and false positives from our initial dataset of 282 predicted GPCR sequences obtained using Pfam search. Phylogenetic analysis reveals that the sponge genome contains members belonging to four of the five major GRAFS families including Glutamate (33), Rhodopsin (126), Adhesion (40) and Frizzled (3). Interestingly, the sponge Rhodopsin family sequences lack orthologous relationships with those found in eumetazoan and bilaterian lineages, since they clustered separately to form sponge specific groups in the phylogenetic analysis. This suggests that sponge Rhodopsins diverged considerably from that found in other basal metazoans. A few sponge Adhesions clustered basal to Adhesion subfamilies commonly found in most vertebrates, suggesting some Adhesion subfamilies may have diverged prior to the emergence of Bilateria. Furthermore, at least eight of the sponge Adhesion members have a hormone binding motif (HRM domain) in their N-termini, although hormones have yet to be identified in sponges. We also phylogenetically clarified that sponge has homologs of metabotropic glutamate (mGluRs) and GABA receptors. Conclusion: Our phylogenetic comparisons of sponge GPCRs with other metazoan genomes suggest that sponge contains a significantly diversified set of GPCRs. This is evident at the family/subfamily level comparisons for most GPCR families, in particular for the Rhodopsin family of GPCRs. In summary, this study provides a framework to perform future experimental and comparative studies to further verify and understand the roles of GPCRs that predates the divergence of bilaterian and eumetazoan lineages.
Towards a comprehensive characterization of durum wheat landraces in Moroccan traditional agrosystems: analysing genetic diversity in the light of geography, farmers¿ taxonomy and tetraploid wheat domestication history
Background: Crop diversity managed by smallholder farmers in traditional agrosystems is the outcome of historical and current processes interacting at various spatial scales, and influenced by factors such as farming practices and environmental pressures. Only recently have studies started to consider the complexity of these processes instead of simply describing diversity for breeding purposes. A first step in that aim is to add multiple references to the collection of genetic data, including the farmers? varietal taxonomy and practices and the historical background of the crop. Results: On the basis of interview data collected in a previous study, we sampled 166 populations of durum wheat varieties in two traditional Moroccan agrosystems, in the Pre-Rif and Atlas Mountains regions. Using a common garden experiment, we detected a high phenotypic variability on traits indicative of taxonomical position and breeding status, namely spike shape and plant height. Populations often combined modern (short) with traditional-like (tall) statures, and classical durum squared spike shape (5 flowers / spikelet) with flat spike shape (3 flowers/ spikelet) representative of primitive domesticated tetraploid wheat (ssp. dicoccum). By contrast, the genetic diversity assessed using 14 microsdatellite markers was relatively limited. When compared to the genetic diversity found in a large collection of tetraploid wheat, it corresponded to free-threshing tetraploid wheat. Within Morocco, the two studied regions differed for both genetic diversity and variety names. Within regions, neither geography nor variety names nor even breeding status constituted strong barriers to gene exchange despite a few significant patterns. Conclusions: This first assessment of morphological and genetic diversity allowed pointing out some important factors that may have influenced the structure and evolutionary dynamics of durum wheat in Morocco: the significance of variety names, the occurrence of mixtures within populations, the relative strength of seed exchange between farmers and local adaptation, as well as the fate of modern varieties once they have been introduced. Further, multidisciplinary studies at different spatial scales are needed to better understand these complex agrosystems of invaluable importance for food security.
Mitochondrial genomes of the Baltic clam Macoma balthica (Bivalvia: Tellinidae): setting the stage for studying mito-nuclear incompatibilities
Background: Allopatric divergence across lineages can lead to post-zygotic reproductive isolation upon secondary contact and disrupt coevolution between mitochondrial and nuclear genomes, promoting emergence of genetic incompatibilities. A previous F ST scan on the transcriptome of the Baltic clam Macoma balthica highlighted several genes potentially involved in mito-nuclear incompatibilities (MNIs). As proteins involved in the mitochondrial oxidative phosphorylation (OXPHO) chain are prone to MNIs and can contribute to the maintenance of genetic barriers, the mitochondrial genomes of six Ma. balthica individuals spanning two secondary contact zones were sequenced using the Illumina MiSeq plateform. Results: The mitogenome has an approximate length of 16,806?bp and encodes 13 protein-coding genes, 2 rRNAs and 22 tRNAs, all located on the same strand. atp8, a gene long reported as rare in bivalves, was detected. It encodes 42 amino acids and is putatively expressed and functional. A large unassigned region was identified between rrnS and tRNA Met and could likely correspond to the Control Region. Replacement and synonymous mutations were mapped on the inferred secondary structure of all protein-coding genes of the OXPHO chain. The atp6 and atp8 genes were characterized by background levels of replacement mutations, relative to synonymous mutations. However, most nad genes (notably nad2 and nad5) were characterized by an elevated proportion of replacement mutations. Conclusions: Six nearly complete mitochondrial genomes were successfully assembled and annotated, providing the necessary roadmap to study MNIs at OXPHO loci. Few replacement mutations were mapped on mitochondrial-encoded ATP synthase subunits, which is in contrast with previous data on nuclear-encoded subunits. Conversely, the high population divergence and the prevalence of non-synonymous mutations at nad genes are congruent with previous observations from the nuclear transcriptome. This further suggest that MNIs between subunits of Complex I of the OXPHO chain, coding for NADH Dehydrogenase, may play a role in maintaining barriers to gene flow in Ma. balthica.
December 20, 2014
High intralocus variability and interlocus recombination promote immunological diversity in a minimal major histocompatibility system
Background: The genes of the major histocompatibility complex (MHC/MH) have attracted considerable scientific interest due to their exceptional levels of variability and important function as part of the adaptive immune system. Despite a large number of studies on MH class II diversity of both model and non-model organisms, most research has focused on patterns of genetic variability at individual loci, failing to capture the functional diversity of the biologically active dimeric molecule. Here, we take a systematic approach to the study of MH variation, analyzing patterns of genetic variation at MH class II? and II? loci of the seahorse, which together form the immunologically active peptide binding cleft of the MH class II molecule. Results: The seahorse carries a minimal class II system, consisting of single copies of both MH class II? and II?, which are physically linked and inherited in a Mendelian fashion. Both genes are ubiquitously expressed and detectible in the brood pouch of male seahorses throughout pregnancy. Genetic variability of the two genes is high, dominated by non-synonymous variation concentrated in their peptide-binding regions. Coding variation outside these regions is negligible, a pattern thought to be driven by intra- and interlocus recombination. Despite the tight physical linkage of MH II? and II? loci, recombination has produced novel composite alleles, increasing functional diversity at sites responsible for antigen recognition. Conclusions: Antigen recognition by the adaptive immune system of the seahorse is enhanced by high variability at both MH class II? and II? loci. Strong positive selection on sites involved in pathogen recognition, coupled with high levels of intra- and interlocus recombination, produce a patchwork pattern of genetic variation driven by genetic hitchhiking. Studies focusing on variation at individual MH loci may unintentionally overlook an important component of ecologically relevant variation.
December 19, 2014
Clock model makes a large difference to age estimates of long-stemmed clades with no internal calibration: a test using Australian grasstrees
Background: Estimating divergence times in phylogenies using a molecular clock depends on accurate modeling of nucleotide substitution rates in DNA sequences. Rate heterogeneity among lineages is likely to affect estimates, especially in lineages with long stems and short crowns (?broom? clades) and no internal calibration. We evaluate the performance of the random local clocks model (RLC) and the more routinely employed uncorrelated lognormal relaxed clock model (UCLN) in situations in which a significant rate shift occurs on the stem branch of a broom clade. We compare the results of simulations to empirical results from analyses of a real rate-heterogeneous taxon ? Australian grass trees (Xanthorrhoea) ? whose substitution rate is slower than in its sister groups, as determined by relative rate tests. Results: In the simulated datasets, the RLC model performed much better than UCLN: RLC correctly estimated the age of the crown node of slow-rate broom clades, whereas UCLN estimates were consistently too young. Similarly, in the Xanthorrhoea dataset, UCLN returned significantly younger crown ages than RLC (mean estimates respectively 3?6 Ma versus 25?35 Ma). In both real and simulated datasets, Bayes Factor tests strongly favored the RLC model over the UCLN model. Conclusions: The choice of an unsuitable molecular clock model can strongly bias divergence time estimates. In particular, for data predicted to have more rate variation among than within clades, dating with RLC is much more likely to be accurate than with UCLN. The choice of clocks should be informed by the biology of the study group (e.g., life-form) or assessed with relative rate tests and post-hoc model comparisons.
Divergent evolutionary rates in vertebrate and mammalian specific conserved non-coding elements (CNEs) in echolocating mammals
Background: The majority of DNA contained within vertebrate genomes is non-coding, with a certain proportion of this thought to play regulatory roles during development. Conserved Non-coding Elements (CNEs) are an abundant group of putative regulatory sequences that are highly conserved across divergent groups and thus assumed to be under strong selective constraint. Many CNEs may contain regulatory factor binding sites, and their frequent spatial association with key developmental genes ? such as those regulating sensory system development ? suggests crucial roles in regulating gene expression and cellular patterning. Yet surprisingly little is known about the molecular evolution of CNEs across diverse mammalian taxa or their role in specific phenotypic adaptations. We examined 3,110 vertebrate-specific and ~82,000 mammalian-specific CNEs across 19 and 9 mammalian orders respectively, and tested for changes in the rate of evolution of CNEs located in the proximity of genes underlying the development or functioning of auditory systems. As we focused on CNEs putatively associated with genes underlying the development/functioning of auditory systems, we incorporated echolocating taxa in our dataset because of their highly specialised and derived auditory systems. Results: Phylogenetic reconstructions of concatenated CNEs broadly recovered accepted mammal relationships despite high levels of sequence conservation. We found that CNE substitution rates were highest in rodents and lowest in primates, consistent with previous findings. Comparisons of CNE substitution rates from several genomic regions containing genes linked to auditory system development and hearing revealed differences between echolocating and non-echolocating taxa. Wider taxonomic sampling of four CNEs associated with the homeobox genes Hmx2 and Hmx3 ? which are required for inner ear development ? revealed family-wise variation across diverse bat species. Specifically within one family of echolocating bats that utilise frequency-modulated echolocation calls varying widely in frequency and intensity high levels of sequence divergence were found. Conclusions: Levels of selective constraint acting on CNEs differed both across genomic locations and taxa, with observed variation in substitution rates of CNEs among bat species. More work is needed to determine whether this variation can be linked to echolocation, and wider taxonomic sampling is necessary to fully document levels of conservation in CNEs across diverse taxa.
December 17, 2014
Background: Synonymous codon usage bias (SCUB) is an inevitable phenomenon in organismic taxa, generally referring to differences in the occurrence frequency of codons across different species or within the genome of the same species. SCUB happens in various degrees under pressure from nature selection, mutation bias and other factors in different ways. It also attaches great significance to gene expression and species evolution, however, a systematic investigation towards the codon usage in Bombyx mori (B. mori) has not been reported yet. Moreover, it is still indistinct about the reasons contributing to the bias or the relationship between the bias and the evolution of B. mori. Results: The comparison of the codon usage pattern between the genomic DNA (gDNA) and the mitochondrial DNA (mtDNA) from B. mori suggests that mtDNA has a higher level of codon bias. Furthermore, the correspondence analysis suggests that natural selection, such as gene length, gene function and translational selection, dominates the codon preference of mtDNA, while the composition constraints for mutation bias only plays a minor role. Additionally, the clustering results of the silkworm superfamily suggest a lack of explicitness in the relationship between the codon usage of mitogenome and species evolution. Conclusions: Among the complicated influence factors leading to codon bias, natural selection is found to play a major role in shaping the high bias in the mtDNA of B. mori from our current data. Although the cluster analysis reveals that codon bias correlates little with the species evolution, furthermore, a detailed analysis of codon usage of mitogenome provides better insight into the evolutionary relationships in Lepidoptera. However, more new methods and data are needed to investigate the relationship between the mtDNA bias and evolution.
Background: There is great interest in understanding the genomic underpinnings of social evolution, in particular, the evolution of eusociality (caste-containing societies with non-reproductives that care for siblings). Subsociality is a key precursor for the evolution of eusociality and characterized by prolonged parental care and parent-offspring interaction. Here, we provide the first transcriptomic data for the small carpenter bee, Ceratina calcarata. This species is of special interest because it is subsocial and in the same family as the highly eusocial honey bee, Apis mellifera. In addition, some C. calcarata females demonstrate alloparental care without reproduction, which provides a unique opportunity to study worker behaviour in a non-eusocial species. Results: We uncovered similar gene expression patterns related to maternal care and sibling care in different groups of females. This agrees with the maternal heterochrony hypothesis, specifically, that changes in timing of offspring care gene expression are related to worker behaviour in incipient insect societies. In addition, we also detected some similarity to caste-related gene expression patterns in highly eusocial honey bees, and uncovered large lifetime changes in gene expression that accompany shifts in reproductive and maternal care behaviour. Conclusions: For Ceratina calcarata, we found that transcript expression profiles were most similar between sibling care and maternal care females. The maternal care behaviour exhibited post-reproductively by Ceratina mothers is concordant in terms of transcript expression with the alloparental care exhibited by workers. In line with theoretical predictions, our data are consistent with the maternal heterochrony hypothesis for the evolutionary development of worker behaviour in subsocial bees.
Background: In bacteria, cell size affects chromosome replication, the assembly of division machinery, cell wall synthesis, membrane synthesis and ultimately growth rate. In addition, cell size can also be a target for Darwinian evolution for protection from predators. This strong coupling of cell size and growth, however, could lead to the introduction of growth defects after size evolution. An important question remains: can bacterial cell size change and/or evolve without imposing a growth burden? Results: The directed evolution of particular cell sizes, without a growth burden, was tested with a laboratory Escherichia coli strain. Cells of defined size ranges were collected by a cell sorter and were subsequently cultured. This selection-propagation cycle was repeated, and significant changes in cell size were detected within 400 generations. In addition, the width of the size distribution was altered. The changes in cell size were unaccompanied by a growth burden. Whole genome sequencing revealed that only a few mutations in genes related to membrane synthesis conferred the size evolution. Conclusions: In conclusion, bacterial cell size could evolve, through a few mutations, without growth reduction. The size evolution without growth reduction suggests a rapid evolutionary change to diverse cell sizes in bacterial survival strategies.
December 14, 2014
Pan-African phylogeny of Mus (subgenus Nannomys ) reveals one of the most successful mammal radiations in Africa
Background: Rodents of the genus Mus represent one of the most valuable biological models for biomedical and evolutionary research. Out of the four currently recognized subgenera, Nannomys (African pygmy mice, including the smallest rodents in the world) comprises the only original African lineage. Species of this subgenus became important models for the study of sex determination in mammals and they are also hosts of potentially dangerous pathogens. Nannomys ancestors colonized Africa from Asia at the end of Miocene and Eastern Africa should be considered as the place of their first radiation. In sharp contrast with this fact and despite the biological importance of Nannomys, the specimens from Eastern Africa were obviously under-represented in previous studies and the phylogenetic and distributional patterns were thus incomplete. Results: We performed comprehensive genetic analysis of 657 individuals of Nannomys collected at approximately 300 localities across the whole sub-Saharan Africa. Phylogenetic reconstructions based on mitochondrial (CYTB) and nuclear (IRBP) genes identified five species groups and three monotypic ancestral lineages. We provide evidence for important cryptic diversity and we defined and mapped the distribution of 27 molecular operational taxonomic units (MOTUs) that may correspond to presumable species. Biogeographical reconstructions based on data spanning all of Africa modified the previous evolutionary scenarios. First divergences occurred in Eastern African mountains soon after the colonization of the continent and the remnants of these old divergences still occur there, represented by long basal branches of M. (previously Muriculus) imberbis and two undescribed species from Ethiopia and Malawi. The radiation in drier lowland habitats associated with the decrease of body size is much younger, occurred mainly in a single lineage (called the minutoides group, and especially within the species M. minutoides), and was probably linked to aridification and climatic fluctuations in middle Pliocene/Pleistocene. Conclusions: We discovered very high cryptic diversity in African pygmy mice making the genus Mus one of the richest genera of African mammals. Our taxon sampling allowed reliable phylogenetic and biogeographic reconstructions that (together with detailed distributional data of individual MOTUs) provide a solid basis for further evolutionary, ecological and epidemiological studies of this important group of rodents.
December 13, 2014
Impacts of genetic correlation on the independent evolution of body mass and skeletal size in mammals
Background: Mammals show a predictable scaling relationship between limb bone size and body mass. This relationship has a genetic basis which likely evolved via natural selection, but it is unclear how much the genetic correlation between these traits in turn impacts their capacity to evolve independently. We selectively bred laboratory mice for increases in tibia length independent of body mass, to test the hypothesis that a genetic correlation with body mass constrains evolutionary change in tibia length. Results: Over 14 generations, we produced mean tibia length increases of 9-13%, while mean body mass was unchanged, in selectively bred mice and random-bred controls. Using evolutionary scenarios with different selection and quantitative genetic parameters, we also found that this genetic correlation impedes the rate of evolutionary change in both traits, slowing increases in tibia length while preventing decreases in body mass, despite the latter’s negative effect on fitness. Conclusions: Overall, results from this ongoing selection experiment suggest that parallel evolution of relatively longer hind limbs among rodents, for example in the context of strong competition for resources and niche partitioning in heterogeneous environments, may have occurred very rapidly on geological timescales, in spite of a moderately strong genetic correlation between tibia length and body mass.
Background: Collaborative tools are of great help in conducting projects involving distant workers. Recent web technologies have helped to build such tools for jointly editing office documents and scientific data, yet none are available for handling phylogenies. Though a large number of studies and projects in evolutionary biology and systematics involve collaborations between scientists of different institutes, current tree comparison visualization software and websites are directed toward single-user access. Moreover, tree comparison functionalities are dispersed between different software that mainly focus on high level single tree visualization but to the detriment of basic tree comparison features. Results: The web platform presented here, named CompPhy, intends to fill this gap by allowing collaborative work on phylogenies and by gathering simple advanced tools dedicated to tree comparison. It offers functionalities for tree edition, tree comparison, supertree inference and data management in a collaborative environment. The latter aspect is a specific feature of the platform, allowing people located in different places to work together at the same time on a common project. CompPhy thus proposes shared tree visualization, both synchronous and asynchronous tree manipulation, data exchange/storage, as well as facilities to keep track of the progress of analyses in working sessions. Specific advanced comparison tools are also available, such as consensus and supertree inference, or automated branch swaps of compared trees. As projects can be readily created and shared, CompPhy is also a tool that can be used easily to interact with students in a educational setting, either in the classroom or for assignments. Conclusions: CompPhy is the first web platform devoted to the comparison of phylogenetic trees allowing real-time distant collaboration on a phylogenetic/phylogenomic project. This application can be accessed freely with a recent browser at the following page of the ATGC bioinformatics platform: http://www.atgc-montpellier.fr/compphy/.
December 12, 2014
Background: Sex chromosomes exhibit many unusual patterns in sequence and gene expression relative to autosomes. Birds have evolved a female heterogametic sex system (male ZZ, female ZW), through stepwise suppression of recombination between chrZ and chrW. To address the broad patterns and complex driving forces of Z chromosome evolution, we analyze here 45 newly available bird genomes and four species’ transcriptomes, over their course of recombination loss between the sex chromosomes. Results: We show Z chromosomes in general have a significantly higher substitution rate in introns and synonymous protein-coding sites than autosomes, driven by the male-to-female mutation bias (‘male-driven evolution’ effect). Our genome-wide estimate reveals that the degree of such a bias ranges from 1.6 to 3.8 among different species. G + C content of third codon positions exhibits the same trend of gradual changes with that of introns, between chrZ and autosomes or regions with increasing ages of becoming Z-linked, therefore codon usage bias in birds is probably driven by the mutational bias. On the other hand, Z chromosomes also evolve significantly faster at nonsynonymous sites relative to autosomes (‘fast-Z’ evolution). And species with a lower level of intronic heterozygosities tend to evolve even faster on the Z chromosome. Further analysis of fast-evolving genes’ enriched functional categories and sex-biased expression patterns support that, fast-Z evolution in birds is mainly driven by genetic drift. Finally, we show in species except for chicken, gene expression becomes more male-biased within Z-linked regions that have became hemizygous in females for a longer time, suggesting a lack of global dosage compensation in birds, and the reported regional dosage compensation in chicken has only evolved very recently. Conclusions: In conclusion, we uncover that the sequence and expression patterns of Z chromosome genes covary with their ages of becoming Z-linked. In contrast to the mammalian X chromosomes, such patterns are mainly driven by mutational bias and genetic drift in birds, due to the opposite sex-biased inheritance of Z vs. X.
December 11, 2014
Dynamic evolution of the alpha (¿) and beta (ß) keratins has accompanied integument diversification and the adaptation of birds into novel lifestyles
Background: Vertebrate skin appendages are constructed of keratins produced by multigene families. Alpha (α) keratins are found in all vertebrates, while beta (β) keratins are found exclusively in reptiles and birds. We have studied the molecular evolution of these gene families in the genomes of 48 phylogenetically diverse birds and their expression in the scales and feathers of the chicken. Results: We found that the total number of α-keratins is lower in birds than mammals and non-avian reptiles, yet two α-keratin genes (KRT42 and KRT75) have expanded in birds. The β-keratins, however, demonstrate a dynamic evolution associated with avian lifestyle. The avian specific feather β-keratins comprise a large majority of the total number of β-keratins, but independently derived lineages of aquatic and predatory birds have smaller proportions of feather β-keratin genes and larger proportions of keratinocyte β-keratin genes. Additionally, birds of prey have a larger proportion of claw β-keratins. Analysis of α- and β-keratin expression during development of chicken scales and feathers demonstrates that while α-keratins are expressed in these tissues, the number and magnitude of expressed β-keratin genes far exceeds that of α-keratins. Conclusions: These results support the view that the number of α- and β-keratin genes expressed, the proportion of the β-keratin subfamily genes expressed and the diversification of the β-keratin genes have been important for the evolution of the feather and the adaptation of birds into multiple ecological niches.
December 9, 2014
The role of aedeagus size and shape in failed mating interactions among recently diverged taxa in the Drosophila mojavensis species cluster
Background: Investigating the evolution of species-specific insect genitalia is central to understanding how morphological diversification contributes to reproductive isolation and lineage divergence. While many studies evoke some form of sexual selection to explain genitalia diversity, the basis of selection and the mechanism of heterospecific mate exclusion remains vague. I conducted reciprocal mate pair trials in the Drosophila mojavensis species cluster to quantify the frequency of failed insemination attempts, historically referred to as pseudocopulation, between lineages with discrete size and shape differences of the male aedeagus. Results: In cross-taxon matings aedeagus size had a significant effect on pseudocopulation frequencies, while aedeagus shape and genetic distance did not. The direction of the size difference was an important factor for successful mating. When females were mated to a cross-taxon male with a larger aedeagus than males from her own species, the pair could not establish a successful mating interaction. Females mated to cross-taxon males with a smaller aedeagus than conspecific males were able to establish the mating interaction but had issues disengaging at the end of the interaction. Conclusions: The results of this study support a role for aedeagus size in the male-female mating interaction, with a secondary role for aedeagus shape. In natural populations, mating failure based on aedeagus size could serve as an important reproductive isolating mechanism resulting in failed insemination attempts after both the male and female show a willingness to mate.
Rapid evolutionary responses of life history traits to different experimentally-induced pollutions in Caenorhabditis elegans
Background: Anthropogenic disturbances can lead to intense selection pressures on traits and very rapid evolutionary changes. Evolutionary responses to environmental changes, in turn, reflect changes in the genetic structure of the traits, accompanied by a reduction of evolutionary potential of the populations under selection. Assessing the effects of pollutants on the evolutionary responses and on the genetic structure of populations is thus important to understanding the mechanisms that entail specialization to novel environmental conditions or resistance to novel stressors. Results: Using an experimental evolution approach we exposed Caenorhabditis elegans populations to uranium, salt and alternating uranium-salt environments over 22 generations. We analyzed the changes in the average values of life history traits and the consequences at the demographic level in these populations. We also estimated the phenotypic and genetic (co)variance structure of these traits at different generations. Compared to populations in salt, populations in uranium showed a reduction of the stability of their trait structure and a higher capacity to respond by acclimation. However, the evolutionary responses of traits were generally lower for uranium compared to salt treatment; and the evolutionary responses to the alternating uranium?salt environment were between those of constant environments. Consequently, at the end of the experiment, the population rate of increase was higher in uranium than in salt and intermediate in the alternating environment. Conclusions: Our multigenerational experiment confirmed that rapid adaptation to different polluted environments may involve different evolutionary responses resulting in demographic consequences. These changes are partly explained by the effects of the pollutants on the genetic (co)variance structure of traits and the capacity of acclimation to novel conditions. Finally, our results in the alternating environment may confirm the selection of a generalist type in this environment.
December 4, 2014
Background: The allele frequency spectrum (AFS) consists of counts of the number of single nucleotide polymorphism (SNP) loci with derived variants present at each given frequency in a sample. Multiple approaches have recently been developed for parameter estimation and calculation of model likelihoods based on the joint AFS from two or more populations. We conducted a simulation study of one of these approaches, implemented in the Python module ?a?i, to compare parameter estimation and model selection accuracy given different sample sizes under one- and two-population models. Results: Our simulations included a variety of demographic models and two parameterizations that differed in the timing of events (divergence or size change). Using a number of SNPs reasonably obtained through next-generation sequencing approaches (10,000 - 50,000), accurate parameter estimates and model selection were possible for models with more ancient demographic events, even given relatively small numbers of sampled individuals. However, for recent events, larger numbers of individuals were required to achieve accuracy and precision in parameter estimates similar to that seen for models with older divergence or population size changes. We quantify i) the uncertainty in model selection, using tools from information theory, and ii) the accuracy and precision of parameter estimates, using the root mean squared error, as a function of the timing of demographic events, sample sizes used in the analysis, and complexity of the simulated models. Conclusions: Here, we illustrate the utility of the genome-wide AFS for estimating demographic history and provide recommendations to guide sampling in population genomics studies that seek to draw inference from the AFS. Our results indicate that larger samples of individuals (and thus larger AFS) provide greater power for model selection and parameter estimation for more recent demographic events.
Body size and allometric shape variation in the molly Poecilia vivipara along a gradient of salinity and predation
Background: Phenotypic diversity among populations may result from divergent natural selection acting directly on traits or via correlated responses to changes in other traits. One of the most frequent patterns of correlated response is the proportional change in the dimensions of anatomical traits associated with changes in growth or absolute size, known as allometry. Livebearing fishes subject to predation gradients have been shown to repeatedly evolve larger caudal peduncles and smaller cranial regions under high predation regimes. Poecilia vivipara is a livebearing fish commonly found in coastal lagoons in the north of the state of Rio de Janeiro, Brazil. Similar to what is observed in other predation gradients, lagoons inhabited by P. vivipara vary in the presence of piscivorous fishes; contrary to other poeciliid systems, populations of P. vivipara vary greatly in body size, which opens the possibility of strong allometric effects on shape variation. Here we investigated body shape diversification among six populations of P. vivipara along a predation gradient and its relationship with allometric trajectories within and among populations. Results: We found substantial body size variation and correlated shape changes among populations. Multivariate regression analysis showed that size variation among populations accounted for 66% of shape variation in females and 38% in males, suggesting that size is the most important dimension underlying shape variation among populations of P. vivipara in this system. Changes in the relative sizes of the caudal peduncle and cranial regions were only partly in line with predictions from divergent natural selection associated with predation regime. Conclusions: Our results suggest the possibility that adaptive shape variation among populations has been partly constrained by allometry in P. vivipara. Processes governing body size changes are therefore important in the diversification of this species. We conclude that in species characterized by substantial among-population differences in body size, ignoring allometric effects when investigating divergent natural selection?s role in phenotypic diversification might not be warranted.
Empirical demonstration of environmental sensing in catalytic RNA: evolution of interpretive behavior at the origins of life
Background: The origins of life on the Earth required chemical entities to interact with their environments in ways that could respond to natural selection. The concept of interpretation, where biotic entities use signs in their environment as proxy for the existence of other items of selective value in their environment, has been proposed on theoretical grounds to be relevant to the origins and early evolution of life. However this concept has not been demonstrated empirically. Results: Here, we present data that certain catalytic RNA sequences have properties that would enable interpretation of divalent cation levels in their environment. By assaying the responsiveness of two variants of the Tetrahymena ribozyme to the Ca2+ ion as a sign for the more catalytically useful Mg2+ ion, we show an empirical proof-of-principle that interpretation can be an evolvable trait in RNA, often suggested as a model system for early life. In particular we demonstrate that in vitro, the wild-type version of the Tetrahymena ribozyme is not interpretive, in that it cannot use Ca2+ as a sign for Mg2+. Yet a variant of this sequence containing five mutations that alter its ability to utilize the Ca2+ ion engenders a strong interpretive characteristic in this RNA. Conclusions: We have shown that RNA molecules in a test tube can meet the minimum criteria for the evolution of interpretive behaviour in regards to their responses to divalent metal ion concentrations in their environment. Interpretation in RNA molecules provides a property entirely dependent on natural physico-chemical interactions, but capable of shaping the evolutionary trajectory of macromolecules, especially in the earliest stages of life’s history.
The Genealogical World of Phylogenetic Networks
BMC Evolutionary Biology