Transcriptional programs of neoantigen-specific TIL in anti-PD-1-treated lung cancers.

PD-1 blockade unleashes CD8 T cells1, including those specific for mutation-associated neoantigens (MANA), but factors in the tumour microenvironment can inhibit these T cell responses. Single-cell transcriptomics have revealed global T cell dysfunction programs in tumour-infiltrating lymphocytes (TIL). However, the majority of TIL do not recognize tumour antigens2, and little is known about transcriptional programs of MANA-specific TIL. Here, we identify MANA-specific T cell clones using the MANA functional expansion of specific T cells assay3 in neoadjuvant anti-PD-1-treated non-small cell lung cancers (NSCLC). We use their T cell receptors as a 'barcode' to track and analyse their transcriptional programs in the tumour microenvironment using coupled single-cell RNA sequencing and T cell receptor sequencing. We find both MANA- and virus-specific clones in TIL, regardless of response, and MANA-, influenza- and Epstein-Barr virus-specific TIL each have unique transcriptional programs. Despite exposure to cognate antigen, MANA-specific TIL express an incompletely activated cytolytic program. MANA-specific CD8 T cells have hallmark transcriptional programs of tissue-resident memory (TRM) cells, but low levels of interleukin-7 receptor (IL-7R) and are functionally less responsive to interleukin-7 (IL-7) compared with influenza-specific TRM cells. Compared with those from responding tumours, MANA-specific clones from non-responding tumours express T cell receptors with markedly lower ligand-dependent signalling, are largely confined to HOBIThigh TRM subsets, and coordinately upregulate checkpoints, killer inhibitory receptors and inhibitors of T cell activation. These findings provide important insights for overcoming resistance to PD-1 blockade.

Mobile interspersed repeats are major structural variants in the human genome.

Characterizing structural variants in the human genome is of great importance, but a genome wide analysis to detect interspersed repeats has not been done. Thus, the degree to which mobile DNAs contribute to genetic diversity, heritable disease, and oncogenesis remains speculative. We perform transposon insertion profiling by microarray (TIP-chip) to map human L1(Ta) retrotransposons (LINE-1 s) genome-wide. This identified numerous novel human L1(Ta) insertional polymorphisms with highly variant allelic frequencies. We also explored TIP-chip's usefulness to identify candidate alleles associated with different phenotypes in clinical cohorts. Our data suggest that the occurrence of new insertions is twice as high as previously estimated, and that these repeats are under-recognized as sources of human genomic and phenotypic diversity. We have just begun to probe the universe of human L1(Ta) polymorphisms, and as TIP-chip is applied to other insertions such as Alu SINEs, it will expand the catalog of genomic variants even further.

Single-cell atlas of epithelial and stromal cell heterogeneity by lobe and strain in the mouse prostate.

BACKGROUND: Evaluating the complex interplay of cell types in the tissue microenvironment is critical to understanding the origin and progression of diseases in the prostate and potential opportunities for intervention. Mouse models are an essential tool to investigate the molecular and cell-type-specific contributions of prostate disease at an organismal level. While there are well-documented differences in the extent, timing, and nature of disease development in various genetically engineered and exposure-based mouse models in different mouse strains and prostate lobes within each mouse strain, the underlying molecular phenotypic differences in cell types across mouse strains and prostate lobes are incompletely understood. METHODS: In this study, we used single-cell RNA-sequencing (scRNA-seq) methods to assess the single-cell transcriptomes of 6-month-old mouse prostates from two commonly used mouse strains, friend virus B/NIH jackson (FVB/NJ) (N = 2) and C57BL/6J (N = 3). For each mouse, the lobes of the prostate were dissected (anterior, dorsal, lateral, and ventral), and individual scRNA-seq libraries were generated. In situ and pathological analyses were used to explore the spatial and anatomical distributions of novel cell types and molecular markers defining these cell types. RESULTS: Data dimensionality reduction and clustering analysis of scRNA-seq data revealed that basal and luminal cells possessed strain-specific transcriptomic differences, with luminal cells also displaying marked lobe-specific differences. Gene set enrichment analysis comparing luminal cells by strain showed enrichment of proto-Oncogene targets in FVB/NJ mice. Additionally, three rare populations of epithelial cells clustered independently of strain and lobe: one population of luminal cells expressing Foxi1 and components of the vacuolar ATPase proton pump (Atp6v0d2 and Atp6v1g3), another population expressing Psca and other stem cell-associated genes (Ly6a/Sca-1, Tacstd2/Trop-2), and a neuroendocrine population expressing Chga, Chgb, and Syp. In contrast, stromal cell clusters, including fibroblasts, smooth muscle cells, endothelial cells, pericytes, and immune cell types, were conserved across strain and lobe, clustering largely by cell type and not by strain or lobe. One notable exception to this was the identification of two distinct fibroblast populations that we term subglandular fibroblasts and interstitial fibroblasts based on their strikingly distinct spatial distribution in the mouse prostate. CONCLUSIONS: Altogether, these data provide a practical reference of the transcriptional profiles of mouse prostate from two commonly used mouse strains and across all four prostate lobes.

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata.

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis.

This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.

De novo genome assembly of Candida glabrata reveals cell wall protein complement and structure of dispersed tandem repeat arrays.

Candida glabratais an opportunistic pathogen in humans, responsible for approximately 20% of disseminated candidiasis. Candida glabrata's ability to adhere to host tissue is mediated by GPI-anchored cell wall proteins (GPI-CWPs); the corresponding genes contain long tandem repeat regions. These repeat regions resulted in assembly errors in the reference genome. Here, we performed a de novo assembly of the C. glabrata type strain CBS138 using long single-molecule real-time reads, with short read sequences (Illumina) for refinement, and constructed telomere-to-telomere assemblies of all 13 chromosomes. Our assembly has excellent agreement overall with the current reference genome, but we made substantial corrections within tandem repeat regions. Specifically, we removed 62 genes of which 45 were scrambled due to misassembly in the reference. We annotated 31 novel ORFs of which 24 ORFs are GPI-CWPs. In addition, we corrected the tandem repeat structure of an additional 21 genes. Our corrections to the genome were substantial, with the length of new genes and tandem repeat corrections amounting to approximately 3.8% of the ORFeome length. As most corrections were within the coding regions of GPI-CWP genes, our genome assembly establishes a high-quality reference set of genes and repeat structures for the functional analysis of these cell surface proteins.

Genome-wide consequences of deleting any single gene.

Loss or duplication of chromosome segments can lead to further genomic changes associated with cancer. However, it is not known whether only a select subset of genes is responsible for driving further changes. To determine whether perturbation of any given gene in a genome suffices to drive subsequent genetic changes, we analyzed the yeast knockout collection for secondary mutations of functional consequence. Unlike wild-type, most gene knockout strains were found to have one additional mutant gene affecting nutrient responses and/or heat-stress-induced cell death. Moreover, independent knockouts of the same gene often evolved mutations in the same secondary gene. Genome sequencing identified acquired mutations in several human tumor suppressor homologs. Thus, mutation of any single gene may cause a genomic imbalance, with consequences sufficient to drive adaptive genetic changes. This complicates genetic analyses but is a logical consequence of losing a functional unit originally acquired under pressure during evolution.

Multi-ancestry fine mapping of interferon lambda and the outcome of acute hepatitis C virus infection.

Clearance of acute infection with hepatitis C virus (HCV) is associated with the chr19q13.13 region containing the rs368234815 (TT/ΔG) polymorphism. We fine-mapped this region to detect possible causal variants that may contribute to HCV clearance. First, we performed sequencing of IFNL1-IFNL4 region in 64 individuals sampled according to rs368234815 genotype: TT/clearance (N = 16) and ΔG/persistent (N = 15) (genotype-outcome concordant) or TT/persistent (N = 19) and ΔG/clearance (N = 14) (discordant). 25 SNPs had a difference in counts of alternative allele >5 between clearance and persistence individuals. Then, we evaluated those markers in an association analysis of HCV clearance conditioning on rs368234815 in two groups of European (692 clearance/1 025 persistence) and African ancestry (320 clearance/1 515 persistence) individuals. 10/25 variants were associated (P < 0.05) in the conditioned analysis leaded by rs4803221 (P value = 4.9 × 10-04) and rs8099917 (P value = 5.5 × 10-04). In the European ancestry group, individuals with the haplotype rs368234815ΔG/rs4803221C were 1.7× more likely to clear than those with the rs368234815ΔG/rs4803221G haplotype (P value = 3.6 × 10-05). For another nearby SNP, the haplotype of rs368234815ΔG/rs8099917T was associated with HCV clearance compared to rs368234815ΔG/rs8099917G (OR: 1.6, P value = 1.8 × 10-04). We identified four possible causal variants: rs368234815, rs12982533, rs10612351 and rs4803221. Our results suggest a main signal of association represented by rs368234815, with contributions from rs4803221, and/or nearby SNPs including rs8099917.

The Terabase Search Engine: a large-scale relational database of short-read sequences.

MOTIVATION: DNA sequencing archives have grown to enormous scales in recent years, and thousands of human genomes have already been sequenced. The size of these data sets has made searching the raw read data infeasible without high-performance data-query technology. Additionally, it is challenging to search a repository of short-read data using relational logic and to apply that logic across samples from multiple whole-genome sequencing samples. RESULTS: We have built a compact, efficiently-indexed database that contains the raw read data for over 250 human genomes, encompassing trillions of bases of DNA, and that allows users to search these data in real-time. The Terabase Search Engine enables retrieval from this database of all the reads for any genomic location in a matter of seconds. Users can search using a range of positions or a specific sequence that is aligned to the genome on the fly. AVAILABILITY AND IMPLEMENTATION: Public access to the Terabase Search Engine database is available at http://tse.idies.jhu.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MYC drives overexpression of telomerase RNA (hTR/TERC) in prostate cancer.

Telomerase consists of at least two essential elements, an RNA component hTR or TERC that contains the template for telomere DNA addition and a catalytic reverse transcriptase (TERT). While expression of TERT has been considered the key rate-limiting component for telomerase activity, increasing evidence suggests an important role for the regulation of TERC in telomere maintenance and perhaps other functions in human cancer. By using three orthogonal methods including RNAseq, RT-qPCR, and an analytically validated chromogenic RNA in situ hybridization assay, we report consistent overexpression of TERC in prostate cancer. This overexpression occurs at the precursor stage (e.g. high-grade prostatic intraepithelial neoplasia or PIN) and persists throughout all stages of disease progression. Levels of TERC correlate with levels of MYC (a known driver of prostate cancer) in clinical samples and we also show the following: forced reductions of MYC result in decreased TERC levels in eight cancer cell lines (prostate, lung, breast, and colorectal); forced overexpression of MYC in PCa cell lines, and in the mouse prostate, results in increased TERC levels; human TERC promoter activity is decreased after MYC silencing; and MYC occupies the TERC locus as assessed by chromatin immunoprecipitation (ChIP). Finally, we show that knockdown of TERC by siRNA results in reduced proliferation of prostate cancer cell lines. These studies indicate that TERC is consistently overexpressed in all stages of prostatic adenocarcinoma and that its expression is regulated by MYC. These findings nominate TERC as a novel prostate cancer biomarker and therapeutic target. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

ETS2 is a prostate basal cell marker and is highly expressed in prostate cancers aberrantly expressing p63.

BACKGROUND: Rare prostate carcinomas aberrantly express p63 and have an immunophenotype intermediate between basal and luminal cells. Here, we performed gene expression profiling on p63-expressing prostatic carcinomas and compared them to usual-type adenocarcinoma. We identify ETS2 as highly expressed in p63-expressing prostatic carcinomas and benign prostate basal cells, with lower expression in luminal cells and primary usual-type adenocarcinomas. METHODS: A total of 8 p63-expressing prostate carcinomas at radical prostatectomy were compared to 358 usual-type adenocarcinomas by gene expression profiling performed on formalin fixed paraffin embedded tumor tissue using Affymetrix 1.0 ST microarrays. Correlation between differentially expressed genes and TP63 expression was performed in 5239 prostate adenocarcinomas available in the Decipher GRID. For validation, ETS2 in situ hybridization was performed on 19 p63-expressing prostate carcinomas and 30 usual-type adenocarcinomas arrayed on tissue microarrays (TMA). RESULTS: By gene expression, p63-expressing prostate carcinomas showed low cell cycle activity and low Decipher prognostic scores, but were predicted to have high Gleason grade compared to usual-type adenocarcinomas by gene expression signatures and morphology. Among the genes over-expressed in p63-expressing carcinoma relative to usual-type adenocarcinoma were known p63-regulated genes, along with ETS2, an ETS family member previously implicated as a prostate cancer tumor suppressor gene. Across several cohorts of prostate samples, ETS2 gene expression was correlated with TP63 expression and was significantly higher in benign prostate compared to usual-type adenocarcinoma. By in situ hybridization, ETS2 gene expression was high in benign basal cells, and low to undetectable in benign luminal cells or usual-type adenocarcinoma. In contrast, ETS2 was highly expressed in 95% (18/19) of p63-expressing prostate carcinomas. CONCLUSIONS: ETS2 is a predominantly basally-expressed gene in the prostate, with low expression in usual-type adenocarcinoma and high expression in p63-expressing carcinomas. Given this pattern, the significance of ETS2 loss by deletion or mutation in usual-type adenocarcinomas is uncertain.

A primary breast cancer with distinct foci of estrogen receptor-alpha positive and negative cells derived from the same clonal origin as revealed by whole exome sequencing.

BACKGROUND/PURPOSE: Tumor heterogeneity is a now well-recognized phenomenon that can affect the classification, prognosis and treatment of human cancers. Heterogeneity is often described in primary breast cancers based upon histologic subtypes, hormone- and HER2-receptor status, and immunolabeling for various markers, which can be seen within a single tumor as mixed cellular populations, or as separate discrete foci. EXPERIMENTAL DESIGN/METHODS: Here, we present a case report of a patient's primary breast cancer that had two separate but adjacent histologic components, one that was estrogen receptor (ER) positive, and the other ER negative. Each component was subjected to whole exome sequencing and compared for gene identity to determine clonal origin. RESULTS: Using prior bioinformatic tools, we demonstrated that both the ER positive and negative components shared many variants, including passenger and driver alterations. Copy number variations also supported the two components were derived from a single common clone. CONCLUSIONS: These analyses strongly suggest that the two ER components of this patient's breast cancer were derived from the same clonal origin. Our results have implications for the evolution of breast cancers with mixed histologies, and how they might be best managed for optimal therapy.

StereoGene: rapid estimation of genome-wide correlation of continuous or interval feature data.

MOTIVATION: Genomics features with similar genome-wide distributions are generally hypothesized to be functionally related, for example, colocalization of histones and transcription start sites indicate chromatin regulation of transcription factor activity. Therefore, statistical algorithms to perform spatial, genome-wide correlation among genomic features are required. RESULTS: Here, we propose a method, StereoGene, that rapidly estimates genome-wide correlation among pairs of genomic features. These features may represent high-throughput data mapped to reference genome or sets of genomic annotations in that reference genome. StereoGene enables correlation of continuous data directly, avoiding the data binarization and subsequent data loss. Correlations are computed among neighboring genomic positions using kernel correlation. Representing the correlation as a function of the genome position, StereoGene outputs the local correlation track as part of the analysis. StereoGene also accounts for confounders such as input DNA by partial correlation. We apply our method to numerous comparisons of ChIP-Seq datasets from the Human Epigenome Atlas and FANTOM CAGE to demonstrate its wide applicability. We observe the changes in the correlation between epigenomic features across developmental trajectories of several tissue types consistent with known biology and find a novel spatial correlation of CAGE clusters with donor splice sites and with poly(A) sites. These analyses provide examples for the broad applicability of StereoGene for regulatory genomics. AVAILABILITY AND IMPLEMENTATION: The StereoGene C ++ source code, program documentation, Galaxy integration scripts and examples are available from the project homepage http://stereogene.bioinf.fbb.msu.ru/. CONTACT: favorov@sensi.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Characterization of functionally active gene fusions in human papillomavirus related oropharyngeal squamous cell carcinoma.

The Cancer Genome Atlas (TCGA) sequencing analysis of head and neck squamous cell carcinoma (HNSCC) recently reported on gene fusions, however, few human papillomavirus (HPV) positive samples were included, and the functional relevance of identified fusions was not explored. We therefore performed an independent analysis of gene fusions in HPV-positive oropharyngeal SCC (OPSCC). RNA sequencing was performed on 47 HPV-positive OPSCC primary tumors and 25 normal mucosal samples from cancer unaffected controls on an Illumina TruSeq platform. MapSplice2 was used for alignment and identification of fusion candidates. Putative fusions with less than five spanning reads, detected in normal tissues, or that mapped to the same gene were filtered out. Selected fusions were validated by RT-PCR and Sanger sequencing. Within 47 HPV-positive OPSCC tumors, 282 gene fusions were identified. Most fusions (85.1%) occurred in a single tumor, and the remaining fusions recurred in 2-16 tumors. Gene fusions were associated with significant up regulation of 16 genes (including EGFR and ERBB4) and down regulation of four genes (PTPRT, ZNF750, DLG2, SLCO5A1). Expression of these genes followed similar patterns of up regulation and down regulation in tumors without these fusions compared to normal tissue. Five of six gene fusions selected for validation were confirmed through RT-PCR and sequencing. This integrative analysis provides a method of prioritizing functionally relevant gene fusions that may be expanded to other tumor types. These results demonstrate that gene fusions may be one mechanism by which functionally relevant genes are altered in HPV-positive OPSCC.

Controlled insertional mutagenesis using a LINE-1 (ORFeus) gene-trap mouse model.

A codon-optimized mouse LINE-1 element, ORFeus, exhibits dramatically higher retrotransposition frequencies compared with its native long interspersed element 1 counterpart. To establish a retrotransposon-mediated mouse model with regulatable and potent mutagenic capabilities, we generated a tetracycline (tet)-regulated ORFeus element harboring a gene-trap cassette. Here, we show that mice expressing tet-ORFeus broadly exhibit robust retrotransposition in somatic tissues when treated with doxycycline. Consistent with a significant mutagenic burden, we observed a reduced number of double transgenic animals when treated with high-level doxycycline during embryogenesis. Transgene induction in skin resulted in a white spotting phenotype due to somatic ORFeus-mediated mutations that likely disrupt melanocyte development. The data suggest a high level of transposition in melanocyte precursors and consequent mutation of genes important for melanoblast proliferation, differentiation, or migration. These findings reveal the utility of a retrotransposon-based mutagenesis system as an alternative to existing DNA transposon systems. Moreover, breeding these mice to different tet-transactivator/reversible tet-transactivator lines supports broad functionality of tet-ORFeus because of the potential for dose-dependent, tissue-specific, and temporal-specific mutagenesis.

piggyBac transposase tools for genome engineering.

The transposon piggyBac is being used increasingly for genetic studies. Here, we describe modified versions of piggyBac transposase that have potentially wide-ranging applications, such as reversible transgenesis and modified targeting of insertions. piggyBac is distinguished by its ability to excise precisely, restoring the donor site to its pretransposon state. This characteristic makes piggyBac useful for reversible transgenesis, a potentially valuable feature when generating induced pluripotent stem cells without permanent alterations to genomic sequence. To avoid further genome modification following piggyBac excision by reintegration, we generated an excision competent/integration defective (Exc(+)Int(-)) transposase. Our findings also suggest the position of a target DNA-transposase interaction. Another goal of genome engineering is to develop reagents that can guide transgenes to preferred genomic regions. Others have shown that piggyBac transposase can be active when fused to a heterologous DNA-binding domain. An Exc(+)Int(-) transposase, the intrinsic targeting of which is defective, might also be a useful intermediate in generating a transposase whose integration activity could be rescued and redirected by fusion to a site-specific DNA-binding domain. We show that fusion to two designed zinc finger proteins rescued the Int(-) phenotype. Successful guided transgene integration into genomic DNA would have broad applications to gene therapy and molecular genetics. Thus, an Exc(+)Int(-) transposase is a potentially useful reagent for genome engineering and provides insight into the mechanism of transposase-target DNA interaction.

Transflip mutations produce deletions in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is driven by the inactivation of the tumor suppressor genes (TSGs), CDKN2A (P16) and SMAD4 (DPC4), commonly by homozygous deletions (HDs). Using a combination of high density single-nucleotide polymorphism (SNP) microarray and whole genome sequencing (WGS), we fine-mapped novel breakpoints surrounding deletions of CDKN2A and SMAD4 and characterized them by their underlying structural variants (SVs). Only one third of CDKN2A and SMAD4 deletions (6 of 18) were simple interstitial deletions, rather, the majority of deletions were caused by complex rearrangements, specifically, a translocation on one side of the TSG in combination with an inversion on the other side. We designate these as "TransFlip" mutations. Characteristics of TransFlip mutations are: (1) a propensity to target the TSGs CDKN2A and SMAD4 (P < 0.005), (2) not present in the germline of the examined samples, (3) non-recurrent breakpoints, (4) relatively small (47 bp to 3.4 kb) inversions, (5) inversions can be either telomeric or centromeric to the TSG, and (6) non-reciprocal, and non-recurrent translocations. TransFlip mutations are novel complex genomic rearrangements with unique breakpoint signatures in pancreatic cancer. We hypothesize that they are a common but poorly understood mechanism of TSG inactivation in human cancer. © 2015 Wiley Periodicals, Inc.

Retrotransposon Ty1 integration targets specifically positioned asymmetric nucleosomal DNA segments in tRNA hotspots.

The Saccharomyces cerevisiae genome contains about 35 copies of dispersed retrotransposons called Ty1 elements. Ty1 elements target regions upstream of tRNA genes and other Pol III-transcribed genes when retrotransposing to new sites. We used deep sequencing of Ty1-flanking sequence amplicons to characterize Ty1 integration. Surprisingly, some insertions were found in mitochondrial DNA sequences, presumably reflecting insertion into mitochondrial DNA segments that had migrated to the nucleus. The overwhelming majority of insertions were associated with the 5' regions of Pol III transcribed genes; alignment of Ty1 insertion sites revealed a strong sequence motif centered on but extending beyond the target site duplication. A strong sequence-independent preference for nucleosomal integration sites was observed, in distinction to the preferences of the Hermes DNA transposon engineered to jump in yeast and the Tf1 retrotransposon of Schizosaccharomyces pombe, both of which prefer nucleosome free regions. Remarkably, an exquisitely specific relationship between Ty1 integration and nucleosomal position was revealed by alignment of hotspot Ty1 insertion position regions to peak nucleosome positions, geographically implicating nucleosomal DNA segments at specific positions on the nucleosome lateral surface as targets, near the "bottom" of the nucleosome. The specificity is observed in the three tRNA 5'-proximal nucleosomes, with insertion frequency dropping off sharply 5' of the tRNA gene. The sites are disposed asymmetrically on the nucleosome relative to its dyad axis, ruling out several simple molecular models for Ty1 targeting, and instead suggesting association with a dynamic or directional process such as nucleosome remodeling associated with these regions.