Accumulation of point mutations and reassortment of genomic RNA segments are involved in the microevolution of Puumala hantavirus in a bank vole (Myodes glareolus) population

"The genetic diversity of Puumala hantavirus (PUUV) was studied in a local population of its natural host, the bank vole (Myodes glareolus). The trapping area (2.5x2.5 km) at Konnevesi, Central Finland, included 14 trapping sites, at least 500 m apart; altogether, 147 voles were captured during May and October 2005. Partial sequences of the S, M and L viral genome segments were recovered from 40 animals. Seven, 12 and 17 variants were detected for the S, M and L sequences, respectively; these represent new wild-type PUUV strains that belong to the Finnish genetic lineage. The genetic diversity of PUUV strains from Konnevesi was 0.2-4.9% for the S segment, 0.2-4.8% for the M segment and 0.2-9.7% for the L segment. Most nucleotide substitutions were synonymous and most deduced amino acid substitutions were conservative, probably due to strong stabilizing selection operating at the protein level. Based on both sequence markers and phylogenetic clustering, the S, M and L sequences could be assigned to two groups, 'A' and 'B'. Notably, not all bank voles carried S, M and L sequences belonging to the same group, i.e. SAMALA or SBMBLB.. A substantial proportion (8/40, 20%) of the newly characterized PUUV strains possessed reassortant genomes such as SBMALA, SAMBLB or SBMALB. These results suggest that at least some of the PUUV reassortants are viable and can survive in the presence of their parental strains."

Analysis of the Genome and Transcriptome of Cryptococcus neoformans var. grubii Reveals Complex RNA Expression and Microevolution Leading to Virulence Attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation.

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Neutral Genomic Microevolution of a Recently Emerged Pathogen, Salmonella enterica Serovar Agona

Salmonella enterica serovar Agona has caused multiple food-borne outbreaks of gastroenteritis since it was first isolated in 1952. We analyzed the genomes of 73 isolates from global sources, comparing five distinct outbreaks with sporadic infections as well as food contamination and the environment. Agona consists of three lineages with minimal mutational diversity: only 846 single nucleotide polymorphisms (SNPs) have accumulated in the non-repetitive, core genome since Agona evolved in 1932 and subsequently underwent a major population expansion in the 1960s. Homologous recombination with other serovars of S. enterica imported 42 recombinational tracts (360 kb) in 5/143 nodes within the genealogy, which resulted in 3,164 additional SNPs. In contrast to this paucity of genetic diversity, Agona is highly diverse according to pulsed-field gel electrophoresis (PFGE), which is used to assign isolates to outbreaks. PFGE diversity reflects a highly dynamic accessory genome associated with the gain or loss (indels) of 51 bacteriophages, 10 plasmids, and 6 integrative conjugational elements (ICE/IMEs), but did not correlate uniquely with outbreaks. Unlike the core genome, indels occurred repeatedly in independent nodes (homoplasies), resulting in inaccurate PFGE genealogies. The accessory genome contained only few cargo genes relevant to infection, other than antibiotic resistance. Thus, most of the genetic diversity within this recently emerged pathogen reflects changes in the accessory genome, or is due to recombination, but these changes seemed to reflect neutral processes rather than Darwinian selection. Each outbreak was caused by an independent clade, without universal, outbreak-associated genomic features, and none of the variable genes in the pan-genome seemed to be associated with an ability to cause outbreaks. © 2013 Achtman et al

The changing pace of insular life: 5000 years of microevolution in the Orkney vole (Microtus arvalis orcadensis)

Island evolution may be expected to involve fast initial morphological divergence followed by stasis. We tested this model using the dental phenotype of modern and ancient common voles (Microtus arvalis), introduced onto the Orkney archipelago (Scotland) from continental Europe some 5000 years ago. First, we investigated phenotypic divergence of Orkney and continental European populations and assessed climatic influences. Second, phenotypic differentiation among Orkney populations was tested against geography, time, and neutral genetic patterns. Finally, we examined evolutionary change along a time series for the Orkney Mainland. Molar gigantism and anterior-lobe hypertrophy evolved rapidly in Orkney voles following introduction, without any transitional forms detected. Founder events and adaptation appear to explain this initial rapid evolution. Idiosyncrasy in dental features among different island populations of Orkney voles is also likely the result of local founder events following Neolithic translocation around the archipelago. However, against our initial expectations, a second marked phenotypic shift occurred between the 4th and 12th centuries AD, associated with increased pastoral farming and introduction of competitors (mice and rats) and terrestrial predators (foxes and cats). These results indicate that human agency can generate a more complex pattern of morphological evolution than might be expected in island rodents.

'Nothing exists for one cause' : putting biology back into evolution

The opening phrase of the title is from Charles Darwin’s notebooks (Schweber 1977). It is a double reminder, firstly that mainstream evolutionary theory is not just about describing nature but is particularly looking for mechanisms or ‘causes’, and secondly, that there will usually be several causes affecting any particular outcome. The second part of the title is our concern at the almost universal rejection of the idea that biological mechanisms are sufficient for macroevolutionary changes, thus rejecting a cornerstone of Darwinian evolutionary theory. Our primary aim here is to consider ways of making it easier to develop and to test hypotheses about evolution. Formalizing hypotheses can help generate tests. In an absolute sense, some of the discussion by scientists about evolution is little better than the lack of reasoning used by those advocating intelligent design. Our discussion here is in a Popperian framework where science is defined by that area of study where it is possible, in principle, to find evidence against hypotheses – they are in principle falsifiable. However, with time, the boundaries of science keep expanding. In the past, some aspects of evolution were outside the current boundaries of falsifiable science, but increasingly new techniques and ideas are expanding the boundaries of science and it is appropriate to re-examine some topics. It often appears that over the last few decades there has been an increasingly strong assumption to look first (and only) for a physical cause. This decision is virtually never formally discussed, just an assumption is made that some physical factor ‘drives’ evolution. It is necessary to examine our assumptions much more carefully. What is meant by physical factors ‘driving’ evolution, or what is an ‘explosive radiation’. Our discussion focuses on two of the six mass extinctions, the fifth being events in the Late Cretaceous, and the sixth starting at least 50,000 years ago (and is ongoing). Cretaceous/Tertiary boundary; the rise of birds and mammals. We have had a long-term interest (Cooper and Penny 1997) in designing tests to help evaluate whether the processes of microevolution are sufficient to explain macroevolution. The real challenge is to formulate hypotheses in a testable way. For example the numbers of lineages of birds and mammals that survive from the Cretaceous to the present is one test. Our first estimate was 22 for birds, and current work is tending to increase this value. This still does not consider lineages that survived into the Tertiary, and then went extinct later. Our initial suggestion was probably too narrow in that it lumped four models from Penny and Phillips (2004) into one model. This reduction is too simplistic in that we need to know about survival and ecological and morphological divergences during the Late Cretaceous, and whether Crown groups of avian or mammalian orders may have existed back into the Cretaceous. More recently (Penny and Phillips 2004) we have formalized hypotheses about dinosaurs and pterosaurs, with the prediction that interactions between mammals (and groundfeeding birds) and dinosaurs would be most likely to affect the smallest dinosaurs, and similarly interactions between birds and pterosaurs would particularly affect the smaller pterosaurs. There is now evidence for both classes of interactions, with the smallest dinosaurs and pterosaurs declining first, as predicted. Thus, testable models are now possible. Mass extinction number six: human impacts. On a broad scale, there is a good correlation between time of human arrival, and increased extinctions (Hurles et al. 2003; Martin 2005; Figure 1). However, it is necessary to distinguish different time scales (Penny 2005) and on a finer scale there are still large numbers of possibilities. In Hurles et al. (2003) we mentioned habitat modification (including the use of Geogenes III July 2006 31 fire), introduced plants and animals (including kiore) in addition to direct predation (the ‘overkill’ hypothesis). We need also to consider prey switching that occurs in early human societies, as evidenced by the results of Wragg (1995) on the middens of different ages on Henderson Island in the Pitcairn group. In addition, the presence of human-wary or humanadapted animals will affect the distribution in the subfossil record. A better understanding of human impacts world-wide, in conjunction with pre-scientific knowledge will make it easier to discuss the issues by removing ‘blame’. While continued spontaneous generation was accepted universally, there was the expectation that animals continued to reappear. New Zealand is one of the very best locations in the world to study many of these issues. Apart from the marine fossil record, some human impact events are extremely recent and the remains less disrupted by time.

4000 years of phenotypic change in an island bird: heterogeneity of selection over three microevolutionary timescales.

Pronounced phenotypic shifts in island populations are typically attributed to natural selection, but reconstructing heterogeneity in long-term selective regimes remains a challenge. We examined a scenario of divergence proposed for species colonizing a new environment, involving directional selection with a rapid shift to a new optimum and subsequent stabilization. We provide some of the first empirical evidence for this model of evolution using morphological data from three timescales in an island bird, Zosterops lateralis chlorocephalus. In less than four millennia since separation from its mainland counterpart, a substantial increase in body size has occurred and was probably achieved in fewer than 500 generations after colonization. Over four recent decades, morphological traits have fluctuated in size but showed no significant directional trends, suggesting maintenance of a relatively stable phenotype. Finally, estimates of contemporary selection gradients indicated generally weak directional selection. These results provide a rare description of heterogeneity in long-term natural regimes, and caution that observations of current selection may be of limited value in inferring mechanisms of past adaptation due to a lack of constancy even over short time-frames.

Recent emergence of dengue virus serotype 4 in French Polynesia results from multiple introductions from other south pacific islands

BACKGROUND: Infection by dengue virus (DENV) is a major public health concern in hundreds of tropical and subtropical countries. French Polynesia (FP) regularly experiences epidemics that initiate, or are consecutive to, DENV circulation in other South Pacific Island Countries (SPICs). In January 2009, after a decade of serotype 1 (DENV-1) circulation, the first cases of DENV-4 infection were reported in FP. Two months later a new epidemic emerged, occurring about 20 years after the previous circulation of DENV-4 in FP. In this study, we investigated the epidemiological and molecular characteristics of the introduction, spread and genetic microevolution of DENV-4 in FP. METHODOLOGY/PRINCIPAL FINDINGS: Epidemiological data suggested that recent transmission of DENV-4 in FP started in the Leeward Islands and this serotype quickly displaced DENV-1 throughout FP. Phylogenetic analyses of the nucleotide sequences of the envelope (E) gene of 64 DENV-4 strains collected in FP in the 1980s and in 2009-2010, and some additional strains from other SPICs showed that DENV-4 strains from the SPICs were distributed into genotypes IIa and IIb. Recent FP strains were distributed into two clusters, each comprising viruses from other but distinct SPICs, suggesting that emergence of DENV-4 in FP in 2009 resulted from multiple introductions. Otherwise, we observed that almost all strains collected in the SPICs in the 1980s exhibit an amino acid (aa) substitution V287I within domain I of the E protein, and all recent South Pacific strains exhibit a T365I substitution within domain III. CONCLUSIONS/SIGNIFICANCE: This study confirmed the cyclic re-emergence and displacement of DENV serotypes in FP. Otherwise, our results showed that specific aa substitutions on the E protein were present on all DENV-4 strains circulating in SPICs. These substitutions probably acquired and subsequently conserved could reflect a founder effect to be associated with epidemiological, geographical, eco-biological and social specificities in SPICs.

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.

Transposable elements in a marginal plant population

Background
How new forms arise in nature has engaged evolutionary biologists since Darwin's seminal treatise on the origin of species. Transposable elements (TEs) may be among the most important internal sources for intraspecific variability. Thus, we aimed to explore the temporal dynamics of several TEs in individual genotypes from a small, marginal population of Aegilops speltoides. A diploid cross-pollinated grass species, it is a wild relative of the various wheat species known for their large genome sizes contributed by an extraordinary number of TEs, particularly long terminal repeat (LTR) retrotransposons. The population is characterized by high heteromorphy and possesses a wide spectrum of chromosomal abnormalities including supernumerary chromosomes, heterozygosity for translocations, and variability in the chromosomal position or number of 45S and 5S ribosomal DNA (rDNA) sites. We propose that variability on the morphological and chromosomal levels may be linked to variability at the molecular level and particularly in TE proliferation.

Results
Significant temporal fluctuation in the copy number of TEs was detected when processes that take place in small, marginal populations were simulated. It is known that under critical external conditions, outcrossing plants very often transit to self-pollination. Thus, three morphologically different genotypes with chromosomal aberrations were taken from a wild population of Ae. speltoides, and the dynamics of the TE complex traced through three rounds of selfing. It was discovered that: (i) various families of TEs vary tremendously in copy number between individuals from the same population and the selfed progenies; (ii) the fluctuations in copy number are TE-family specific; (iii) there is a great difference in TE copy number expansion or contraction between gametophytes and sporophytes; and (iv) a small percentage of TEs that increase in copy number can actually insert at novel locations and could serve as a bona fide mutagen.

Conclusions
We hypothesize that TE dynamics could promote or intensify morphological and karyotypical changes, some of which may be potentially important for the process of microevolution, and allow species with plastic genomes to survive as new forms or even species in times of rapid climatic change.

Mechanisms of azole resistance and carcinogenic acetaldehyde production in chronic mucosal candidosis

Autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy (APECED, APS1) is an autoimmune disease caused by a loss-of function mutation in the autoregulator gene (AIRE). Patients with APECED suffer from chronic mucocutaneous candidosis (CMC) of the oral cavity and oesophagus often since early childhood. The patients are mainly colonized with Candida albicans and decades of exposure to antifungal agents have lead to the development of clinical and microbiological resistance in the treatment of CMC in the APECED patient population in Finland. A high incidence of oral squamous cell carcinoma is associated with oral CMC lesions in the APECED patients over the age of 25. The overall aim of this study was firstly, to investigate the effect of long-term azole exposure on the metabolism of oral C. albicans isolates from APECED patients with CMC and secondly, to analyse the specific molecular mechanisms that are responsible for these changes. The aim of the first study was to examine C. albicans strains from APECED patients and the level of cross-resistance to miconazole, the recommended topical compound for the treatment of oral candidosis. A total of 16% of the strains had decreased susceptibility to miconazole and all of these isolates had decreased susceptibility to fluconazole. Miconazole MICs also correlated with MICs to voriconazole and posaconazole. A significant positive correlation between the years of miconazole exposure and the MICs to azole antifungal agents was also found. These included azoles the patients had not been exposed to. The aim of our second study was to determine if the APECED patients are continuously colonized with the same C. albicans strains despite extensive antifungal treatment and to gain a deeper insight into the genetic changes leading to azole resistance. The strains were typed using MLST and our results confirmed that all patients were persistently colonized with the same or a genetically related strain despite antifungal treatment between isolations. No epidemic strains were found. mRNA expression was analysed by Northern blotting, protein level by western blotting, and TAC1 and ERG11 genes were sequenced. The main molecular mechanisms resulting in azole resistance were gain-of-function mutations in TAC1 leading to over expression of CDR1 and CDR2, genes linked to azole resistance. Several strains had also developed point mutations in ERG11, another gene linked to azole resistance. In the third study we used gas chromatography to test whether the level of carcinogenic acetaldehyde produced by C. albicans strains isolated from APECED patients were different from the levels produced by strains isolated from healthy controls and oral carcinoma patients. Acetaldehyde is a carcinogenic product of alcohol fermentation and metabolism in microbes associated with cancers of the upper digestive tract. In yeast, acetaldehyde is a by-product of the pyruvate bypass that converts pyruvate into acetyl-CoA during fermentation. Our results showed that strains isolated from APECED patients produced mutagenic levels of acetaldehyde in the presence of glucose (100mM, 18g/l) and the levels produced were significantly higher than those from strains isolated from controls and oral carcinoma patients. All strains in the study, however, were found to produce mutagenic levels of acetaldehyde in the presence of ethanol (11mM). The glucose and ethanol levels used in this study are equivalent to those found in food and beverages and our results highlight the role of dietary sugars and ethanol on carcinogenesis. The aims of our fourth study were to research the effect of growth conditions in the levels of acetaldehyde produced by C. albicans and to gain deeper insight into the role of different genes in the pyruvate-bypass in the production of high acetaldehyde levels. Acetaldehyde production in the presence of glucose increased by 17-fold under moderately hypoxic conditions compared to the levels produced under normoxic conditions. Under moderately hypoxic conditions acetaldehyde levels did not correlate with the expression of ADH1 and ADH2, genes catalyzing the oxidation of ethanol to acetaldehyde, or PDC11, the gene catalyzing the oxidation of pyruvate to acetaldehyde but correlated with the expression of down-stream genes ALD6 and ACS1. Our results highlight a problem where indiscriminate use of azoles may influence azole susceptibility and lead to the development of cross-resistance. Despite clinically successful treatment leading to relief of symptoms, colonization by C. albicans strains is persistent within APECED patients. Microevolution and point mutations that occur in strains may lead to the development of azole-resistant isolates and metabolic changes leading to increased production of carcinogenic acetaldehyde.

Temporal Stability of Molecular Diversity Measures in Natural Populations of Drosophila pseudoobscura and Drosophila persimilis

Many molecular ecological and evolutionary studies sample wild populations at a single point in time, failing to consider that data they collect represents genetic variation from a potentially unrepresentative snapshot in time. Variation across time in genetic parameters may occur quickly in species that produce multiple generations of offspring per year. However, many studies of rapid contemporary microevolution examine phenotypic trait divergence as opposed to molecular evolutionary divergence. Here, we compare genetic diversity in wild caught populations of Drosophila persimilis and D. pseudoobscura collected 16 years apart at the same time of year and same site at four X-linked and two mitochondrial loci to assess genetic stability. We found no major changes in nucleotide diversity in either species, but we observed a drastic shift in Tajima’s D between D. pseudoobscura timepoints at one locus associated with the increased abundance of a set of related haplotypes. Our data also suggests that D. persimilis may have recently accelerated its demographic expansion. While the changes we observed were modest, this study reinforces the importance of considering potential temporal variation in genetic parameters within single populations over short evolutionary timescales.

Strain variation in ppGpp concentration and RpoS levels in laboratory strains of Escherichia coli K-12

Laboratory strains and natural isolates of Escherichia coli differ in their level of stress resistance due to strain variation in the level of the sigma factor sigma(S) (or RpoS), the transcriptional master controller of the general stress response. We found that the high level of RpoS in one laboratory strain (MC4100) was partially dependent on an elevated basal level of ppGpp, an alarmone responding to stress and starvation. The elevated ppGpp was caused by two mutations in spoT, a gene associated with ppGpp synthesis and degradation. The nature of the spoT allele influenced the level of ppGpp in both MC4100 and another commonly used K-12 strain, MG1655. Introduction of the spoT mutation into MG1655 also resulted in an increased level of RpoS, but the amount of RpoS was lower in MG1655 than in MC4100 with either the wild-type or mutant spoT allele. In both MC4100 and MG1655, high ppGpp concentration increased RpoS levels, which in turn reduced growth with poor carbon sources like acetate. The growth inhibition resulting from elevated ppGpp was relieved by rpoS mutations. The extent of the growth inhibition by ppGpp, as well as the magnitude of the relief by rpoS mutations, differed between MG1655 and MC4100. These results together suggest that spoT mutations represent one of several polymorphisms influencing the strain variation of RpoS levels. Stress resistance was higher in strains with the spoT mutation, which is consistent with the conclusion that microevolution affecting either or both ppGpp and RpoS can reset the balance between self-protection and nutritional capability, the SPANC balance, in individual strains of E coli.