Search

Cystic fibrosis (CF) is a life‐limiting autosomal recessive genetic disease caused by variants in the CFTR gene, most commonly by the [F508del] variant. Although CF is a classical Mendelian disease, genetic variants in several modifier genes have been associated with variation of the clinical phenotype for pulmonary and gastrointestinal function and urogenital development. We hypothesized that whole genome sequencing of a well‐phenotyped CF populations might identify novel variants in known, or hitherto unknown, modifier genes. Whole genome sequencing was performed on the Illumina HiSeq X platform for 98 clinically diagnosed cystic fibrosis patient samples from the Adult CF Clinic at the University of California San Diego (UCSD). We compared protein‐coding, non‐silent variants genome wide between CFTR [F508del] homozygotes vs CFTR compound heterozygotes. Based on a single variant score test, we found 3 SNPs in common variants (MAF >5%) that occurred at significantly different rates between homozygous [F508del]CFTR and compound heterozygous [F508del]CFTR patients. The 3 SNPs were all located in one gene on chromosome 2: Tensin 1 (TNS1: rs3796028; rs2571445: and rs918949). We observed significantly lower BMIs in homozygous [F508del]CFTR patients who were also homozygous for Tensin 1 rs918949 (T/T) (p = 0.023) or rs2571445 (G/G) (p = 0.02) variants. The Tensin 1 gene is thus a potential modifier gene for low BMI in CF patients homozygous for the [F508del]CFTR variant., We hypothesized that whole genome sequencing of a well‐phenotyped CF populations might identify novel variants in known, or hitherto unknown, modifier genes. We compared protein‐coding, non‐silent variants genomewide between CFTR [F508del] homozygotes vs CFTR compound heterozygotes. We found that the Tensin 1 gene is a potential modifier gene for low BMI in CF patients homozygous for the [F508del]CFTR variant.

Background: Genomic evaluation of structural and molecular alterations is essential for the classification and risk stratification of patients with acute myeloid leukemia (AML). Cytogenetics (CG) has traditionally been used for the identification of large structural or copy number changes and molecular studies for small sequence changes. Rapid technical advances and decreased sequencing (seq) costs have made whole-genome sequencing (WGS) more feasible and recent efforts have suggested it can match the analytic performance of traditional CG while providing additional prognostic information. However, the use of WGS in the clinic is still in its infancy and questions remain regarding the utility of many factors including seq depth, structural events not identifiable by traditional CG, normal samples, AML enrichment, paired RNA-seq, and single-cell DNA-seq. Methods: Baseline bone marrow (BM) samples from 10 AML patients underwent evaluation using standard clinical CG and targeted next generation seq tests. Additionally, samples were extensively characterized using: 200X WGS (bulk and enriched BM), 50X germline WGS (buccal swab and remission blood), 50X RNA-seq, targeted error-corrected (EC) DNA-seq, and single-cell DNA-seq with oligonucleotide-conjugated antibodies (scDNA-seq+ab). Results: The AML patients analyzed had a mean age of 71 (30-81) and a range of BM blast percentage by clinical flow cytometry (0.4-33%). Standard clinical assessment revealed 14 CG abnormalities and 20 mutations across all patients. 50X WGS of the bulk BM utilizing previously published conditions (PMID:33704937) successfully identified 89% of the alterations. Increasing the WGS depth to 100X recovered 2 additional abnormalities, increasing the capture rate to 94%. Enrichment of AML using the primary cell surface marker prior to WGS achieved a 94% capture rate at only 50X depth, with a depth of only 10X required for the CG abnormalities. Integration of 50X WGS on buccal swab samples reduced the large number of structural variant calls and private mutations without the need for aggressive filtering conditions required for tumor-only samples. However, this also resulted in the removal of a few true mutations due to buccal contamination with immune cells. 50X WGS on remission blood samples from 2 patients mitigated this problem. Next, we sought to determine what additional information WGS could provide that was missed by standard CG (>5Mb) and seq (>5% variant allele fraction, VAF) analyses. Firstly, 3 of the deletions reported by CG were found to be complex translocation events, resulting in different losses and gains than anticipated. Most importantly, decreasing the CNV size threshold ( Targeted EC DNA-seq was used to detect of variants Finally, personalized scDNA-seq+ab provided further resolution of each patient's AML through confirmation of new structural and mutational findings, identification of mutations associated with clonal hematopoiesis versus AML, and definition of AML clonal structure (Figure 1). Conclusions: We confirmed that 50X WGS can recapitulate standard clinical CG and mutational profiling in AML patients while providing improved resolution. Deeper seq, AML pre-enrichment, and integrated RNA-seq further aided in identifying low level events. We also demonstrated that WGS can identify additional important structural events missed by CG. Finally, personalized scDNA-seq+ab was able to define clonal architecture important for prognostic and residual disease tracking. Together, these results provide a framework for future integration of WGS into the clinic and personalized patient care. Figure 1 Figure 1. Disclosures Sciambi: Mission Bio Inc.: Current Employment. Durruthy-Durruthy: Mission Bio Inc.: Current Employment. Hourigan: Sellas: Research Funding.

Mounting evidence suggests that copy number variations (CNVs) can contribute to cancer susceptibility. The main goal of this study was to evaluate the role of germline CNVs in melanoma predisposition in high-risk melanoma families. We used genome-wide tiling comparative genomic hybridization and single nucleotide polymorphism arrays to characterize CNVs in 335 individuals (240 melanoma cases) from American melanoma-prone families (22 with germline CDKN2A or CDK4 mutations). We found that the global burden of overall CNVs (or deletions or duplications separately) was not significantly associated with case-control or CDKN2A/CDK4 mutation status after accounting for the familial dependence. However, we identified several rare CNVs that either involved known melanoma genes (e.g., PARP1, CDKN2A) or cosegregated with melanoma (duplication on 10q23.23, 3p12.2 and deletions on 8q424.3, 2q22.1) in families without mutations in known melanoma high-risk genes. Some of these CNVs were correlated with expression changes in disrupted genes based on RNASeq data from a subset of melanoma cases included in the CNV study. These results suggest that rare cosegregating CNVs may influence melanoma susceptibility in some melanoma-prone families and genes found in our study warrant further evaluation in future genetic analyses of melanoma.

Castration Resistant Prostate Cancer (CRPC) is thought to be driven by a collaborative mechanism between TNFα/NFκB and TGFβ signaling, leading to inflammation, Epithelial-to-Mesenchymal-Transition (EMT), and metastasis. Initially, TGFβ is a tumor suppressor, but in advanced metastatic disease it switches to being a tumor promoter. TGFBR2 may play a critical role in this collaboration, as its expression is driven by NFκB and it is the primary receptor for TGFβ. We have previously reported that the cardenolide drug digitoxin blocks TNFα/NFκB-driven proinflammatory signaling. We therefore hypothesized that digitoxin might break the collaborative process between NFκB and TGFβ by also inhibiting expression of TGFBR2. We therefore tested whether TGFβ-driven EMT and resulting metastases would be suppressed. Here we show, in vitro, that digitoxin inhibits NFκB-driven TGFBR2 expression, as well as Vimentin, while elevating E-cadherin expression. Digitoxin also significantly reduces HSPB1 mRNA and the HSPB1/RBFOX2 mRNA ratio in PC3 cells. In vivo, in a syngeneic, immune competent rat model of metastatic CRPC, we show that digitoxin also suppresses Tgfbr2 expression, as well as expression of other genes classically driven by NFκB, and of multiple EMT genes associated with metastasis. Concurrently, digitoxin suppresses tumor growth and metastasis in these animals, and prolongs survival. Gross tumor recurrence following tumor resection also appears prevented in ca 30% of cases. While the existence of a collaboration between NFκB and TGFβ to drive EMT and metastasis has previously been appreciated, we show here, for the first time, that chronic, low concentrations of digitoxin are able to block CRPC tumor progression, EMT and the ensuing metastatic disease.

Posttraumatic stress disorder is developed by exposure to a threatening and/or a horrifying event and characterized by the presence of anxiety, hyperarousal, avoidance, and sleep abnormality for a prolonged period of time. To elucidate the potential molecular mechanisms, we constructed a mouse model by electric foot shock followed by situational reminders and performed transcriptome analysis in brain tissues. The stressed mice acquired anxiety-like behavior after 2 weeks and exaggerated startle response after 4 weeks. Avoidance latency and freezing behavior were sustained up to 5 weeks post stress and abnormal static behavior was observed during the sleep period. RNA sequencing was performed in two of the emotional regulatory regions, anterior cingulate cortex and amygdala, at 2 and 5 weeks post stress. More than 1000 differentially expressed genes were identified at 2 weeks in both regions. The number of the regulated genes remained constant in amygdala at 5 weeks post stress, whereas those in anterior cingulate cortex were plummeted. Although synaptic remodeling and endocrine system were the most enriched signaling pathways in both anterior cingulate cortex and amygdala, the individual gene expression profile was regulated in a region- and time-dependent manner. In addition, several genes associated with PTSD involved in Hypothalamic-Pituitary-Adrenal axis were differentially regulated. These findings suggested that global gene expression profile was dynamically regulated in accordance with the disease development stage, and therefore targeting the distinct signaling molecules in different region and development stage might be critical for effective treatment to PTSD. Electronic supplementary material The online version of this article (10.1186/s13041-019-0449-0) contains supplementary material, which is available to authorized users.

Hodgkin lymphoma shows strong familial aggregation but no major susceptibility genes have been identified to date. The goal of this study was to identify high-penetrance variants using whole exome sequencing in 17 Hodgkin lymphoma prone families with three or more affected cases or obligate carriers (69 individuals), followed by targeted sequencing in an additional 48 smaller HL families (80 individuals). Alignment and variant calling were performed using standard methods. Dominantly segregating, rare, coding or potentially functional variants were further prioritized based on predicted deleteriousness, conservation, and potential importance in lymphoid malignancy pathways. We selected 23 genes for targeted sequencing. Only the p.A1065T variant in KDR (kinase insert domain receptor) also known as VEGFR2 (vascular endothelial growth factor receptor 2) was replicated in two independent Hodgkin lymphoma families. KDR is a type III receptor tyrosine kinase, the main mediator of vascular endothelial growth factor induced proliferation, survival, and migration. Its activity is associated with several diseases including lymphoma. Functional experiments have shown that p.A1065T, located in the activation loop, can promote constitutive autophosphorylation on tyrosine in the absence of vascular endothelial growth factor and that the kinase activity was abrogated after exposure to kinase inhibitors. A few other promising mutations were identified but appear to be "private". In conclusion, in the largest sequenced cohort of Hodgkin lymphoma families to date, we identified a causal mutation in the KDR gene. While independent validation is needed, this mutation may increase downstream tumor cell proliferation activity and might be a candidate for targeted therapy.

For unknown reasons, there is huge variability in risk conferred by different HPV types and, remarkably, strong differences even between closely related variant lineages within each type. HPV16 is a uniquely powerful carcinogenic type, causing approximately half of cervical cancer and most other HPV-related cancers. To permit the large-scale study of HPV genome variability and precancer/cancer, starting with HPV16 and cervical cancer, we developed a high-throughput next-generation sequencing (NGS) whole-genome method. We designed a custom HPV16 AmpliSeq™ panel that generated 47 overlapping amplicons covering 99% of the genome sequenced on the Ion Torrent Proton platform. After validating with Sanger, the current “gold standard” of sequencing, in 89 specimens with concordance of 99.9%, we used our NGS method and custom annotation pipeline to sequence 796 HPV16-positive exfoliated cervical cell specimens. The median completion rate per sample was 98.0%.Our method enabled us to discover novel SNPs, large contiguous deletions suggestive of viral integration (OR of 27.3, 95% CI 3.3–222, P=0.002), and the sensitive detection of variant lineage coinfections. This method represents an innovative high-throughput, ultra-deep coverage technique for HPV genomic sequencing, which, in turn, enables the investigation of the role of genetic variation in HPV epidemiology and carcinogenesis. Keywords: HPV16, HPV epidemiology, HPV genomics

Lung adenocarcinoma is a highly lethal tumor that displays extensive molecular heterogeneity of which deep characterization may drive therapeutic development and improve clinical outcomes. Through the Applied Proteogenomics Organizational Learning and Outcomes (APOLLO) research network, we utilized five molecular profiling technologies (DNA whole genome sequencing, RNA sequencing, total and phospho-proteomics by mass spectrometry, and reverse phase protein arrays [RPPA]) to characterize a longitudinally-annotated cohort of 87 lung adenocarcinomas. Through whole genome sequencing, we identified molecular signatures from patterns of somatic SNVs, indels, and large structural alterations that stratified tumors into three groups associated with patient smoke exposures. We also identified TP53, EGFR, KRAS, and STK11 as recurrently mutated genes, which together represent 80% of the cohort, in addition to genes mutated in smaller cohort subsets (e.g. RBM10), fusion genes, and pathogenic germline variants. To characterize tumor proteomes, we quantified >7,000 proteins and >10,000 phosphopeptides by mass spectrometry and >300 species by RPPA. Matched RNAs and proteins were typically positively correlated across samples (median ρ = 0.49). Through quantitative trait loci analyses, we identified genes whose RNA and protein expression levels were significantly modified by somatic mutations. We then classified tumors into RNA expression subtypes and found coordinated proteogenomic alterations and distinct clinical associations: terminal respiratory unit subtype – EGFR mutations and RNA/protein overexpression, acinar histology, non-smokers; proximal-proliferative subtype – STK11 mutations and RNA/protein underexpression, high smoking signature; proximal-inflammatory subtype – high tumor mutational burden. We also identified phospho-peptide signatures associated with these subtypes, including downregulation of CDK1/2 targets in terminal respiratory unit tumors. Protein co-expression network analysis discovered biologically-diverse pathway activities of the RNA expression subtypes. To interrogate somatic mutations in the context of molecular pathways, we projected DNA alterations onto known interaction networks and identified four subtypes with markedly distinct proteomic and microenvironment characteristics. Finally, several molecular characteristics were found to significantly predict patient outcomes, including RNA expression subtype classification against metastasis-free survival. Thus, our integrative, proteogenomic characterization of lung adenocarcinoma uncovered novel tumor biology and identified potential molecular markers for predicting patient outcomes. The views expressed in this abstract are solely of the authors and do not reflect the official policy of the Departments of Army/Navy/Air Force, Department of Defense, USUHS, HJF, or U.S. Government. Citation Format: Anthony R. Soltis, Nicholas W. Bateman, Thomas P. Conrads, Clifton L. Dalgard, Hai Hu, Teri J. Franks, Jianfang Liu, Daoud Meerzaman, Emanuel F. Petricoin, Qingrong Chen, Chunhua Yan, Xijun Zhang, Clesson E. Turner, The APOLLO Research Network, Craig D. Shriver, Christopher A. Moskaluk, Robert F. Browning, Matthew D. Wilkerson. Comprehensive proteogenomic analysis and classification of lung adenocarcinoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5893.

Background: DNA damage repair genes (DDRGs) play a critical role in protecting genome integrity and have been implicated in several cancer types. In the context of prostate cancer (CaP) emerging data provide potential role of this pathway in aggressive disease. It has also been recently demonstrated that PARP inhibitors can extend overall survival in metastatic patients with DDRG mutations. However, the association and therapeutic stratification based on inherited mutations of DDRGs remains to be defined in African American (AA) CaP patients. Our objective was to assess the frequency and association with disease aggressiveness of all known DDRGs in blood derived genomic DNAs of AA and Caucasian American (CA) CaP patients archived at the DOD Center for Prostate Disease Research and assess how this information can refine patient stratification for specific targeted therapeutic options. Method: Germline mutations in all DDRGs was evaluated by whole genome sequence (WGS) analysis of archived blood DNA samples from 600 CaP patients (300 AA and 300 CA) who underwent primary treatment at Walter Reed National Military Medical Center (WRNMMC) over the past 20 years. These patients had equal access Department of Defense healthcare system with up to 20 years of follow-up time. Following quality control steps to assess DNA quantity by Qubit assay and DNA quality by Bioanalyzer assay, DNA samples were used to generate PCR-free libraries for WGS using the NovaSeq (Illumina) platform. Out of the 600 libraries we generated, 14 dropped out, achieving a success rate of 97.6%. The successful libraries had an excellent quality based on DNA library metrics (yield and fragment length). Whole genome sequencing depth of the samples exceeded 37x on average and about 4 million SNPs were identified in the samples. Patient genotypes were projected onto principal components from reference populations. Sample ancestries were predicted by using the “Peddy” program, which uses a machine learning model trained on individuals of diverse ancestries from the 1000 Genomes Project reference panel. Due to mismatched ancestry, 33 samples were excluded from further analyses. An additional 17 samples were excluded due to higher than minimal noise level based on analyses using “ContEst” tool from Broad GATK package and Illumina noise percent values. Results: Interrogation of an inclusive DDRG set of 180 genes predicted to have non-silent effects on the protein sequence (e.g. missense, nonsense, frameshift) for germline mutations in this cohort revealed several known and novel mutations. The analysis has not been finalized yet. Mutation data will be assessed for association with the extensive available clinical and pathological data, including disease progression (metastasis), family history and African ancestry. Novel mutations, anticipated especially in the understudied AA context, will be analyzed for functional impact. Citation Format: Kevin Babcock, Xijun Zhang, Matthew Wilkerson, Clifton L. Dalgard, Shyh-Han Tan, Lakshmi Ravindranath, Yongmei Chen, Jennifer Cullen, Shiv Srivastava, Inger L. Rosner, Gyorgy Petrovics. Defining germline mutations of DNA damage repair genes in African American prostate cancer patients [abstract]. In: Proceedings of the Twelfth AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2019 Sep 20-23; San Francisco, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(6 Suppl_2):Abstract nr B075.

ObjectiveTo elucidate the genetic architecture of gene expression in pancreatic tissues.DesignWe performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison.ResultsWe identified 38 615 cis-eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis-eQTL (in 237 genes) in tumour-derived tissues (false discovery rate cis-eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10−8) and tumour-derived (p=8.3×10−5) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the ‘O’ mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO ‘O’ mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL.ConclusionsWe have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

Juvenile myelomonocytic leukemia (JMML) is a pediatric myeloproliferative neoplasm that arises from malignant transformation of the stem cell compartment and results in increased production of myeloid cells. Somatic and germline variants in CBL (Casitas B-lineage lymphoma proto-oncogene) have been associated with JMML. We report an incompletely penetrant CBL Y371C mutation discovered by whole-exome sequencing in three individuals with JMML in a large pedigree with 35 years of follow-up. The Y371 residue is highly evolutionarily conserved among CBL orthologs and paralogs. In silico bioinformatics prediction programs suggested that the Y371C mutation is highly deleterious. Protein structural modeling revealed that the Y371C mutation abrogated the ability of the CBL protein to adopt a conformation that is required for ubiquitination. Clinically, the three mutation-positive JMML individuals exhibited variable clinical courses; in two out of three, primary hematologic abnormalities persisted into adulthood with minimal clinical symptoms. The penetrance of the CBL Y371C mutation was 30 % for JMML and 40 % for all leukemia. Of the 8 mutation carriers in the family with available photographs, only one had significant dysmorphic features; we found no evidence of a clinical phenotype consistent with a “CBL syndrome”. Although CBL Y371C has been previously reported in familial JMML, we are the first group to follow a complete pedigree harboring this mutation for an extended period, revealing additional information about this variant’s penetrance, function and natural history.

Although CDKN2A is the most frequent high-risk melanoma susceptibility gene, the underlying genetic factors for most melanoma-prone families remain unknown. Using whole-exome sequencing, we identified a rare variant that arose as a founder mutation in the telomere shelterin gene POT1 (chromosome 7, g.124493086C>T; p.Ser270Asn) in five unrelated melanoma-prone families from Romagna, Italy. Carriers of this variant had increased telomere lengths and numbers of fragile telomeres, suggesting that this variant perturbs telomere maintenance. Two additional rare POT1 variants were identified in all cases sequenced in two separate Italian families, one variant per family, yielding a frequency for POT1 variants comparable to that for CDKN2A mutations in this population. These variants were not found in public databases or in 2,038 genotyped Italian controls. We also identified two rare recurrent POT1 variants in US and French familial melanoma cases. Our findings suggest that POT1 is a major susceptibility gene for familial melanoma in several populations.

Objective The objective of this study was to fine-map common pancreatic cancer susceptibility regions. Methods We conducted targeted Roche-454 resequencing across 428 kb in 3 genomic regions identified in genome-wide association studies (GWAS) of pancreatic cancer, on chromosomes 1q32.1, 5p15.33, and 13q22.1. Results An analytical pipeline for calling genotypes was developed using HapMap samples sequenced on chr5p15.33. Concordance to 1000 Genomes data for chr5p15.33 was greater than 96%. The concordance for chr1q32.1 and chr13q22.1 with pancreatic cancer GWAS data was greater than 99%. Between 9.2% and 19.0% of variants detected were not present in 1000 Genomes for the respective continental population. The majority of completely novel single-nucleotide polymorphisms (SNPs) were less common (minor allele frequency [MAF], ≤5%) or rare (MAF, ≤2%), illustrating the value of enlarging test sets for discovery of less common variants. Using the data set, we examined haplotype blocks across each region using a tag SNP analysis (r2 > 0.8 for MAF of ≥5%) and determined that at least 196, 243, and 63 SNPs are required for fine-mapping chr1q32.1, chr5p15.33, and chr13q22.1, respectively, in European populations. Conclusions We have characterized germline variation in 3 regions associated with pancreatic cancer risk and show that targeted resequencing leads to the discovery of novel variants and improves the completeness of germline sequence variants for fine-mapping GWAS susceptibility loci.

Germline inactivating mutations of isoform 4 of phosphodiesterase (PDE) 11A (coded by the PDE11A gene) have been associated with familial adrenocortical tumors and familial testicular cancer. Testicular tissue is unique in expressing all four isoforms of PDE11A. In a prior candidate gene study of 94 familial testicular germ cell tumor (TGCT) subjects, we identified a significant association between the presence of functionally abnormal variants in PDE11A and familial TGCT risk. To validate this novel observation, we sequenced the PDE11A coding region in 259 additional TGCT patients (both familial and sporadic) and 363 controls. We identified 55 PDE11A variants: 20 missense, four splice-site, two nonsense, seven synonymous, and 22 intronic. Ten missense variants were novel; nine occurred in transcript variant 4 and one in transcript variant 3. Five rare mutations (p.F258Y, p.G291R, p.V820M, p.R545X, and p.K568R) were present only in cases and were significantly more common in cases vs controls (P=0.0037). The latter two novel variants were functionally characterized and shown to be functionally inactivating, resulting in reduced PDE activity and increased cAMP levels. In further analysis of this cohort, we focused on white participants only to minimize confounding due to population stratification. This study builds upon our prior reports implicating PDE11A variants in familial TGCT, provides the first independent validation of those findings, extends that work to sporadic testicular cancer, demonstrates that these variants are uncommonly but reproducibly associated with TGCT, and refines our understanding regarding which specific inactivating PDE11A variants are most likely to be associated with TGCT risk.

BACKGROUND Genome-wide association studies of prostate cancer have identified single nucleotide polymorphism (SNP) markers in a region of chromosome 11q13.3 in men of European decent. A fine-mapping analysis with tag SNPs in the cancer genetic markers of susceptibility study identified three independent loci, marked by rs10896438, rs12793759, and rs10896449. This study further annotates common and uncommon variation across this region. METHODS A next generation resequence analysis of a 122.9-kb region of 11q13.3 (68,642,755-68,765,690) was conducted in 78 unrelated individuals of European background, 1 CEPH trio, and 1 YRI trio. RESULTS In total, 644 polymorphic loci were identified by our sequence analysis. Of these, 166 variants—118 SNPs and 48 insertion-deletion polymorphisms (indels)—were novel, namely not present in the 1000 Genomes or International HapMap Projects. We identified 22, 25, 6, and 4 variants strongly correlated (r2 ≥ 0.8) with rs10896438, rs10896449, rs12793759, and rs11228565, respectively. HapMap SNPs were in linkage disequilibrium (r2 ≥ 0.8) with 48%, 69%, 14%, and 60% of SNPs marking bins by rs10896438, rs10896449, rs12793759, and rs11228565, respectively. CONCLUSIONS Our next generation resequence analysis compliments publicly available datasets of European descent (HapMap, build 28 and 1000 Genome, Pilot 1, October 2010), underscoring the value of targeted resequence analysis prior to initiating functional studies based on public databases alone. Increasing the number of common variants enables investigators to better prioritize variants for functional studies designed to uncover the biological basis of the direct association(s) in the region. Prostate 72:476–486, 2012. © 2011 Wiley Periodicals, Inc.

The relationship between chromosomal instability and integration of the human papillomavirus (HPV) DNA into the host genome is not well understood. We found that chromosomal instability increases with HPV integration and diagnosis. However, chromosomal instability could occur without HPV integration.

The shelterin complex protects chromosomal ends by regulating how the telomerase complex interacts with telomeres. Following the recent finding in familial melanoma of inactivating germline mutations in POT1, encoding a member of the shelterin complex, we searched for mutations in the other five components of the shelterin complex in melanoma families.Next-generation sequencing techniques were used to screen 510 melanoma families (with unknown genetic etiology) and control cohorts for mutations in shelterin complex encoding genes: ACD, TERF2IP, TERF1, TERF2, and TINF 2. Maximum likelihood and LOD [logarithm (base 10) of odds] analyses were used. Mutation clustering was assessed with χ(2) and Fisher's exact tests. P values under .05 were considered statistically significant (one-tailed with Yates' correction).Six families had mutations in ACD and four families carried TERF2IP variants, which included nonsense mutations in both genes (p.Q320X and p.R364X, respectively) and point mutations that cosegregated with melanoma. Of five distinct mutations in ACD, four clustered in the POT1 binding domain, including p.Q320X. This clustering of novel mutations in the POT1 binding domain of ACD was statistically higher (P = .005) in melanoma probands compared with population control individuals (n = 6785), as were all novel and rare variants in both ACD (P = .040) and TERF2IP (P = .022). Families carrying ACD and TERF2IP mutations were also enriched with other cancer types, suggesting that these variants also predispose to a broader spectrum of cancers than just melanoma. Novel mutations were also observed in TERF1, TERF2, and TINF2, but these were not convincingly associated with melanoma.Our findings add to the growing support for telomere dysregulation as a key process associated with melanoma susceptibility.

The completion of the 5,373,180-bp genome sequence of the marine psychrophilic bacterium Colwellia psychrerythraea 34H, a model for the study of life in permanently cold environments, reveals capabilities important to carbon and nutrient cycling, bioremediation, production of secondary metabolites, and cold-adapted enzymes. From a genomic perspective, cold adaptation is suggested in several broad categories involving changes to the cell membrane fluidity, uptake and synthesis of compounds conferring cryotolerance, and strategies to overcome temperature-dependent barriers to carbon uptake. Modeling of three-dimensional protein homology from bacteria representing a range of optimal growth temperatures suggests changes to proteome composition that may enhance enzyme effectiveness at low temperatures. Comparative genome analyses suggest that the psychrophilic lifestyle is most likely conferred not by a unique set of genes but by a collection of synergistic changes in overall genome content and amino acid composition.

Background: Accurate detection of somatic mutations is critical for informing targeted therapy options. Prevalent non-cancer cell admixture complicates this detection in breast cancer. Conventional mutation detection relies on DNA sequencing; however in prior work, we demonstrated that combining RNA and DNA sequencing increases mutation signal strength, or mutant allele fraction (MAF). The ratio of RNA MAF versus DNA MAF (RNA:DNA MAF) was greatest in low purity breast tumors. We hypothesized that this elevation is biologically driven and would be conserved in a second, distinct tissue specimen of the same tumors. Here, we compare mutation characteristics between two tissue blocks in a cohort of breast tumors (n = 8) to evaluate possible preservation of RNA versus DNA mutation signal throughout the tumor. Methods: We selected four high purity and four low purity breast tumors (“Block1”) from The Cancer Genome Atlas (TCGA) cohort and associated ABSOLUTE purity analysis. For these tumors, we acquired a second tissue block (“Block2”) not analyzed by TCGA, cut sections, analyzed sections by H&E stains, and extracted nucleic acids. Whole genome DNA sequencing and mRNA sequencing was performed for Block2 specimens using Illumina X and NextSeq 500 sequencers, respectively. Somatic mutations in Block2 were detected using UNCeqR and compared to published UNCeqR somatic mutations from TCGA. We then evaluate MAF characteristics in the entire TCGA breast tumor cohort (n = 695). Results: Tumor purity estimates, determined by histology and by sequencing, were reduced in Block2 of the low purity tumor set versus the high purity tumor set, consistent with Block1 analysis. Molecular properties of genome-wide gene expression and somatic DNA copy number were highly similar between block-mated specimens (p < 0.01). We then identified expressed mutations present in Block1 and Block2 of the same tumor and compared the MAFs on these common mutations. DNA MAF and RNA MAF were each significantly correlated between Block1 and Block2 (p < 1e-12 in both cases). The average RNA:DNA MAF was 2.5 for the cohort, indicating that RNA mutation signal is greater than DNA in general. In Block2 specimens, the RNA:DNA MAFs were significantly greater in the low purity tumor set than the high purity tumor set (mean 2.7 versus 2.1, p < 6e-5), reflecting the same trend observed in Block1 specimens. Analyzing the entire TCGA cohort, RNA:DNA MAF was positively correlated with proliferation pathway gene expression (p < 3e-16 ) and was greatest in the Basal subtype versus other subtypes (p < 2e-9). Conclusion: Mutant allele fraction both of DNA and of RNA was conserved across breast tumor subsections. Low purity and basal subtype breast tumors had elevated RNA:DNA MAF supporting a relationship to underlying biology and identifying classes of tumors with pronounced benefit for DNA and RNA integrated mutation analysis. Citation Format: Jerez Te, Coralie Viollet, Xijun Zhang, Jatinder Singh, Jeffrey A. Hooke, Harvey B. Pollard, Hai Hu, Craig D. Shriver, Clifton L. Dalgard, Matthew D. Wilkerson. Reproducible elevation of RNA versus DNA mutation signal in low purity breast tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5403. doi:10.1158/1538-7445.AM2017-5403

Neurofibromatosis type 1 (NF1) is an autosomal dominant, monogenic disorder of dysregulated neurocutaneous tissue growth. Pleiotropy, variable expressivity and few NF1 genotype-phenotype correlates limit clinical prognostication in NF1. Phenotype complexity in NF1 is hypothesized to derive in part from genetic modifiers unlinked to the NF1 locus. In this study, we hypothesized that normal variation in germline gene expression confers risk for certain phenotypes in NF1. In a set of 79 individuals with NF1, we examined the association between gene expression in lymphoblastoid cell lines with NF1-associated phenotypes and sequenced select genes with significant phenotype/expression correlations. In a discovery cohort of 89 self-reported European-Americans with NF1 we examined the association between germline sequence variants of these genes with café-au-lait macule (CALM) count, a tractable, tumor-like phenotype in NF1. Two correlated, common SNPs (rs4660761 and rs7161) between DPH2 and ATP6V0B were significantly associated with the CALM count. Analysis with tiled regression also identified SNP rs4660761 as significantly associated with CALM count. SNP rs1800934 and 12 rare variants in the mismatch repair gene MSH6 were also associated with CALM count. Both SNPs rs7161 and rs4660761 (DPH2 and ATP6V0B) were highly significant in a mega-analysis in a combined cohort of 180 self-reported European-Americans; SNP rs1800934 (MSH6) was near-significant in a meta-analysis assuming dominant effect of the minor allele. SNP rs4660761 is predicted to regulate ATP6V0B, a gene associated with melanosome biology. Individuals with homozygous mutations in MSH6 can develop an NF1-like phenotype, including multiple CALMs. Through a multi-platform approach, we identified variants that influence NF1 CALM count., Author Summary Neurofibromatosis type 1 (NF1) is a relatively common genetic disease that increases the chance to develop a variety of benign and malignant tumors. People with NF1 also typically feature a large number of birthmarks called café-au-lait macules. It is difficult to predict severity or specific problems in NF1. We sought to identify genes (other than NF1, the gene that causes the disease) that influence severity in NF1. We determined the number of café-au-lait macules in two groups of people with NF1. We measured the gene expression of about 10,000 genes in the cultured white blood cells from one group of people. We then sequenced a group of genes whose expression level was increased in people with higher numbers of café-au-lait macules. In the first group, we found common variants in genes MSH6 and near DPH2 and ATP6V0B that were significantly associated with the number of café-au-lait macules. Some of these variants were close to significant in the second group of people. The two variants near DPH2 and ATP6V0B were very significant when analysed in both groups combined. Our work is among the first to identify genetic variants that influence the severity of NF1.

We assessed the performance of the new Life Technologies Proton sequencer by comparing whole-exome sequence data in a Centre d’Etude du Polymorphisme Humain trio (family 1463) to the Illumina HiSeq instrument. To simulate a typical user’s results, we utilized the standard capture, alignment and variant calling methods specific to each platform. We restricted data analysis to include the capture region common to both methods. The Proton produced high quality data at a comparable average depth and read length, and the Ion Reporter variant caller identified 96 % of single nucleotide polymorphisms (SNPs) detected by the HiSeq and GATK pipeline. However, only 40 % of small insertion and deletion variants (indels) were identified by both methods. Usage of the trio structure and segregation of platform-specific alleles supported this result. Further comparison of the trio data with Complete Genomics sequence data and Illumina SNP microarray genotypes documented high concordance and accurate SNP genotyping of both Proton and Illumina platforms. However, our study underscored the problem of accurate detection of indels for both the Proton and HiSeq platforms.

This study examined the effect of P85 (a pluronic block copolymer) and microbubble (MB) ultrasound contrast agents under ultrasound irradiation on gene transfection and expression. The pEGFP plasmids that can encode enhanced green fluorescent protein (pEGFP) served as a report gene and were mixed with different concentrations of MB/0.05% (w/v) P85. Then the plasmids were transfected into human hepatoma G2 (HepG2) cells. The HepG2 cells treated with MB/P85 or without treatment were exposed to ultrasound (US parameters: 1 MHz, 1.0 W/cm(2), 20 s, 20% duty cycle). Twenty-four hours later, the transfection efficiency was assessed by fluorescence microscopy and fluorescence activated cell sorting (FACS) analysis. The cell viability was evaluated by Trypan blue exclusion test. The results showed that the gene transfection efficiency in HepG2 cells under ultrasound irradiation was significantly higher than that without ultrasound irradiation. HepG2 cells in the MB or P85 group in the absence of ultrasound expressed less amount of green fluorescent protein. The expression efficiency reached (22.14 ± 3.06)% and the survival rate was as high as (55.73 ± 3.32)% in the 30% MB plus P85 group. It was concluded that MB and P85 in the presence of ultrasound can enhance gene transfection and expression.

Wistar Furth (WF) intestinal allografts were transplanted into Sprague-Dawley (SD) rats. Recipients were randomly allocated into the following groups: (1) no treatment; (2) cyclosporine (CsA) 6 mg/kg/day, daily, subcutaneously (s.c.; full-dose therapy); (3) CsA 3 mg/kg/day, daily, s.c. (half-dose therapy); and (4) CsA 3 mg/kg/day, daily, s.c. + Tripterygium Wilfordii Hook. WF (TW) 3 mg/kg/day, daily, s.c. WF rats with intestinal autografts were used as controls. CD45RO intragraft expression and its index (CD45RO/CD45), measured by reverse transcription polymerase chain reaction (RT-PCR), were significantly elevated in untreated and half-dose CsA-treated allografts as early as postoperative day (POD) 4, when rejection of intestinal allografts was not detected by routine pathology. Intestinal permeability measured by Tc-DTPA radioassay was significantly elevated in untreated allografts on POD 6. Histology showed that there was severe rejection in untreated intestinal allografts and mild rejection in allografts treated with a half dose of CsA on POD 6. There was a normal CD45RO expression, permeability, and histology in intestinal allografts treated with either a full dose of CsA or a half dose of CsA + TW. These data indicate that CD45RO intragraft gene expression is an early marker to detect intestinal allograft rejection in rats. © 1999 Wiley-Liss, Inc. MICROSURGERY 19:348–350 1999