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

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 (T_MRCA) 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 T_MRCA 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. 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.

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.

Background: The genusElaeishas two species of economic importance for the oil palm agroindustry:Elaeis oleifera(O), native to the Americas, andElaeis 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.

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.

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