Efficacy of National Comprehensive Cancer Network Guidelines in Identifying Pathogenic Germline Variants Among Unselected Patients with Prostate Cancer: The PROCLAIM Trial

BACKGROUND

Prostate cancer (PCa) patients with pathogenic/likely pathogenic germline variants (PGVs) in cancer predisposition genes may be eligible for U.S. Food and Drug Administration-approved targeted therapies, clinical trials, or enhanced screening. Studies suggest that eligible patients are missing genetics-informed care due to restrictive testing criteria.

OBJECTIVE

To establish the prevalence of actionable PGVs among prospectively accrued, unselected PCa patients, stratified by their guideline eligibility.

DESIGN, SETTING, AND PARTICIPANTS

Consecutive, unselected PCa patients were enrolled at 15 sites in the USA from October 2019 to August 2021, and had multigene cancer panel testing.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Correlates between the prevalence of PGVs and clinician-reported demographic and clinical characteristics were examined.

RESULTS AND LIMITATIONS

Among 958 patients (median [quartiles] age at diagnosis 65 [60, 71] yr), 627 (65%) had low- or intermediate-risk disease (grade group 1, 2, or 3). A total of 77 PGVs in 17 genes were identified in 74 patients (7.7%, 95% confidence interval [CI] 6.2-9.6%). No significant difference was found in the prevalence of PGVs among patients who met the 2019 National Comprehensive Cancer Network Prostate criteria (8.8%, 43/486, 95% CI 6.6-12%) versus those who did not (6.6%, 31/472, 95% CI 4.6-9.2%; odds ratio 1.38, 95% CI 0.85-2.23), indicating that these criteria would miss 42% of patients (31/74, 95% CI 31-53%) with PGVs. The criteria were less effective at predicting PGVs in patients from under-represented populations. Most PGVs (81%, 60/74) were potentially clinically actionable. Limitations include the inability to stratify analyses based on individual ethnicity due to low numbers of non-White patients with PGVs.

CONCLUSIONS

Our results indicate that almost half of PCa patients with PGVs are missed by current testing guidelines. Comprehensive germline genetic testing should be offered to all patients with PCa.

PATIENT SUMMARY

One in 13 patients with prostate cancer carries an inherited variant that may be actionable for the patient's current care or prevention of future cancer, and could benefit from expanded testing criteria.

The hydrocortisone-responsive urinary metabolome of premature infants

BACKGROUND

Extremely premature infants are at risk for circulatory collapse or respiratory failure that are often treated with hydrocortisone (HC); however, there is no information on the metabolic consequences of this therapy.

METHODS

Longitudinal urine samples from infants <28 weeks gestation in the Trial of Late Surfactant were analyzed by untargeted UHPLC:MS/MS. Fourteen infants who received a tapering course of HC beginning at 3 mg/kg/day for ≥9 days were compared to 14 matched control infants. A secondary cross-sectional analysis by logistic regression used urines from 314 infants.

RESULTS

Of 1145 urinary metabolites detected, abundance of 219, representing all the major biochemical pathways, changed at p < 0.05 in the HC-treated group with 90% decreasing; 3 cortisol derivatives increased ~2-fold with HC therapy. Only 11% of regulated metabolites remained responsive at the lowest HC dose. Regulated metabolites included two steroids and thiamin that are associated with lung inflammation in infants. HC responsiveness was confirmed in 57% of metabolites by cross-sectional analysis.

CONCLUSIONS

HC treatment of premature infants influenced in a dose-dependent manner abundance of 19% of identified urinary metabolites of diverse biochemical systems, primarily reducing concentrations. These findings indicate that exposure to HC reversibly impacts the nutritional status of premature infants.

IMPACT

Hydrocortisone treatment of premature infants with respiratory failure or circulatory collapse alters levels of a subset of urinary metabolites representing all major biochemical pathways. This is the first description of the scope, magnitude, timing and reversibility of metabolomic changes in infants in response to hydrocortisone, and it confirms corticosteroid regulation of three biochemicals that are associated with lung inflammatory status. The findings indicate a dose-dependency of hydrocortisone for metabolomic and anti-inflammatory effects, that prolonged therapy may lower the supply of many nutrients, and that monitoring concentrations of cortisol and inflammation markers may be a useful clinical approach during corticosteroid therapy.

Effect of different ankle joint positions on medial gastrocnemius muscle fiber strains during isometric plantarflexion

Muscle force production is influenced by muscle fiber and aponeurosis architecture. This prospective cohort study utilizes special MR imaging sequences to examine the structure-function in-vivo in the Medial Gastrocnemius (MG) at three-ankle angles (dorsiflexion, plantar flexion-low and high) and two sub-maximal levels of maximum voluntary contraction (25% and 50%MVC). The study was performed on 6 young male participants. Muscle fiber and aponeurosis strain, fiber strain normalized to force, fiber length and pennation angle (at rest and peak contraction) were analyzed for statistical differences between ankle positions and %MVC. A two-way repeated measures ANOVA and post hoc Bonferroni-adjusted tests were conducted for normal data. A related samples test with Friedman's 2-way ANOVA by ranks with corrections for multiple comparisons was conducted for non-normal data. The dorsiflexed ankle position generated significantly higher force with lower fiber strain than the plantarflexed positions. Sarcomere length extracted from muscle fiber length at each ankle angle was used to track the location on the Force-Length curve and showed the MG operates on the curve's ascending limb. Muscle force changes predicted from the F-L curve going from dorsi- to plantarflexion was less than that experimentally observed suggesting other determinants of force changes with ankle position.

The Urinary Metabolomic Fingerprint in Extremely Preterm Infants on Total Parenteral Nutrition vs. Enteral Feeds

Total Parenteral Nutrition (TPN), which uses intravenous administration of nutrients, minerals and vitamins, is essential for sustaining premature infants until they transition to enteral feeds, but there is limited information on metabolomic differences between infants on TPN and enteral feeds. We performed untargeted global metabolomics on urine samples collected between 23-30 days of life from 314 infants born <29 weeks gestational age from the TOLSURF and PROP cohorts. Principal component analysis across all metabolites showed a separation of infants solely on TPN compared to infants who had transitioned to enteral feeds, indicating global metabolomic differences between infants based on feeding status. Among 913 metabolites that passed quality control filters, 609 varied in abundance between infants on TPN vs. enteral feeds at p < 0.05. Of these, 88% were in the direction of higher abundance in the urine of infants on enteral feeds. In a subset of infants in a longitudinal analysis, both concurrent and delayed changes in metabolite levels were observed with the initiation of enteral feeds. These infants had higher concentrations of essential amino acids, lipids, and vitamins, which are necessary for growth and development, suggesting the nutritional benefit of an enteral feeding regimen.

Effect of Ocrelizumab on B- and T-Cell Receptor Repertoire Diversity in Patients With Relapsing Multiple Sclerosis From the Randomized Phase III OPERA Trial

BACKGROUND AND OBJECTIVES

The B cell-depleting anti-CD20 antibody ocrelizumab (OCR) effectively reduces MS disease activity and slows disability progression. Given the role of B cells as antigen-presenting cells, the primary goal of this study was to evaluate the effect of OCR on the T-cell receptor repertoire diversity.

METHODS

To examine whether OCR substantially alters the molecular diversity of the T-cell receptor repertoire, deep immune repertoire sequencing (RepSeq) of CD4⁺ and CD8⁺ T-cell receptor β-chain variable regions was performed on longitudinal blood samples. The IgM and IgG heavy chain variable region repertoire was also analyzed to characterize the residual B-cell repertoire under OCR treatment.

RESULTS

Peripheral blood samples for RepSeq were obtained from 8 patients with relapsing MS enrolled in the OPERA I trial over a period of up to 39 months. Four patients each were treated with OCR or interferon β1-a during the double-blind period of OPERA I. All patients received OCR during the open-label extension. The diversity of the CD4⁺/CD8⁺ T-cell repertoires remained unaffected in OCR-treated patients. The expected OCR-associated B-cell depletion was mirrored by reduced B-cell receptor diversity in peripheral blood and a shift in immunoglobulin gene usage. Despite deep B-cell depletion, longitudinal persistence of clonally related B-cells was observed.

DISCUSSION

Our data illustrate that the diversity of CD4⁺/CD8⁺ T-cell receptor repertoires remained unaltered in OCR-treated patients with relapsing MS. Persistence of a highly diverse T-cell repertoire suggests that aspects of adaptive immunity remain intact despite extended anti-CD20 therapy.

TRIAL REGISTRATION INFORMATION

This is a substudy (BE29353) of the OPERA I (WA21092; NCT01247324) trial. Date of registration, November 23, 2010; first patient enrollment, August 31, 2011.

Assessing the contribution of rare variants to complex trait heritability from whole-genome sequence data

Analyses of data from genome-wide association studies on unrelated individuals have shown that, for human traits and diseases, approximately one-third to two-thirds of heritability is captured by common SNPs. However, it is not known whether the remaining heritability is due to the imperfect tagging of causal variants by common SNPs, in particular whether the causal variants are rare, or whether it is overestimated due to bias in inference from pedigree data. Here we estimated heritability for height and body mass index (BMI) from whole-genome sequence data on 25,465 unrelated individuals of European ancestry. The estimated heritability was 0.68 (standard error 0.10) for height and 0.30 (standard error 0.10) for body mass index. Low minor allele frequency variants in low linkage disequilibrium (LD) with neighboring variants were enriched for heritability, to a greater extent for protein-altering variants, consistent with negative selection. Our results imply that rare variants, in particular those in regions of low linkage disequilibrium, are a major source of the still missing heritability of complex traits and disease.

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.

A Cutibacterium acnes antibiotic modulates human skin microbiota composition in hair follicles

The composition of the skin microbiota varies widely among individuals when sampled at the same body site. A key question is which molecular factors determine strain-level variability within sub-ecosystems of the skin microbiota. Here, we used a genomics-guided approach to identify an antibacterial biosynthetic gene cluster in Cutibacterium acnes (formerly Propionibacterium acnes), a human skin commensal bacterium that is widely distributed across individuals and skin sites. Experimental characterization of this biosynthetic gene cluster resulted in identification of a new thiopeptide antibiotic, cutimycin. Analysis of individual human skin hair follicles revealed that cutimycin contributed to the ecology of the skin hair follicle microbiota and helped to reduce colonization of skin hair follicles by Staphylococcus species.

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

Recent shifts in the genomic ancestry of Mexican Americans may alter the genetic architecture of biomedical traits

People in the Americas represent a diverse continuum of populations with varying degrees of admixture among African, European, and Amerindigenous ancestries. In the United States, populations with non-European ancestry remain understudied, and thus little is known about the genetic architecture of phenotypic variation in these populations. Using genotype data from the Hispanic Community Health Study/Study of Latinos, we find that Amerindigenous ancestry increased by an average of ~20% spanning 1940s-1990s in Mexican Americans. These patterns result from complex interactions between several population and cultural factors which shaped patterns of genetic variation and influenced the genetic architecture of complex traits in Mexican Americans. We show for height how polygenic risk scores based on summary statistics from a European-based genome-wide association study perform poorly in Mexican Americans. Our findings reveal temporal changes in population structure within Hispanics/Latinos that may influence biomedical traits, demonstrating a need to improve our understanding of admixed populations.

Population genetic simulation study of power in association testing across genetic architectures and study designs

While it is well established that genetics can be a major contributor to population variation of complex traits, the relative contributions of rare and common variants to phenotypic variation remains a matter of considerable debate. Here, we simulate genetic and phenotypic data across different case/control panel sampling strategies, sequencing methods, and genetic architecture models based on evolutionary forces to determine the statistical performance of rare variant association tests (RVATs) widely in use. We find that the highest statistical power of RVATs is achieved by sampling case/control individuals from the extremes of an underlying quantitative trait distribution. We also demonstrate that the use of genotyping arrays, in conjunction with imputation from a whole-genome sequenced (WGS) reference panel, recovers the vast majority (90%) of the power that could be achieved by sequencing the case/control panel using current tools. Finally, we show that for dichotomous traits, the statistical performance of RVATs decreases as rare variants become more important in the trait architecture. Our results extend previous work to show that RVATs are insufficiently powered to make generalizable conclusions about the role of rare variants in dichotomous complex traits.

Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.

Ultrarare variants drive substantial cis heritability of human gene expression

The vast majority of human mutations have minor allele frequencies under 1%, with the plurality observed only once (that is, 'singletons'). While Mendelian diseases are predominantly caused by rare alleles, their cumulative contribution to complex phenotypes is largely unknown. We develop and rigorously validate an approach to jointly estimate the contribution of all alleles, including singletons, to phenotypic variation. We apply our approach to transcriptional regulation, an intermediate between genetic variation and complex disease. Using whole-genome DNA and lymphoblastoid cell line RNA sequencing data from 360 European individuals, we conservatively estimate that singletons contribute approximately 25% of cis heritability across genes (dwarfing the contributions of other frequencies). The majority (approximately 76%) of singleton heritability derives from ultrarare variants absent from thousands of additional samples. We develop an inference procedure to demonstrate that our results are consistent with pervasive purifying selection shaping the regulatory architecture of most human genes.

An admixture mapping meta-analysis implicates genetic variation at 18q21 with asthma susceptibility in Latinos

BACKGROUND

Asthma is a common but complex disease with racial/ethnic differences in prevalence, morbidity, and response to therapies.

OBJECTIVE

We sought to perform an analysis of genetic ancestry to identify new loci that contribute to asthma susceptibility.

METHODS

We leveraged the mixed ancestry of 3902 Latinos and performed an admixture mapping meta-analysis for asthma susceptibility. We replicated associations in an independent study of 3774 Latinos, performed targeted sequencing for fine mapping, and tested for disease correlations with gene expression in the whole blood of more than 500 subjects from 3 racial/ethnic groups.

RESULTS

We identified a genome-wide significant admixture mapping peak at 18q21 in Latinos (P = 6.8 × 10-6), where Native American ancestry was associated with increased risk of asthma (odds ratio [OR], 1.20; 95% CI, 1.07-1.34; P = .002) and European ancestry was associated with protection (OR, 0.86; 95% CI, 0.77-0.96; P = .008). Our findings were replicated in an independent childhood asthma study in Latinos (P = 5.3 × 10-3, combined P = 2.6 × 10-7). Fine mapping of 18q21 in 1978 Latinos identified a significant association with multiple variants 5' of SMAD family member 2 (SMAD2) in Mexicans, whereas a single rare variant in the same window was the top association in Puerto Ricans. Low versus high SMAD2 blood expression was correlated with case status (13.4% lower expression; OR, 3.93; 95% CI, 2.12-7.28; P < .001). In addition, lower expression of SMAD2 was associated with more frequent exacerbations among Puerto Ricans with asthma.

CONCLUSION

Ancestry at 18q21 was significantly associated with asthma in Latinos and implicated multiple ancestry-informative noncoding variants upstream of SMAD2 with asthma susceptibility. Furthermore, decreased SMAD2 expression in blood was strongly associated with increased asthma risk and increased exacerbations.

A genome-wide association and admixture mapping study of bronchodilator drug response in African Americans with asthma

Short-acting β₂-adrenergic receptor agonists (SABAs) are the most commonly prescribed asthma medications worldwide. Response to SABAs is measured as bronchodilator drug response (BDR), which varies among racial/ethnic groups in the U.S^{1, 2}. However, the genetic variation that contributes to BDR is largely undefined in African Americans with asthma³. To identify genetic variants that may contribute to differences in BDR in African Americans with asthma, we performed a genome-wide association study (GWAS) of BDR in 949 African American children with asthma, genotyped with the Axiom World Array 4 (Affymetrix, Santa Clara, CA) followed by imputation using 1000 Genomes phase III genotypes. We used linear regression models adjusting for age, sex, body mass index (BMI) and genetic ancestry to test for an association between BDR and genotype at single nucleotide polymorphisms (SNPs). To increase power and distinguish between shared vs. population-specific associations with BDR in children with asthma, we performed a meta-analysis across 949 African Americans and 1,830 Latinos (Total=2,779). Lastly, we performed genome-wide admixture mapping to identify regions whereby local African or European ancestry is associated with BDR in African Americans. We identified a population-specific association with an intergenic SNP on chromosome 9q21 that was significantly associated with BDR (rs73650726, p=7.69×10⁻⁹). A trans-ethnic meta-analysis across African Americans and Latinos identified three additional SNPs within the intron of PRKG1 that were significantly associated with BDR (rs7903366, rs7070958, and rs7081864, p≤5×10⁻⁸). Our results failed to replicate in three additional populations of 416 Latinos and 1,615 African Americans. Our findings indicate that both population specific and shared genetic variation contributes to differences in BDR in minority children with asthma, and that the genetic underpinnings of BDR may differ between racial/ethnic groups.

Ebola virus, but not Marburg virus, replicates efficiently and without required adaptation in snake cells

Ebola virus (EBOV) disease is a viral hemorrhagic fever with a high case-fatality rate in humans. This disease is caused by four members of the filoviral genus Ebolavirus, including EBOV. The natural hosts reservoirs of ebolaviruses remain to be identified. Glycoprotein 2 of reptarenaviruses, known to infect only boa constrictors and pythons, is similar in sequence and structure to ebolaviral glycoprotein 2, suggesting that EBOV may be able to infect reptilian cells. Therefore, we serially passaged EBOV and a distantly related filovirus, Marburg virus (MARV), in boa constrictor JK cells and characterized viral infection/replication and mutational frequency by confocal imaging and sequencing. We observed that EBOV efficiently infected and replicated in JK cells, but MARV did not. In contrast to most cell lines, EBOV-infected JK cells did not result in an obvious cytopathic effect. Surprisingly, genomic characterization of serial-passaged EBOV in JK cells revealed that genomic adaptation was not required for infection. Deep sequencing coverage (>10,000×) demonstrated the existence of only a single nonsynonymous variant (EBOV glycoprotein precursor pre-GP T544I) of unknown significance within the viral population that exhibited a shift in frequency of at least 10 per cent over six serial passages. In summary, we present the first reptilian cell line that replicates a filovirus at high titers, and for the first time demonstrate a filovirus genus-specific restriction to MARV in a cell line. Our data suggest the possibility that there may be differences between the natural host spectra of ebolaviruses and marburgviruses.

Human demographic history has amplified the effects of background selection across the genome

Natural populations often grow, shrink, and migrate over time. Such demographic processes can affect genome-wide levels of genetic diversity. Additionally, genetic variation in functional regions of the genome can be altered by natural selection, which drives adaptive mutations to higher frequencies or purges deleterious ones. Such selective processes affect not only the sites directly under selection but also nearby neutral variation through genetic linkage via processes referred to as genetic hitchhiking in the context of positive selection and background selection (BGS) in the context of purifying selection. While there is extensive literature examining the consequences of selection at linked sites at demographic equilibrium, less is known about how non-equilibrium demographic processes influence the effects of hitchhiking and BGS. Utilizing a global sample of human whole-genome sequences from the Thousand Genomes Project and extensive simulations, we investigate how non-equilibrium demographic processes magnify and dampen the consequences of selection at linked sites across the human genome. When binning the genome by inferred strength of BGS, we observe that, compared to Africans, non-African populations have experienced larger proportional decreases in neutral genetic diversity in strong BGS regions. We replicate these findings in admixed populations by showing that non-African ancestral components of the genome have also been affected more severely in these regions. We attribute these differences to the strong, sustained/recurrent population bottlenecks that non-Africans experienced as they migrated out of Africa and throughout the globe. Furthermore, we observe a strong correlation between FST and the inferred strength of BGS, suggesting a stronger rate of genetic drift. Forward simulations of human demographic history with a model of BGS support these observations. Our results show that non-equilibrium demography significantly alters the consequences of selection at linked sites and support the need for more work investigating the dynamic process of multiple evolutionary forces operating in concert.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma

RATIONALE

Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response.

OBJECTIVES

To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children.

METHODS

We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR.

MEASUREMENTS AND MAIN RESULTS

We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and β-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings.

CONCLUSIONS

The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

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.

ROP: dumpster diving in RNA-sequencing to find the source of 1 trillion reads across diverse adult human tissues

High-throughput RNA-sequencing (RNA-seq) technologies provide an unprecedented opportunity to explore the individual transcriptome. Unmapped reads are a large and often overlooked output of standard RNA-seq analyses. Here, we present Read Origin Protocol (ROP), a tool for discovering the source of all reads originating from complex RNA molecules. We apply ROP to samples across 2630 individuals from 54 diverse human tissues. Our approach can account for 99.9% of 1 trillion reads of various read length. Additionally, we use ROP to investigate the functional mechanisms underlying connections between the immune system, microbiome, and disease. ROP is freely available at https://github.com/smangul1/rop/wiki.

Cancer-associated arginine-to-histidine mutations confer a gain in pH sensing to mutant proteins

The intracellular pH (pHi) of most cancers is constitutively higher than that of normal cells and enhances proliferation and cell survival. We found that increased pHi enabled the tumorigenic behaviors caused by somatic arginine-to-histidine mutations, which are frequent in cancer and confer pH sensing not seen with wild-type proteins. Experimentally raising the pHi increased the activity of R776H mutant epidermal growth factor receptor (EGFR-R776H), thereby increasing proliferation and causing transformation in fibroblasts. An Arg-to-Gly mutation did not confer these effects. Molecular dynamics simulations of EGFR suggested that decreased protonation of His⁷⁷⁶ at high pH causes conformational changes in the αC helix that may stabilize the active form of the kinase. An Arg-to-His, but not Arg-to-Lys, mutation in the transcription factor p53 (p53-R273H) decreased its transcriptional activity and attenuated the DNA damage response in fibroblasts and breast cancer cells with high pHi. Lowering pHi attenuated the tumorigenic effects of both EGFR-R776H and p53-R273H. Our data suggest that some somatic mutations may confer a fitness advantage to the higher pHi of cancer cells.

Prominent features of the amino acid mutation landscape in cancer

Cancer can be viewed as a set of different diseases with distinctions based on tissue origin, driver mutations, and genetic signatures. Accordingly, each of these distinctions have been used to classify cancer subtypes and to reveal common features. Here, we present a different analysis of cancer based on amino acid mutation signatures. Non-negative Matrix Factorization and principal component analysis of 29 cancers revealed six amino acid mutation signatures, including four signatures that were dominated by either arginine to histidine (Arg>His) or glutamate to lysine (Glu>Lys) mutations. Sample-level analyses reveal that while some cancers are heterogeneous, others are largely dominated by one type of mutation. Using a non-overlapping set of samples from the COSMIC somatic mutation database, we validate five of six mutation signatures, including signatures with prominent arginine to histidine (Arg>His) or glutamate to lysine (Glu>Lys) mutations. This suggests that our classification of cancers based on amino acid mutation patterns may provide avenues of inquiry pertaining to specific protein mutations that may generate novel insights into cancer biology.

Functional Segregation of Overlapping Genes in HIV

Overlapping genes pose an evolutionary dilemma as one DNA sequence evolves under the selection pressures of multiple proteins. Here, we perform systematic statistical and mutational analyses of the overlapping HIV-1 genes tat and rev and engineer exhaustive libraries of non-overlapped viruses to perform deep mutational scanning of each gene independently. We find a "segregated" organization in which overlapped sites encode functional residues of one gene or the other, but never both. Furthermore, this organization eliminates unfit genotypes, providing a fitness advantage to the population. Our comprehensive analysis reveals the extraordinary manner in which HIV minimizes the constraint of overlapping genes and repurposes that constraint to its own advantage. Thus, overlaps are not just consequences of evolutionary constraints, but rather can provide population fitness advantages.

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.

Statistical inference of a convergent antibody repertoire response to influenza vaccine

Background

Vaccines dramatically affect an individual’s adaptive immune system and thus provide an excellent means to study human immunity. Upon vaccination, the B cells that express antibodies (Abs) that happen to bind the vaccine are stimulated to proliferate and undergo mutagenesis at their Ab locus. This process may alter the composition of B cell lineages within an individual, which are known collectively as the antibody repertoire (AbR). Antibodies are also highly expressed in whole blood, potentially enabling RNA sequencing (RNA-seq) technologies to query this diversity. Less is known about the diversity of AbR responses across individuals to a given vaccine and if individuals tend to yield a similar response to the same antigenic stimulus.

Methods

Here we implement a bioinformatic pipeline that extracts the AbR information from a time-series RNA-seq dataset of five patients who were administered a seasonal trivalent influenza vaccine (TIV). We harness the detailed time-series nature of this dataset and use methods based in functional data analysis (FDA) to identify the Abs that respond to the vaccine. We then design and implement rigorous statistical tests in order to ask whether or not these patients exhibit a convergent AbR response to the same TIV.

Results

We find that high-resolution time-series data can be used to help identify the Abs that respond to an antigenic stimulus and that this response can exhibit a convergent nature across patients inoculated with the same vaccine. However, correlations in AbR diversity among individuals prior to inoculation can confound inference of a convergent signal unless it is taken into account.

Conclusions

We developed a framework to identify the elements of an AbR that respond to an antigen. This information could be used to understand the diversity of different immune responses in different individuals, as well as to gauge the effectiveness of the immune response to a given stimulus within an individual. We also present a framework for testing a convergent hypothesis between AbRs; a hypothesis that is more difficult to test than previously appreciated. Our discovery of a convergent signal suggests that similar epitopes do select for antibodies with similar sequence characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0314-z) contains supplementary material, which is available to authorized users.

Selection and explosive growth alter genetic architecture and hamper the detection of causal rare variants

The role of rare alleles in complex phenotypes has been hotly debated, but most rare variant association tests (RVATs) do not account for the evolutionary forces that affect genetic architecture. Here, we use simulation and numerical algorithms to show that explosive population growth, as experienced by human populations, can dramatically increase the impact of very rare alleles on trait variance. We then assess the ability of RVATs to detect causal loci using simulations and human RNA-seq data. Surprisingly, we find that statistical performance is worst for phenotypes in which genetic variance is due mainly to rare alleles, and explosive population growth decreases power. Although many studies have attempted to identify causal rare variants, few have reported novel associations. This has sometimes been interpreted to mean that rare variants make negligible contributions to complex trait heritability. Our work shows that RVATs are not robust to realistic human evolutionary forces, so general conclusions about the impact of rare variants on complex traits may be premature.

Pooled Sequencing of Candidate Genes Implicates Rare Variants in the Development of Asthma Following Severe RSV Bronchiolitis in Infancy

Severe infection with respiratory syncytial virus (RSV) during infancy is strongly associated with the development of asthma. To identify genetic variation that contributes to asthma following severe RSV bronchiolitis during infancy, we sequenced the coding exons of 131 asthma candidate genes in 182 European and African American children with severe RSV bronchiolitis in infancy using anonymous pools for variant discovery, and then directly genotyped a set of 190 nonsynonymous variants. Association testing was performed for physician-diagnosed asthma before the 7^th birthday (asthma) using genotypes from 6,500 individuals from the Exome Sequencing Project (ESP) as controls to gain statistical power. In addition, among patients with severe RSV bronchiolitis during infancy, we examined genetic associations with asthma, active asthma, persistent wheeze, and bronchial hyperreactivity (methacholine PC₂₀) at age 6 years. We identified four rare nonsynonymous variants that were significantly associated with asthma following severe RSV bronchiolitis, including single variants in ADRB2, FLG and NCAM1 in European Americans (p = 4.6x10^-4, 1.9x10^-13 and 5.0x10^-5, respectively), and NOS1 in African Americans (p = 2.3x10^-11). One of the variants was a highly functional nonsynonymous variant in ADRB2 (rs1800888), which was also nominally associated with asthma (p = 0.027) and active asthma (p = 0.013) among European Americans with severe RSV bronchiolitis without including the ESP. Our results suggest that rare nonsynonymous variants contribute to the development of asthma following severe RSV bronchiolitis in infancy, notably in ADRB2. Additional studies are required to explore the role of rare variants in the etiology of asthma and asthma-related traits following severe RSV bronchiolitis.

CauseMap: fast inference of causality from complex time series

Background. Establishing health-related causal relationships is a central pursuit in biomedical research. Yet, the interdependent non-linearity of biological systems renders causal dynamics laborious and at times impractical to disentangle. This pursuit is further impeded by the dearth of time series that are sufficiently long to observe and understand recurrent patterns of flux. However, as data generation costs plummet and technologies like wearable devices democratize data collection, we anticipate a coming surge in the availability of biomedically-relevant time series data. Given the life-saving potential of these burgeoning resources, it is critical to invest in the development of open source software tools that are capable of drawing meaningful insight from vast amounts of time series data. Results. Here we present CauseMap, the first open source implementation of convergent cross mapping (CCM), a method for establishing causality from long time series data (≳25 observations). Compared to existing time series methods, CCM has the advantage of being model-free and robust to unmeasured confounding that could otherwise induce spurious associations. CCM builds on Takens' Theorem, a well-established result from dynamical systems theory that requires only mild assumptions. This theorem allows us to reconstruct high dimensional system dynamics using a time series of only a single variable. These reconstructions can be thought of as shadows of the true causal system. If reconstructed shadows can predict points from opposing time series, we can infer that the corresponding variables are providing views of the same causal system, and so are causally related. Unlike traditional metrics, this test can establish the directionality of causation, even in the presence of feedback loops. Furthermore, since CCM can extract causal relationships from times series of, e.g., a single individual, it may be a valuable tool to personalized medicine. We implement CCM in Julia, a high-performance programming language designed for facile technical computing. Our software package, CauseMap, is platform-independent and freely available as an official Julia package. Conclusions. CauseMap is an efficient implementation of a state-of-the-art algorithm for detecting causality from time series data. We believe this tool will be a valuable resource for biomedical research and personalized medicine.

Global mapping of herpesvirus-host protein complexes reveals a transcription strategy for late genes

Mapping host-pathogen interactions has proven instrumental for understanding how viruses manipulate host machinery and how numerous cellular processes are regulated. DNA viruses such as herpesviruses have relatively large coding capacity and thus can target an extensive network of cellular proteins. To identify the host proteins hijacked by this pathogen, we systematically affinity tagged and purified all 89 proteins of Kaposi's sarcoma-associated herpesvirus (KSHV) from human cells. Mass spectrometry of this material identified over 500 virus-host interactions. KSHV causes AIDS-associated cancers, and its interaction network is enriched for proteins linked to cancer and overlaps with proteins that are also targeted by HIV-1. We found that the conserved KSHV protein ORF24 binds to RNA polymerase II and brings it to viral late promoters by mimicking and replacing cellular TATA-box-binding protein (TBP). This is required for herpesviral late gene expression, a complex and poorly understood phase of the viral lifecycle.

Genome-wide association study and admixture mapping reveal new loci associated with total IgE levels in Latinos

BACKGROUND

IgE is a key mediator of allergic inflammation, and its levels are frequently increased in patients with allergic disorders.

OBJECTIVE

We sought to identify genetic variants associated with IgE levels in Latinos.

METHODS

We performed a genome-wide association study and admixture mapping of total IgE levels in 3334 Latinos from the Genes-environments & Admixture in Latino Americans (GALA II) study. Replication was evaluated in 454 Latinos, 1564 European Americans, and 3187 African Americans from independent studies.

RESULTS

We confirmed associations of 6 genes identified by means of previous genome-wide association studies and identified a novel genome-wide significant association of a polymorphism in the zinc finger protein 365 gene (ZNF365) with total IgE levels (rs200076616, P = 2.3 × 10(-8)). We next identified 4 admixture mapping peaks (6p21.32-p22.1, 13p22-31, 14q23.2, and 22q13.1) at which local African, European, and/or Native American ancestry was significantly associated with IgE levels. The most significant peak was 6p21.32-p22.1, where Native American ancestry was associated with lower IgE levels (P = 4.95 × 10(-8)). All but 22q13.1 were replicated in an independent sample of Latinos, and 2 of the peaks were replicated in African Americans (6p21.32-p22.1 and 14q23.2). Fine mapping of 6p21.32-p22.1 identified 6 genome-wide significant single nucleotide polymorphisms in Latinos, 2 of which replicated in European Americans. Another single nucleotide polymorphism was peak-wide significant within 14q23.2 in African Americans (rs1741099, P = 3.7 × 10(-6)) and replicated in non-African American samples (P = .011).

CONCLUSION

We confirmed genetic associations at 6 genes and identified novel associations within ZNF365, HLA-DQA1, and 14q23.2. Our results highlight the importance of studying diverse multiethnic populations to uncover novel loci associated with total IgE levels.

Population genetic simulations of complex phenotypes with implications for rare variant association tests

Demographic events and natural selection alter patterns of genetic variation within populations and may play a substantial role in shaping the genetic architecture of complex phenotypes and disease. However, the joint impact of these basic evolutionary forces is often ignored in the assessment of statistical tests of association. Here, we provide a simulation-based framework for generating DNA sequences that incorporates selection and demography with flexible models for simulating phenotypic variation (sfs_coder). This tool also allows the user to perform locus-specific simulations by automatically querying annotated genomic functional elements and genetic maps. We demonstrate the effects of evolutionary forces on patterns of genetic variation by simulating recently inferred models of human selection and demography. We use these simulations to show that the demographic model and locus-specific features, such as the proportion of sites under selection, may have practical implications for estimating the statistical power of sequencing-based rare variant association tests. In particular, for some phenotype models, there may be higher power to detect rare variant associations in African populations compared to non-Africans, but power is considerably reduced in regions of the genome with rampant negative selection. Furthermore, we show that existing methods for simulating large samples based on resampling from a small set of observed haplotypes fail to recapitulate the distribution of rare variants in the presence of rapid population growth (as has been observed in several human populations).

Genetic simulation tools for post-genome wide association studies of complex diseases

Genetic simulation programs are used to model data under specified assumptions to facilitate the understanding and study of complex genetic systems. Standardized data sets generated using genetic simulation are essential for the development and application of novel analytical tools in genetic epidemiology studies. With continuing advances in high-throughput genomic technologies and generation and analysis of larger, more complex data sets, there is a need for updating current approaches in genetic simulation modeling. To provide a forum to address current and emerging challenges in this area, the National Cancer Institute (NCI) sponsored a workshop, entitled "Genetic Simulation Tools for Post-Genome Wide Association Studies of Complex Diseases" at the National Institutes of Health (NIH) in Bethesda, Maryland on March 11-12, 2014. The goals of the workshop were to (1) identify opportunities, challenges, and resource needs for the development and application of genetic simulation models; (2) improve the integration of tools for modeling and analysis of simulated data; and (3) foster collaborations to facilitate development and applications of genetic simulation. During the course of the meeting, the group identified challenges and opportunities for the science of simulation, software and methods development, and collaboration. This paper summarizes key discussions at the meeting, and highlights important challenges and opportunities to advance the field of genetic simulation.

selscan: an efficient multithreaded program to perform EHH-based scans for positive selection

Haplotype-based scans to detect natural selection are useful to identify recent or ongoing positive selection in genomes. As both real and simulated genomic data sets grow larger, spanning thousands of samples and millions of markers, there is a need for a fast and efficient implementation of these scans for general use. Here, we present selscan, an efficient multithreaded application that implements Extended Haplotype Homozygosity (EHH), Integrated Haplotype Score (iHS), and Cross-population EHH (XPEHH). selscan accepts phased genotypes in multiple formats, including TPED, and performs extremely well on both simulated and real data and over an order of magnitude faster than existing available implementations. It calculates iHS on chromosome 22 (22,147 loci) across 204 CEU haplotypes in 353 s on one thread (33 s on 16 threads) and calculates XPEHH for the same data relative to 210 YRI haplotypes in 578 s on one thread (52 s on 16 threads). Source code and binaries (Windows, OSX, and Linux) are available at https://github.com/szpiech/selscan.

Dysbiosis of the gut microbiota is associated with HIV disease progression and tryptophan catabolism

Progressive HIV infection is characterized by dysregulation of the intestinal immune barrier, translocation of immunostimulatory microbial products, and chronic systemic inflammation that is thought to drive progression of disease to AIDS. Elements of this pathologic process persist despite viral suppression during highly active antiretroviral therapy (HAART), and drivers of these phenomena remain poorly understood. Disrupted intestinal immunity can precipitate dysbiosis that induces chronic inflammation in the mucosa and periphery of mice. However, putative microbial drivers of HIV-associated immunopathology versus recovery have not been identified in humans. Using high-resolution bacterial community profiling, we identified a dysbiotic mucosal-adherent community enriched in Proteobacteria and depleted of Bacteroidia members that was associated with markers of mucosal immune disruption, T cell activation, and chronic inflammation in HIV-infected subjects. Furthermore, this dysbiosis was evident among HIV-infected subjects undergoing HAART, and the extent of dysbiosis correlated with activity of the kynurenine pathway of tryptophan catabolism and plasma concentrations of the inflammatory cytokine interleukin-6 (IL-6), two established markers of disease progression. Gut-resident bacteria with capacity to catabolize tryptophan through the kynurenine pathway were found to be enriched in HIV-infected subjects, strongly correlated with kynurenine levels in HIV-infected subjects, and capable of kynurenine production in vitro. These observations demonstrate a link between mucosal-adherent colonic bacteria and immunopathogenesis during progressive HIV infection that is apparent even in the setting of viral suppression during HAART. This link suggests that gut-resident microbial populations may influence intestinal homeostasis during HIV disease.

A genome-wide association study of bronchodilator response in Latinos implicates rare variants

BACKGROUND

The primary rescue medication to treat acute asthma exacerbation is the short-acting β₂-adrenergic receptor agonist; however, there is variation in how well a patient responds to treatment. Although these differences might be due to environmental factors, there is mounting evidence for a genetic contribution to variability in bronchodilator response (BDR).

OBJECTIVE

To identify genetic variation associated with bronchodilator drug response in Latino children with asthma.

METHODS

We performed a genome-wide association study (GWAS) for BDR in 1782 Latino children with asthma using standard linear regression, adjusting for genetic ancestry and ethnicity, and performed replication studies in an additional 531 Latinos. We also performed admixture mapping across the genome by testing for an association between local European, African, and Native American ancestry and BDR, adjusting for genomic ancestry and ethnicity.

RESULTS

We identified 7 genetic variants associated with BDR at a genome-wide significant threshold (P < 5 × 10(-8)), all of which had frequencies of less than 5%. Furthermore, we observed an excess of small P values driven by rare variants (frequency, <5%) and by variants in the proximity of solute carrier (SLC) genes. Admixture mapping identified 5 significant peaks; fine mapping within these peaks identified 2 rare variants in SLC22A15 as being associated with increased BDR in Mexicans. Quantitative PCR and immunohistochemistry identified SLC22A15 as being expressed in the lung and bronchial epithelial cells.

CONCLUSION

Our results suggest that rare variation contributes to individual differences in response to albuterol in Latinos, notably in SLC genes that include membrane transport proteins involved in the transport of endogenous metabolites and xenobiotics. Resequencing in larger, multiethnic population samples and additional functional studies are required to further understand the role of rare variation in BDR.

Population genetics of rare variants and complex diseases

OBJECTIVES

Identifying drivers of complex traits from the noisy signals of genetic variation obtained from high-throughput genome sequencing technologies is a central challenge faced by human geneticists today. We hypothesize that the variants involved in complex diseases are likely to exhibit non-neutral evolutionary signatures. Uncovering the evolutionary history of all variants is therefore of intrinsic interest for complex disease research. However, doing so necessitates the simultaneous elucidation of the targets of natural selection and population-specific demographic history.

METHODS

Here we characterize the action of natural selection operating across complex disease categories, and use population genetic simulations to evaluate the expected patterns of genetic variation in large samples. We focus on populations that have experienced historical bottlenecks followed by explosive growth (consistent with many human populations), and describe the differences between evolutionarily deleterious mutations and those that are neutral.

RESULTS

Genes associated with several complex disease categories exhibit stronger signatures of purifying selection than non-disease genes. In addition, loci identified through genome-wide association studies of complex traits also exhibit signatures consistent with being in regions recurrently targeted by purifying selection. Through simulations, we show that population bottlenecks and rapid growth enable deleterious rare variants to persist at low frequencies just as long as neutral variants, but low-frequency and common variants tend to be much younger than neutral variants. This has resulted in a large proportion of modern-day rare alleles that have a deleterious effect on function and that potentially contribute to disease susceptibility.

CONCLUSIONS

The key question for sequencing-based association studies of complex traits is how to distinguish between deleterious and benign genetic variation. We used population genetic simulations to uncover patterns of genetic variation that distinguish these two categories, especially derived allele age, thereby providing inroads into novel methods for characterizing rare genetic variation driving complex diseases.

Case-control admixture mapping in Latino populations enriches for known asthma-associated genes

BACKGROUND

Polymorphisms in more than 100 genes have been associated with asthma susceptibility, yet much of the heritability remains to be explained. Asthma disproportionately affects different racial and ethnic groups in the United States, suggesting that admixture mapping is a useful strategy to identify novel asthma-associated loci.

OBJECTIVE

We sought to identify novel asthma-associated loci in Latino populations using case-control admixture mapping.

METHODS

We performed genome-wide admixture mapping by comparing levels of local Native American, European, and African ancestry between children with asthma and nonasthmatic control subjects in Puerto Rican and Mexican populations. Within candidate peaks, we performed allelic tests of association, controlling for differences in local ancestry.

RESULTS

Between the 2 populations, we identified a total of 62 admixture mapping peaks at a P value of less than 10(-3) that were significantly enriched for previously identified asthma-associated genes (P= .0051). One of the peaks was statistically significant based on 100 permutations in the Mexican sample (6q15); however, it was not significant in Puerto Rican subjects. Another peak was identified at nominal significance in both populations (8q12); however, the association was observed with different ancestries.

CONCLUSION

Case-control admixture mapping is a promising strategy for identifying novel asthma-associated loci in Latino populations and implicates genetic variation at 6q15 and 8q12 regions with asthma susceptibility. This approach might be useful for identifying regions that contribute to both shared and population-specific differences in asthma susceptibility.

A fine-scale chimpanzee genetic map from population sequencing

To study the evolution of recombination rates in apes, we developed methodology to construct a fine-scale genetic map from high-throughput sequence data from 10 Western chimpanzees, Pan troglodytes verus. Compared to the human genetic map, broad-scale recombination rates tend to be conserved, but with exceptions, particularly in regions of chromosomal rearrangements and around the site of ancestral fusion in human chromosome 2. At fine scales, chimpanzee recombination is dominated by hotspots, which show no overlap with those of humans even though rates are similarly elevated around CpG islands and decreased within genes. The hotspot-specifying protein PRDM9 shows extensive variation among Western chimpanzees, and there is little evidence that any sequence motifs are enriched in hotspots. The contrasting locations of hotspots provide a natural experiment, which demonstrates the impact of recombination on base composition.

Global landscape of HIV-human protein complexes

Human immunodeficiency virus (HIV) has a small genome and therefore relies heavily on the host cellular machinery to replicate. Identifying which host proteins and complexes come into physical contact with the viral proteins is crucial for a comprehensive understanding of how HIV rewires the host's cellular machinery during the course of infection. Here we report the use of affinity tagging and purification mass spectrometry to determine systematically the physical interactions of all 18 HIV-1 proteins and polyproteins with host proteins in two different human cell lines (HEK293 and Jurkat). Using a quantitative scoring system that we call MiST, we identified with high confidence 497 HIV-human protein-protein interactions involving 435 individual human proteins, with ∼40% of the interactions being identified in both cell types. We found that the host proteins hijacked by HIV, especially those found interacting in both cell types, are highly conserved across primates. We uncovered a number of host complexes targeted by viral proteins, including the finding that HIV protease cleaves eIF3d, a subunit of eukaryotic translation initiation factor 3. This host protein is one of eleven identified in this analysis that act to inhibit HIV replication. This data set facilitates a more comprehensive and detailed understanding of how the host machinery is manipulated during the course of HIV infection.

A population genetics-phylogenetics approach to inferring natural selection in coding sequences

Through an analysis of polymorphism within and divergence between species, we can hope to learn about the distribution of selective effects of mutations in the genome, changes in the fitness landscape that occur over time, and the location of sites involved in key adaptations that distinguish modern-day species. We introduce a novel method for the analysis of variation in selection pressures within and between species, spatially along the genome and temporally between lineages. We model codon evolution explicitly using a joint population genetics-phylogenetics approach that we developed for the construction of multiallelic models with mutation, selection, and drift. Our approach has the advantage of performing direct inference on coding sequences, inferring ancestral states probabilistically, utilizing allele frequency information, and generalizing to multiple species. We use a Bayesian sliding window model for intragenic variation in selection coefficients that efficiently combines information across sites and captures spatial clustering within the genome. To demonstrate the utility of the method, we infer selective pressures acting in Drosophila melanogaster and D. simulans from polymorphism and divergence data for 100 X-linked coding regions.

Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations

Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies of asthma in 5,416 individuals with asthma (cases) including individuals of European American, African American or African Caribbean, and Latino ancestry, with replication in an additional 12,649 individuals from the same ethnic groups. We identified five susceptibility loci. Four were at previously reported loci on 17q21, near IL1RL1, TSLP and IL33, but we report for the first time, to our knowledge, that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a new asthma susceptibility locus at PYHIN1, with the association being specific to individuals of African descent (P = 3.9 × 10(-9)). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.

SNP discovery, expression and cis-regulatory variation in the UGT2B genes

UGT2B enzymes metabolize multiple endogenous and exogenous molecules, including steroid hormones and clinical drugs. However, little is known about the inter-individual variation in gene expression and its determinants. We re-sequenced candidate regulatory regions and the partial coding regions (41.1 kb) of UGT2B genes and identified 332 genetic variants. We measured gene expression in normal breast and liver samples and observed different patterns. The expression levels varied greatly across individuals in both tissues and were significantly correlated with each other in liver. Genotyping of tagging single-nucleotide polymorphisms (SNPs) in the same samples and association tests between genotype and transcript levels identified 62 variants that were associated with at least one UGT2B mRNA levels in either tissue. Most of these cis-regulatory SNPs were not shared between tissues, suggesting that this gene family is regulated in a tissue-specific manner. Our results provide insight into studying the role of UGT2B variation in hormone-dependent cancers and drug response.

Classic selective sweeps were rare in recent human evolution

Efforts to identify the genetic basis of human adaptations from polymorphism data have sought footprints of "classic selective sweeps" (in which a beneficial mutation arises and rapidly fixes in the population).Yet it remains unknown whether this form of natural selection was common in our evolution. We examined the evidence for classic sweeps in resequencing data from 179 human genomes. As expected under a recurrent-sweep model, we found that diversity levels decrease near exons and conserved noncoding regions. In contrast to expectation, however, the trough in diversity around human-specific amino acid substitutions is no more pronounced than around synonymous substitutions. Moreover, relative to the genome background, amino acid and putative regulatory sites are not significantly enriched in alleles that are highly differentiated between populations. These findings indicate that classic sweeps were not a dominant mode of human adaptation over the past ~250,000 years.

Inferring the joint demographic history of multiple populations from multidimensional SNP frequency data

Demographic models built from genetic data play important roles in illuminating prehistorical events and serving as null models in genome scans for selection. We introduce an inference method based on the joint frequency spectrum of genetic variants within and between populations. For candidate models we numerically compute the expected spectrum using a diffusion approximation to the one-locus, two-allele Wright-Fisher process, involving up to three simultaneous populations. Our approach is a composite likelihood scheme, since linkage between neutral loci alters the variance but not the expectation of the frequency spectrum. We thus use bootstraps incorporating linkage to estimate uncertainties for parameters and significance values for hypothesis tests. Our method can also incorporate selection on single sites, predicting the joint distribution of selected alleles among populations experiencing a bevy of evolutionary forces, including expansions, contractions, migrations, and admixture. We model human expansion out of Africa and the settlement of the New World, using 5 Mb of noncoding DNA resequenced in 68 individuals from 4 populations (YRI, CHB, CEU, and MXL) by the Environmental Genome Project. We infer divergence between West African and Eurasian populations 140 thousand years ago (95% confidence interval: 40–270 kya). This is earlier than other genetic studies, in part because we incorporate migration. We estimate the European (CEU) and East Asian (CHB) divergence time to be 23 kya (95% c.i.: 17–43 kya), long after archeological evidence places modern humans in Europe. Finally, we estimate divergence between East Asians (CHB) and Mexican-Americans (MXL) of 22 kya (95% c.i.: 16.3–26.9 kya), and our analysis yields no evidence for subsequent migration. Furthermore, combining our demographic model with a previously estimated distribution of selective effects among newly arising amino acid mutations accurately predicts the frequency spectrum of nonsynonymous variants across three continental populations (YRI, CHB, CEU).

The demographic history of our species is reflected in patterns of genetic variation within and among populations. We developed an efficient method for calculating the expected distribution of genetic variation, given a demographic model including such events as population size changes, population splits and joins, and migration. We applied our approach to publicly available human sequencing data, searching for models that best reproduce the observed patterns. Our joint analysis of data from African, European, and Asian populations yielded new dates for when these populations diverged. In particular, we found that African and Eurasian populations diverged around 100,000 years ago. This is earlier than other genetic studies suggest, because our model includes the effects of migration, which we found to be important for reproducing observed patterns of variation in the data. We also analyzed data from European, Asian, and Mexican populations to model the peopling of the Americas. Here, we find no evidence for recurrent migration after East Asian and Native American populations diverged. Our methods are not limited to studying humans, and we hope that future sequencing projects will offer more insights into the history of both our own species and others.

Evolutionary processes acting on candidate cis-regulatory regions in humans inferred from patterns of polymorphism and divergence

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.

It has been suggested that changes in gene expression may have played a more important role in the evolution of modern humans than changes in protein-coding sequences. In order to identify signatures of natural selection on candidate cis-regulatory regions, we examined single nucleotide polymorphisms obtained from the complete re-sequencing of conserved non-coding sites (CNCs) in the flanking regions of over 15,000 genes in 35 humans. Patterns of allele frequencies in CNCs indicate the presence of both positive and negative selection acting on standing variation within these candidate cis-regulatory regions, particularly for the 5′ and 3′ UTRs of genes. Gene-specific tests comparing levels of polymorphism and divergence identify several genes with strong signatures of selection on candidate cis-regulatory regions and suggest that the biological characteristics of genes subject to selection are different between coding and candidate cis-regulatory regions with respect to gene expression and function. For example, we find stronger signatures of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, which we do not observe in a concurrent analysis on protein-coding regions. Our results suggest that both positive and negative selection have acted on candidate cis-regulatory regions and that the evolution of non-coding DNA has played an important role throughout hominid evolution.

Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds

The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.

A flexible forward simulator for populations subject to selection and demography

This article introduces a new forward population genetic simulation program that can efficiently generate samples from populations with complex demographic histories under various models of natural selection. The program (SFS_CODE) is highly flexible, allowing the user to simulate realistic genomic regions with several loci evolving according to a variety of mutation models (from simple to context-dependent), and allows for insertions and deletions. Each locus can be annotated as either coding or non-coding, sex-linked or autosomal, selected or neutral, and have an arbitrary linkage structure (from completely linked to independent).

AVAILABILITY

The source code (written in the C programming language) is available at http://sfscode.sourceforge.net, and a web server (http://cbsuapps.tc.cornell.edu/sfscode.aspx) allows the user to perform simulations using the high-performance computing cluster hosted by the Cornell University Computational Biology Service Unit.