Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Parallel evolution of facial stripe patterns in the Neolamprologus brichardilpulcher species complex endemic to Lake Tanganyika

Colour pattern diversity can be due to random processes or to natural or sexual selection. Consequently, similarities in colour patterns are not always correlated with common ancestry, but may result from convergent evolution under shared selection pressures or drift. Neolamprologus brichardi and Neolamprologus pulcher have been described as two distinct species based on differences in the arrangement of two dark bars on the operculum. Our study uses DNA sequences of the mitochondrial control region to show that relatedness of haplotypes disagrees with species assignment based on head colour pattern. This suggests repeated parallel evolution of particular stripe patterns. The complete lack of shared haplotypes between populations of the same or different phenotypes reflects strong philopatric behaviour, possibly induced by the cooperative breeding mode in which offspring remain in their natal territory and serve as helpers until they disperse to nearby territories or take over a breeding position. Concordant phylogeographic patterns between N. brichardi/N. pulcher populations and other rock-dwelling cichlids suggest that the same colonization routes have been taken by sympatric species and that these routes were affected by lake level fluctuations in the past. (C) 2007 Elsevier Inc. All rights reserved.

Single-Molecule Conductance of Functionalized Oligoynes: Length Dependence and Junction Evolution

We report a combined experimental and theoretical investigation of the length dependence and anchor group dependence of the electrical conductance of a series of oligoyne molecular wires in single-molecule junctions with gold contacts. Experimentally, we focus on the synthesis and properties of diaryloligoynes with n = 1, 2, and 4 triple bonds and the anchor dihydrobenzo[b]thiophene (BT). For comparison, we also explored the aurophilic anchor group cyano (CN), amino (NH2), thiol (SH), and 4-pyridyl (PY). Scanning tunneling microscopy break junction (STM-BJ) and mechanically controllable break junction (MCBJ) techniques are employed to investigate single-molecule conductance characteristics. The BT moiety is superior as compared to traditional anchoring groups investigated so far. BT-terminated oligoynes display a 100% probability of junction formation and possess conductance values which are the highest of the oligoynes studied and, moreover, are higher than other conjugated molecular wires of similar length. Density functional theory (DFT)-based calculations are reported for oligoynes with n = 1−4 triple bonds. Complete conductance traces and conductance distributions are computed for each family of molecules. The sliding of the anchor groups leads to oscillations in both the electrical conductance and the binding energies of the studied molecular wires. In agreement with experimental results, BT-terminated oligoynes are predicted to have a high electrical conductance. The experimental attenuation constants βH range between 1.7 nm−1 (CN) and 3.2 nm−1 (SH) and show the following trend: βH(CN) < βH(NH2) < βH(BT) < βH(PY) ≈ βH(SH). DFT-based calculations yield lower values, which range between 0.4 nm−1 (CN) and 2.2 nm−1 (PY).

Molecular epidemiology of Flavobacterium psychrophilum from Swiss fish farms

Bacterial cold water disease (BCWD) and rainbow trout fry syndrome (RTFS) caused by Flavobacterium psychrophilum are 2 of the major diseases causing high fish mortality in salmonid fish farms. The molecular epidemiology of F. psychrophilum is still largely unknown. Multilocus sequence typing (MLST) has been previously used for this pathogen and underscored a correlation between clonal complexes and host fish species. Here we used MLST to study the relationships among 112 F. psychrophilum isolates from rainbow trout Oncorhynchus mykiss and brown trout Salmo trutta fario and S. t. lacustris in Swiss fish farms between 1993 and 2012. The isolates belonged to 27 different sequence types (STs). Most of the Swiss outbreaks were associated with strains belonging to clonal complexes CC-ST2 and CC-ST90, found in both rainbow trout and brown trout and represented by several STs. Eight ST singletons could not be connected to any known clonal complex. Already reported from other parts of Europe and North America, CC-ST2 was the most frequent clonal complex observed, and it caused the majority of outbreaks in Switzerland, with CC-ST90 being the second most important type. In the tightly interconnected Swiss fish farms, no association between clonal complex and host fish was detected, but a temporal evolution of the frequency of some STs was observed. The occurrence of sporadic STs suggests high F. psychrophilum diversity and may reflect the presence of different sequence types in the environment.

Adaptive evolution during an ongoing range expansion: the invasive bank vole (Myodes glareolus) in Ireland

Range expansions are extremely common, but have only recently begun to attract attention in terms of their genetic consequences. As populations expand, demes at the wave front experience strong genetic drift, which is expected to reduce genetic diversity and potentially cause ‘allele surfing’, where alleles may become fixed over a wide geographical area even if their effects are deleterious. Previous simulation models show that range expansions can generate very strong selective gradients on dispersal, reproduction, competition and immunity. To investigate the effects of range expansion on genetic diversity and adaptation, we studied the population genomics of the bank vole (Myodes glareolus) in Ireland. The bank vole was likely introduced in the late 1920s and is expanding its range at a rate of ~2.5 km/year. Using genotyping-by-sequencing, we genotyped 281 bank voles at 5979 SNP loci. Fourteen sample sites were arranged in three transects running from the introduction site to the wave front of the expansion. We found significant declines in genetic diversity along all three transects. However, there was no evidence that sites at the wave front had accumulated more deleterious mutations. We looked for outlier loci with strong correlations between allele frequency and distance from the introduction site, where the direction of correlation was the same in all three transects. Amongst these outliers, we found significant enrichment for genic SNPs, suggesting the action of selection. Candidates for selection included several genes with immunological functions and several genes that could influence behaviour.

The evolution of host associations in the parasitic wasp genus Ichneumon (Hymenoptera: Ichneumonidae): convergent adaptations to host pupation sites

Background: The diversification of organisms with a parasitic lifestyle is often tightly linked to the evolution of their host associations. If a tight host association exists, closely related species tend to attack closely related hosts; host associations are less stable if associations are determined by more plastic traits like parasitoid searching and oviposition behaviour. The pupal-parasitoids of the genus Ichneumon attack a variety of macrolepidopteran hosts.They are either monophagous or polyphagous, and therefore offer a promissing system to investigate the evolution of host associations. Ichneumon was previously divided into two groups based on general body shape; however, a stout shape has been suggested as an adaptation to buried host pupation sites, and might thus not represent a reliable phylogenetic character. Results: We here reconstruct the first molecular phylogeny of the genus Ichneumon using two mitochondrial (CO1 and NADH1) and one nuclear marker (28S). The resulting phylogeny only supports monophyly of Ichneumon when Ichneumon lugens Gravenhorst, 1829 (formerly in Chasmias, stat. rev.) and Ichneumon deliratorius Linnaeus, 1758 (formerly Coelichneumon) are included. Neither parasitoid species that attack hosts belonging to one family nor those attacking butterflies (Rhopalocera) form monophyletic clades. Ancestral state reconstructions suggest multiple transitions between searching for hosts above versus below ground and between a stout versus elongated body shape. A model assuming correlated evolution between the two characters was preferred over independent evolution of host-searching niche and body shape. Conclusions: Host relations, both in terms of phylogeny and ecology, evolved at a high pace in the genus Ichneumon. Numerous switches between hosts of different lepidopteran families have occurred, a pattern that seems to be the rule among idiobiont parasitoids. A stout body and antennal shape in the parasitoid female is confirmed as an ecological adaptation to host pupation sites below ground and has evolved convergently several times. Morphological characters that might be involved in adaptation to hosts should be avoided as diagnostic characters for phylogeny and classification, as they can be expected to show high levels of homoplasy.

BCL2 mutation spectrum in B-cell non-Hodgkin lymphomas and patterns associated with evolution of follicular lymphoma

BCL2 is a target of somatic hypermutation in t(14;18) positive and also in a small fraction of t(14;18) negative diffuse large B-cell lymphoma (DLBCL), suggesting an aberrant role of somatic hypermutation (ASHM). To elucidate the prevalence of BCL2 mutations in lymphomas other than DLBCL, we Sanger-sequenced the hypermutable region of the BCL2 gene in a panel of 69 mature B-cell lymphomas, including Richter's syndrome DLBCL, marginal-zone lymphomas, post-transplant lymphoproliferative disorders, HIV-associated and common-variable immunodeficiency-associated DLBCL, all known to harbour ASHM-dependent mutations in other genes, as well as 16 t(14,18) negative and 21 t(14;18) positive follicular lymphomas (FLs). We also investigated the pattern of BCL2 mutations in longitudinal samples from 10 FL patients relapsing to FL or transforming to DLBCL (tFL). By direct sequencing, we found clonally represented BCL2 mutations in 2/16 (13%) of t(14;18) negative FLs, 2/16 (13%) HIV-DLBCLs, 1/9 (11%) of Richter's syndrome DLBCL, 1/17 (6%) of post-transplant lymphoproliferative disorders and 1/2 (50%) common-variable immunodeficiency-associated DLBCL. The proportion of mutated cases was significantly lower than in FLs carrying the t(14;18) translocation (15/21, 71%). However, the absence of t(14;18) by FISH or PCR and the molecular features of the mutations strongly suggest that BCL2 represents an additional target of ASHM in these entities. Analysis of the BCL2 mutation pattern in clonally related FL/FL and FL/tFL samples revealed two distinct scenarios of genomic evolution: (i) direct evolution from the antecedent FL clone, with few novel clonal mutations acquired by the tFL major clone, and (ii) evolution from a common mutated long-lived progenitor cell, which subsequently acquired distinct mutations in the FL and in the relapsed or transformed counterpart. Copyright © 2014 John Wiley & Sons, Ltd.

Genome profiling of sterol synthesis shows convergent evolution in parasites and guides chemotherapeutic attack

Sterols are an essential class of lipids in eukaryotes, where they serve as structural components of membranes and play important roles as signaling molecules. Sterols are also of high pharmacological significance: cholesterol-lowering drugs are blockbusters in human health, and inhibitors of ergosterol biosynthesis are widely used as antifungals. Inhibitors of ergosterol synthesis are also being developed for Chagas's disease, caused by Trypanosoma cruzi. Here we develop an in silico pipeline to globally evaluate sterol metabolism and perform comparative genomics. We generate a library of hidden Markov model-based profiles for 42 sterol biosynthetic enzymes, which allows expressing the genomic makeup of a given species as a numerical vector. Hierarchical clustering of these vectors functionally groups eukaryote proteomes and reveals convergent evolution, in particular metabolic reduction in obligate endoparasites. We experimentally explore sterol metabolism by testing a set of sterol biosynthesis inhibitors against trypanosomatids, Plasmodium falciparum, Giardia, and mammalian cells, and by quantifying the expression levels of sterol biosynthetic genes during the different life stages of T. cruzi and Trypanosoma brucei. The phenotypic data correlate with genomic makeup for simvastatin, which showed activity against trypanosomatids. Other findings, such as the activity of terbinafine against Giardia, are not in agreement with the genotypic profile.

Evolution of colour patterns in East African cichlid fish

African cichlid fishes have undergone outbursts of explosive speciation in several lakes, accompanied by rapid radiations in coloration and ecology. Little is known about the evolutionary forces that triggered these events but a hypothesis, published by Wallace Dominey in 1984, has figured prominently. It states that the evolution of colour patterns is driven by sexual selection and that these colour patterns are important in interspecific mate choice, a combination which holds the potential for rapid speciation. Here we present phylogenetic analyses that describe major events in colour evolution and test predictions yielded by Dominey's hypothesis. We assembled information on stripe patterns and the presence or absence of nuptial coloration from more than 700 cichlid species representing more than 90 taxa for which molecular phylogenetic hypotheses were available. We show that sexual selection is most likely the selection force that made male nuptial coloration arise and evolve quickly. In contrast, stripe patterns, though phylogenetically not conserved either, are constrained ecologically. The evolution of vertical bar patterns is associated with structurally complex habitats, such as rocky substrates or vegetation. The evolution of a horizontal stripe is associated with a piscivorous feeding mode. Horizontal stripes are also associated with shoaling behaviour. Strength of sexual selection, measured in terms of the mating system (weak in monogamous, strong in promiscuous species), has no detectable effects on stripe pattern evolution. In promiscuous species the frequency of difference between sister species in nuptial hue is higher than in pair bonding and harem forming species, but the frequency of difference in stripe pattern is lower. We argue that differences between the two components of coloration in their exposure to natural selection explain their very different evolutionary behaviour. Finally, we suggest that habitat-mediated selection upon chromomotor flexibility, a special form of phenotypic plasticity found in the river-dwelling outgroups of the lake-dwelling cichlids, explains the rapid and recurrent ecology-associated radiation of stripe patterns in lake environments, a new hypothesis that yields experimentally testable predictions.

In vivo Evolution of CMY-2 to CMY-33 β-Lactamase in Escherichia coli ST131: Characterization of an Acquired Extended-Spectrum AmpC (ESAC) Conferring Resistance to Cefepime.

Cefepime is frequently prescribed to treat infections caused by AmpC-producing Gram-negative bacteria. CMY-2 is the most common plasmid-mediated AmpC (pAmpC) β-lactamase. Unfortunately, CMY variants conferring enhanced cefepime resistance are reported. Here, we describe the evolution of CMY-2 to an extended-spectrum AmpC (ESAC) in clonally identical E. coli isolates obtained from a patient. The CMY-2-producing E. coli (CMY-2-Ec) was isolated from a wound. Thirty days later, one CMY-33-producing E. coli (CMY-33-Ec) was detected in bronchoalveolar lavage. Two weeks before the isolation of CMY-33-Ec, the patient received cefepime.CMY-33-Ec and CMY-2-Ec were identical by rep-PCR, being of hyperepidemic ST131, but showed different β-lactam MICs (e.g., cefepime 16 vs. ≤0.5 μg/ml). Identical CMY-2-Ec isolates were also found in a rectal swab. CMY-33 differs from CMY-2 by a Leu293-Ala294 deletion. Expressed in E. coli DH10B, both CMYs conferred resistance to ceftazidime (≥256 μg/ml), but cefepime MICs were higher for CMY-33 than CMY-2 (8 vs. 0.25 μg/ml). The kcat/Km or kinact/KI (μM(-1) s(-1)) indicated that CMY-33 possesses an ESBL-like spectrum compared to CMY-2 (cefoxitin: 0.2 vs. 0.4; ceftazidime: 0.2 vs. not measurable; cefepime: 0.2 vs. not measurable; tazobactam 0.0018 vs. 0.0009). Using molecular modeling, we show that a widened active site (∼4 Å shift) may play a significant role in enhancing cefepime hydrolysis. This is the first in vivo demonstration of a pAmpC that under cephalosporin treatment expands its substrate spectrum resembling an ESBL. The prevalence of CMY-2-Ec isolates is rapidly increasing worldwide, therefore awareness that cefepime treatment may select for resistant isolates is critical.

Disease evolution and outcomes in familial AML with germline CEBPA mutations

In-depth molecular investigation of familial leukemia has been limited by the rarity of recognized cases. This study examines the genetic events initiating leukemia and details the clinical progression of disease across multiple families harboring germ-line CEBPA mutations. Clinical data were collected from 10 CEBPA-mutated families, representing 24 members with acute myeloid leukemia (AML). Whole-exome (WES) and deep sequencing were performed to genetically profile tumors and define patterns of clonal evolution. Germline CEBPA mutations clustered within the N-terminal and were highly penetrant, with AML presenting at a median age of 24.5 years (range, 1.75-46 years). In all diagnostic tumors tested (n = 18), double CEBPA mutations (CEBPAdm) were detected, with acquired (somatic) mutations preferentially targeting the C-terminal. Somatic CEBPA mutations were unstable throughout the disease course, with different mutations identified at recurrence. Deep sequencing of diagnostic and relapse paired samples confirmed that relapse-associated CEBPA mutations were absent at diagnosis, suggesting recurrence was triggered by novel, independent clones. Integrated WES and deep sequencing subsequently revealed an entirely new complement of mutations at relapse, verifying the presentation of a de novo leukemic episode. The cumulative incidence of relapse in familial AML was 56% at 10 years (n = 11), and 3 patients experienced ≥3 disease episodes over a period of 17 to 20 years. Durable responses to secondary therapies were observed, with prolonged median survival after relapse (8 years) and long-term overall survival (10-year overall survival, 67%). Our data reveal that familial CEBPA-mutated AML exhibits a unique model of disease progression, associated with favorable long-term outcomes.

Mechanisms and evolution of plant resistance to aphids

Aphids are important herbivores of both wild and cultivated plants. Plants rely on unique mechanisms of recognition, signalling and defence to cope with the specialized mode of phloem feeding by aphids. Aspects of the molecular mechanisms underlying aphid-plant interactions are beginning to be understood. Recent advances include the identification of aphid salivary proteins involved in host plant manipulation, and plant receptors involved in aphid recognition. However, a complete picture of aphid-plant interactions requires consideration of the ecological outcome of these mechanisms in nature, and the evolutionary processes that shaped them. Here we identify general patterns of resistance, with a special focus on recognition, phytohormonal signalling, secondary metabolites and induction of plant resistance. We discuss how host specialization can enable aphids to co-opt both the phytohormonal responses and defensive compounds of plants for their own benefit at a local scale. In response, systemically induced resistance in plants is common and often involves targeted responses to specific aphid species or even genotypes. As co-evolutionary adaptation between plants and aphids is ongoing, the stealthy nature of aphid feeding makes both the mechanisms and outcomes of these interactions highly distinct from those of other herbivore-plant interactions. © 2016 Macmillan Publishers Limited.

Molecular diagnostics of myeloproliferative neoplasms

Since the discovery of the JAK2 V617F mutation in the majority of the myeloproliferative neoplasms (MPN) of polycythemia vera, essential thrombocythemia and primary myelofibrosis ten years ago, further MPN-specific mutational events, notably in JAK2 exon 12, MPL exon 10 and CALR exon 9 have been identified. These discoveries have been rapidly incorporated into evolving molecular diagnostic algorithms. Whilst many of these mutations appear to have prognostic implications, establishing MPN diagnosis is of immediate clinical importance with selection, implementation and the continual evaluation of the appropriate laboratory methodology to achieve this diagnosis similarly vital. The advantages and limitations of these approaches in identifying and quantitating the common MPN-associated mutations are considered herein with particular regard to their clinical utility. The evolution of molecular diagnostic applications and platforms has occurred in parallel with the discovery of MPN-associated mutations, and it therefore appears likely that emerging technologies such as next-generation sequencing and digital PCR will in the future play an increasing role in the molecular diagnosis of MPN. Accepted for publication 30 April 2015 doi:10.1111/ejh.12578