Search

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 individuals1. 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.

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.

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.

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.

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.

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.

Anthurium andraeanum(Hort.) is an important ornamental in the tropical cut-flower industry. However, there is currently not enough 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 ofA. andraeanumcut-flowers were characterized by transcriptome sequencing technology. For thede novotranscriptome, 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 varieties that differ in spathe color (‘Sasha’ (white), ‘Honduras’ (red), and ‘Rapido’ (purple)). A large number of differentially expressed genes (8,324) that were potentially involved in multiple biological and metabolic pathways were identified, including the flavonoid and anthocyanin biosynthetic pathways. Our results showed that chalcone synthase (CHS) and flavonoid 3’-hydroxylase (F3’H) were the main genes differentially expressed in the white/red/purple comparison. We also identified a differentially expressed cytochromeP450in the late developmental stage of the purple spathe that appeared to determine the difference between the red- and purple-colored spathes. Additionally, putative MYB-domain protein candidates that could be responsible for the control of the biosynthetic pathway were identified. The results provided basic sequence information for future research on spathe color, which have important implications for breeding strategies in this ornamental.Core ideasRNA-seq was performed on three anthurium varieties.Gene expression was compared for developmental stage and spathe color.Differentially expressed unigenes were identified.Putative MYB-domain protein candidates of the anthocyanin biosynthetic pathway were identified.

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 (

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.

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)., Omar Cornejo et al. report a genomic analysis of 200 cacao plants (Theobroma cacao L.) representing more than 10 genetically distinct populations. They identify metabolic and disease resistance genes as contributing to the domestication of cacao and show that domesticated populations maintain a high proportion of deleterious mutations.

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.

Abstract 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. The work presented herein provides, to our knowledge, the first association mapping study in an interspecific OxG hybrid population of oil palm which presents tolerance to pests and diseases, high oil quality, and acceptable fruit bunch production. Results Using genotyping-by-sequencing (GBS), we identified a total of 3,776 single nucleotide polymorphisms (SNPs) that were used to perform a genome-wide association analysis (GWAS) in 378 OxG hybrids for 10 agronomic traits. Twelve genomic regions 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 candidate genes that mapped on genomic regions involved in plant architecture and yield; however, these potential candidate genes need to be confirmed for future targeted functional analysis. The associated markers may be valuable resources for the development of marker-assisted selection in oil palm breeding. Keywords: Association mapping, Elaeis guineensis, Elaeis oleifera, genotyping-by-sequencing, plant architecture, yield.

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. The breeding program in Colombia relies on interspecific OxG crossing populations with tolerance to pests and diseases, high oil quality, and acceptable fruit bunch production. The identification of loci associated to morphological and yield-related traits and the dissection of their genetic architecture will provide essential insights for oil palm breeding strategies. Results The genotypes of 471 oil palms, including 62 E. oleifera (O), 31 E. guineensis (G) and 378 OxG samples were analyzed in this study. A total of 3,776 single nucleotide polymorphisms (SNP) were detected across the 16 oil palm chromosomes using the genotyping-by-sequencing (GBS) technique. The genetic variation and population structure analyses grouped the samples into two clades according to the parental relatedness. A genome wide association analysis (GWAS) was conducting using the OxG hybrid population, resulting in 12 SNPs significantly associated with ten different morphological and yield-related traits. Conclusions The work presented herein provides to our knowledge the first association mapping study in an interspecific OxG hybrid population of oil palm. We provide new insights on candidate genes involved in tissue development and plant architecture associated to traits such as: rachis length, trunk diameter, bunch number, and bunch weight. The genes identified in our analysis are putative candidates for future targeted functional analysis. They are valuable resources for the development of marker-assisted selection in oil palm breeding. Keywords: Association mapping, Elaeis guineensis, Elaeis oleifera, genotyping-by-sequencing, plant architecture, yield.

Identifying loci that are under selection versus those that are evolving neutrally is a common challenge in evolutionary genetics. Moreover, with the increase in sequence data, genomic studies have begun to incorporate the use of multiple methods to identify candidate loci under selection. Composite methods are usually implemented to transform the data into a multi-dimensional scatter where outliers are identified using a distance metric, the most common being Mahalanobis distance. However, studies have shown that the power of Mahalanobis distance reduces as the number of dimensions increases. Because the number of methods for detecting selection continue to grow, this is an undesirable feature of Mahalanobis distance. Other composite methods such as invariant coordinate selection (ICS) have proven to be a robust method for identifying outliers in multi-dimensional space; though, this method has not been implemented for genomic data. Here we use simulated genomic data to test the performance of ICS in identifying outlier loci from multiple selection scans and compare the results to the performance of Mahalanobis distances. We show that the ICS outperformed Mahalanobis distance in all aspects including false positives, false negatives, and recall. Furthermore, ICS also performed better when identifying loci with weaker selection coefficients. We also introduce a pipeline in a R-Shiny smart wrapper environment that implements the ICS on multiple scans of selection. Importantly, we show that the ICS is a robust method for identifying outliers in multi-dimensional space and recommend its use for studies aimed at identifying loci under selection in the genome.

Background Carbonic anhydrase (CA) catalyzes the hydration of CO2 in the first biochemical step of C4 photosynthesis, and has been considered a potentially rate-limiting step when CO2 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 CO2concentrations. However, the ca1ca2 mutants exhibited reduced rates of photosynthesis at sub-ambient CO2, and accumulated less biomass when grown under sub-ambient CO2 (9.2 Pa). To further clarify the importance of CA for C4 photosynthesis, we assessed gene expression changes in wild-type, ca1 and ca1ca2 mutants in response to changes in pCO2 from 920 to 9.2 Pa. Results Leaf samples from each genotype were collected for RNA-seq analysis at high CO2 and at two time points after the low CO2 transition, in order to identify early and longer-term responses to CO2 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 CO2, 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 CO2 stomatal signaling pathway were differentially expressed in the CA mutants under low CO2. 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 CO2. Our results also support a role for CA in the CO2 stomatal signaling pathway. This study provides insight into the importance of CA for C4 photosynthesis and its role in stomatal signaling.

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.

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 commonlyfound in healthy asymptomatic women of all ages. This raisesthe question if genetic differences among strains mightdistinguish potentially pathogenic from commensal strains. To disentangle the diversity ofG.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 comparativegenomic analysis of 35strains 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 inG.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.

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.

Domestication has had a strong impact on the development of modern societies. We sequenced 200 genomes of the chocolate plantTheobroma cacaoL. to show for the first time that a single population underwent strong domestication approximately 3,600 years (95% CI: 2481 – 10,903 years ago) ago, the Criollo population. 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 ofT. cacao(Criollo) maintain a higher proportion of high frequency deleterious mutations. We also show for the first time the negative consequences the increase accumulation of deleterious mutations during domestication on the fitness of individuals (significant negative correlation between Criollo ancestry and Kg of beans per hectare per year, P = 0.000425).

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.

Aneuploidy is prevalent in human preimplantation embryos and is the leading cause of pregnancy loss. Many aneuploidies arise during oogenesis, increasing in frequency with maternal age. Superimposed on these meiotic aneuploidies are a range of errors occurring during early mitotic divisions of the embryo, 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 to quantify their incidences and gain insight into their mechanisms of formation. In addition to well-described patterns such as triploidy and haploidy, we identified 4.7% of day-3 blastomeres possessing karyotypes suggestive of tripolar mitosis in normally-fertilized diploid zygotes or descendant diploid cells. We further supported this hypothesis using time-lapse data from an intersecting set of 77 cleavage-stage embryos. The diploid tripolar signature was rare among day-5 blastocyst biopsies (0.5%), suggesting that complex aneuploidy generated by tripolar mitosis impairs cellular and/or early embryonic survival. Strikingly, we found that the tripolar mitosis mechanism is responsible for the previously described association with common maternal genetic variants spanning PLK4. Our findings are consistent with the role of PLK4 as the master regulator of centriole duplication with a known capacity to induce tripolar mitosis when mutated or mis-expressed. Taken together, we propose that tripolar mitosis is a key mechanism generating karyotype-wide aneuploidy in cleavage-stage embryos and implicate PLK4-mediated centrosome abnormality as a factor influencing its occurrence.

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., Author Summary Infectious diseases have undoubtedly played an important role in ancient and modern human history. Yet, there are relatively few regions of the genome involved in resistance to pathogens that show a strong selection signal in current genome-wide searches for this kind of signal. We revisit the evolutionary history of a gene associated with resistance to the most common malaria-causing parasite, Plasmodium vivax, and show that it is one of regions of the human genome that has been under strongest selective pressure in our evolutionary history (selection coefficient: 4.3%). Our results are consistent with a complex evolutionary history of the locus involving selection on a mutation that was at a very low frequency in the ancestral African population (standing variation) and subsequent differentiation between European, Asian and African populations.

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.

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 C-14 years before present (BP) (13,000 to 12,600 calendar years BP)(1,2). 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(3). 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(2). An alternative, Solutrean, hypothesis posits that the Clovis predecessors emigrated from southwestern Europe during the Last Glacial Maximum(4). 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 C-14 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.4x and show that the gene flow from the Siberian Upper Palaeolithic Mal'ta population(5) 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.

Natural transformation in the Gram-positive pathogenStreptococcus pneumoniaeoccurs when cells become “competent”, a state that is induced in response to high extracellular concentrations of a secreted peptide signal called CSP (Competence Stimulating Peptide) encoded by thecomClocus. Two main CSP signal types (pherotypes) are known to dominate the pherotype diversity across strains. Using thousands of 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 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 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 marginally facilitate, recombination between strains. Furthermore, our results suggest that the CSP balanced polymorphism does not causally underlie population differentiation. Therefore, when strains carrying different pherotypes encounter one another during co-colonization, genetic exchange can freely occur.

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., Facultad de Ciencias Naturales y Museo, Instituto Multidisciplinario de Biología Celular

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 (IA 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.

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.

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 ofStreptococcus pneumoniaeare 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 withS. 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. pneumoniaegenetic constructs knocked out forlytAand other genes coding for known candidates for this agent oscillate in chemostat culture. Since this toxin lyses different strains ofS. pneumoniaeas well as other closely related species ofStreptococcus, 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.

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 find 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 (KYA), and after no more than 8,000-year isolation period in Beringia. Following their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 KYA, one that is now dispersed across North and South America and the other is 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.

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 (β-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.

Arctic genetics comes in from the coldDespite a well-characterized archaeological record, the genetics of the people who inhabit the Arctic have been unexplored. Raghavanet al.sequenced ancient and modern genomes of individuals from the North American Arctic (see the Perspective by Park). Analyses of these genomes indicate that the Arctic was colonized 6000 years ago by a migration separate from the one that gave rise to other Native American populations. Furthermore, the original paleo-inhabitants of the Arctic appear to have been completely replaced approximately 700 years ago.Science, this issue10.1126/science.1255832; see also p.1004