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

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

EDENetworks: A user‐friendly software to build and analyse networks in biogeography,ecology and population genetics

The recent application of graph‐based network theory analysis to biogeography, community ecology and population genetics has created a need for user‐friendly software, which would allow a wider accessibility to and adaptation of these methods. EDENetworks aims to fill this void by providing an easy‐to‐use interface for the whole analysis pipeline of ecological and evolutionary networks starting from matrices of species distributions, genotypes, bacterial OTUs or populations characterized genetically. The user can choose between several different ecological distance metrics, such as Bray‐Curtis or Sorensen distance, or population genetic metrics such as F_ST or Goldstein distances, to turn the raw data into a distance/dissimilarity matrix. This matrix is then transformed into a network by manual or automatic thresholding based on percolation theory or by building the minimum spanning tree. The networks can be visualized along with auxiliary data and analysed with various metrics such as degree, clustering coefficient, assortativity and betweenness centrality. The statistical significance of the results can be estimated either by resampling the original biological data or by null models based on permutations of the data.     

Evolution and biodiversity of Antarctic organisms: a molecular perspective

The Antarctic biota is highly endemic, and the diversity and abundance of taxonomic groups differ from elsewhere in the world. Such characteristics have resulted from evolution in isolation in an increasingly extreme environment over the last 100 Myr. Studies on Antarctic species represent some of the best examples of natural selection at the molecular, structural and physiological levels. Analyses of molecular genetics data are consistent with the diversity and distribution of marine and terrestrial taxa having been strongly influenced by geological and climatic cooling events over the last 70 Myr. Such events have resulted in vicariance driven by continental drift and thermal isolation of the Antarctic, and in pulses of species range contraction into refugia and subsequent expansion and secondary contact of genetically distinct populations or sister species during cycles of glaciation. Limited habitat availability has played a major role in structuring populations of species both in the past and in the present day. For these reasons, despite the apparent simplicity or homogeneity of Antarctic terrestrial and marine environments, populations of species are often geographically structured into genetically distinct lineages. In some cases, genetic studies have revealed that species defined by morphological characters are complexes of cryptic or sibling species. Climate change will cause changes in the distribution of many Antarctic and sub-Antarctic species through affecting population-level processes such as life history and dispersal.     

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae)

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.     

Hybridization effects and genetic diversity of the common and black‐tufted marmoset (Callithrix jacchus and Callithrix penicillata) mitochondrial control region

Hybridization is continually documented in primates, but effects of natural and anthropogenic hybridization on biodiversity are still unclear and differentiating between these contexts remains challenging in regards to primate evolution and conservation. Here, we examine hybridization effects on the mitochondrial DNA (mtDNA) control region of Callithrix marmosets, which provide a unique glimpse into interspecific mating under distinct anthropogenic and natural conditions. DNA was sampled from 40 marmosets along a 50‐km transect from a previously uncharacterized hybrid zone in NE Brazil between the ranges of Callithrix jacchus and Callithrix penicillata. DNA was also collected from 46 marmosets along a 30‐km transect in a hybrid zone in Rio de Janeiro state, Brazil, where exotic marmosets appeared in the 1980s. Combining Callithrix DNA sampled inside and outside of these hybrid zones, phylogenetic and network analyses show C. jacchus and C. penicillata being parental species to sampled hybrids. We expand limited Callithrix population genetics work by describing mtDNA diversity and demographic history of these parental species. We show ancient population expansion in C. jacchus and historically constant population size in C. penicillata, with the latter being more genetically diverse than the former. The natural hybrid zone contained higher genetic diversity relative to the anthropogenic zone. While our data suggest hybrid swarm formation within the anthropogenic zone due to removed physical reproductive barriers, this pattern is not seen in the natural hybrid zone. These results suggest different genetic dynamics within natural and anthropogenic hybridization contexts that carry important implications for primate evolution and conservation. Am J Phys Anthropol 155:522–536, 2014. © 2014 Wiley Periodicals, Inc.     

Into the wild—WAMBAM goes to Canada

The sixth Wild Animal Models Bi‐Annual Meeting was held in July 2017 in Québec, with 42 participants. This report documents the evolution of questions asked and approaches used in evolutionary quantitative genetic studies of wild populations in recent decades, and how these questions and approaches were represented at the recent meeting. We explore how ideas from previous meetings in this series have developed to their present states, and consider how the format of the meetings may be particularly useful at fostering the rapid development and proliferation of ideas and approaches.     

A (not-so-radical) solution to the species problem

What are species? One popular answer is that species are individuals. Here I develop another approach to thinking about species, an approach based on the notion of a lineage. A lineage is a sequence of reproducing entities, individuated in terms of its components. I argue that one can conceive of species as groups of lineages, either organism lineages or population lineages. Conceiving of species as groups of lineages resolves the problems that the individual conception of species is supposed to resolve. It has added the virtue of focusing attention on the characteristic of species that is most relevant to understanding their role in evolutionary processes, namely, the lineage structure of species.     

Population Genetics and Molecular Evolution of DNA Sequences in Transposable Elements. I. A Simulation Framework

A population genetic simulation framework is developed to understand the behavior and molecular evolution of DNA sequences of transposable elements. Our model incorporates random transposition and excision of transposable element (TE) copies, two modes of selection against TEs, and degeneration of transpositional activity by point mutations. We first investigated the relationships between the behavior of the copy number of TEs and these parameters. Our results show that when selection is weak, the genome can maintain a relatively large number of TEs, but most of them are less active. In contrast, with strong selection, the genome can maintain only a limited number of TEs but the proportion of active copies is large. In such a case, there could be substantial fluctuations of the copy number over generations. We also explored how DNA sequences of TEs evolve through the simulations. In general, active copies form clusters around the original sequence, while less active copies have long branches specific to themselves, exhibiting a star-shaped phylogeny. It is demonstrated that the phylogeny of TE sequences could be informative to understand the dynamics of TE evolution.     

Biogeography of the grasses (Poaceae): a phylogenetic approach to reveal evolutionary history in geographical space and geological time

The grasses (Poaceae) are the fifth most diverse family of angiosperms, including 800 genera and more than 10 000 species. Few phylogenetic studies have tried to investigate palaeo‐biogeographical and palaeo‐ecological scenarios that may have led to present‐day distribution and diversity of grasses at the family level. We produced a dated phylogenetic tree based on combined plastid DNA sequences and a comprehensive sample of Poaceae. Furthermore, we produced an additional tree using a supermatrix of morphological and molecular data that included all 800 grass genera so that ancestral biogeography and ecological habitats could be inferred. We used a likelihood‐based method, which allows the estimation of ancestral polymorphism in both biogeographical and ecological analyses for large data sets. The origin of Poaceae was retrieved as African and shade adapted. The crown node of the BEP + PACCMAD clade was dated at 57 Mya, in the early Eocene. Grasses dispersed to all continents by approximately 60 million years after their Gondwanan origin in the late Cretaceous. PACCMAD taxa adapted to open habitats as early as the late Eocene, a date consistent with recent phytolith fossil data for North America. C₄ photosynthesis first originated in Africa, at least for Chloridoideae in the Eocene at c. 30 Mya. The BEP clade members adapted to open habitats later than PACCMAD members; this was inferred to occur in Eurasia in the Oligocene. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 162 , 543–557.     

Stereoselectivity of Mucorales lipases toward triradylglycerols--a simple solution to a complex problem

The lipases from Rhizopus and Rhizomucor are members of the family of Mucorales lipases. Although they display high sequence homology, their stereoselectivity toward triradylglycerols (sn-2 substituted triacylglycerols) varies. Four different triradylglycerols were investigated, which were classified into two groups: flexible substrates with rotatable O'-C1' ether or ester bonds adjacent to C2 of glycerol and rigid substrates with a rigid N'-C1' amide bond or a phenyl ring in sn-2. Although Rhizopus lipase shows opposite stereopreference for flexible and rigid substrates (hydrolysis in sn-1 and sn-3, respectively), Rhizomucor lipase hydrolyzes both groups of triradylglycerols preferably in sn-1. To explain these experimental observations, computer-aided molecular modeling was applied to study the molecular basis of stereoselectivity. A generalized model for both lipases of the Mucorales family highlights the residues mediating stereoselectivity: (1) L258, the C-terminal neighbor of the catalytic histidine, and (2) G266, which is located in a loop contacting the glycerol backbone of a bound substrate. Interactions with triradylglycerol substrates are dominated by van der Waals contacts. Stereoselectivity can be predicted by analyzing the value of a single substrate torsion angle that discriminates between sn-1 and sn-3 stereopreference for all substrates and lipases investigated here. This simple model can be easily applied in enzyme and substrate engineering to predict Mucorales lipase variants and synthetic substrates with desired stereoselectivity.     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.     

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill‐coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function‐altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal‐level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.     

THE INTENSITY OF SELECTION ACTING ON THE COUCH POTATO GENE—SPATIAL–TEMPORAL VARIATION IN A DIAPAUSE CLINE

Cosmopolitan populations of Drosophila melanogaster have co‐opted a form of reproductive diapause to overwinter in northern populations. Polymorphism in the couch potato gene has been implicated in genetic variation for this diapause trait. Using a collection of 20 populations from Florida to Canada and 11 collections from 3 years in a Pennsylvania orchard, we estimated the allele frequencies for 15 single nucleotide polymorphisms (SNPs) in the couch potato gene. These include the specific polymorphism associated with diapause inducability. We find that the SNP polymorphism, 48034(A/T), is correlated with latitude and its frequencies are predicted by the incidence of diapause trait. We find that the clinal patterns for cpo SNPs sampled in 1997 are similar to the same SNPs sampled in 2009–2010. SNPs that show apparent associations with cpo expression are also clinal with the low‐expression allele increasing in frequency, as would be predicted from functional knockout studies of cpo. Finally, we see a significant pattern where the frequency of the diapause‐causing allele drops in frequency during the summer season, consistent with the drop in the incidence of the diapause trait. The selection required to drive this response is large, roughly 24% to 59% per generation depending on the degree of dominance.     

Land plants and DNA barcodes: short-term and long-term goals

Land plants have had the reputation of being problematic for DNA barcoding for two general reasons: (i) the standard DNA regions used in algae, animals and fungi have exceedingly low levels of variability and (ii) the typically used land plant plastid phylogenetic markers (e.g. rbcL, trnL-F, etc.) appear to have too little variation. However, no one has assessed how well current phylogenetic resources might work in the context of identification (versus phylogeny reconstruction). In this paper, we make such an assessment, particularly with two of the markers commonly sequenced in land plant phylogenetic studies, plastid rbcL and internal transcribed spacers of the large subunits of nuclear ribosomal DNA (ITS), and find that both of these DNA regions perform well even though the data currently available in GenBank/EBI were not produced to be used as barcodes and BLAST searches are not an ideal tool for this purpose. These results bode well for the use of even more variable regions of plastid DNA (such as, for example, psbA-trnH) as barcodes, once they have been widely sequenced. In the short term, efforts to bring land plant barcoding up to the standards being used now in other organisms should make swift progress. There are two categories of DNA barcode users, scientists in fields other than taxonomy and taxonomists. For the former, the use of mitochondrial and plastid DNA, the two most easily assessed genomes, is at least in the short term a useful tool that permits them to get on with their studies, which depend on knowing roughly which species or species groups they are dealing with, but these same DNA regions have important drawbacks for use in taxonomic studies (i.e. studies designed to elucidate species limits). For these purposes, DNA markers from uniparentally (usually maternally) inherited genomes can only provide half of the story required to improve taxonomic standards being used in DNA barcoding. In the long term, we will need to develop more sophisticated barcoding tools, which would be multiple, low-copy nuclear markers with sufficient genetic variability and PCR-reliability; these would permit the detection of hybrids and permit researchers to identify the 'genetic gaps' that are useful in assessing species limits.     

These aren’t the loci you’e looking for: Principles of effective SNP filtering for molecular ecologists

Sequencing reduced‐representation libraries of restriction site‐associated DNA (RADseq) to identify single nucleotide polymorphisms (SNPs) is quickly becoming a standard methodology for molecular ecologists. Because of the scale of RADseq data sets, putative loci cannot be assessed individually, making the process of filtering noise and correctly identifying biologically meaningful signal more difficult. Artefacts introduced during library preparation and/or bioinformatic processing of SNP data can create patterns that are incorrectly interpreted as indicative of population structure or natural selection. Therefore, it is crucial to carefully consider types of errors that may be introduced during laboratory work and data processing, and how to minimize, detect and remove these errors. Here, we discuss issues inherent to RADseq methodologies that can result in artefacts during library preparation and locus reconstruction resulting in erroneous SNP calls and, ultimately, genotyping error. Further, we describe steps that can be implemented to create a rigorously filtered data set consisting of markers accurately representing independent loci and compare the effect of different combinations of filters on four RAD data sets. At last, we stress the importance of publishing raw sequence data along with final filtered data sets in addition to detailed documentation of filtering steps and quality control measures.     

Speciation and mitochondrial DNA diversification of the diving beetles Agabus bipustulatus and A. wollastoni (Coleoptera, Dytiscidae) within Macaronesia

Diversification of populations of two morphologically similar diving beetles within the Agabus tristis group, A. wollastoni and A. bipustulatus, was investigated, with partial mtDNA cytochrome b (Cyt b) sequences, allozymes and landmark‐based morphometrics. The Madeiran endemic A. wollastoni was collected from 11 localities. Population genetic and morphological variation was compared to Scandinavian localities of the widespread west Palearctic A. bipustulatus, recorded also from the Azores. Agabus wollastoni and European A. bipustulatus specimens representing eight and 13 localities respectively, were used in evaluating their phylogenetic relationship. Maximum parsimony analysis of the Cyt b sequences showed that both the A. bipustulatus and A. wollastoni specimens form well‐supported monophyletic groups. Three lines of evidence suggest that Agabus wollastoni has speciated through a few founders: (1) a well‐supported mtDNA line; (2) the mean heterozygosity of A. wollastoni is lower when compared to A. bipustulatus on the mainland; and (3) several uncommon alleles of A. bipustulatus are missing in A. wollastoni. The Azorean A. bipustulatus population was drastically affected by the colonization event, since several loci have become fixed with a resulting lower mean heterozygosity. The colonization was relatively recent, as the mtDNA lineage represented in the Azores is deeply nested within the A. bipustulatus clade. Population structure shows moderate inbreeding of A. wollastoni, and extensive substructuring at all localities with moderate gene‐flow between them. Morphological variation in A. wollastoni showed significant differentiation among several populations. Island colonizations, population structure of A. wollastoni, and an observed pattern of variation of the α‐glycerophosphate dehydrogenase locus are discussed. © 2003 The Linnean Society of London, Biological Journal of the Linnean Society, 2003, 79, 653–666.     

Mitochondrial DNA phylogenetics of the Galaxias vulgaris complex from South Island, New Zealand: rapid radiation of a species flock

Mitochondrial control region sequence variation was examined in the Galaxias vulgaris complex, a group of freshwater-limited galaxiid fishes endemic to South Island, New Zealand. Phylogenetic analyses were used to test the monophyly of seven non-migratory ( G. vulgaris complex) and one migratory ( G. brevipinnis ) species. Newly-described taxa, some diagnosed on the basis of subtle differences in morphology, are associated with strongly monophyletic mtDNA lineages. Although the current taxonomy is supported generally, sequence data suggest that Southland G. anomalus should be attributed to G. gollumoides , whereas Southland and Stewart Island G. depressiceps represent a new southern lineage. Molecular clock calibrations suggest that migratory and non-migratory forms, separated by a maximum of 6·4% sequence divergence, diverged no earlier than the mid-late Pliocene. Biogeographical implications of high diversity (eight lineages) in the south (Otago-Southland) and low diversity (one lineage) in central South Island (Canterbury) are discussed. The unresolved relationships among species may reflect a rapid evolutionary radiation, essentially a star phylogeny.     

The First Internal Transcribed Spacer (ITS-1) of Ribosomal DNA as a Molecular Marker for Phylogenetic and Population Analyses in Crustacea

The objective of the present study is to explore the feasibility of using the first internal transcribed spacer (ITS-1) of ribosomal DNA as a molecular marker for studying the interspecific and intraspecific genetic variations among crustaceans. We designed primers that could amplify ITS-1 from a majority of taxonomic groups of crustaceans. The gene was found to exhibit a high degree of length polymorphism among different groups, ranging from 182 bp in the barnacle Balanus amphitrite to approximately 820 bp in the spiny lobster Panulirus japonicus. With respect to differences between congeneric species, it was found that the ITS-1 sequences of 3 mitten crabs, Eriocheir sinensis, Eriocheir leptognathus, and Eriocheir formosa, exhibit 5.4% to 16.3% nucleotide divergence, suggesting that ITS-1 is informative for phylogenetic analysis at the species level. Yet there are extensive (0.9%–2.3%) variations within individual E. formosa, so that phylogenetic analyses could be obscured. ITS-1 was found to vary between 2 geographical populations of the shrimp Penaeus japonicus. The variations involved substitutions as well as insertions/deletions between shrimp from Australia and South China Sea. These results show that ITS-1 is highly divergent among different crustaceans and could be an appropriate marker for molecular systematic studies at the species and population levels, although the presence of intragenomic variation needs to be taken into consideration. Received August 15, 2000; accepted February 9, 2001