Strong positive selection biases identity-by-descent-based inferences of recent demography and population structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD), yet strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we use simulations, a true IBD inference algorithm, and empirical data sets from different malaria transmission settings to investigate the extent of this bias and explore potential correction strategies. We analyze whole genome sequence data generated from 640 new and 3089 publicly available Plasmodium falciparum clinical isolates. We demonstrate that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discover that the removal of IBD peak regions partially restores the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness and extent of inbreeding. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Survival disparities in non-Hispanic Black and White cervical cancer patients vary by histology and are largely explained by modifiable factors

PURPOSE

We investigated racial disparities in survival by histology in cervical cancer and examined the factors contributing to these disparities.

METHODS

Non-Hispanic Black and non-Hispanic White (hereafter known as Black and White) patients with stage I-IV cervical carcinoma diagnosed between 2004 and 2017 in the National Cancer Database were studied. Survival differences were compared using Cox modeling to estimate hazard ratio (HR) or adjusted HR (AHR) and 95% confidence interval (CI). The contribution of demographic, socioeconomic and clinical factors to the Black vs White differences in survival was estimated after applying propensity score weighting in patients with squamous cell carcinoma (SCC) or adenocarcinoma (AC).

RESULTS

This study included 10,111 Black and 43,252 White patients with cervical cancer. Black patients had worse survival than White cervical cancer patients (HR = 1.40, 95% CI = 1.35-1.45). Survival disparities between Black and White patients varied significantly by histology (HR = 1.20, 95% CI = 1.15-1.24 for SCC; HR = 2.32, 95% CI = 2.12-2.54 for AC, interaction p < 0.0001). After balancing the selected demographic, socioeconomic and clinical factors, survival in Black vs. White patients was no longer different in those with SCC (AHR = 1.01, 95% CI 0.97-1.06) or AC (AHR = 1.09, 95% CI = 0.96-1.24). In SCC, the largest contributors to survival disparities were neighborhood income and insurance. In AC, age was the most significant contributor followed by neighborhood income, insurance, and stage. Diagnosis of AC (but not SCC) at ≥65 years old was more common in Black vs. White patients (26% vs. 13%, respectively).

CONCLUSIONS

Histology matters in survival disparities and diagnosis at ≥65 years old between Black and White cervical cancer patients. These disparities were largely explained by modifiable factors.

Racial, ethnic and country of origin disparities in aggressive endometrial cancer histologic subtypes

OBJECTIVE

This study investigated the risk of an aggressive endometrial cancer (EC) diagnosis by race, ethnicity, and country of origin to further elucidate histologic disparities in non-Hispanic Black (NHB), Hispanic, Asian/Pacific Islander (API), American Indian/Alaskan Native (AIAN) vs. non-Hispanic White (NHW) patients, particularly in Hispanic or API subgroups.

METHODS

Patient diagnosed between 2004 and 2020 with low grade (LG)-endometrioid endometrial cancer (ECC) or an aggressive EC including grade 3 EEC, serous carcinoma, clear cell carcinoma, mixed epithelial carcinoma, or carcinosarcoma in the National Cancer Database were studied. The odds ratio (OR) and 95% confidence interval (CI) for diagnosis of an aggressive EC histology was estimated using logistic modeling.

RESULTS

There were 343,868 NHW, 48,897 NHB, 30,013 Hispanic, 15,015 API and 1646 AIAN patients. The OR (95% CI) for an aggressive EC diagnosis was 3.07 (3.01-3.13) for NHB, 1.08 (1.06-1.11) for Hispanic, 1.17 (1.13-1.21) for API and 1.07 (0.96-1.19) for AIAN, relative to NHW patients. Subset analyses by country of origin illustrated the diversity in the OR for an aggressive EC diagnosis among Hispanic (1.18 for Mexican to 1.87 for Dominican), Asian (1.14 Asian Indian-Pakistani to 1.48 Korean) and Pacific Islander (1.00 for Hawaiian to 1.33 for Samoan) descendants. Hispanic, API and AIAN patients were diagnosed 5-years younger that NHW patients, and the risk for an aggressive EC histology were all significantly higher than NHW patients after correcting for age. Insurance status was another independent risk factor for aggressive histology.

CONCLUSIONS

Risk of an aggressive EC diagnosis varied by race, ethnicity, and country of origin. NHB patients had the highest risk, followed by Dominican, South/Central American, Cuban, Korean, Thai, Vietnamese, and Filipino descendants.

Genome-wide association studies and fine-mapping identify genomic loci for n-3 and n-6 polyunsaturated fatty acids in Hispanic American and African American cohorts

Omega-3 (n-3) and omega-6 (n-6) polyunsaturated fatty acids (PUFAs) play critical roles in human health. Prior genome-wide association studies (GWAS) of n-3 and n-6 PUFAs in European Americans from the CHARGE Consortium have documented strong genetic signals in/near the FADS locus on chromosome 11. We performed a GWAS of four n-3 and four n-6 PUFAs in Hispanic American (n = 1454) and African American (n = 2278) participants from three CHARGE cohorts. Applying a genome-wide significance threshold of P < 5 × 10-8, we confirmed association of the FADS signal and found evidence of two additional signals (in DAGLA and BEST1) within 200 kb of the originally reported FADS signal. Outside of the FADS region, we identified novel signals for arachidonic acid (AA) in Hispanic Americans located in/near genes including TMX2, SLC29A2, ANKRD13D and POLD4, and spanning a > 9 Mb region on chromosome 11 (57.5 Mb ~ 67.1 Mb). Among these novel signals, we found associations unique to Hispanic Americans, including rs28364240, a POLD4 missense variant for AA that is common in CHARGE Hispanic Americans but absent in other race/ancestry groups. Our study sheds light on the genetics of PUFAs and the value of investigating complex trait genetics across diverse ancestry populations.

Strong Positive Selection Biases Identity-By-Descent-Based Inferences of Recent Demography and Population Structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD). Yet, strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we utilized simulations, a true IBD inference algorithm, and empirical datasets from different malaria transmission settings to investigate the extent of such bias and explore potential correction strategies. We analyzed whole genome sequence data generated from 640 new and 4,026 publicly available Plasmodium falciparum clinical isolates. Our findings demonstrated that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discovered that the removal of IBD peak regions partially restored the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Contribution of Common Genetic Variants to Risk of Early-Onset Ischemic Stroke

BACKGROUND AND OBJECTIVES

Current genome-wide association studies of ischemic stroke have focused primarily on late-onset disease. As a complement to these studies, we sought to identify the contribution of common genetic variants to risk of early-onset ischemic stroke.

METHODS

We performed a meta-analysis of genome-wide association studies of early-onset stroke (EOS), ages 18-59 years, using individual-level data or summary statistics in 16,730 cases and 599,237 nonstroke controls obtained across 48 different studies. We further compared effect sizes at associated loci between EOS and late-onset stroke (LOS) and compared polygenic risk scores (PRS) for venous thromboembolism (VTE) between EOS and LOS.

RESULTS

We observed genome-wide significant associations of EOS with 2 variants in ABO, a known stroke locus. These variants tag blood subgroups O1 and A1, and the effect sizes of both variants were significantly larger in EOS compared with LOS. The odds ratio (OR) for rs529565, tagging O1, was 0.88 (95% confidence interval [CI]: 0.85-0.91) in EOS vs 0.96 (95% CI: 0.92-1.00) in LOS, and the OR for rs635634, tagging A1, was 1.16 (1.11-1.21) for EOS vs 1.05 (0.99-1.11) in LOS; p-values for interaction = 0.001 and 0.005, respectively. Using PRSs, we observed that greater genetic risk for VTE, another prothrombotic condition, was more strongly associated with EOS compared with LOS (p = 0.008).

DISCUSSION

The ABO locus, genetically predicted blood group A, and higher genetic propensity for venous thrombosis are more strongly associated with EOS than with LOS, supporting a stronger role of prothrombotic factors in EOS.

Polygenic risk prediction and SNCA haplotype analysis in a Latino Parkinson's disease cohort

BACKGROUND

Large-scale Parkinson's disease (PD) genome-wide association studies (GWAS) have, until recently, only been conducted on subjects with European-ancestry. Consequently, polygenic risk scores (PRS) constructed using PD GWAS data are likely to be less predictive when applied to non-European cohorts.

METHODS

Using GWAS data from the largest study to date, we constructed a PD PRS for a Latino PD cohort (1497 subjects from LARGE-PD) and tested it for association with PD status and age at onset. We validated the PRS performance by testing it in an independent Latino cohort (448 subjects) and by repeating the analysis in LARGE-PD with the addition of 440 external Peruvian controls. We also tested SNCA haplotypes for association with PD risk in LARGE-PD and a European-ancestry PD cohort.

RESULTS

The GWAS-significant PD PRS had an area under the receiver-operator curve (AUC) of 0.668 (95% CI: 0.640-0.695) in LARGE-PD. The inclusion of external Peruvian controls mitigated this result, dropping the AUC 0.632 (95% CI: 0.607-0.657). At the SNCA locus, haplotypes differ by ancestry. Ancestry-specific SNCA haplotypes were associated with PD status in both LARGE-PD and the European-ancestry cohort (p-value < 0.05). These haplotypes both include the rs356182 G-allele, but only share 14% of their variants overall.

CONCLUSION

The PD PRS has potential for PD risk prediction in Latinos, but variability caused by admixture patterns and bias in a European-ancestry PD PRS data limits its utility. The inclusion of diverse subjects can help elucidate PD risk loci and improve risk prediction in non-European cohorts.

Population genetics of PDE4B (Phosphodiesterase-4B) in neglected native americans: implications for cancer pharmacogenetics

PDE4B (phosphodiesterase‐4B) has an important role in cancer and in pharmacology of some disorders, such as inflammatory diseases. Remarkably in Native Americans, PDE4B variants are associated with acute lymphoblastic leukemia (ALL) relapse, as this gene modulates sensitivity of glucocorticoids used in ALL chemotherapy. PDE4B allele rs6683977.G, associated with genomic regions of Native American origin in US‐Hispanics (admixed among Native Americans, Europeans, and Africans), increases ALL relapse risk, contributing to an association between Native American ancestry and ALL relapse that disappeared with an extra‐phase of chemotherapy. This result insinuates that indigenous populations along the Americas may have high frequencies of rs6683977.G, but this has never been corroborated. We studied ancestry and PDE4B diversity in 951 healthy individuals from nine Latin American populations. In non‐admixed Native American populations rs6683977.G has frequencies greater than 90%, is in linkage disequilibrium with other ALL relapse associated and regulatory variants in PDE4B‐intron‐7, conforming haplotypes showing their highest worldwide frequencies in Native Americans (>0.82). Our findings inform the discussion on the pertinence of an extra‐phase of chemotherapy in Native American populations, and exemplifies how knowledge generated in US‐Hispanics is relevant for their even more neglected and vulnerable Native American ancestors along the American continent.

Interpreting Clinical Trials With Omega-3 Supplements in the Context of Ancestry and FADS Genetic Variation

Human diets in developed countries such as the US have changed dramatically over the past 75 years, leading to increased obesity, inflammation, and cardiometabolic dysfunction. Evidence over the past decade indicates that the interaction of genetic variation with changes in the intake of 18-carbon essential dietary omega-6 (n-6) and omega-3 (n-3) polyunsaturated fatty acids (PUFA), linoleic acid (LA) and α-linolenic acid (ALA), respectively, has impacted numerous molecular and clinical phenotypes. Interactions are particularly relevant with the FADS1 and FADS2 genes, which encode key fatty acid desaturases in the pathway that converts LA and ALA to their long chain (≥20 carbons), highly unsaturated fatty acid (HUFA) counterparts. These gene by nutrient interactions affect the levels and balance of n-6 and n-3 HUFA that in turn are converted to a wide array of lipids with signaling roles, including eicosanoids, docosanoids, other oxylipins and endocannabinoids. With few exceptions, n-6 HUFA are precursors of pro-inflammatory/pro-thrombotic signaling lipids, and n-3 HUFA are generally anti-inflammatory/anti-thrombotic. We and others have demonstrated that African ancestry populations have much higher frequencies (vs. European-, Asian- or indigenous Americas-ancestry populations) of a FADS “derived” haplotype that is associated with the efficient conversion of high levels of dietary n-6 PUFA to pro-inflammatory n-6 HUFA. By contrast, an “ancestral” haplotype, carrying alleles associated with a limited capacity to synthesize HUFA, which can lead to n-3 HUFA deficiency, is found at high frequency in certain Hispanic populations and is nearly fixed in several indigenous populations from the Americas. Based on these observations, a focused secondary subgroup analysis of the VITAL n-3 HUFA supplementation trial stratifying the data based on self-reported ancestry revealed that African Americans may benefit from n-3 HUFA supplementation, and both ancestry and FADS variability should be factored into future clinical trials design.

STRIDE: a command-line HMM-based identifier and sub-classifier of Plasmodium falciparum RIFIN and STEVOR variant surface antigen families

Background

RIFINs and STEVORs are variant surface antigens expressed by P. falciparum that play roles in severe malaria pathogenesis and immune evasion. These two highly diverse multigene families feature multiple paralogs, making their classification challenging using traditional bioinformatic methods.

Results

STRIDE (STevor and RIfin iDEntifier) is an HMM-based, command-line program that automates the identification and classification of RIFIN and STEVOR protein sequences in the malaria parasite Plasmodium falciparum. STRIDE is more sensitive in detecting RIFINs and STEVORs than available PFAM and TIGRFAM tools and reports RIFIN subtypes and the number of sequences with a FHEYDER amino acid motif, which has been associated with severe malaria pathogenesis.

Conclusions

STRIDE will be beneficial to malaria research groups analyzing genome sequences and transcripts of clinical field isolates, providing insight into parasite biology and virulence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04515-8.

Peptide ancestry informative markers in uterine neoplasms from women of European, African, and Asian ancestry

Characterization of ancestry-linked peptide variants in disease-relevant patient tissues represents a foundational step to connect patient ancestry with disease pathogenesis. Nonsynonymous single-nucleotide polymorphisms encoding missense substitutions within tryptic peptides exhibiting high allele frequencies in European, African, and East Asian populations, termed peptide ancestry informative markers (pAIMs), were prioritized from 1000 genomes. In silico analysis identified that as few as 20 pAIMs can determine ancestry proportions similarly to >260K SNPs (R² = 0.99). Multiplexed proteomic analysis of >100 human endometrial cancer cell lines and uterine leiomyoma tissues combined resulted in the quantitation of 62 pAIMs that correlate with patient race and genotype-confirmed ancestry. Candidates include a D451E substitution in GC vitamin D-binding protein previously associated with altered vitamin D levels in African and European populations. pAIMs will support generalized proteoancestry assessment as well as efforts investigating the impact of ancestry on the human proteome and how this relates to the pathogenesis of uterine neoplasms.

Population sequencing data reveal a compendium of mutational processes in the human germ line

Biological mechanisms underlying human germline mutations remain largely unknown. We statistically decompose variation in the rate and spectra of mutations along the genome using volume-regularized nonnegative matrix factorization. The analysis of a sequencing dataset (TOPMed) reveals nine processes that explain the variation in mutation properties between loci. We provide a biological interpretation for seven of these processes. We associate one process with bulky DNA lesions that are resolved asymmetrically with respect to transcription and replication. Two processes track direction of replication fork and replication timing, respectively. We identify a mutagenic effect of active demethylation primarily acting in regulatory regions and a mutagenic effect of long interspersed nuclear elements. We localize a mutagenic process specific to oocytes from population sequencing data. This process appears transcriptionally asymmetric.

Characterizing the Genetic Architecture of Parkinson's Disease in Latinos

OBJECTIVE

This work was undertaken in order to identify Parkinson's disease (PD) risk variants in a Latino cohort, to describe the overlap in the genetic architecture of PD in Latinos compared to European-ancestry subjects, and to increase the diversity in PD genome-wide association (GWAS) data.

METHODS

We genotyped and imputed 1,497 PD cases and controls recruited from nine clinical sites across South America. We performed a GWAS using logistic mixed models; variants with a p-value <1 × 10^-5 were tested in a replication cohort of 1,234 self-reported Latino PD cases and 439,522 Latino controls from 23andMe, Inc. We also performed an admixture mapping analysis where local ancestry blocks were tested for association with PD status.

RESULTS

One locus, SNCA, achieved genome-wide significance (p-value <5 × 10^-8 ); rs356182 achieved genome-wide significance in both the discovery and the replication cohorts (discovery, G allele: 1.58 OR, 95% CI 1.35-1.86, p-value 2.48 × 10^-8 ; 23andMe, G allele: 1.26 OR, 95% CI 1.16-1.37, p-value 4.55 × 10^-8 ). In our admixture mapping analysis, a locus on chromosome 14, containing the gene STXBP6, achieved significance in a joint test of ancestries and in the Native American single-ancestry test (p-value <5 × 10^-5 ). A second locus on chromosome 6, containing the gene RPS6KA2, achieved significance in the African single-ancestry test (p-value <5 × 10^-5 ).

INTERPRETATION

This study demonstrated the importance of the SNCA locus for the etiology of PD in Latinos. By leveraging the demographic history of our cohort via admixture mapping, we identified two potential PD risk loci that merit further study. ANN NEUROL 2021.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Evolutionary history of modern Samoans

Archaeological studies estimate the initial settlement of Samoa at 2,750 to 2,880 y ago and identify only limited settlement and human modification to the landscape until about 1,000 to 1,500 y ago. At this point, a complex history of migration is thought to have begun with the arrival of people sharing ancestry with Near Oceanic groups (i.e., Austronesian-speaking and Papuan-speaking groups), and was then followed by the arrival of non-Oceanic groups during European colonialism. However, the specifics of this peopling are not entirely clear from the archaeological and anthropological records, and is therefore a focus of continued debate. To shed additional light on the Samoan population history that this peopling reflects, we employ a population genetic approach to analyze 1,197 Samoan high-coverage whole genomes. We identify population splits between the major Samoan islands and detect asymmetrical gene flow to the capital city. We also find an extreme bottleneck until about 1,000 y ago, which is followed by distinct expansions across the islands and subsequent bottlenecks consistent with European colonization. These results provide for an increased understanding of Samoan population history and the dynamics that inform it, and also demonstrate how rapid demographic processes can shape modern genomes.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.

Inferring the "Primordial Immune Complex": Origins of MHC Class I and Antigen Receptors Revealed by Comparative Genomics

Comparative analyses suggest that the MHC was derived from a prevertebrate "primordial immune complex" (PIC). PIC duplicated twice in the well-studied two rounds of genome-wide duplications (2R) early in vertebrate evolution, generating four MHC paralogous regions (predominantly on human chromosomes [chr] 1, 6, 9, 19). Examining chiefly the amphibian Xenopus laevis, but also other vertebrates, we identified their MHC paralogues and mapped MHC class I, AgR, and "framework" genes. Most class I genes mapped to MHC paralogues, but a cluster of Xenopus MHC class Ib genes (xnc), which previously was mapped outside of the MHC paralogues, was surrounded by genes syntenic to mammalian CD1 genes, a region previously proposed as an MHC paralogue on human chr 1. Thus, this gene block is instead the result of a translocation that we call the translocated part of the MHC paralogous region (MHCtrans) Analyses of Xenopus class I genes, as well as MHCtrans, suggest that class I arose at 1R on the chr 6/19 ancestor. Of great interest are nonrearranging AgR-like genes mapping to three MHC paralogues; thus, PIC clearly contained several AgR precursor loci, predating MHC class I/II. However, all rearranging AgR genes were found on paralogues derived from the chr 19 precursor, suggesting that invasion of a variable (V) exon by the RAG transposon occurred after 2R. We propose models for the evolutionary history of MHC/TCR/Ig and speculate on the dichotomy between the jawless (lamprey and hagfish) and jawed vertebrate adaptive immune systems, as we found genes related to variable lymphocyte receptors also map to MHC paralogues.

Evolution of Hominin Polyunsaturated Fatty Acid Metabolism: From Africa to the New World

The metabolic conversion of dietary omega-3 and omega-6 18 carbon (18C) to long chain (>20 carbon) polyunsaturated fatty acids (LC-PUFAs) is vital for human life. The rate-limiting steps of this process are catalyzed by fatty acid desaturase (FADS) 1 and 2. Therefore, understanding the evolutionary history of the FADS genes is essential to our understanding of hominin evolution. The FADS genes have two haplogroups, ancestral and derived, with the derived haplogroup being associated with more efficient LC-PUFA biosynthesis than the ancestral haplogroup. In addition, there is a complex global distribution of these haplogroups that is suggestive of Neanderthal introgression. We confirm that Native American ancestry is nearly fixed for the ancestral haplogroup, and replicate a positive selection signal in Native Americans. This positive selection potentially continued after the founding of the Americas, although simulations suggest that the timing is dependent on the allele frequency of the ancestral Beringian population. We also find that the Neanderthal FADS haplotype is more closely related to the derived haplogroup and the Denisovan clusters closer to the ancestral haplogroup. Furthermore, the derived haplogroup has a time to the most recent common ancestor of 688,474 years before present. These results support an ancient polymorphism, as opposed to Neanderthal introgression, forming in the FADS region during the Pleistocene with possibly differential selection pressures on both haplogroups. The near fixation of the ancestral haplogroup in Native American ancestry calls for future studies to explore the potential health risk of associated low LC-PUFA levels in these populations.

Genomic structure and diversity of Plasmodium falciparum in Southeast Asia reveal recent parasite migration patterns

Estimates of Plasmodium falciparum migration may inform strategies for malaria elimination. Here we elucidate fine-scale parasite population structure and infer recent migration across Southeast Asia using identity-by-descent (IBD) approaches based on genome-wide single nucleotide polymorphisms called in 1722 samples from 54 districts. IBD estimates are consistent with isolation-by-distance. We observe greater sharing of larger IBD segments between artemisinin-resistant parasites versus sensitive parasites, which is consistent with the recent spread of drug resistance. Our IBD analyses reveal actionable patterns, including isolated parasite populations, which may be prioritized for malaria elimination, as well as asymmetrical migration identifying potential sources and sinks of migrating parasites.

Native American Genomic Diversity through Ancient DNA

Ancient DNA is a powerful tool to understand the evolutionary dynamics of both current and ancestral populations. Posth et al. use ancient DNA to elucidate important questions surrounding the peopling of Central and South America, giving us greater insights into the ancestry of genetically understudied populations.

Improving Cancer Detection and Treatment with Liquid Biopsies and ptDNA

Liquid biopsy, or the capacity to noninvasively isolate and analyze plasma tumor DNA (ptDNA) using blood samples, represents an important tool for modern oncology that enables increasingly safe, personalized, and robust cancer diagnosis and treatment. Here, we review advances in the development and implementation of liquid biopsy approaches, and we focus on the capacity of liquid biopsy to noninvasively detect oncological disease and enhance early detection strategies. In addition to noting the distinctions between mutation-targeted and mutation-agnostic approaches, we discuss the potential for genomic analysis and longitudinal testing to identify somatic lesions early and to guide intervention at more manageable disease stages.

Using genotype array data to compare multi- and single-sample variant calls and improve variant call sets from deep coverage whole-genome sequencing data

Motivation

Variant calling from next-generation sequencing (NGS) data is susceptible to false positive calls due to sequencing, mapping and other errors. To better distinguish true from false positive calls, we present a method that uses genotype array data from the sequenced samples, rather than public data such as HapMap or dbSNP, to train an accurate classifier using Random Forests. We demonstrate our method on a set of variant calls obtained from 642 African-ancestry genomes from the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA), sequenced to high depth (30X).

Results

We have applied our classifier to compare call sets generated with different calling methods, including both single-sample and multi-sample callers. At a False Positive Rate of 5%, our method determines true positive rates of 97.5%, 95% and 99% on variant calls obtained using Illuminas single-sample caller CASAVA, Real Time Genomics multisample variant caller, and the GATK UnifiedGenotyper, respectively. Since NGS sequencing data may be accompanied by genotype data for the same samples, either collected concurrent to sequencing or from a previous study, our method can be trained on each dataset to provide a more accurate computational validation of site calls compared to generic methods. Moreover, our method allows for adjustment based on allele frequency (e.g. a different set of criteria to determine quality for rare versus common variants) and thereby provides insight into sequencing characteristics that indicate call quality for variants of different frequencies.

Availability and Implementation

Code is available on Github at: https://github.com/suyashss/variant_validation.

Contacts

suyashs@stanford.edu or mtaub@jhsph.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Accurate and equitable medical genomic analysis requires an understanding of demography and its influence on sample size and ratio

In a recent study, Petrovski and Goldstein reported that (non-Finnish) Europeans have significantly fewer nonsynonymous singletons in Online Mendelian Inheritance in Man (OMIM) disease genes compared with Africans, Latinos, South Asians, East Asians, and other unassigned non-Europeans. We use simulations of Exome Aggregation Consortium (ExAC) data to show that sample size and ratio interact to influence the number of these singletons identified in a cohort. These interactions are different across ancestries and can lead to the same number of identified singletons in both Europeans and non-Europeans without an equal number of samples. We conclude that there is a need to account for the ancestry-specific influence of demography on genomic architecture and rare variant analysis in order to address inequalities in medical genomic analysis.The authors of the original article were invited to submit a response, but declined to do so. Please see related Open Letter: http://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-1016-y.

A continuum of admixture in the Western Hemisphere revealed by the African Diaspora genome

The African Diaspora in the Western Hemisphere represents one of the largest forced migrations in history and had a profound impact on genetic diversity in modern populations. To date, the fine-scale population structure of descendants of the African Diaspora remains largely uncharacterized. Here we present genetic variation from deeply sequenced genomes of 642 individuals from North and South American, Caribbean and West African populations, substantially increasing the lexicon of human genomic variation and suggesting much variation remains to be discovered in African-admixed populations in the Americas. We summarize genetic variation in these populations, quantifying the postcolonial sex-biased European gene flow across multiple regions. Moreover, we refine estimates on the burden of deleterious variants carried across populations and how this varies with African ancestry. Our data are an important resource for empowering disease mapping studies in African-admixed individuals and will facilitate gene discovery for diseases disproportionately affecting individuals of African ancestry.

An integrated genomic and transcriptomic survey of mucormycosis-causing fungi

Mucormycosis is a life-threatening infection caused by Mucorales fungi. Here we sequence 30 fungal genomes, and perform transcriptomics with three representative Rhizopus and Mucor strains and with human airway epithelial cells during fungal invasion, to reveal key host and fungal determinants contributing to pathogenesis. Analysis of the host transcriptional response to Mucorales reveals platelet-derived growth factor receptor B (PDGFRB) signaling as part of a core response to divergent pathogenic fungi; inhibition of PDGFRB reduces Mucorales-induced damage to host cells. The unique presence of CotH invasins in all invasive Mucorales, and the correlation between CotH gene copy number and clinical prevalence, are consistent with an important role for these proteins in mucormycosis pathogenesis. Our work provides insight into the evolution of this medically and economically important group of fungi, and identifies several molecular pathways that might be exploited as potential therapeutic targets.

AdmixKJump: identifying population structure in recently diverged groups

Motivation

Correctly modeling population structure is important for understanding recent evolution and for association studies in humans. While pre-existing knowledge of population history can be used to specify expected levels of subdivision, objective metrics to detect population structure are important and may even be preferable for identifying groups in some situations. One such metric for genomic scale data is implemented in the cross-validation procedure of the program ADMIXTURE, but it has not been evaluated on recently diverged and potentially cryptic levels of population structure. Here, I develop a new method, AdmixKJump, and test both metrics under this scenario.

Findings

I show that AdmixKJump is more sensitive to recent population divisions compared to the cross-validation metric using both realistic simulations, as well as 1000 Genomes Project European genomic data. With two populations of 50 individuals each, AdmixKJump is able to detect two populations with 100% accuracy that split at least 10KYA, whereas cross-validation obtains this 100% level at 14KYA. I also show that AdmixKJump is more accurate with fewer samples per population. Furthermore, in contrast to the cross-validation approach, AdmixKJump is able to detect the population split between the Finnish and Tuscan populations of the 1000 Genomes Project.

Conclusion

AdmixKJump has more power to detect the number of populations in a cohort of samples with smaller sample sizes and shorter divergence times.

Availability

A java implementation can be found at https://sites.google.com/site/igsevolgenomicslab/home/downloads

Rare variation facilitates inferences of fine-scale population structure in humans

Understanding the genetic structure of human populations has important implications for the design and interpretation of disease mapping studies and reconstructing human evolutionary history. To date, inferences of human population structure have primarily been made with common variants. However, recent large-scale resequencing studies have shown an abundance of rare variation in humans, which may be particularly useful for making inferences of fine-scale population structure. To this end, we used an information theory framework and extensive coalescent simulations to rigorously quantify the informativeness of rare and common variation to detect signatures of fine-scale population structure. We show that rare variation affords unique insights into patterns of recent population structure. Furthermore, to empirically assess our theoretical findings, we analyzed high-coverage exome sequences in 6,515 European and African American individuals. As predicted, rare variants are more informative than common polymorphisms in revealing a distinct cluster of European–American individuals, and subsequent analyses demonstrate that these individuals are likely of Ashkenazi Jewish ancestry. Our results provide new insights into the population structure using rare variation, which will be an important factor to account for in rare variant association studies.

Genetic architecture of quantitative traits and complex diseases

More than 150 years after Mendel discovered the laws of heredity, the genetic architecture of phenotypic variation remains elusive. Here, we discuss recent progress in deciphering how genotypes map onto phenotypes, sources of genetic complexity, and how model organisms are illuminating general principles about the relationship between genetic and phenotypic variation. Moreover, we highlight insights gleaned from large-scale sequencing studies in humans, and how this knowledge informs outstanding questions about the genetic architecture of quantitative traits and complex diseases. Finally, we articulate how the confluence of technologies enabling whole-genome sequencing, comprehensive phenotyping, and high-throughput functional assays of polymorphisms will facilitate a more principled and mechanistic understanding of the genetic architecture of phenotypic variation.

Evolution and functional impact of rare coding variation from deep sequencing of human exomes

As a first step toward understanding how rare variants contribute to risk for complex diseases, we sequenced 15,585 human protein-coding genes to an average median depth of 111× in 2440 individuals of European (n = 1351) and African (n = 1088) ancestry. We identified over 500,000 single-nucleotide variants (SNVs), the majority of which were rare (86% with a minor allele frequency less than 0.5%), previously unknown (82%), and population-specific (82%). On average, 2.3% of the 13,595 SNVs each person carried were predicted to affect protein function of ~313 genes per genome, and ~95.7% of SNVs predicted to be functionally important were rare. This excess of rare functional variants is due to the combined effects of explosive, recent accelerated population growth and weak purifying selection. Furthermore, we show that large sample sizes will be required to associate rare variants with complex traits.

Insights into hominid evolution from the gorilla genome sequence

Gorillas are humans’ closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago (Mya). In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.

The promise and limitations of population exomics for human evolution studies

Exome sequencing is poised to yield substantial insights into human genetic variation and evolutionary history, but there are significant challenges to overcome before this becomes a reality.

Genotype-phenotype associations: substitution models to detect evolutionary associations between phenotypic variables and genotypic evolutionary rate

Motivation: Mapping between genotype and phenotype is one of the primary goals of evolutionary genetics but one that has received little attention at the interspecies level. Recent developments in phylogenetics and statistical modelling have typically been used to examine molecular and phenotypic evolution separately. We have used this background to develop phylogenetic substitution models to test for associations between evolutionary rate of genotype and phenotype. We do this by creating hybrid rate matrices between genotype and phenotype.

Results: Simulation results show our models to be accurate in detecting genotype–phenotype associations and robust for various factors that typically affect maximum likelihood methods, such as number of taxa, level of relevant signal, proportion of sites affected and length of evolutionary divergence. Further, simulations show that our method is robust to homogeneity assumptions. We apply the models to datasets of male reproductive system genes in relation to mating systems of primates. We show that evolution of semenogelin II is significantly associated with mating systems whereas two negative control genes (cytochrome b and peptidase inhibitor 3) show no significant association. This provides the first hybrid substitution model of which we are aware to directly test the association between genotype and phenotype using a phylogenetic framework.

Availability: Perl and HYPHY scripts are available upon request from the authors.

Contact:to252@cam.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Pharmacogenomics: analysing SNPs in the CYP2D6 gene using amino acid properties

The CYP2D6 gene is responsible for metabolising a large portion of the commonly prescribed drugs. Because of its importance, various approaches have been taken to analyse CYP2D6 and Single Nucleotide Polymorphisms (SNPs) throughout its sequence. This study introduces a novel method to analyse the effects of SNPs on encoded protein complexes by focusing on the biochemical properties of each non-synonymous substitution using the program TreeSAAP. Our results show four SNPs in CYP2D6 that exhibit radical changes in amino acid properties which may cause a lack of functionality in the CYP2D6 gene and contribute to a person's inability to metabolise specific drugs.

Effects of position on respiratory muscle function during CO2 rebreathing

The effects of changing from the sitting to supine position on respiratory muscle function was assessed during CO2 rebreathing. Gastric (Pg), pleural (Ppl) and transdiaphragmatic (Pdi) pressures and thoracoabdominal motion were monitored. Diaphragmatic EMG was measured by a bipolar esophageal electrode and quantitated as a moving time average (EMGdi). From sitting to supine, in only 2 of 7 subjects (group A) the diaphragm gained a mechanical advantage as evident by an increased slope of the Pdi versus EMGdi relationship not present in the other 5 subjects (group B). At high levels of ventilation while sitting, only group B increased expiratory abdominal muscle activity leading to a more favorable diaphragm length and a passive descent of the abdomen-diaphragm on inspiration. In the supine position functional residual capacity progressively increased in all subjects and the above abdominal pattern was not seen. We conclude that during upright CO2 rebreathing the recruitment of the expiratory abdominal muscles assists diaphragmatic function by placing the diaphragm in an advantageous pressure generating configuration.

Diaphragmatic EMG and transdiaphragmatic pressure measurements with a single catheter

A single gastroesophageal catheter was used for simultaneous measurements of diaphragmatic electromyogram (EMGdi), esophageal, and gastric pressures in 10 normal volunteers. The catheter consisted of 2 polyethylene tubes, each with an outer diameter of 1.70 mm and an inner diameter of 1.19 mm, 2 platinum wire coils, and esophageal and gastric latex balloons. In all subjects studied, the increase in EMGdi, quantified as the average rate of rise of inspiratory moving average activity, and transdiaphragmatic pressure (Pdi) were linearly related to the increase in end-tidal PCO2 during CO2 rebreathing (range of r, 0.88 to 0.99 and 0.85 to 0.99, respectively). The relation between changes in EMGdi and that of Pdi, mean inspiratory flow, and occlusion pressure were also linear (range of r, 0.83-0.99, 0.73-0.99, and 0.84-0.99, respectively), indicating reliable recordings of EMGdi and Pdi during CO2 rebreathing in upright normal humans.

Diaphragmatic and genioglossal electromyogram responses to isocapnic hypoxia in humans

In order to define the relationship between central control of upper airway and respiratory muscle function, diaphragmatic electromyogram (EMGdi) and genioglossal EMG (EMGge) responses to isocapnic hypoxia were studied in 6 awake supine volunteers. Both EMGs were processed and quantitated as moving time average activity. In all subjects, EMGge showed phasic inspiratory activity synchronous with EMGdi. Increases seen in EMGdi and EMGge were linearly related to the decrease in oxygen saturation (r = 0.89 +/- 0.08 and 0.89 +/- 0.08, respectively). There was also a linear relationship between the relative responses of both EMGs to hypoxia such that a low EMGdi response was associated with a low EMGge response and vice versa (r = 0.92, p less than 0.001). These results indicated that the genioglossus muscle behaves like a respiratory muscle and suggested that central control of upper airway and respiratory muscles in humans are intimately related.

Diaphragmatic and genioglossal electromyogram responses to CO2 rebreathing in humans

To assess the relationship between central control of upper airway and respiratory muscle, simultaneously recorded diaphragmatic electromyogram (EMGdi) and genioglossal EMG (EMG ge) responses to CO2 rebreathing were compared in five supine volunteers. Both EMGs were quantitated in terms of inspiratory peak moving time-average activity. In all subjects both EMGdi and EMGge increased linearly with increasing alveolar CO2 pressure (r = 0.93 +/- 0.04 and 0.87 +/- 0.07, respectively), resulting in a significantly linear EMGge vs. EMGdi relationship (r = 0.91 +/- 0.04). CO2 response slopes of both EMGs were similar and linearly related (r = 0.96, P less than 0.001) such that subjects with low EMGdi response also had a low EMGge response and vice versa. Although the onset of EMGge activity preceded that of EMGdi, and the pattern of both EMGs were different, inspiration duration of both EMGs were similar. These data indicate that in humans both diaphragm and genioglossus muscle share similar control mechanisms and suggest that upper airway function is intimately related to the regulation of breathing.