Genomic analyses inform on migration events during the peopling of Eurasia.

High-coverage whole-genome sequence studies have so far focused on a limited number of geographically restricted populations, or been targeted at specific diseases, such as cancer. Nevertheless, the availability of high-resolution genomic data has led to the development of new methodologies for inferring population history and refuelled the debate on the mutation rate in humans. Here we present the Estonian Biocentre Human Genome Diversity Panel (EGDP), a dataset of 483 high-coverage human genomes from 148 populations worldwide, including 379 new genomes from 125 populations, which we group into diversity and selection sets. We analyse this dataset to refine estimates of continent-wide patterns of heterozygosity, long- and short-distance gene flow, archaic admixture, and changes in effective population size through time as well as for signals of positive or balancing selection. We find a genetic signature in present-day Papuans that suggests that at least 2% of their genome originates from an early and largely extinct expansion of anatomically modern humans (AMHs) out of Africa. Together with evidence from the western Asian fossil record, and admixture between AMHs and Neanderthals predating the main Eurasian expansion, our results contribute to the mounting evidence for the presence of AMHs out of Africa earlier than 75,000 years ago.

The association of genomic lesions and PD-1/PD-L1 expression in resected triple-negative breast cancers.

BACKGROUND: Elevated PD-L1 expression on tumor cells, a context associated with an adaptive immune response, has been linked to the total burden of copy number variants (CNVs) in aneuploid tumors, to microsatellite instability (MSI), and to specific genomic driver lesions, including loss of PTEN, MYC amplification, and activating mutations in driver oncogenes such as KRAS and PIK3CA. Triple-negative breast cancers (TNBCs) typically have high levels of CNVs and diverse driver lesions in their genomes. Thus, there is significant interest in exploiting genomic data to develop predictive immunotherapy biomarkers for patients with TNBC. METHODS: Whole tissue samples from 55 resected TNBCs were screened by immunohistochemistry (IHC) for PD-1 and PD-L1 by using validated antibodies and established scoring methods for staining of tumor and non-tumor cells. In parallel, we interrogated biopsies from each resection with DNA content flow cytometry and sorted the nuclei of diploid, tetraploid, and aneuploid cell populations. CNVs were mapped with CNV oligonucleotide arrays by using purified (>95%) tumor populations. We generated whole exome data for 12 sorted tumor samples to increase the resolution within loci of interest and to incorporate somatic mutations into our genomic signatures. RESULTS AND CONCLUSIONS: PD-L1 staining was detected on tumor cells in 29 out of 54 (54%) evaluable cases and was associated with increased overall survival (P = 0.0024). High levels of PD-1 and PD-L1 (IHC ≥4) were present in 11 out of 54 (20%) and 20 out of 54 (37%) cases with staining of PD-L1 primarily on tumor cells for 17 out of 20 (85%) cases. The latter included tumors with both high (>50) and low (<20) numbers of CNVs. Notably, homozygous deletion of PTEN (n = 6) or activating mutation in PIK3CA (n = 1) was not associated with increased expression of either immune checkpoint activator in TNBC. In contrast, two treatment-naïve cases with EGFR driver amplicons had high PD-L1 tumor staining. High mutational load and predicted neoepitopes were observed in MSI+ and high CNV burden TNBCs but were not associated with high PD-L1 expression on tumor cells. Our results challenge current models of genomic-based immunotherapy signatures yet suggest that discrete genomic lesions may complement existing biomarkers to advance immune checkpoint therapies for patients with TNBC.

A recent bottleneck of Y chromosome diversity coincides with a global change in culture.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Natural selection reduced diversity on human y chromosomes.

The human Y chromosome exhibits surprisingly low levels of genetic diversity. This could result from neutral processes if the effective population size of males is reduced relative to females due to a higher variance in the number of offspring from males than from females. Alternatively, selection acting on new mutations, and affecting linked neutral sites, could reduce variability on the Y chromosome. Here, using genome-wide analyses of X, Y, autosomal and mitochondrial DNA, in combination with extensive population genetic simulations, we show that low observed Y chromosome variability is not consistent with a purely neutral model. Instead, we show that models of purifying selection are consistent with observed Y diversity. Further, the number of sites estimated to be under purifying selection greatly exceeds the number of Y-linked coding sites, suggesting the importance of the highly repetitive ampliconic regions. While we show that purifying selection removing deleterious mutations can explain the low diversity on the Y chromosome, we cannot exclude the possibility that positive selection acting on beneficial mutations could have also reduced diversity in linked neutral regions, and may have contributed to lowering human Y chromosome diversity. Because the functional significance of the ampliconic regions is poorly understood, our findings should motivate future research in this area.

Gene survival and death on the human Y chromosome.

Y chromosomes have long been dismissed as "graveyards of genes," but there is still much to be learned from the genetic relics of genes that were once functional on the human Y. We identified human X-linked genes whose gametologs have been pseudogenized or completely lost from the Y chromosome and inferred which evolutionary forces may be acting to retain genes on the Y. Although gene loss appears to be largely correlated with the suppression of recombination, we observe that X-linked genes with functional Y homologs evolve under stronger purifying selection and are expressed at higher levels than X-linked genes with nonfunctional Y homologs. Additionally, we support and expand upon the hypothesis that X inactivation is primarily driven by gene loss on the Y. Using linear discriminant analysis, we show that X-inactivation status can successfully classify 90% of X-linked genes into those with functional or nonfunctional Y homologs.

A scan for human-specific relaxation of negative selection reveals unexpected polymorphism in proteasome genes.

Environmental or genomic changes during evolution can relax negative selection pressure on specific loci, permitting high frequency polymorphisms at previously conserved sites. Here, we jointly analyze population genomic and comparative genomic data to search for functional processes showing relaxed negative selection specifically in the human lineage, whereas remaining evolutionarily conserved in other mammals. Consistent with previous studies, we find that olfactory receptor genes display such a signature of relaxation in humans. Intriguingly, proteasome genes also show a prominent signal of human-specific relaxation: multiple proteasome subunits, including four members of the catalytic core particle, contain high frequency nonsynonymous polymorphisms at sites conserved across mammals. Chimpanzee proteasome genes do not display a similar trend. Human proteasome genes also bear no evidence of recent positive or balancing selection. These results suggest human-specific relaxation of negative selection in proteasome subunits; the exact biological causes, however, remain unknown.

Genome analyses substantiate male mutation bias in many species.

In many species the mutation rate is higher in males than in females, a phenomenon denoted as male mutation bias. This is often observed in animals where males produce many more sperm than females produce eggs, and is thought to result from differences in the number of replication-associated mutations accumulated in each sex. Thus, studies of male mutation bias have the capacity to reveal information about the replication-dependent or replication-independent nature of different mutations. The availability of whole genome sequences for many species, as well as for multiple individuals within a species, has opened the door to studying factors, both sequence-specific and those acting on the genome globally, that affect differences in mutation rates between males and females. Here, we assess the advantages that genomic sequences provide for studies of male mutation bias and general mutation mechanisms, discuss major challenges left unresolved, and speculate about the direction of future studies.

HGV2011: personalized genomic medicine meets the incidentalome.

The 12th International Meeting on Human Genome Variation and Complex Genome Analysis (HGV2011: Berkeley, California, USA, 8th-10th September 2011) was a stimulating workshop where researchers from academia and industry explored the latest progress, challenges, and opportunities in genome variation research. Key themes included progress beyond GWAS, variation in human populations, use of sequence data in medical settings, large-scale sequencing data analysis, and bioinformatics approaches to large datasets.

Genetic Diversity on the Sex Chromosomes.

Levels and patterns of genetic diversity can provide insights into a population's history. In species with sex chromosomes, differences between genomic regions with unique inheritance patterns can be used to distinguish between different sets of possible demographic and selective events. This review introduces the differences in population history for sex chromosomes and autosomes, provides the expectations for genetic diversity across the genome under different evolutionary scenarios, and gives an introductory description for how deviations in these expectations are calculated and can be interpreted. Predominantly, diversity on the sex chromosomes has been used to explore and address three research areas: 1) Mating patterns and sex-biased variance in reproductive success, 2) signatures of selection, and 3) evidence for modes of speciation and introgression. After introducing the theory, this review catalogs recent studies of genetic diversity on the sex chromosomes across species within the major research areas that sex chromosomes are typically applied to, arguing that there are broad similarities not only between male-heterogametic (XX/XY) and female-heterogametic (ZZ/ZW) sex determination systems but also any mating system with reduced recombination in a sex-determining region. Further, general patterns of reduced diversity in nonrecombining regions are shared across plants and animals. There are unique patterns across populations with vastly different patterns of mating and speciation, but these do not tend to cluster by taxa or sex determination system.

Genome Sequence of a Single-Stranded DNA Virus Identified in Gila Monster Feces.

The Gila monster (Heloderma suspectum) is native to the Sonoran Desert. Metagenomic analyses of a Gila monster fecal sample revealed the presence of a small, circular, single-stranded DNA virus that is most closely related to a gemykrogvirus (family Genomoviridae) genome from caribou feces sharing 88% genome-wide pairwise identity.

The synthetic histone-binding regulator protein PcTF activates interferon genes in breast cancer cells.

BACKGROUND: Mounting evidence from genome-wide studies of cancer shows that chromatin-mediated epigenetic silencing at large cohorts of genes is strongly linked to a poor prognosis. This mechanism is thought to prevent cell differentiation and enable evasion of the immune system. Drugging the cancer epigenome with small molecule inhibitors to release silenced genes from the repressed state has emerged as a powerful approach for cancer research and drug development. Targets of these inhibitors include chromatin-modifying enzymes that can acquire drug-resistant mutations. In order to directly target a generally conserved feature, elevated trimethyl-lysine 27 on histone H3 (H3K27me3), we developed the Polycomb-based Transcription Factor (PcTF), a fusion activator that targets methyl-histone marks via its N-terminal H3K27me3-binding motif, and co-regulates sets of silenced genes. RESULTS: Here, we report transcriptome profiling analyses of PcTF-treated breast cancer model cell lines. We identified a set of 19 PcTF-upregulated genes, or PUGs, that were consistent across three distinct breast cancer cell lines. These genes are associated with the interferon response pathway. CONCLUSIONS: Our results demonstrate for the first time a chromatin-mediated interferon-related transcriptional response driven by an engineered fusion protein that physically links repressive histone marks with active transcription.

Human Papilloma Virus Specific Immunogenicity and Dysfunction of CD8+ T Cells in Head and Neck Cancer.

Human papillomavirus subtype 16 (HPV16) is the primary cause of an increasing number of head and neck squamous cell carcinomas (HNSCC), providing strong rationale for T-cell immune therapies against HPV+ HNSCC. Here we assess immunogenicity of HPV16-specific CD8+ T cells (CTL) and characterize HPV-specific mechanisms of T-cell dysfunction. We identified 16 strong and 29 moderately immunogenic CTL-epitopes from HPV16 E2, E6, and E7 antigens restricted by 12 common HLA class I alleles. E2-specific CTL-reactivity was higher in patients with HPV+ HNSCC than in healthy controls (>3-fold; P = 0.026). Patient-derived E2, E6, and E7 peripheral CTLs exhibited heterogeneity in dysfunctional phenotypes. Immunogenomic analyses of 119 HNSCC transcriptomes revealed high T-cell infiltration and dysfunction in HPV+ HNSCC and correlation of HPV antigen expression with T-cell exhaustion gene signatures. Indoleamine 2,3-dioxygenase (IDO-1) was strongly expressed in HPV+ HNSCC versus HPV- HNSCC (P = 0.001) and correlated with E7 expression (R 2 = 0.84; P = 0.033). Combination treatment with PD-1 blockade and IDO-1 inhibition overcame profound CTL-dysfunction, enhancing HPV+ HNSCC sensitivity to CTL-cytotoxicity in vitro (up to 10-fold in E7-CTLs, P = 0.011). Our findings implicate mechanisms of T-cell escape in HPV+ HNSCC, wherein high tumoral HPV-antigen load results in high expression of immune dysfunction genes on tumor cells (e.g., IDO-1), and dysfunction of HPV-specific CTLs (e.g., E7, E2-CTLs). The HPV16 CTL-epitopes identified in this study, in combination with blockade of HPV+ HNSCC-specific PD-1/IDO-1 checkpoints, may be useful for targeted immunotherapy.Significance: This study evaluates the HPV antigen T-cell immunogenicity role of inhibitory receptors and other exhaustion markers in the cytotoxic function of HPV antigen-specific CTLs and identifies combined inhibition of PD-1/IDO-1 as a strategy to enhance CTL targeting of HPV+ HNSCC. Cancer Res; 78(21); 6159-70. ©2018 AACR.

Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards.

Squamates include all lizards and snakes, and display some of the most diverse and extreme morphological adaptations among vertebrates. However, compared with birds and mammals, relatively few resources exist for comparative genomic analyses of squamates, hampering efforts to understand the molecular bases of phenotypic diversification in such a speciose clade. In particular, the ∼400 species of anole lizard represent an extensive squamate radiation. Here, we sequence and assemble the draft genomes of three anole species-Anolis frenatus, Anolis auratus, and Anolis apletophallus-for comparison with the available reference genome of Anolis carolinensis. Comparative analyses reveal a rapid background rate of molecular evolution consistent with a model of punctuated equilibrium, and strong purifying selection on functional genomic elements in anoles. We find evidence for accelerated evolution in genes involved in behavior, sensory perception, and reproduction, as well as in genes regulating limb bud development and hindlimb specification. Morphometric analyses of anole fore and hindlimbs corroborated these findings. We detect signatures of positive selection across several genes related to the development and regulation of the forebrain, hormones, and the iguanian lizard dewlap, suggesting molecular changes underlying behavioral adaptations known to reinforce species boundaries were a key component in the diversification of anole lizards.

Bioinformatics core competencies for undergraduate life sciences education.

Although bioinformatics is becoming increasingly central to research in the life sciences, bioinformatics skills and knowledge are not well integrated into undergraduate biology education. This curricular gap prevents biology students from harnessing the full potential of their education, limiting their career opportunities and slowing research innovation. To advance the integration of bioinformatics into life sciences education, a framework of core bioinformatics competencies is needed. To that end, we here report the results of a survey of biology faculty in the United States about teaching bioinformatics to undergraduate life scientists. Responses were received from 1,260 faculty representing institutions in all fifty states with a combined capacity to educate hundreds of thousands of students every year. Results indicate strong, widespread agreement that bioinformatics knowledge and skills are critical for undergraduate life scientists as well as considerable agreement about which skills are necessary. Perceptions of the importance of some skills varied with the respondent's degree of training, time since degree earned, and/or the Carnegie Classification of the respondent's institution. To assess which skills are currently being taught, we analyzed syllabi of courses with bioinformatics content submitted by survey respondents. Finally, we used the survey results, the analysis of the syllabi, and our collective research and teaching expertise to develop a set of bioinformatics core competencies for undergraduate biology students. These core competencies are intended to serve as a guide for institutions as they work to integrate bioinformatics into their life sciences curricula.

Evolution of Dosage Compensation in Anolis carolinensis, a Reptile with XX/XY Chromosomal Sex Determination.

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.

Variable Autosomal and X Divergence Near and Far from Genes Affects Estimates of Male Mutation Bias in Great Apes.

Male mutation bias, when more mutations are passed on via the male germline than via the female germline, is observed across mammals. One common way to infer the magnitude of male mutation bias, α, is to compare levels of neutral sequence divergence between genomic regions that spend different amounts of time in the male and female germline. For great apes, including human, we show that estimates of divergence are reduced in putatively unconstrained regions near genes relative to unconstrained regions far from genes. Divergence increases with increasing distance from genes on both the X chromosome and autosomes, but increases faster on the X chromosome than autosomes. As a result, ratios of X/A divergence increase with increasing distance from genes and corresponding estimates of male mutation bias are significantly higher in intergenic regions near genes versus far from genes. Future studies in other species will need to carefully consider the effect that genomic location will have on estimates of male mutation bias.

Genomic signatures of sex-biased demography: progress and prospects.

Sex-biased demographic events have played a crucial role in shaping human history. Many of these processes affect genetic variation and can therefore leave detectable signatures in the genome because autosomal, X-linked, Y-linked, and mitochondrial DNA inheritance differ between sexes. Here, we discuss how sex-biased processes shape patterns of genetic diversity across the genome, review recent genomic evidence for sex-biased demography in modern human populations, and suggest directions for future research.

Genetic Diversity on the Human X Chromosome Does Not Support a Strict Pseudoautosomal Boundary.

Unlike the autosomes, recombination between the X chromosome and the Y chromosome is often thought to be constrained to two small pseudoautosomal regions (PARs) at the tips of each sex chromosome. PAR1 spans the first 2.7 Mb of the proximal arm of the human sex chromosomes, whereas the much smaller PAR2 encompasses the distal 320 kb of the long arm of each sex chromosome. In addition to PAR1 and PAR2, there is a human-specific X-transposed region that was duplicated from the X to the Y chromosome. The X-transposed region is often not excluded from X-specific analyses, unlike the PARs, because it is not thought to routinely recombine. Genetic diversity is expected to be higher in recombining regions than in nonrecombining regions because recombination reduces the effect of linked selection. In this study, we investigated patterns of genetic diversity in noncoding regions across the entire X chromosome of a global sample of 26 unrelated genetic females. We found that genetic diversity in PAR1 is significantly greater than in the nonrecombining regions (nonPARs). However, rather than an abrupt drop in diversity at the pseudoautosomal boundary, there is a gradual reduction in diversity from the recombining through the nonrecombining regions, suggesting that recombination between the human sex chromosomes spans across the currently defined pseudoautosomal boundary. A consequence of recombination spanning this boundary potentially includes increasing the rate of sex-linked disorders (e.g., de la Chapelle) and sex chromosome aneuploidies. In contrast, diversity in PAR2 is not significantly elevated compared to the nonPARs, suggesting that recombination is not obligatory in PAR2. Finally, diversity in the X-transposed region is higher than in the surrounding nonPARs, providing evidence that recombination may occur with some frequency between the X and Y chromosomes in the X-transposed region.

Punctuated bursts in human male demography inferred from 1,244 worldwide Y-chromosome sequences.

We report the sequences of 1,244 human Y chromosomes randomly ascertained from 26 worldwide populations by the 1000 Genomes Project. We discovered more than 65,000 variants, including single-nucleotide variants, multiple-nucleotide variants, insertions and deletions, short tandem repeats, and copy number variants. Of these, copy number variants contribute the greatest predicted functional impact. We constructed a calibrated phylogenetic tree on the basis of binary single-nucleotide variants and projected the more complex variants onto it, estimating the number of mutations for each class. Our phylogeny shows bursts of extreme expansion in male numbers that have occurred independently among each of the five continental superpopulations examined, at times of known migrations and technological innovations.