Search

The relative importance of genetic drift and local adaptation in facilitating speciation remains unclear. This is particularly true for seabirds, which can disperse over large geographic distances, providing opportunities for intermittent gene flow among distant colonies that span the temperature and salinity gradients of the oceans. Here, we delve into the genomic basis of adaptation and speciation of banded penguins, Galápagos (Spheniscus mendiculus), Humboldt (Spheniscus humboldti), Magellanic (Spheniscus magellanicus), and African penguins (Spheniscus demersus), by analyzing 114 genomes from the main 16 breeding colonies. We aim to identify the molecular mechanism and genomic adaptive traits that have facilitated their diversifications. Through positive selection and gene family expansion analyses, we identified candidate genes that may be related to reproductive isolation processes mediated by ecological thermal niche divergence. We recover signals of positive selection on key loci associated with spermatogenesis, especially during the recent peripatric divergence of the Galápagos penguin from the Humboldt penguin. High temperatures in tropical habitats may have favored selection on loci associated with spermatogenesis to maintain sperm viability, leading to reproductive isolation among young species. Our results suggest that genome-wide selection on loci associated with molecular pathways that underpin thermoregulation, osmoregulation, hypoxia, and social behavior appears to have been crucial in local adaptation of banded penguins. Overall, these results contribute to our understanding of how the complexity of biotic, but especially abiotic, factors, along with the high dispersal capabilities of these marine species, may promote both neutral and adaptive lineage divergence even in the presence of gene flow., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Malaria parasites are known to infect a variety of vertebrate hosts, including ungulates. However, ungulates of Amazonia have not been investigated. We report for the first time, the presence of parasite lineages closely related to Plasmodium odocoilei clade 1 and clade 2 in free-ranging South American red-brocket deer (Mazama americana; 44.4%, 4/9) and gray-brocket deer (Mazama nemorivaga; 50.0%, 1/2). We performed PCR-based analysis of blood samples from 47 ungulates of five different species collected during subsistence hunting by an indigenous community in the Peruvian Amazon. We detected Plasmodium malariae/brasilianum lineage in a sample from red-brocket deer. However, no parasite DNA was detected in collared peccary (Pecari tajacu; 0.0%, 0/10), white-lipped peccary (Tayassu pecari; 0.0%, 0/15), and tapir (Tapirus terrestris; 0.0%, 0/11). Concordant phylogenetic analyses suggested a possible co-evolutionary relationship between the Plasmodium lineages found in American deer and their hosts., Competing Interests: Declaration of competing interest Nothing to declare., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Hibernation in bears involves a suite of metabolical and physiological changes, including the onset of insulin resistance, that are driven in part by sweeping changes in gene expression in multiple tissues. Feeding bears glucose during hibernation partially restores active season physiological phenotypes, including partial resensitization to insulin, but the molecular mechanisms underlying this transition remain poorly understood. Here, we analyze tissue-level gene expression in adipose, liver, and muscle to identify genes that respond to midhibernation glucose feeding and thus potentially drive postfeeding metabolical and physiological shifts. We show that midhibernation feeding stimulates differential expression in all analyzed tissues of hibernating bears and that a subset of these genes responds specifically by shifting expression toward levels typical of the active season. Inferences of upstream regulatory molecules potentially driving these postfeeding responses implicate peroxisome proliferator-activated receptor gamma (PPARG) and other known regulators of insulin sensitivity, providing new insight into high-level regulatory mechanisms involved in shifting metabolic phenotypes between hibernation and active states.

Archaic admixture has had a substantial impact on human evolution with multiple events across different clades, including from extinct hominins such as Neanderthals and Denisovans into modern humans. In great apes, archaic admixture has been identified in chimpanzees and bonobos but the possibility of such events has not been explored in other species. Here, we address this question using high-coverage whole-genome sequences from all four extant gorilla subspecies, including six newly sequenced eastern gorillas from previously unsampled geographic regions. Using approximate Bayesian computation with neural networks to model the demographic history of gorillas, we find a signature of admixture from an archaic 'ghost' lineage into the common ancestor of eastern gorillas but not western gorillas. We infer that up to 3% of the genome of these individuals is introgressed from an archaic lineage that diverged more than 3 million years ago from the common ancestor of all extant gorillas. This introgression event took place before the split of mountain and eastern lowland gorillas, probably more than 40 thousand years ago and may have influenced perception of bitter taste in eastern gorillas. When comparing the introgression landscapes of gorillas, humans and bonobos, we find a consistent depletion of introgressed fragments on the X chromosome across these species. However, depletion in protein-coding content is not detectable in eastern gorillas, possibly as a consequence of stronger genetic drift in this species., (© 2023. The Author(s).)

Objectives: Complex physiological adaptations often involve the coordination of molecular responses across multiple tissues. Establishing transcriptomic resources for non-traditional model organisms with phenotypes of interest can provide a foundation for understanding the genomic basis of these phenotypes, and the degree to which these resemble, or contrast, those of traditional model organisms. Here, we present a one-of-a-kind gene expression dataset generated from multiple tissues of two hibernating brown bears (Ursus arctos)., Data Description: This dataset is comprised of 26 samples collected from 13 tissues of two hibernating brown bears. These samples were collected opportunistically and are typically not possible to attain, resulting in a highly unique and valuable gene expression dataset. In combination with previously published datasets, this new transcriptomic resource will facilitate detailed investigation of hibernation physiology in bears, and the potential to translate aspects of this biology to treat human disease., (© 2023. The Author(s).)

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals 1 . These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample., (© 2023. The Author(s).)

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology., (© 2023. Springer Nature Limited.)

Understanding how metabolic reprogramming happens in cells will aid the progress in the treatment of a variety of metabolic disorders. Brown bears undergo seasonal shifts in insulin sensitivity, including reversible insulin resistance in hibernation. We performed RNA-sequencing on brown bear adipocytes and proteomics on serum to identify changes possibly responsible for reversible insulin resistance. We observed dramatic transcriptional changes, which depended on both the cell and serum season of origin. Despite large changes in adipocyte gene expression, only changes in eight circulating proteins were identified as related to the seasonal shifts in insulin sensitivity, including some that have not previously been associated with glucose homeostasis. The identified serum proteins may be sufficient for shifting hibernation adipocytes to an active-like state., Competing Interests: The MacCoss Lab at the University of Washington (members G.E.M., J.P., and M.J.M.) has a sponsored research agreement with Thermo Fisher Scientific, the manufacturer of the instrumentation used in this research. Moreover, M.J.M. is a paid consultant for Thermo Fisher Scientific. The remaining authors declare no competing interests., (© 2022 The Author(s).)

The SARS-CoV-2 pandemic has differentially impacted populations across race and ethnicity. A multi-omic approach represents a powerful tool to examine risk across multi-ancestry genomes. We leverage a pandemic tracking strategy in which we sequence viral and host genomes and transcriptomes from nasopharyngeal swabs of 1049 individuals (736 SARS-CoV-2 positive and 313 SARS-CoV-2 negative) and integrate them with digital phenotypes from electronic health records from a diverse catchment area in Northern California. Genome-wide association disaggregated by admixture mapping reveals novel COVID-19-severity-associated regions containing previously reported markers of neurologic, pulmonary and viral disease susceptibility. Phylodynamic tracking of consensus viral genomes reveals no association with disease severity or inferred ancestry. Summary data from multiomic investigation reveals metagenomic and HLA associations with severe COVID-19. The wealth of data available from residual nasopharyngeal swabs in combination with clinical data abstracted automatically at scale highlights a powerful strategy for pandemic tracking, and reveals distinct epidemiologic, genetic, and biological associations for those at the highest risk., (© 2022. The Author(s).)

PRDM9 drives recombination hotspots in some mammals, including mice and apes. Non-functional orthologs of PRDM9 are present in a wide variety of vertebrates, but why it is functionally maintained in some lineages is not clear. One possible explanation is that PRDM9 plays a role in ensuring that meiosis is successful. During meiosis, available DNA may be a limiting resource given the tight packaging of chromosomes and could lead to competition between two key processes: meiotic transcription and recombination. Here we explore this potential competition and the role that PRDM9 might play in their interaction. Leveraging existing mouse genomic data, we use resampling schemes that simulate shuffled features along the genome and models that account for the rarity of features in the genome, to test if PRDM9 influences interactions between recombination hotspots and meiotic transcription in a whole genome framework. We also explored possible DNA sequence motifs associated to clusters of hotspots not tied to transcription or PRDM9. We find evidence of competition between meiotic transcription and recombination, with PRDM9 appearing to relocate recombination to avoid said conflict. We also find that retrotransposons may be playing a role in directing hotspots in the absence of other factors., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Anthurium andraeanum (Hort.) is an important ornamental in the tropical cut-flower industry. However, there is currently insufficient information to establish a clear connection between the genetic model(s) proposed and the putative genes involved in the differentiation between colors. In this study, 18 cDNA libraries related to the spathe color and developmental stages of A. andraeanum were characterized by transcriptome sequencing (RNA-seq). For the de novo transcriptome, a total of 114,334,082 primary sequence reads were obtained from the Illumina sequencer and were assembled into 151,652 unigenes. Approximately 58,476 transcripts were generated and used for comparative transcriptome analysis between three cultivars that differ in spathe color ('Sasha' (white), 'Honduras' (red), and 'Rapido' (purple)). A large number of differentially expressed genes (8,324), potentially involved in multiple biological and metabolic pathways, were identified, including genes in the flavonoid and anthocyanin biosynthetic pathways. Our results showed that the chalcone isomerase (CHI) gene presented the strongest evidence for an association with differences in color and the highest correlation with other key genes (flavanone 3-hydroxylase (F3H), flavonoid 3'5' hydroxylase (F3'5'H)/ flavonoid 3'-hydroxylase (F3'H), and leucoanthocyanidin dioxygenase (LDOX)) in the anthocyanin pathway. We also identified a differentially expressed cytochrome P450 gene in the late developmental stage of the purple spathe that appeared to determine the difference between the red- and purple-colored spathes. Furthermore, transcription factors related to putative MYB-domain protein that may control anthocyanin pathway were identified through a weighted gene co-expression network analysis (WGCNA). The results provided basic sequence information for future research on spathe color, which have important implications for this ornamental breeding strategies., Competing Interests: The authors have declared that no competing interests exist.

Predicting the potential fate of a species in the face of climate change requires knowing the distribution of molecular adaptations across the geographic range of the species. In this work, we analysed 79 genomes of Theobroma cacao, an Amazonian tree known for the fruit from which chocolate is produced, to evaluate how local and regional molecular signatures of adaptation are distributed across the natural range of the species. We implemented novel techniques that incorporate summary statistics from multiple selection scans to infer selective sweeps. The majority of the molecular adaptations in the genome are not shared among populations. We show that ~71.5% of genes under selection also show significant associations with changes in environmental variables. Our results support the interpretation that these genes contribute to local adaptation of the populations in response to abiotic factors. We also found strong patterns of molecular adaptation in a diverse array of disease resistance genes (6.5% of selective sweeps), suggesting that differential adaptation to pathogens also contributes significantly to local adaptations. Our results are consistent with the interpretation that local selective pressures are more important than regional selective pressures in explaining adaptation across the range of a species., (© 2020 John Wiley & Sons Ltd.)

Background: Recombination plays an important evolutionary role by breaking up haplotypes and shuffling genetic variation. This process impacts the ability of selection to eliminate deleterious mutations or increase the frequency of beneficial mutations in a population. To understand the role of recombination generating and maintaining haplotypic variation in a population, we can construct fine-scale recombination maps. Such maps have been used to study a variety of model organisms and proven to be informative of how selection and demographics shape species-wide variation. Here we present a fine-scale recombination map for ten populations of Theobroma cacao - a non-model, long-lived, woody crop. We use this map to elucidate the dynamics of recombination rates in distinct populations of the same species, one of which is domesticated., Results: Mean recombination rates in range between 2.5 and 8.6 cM/Mb for most populations of T. cacao with the exception of the domesticated Criollo (525 cM/Mb) and Guianna, a more recently established population (46.5 cM/Mb). We found little overlap in the location of hotspots of recombination across populations. We also found that hotspot regions contained fewer known retroelement sequences than expected and were overrepresented near transcription start and termination sites. We find mutations in FIGL-1, a protein shown to downregulate cross-over frequency in Arabidopsis, statistically associated to higher recombination rates in domesticated Criollo., Conclusions: We generated fine-scale recombination maps for ten populations of Theobroma cacao and used them to understand what processes are associated with population-level variation in this species. Our results provide support to the hypothesis of increased recombination rates in domesticated plants (Criollo population). We propose a testable mechanistic hypothesis for the change in recombination rate in domesticated populations in the form of mutations to a previously identified recombination-suppressing protein. Finally, we establish a number of possible correlates of recombination hotspots that help explain general patterns of recombination in this species.

Background: The genus Elaeis has two species of economic importance for the oil palm agroindustry: Elaeis oleifera (O), native to the Americas, and Elaeis guineensis (G), native to Africa. This work provides to our knowledge, the first association mapping study in an interspecific OxG oil palm population, which shows tolerance to pests and diseases, high oil quality, and acceptable fruit bunch production., Results: Using genotyping-by-sequencing (GBS), we identified a total of 3776 single nucleotide polymorphisms (SNPs) that were used to perform a genome-wide association analysis (GWAS) in 378 OxG hybrid population for 10 agronomic traits. Twelve genomic regions (SNPs) were located near candidate genes implicated in multiple functional categories, such as tissue growth, cellular trafficking, and physiological processes., Conclusions: We provide new insights on genomic regions that mapped on candidate genes involved in plant architecture and yield. These potential candidate genes need to be confirmed for future targeted functional analyses. Associated markers to the traits of interest may be valuable resources for the development of marker-assisted selection in oil palm breeding.

Revealing the mechanisms underlying the reversible physiology of hibernation could have applications to both human and animal health as hibernation is often associated with disease-like states. The present study uses RNA-sequencing to reveal the tissue and seasonal transcriptional changes occurring in grizzly bears ( Ursus arctos horribilis ). Comparing hibernation to other seasons, bear adipose has a greater number of differentially expressed genes than liver and skeletal muscle. During hyperphagia, adipose has more than 900 differentially expressed genes compared to active season. Hibernation is characterized by reduced expression of genes associated with insulin signaling, muscle protein degradation, and urea production, and increased expression within muscle protein anabolic pathways. Across all three tissues we find a subset of shared differentially expressed genes, some of which are uncharacterized, that together may reflect a common regulatory mechanism. The identified gene families could be useful for developing novel therapeutics to treat human and animal diseases., Competing Interests: Competing interestsThe authors declare no competing interests.

Background: Carbonic anhydrase (CA) catalyzes the hydration of CO 2 in the first biochemical step of C 4 photosynthesis, and has been considered a potentially rate-limiting step when CO 2 availability within a leaf is low. Previous work in Zea mays (maize) with a double knockout of the two highest-expressed β-CA genes, CA1 and CA2, reduced total leaf CA activity to less than 3% of wild-type. Surprisingly, this did not limit photosynthesis in maize at ambient or higher CO 2 concentrations. However, the ca1ca2 mutants exhibited reduced rates of photosynthesis at sub-ambient CO 2 , and accumulated less biomass when grown under sub-ambient CO 2 (9.2 Pa). To further clarify the importance of CA for C 4 photosynthesis, we assessed gene expression changes in wild-type, ca1 and ca1ca2 mutants in response to changes in pCO 2 from 920 to 9.2 Pa., Results: Leaf samples from each genotype were collected for RNA-seq analysis at high CO 2 and at two time points after the low CO 2 transition, in order to identify early and longer-term responses to CO 2 deprivation. Despite the existence of multiple isoforms of CA, no other CA genes were upregulated in CA mutants. Although photosynthetic genes were downregulated in response to low CO 2 , differential expression was not observed between genotypes. However, multiple indicators of carbon starvation were present in the mutants, including amino acid synthesis, carbohydrate metabolism, and sugar signaling. In particular, multiple genes previously implicated in low carbon stress such as asparagine synthetase, amino acid transporters, trehalose-6-phosphate synthase, as well as many transcription factors, were strongly upregulated. Furthermore, genes in the CO 2 stomatal signaling pathway were differentially expressed in the CA mutants under low CO 2 ., Conclusions: Using a transcriptomic approach, we showed that carbonic anhydrase mutants do not compensate for the lack of CA activity by upregulating other CA or photosynthetic genes, but rather experienced extreme carbon stress when grown under low CO 2 . Our results also support a role for CA in the CO 2 stomatal signaling pathway. This study provides insight into the importance of CA for C 4 photosynthesis and its role in stomatal signaling.

Puccinia striiformis f. sp. tritici causes devastating stripe (yellow) rust on wheat and P. striiformis f. sp. hordei causes stripe rust on barley. Several P. striiformis f. sp. tritici genomes are available, but no P. striiformis f. sp. hordei genome is available. More genomes of P. striiformis f. sp. tritici and P. striiformis f. sp. hordei are needed to understand the genome evolution and molecular mechanisms of their pathogenicity. We sequenced P. striiformis f. sp. tritici isolate 93-210 and P. striiformis f. sp. hordei isolate 93TX-2, using PacBio and Illumina technologies and RNA sequencing. Their genomic sequences were assembled to contigs with high continuity and showed significant structural differences. The circular mitochondria genomes of both were complete. These genomes provide high-quality resources for deciphering the genomic basis of rapid evolution and host adaptation, identifying genes for avirulence and other important traits, and studying host-pathogen interactions.

As the source of chocolate, cacao has become one of the most important crops in the world. The identification of molecular markers to understand the demographic history, genetic diversity and population structure plays a pivotal role in cacao breeding programs. Here, we report the use of a modified genotyping-by-sequencing (GBS) approach for large-scale single nucleotide polymorphism (SNP) discovery and allele ancestry mapping. We identified 12,357 bi-allelic SNPs after filtering, of which, 7,009 variants were ancestry informative. The GBS approach proved to be rapid, cost-effective, and highly informative for ancestry assignment in this species., Competing Interests: Jaime Osorio is an employee of Agrosavia (Colombia) that provided with some of the funds necessary for the completion of this work contribution of Agrosavia does adhere to PLOS policies of data sharing and materials.

Domestication has had a strong impact on the development of modern societies. We sequenced 200 genomes of the chocolate plant Theobroma cacao L. to show for the first time to our knowledge that a single population, the Criollo population, underwent strong domestication ~3600 years ago (95% CI: 2481-13,806 years ago). We also show that during the process of domestication, there was strong selection for genes involved in the metabolism of the colored protectants anthocyanins and the stimulant theobromine, as well as disease resistance genes. Our analyses show that domesticated populations of T. cacao (Criollo) maintain a higher proportion of high-frequency deleterious mutations. We also show for the first time the negative consequences of the increased accumulation of deleterious mutations during domestication on the fitness of individuals (significant reduction in kilograms of beans per hectare per year as Criollo ancestry increases, as estimated from a GLM, P  = 0.000425)., Competing Interests: Co-authors Donald Livingstone III, Conrad Stack, Alberto Romero, Stefan Royaert, Osman Gutierrez, and Juan C. Motamayor are or have been employees of the funding agency MARS Inc. The authors declare no other competing financial or non-financial interests.

Background: Plant fungal pathogens can rapidly evolve and adapt to new environmental conditions in response to sudden changes of host populations in agro-ecosystems. However, the genomic basis of their host adaptation, especially at the forma specialis level, remains unclear., Results: We sequenced two isolates each representing Puccinia striiformis f. sp. tritici (Pst) and P. striiformis f. sp. hordei (Psh), different formae speciales of the stripe rust fungus P. striiformis highly adapted to wheat and barley, respectively. The divergence of Pst and Psh, estimated to start 8.12 million years ago, has been driven by high nucleotide mutation rates. The high genomic variation within dikaryotic urediniospores of P. striiformis has provided raw genetic materials for genome evolution. No specific gene families have enriched in either isolate, but extensive gene loss events have occurred in both Pst and Psh after the divergence from their most recent common ancestor. A large number of isolate-specific genes were identified, with unique genomic features compared to the conserved genes, including 1) significantly shorter in length; 2) significantly less expressed; 3) significantly closer to transposable elements; and 4) redundant in pathways. The presence of specific genes in one isolate (or forma specialis) was resulted from the loss of the homologues in the other isolate (or forma specialis) by the replacements of transposable elements or losses of genomic fragments. In addition, different patterns and numbers of telomeric repeats were observed between the isolates., Conclusions: Host adaptation of P. striiformis at the forma specialis level is a complex pathogenic trait, involving not only virulence-related genes but also other genes. Gene loss, which might be adaptive and driven by transposable element activities, provides genomic basis for host adaptation of different formae speciales of P. striiformis.

Aneuploidy is prevalent in human embryos and is the leading cause of pregnancy loss. Many aneuploidies arise during oogenesis, increasing with maternal age. Superimposed on these meiotic aneuploidies are frequent errors occurring during early mitotic divisions, contributing to widespread chromosomal mosaicism. Here we reanalyzed a published dataset comprising preimplantation genetic testing for aneuploidy in 24 653 blastomere biopsies from day-3 cleavage-stage embryos, as well as 17 051 trophectoderm biopsies from day-5 blastocysts. We focused on complex abnormalities that affected multiple chromosomes simultaneously, seeking insights into their formation. In addition to well-described patterns such as triploidy and haploidy, we identified 4.7% of blastomeres possessing characteristic hypodiploid karyotypes. We inferred this signature to have arisen from tripolar chromosome segregation in normally fertilized diploid zygotes or their descendant diploid cells. This could occur via segregation on a tripolar mitotic spindle or by rapid sequential bipolar mitoses without an intervening S-phase. Both models are consistent with time-lapse data from an intersecting set of 77 cleavage-stage embryos, which were enriched for the tripolar signature among embryos exhibiting abnormal cleavage. The tripolar signature was strongly associated with common maternal genetic variants spanning the centrosomal regulator PLK4, driving the association we previously reported with overall mitotic errors. Our findings are consistent with the known capacity of PLK4 to induce tripolar mitosis or precocious M-phase upon dysregulation. Together, our data support tripolar chromosome segregation as a key mechanism generating complex aneuploidy in cleavage-stage embryos and implicate maternal genotype at a quantitative trait locus spanning PLK4 as a factor influencing its occurrence.

Stripe (yellow) rust, caused by Puccinia striiformis f. sp. tritici ( Pst ), is one of the most destructive diseases of wheat worldwide. Planting resistant cultivars is an effective way to control this disease, but race-specific resistance can be overcome quickly due to the rapid evolving Pst population. Studying the pathogenicity mechanisms is critical for understanding how Pst virulence changes and how to develop wheat cultivars with durable resistance to stripe rust. We re-sequenced 7 Pst isolates and included additional 7 previously sequenced isolates to represent balanced virulence/avirulence profiles for several avirulence loci in seretome analyses. We observed an uneven distribution of heterozygosity among the isolates. Secretome comparison of Pst with other rust fungi identified a large portion of species-specific secreted proteins, suggesting that they may have specific roles when interacting with the wheat host. Thirty-two effectors of Pst were identified from its secretome. We identified candidates for Avr genes corresponding to six Yr genes by correlating polymorphisms for effector genes to the virulence/avirulence profiles of the 14 Pst isolates. The putative AvYr76 was present in the avirulent isolates, but absent in the virulent isolates, suggesting that deleting the coding region of the candidate avirulence gene has produced races virulent to resistance gene Yr76 . We conclude that incorporating avirulence/virulence phenotypes into correlation analysis with variations in genomic structure and secretome, particularly presence/absence polymorphisms of effectors, is an efficient way to identify candidate Avr genes in Pst . The candidate effector genes provide a rich resource for further studies to determine the evolutionary history of Pst populations and the co-evolutionary arms race between Pst and wheat. The Avr candidates identified in this study will lead to cloning avirulence genes in Pst , which will enable us to understand molecular mechanisms underlying Pst -wheat interactions, to determine the effectiveness of resistance genes and further to develop durable resistance to stripe rust.

Cacao ( Theobroma cacao L.) is an important cash crop in tropical regions around the world and has a rich agronomic history in South America. As a key component in the cosmetic and confectionary industries, millions of people worldwide use products made from cacao, ranging from shampoo to chocolate. An Illumina Infinity II array was created using 13,530 SNPs identified within a small diversity panel of cacao. Of these SNPs, 12,643 derive from variation within annotated cacao genes. The genotypes of 3,072 trees were obtained, including two mapping populations from Ecuador. High-density linkage maps for these two populations were generated and compared to the cacao genome assembly. Phenotypic data from these populations were combined with the linkage maps to identify the QTLs for yield and disease resistance.

Gardnerella vaginalis has long been associated with bacterial vaginosis, a condition that increases the risk of women to preterm birth, sexually transmitted infections, and other adverse sequelae. However, G. vaginalis is also commonly found in healthy asymptomatic women of all ages. This raises the question if genetic differences among strains might distinguish potentially pathogenic from commensal strains. To disentangle the diversity of  G. vaginalis, we invoked the concept of ecotypes-lineages of genetically and ecologically distinct strains within a named species-to better understand their evolutionary history and identify functional characteristics. We compared the genomes of G. vaginalis to six species in the closely related Bifidobacterium genus and found that G. vaginalis has a large accessory genome relative to Bifidobacterium , including many unique genes possibly involved in metabolism, drug resistance, and virulence. We then performed a comparative genomic analysis of 35 strains of G. vaginalis to infer a phylogeny based on the combined analysis of the core genome, using nucleotide substitution models, and the accessory genome, using gene gain/loss models. With the inferred tree topology, we performed comparisons of functional gene content among lineages that diverged at varying depths in the phylogeny and found significant differences in the representation of genes putatively involved in pathogenicity. Our functional enrichment analysis suggests that some lineages of G. vaginalis may possess enhanced pathogenic capabilities, including genes involved in mucus degradation like sialidases, while others may be commensal strains, lacking many of these pathogenic capabilities. The combined phylogenetic evidence and functional enrichment analysis allowed us to identify distinct ecotypes that have evolved in  G. vaginalis as the result of the differential gene gain/loss for specific functions, including the capability to cause disease. We finally discuss how this analysis framework could be used to gain insight into the etiology of bacterial vaginosis and improve diagnosis.

Natural transformation in the Gram-positive pathogen Streptococcus pneumoniae occurs when cells become "competent," a state that is induced in response to high extracellular concentrations of a secreted peptide signal called competence stimulating peptide (CSP) encoded by the comC locus. Two main CSP signal types (pherotypes) are known to dominate the pherotype diversity across strains. Using 4,089 fully sequenced pneumococcal genomes, we confirm that pneumococcal populations are highly genetically structured and that there is significant variation among diverged populations in pherotype frequencies; most carry only a single pherotype. Moreover, we find that the relative frequencies of the two dominant pherotypes significantly vary within a small range across geographical sites. It has been variously proposed that pherotypes either promote genetic exchange among cells expressing the same pherotype, or conversely that they promote recombination between strains bearing different pherotypes. We attempt to distinguish these hypotheses using a bioinformatics approach by estimating recombination frequencies within and between pherotypes across 4,089 full genomes. Despite underlying population structure, we observe extensive recombination between populations; additionally, we found significantly higher (although marginal) rates of genetic exchange between strains expressing different pherotypes than among isolates carrying the same pherotype. Our results indicate that pherotypes do not restrict, and may even slightly facilitate, recombination between strains; however, these marginal effects suggest the more likely possibility that the cause of CSP polymorphism lies outside of its effects on transformation. Our results suggest that the CSP balanced polymorphism does not causally underlie population differentiation. Therefore, when strains carrying different pherotypes encounter one another during cocolonization, genetic exchange can occur without restriction., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The human DARC (Duffy antigen receptor for chemokines) gene encodes a membrane-bound chemokine receptor crucial for the infection of red blood cells by Plasmodium vivax, a major causative agent of malaria. Of the three major allelic classes segregating in human populations, the FY*O allele has been shown to protect against P. vivax infection and is at near fixation in sub-Saharan Africa, while FY*B and FY*A are common in Europe and Asia, respectively. Due to the combination of strong geographic differentiation and association with malaria resistance, DARC is considered a canonical example of positive selection in humans. Despite this, details of the timing and mode of selection at DARC remain poorly understood. Here, we use sequencing data from over 1,000 individuals in twenty-one human populations, as well as ancient human genomes, to perform a fine-scale investigation of the evolutionary history of DARC. We estimate the time to most recent common ancestor (TMRCA) of the most common FY*O haplotype to be 42 kya (95% CI: 34-49 kya). We infer the FY*O null mutation swept to fixation in Africa from standing variation with very low initial frequency (0.1%) and a selection coefficient of 0.043 (95% CI:0.011-0.18), which is among the strongest estimated in the human genome. We estimate the TMRCA of the FY*A mutation in non-Africans to be 57 kya (95% CI: 48-65 kya) and infer that, prior to the sweep of FY*O, all three alleles were segregating in Africa, as highly diverged populations from Asia and ≠Khomani San hunter-gatherers share the same FY*A haplotypes. We test multiple models of admixture that may account for this observation and reject recent Asian or European admixture as the cause.

Frosty pod rot (FPR) disease on cocoa, caused by Moniliophthora roreri, is one of the most devastating cocoa disease in the Western Hemisphere. In Colombia, the disease is particularly severe in the Magdalena Valley, which is considered the possible center of origin for the pathogen species. We analyzed the genetic diversity of isolates from the departments of Santander, Antioquia, Tolima, and Huila in Colombia using 23 simple-sequence repeats (SSR) markers. In total, 117 different multilocus genotypes were found among 120 isolates, each one representing a unique haplotype. High mutation rates in the SSR and gene flow can explain the high levels of diversity. Also, the observed and standardized indexes of association (I A and řd) indicate that the populations of M. roreri are clonal. Furthermore, given the high haplotype diversity and the significant linkage disequilibrium observed, we hypothesize that M. roreri could be a primarily asexual species undergoing sporadic recombination or partial recombination through parasexuality. A Bayesian clustering analysis implemented by STRUCTURE showed that the most probable number of genetic groups in the data was three, confirming the geographical differentiation among isolates. Similar results were obtained by a discriminant analysis of principal components, a principal coordinate analysis, and a neighbor-joining tree from microsatellite loci base on Nei distance. Cacao genotypes and environmental variables did contribute to the genetic differentiation of the groups. We discuss how this information could be used to improve the management of FPR at the regional level.

Transmission-blocking (TB) vaccines are considered an important tool for malaria control and elimination. Among all the antigens characterized as TB vaccines against Plasmodium vivax, the ookinete surface proteins Pvs28 and Pvs25 are leading candidates. These proteins likely originated by a gene duplication event that took place before the radiation of the known Plasmodium species to primates. We report an evolutionary genetic analysis of a worldwide sample of pvs28 and pvs25 alleles. Our results show that both genes display low levels of genetic polymorphism when compared to the merozoite surface antigens AMA-1 and MSP-1; however, both ookinete antigens can be as polymorphic as other merozoite antigens such as MSP-8 and MSP-10. We found that parasite populations in Asia and the Americas are geographically differentiated with comparable levels of genetic diversity and specific amino acid replacements found only in the Americas. Furthermore, the observed variation was mainly accumulated in the EGF2- and EGF3-like domains for P. vivax in both proteins. This pattern was shared by other closely related non-human primate parasites such as Plasmodium cynomolgi, suggesting that it could be functionally important. In addition, examination with a suite of evolutionary genetic analyses indicated that the observed patterns are consistent with positive natural selection acting on Pvs28 and Pvs25 polymorphisms. The geographic pattern of genetic differentiation and the evidence for positive selection strongly suggest that the functional consequences of the observed polymorphism should be evaluated during development of TBVs that include Pvs25 and Pvs28.

Reduced representation sequencing methods such as genotyping-by-sequencing (GBS) enable low-cost measurement of genetic variation without the need for a reference genome assembly. These methods are widely used in genetic mapping and population genetics studies, especially with non-model organisms. Variant calling error rates, however, are higher in GBS than in standard sequencing, in particular due to restriction site polymorphisms, and few computational tools exist that specifically model and correct these errors. We developed a statistical method to remove errors caused by restriction site polymorphisms, implemented in the software package GBStools. We evaluated it in several simulated data sets, varying in number of samples, mean coverage and population mutation rate, and in two empirical human data sets (N = 8 and N = 63 samples). In our simulations, GBStools improved genotype accuracy more than commonly used filters such as Hardy-Weinberg equilibrium p-values. GBStools is most effective at removing genotype errors in data sets over 100 samples when coverage is 40X or higher, and the improvement is most pronounced in species with high genomic diversity. We also demonstrate the utility of GBS and GBStools for human population genetic inference in Argentine populations and reveal widely varying individual ancestry proportions and an excess of singletons, consistent with recent population growth.

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model., (Copyright © 2015, American Association for the Advancement of Science.)

Plasmodium vivax is the most prevalent human malaria parasite outside of Africa. Yet, studies aimed to identify genes with signatures consistent with natural selection are rare. Here, we present a comparative analysis of the pattern of genetic variation of five sequenced isolates of P. vivax and its divergence with two closely related species, Plasmodium cynomolgi and Plasmodium knowlesi, using a set of orthologous genes. In contrast to Plasmodium falciparum, the parasite that causes the most lethal form of human malaria, we did not find significant constraints on the evolution of synonymous sites genome wide in P. vivax. The comparative analysis of polymorphism and divergence across loci allowed us to identify 87 genes with patterns consistent with positive selection, including genes involved in the "exportome" of P. vivax, which are potentially involved in evasion of the host immune system. Nevertheless, we have found a pattern of polymorphism genome wide that is consistent with a significant amount of constraint on the replacement changes and prevalent negative selection. Our analyses also show that silent polymorphism tends to be larger toward the ends of the chromosomes, where many genes involved in antigenicity are located, suggesting that natural selection acts not only by shaping the patterns of variation within the genes but it also affects genome organization., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The New World Arctic, the last region of the Americas to be populated by humans, has a relatively well-researched archaeology, but an understanding of its genetic history is lacking. We present genome-wide sequence data from ancient and present-day humans from Greenland, Arctic Canada, Alaska, Aleutian Islands, and Siberia. We show that Paleo-Eskimos (~3000 BCE to 1300 CE) represent a migration pulse into the Americas independent of both Native American and Inuit expansions. Furthermore, the genetic continuity characterizing the Paleo-Eskimo period was interrupted by the arrival of a new population, representing the ancestors of present-day Inuit, with evidence of past gene flow between these lineages. Despite periodic abandonment of major Arctic regions, a single Paleo-Eskimo metapopulation likely survived in near-isolation for more than 4000 years, only to vanish around 700 years ago., (Copyright © 2014, American Association for the Advancement of Science.)

Clovis, with its distinctive biface, blade and osseous technologies, is the oldest widespread archaeological complex defined in North America, dating from 11,100 to 10,700 (14)C years before present (bp) (13,000 to 12,600 calendar years bp). Nearly 50 years of archaeological research point to the Clovis complex as having developed south of the North American ice sheets from an ancestral technology. However, both the origins and the genetic legacy of the people who manufactured Clovis tools remain under debate. It is generally believed that these people ultimately derived from Asia and were directly related to contemporary Native Americans. An alternative, Solutrean, hypothesis posits that the Clovis predecessors emigrated from southwestern Europe during the Last Glacial Maximum. Here we report the genome sequence of a male infant (Anzick-1) recovered from the Anzick burial site in western Montana. The human bones date to 10,705 ± 35 (14)C years bp (approximately 12,707-12,556 calendar years bp) and were directly associated with Clovis tools. We sequenced the genome to an average depth of 14.4× and show that the gene flow from the Siberian Upper Palaeolithic Mal'ta population into Native American ancestors is also shared by the Anzick-1 individual and thus happened before 12,600 years bp. We also show that the Anzick-1 individual is more closely related to all indigenous American populations than to any other group. Our data are compatible with the hypothesis that Anzick-1 belonged to a population directly ancestral to many contemporary Native Americans. Finally, we find evidence of a deep divergence in Native American populations that predates the Anzick-1 individual.

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to -9.5 × 10(-4) (95% CI -1.1 × 10(-3) to -6.8 × 10(-4)); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620-1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Background: We address the task of extracting accurate haplotypes from genotype data of individuals of large F1 populations for mapping studies. While methods for inferring parental haplotype assignments on large F1 populations exist in theory, these approaches do not work in practice at high levels of accuracy., Results: We have designed iXora (Identifying crossovers and recombining alleles), a robust method for extracting reliable haplotypes of a mapping population, as well as parental haplotypes, that runs in linear time. Each allele in the progeny is assigned not just to a parent, but more precisely to a haplotype inherited from the parent. iXora shows an improvement of at least 15% in accuracy over similar systems in literature. Furthermore, iXora provides an easy-to-use, comprehensive environment for association studies and hypothesis checking in populations of related individuals., Conclusions: iXora provides detailed resolution in parental inheritance, along with the capability of handling very large populations, which allows for accurate haplotype extraction and trait association. iXora is available for non-commercial use from http://researcher.ibm.com/project/3430.

Streptococcus mutans is widely recognized as one of the key etiological agents of human dental caries. Despite its role in this important disease, our present knowledge of gene content variability across the species and its relationship to adaptation is minimal. Estimates of its demographic history are not available. In this study, we generated genome sequences of 57 S. mutans isolates, as well as representative strains of the most closely related species to S. mutans (S. ratti, S. macaccae, and S. criceti), to identify the overall structure and potential adaptive features of the dispensable and core components of the genome. We also performed population genetic analyses on the core genome of the species aimed at understanding the demographic history, and impact of selection shaping its genetic variation. The maximum gene content divergence among strains was approximately 23%, with the majority of strains diverging by 5-15%. The core genome consisted of 1,490 genes and the pan-genome approximately 3,296. Maximum likelihood analysis of the synonymous site frequency spectrum (SFS) suggested that the S. mutans population started expanding exponentially approximately 10,000 years ago (95% confidence interval [CI]: 3,268-14,344 years ago), coincidental with the onset of human agriculture. Analysis of the replacement SFS indicated that a majority of these substitutions are under strong negative selection, and the remainder evolved neutrally. A set of 14 genes was identified as being under positive selection, most of which were involved in either sugar metabolism or acid tolerance. Analysis of the core genome suggested that among 73 genes present in all isolates of S. mutans but absent in other species of the mutans taxonomic group, the majority can be associated with metabolic processes that could have contributed to the successful adaptation of S. mutans to its new niche, the human mouth, and with the dietary changes that accompanied the origin of agriculture.

Background: Ancestral recombinations graph (ARG) is a topological structure that captures the relationship between the extant genomic sequences in terms of genetic events including recombinations. IRiS is a system that estimates the ARG on sequences of individuals, at genomic scales, capturing the relationship between these individuals of the species. Recently, this system was used to estimate the ARG of the recombining X Chromosome of a collection of human populations using relatively dense, bi-allelic SNP data., Results: While the ARG is a natural model for capturing the inter-relationship between a single chromosome of the individuals of a species, it is not immediately apparent how the model can utilize whole-genome (across chromosomes) diploid data. Also, the sheer complexity of an ARG structure presents a challenge to graph visualization techniques. In this paper we examine the ARG reconstruction for (1) genome-wide or multiple chromosomes, (2) multi-allelic and (3) extremely sparse data. To aid in the visualization of the results of the reconstructed ARG, we additionally construct a much simplified topology, a classification tree, suggested by the ARG.As the test case, we study the problem of extracting the relationship between populations of Theobroma cacao. The chocolate tree is an outcrossing species in the wild, due to self-incompatibility mechanisms at play. Thus a principled approach to understanding the inter-relationships between the different populations must take the shuffling of the genomic segments into account. The polymorphisms in the test data are short tandem repeats (STR) and are multi-allelic (sometimes as high as 30 distinct possible values at a locus). Each is at a genomic location that is bilaterally transmitted, hence the ARG is a natural model for this data. Another characteristic of this plant data set is that while it is genome-wide, across 10 linkage groups or chromosomes, it is very sparse, i.e., only 96 loci from a genome of approximately 400 megabases. The results are visualized both as MDS plots and as classification trees. To evaluate the accuracy of the ARG approach, we compare the results with those available in literature., Conclusions: We have extended the ARG model to incorporate genome-wide (ensemble of multiple chromosomes) data in a natural way. We present a simple scheme to implement this in practice. Finally, this is the first time that a plant population data set is being studied by estimating its underlying ARG. We demonstrate an overall precision of 0.92 and an overall recall of 0.93 of the ARG-based classification, with respect to the gold standard. While we have corroborated the classification of the samples with that in literature, this opens the door to other potential studies that can be made on the ARG.

Background: Timing the origin of human malarias has been a focus of great interest. Previous studies on the mitochondrial genome concluded that Plasmodium in primates, including those parasitic to humans, radiated relatively recently during a process where host switches were common. Those investigations, however, assumed constant rate of evolution and tightly bound (fixed) calibration points based on host fossils or host distribution. We investigate the effect of such assumptions using different molecular dating methods. We include parasites from Lemuroidea since their distribution provides an external validation to time estimates allowing us to disregard scenarios that cannot explain their introduction in Madagascar., Results: We reject the assumption that the Plasmodium mitochondrial genome, as a unit or each gene separately, evolves at a constant rate. Our analyses show that Lemuroidea parasites are a monophyletic group that shares a common ancestor with all Catarrhini malarias except those related to P. falciparum. However, we found no evidence that this group of parasites branched with their hosts early in the evolution of primates. We applied relaxed clock methods and different calibrations points to explore the origin of primate malarias including those found in African apes. We showed that previous studies likely underestimated the origin of malarial parasites in primates., Conclusions: The use of fossils from the host as absolute calibration and the assumption of a strict clock likely underestimate time when performing molecular dating analyses on malarial parasites. Indeed, by exploring different calibration points, we found that the time for the radiation of primate parasites may have taken place in the Eocene, a time consistent with the radiation of African anthropoids. The radiation of the four human parasite lineages was part of such events. The time frame estimated in this investigation, together with our phylogenetic analyses, made plausible a scenario where gorillas and humans acquired malaria from a Pan lineage.

Whole-genome sequencing harbors unprecedented potential for characterization of individual and family genetic variation. Here, we develop a novel synthetic human reference sequence that is ethnically concordant and use it for the analysis of genomes from a nuclear family with history of familial thrombophilia. We demonstrate that the use of the major allele reference sequence results in improved genotype accuracy for disease-associated variant loci. We infer recombination sites to the lowest median resolution demonstrated to date (< 1,000 base pairs). We use family inheritance state analysis to control sequencing error and inform family-wide haplotype phasing, allowing quantification of genome-wide compound heterozygosity. We develop a sequence-based methodology for Human Leukocyte Antigen typing that contributes to disease risk prediction. Finally, we advance methods for analysis of disease and pharmacogenomic risk across the coding and non-coding genome that incorporate phased variant data. We show these methods are capable of identifying multigenic risk for inherited thrombophilia and informing the appropriate pharmacological therapy. These ethnicity-specific, family-based approaches to interpretation of genetic variation are emblematic of the next generation of genetic risk assessment using whole-genome sequencing., Competing Interests: JVT and AWZ are founders, consultants, and equity holders in Clinical Future; GMC has advisory roles in and research sponsorships from several companies involved in genome sequencing technology and personal genomics (see http://arep.med.harvard.edu/gmc/tech.html); MS is on the scientific advisory board of DNA Nexus and holds stock in Personalis; RBA has received consultancy fees from Novartis and 23andMe and holds stock in Personalis; AJB is a scientific advisory board member and founder for NuMedii and Genstruct, a scientific advisory board member for Johnson and Johnson, has received consultancy fees from Lilly, NuMedii, Johnson and Johnson, Genstruct, Tercica, and Prevendia and honoraria from Lilly and Siemens, and holds stock in NuMedii, Genstruct, and Personalis. EAA holds stock in Personalis.

Understanding the factors that limit recombination in bacteria is critical in order to better understand and assess its effects on genetic variation and bacterial population genetic structure. Transformation in the naturally competent bacterium, Streptococcus pneumoniae, is regulated by a polymorphic competence (com) apparatus. It has been suggested that polymorphic types, called pherotypes, generate and maintain subpopulation genetic structure within this species. We test predictions stemming from this hypothesis using a cosmopolitan sample of clinical pneumococcal isolates. We sequenced the locus encoding the peptide that induces competence (comC) to assign clones to each known pherotype class and then used multilocus sequence typing to determine whether there is significant genetic differentiation between pherotypes subgroups. We find two dominant pherotypes within our sample, and both are maintained at high frequencies (CSP1 74%, CSP2 26%). Our analyses fail to detect significant genetic differentiation between pherotype groups and find strong evidence, from a coalescent analysis, for interpherotype recombination. In addition, our analyses indicate that positive selection may account for the maintenance of the fixed polymorphism in this locus (comC). Altogether, these results fail to support the prediction that the polymorphism in the competence system acts to limit recombination within S. pneumoniae populations. We discuss why this result is expected given the mechanism underlying transformation and outline a scenario to explain the evolution of polymorphism in the competence system.

The origin of Plasmodium falciparum, the etiological agent of the most dangerous forms of human malaria, remains controversial. Although investigations of homologous parasites in African Apes are crucial to resolve this issue, studies have been restricted to a chimpanzee parasite related to P. falciparum, P. reichenowi, for which a single isolate was available until very recently. Using PCR amplification, we detected Plasmodium parasites in blood samples from 18 of 91 individuals of the genus Pan, including six chimpanzees (three Pan troglodytes troglodytes, three Pan t. schweinfurthii) and twelve bonobos (Pan paniscus). We obtained sequences of the parasites' mitochondrial genomes and/or from two nuclear genes from 14 samples. In addition to P. reichenowi, three other hitherto unknown lineages were found in the chimpanzees. One is related to P. vivax and two to P. falciparum that are likely to belong to distinct species. In the bonobos we found P. falciparum parasites whose mitochondrial genomes indicated that they were distinct from those present in humans, and another parasite lineage related to P. malariae. Phylogenetic analyses based on this diverse set of Plasmodium parasites in African Apes shed new light on the evolutionary history of P. falciparum. The data suggested that P. falciparum did not originate from P. reichenowi of chimpanzees (Pan troglodytes), but rather evolved in bonobos (Pan paniscus), from which it subsequently colonized humans by a host-switch. Finally, our data and that of others indicated that chimpanzees and bonobos maintain malaria parasites, to which humans are susceptible, a factor of some relevance to the renewed efforts to eradicate malaria.

In bacteria, recombination is a rare event, not a part of the reproductive process. Nevertheless, recombination -- broadly defined to include the acquisition of genes from external sources, i.e., horizontal gene transfer (HGT) -- plays a central role as a source of variation for adaptive evolution in many species of bacteria. Much of niche expansion, resistance to antibiotics and other environmental stresses, virulence, and other characteristics that make bacteria interesting and problematic, is achieved through the expression of genes and genetic elements obtained from other populations of bacteria of the same and different species, as well as from eukaryotes and archaea. While recombination of homologous genes among members of the same species has played a central role in the development of the genetics and molecular biology of bacteria, the contribution of homologous gene recombination (HGR) to bacterial evolution is not at all clear. Also, not so clear are the selective pressures responsible for the evolution and maintenance of transformation, the only bacteria-encoded form of HGR. Using a semi-stochastic simulation of mutation, recombination, and selection within bacterial populations and competition between populations, we explore (1) the contribution of HGR to the rate of adaptive evolution in these populations and (2) the conditions under which HGR will provide a bacterial population a selective advantage over non-recombining or more slowly recombining populations. The results of our simulation indicate that, under broad conditions: (1) HGR occurring at rates in the range anticipated for bacteria like Streptococcus pneumoniae, Escherichia coli, Haemophilus influenzae, and Bacillus subtilis will accelerate the rate at which a population adapts to environmental conditions; (2) once established in a population, selection for this capacity to increase rates of adaptive evolution can maintain bacteria-encoded mechanisms of recombination and prevent invasion of non-recombining populations, even when recombination engenders a modest fitness cost; and (3) because of the density- and frequency-dependent nature of HGR in bacteria, this capacity to increase rates of adaptive evolution is not sufficient as a selective force to provide a recombining population a selective advantage when it is rare. Under realistic conditions, homologous gene recombination will increase the rate of adaptive evolution in bacterial populations and, once established, selection for higher rates of evolution will promote the maintenance of bacteria-encoded mechanisms for HGR. On the other hand, increasing rates of adaptive evolution by HGR is unlikely to be the sole or even a dominant selective pressure responsible for the original evolution of transformation., Competing Interests: The authors have declared that no competing interests exist.

The human ZC3H14 gene encodes an evolutionarily conserved Cys(3)His zinc finger protein that binds specifically to polyadenosine RNA and is thus postulated to modulate post-transcriptional gene expression. Expressed sequence tag (EST) data predicts multiple splice variants of both human and mouse ZC3H14. Analysis of ZC3H14 expression in both human cell lines and mouse tissues confirms the presence of multiple alternatively spliced transcripts. Although all of these transcripts encode protein isoforms that contain the conserved C-terminal zinc finger domain, suggesting that they could all bind to polyadenosine RNA, they differ in other functionally important domains. Most of the alternative transcripts encode closely related proteins (termed isoforms 1, 2, 3, and 3 short) that differ primarily in the inclusion of three small exons, 9, 10, and 11, resulting in predicted protein isoforms ranging from 82 to 64 kDa. Each of these closely related isoforms contains predicted classical nuclear localization signals (cNLS) within exons 7 and 11. Consistent with the presence of these putative nuclear targeting signals, these ZC3H14 isoforms are all localized to the nucleus. In contrast, an additional transcript encodes a smaller protein (34 kDa) with an alternative first exon (isoform 4). Consistent with the absence of the predicted cNLS motifs located in exons 7 and 11, ZC3H14 isoform 4 is localized to the cytoplasm. Both EST data and experimental data suggest that this variant is enriched in testes and brain. Using an antibody that detects endogenous ZC3H14 isoforms 1-3 reveals localization of these isoforms to nuclear speckles. These speckles co-localize with the splicing factor, SC35, suggesting a role for nuclear ZC3H14 in mRNA processing. Taken together, these results demonstrate that multiple transcripts encoding several ZC3H14 isoforms exist in vivo. Both nuclear and cytoplasmic ZC3H14 isoforms could have distinct effects on gene expression mediated by the common Cys(3)His zinc finger polyadenosine RNA binding domain.

Agents that kill or induce suicide in the organisms that produce them or other individuals of the same genotype are intriguing puzzles for ecologists and evolutionary biologists. When those organisms are pathogenic bacteria, these suicidal toxins have the added appeal as candidates for the development of narrow spectrum antibiotics to kill the pathogens that produce them.We show that when clinical as well as laboratory strains of Streptococcus pneumoniae are maintained in continuous culture (chemostats), their densities oscillate by as much as five orders of magnitude with an apparently constant period. This dynamic, which is unanticipated for single clones of bacteria in chemostats, can be attributed to population-wide die-offs and recoveries. Using a combination of mathematical models and experiments with S. pneumoniae, we present evidence that these die-offs can be attributed to the autocatalytic production of a toxin that lyses or induces autolysis in members of the clone that produces it. This toxin, which our evidence indicates is a protein, appears to be novel; S. pneumoniae genetic constructs knocked out for lytA and other genes coding for known candidates for this agent oscillate in chemostat culture. Since this toxin lyses different strains of S. pneumoniae as well as other closely related species of Streptococcus, we propose that its ecological role is as an allelopathic agent. Using a mathematical model, we explore the conditions under which toxins that kill members of the same clone that produces them can prevent established populations from invasion by different strains of the same or other species. We postulate that the production of the toxin observed here as well as other bacteria-produced toxins that kill members of the same genotype, 'clonal suicide', evolved and are maintained to prevent colonization of established populations by different strains of the same and closely related species.

The evolutionary history of Plasmodium vivax has recently been addressed in terms of its origin as a parasite of humans and the age of extant populations. The consensus is that P. vivax originated as a result of a host switch from a non-human primate to hominids and that the extant populations did not originate as recently as previously proposed. Here, we show that, in a comparison of parasite isolates from across the world, Asian populations of P. vivax are the oldest. We discuss how this result, together with the phylogenetic evidence that P. vivax derived from Plasmodium found in Southeast Asian macaques, is most simply explained by assuming an Asian origin of this parasite. Nevertheless, the available data show only the tip of the iceberg. We discuss how sampling might affect time estimates to the most recent common ancestor for P. vivax populations and suggest that spatially explicit estimates are needed to understand the demographic history of this parasite better.

The high prevalence of Duffy negativity (lack of the Duffy blood group antigen) among human populations in sub-Saharan Africa has been used to argue that Plasmodium vivax originated on that continent. Here, we investigate the phylogenetic relationships among 10 species of Plasmodium that infect primates by using three genes, two nuclear (beta-tubulin and cell division cycle 2) and a gene from the plastid genome (the elongation factor Tu). We find compelling evidence that P. vivax is derived from a species that inhabited macaques in Southeast Asia. Specifically, those phylogenies that include P. vivax as an ancient lineage from which all of the macaque parasites could originate are significantly less likely to explain the data. We estimate the time to the most recent common ancestor at four neutral gene loci from Asian and South American isolates (a minimum sample of seven isolates per locus). Our analysis estimates that the extant populations of P. vivax originated between 45,680 and 81,607 years ago. The phylogeny and the estimated time frame for the origination of current P. vivax populations are consistent with an "out of Asia" origin for P. vivax as hominoid parasite. The current debate regarding how the Duffy negative trait became fixed in Africa needs to be revisited, taking into account not only human genetic data but also the genetic diversity observed in the extant P. vivax populations and the phylogeny of the genus Plasmodium.

There are few concepts that have been used across disciplines; one of them is natural selection. The impact that this process has on parasite genetic diversity is reviewed here by discussing examples on drug resistance and vaccine antigens. Emphasis is made on how mechanisms need to be addressed rather than associations, and how such investigations were out of reach of biomedical researchers only a decade ago.