A tale of two lineages: unexpected,long‐term persistence of the amphibian‐killing fungus in Brazil

For the past 17 years, scientists have been compiling a list of amphibian species susceptible to infection by the amphibian‐killing chytrid fungus, Batrachochytrium dendrobatidis (Bd), all over the world, with >500 species infected on every continent except Antarctica (Olson et al. 2013 ). Where Bd has been found, the impacts on amphibians has been one of two types: either Bd arrives into a naïve amphibian population followed by a mass die‐off and population declines (e.g. Lips et al. 2006 ), or Bd is present at some moderate prevalence, usually infecting many species but at apparently nonlethal intensities for a long time. In this issue of Molecular Ecology, Rodriguez et al. ( 2014 ) discover that the Atlantic Coastal Forest of Brazil is home to two Bd lineages: the Global Pandemic Lineage (Bd‐GPL) – the strain responsible for mass die‐offs and population declines – and a lineage endemic to Brazil (Bd‐Bz). Even more surprising was that both lineages have been present in this area for the past 100 years, making these the oldest records of Bd infecting amphibians. The team also described a moderate but steady prevalence of ~20% across all sampled anuran families for over 100 years, indicating that Brazil has been in an enzootic disease state for over a century. Most amphibians were infected with Bd‐GPL, suggesting this lineage may be a better competitor than Bd‐Bz or may be replacing the Bd‐Bz lineage. Rodriguez et al. ( 2014 ) also detected likely hybridization of the two Bd lineages, as originally described by Schloegel et al. ( 2012 ).     

Demystifying the RAD fad

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart ( 2014 ) entitled ‘Recent novel approaches for population genomics data analysis’. Restriction‐site‐associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. ( 2007 ) and Baird et al. ( 2008 ) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. 2011 ), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart 2014 ). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. 2007 ; Baird et al. 2008 ), double digest RAD (ddRAD, Peterson et al. 2012 ), ezRAD (Toonen et al. 2013 ) and 2bRAD (Wang et al. 2012 ). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

Homomorphic plant sex chromosomes are coming of age

Sex chromosomes are a very peculiar part of the genome that have evolved independently in many groups of animals and plants (Bull 1983 ). Major research efforts have so far been focused on large heteromorphic sex chromosomes in a few animal and plant species (Chibalina & Filatov 2011 ; Zhou & Bachtrog 2012 ; Bellott et al. 2014 ; Hough et al. 2014 ; Zhou et al. 2014 ), while homomorphic (cytologically indistinguishable) sex chromosomes have largely been neglected. However, this situation is starting to change. In this issue, Geraldes et al. ( 2015 ) describe a small (~100 kb long) sex‐determining region on the homomorphic sex chromosomes of poplars (Populus trichocarpa and related species, Fig.  1 ). All species in Populus and its sister genus Salix are dioecious, suggesting that dioecy and the sex chromosomes, if any, should be relatively old. Contrary to this expectation, Geraldes et al. ( 2015 ) demonstrate that the sex‐determining region in poplars is of very recent origin and probably evolved within the genus Populus only a few million years ago.     

Next‐generation sequencing for molecular ecology: a caveat regarding pooled samples

We develop a model based on the Dirichlet‐compound multinomial distribution (CMD) and Ewens sampling formula to predict the fraction of SNP loci that will appear fixed for alternate alleles between two pooled samples drawn from the same underlying population. We apply this model to next‐generation sequencing (NGS) data from Baltic Sea herring recently published by (Corander et al., 2013 , Molecular Ecology, 2931 –2940), and show that there are many more fixed loci than expected in the absence of genetic structure. However, we show through coalescent simulations that the degree of population structure required to explain the fraction of alternatively fixed SNPs is extraordinarily high and that the surplus of fixed loci is more likely a consequence of limited representation of individual gene copies in the pooled samples, than it is of population structure. Our analysis signals that the use of NGS on pooled samples to identify divergent SNPs warrants caution. With pooled samples, it is hard to diagnose when an NGS experiment has gone awry; especially when NGS data on pooled samples are of low read depth with a limited number of individuals, it may be worthwhile to temper claims of unexpected population differentiation from pooled samples, pending verification with more reliable methods or stricter adherence to recommended sampling designs for pooled sequencing e.g. Futschik & Schlötterer 2010 , Genetics, 186 , 207; Gautier et al., 2013a , Molecular Ecology, 3766 –3779). Analysis of the data and diagnosis of problems is easier and more reliable (and can be less costly) with individually barcoded samples. Consequently, for some scenarios, individual barcoding may be preferable to pooling of samples.     

Don't throw the baby out with the bathwater: identifying and mapping paralogs in salmonids

Many eukaryotic genomes contain a large fraction of gene duplicates (or paralogs) as a result of ancient or recent whole‐genome duplications (Ohno 1970 ; Jaillon et al. 2004 ; Kellis et al. 2004 ). Identifying paralogs with NGS data is a pervasive problem in both ancient polyploids and neopolyploids. Likewise, paralogs are often treated as a nuisance that has to be detected and removed (Everett et al. 2012 ). In this issue of Molecular Ecology Resources, Waples et al. ( 2015 ) show that exclusion might not be necessary and how we may miss out on important genomic information in doing so. They present a novel statistical approach to detect paralogs based on the segregation of RAD loci in haploid offspring and test their method by constructing linkage maps with and without these duplicated loci in chum salmon, Oncorhynchus keta (Fig.  1 ). Their linkage map including the resolved paralogs shows that these are mostly located in the distal regions of several linkage groups. Particularly intriguing is their finding that these homoeologous regions appear impoverished in transposable elements (TE). Given the role that TE play in genome remodelling, it is noteworthy that these elements are of low abundance in regions showing residual tetrasomic inheritance. This raises the question whether re‐diploidization is constrained in these regions and whether they might have a role to play in salmonid speciation. This study provides an original approach to identifying duplicated loci in species with a pedigree, as well as providing a dense linkage map for chum salmon, and interesting insights into the retention of gene duplicates in an ancient polyploid.     

Population genomic analyses reveal possible drivers of population divergence

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. 2011 ), Saccharomyces uvarum (Almeida et al. 2014 ), Coccidioides spp. (Neafsey et al. 2010 ) and Cryptococcus gatti (Engelthaler et al. 2014 ). In this issue of Molecular Ecology, Branco et al. ( 2015 ) apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.     

The number of markers and samples needed for detecting bottlenecks under realistic scenarios,with and without recovery: a simulation‐based study

Detecting bottlenecks is a common task in molecular ecology. While several bottleneck detection methods exist, evaluations of their power have focused only on severe bottlenecks (e.g. to Ne ~10). As a component of a recent review, Peery et al. ( 2012 ) analysed the power of two approaches, the M‐ratio and heterozygote excess tests, to detect moderate bottlenecks (e.g. to Ne ~100), which is realistic for many conservation situations. In this Comment, we address three important points relevant to but not considered in Peery et al. Under moderate bottleneck scenarios, we test the (i) relative advantage of sampling more markers vs. more individuals, (ii) potential power to detect the bottleneck when utilizing dozens of microsatellites (a realistic possibility for contemporary studies) and (iii) reduction in power when postbottleneck recovery has occurred. For the realistic situations examined, we show that (i) doubling the number of loci shows equal or better power than tripling the number of individuals, (ii) increasing the number of markers (up to 100) results in continued additive gains in power, and (iii) recovery after a moderate amount of time or gradual change in size reduces power, by up to one‐half. Our results provide a practical supplement to Peery et al. and encourage the continued use of bottleneck detection methods in the genomic age, but also emphasize that the power under different sampling schemes should be estimated, using simulation modelling, as a routine component of molecular ecology studies.     

Coalescing molecular evolution and DNA barcoding

The DNA barcoding concept (Woese et al. 1990 ; Hebert et al. 2003 ) has considerably boosted taxonomy research by facilitating the identification of specimens and discovery of new species. Used alone or in combination with DNA metabarcoding on environmental samples (Taberlet et al. 2012 ), the approach is becoming a standard for basic and applied research in ecology, evolution and conservation across taxa, communities and ecosystems (Scheffers et al. 2012 ; Kress et al. 2015 ). However, DNA barcoding suffers from several shortcomings that still remain overlooked, especially when it comes to species delineation (Collins & Cruickshank 2012 ). In this issue of Molecular Ecology, Barley & Thomson ( 2016 ) demonstrate that the choice of models of sequence evolution has substantial impacts on inferred genetic distances, with a propensity of the widely used Kimura 2‐parameter model to lead to underestimated species richness. While DNA barcoding has been and will continue to be a powerful tool for specimen identification and preliminary taxonomic sorting, this work calls for a systematic assessment of substitution models fit on barcoding data used for species delineation and reopens the debate on the limitation of this approach.     

Mock communities highlight the diversity of host‐associated eukaryotes

Host‐associated microbes are ubiquitous. Every multicellular eukaryote, and even many unicellular eukaryotes (protists), hosts a diverse community of microbes. High‐throughput sequencing (HTS) tools have illuminated the vast diversity of host‐associated microbes and shown that they have widespread influence on host biology, ecology and evolution (McFall‐Ngai et al. 2013 ). Bacteria receive most of the attention, but protists are also important components of microbial communities associated with humans (Parfrey et al. 2011 ) and other hosts. As HTS tools are increasingly used to study eukaryotes, the presence of numerous and diverse host‐associated eukaryotes is emerging as a common theme across ecosystems. Indeed, HTS studies demonstrate that host‐associated lineages account for between 2 and 12% of overall eukaryotic sequences detected in soil, marine and freshwater data sets, with much higher relative abundances observed in some samples (Ramirez et al. 2014 ; Simon et al. 2015 ; de Vargas et al. 2015 ). Previous studies in soil detected large numbers of predominantly parasitic lineages such as Apicomplexa, but did not delve into their origin [e.g. (Ramirez et al. 2014 )]. In this issue of Molecular Ecology, Geisen et al. ( 2015 ) use mock communities to show that many of the eukaryotic organisms detected by environmental sequencing in soils are potentially associated with animal hosts rather than free‐living. By isolating the host‐associated fraction of soil microbial communities, Geisen and colleagues help explain the surprisingly high diversity of parasitic eukaryotic lineages often detected in soil/terrestrial studies using high‐throughput sequencing (HTS) and reinforce the ubiquity of these host‐associated microbes. It is clear that we can no longer assume that organisms detected in bulk environmental sequencing are free‐living, but instead need to design studies that specifically enumerate the diversity and function of host‐associated eukaryotes. Doing so will allow the field to determine the role host‐associated eukaryotes play in soils and other environments and to evaluate hypotheses on assembly of host‐associated communities, disease ecology and more.     

Massive parallel MHC genotyping: titanium that shines

The power of population genetic analyses is often limited by sample size resulting from constraints in financial resources and time to genotype large numbers of individuals. This particularly applies to nonmodel species where detailed genomic knowledge is lacking. Next‐generation sequencing technology using primers ‘tagged’ with an individual barcode of a few nucleotides offers the opportunity to genotype hundreds of individuals at several loci in parallel ( Binladen et al. 2007 ; Meyer et al. 2008 ). The large number of sequence reads can also be used to identify artefacts by frequency distribution thresholds intrinsically determined for each run and data set. In Babik et al. (2009 ), next‐generation deep sequencing was used to genotype several major histocompatibility complex (MHC) class IIB loci of the European bank vole ( Fig. 1 ). Their approach can be useful for many researchers working with complex multiallelic templates and large sample sizes.

Figure 1 Open in figure viewer PowerPoint Hypothetical example of parallel genotyping of two individuals using individually bar‐coded primers. Polymerase chain reactions (PCRs) are performed separately for each individual using a forward primer with a unique Tag‐sequence of four nucleotides. After sequencing of pooled PCR products, sequences can be sorted by their forward primer Tag (Tag‐sorting error rate was estimated < 0.1%). Rare sequences most likely represent artefacts and due to the large amount of sequences obtained (up to 10⁶) the artefact threshold can be determined intrinsically for each data set and was estimated to be around 3% in the case of bank vole MHC class IIB genes ( Babik et al. 2009 ). Photos by Gabriela Bydlon.     

Bacteriophage exclusion,a new defense system

The ability to withstand viral predation is critical for survival of most microbes. Accordingly, a plethora of phage resistance systems has been identified in bacterial genomes (Labrie et al, 2010 ), including restriction‐modification systems (R‐M) (Tock & Dryden, 2005 ), abortive infection (Abi) (Chopin et al, 2005 ), Argonaute‐based interference (Swarts et al, 2014 ), as well as clustered regularly interspaced short palindromic repeats (CRISPR) and associated protein (Cas) adaptive immune system (CRISPR‐Cas) (Barrangou & Marraffini, 2014 ; Van der Oost et al, 2014 ). Predictably, the dark matter of bacterial genomes contains a wealth of genetic gold. A study published in this issue of The EMBO Journal by Goldfarb et al ( 2015 ) unveils bacteriophage exclusion (BREX) as a novel, widespread bacteriophage resistance system that provides innate immunity against virulent and temperate phage in bacteria.     

It's a family act: the geminin triplets take center stage in motile ciliogenesis

The balance between proliferation and differentiation is a fundamental aspect of multicellular life. Perhaps nowhere is this delicate balance more palpable than in the multiciliated cells (MCCs) that line the respiratory tract, the ependyma, and the oviduct. These cells contain dozens to hundreds of motile cilia that beat in a concerted fashion to generate directed fluid flow over the tissue surface. Although MCCs have exited the cell cycle, remarkably, they retain the ability to duplicate their centrioles and to mature those centrioles into ciliary basal bodies—two features, which are known to be normally under strict cell cycle control (Firat‐Karalar & Stearns, 2014 ). How post‐mitotic MCCs retain this ability, remains unclear. In the past several months, four research articles, including one from Terré et al in this issue of The EMBO Journal, have described a vital role for the geminin coiled‐coil domain‐containing protein (Gemc1) in the MCC gene expression program in multiple tissues and organisms, that bring us closer to understanding this question (Kyrousi et al, 2015 ; Zhou et al, 2015 ; Arbi et al, 2016 ; Terré et al, 2016 ).     

Genomics and museum specimens

Nearly 25 years ago, Allan Wilson and colleagues isolated DNA sequences from museum specimens of kangaroo rats (Dipodomys panamintinus) and compared these sequences with those from freshly collected animals (Thomas et al. 1990 ). The museum specimens had been collected up to 78 years earlier, so the two samples provided a direct temporal comparison of patterns of genetic variation. This was not the first time DNA sequences had been isolated from preserved material, but it was the first time it had been carried out with a population sample. Population geneticists often try to make inferences about the influence of historical processes such as selection, drift, mutation and migration on patterns of genetic variation in the present. The work of Wilson and colleagues was important in part because it suggested a way in which population geneticists could actually study genetic change in natural populations through time, much the same way that experimentalists can do with artificial populations in the laboratory. Indeed, the work of Thomas et al. ( 1990 ) spawned dozens of studies in which museum specimens were used to compare historical and present‐day genetic diversity (reviewed in Wandeler et al. 2007 ). All of these studies, however, were limited by the same fundamental problem: old DNA is degraded into short fragments. As a consequence, these studies mostly involved PCR amplification of short templates, usually short stretches of mitochondrial DNA or microsatellites. In this issue, Bi et al. ( 2013 ) report a breakthrough that should open the door to studies of genomic variation in museum specimens. They used target enrichment (exon capture) and next‐generation (Illumina) sequencing to compare patterns of genetic variation in historic and present‐day population samples of alpine chipmunks (Tamias alpinus) (Fig. 1). The historic samples came from specimens collected in 1915, so the temporal span of this comparison is nearly 100 years.     

Aquatic biodiversity assessment for the lazy

The world is covered in DNA. In any ecosystem, extracellular DNA fragments can be found that once formed the genomes of a variety of micro‐ and macroorganisms. A few years ago, it was proposed to use this environmental DNA (eDNA) as a source of information on local vertebrate biodiversity (Ficetola et al. 2008 ; Taberlet et al. 2012 ). This idea offered an elegant solution to take up the gauntlet of rapidly increasing monitoring needs. Coupled with barcoding efforts, it promised to be cost‐efficient in many respects, for example man‐hours and taxonomic expertise. Ecologists and conservation biologists with an interest in aquatic ecosystems have enthusiastically adopted and pioneered this new method, producing dozens of eDNA studies. Most of these studies have, however, focused on a single or a few aquatic species. In this issue of Molecular Ecology, Valentini et al. ( 2016 ) move the field a step further by demonstrating that metabarcoding approaches – which simultaneously target large groups of organisms such as amphibians or fish – can match and sometimes even outperform other inventory methods.     

Ecological genomics in full colour

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning et al. ( 2014 ) detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson et al. ( 2014 ) identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson et al. ( 2014 ) go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics – efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.     

Ready,Set…Poised!: Polycomb target genes are bound by poised RNA polymerase II throughout differentiation

In embryonic stem cells (ESCs), silent genes with major developmental functions display a unique epigenetic state in which strong and broad binding by Polycomb repressive complexes (PRCs) is accompanied by the presence of poised RNA polymerase II (RNAPII) and activating histone marks (e.g. H3K4me3) (Azuara et al, 2006 ; Bernstein et al, 2006 ; Stock et al, 2007 ; Brookes et al, 2012 ). It has been suggested that the plasticity and broad differentiation potential of pluripotent cells might rely, at least partly, on this unique epigenetic state (Bernstein et al, 2006 ; Stock et al, 2007 ). In their recent study, Pombo and colleagues (Ferrai et al, 2017 ) show that a similar epigenetic state can be found at a subset of major developmental genes throughout the differentiation of ESCs into neurons, providing novel and exciting insights into the molecular basis of cellular plasticity in differentiated cells.     

Trade‐offs and utility of alternative RADseq methods: Reply to Puritz et al.

Puritz et al. provide a review of several RADseq methodological approaches in response to our ‘Population Genomic Data Analysis’ workshop (Sept 2013) review (Andrews & Luikart 2014). We agree with Puritz et al. on the importance for researchers to thoroughly understand RADseq library preparation and data analysis when choosing an approach for answering their research questions. Some of us are currently using multiple RADseq protocols, and we agree that the different methods may offer advantages in different cases. Our workshop review did not intend to provide a thorough review of RADseq because the workshop covered a broad range of topics within the field of population genomics. Similarly, neither the response of Puritz et al. nor our comments here provide sufficient space to thoroughly review RADseq. Nonetheless, here we address some key points that we find unclear or potentially misleading in their evaluation of techniques.     

A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour‐intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra 2008 ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.