Exome and targeted sequencing studies have identified potential driver mutations for a variety of tumor types. Cutaneous squamous cell carcinoma (cSCC) is one of the most highly mutated cancers but is typically correlated with high survival rates and low rates of metastasis. Nevertheless, metastatic cSCC is a significant health threat; up to 8800 individuals are estimated to die yearly from this disease. As it is difficult to predict which cSCCs are more likely to metastasize, and because there are no targeted therapies specifically designated for metastatic cSCC, we performed exome and targeted sequencing of 18 metastatic and 10 primary cSCCs to identify mutations with potential therapeutic benefit. Genes previously shown to be mutated in primary cSCC as well as aggressive tumors such as TP53 and NOTCH pathway genes were mutated at high rates in both metastatic and primary tumors. We compared our results to published sequencing results of an additional 223 primary tumors and 68 aggressive cSCCs. We identified several genes showing higher mutation frequencies in metastatic cSCC relative to primary tumors including the chromatin remodeling gene KMT2D, AKT/PI3K pathway gene PIK3CG, and the classic skin tumor suppressor TP53. These studies uncover potential pathways important in metastatic cSCC that may lead to new therapeutic strategies and understanding of the biology leading to more aggressive tumor behavior.

The statistical power of genome-wide association (GWA) studies to detect risk alleles for human diseases is limited by the unfavorable ratio of SNPs to study subjects. This multiple testing problem can be surmounted with very large population sizes when common alleles of large effects give rise to disease status. However, GWA approaches fall short when many rare alleles may give rise to a common disease, or when the number of subjects that can be recruited is limited. Here, we demonstrate that this multiple testing problem can be overcome by a comparative genomics approach in which an initial genome-wide screen in a genetically amenable model organism is used to identify human orthologues that may harbor risk alleles for adult-onset primary open angle glaucoma (POAG). Glaucoma is a major cause of blindness, which affects over 60 million people worldwide. Several genes have been associated with juvenile onset glaucoma, but genetic factors that predispose to adult onset primary open angle glaucoma (POAG) remain largely unknown. Previous genome-wide analysis in a Drosophila ocular hypertension model identified transcripts with altered regulation and showed induction of the unfolded protein response (UPR) upon overexpression of transgenic human glaucoma-associated myocilin (MYOC). We selected 16 orthologous genes with 62 polymorphic markers and identified in two independent human populations two genes of the UPR that harbor POAG risk alleles, BIRC6 and PDIA5. Thus, effectiveness of the UPR in response to accumulation of misfolded or aggregated proteins may contribute to the pathogenesis of POAG and provide targets for early therapeutic intervention.

The extent to which heritable genetic variants can affect tumor development has yet to be fully elucidated. Tumor selection of single nucleotide polymorphism (SNP) risk alleles, a phenomenon called preferential allelic imbalance (PAI), has been demonstrated in some cancer types. We developed a novel application of digital PCR termed Somatic Mutation Allelic Ratio Test using Droplet Digital PCR (SMART-ddPCR) for accurate assessment of tumor PAI, and have applied this method to test the hypothesis that heritable SNPs associated with childhood acute lymphoblastic leukemia (ALL) may demonstrate tumor PAI. These SNPs are located at CDKN2A (rs3731217) and IKZF1 (rs4132601), genes frequently lost in ALL, and at CEBPE (rs2239633), ARID5B (rs7089424), PIP4K2A (rs10764338), and GATA3 (rs3824662), genes located on chromosomes gained in high-hyperdiploid ALL. We established thresholds of AI using constitutional DNA from SNP heterozygotes, and subsequently measured allelic copy number in tumor DNA from 19-142 heterozygote samples per SNP locus. We did not find significant tumor PAI at these loci, though CDKN2A and IKZF1 SNPs showed a trend towards preferential selection of the risk allele (p = 0.17 and p = 0.23, respectively). Using a genomic copy number control ddPCR assay, we investigated somatic copy number alterations (SCNA) underlying AI at CDKN2A and IKZF1, revealing a complex range of alterations including homozygous and hemizygous deletions and copy-neutral loss of heterozygosity, with varying degrees of clonality. Copy number estimates from ddPCR showed high agreement with those from multiplex ligation-dependent probe amplification (MLPA) assays. We demonstrate that SMART-ddPCR is a highly accurate method for investigation of tumor PAI and for assessment of the somatic alterations underlying AI. Furthermore, analysis of publicly available data from The Cancer Genome Atlas identified 16 recurrent SCNA loci that contain heritable cancer risk SNPs associated with a matching tumor type, and which represent candidate PAI regions warranting further investigation.

Background

Mus spretus diverged from Mus musculus over one million years ago. These mice are genetically and phenotypically divergent. Despite the value of utilizing M. musculus and M. spretus for quantitative trait locus (QTL) mapping, relatively little genomic information on M. spretus exists, and most of the available sequence and polymorphic data is for one strain of M. spretus, Spret/Ei. In previous work, we mapped fifteen loci for skin cancer susceptibility using four different M. spretus by M. musculus F1 backcrosses. One locus, skin tumor susceptibility 5 (Skts5) on chromosome 12, shows strong linkage in one cross.

Results

To identify potential candidate genes for Skts5, we sequenced 65 named and unnamed genes and coding elements mapping to the peak linkage area in outbred spretus, Spret/EiJ, FVB/NJ, and NIH/Ola. We identified polymorphisms in 62 of 65 genes including 122 amino acid substitutions. To look for polymorphisms consistent with the linkage data, we sequenced exons with amino acid polymorphisms in two additional M. spretus strains and one additional M. musculus strain generating 40.1 kb of sequence data. Eight candidate variants were identified that fit with the linkage data. To determine the degree of variation across M. spretus, we conducted phylogenetic analyses. The relatedness of the M. spretus strains at this locus is consistent with the proximity of region of ascertainment of the ancestral mice.

Conclusion

Our analyses suggest that, if Skts5 on chromosome 12 is representative of other regions in the genome, then published genomic data for Spret/EiJ are likely to be of high utility for genomic studies in other M. spretus strains.

Although genetics is the most significant known determinant of human intelligence, specific gene contributions remain largely unknown. To accelerate understanding in this area, we have taken a new approach by studying the relationship between quantitative gene expression and intelligence in a cohort of 65 patients with Williams Syndrome (WS), a neurodevelopmental disorder caused by a 1.5 Mb deletion on chromosome 7q11.23. We find that variation in the transcript levels of the brain gene STX1A correlates significantly with intelligence in WS patients measured by principal component analysis (PCA) of standardized WAIS-R subtests, r = 0.40 (Pearson correlation, Bonferroni corrected p-value = 0.007), accounting for 15.6% of the cognitive variation. These results suggest that syntaxin 1A, a neuronal regulator of presynaptic vesicle release, may play a role in WS and be a component of the cellular pathway determining human intelligence.

HIV infection disturbs the central nervous system (CNS) through inflammation and glial activation. Evidence suggests roles for microRNA (miRNA) in host defense and neuronal homeostasis, though little is known about miRNAs' role in HIV CNS infection. MiRNAs are non-coding RNAs that regulate gene translation through post-transcriptional mechanisms. Messenger-RNA profiling alone is insufficient to elucidate the dynamic dance of molecular expression of the genome. We sought to clarify RNA alterations in the frontal cortex (FC) of HIV-infected individuals and those concurrently infected and diagnosed with major depressive disorder (MDD). This report is the first published study of large-scale miRNA profiling from human HIV-infected FC. The goals of this study were to: 1. Identify changes in miRNA expression that occurred in the frontal cortex (FC) of HIV individuals, 2. Determine whether miRNA expression profiles of the FC could differentiate HIV from HIV/MDD, and 3. Adapt a method to meaningfully integrate gene expression data and miRNA expression data in clinical samples. We isolated RNA from the FC (n = 3) of three separate groups (uninfected controls, HIV, and HIV/MDD) and then pooled the RNA within each group for use in large-scale miRNA profiling. RNA from HIV and HIV/MDD patients (n = 4 per group) were also used for non-pooled mRNA analysis on Affymetrix U133 Plus 2.0 arrays. We then utilized a method for integrating the two datasets in a Target Bias Analysis. We found miRNAs of three types: A) Those with many dysregulated mRNA targets of less stringent statistical significance, B) Fewer dysregulated target-genes of highly stringent statistical significance, and C) unclear bias. In HIV/MDD, more miRNAs were downregulated than in HIV alone. Specific miRNA families at targeted chromosomal loci were dysregulated. The dysregulated miRNAs clustered on Chromosomes 14, 17, 19, and X. A small subset of dysregulated genes had many 3' untranslated region (3'UTR) target-sites for dysregulated miRNAs. We provide evidence that certain miRNAs serve as key elements in gene regulatory networks in HIV-infected FC and may be implicated in neurobehavioral disorder. Finally, our data indicates that some genes may serve as hubs of miRNA activity.

More than 3.5 million nonmelanoma skin cancers were treated in 2006; of these 700,000 were cutaneous squamous cell carcinomas (cSCCs). Despite clear environmental causes for cSCC, studies also suggest genetic risk factors. A cSCC susceptibility locus, Skts5, was identified on mouse chromosome 12 by linkage analysis. The orthologous locus to Skts5 in humans maps to 7p21 and 7q31. These loci show copy number increases in ∼10% of cSCC tumors. Here, we show that an additional 15-22% of tumors exhibit copy-neutral loss of heterozygosity. Furthermore, our previous data identified microsatellite markers on 7p21 and 7q31 that demonstrate preferential allelic imbalance (PAI) in cSCC tumors. On the basis of these results, we hypothesized that the human orthologous locus to Skts5 would house a gene important in human cSCC development and that tumors would demonstrate allele-specific somatic alterations. To test this hypothesis, we performed quantitative genotyping of 108 single nucleotide polymorphisms (SNPs) mapping to candidate genes at human SKTS5 in paired normal and tumor DNAs. Nine SNPs in HDAC9 (rs801540, rs1178108, rs1178112, rs1726610, rs10243618, rs11764116, rs1178355, rs10269422 and rs12540872) showed PAI in tumors. These data suggest that HDAC9 variants may be selected for during cSCC tumorigenesis.

PALB2 has been identified as a breast and pancreatic cancer susceptibility gene. Utilizing a targeted sequencing approach, we discovered two novel germline missense PALB2 variants c.191C>T and c.311C>T, encoding p.Ser64Leu and p.Pro104Leu, respectively, in individuals in a pancreatic cancer registry. No missense PALB2 variants from familial pancreatic cancer patients, and few PALB2 variants overall, have been functionally characterized. Given the known role of PALB2, we tested the impact of p.Ser64Leu and p.Pro104Leu variants on DNA damage responses. Neither p.Ser64Leu nor p.Pro104Leu have clear effects on interactions with BRCA1 and KEAP1, which are mediated by adjacent motifs in PALB2. However, both variants are associated with defective recruitment of PALB2, and the RAD51 recombinase downstream, to DNA damage foci. Furthermore, p.Ser64Leu and p.Pro104Leu both largely compromise DNA double-strand break-initiated homologous recombination, and confer increased cellular sensitivity to ionizing radiation (IR) and the poly (ADP-ribose) polymerase (PARP) inhibitor Olaparib. Taken together, our results represent the first demonstration of functionally deleterious PALB2 missense variants associated with familial pancreatic cancer and of deleterious variants in the N-terminus outside of the coiled-coil domain. Furthermore, our results suggest the possibility of personalized treatments, using IR or PARP inhibitor, of pancreatic and other cancers that carry a deleterious PALB2 variant.

Background

Hispanic/Latino populations are a genetically admixed and heterogeneous group, with variable fractions of European, Indigenous American and African ancestries. The molecular profile of breast cancer has been widely described in non-Hispanic Whites but equivalent knowledge is lacking in Hispanic/Latinas. We have previously reported that the most prevalent breast cancer intrinsic subtype in Colombian women was Luminal B as defined by St. Gallen 2013 criteria. In this study we explored ancestry-associated differences in molecular profiles of Luminal B tumors among these highly admixed women.

Methods

We performed whole-transcriptome RNA-seq analysis in 42 Luminal tumors (21 Luminal A and 21 Luminal B) from Colombian women. Genetic ancestry was estimated from a panel of 80 ancestry-informative markers (AIM). We categorized patients according to Luminal subtype and to the proportion of European and Indigenous American ancestry and performed differential expression analysis comparing Luminal B against Luminal A tumors according to the assigned ancestry groups.

Results

We found 5 genes potentially modulated by genetic ancestry: ERBB2 (log2FC = 2.367, padj<0.01), GRB7 (log2FC = 2.327, padj<0.01), GSDMB (log2FC = 1.723, padj<0.01, MIEN1 (log2FC = 2.195, padj<0.01 and ONECUT2 (log2FC = 2.204, padj<0.01). In the replication set we found a statistical significant association between ERBB2 expression with Indigenous American ancestry (p = 0.02, B = 3.11). This association was not biased by the distribution of HER2+ tumors among the groups analyzed.

Conclusions

Our results suggest that genetic ancestry in Hispanic/Latina women might modify ERBB2 gene expression in Luminal tumors. Further analyses are needed to confirm these findings and explore their prognostic value.

Background

Hispanic/Latino populations are a genetically admixed and heterogeneous group, with variable fractions of European, Indigenous American and African ancestries. The molecular profile of breast cancer has been widely described in non-Hispanic Whites but equivalent knowledge is lacking in Hispanic/Latinas. We have previously reported that the most prevalent breast cancer intrinsic subtype in Colombian women was Luminal B as defined by St. Gallen 2013 criteria. In this study we explored ancestry-associated differences in molecular profiles of Luminal B tumors among these highly admixed women.

Methods

We performed whole-transcriptome RNA-seq analysis in 42 Luminal tumors (21 Luminal A and 21 Luminal B) from Colombian women. Genetic ancestry was estimated from a panel of 80 ancestry-informative markers (AIM). We categorized patients according to Luminal subtype and to the proportion of European and Indigenous American ancestry and performed differential expression analysis comparing Luminal B against Luminal A tumors according to the assigned ancestry groups.

Results

We found 5 genes potentially modulated by genetic ancestry: ERBB2 (log2FC = 2.367, padj<0.01), GRB7 (log2FC = 2.327, padj<0.01), GSDMB (log2FC = 1.723, padj<0.01, MIEN1 (log2FC = 2.195, padj<0.01 and ONECUT2 (log2FC = 2.204, padj<0.01). In the replication set we found a statistical significant association between ERBB2 expression with Indigenous American ancestry (p = 0.02, B = 3.11). This association was not biased by the distribution of HER2+ tumors among the groups analyzed.

Conclusions

Our results suggest that genetic ancestry in Hispanic/Latina women might modify ERBB2 gene expression in Luminal tumors. Further analyses are needed to confirm these findings and explore their prognostic value.