Heterozygous SNPs Research Articles

Allele-specific expression in the human heart and its application to postoperative atrial fibrillation and myocardial ischemia.

BackgroundAllele-specific expression (ASE) is differential expression of each of the two chromosomal alleles of an autosomal gene. We assessed ASE patterns in the human left atrium (LA, n = 62) and paired samples from the left ventricle (LV, n = 76) before and after ischemia, and tested the utility of differential ASE to identify genes associated with postoperative atrial fibrillation (poAF) and myocardial ischemia.MethodsFollowing genotyping from whole blood and whole-genome sequencing of LA and LV samples, we called ASE using sequences overlapping heterozygous SNPs using rigorous quality control to minimize false ASE calling. ASE patterns were compared between cardiac chambers and with a validation cohort from cadaveric tissue. ASE patterns in the LA were compared between patients who had poAF and those who did not. Changes in ASE in the LV were compared between paired baseline and post-ischemia samples.ResultsASE was found for 3404 (5.1%) and 8642 (4.0%) of SNPs analyzed in the LA and LV, respectively. Out of 6157 SNPs with ASE analyzed in both chambers, 2078 had evidence of ASE in both LA and LV (p < 0.0001). The SNP with the greatest ASE difference in the LA of patients with and without postoperative atrial fibrillation was within the gelsolin (GSN) gene, previously associated with atrial fibrillation in mice. The genes with differential ASE in poAF were enriched for myocardial structure genes, indicating the importance of atrial remodeling in the pathophysiology of AF. The greatest change in ASE between paired post-ischemic and baseline samples of the LV was in the zinc finger and homeodomain protein 2 (ZHX2) gene, a modulator of plasma lipids. Genes with differential ASE in ischemia were enriched in the ubiquitin ligase complex pathway associated with the ischemia-reperfusion response.ConclusionsOur results establish a pattern of ASE in the human heart, with a high degree of shared ASE between cardiac chambers as well as chamber-specific ASE. Furthermore, ASE analysis can be used to identify novel genes associated with (poAF) and myocardial ischemia.Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-016-0381-1) contains supplementary material, which is available to authorized users.

Second generation noninvasive fetal genome analysis reveals de novo mutations, single-base parental inheritance, and preferred DNA ends

Plasma DNA obtained from a pregnant woman was sequenced to a depth of 270× haploid genome coverage. Comparing the maternal plasma DNA sequencing data with the parental genomic DNA data and using a series of bioinformatics filters, fetal de novo mutations were detected at a sensitivity of 85% and a positive predictive value of 74%. These results represent a 169-fold improvement in the positive predictive value over previous attempts. Improvements in the interpretation of the sequence information of every base position in the genome allowed us to interrogate the maternal inheritance of the fetus for 618,271 of 656,676 (94.2%) heterozygous SNPs within the maternal genome. The fetal genotype at each of these sites was deduced individually, unlike previously, where the inheritance was determined for a collection of sites within a haplotype. These results represent a 90-fold enhancement in the resolution in determining the fetus's maternal inheritance. Selected genomic locations were more likely to be found at the ends of plasma DNA molecules. We found that a subset of such preferred ends exhibited selectivity for fetal- or maternal-derived DNA in maternal plasma. The ratio of the number of maternal plasma DNA molecules with fetal preferred ends to those with maternal preferred ends showed a correlation with the fetal DNA fraction. Finally, this second generation approach for noninvasive fetal whole-genome analysis was validated in a pregnancy diagnosed with cardiofaciocutaneous syndrome with maternal plasma DNA sequenced to 195× coverage. The causative de novo BRAF mutation was successfully detected through the maternal plasma DNA analysis.

Somatic Variants in the Human Lens Epithelium: A Preliminary Assessment.

PurposeWe hypothesize that somatic mutations accumulate in cells of the human lens and may contribute to the development of cortical or posterior sub-capsular cataracts. Here, we used a Next-generation sequencing (NGS) strategy to screen for low-allelic frequency variants in DNA extracted from human lens epithelial samples.MethodsNext-Generation sequencing of 151 cancer-related genes (WUCaMP2 panel) was performed on DNA extracted from post-mortem or surgical specimens obtained from 24 individuals. Usually, pairwise comparisons were made between two or more ocular samples from the same individual, allowing putative somatic variants detected in lens samples to be differentiated from germline variants.ResultsUse of a targeted hybridization approach enabled high sequence coverage (>1000-fold) of the WUCaMP2 genes. In addition to high-frequency variants (corresponding to homozygous or heterozygous SNPs and Indels), somatic variants with allelic frequencies of 1-4% were detected in the lens epithelial samples. The presence of one such variant, a T > C point substitution at position 32907082 in BRCA2, was verified subsequently using droplet digital PCR.ConclusionsLow-allelic fraction variants are present in the human lens epithelium, at frequencies consistent with the presence of millimeter-sized clones.

Abstract 5271: Optimization and detection of focal somatic copy number variants in whole genome, whole exome and panel sequencing for tumor/normal matched pairs and tumor only analysis

Abstract Often in the clinical setting, tumor samples may be derived from historical tissue or from fresh frozen tissue and it may not be feasible to obtain normal tissue from the individual. Comparing unmatched tissues from different individuals possess a unique challenge of addressing common copy number variations not faced when comparing matched samples. In addition to the challenges of comparing unmatched samples, DNA extracted from samples preserved as formalin-fixed, paraffin-embedded (FFPE) tissue is often highly degraded and when compared to more intact DNA extracted from different source (e.g. blood) can increase the inherent noise in the system. Here we analytically characterize the accuracy and performance of algorithmic approaches for separating germline inherited copy number variation from somatic copy number changes, focusing on both filtering approaches and use of pooled reference samples sequenced under similar conditions. Example approaches include filtering known common copy number variation within 1000 Genomes Phase 3 and Database of Genomic Variants (DGV) Gold Standard. Additionally, we utilized a tumor/reference pool based analysis where a reference pool was constructed by equimolar pooling of multiple individuals. Determination of fold changes between tumor and reference was calculated by determining physical coverage of read pair fragments in 100 bases increments. Next, to address differences in sequencing performance between tumor and reference, the read depth data for each sample is collapsed/averaged into a lower resolution according to user-selected parameters (e.g. distance between points and read depth). Normalized log2 fold-changes between tumor and reference samples are then calculated and an adjustable smoothing window is applied. In addition, we utilize tumor allele frequencies of known heterozygous germline SNPs identified within the normal to both evaluate potential false positives and correct biases. Lastly, a segmentation algorithm is applied to summarize the individual log2 fold-changes into intervals with a constant copy number state. We will present the advantages and limitations of these approaches both when a germline normal is available and when tumor only analysis is necessary. Citation Format: Jessica Aldrich, Jonathan J. Keats, Austin Christofferson, Winnie S. Liang, John D. Carpten, Lisa Baumbach-Reardon, David W. Craig. Optimization and detection of focal somatic copy number variants in whole genome, whole exome and panel sequencing for tumor/normal matched pairs and tumor only analysis. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 5271.

Whole-genome resequencing of honeybee drones to detect genomic selection in a population managed for royal jelly.

Four main evolutionary lineages of A. mellifera have been described including eastern Europe (C) and western and northern Europe (M). Many apiculturists prefer bees from the C lineage due to their docility and high productivity. In France, the routine importation of bees from the C lineage has resulted in the widespread admixture of bees from the M lineage. The haplodiploid nature of the honeybee Apis mellifera, and its small genome size, permits affordable and extensive genomics studies. As a pilot study of a larger project to characterise French honeybee populations, we sequenced 60 drones sampled from two commercial populations managed for the production of honey and royal jelly. Results indicate a C lineage origin, whilst mitochondrial analysis suggests two drones originated from the O lineage. Analysis of heterozygous SNPs identified potential copy number variants near to genes encoding odorant binding proteins and several cytochrome P450 genes. Signatures of selection were detected using the hapFLK haplotype-based method, revealing several regions under putative selection for royal jelly production. The framework developed during this study will be applied to a broader sampling regime, allowing the genetic diversity of French honeybees to be characterised in detail.

Elevated variant density around SV breakpoints in germline lineage lends support to error-prone replication hypothesis.

Copy number variants (CNVs) are a class of structural variants that may involve complex genomic rearrangements (CGRs) and are hypothesized to have additional mutations around their breakpoints. Understanding the mechanisms underlying CNV formation is fundamental for understanding the repair and mutation mechanisms in cells, thereby shedding light on evolution, genomic disorders, cancer, and complex human traits. In this study, we used data from the 1000 Genomes Project to analyze hundreds of loci harboring heterozygous germline deletions in the subjects NA12878 and NA19240. By utilizing synthetic long-read data (longer than 2 kbp) in combination with high coverage short-read data and, in parallel, by comparing with parental genomes, we interrogated the phasing of these deletions with the flanking tens of thousands of heterozygous SNPs and indels. We found that the density of SNPs/indels flanking the breakpoints of deletions (in-phase variants) is approximately twice as high as the corresponding density for the variants on the haplotype without deletion (out-of-phase variants). This fold change was even larger for the subset of deletions with signatures of replication-based mechanism of formation. The allele frequency (AF) spectrum for deletions is enriched for rare events; and the AF spectrum for in-phase SNPs is shifted toward this deletion spectrum, thus offering evidence consistent with the concomitance of the in-phase SNPs/indels with the deletion events. These findings therefore lend support to the hypothesis that the mutational mechanisms underlying CNV formation are error prone. Our results could also be relevant for resolving mutation-rate discrepancies in human and to explain kataegis.

IntSIM: An Integrated Simulator of Next-Generation Sequencing Data.

Next-generation sequencing data has been widely used for DNA variant discovery and tumor study through computational tools. Effective simulation of such data with many realistic features is very necessary for testing existing tools and guiding the development of new tools. We present an integrated simulation system, IntSIM, to simulate common DNA variants and to generate sequencing reads for mixture genomes. IntSIM has three novel features in comparison with other simulation programs: 1) it is able to simulate both germline and somatic variants in the same sequence, 2) it deals with tumor purity so as to generate reads corresponding to heterogeneous genomes and also produce tumor-normal matched samples, and 3) it simulates correlations among SNPs, among CNVs/CNAs based on HMM models trained from real sequencing genomes, and can simulates broad and focal CNV/CNA events. The simulation data of IntSIM can reflect characteristics observed from real data and are consistent with input parameters. The IntSIM software package is freely available at http://intsim.sourceforge.net/. Based on a great number of experiments, IntSIM performs better than other program for some scenarios, such as simulation of heterozygous SNPs, CNVs/CNAs, and can achieve some functions that other programs cannot achieve. Simulation with IntSIM can be expected to evaluate performance of methods in detecting various types of variants, analyzing tumor samples, and especially providing a realistic assessment of effect of tumor purity on identification of somatic mutations.

Abstract 45: Detection of focal somatic copy number variants in whole genome, whole exome, and targeted next-generation sequencing data of tumor/normal pairs

Abstract We describe the development and validation of an analysis pipeline for detection of focal somatic copy number aberrations from next-generation sequencing (NGS) data. In the clinical setting, tumor samples may be derived from historical tissue or from recently collected fresh frozen biopsies. DNA extracted from samples preserved as formalin-fixed, paraffin-embedded (FFPE) tissue are often highly degraded. We describe challenges inherent to sequencing different tissue sources and algorithmic approaches for optimizing detection of focal alterations in these samples with a high degree of accuracy. The first area of optimization takes into account clonal or fragment coverage based on user-defined intervals by accounting for coverage between paired-end reads. Normalized log2 fold-changes between tumor and normal samples are then calculated and a smoothing window is applied. In addition, we utilize tumor allele frequencies of known heterozygous germline SNPs identified within the normal to both evaluate potential false positives and correct biases. Lastly, a segmentation algorithm (circular binary segmentation as implemented in DNAcopy package1) is applied to summarize the individual log2 fold-changes into intervals with a constant copy number state. The results from the segmentation algorithm are then annotated and reformatted into the VCF format. We have analytically optimized and validated our algorithm across a series of samples with known alterations for implementation on both clinical and research samples. 1. Seshan VE and Olshen A. DNAcopy: DNA copy number data analysis. R package version 1.40.0. Citation Format: Jessica Aldrich, Jonathan J. Keats, Winnie S. Liang, John D. Carpten, David W. Craig. Detection of focal somatic copy number variants in whole genome, whole exome, and targeted next-generation sequencing data of tumor/normal pairs. [abstract]. In: Proceedings of the AACR Precision Medicine Series: Integrating Clinical Genomics and Cancer Therapy; Jun 13-16, 2015; Salt Lake City, UT. Philadelphia (PA): AACR; Clin Cancer Res 2016;22(1_Suppl):Abstract nr 45.

Implementing Non-Invasive Prenatal Diagnosis (NIPD) in a National Health Service Laboratory; From Dominant to Recessive Disorders.

Our UK National Health Service regional genetics laboratory offers NIPD for autosomal dominant and de novo conditions (achondroplasia, thanataphoric dysplasia, Apert syndrome), paternal mutation exclusion for cystic fibrosis and a range of bespoke tests. NIPD avoids the risks associated with invasive testing, making prenatal diagnosis more accessible to families at high genetic risk. However, the challenge remains in offering definitive diagnosis for autosomal recessive diseases, which is complicated by the predominance of the maternal mutant allele in the cell-free DNA sample and thus requires a variety of different approaches. Validation and diagnostic implementation for NIPD of congenital adrenal hyperplasia (CAH) is further complicated by presence of a pseudogene that requires a different approach. We have used an assay targeting approximately 6700 heterozygous SNPs around the CAH gene (CYP21A2) to construct the high-risk parental haplotypes and tested this approach in five cases, showing that inheritance of the parental alleles can be correctly identified using NIPD. We are evaluating various measures of the fetal fraction to help determine inheritance of parental mutations. We are currently exploring the utility of an NIPD multi-disorder panel for autosomal recessive disease, to make testing more widely applicable to families with a variety of serious genetic conditions.

Abstract 4813: High-depth sequencing of 800 cancer-associated genes in 48 matched tumor-normal tissues to identify somatic alterations including low-frequency variants in colorectal cancer

Abstract Purpose: Colorectal cancer is the 3rd most common cancer globally accounting for 10% of the global cancer burden. Catalogues of somatic mutations in cancer (e.g. The Cancer Genome Atlas (TCGA)) are based upon exome sequencing, typically performed at depths of ∼ 80x, which will miss sub-clonal low-frequency somatic variants. We performed high depth sequencing of 800 cancer-associated genes across 48 tumor/normal pairs of CRC to catalogue somatic variants including those present at low variant allele frequencies. Method: We used the Agilent SureSelect Target Enrichment system to design a customized enrichment panel comprising all coding exons of ∼800 cancer-associated genes. From genomic DNA extracted from frozen tissue specimens (matched primary tumor and adjacent normal), we performed library construction, target enrichment and illumina next generation sequencing. We used 3 variant-calling algorithms - a Genome Analyzer Toolkit (GATK) based pipeline, LoFreq and MuTect to identify somatic variants including sub-clonal variants present at low variant allele frequencies. Results: Sequencing was performed to a median of 302-fold and 310-fold coverage across target regions in the tumor tissues and normal tissues respectively. The use of 3 bioinformatics variant-calling algorithm allowed sensitive and comprehensive identification of variants down to a variant allele frequency of 2.5%; selected variants were validated by an orthogonal sequencing platform (Ion Torrent). The validation rate was &amp;gt; 97%. Each patient had a variant in a median of 19 genes (range 8-100) with a different complement of clonal and sub-clonal variants. Recurrent alterations (e.g. APC, KRAS, TP53) tended to be clonal. Broad regions of allelic imbalance (based on heterozygous SNPs identified in each patient's normal tissue) could also be mapped in each patient. Conclusion: High Depth sequencing of 800 cancer-associated genes in 48 matched tumor normal tissues permits identification of somatic alterations including low-frequency variants and sub-clonal alterations in patients with colorectal cancer. Citation Format: Clarinda Chua, Dan Liang Ho, Yuka Suzuki, Simeen Malik, John McPherson, Anna Gan, Dennis Koh, Choon Leong Tang, Sarah Boon Hsi Ng, Patrick Tan, Steve Rozen, Iain Beehuat Tan. High-depth sequencing of 800 cancer-associated genes in 48 matched tumor-normal tissues to identify somatic alterations including low-frequency variants in colorectal cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4813. doi:10.1158/1538-7445.AM2015-4813

Abstract 1552: Image-based microchip sorting of pure, immuno-phenotypically defined subpopulations of tumor cells from tiny formalin-fixed paraffin embedded (FFPE) samples reveals their distinct genetic features

Abstract Background: We provide a solution of pressing needs in preparation of FFPE samples for genomic analysis: small sample size, unwanted admixture of normal cells, analysis of tumor rare-cell subpopulations present at low percentages in the tumor fraction. Methods: We disaggregated into cell suspensions archival FFPE samples from 12 ovarian, pancreatic and lung cancer patients, staining for Vimentin, Keratin and DNA. We sorted by DEPArray™ precise numbers (mean = 107, median 58, range = 5-600) of pure homogenous cells from the major population of tumor cells, the contaminant diploid stromal cells, and other minority tumor cell types indicative of epithelial-to-mesenchymal transition (EMT). Using IonTorrent AmpliSeq CHPv2, we generated sequencing libraries, after lysis of the pure cells recovered by DEPArray™ (n = 54), or unsorted samples (either QIAmp DNA columns or disaggregated cells). Libraries were sequenced with IonTorrent PGM (mean depth&amp;gt;2,000x), and analyzed using IonTorrent software. Results: On several loci, we detected somatic mutations with 100% variant frequency, only observable as heterozygous in the unsorted samples and as wild-type in stromal cells of same patient, confirming 100% purity of sorted cells. Moreover, in the EMT-phenotype subpopulations we identified clear somatic mutations, different from tumor cells majority and undetectable in unsorted samples. Frequently, for loci harboring germ-line heterozygous SNPs with variant frequency around 50% for pure stromal cells, we readily detected loss-of-heterozygosis in tumor cells subpopulations as binary (0%/100%) variants. Quantitative traits such as copy number gains and losses were also reproducibly identified in tumor cell replicates as deviations from the 50% variant frequency of germline SNPs of pure stromal cells. Furthermore, we observed an excellent coverage uniformity (mean = 96%) for recoveries (n = 27) in the range of 81-600 cells, even higher than the uniformity obtained with (n = 2) QIAmp-purified DNA (92%). Mean uniformity gradually decreased to 89% for cell recoveries (n = 13) in the range 21-80, and further decreased to 70% for lower cell numbers (n = 14). Highlights: Sorting tumor rare-cell subpopulations reveals their genetic characteristics, undetectable in unsorted samples. Analyzing homogenous cell subpopulations boosts signal-to-noise ratio working around inherent sensitivity/specifitiy trade-offs of rare-variant calls. The proposed workflow further enables reliable detection of quantitative traits such as CNVs. Sorting pure stromal cells yields internal controls for archival samples. Citation Format: Chiara Bolognesi, Anna Doffini, Genny Buson, Rossana Lanzellotto, Giulio Signorini, Valeria Sero, Alex Calanca, Francesca Fontana, Rita Romano, Stefano Gianni, Giulia Bregola, Gianni Medoro, Raimo Tanzi, Giuseppe Giorgini, Hans Morreau, Massimo Barberis, Willem E. Corver, Nicolo Manaresi. Image-based microchip sorting of pure, immuno-phenotypically defined subpopulations of tumor cells from tiny formalin-fixed paraffin embedded (FFPE) samples reveals their distinct genetic features. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 1552. doi:10.1158/1538-7445.AM2015-1552

Abstract P4-02-03: Detection of single nucleotide variations and copy number variations in breast cancer tissue and ctDNA samples using single-nucleotide polymorphism-targeted massively multiplexed PCR

Abstract Genomic instability, the hallmark of cancer, presents with a variety of mutation types, most commonly single nucleotide variations (SNVs) and copy number variations (CNVs), which traditionally have required different methods for identification. It has proven challenging to simultaneously achieve sufficient breadth to detect CNVs and depth to detect SNVs on samples of limited input amount. The objective of this study was to validate a new methodology for detection of SNVs and CNVs in a single assay. We used a massively multiplex PCR/NGS approach combining an SNV panel covering 585 point mutation hotspots in breast cancer (Cosmic) and a CNV panel targeting 28,000 SNPs designed to detect copy number at chromosomes 1, 2, 13, 18, 21, and X, and focal regions 4p16, 5p15, 7q11, 15q, 17p, 22q11, and 22q13. We applied these panels to breast cancer cell lines and fresh frozen (FF) breast tumor samples; the presence of CNVs in circulating cell-free tumor DNA (ctDNA) in the plasma of breast cancer patients was also investigated. The CNV assay methodology was validated using genomic DNA isolated from 96 human samples with known karyotype; sensitivity to single region deletions or duplications was 100% (71/71) and specificity was 100% for normal regions in the same samples. Single-molecule sensitivity for the detection of CNVs was established by analyzing isolated single cells. Performance of the mutation assay was demonstrated with the analysis of 5 matched tumor and normal cell lines, with 24 out of 27 SNVs known to be present in these cell lines detected. The 3 undetected SNVs were determined to be a result of assay design failure. Also, multiple somatic CNVs (median: 13) were detected in all 5 tumor cell lines. Analysis of the normal cell lines found no cancer related SNVs or CNVs. In 32 FF tumor samples, 78.1% (25/32) had SNVs detected; of samples with SNVs, 88% (22/25) had SNVs in TP53 or PIK3CA. Of the same 32 FF breast tumor samples, 96.9% (31/32) showed full or partial CNVs in at least 1 and up to 15 regions; of the 31 samples with detected CNVs, 93.5% had a CNV of either 1q or 17p, two of the three most prevalent breast cancer CNVs (the 16q region was not represented in this panel). Overall, a combination of SNV and CNV testing allowed identification of genetic changes in 100% of the breast tumor samples, a significant improvement in diagnostic yield than using SNV detection alone. Of the 12 breast cancer patients with matched tumor tissue and plasma samples, 83.3% (10/12) had CNVs detected in tissue. The CNVs present in each primary tumor sample were identified in corresponding plasma ctDNA samples (1 stage IIa, 7 stage IIb, and 2 stage III). The ctDNA fractions in these samples ranged from 0.58 to 4.33%; detection required as few as 86 heterozygous SNPs per CNV. Analysis of ctDNA for cancer-associated mutations may allow earlier, safer and more accurate profiling and monitoring of breast cancer. Thus, this targeted PCR approach offers the promise of an assay able to detect both cancer-associated SNVs and CNVs in the same sample with good sensitivity and specificity, and improved detection rates compared to assays that only detect SNVs. Citation Format: Robert J Pelham, Bernhard G Zimmermann, Eser Kirkizlar, Ryan K Swenerton, Bin Hoang, Onur Sakarya, Joshua E Babiarz, Nicholas Wayham, Tudor Constantin, Styrmir Sigurjonsson, Matthew Rabinowitz, Matthew Hill. Detection of single nucleotide variations and copy number variations in breast cancer tissue and ctDNA samples using single-nucleotide polymorphism-targeted massively multiplexed PCR [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P4-02-03.

A hidden Markov approach for ascertaining cSNP genotypes from RNA sequence data in the presence of allelic imbalance by exploiting linkage disequilibrium.

BackgroundAllelic specific expression (ASE) increases our understanding of the genetic control of gene expression and its links to phenotypic variation. ASE testing is implemented through binomial or beta-binomial tests of sequence read counts of alternative alleles at a cSNP of interest in heterozygous individuals. This requires prior ascertainment of the cSNP genotypes for all individuals. To meet the needs, we propose hidden Markov methods to call SNPs from next generation RNA sequence data when ASE possibly exists.ResultsWe propose two hidden Markov models (HMMs), HMM-ASE and HMM-NASE that consider or do not consider ASE, respectively, in order to improve genotyping accuracy. Both HMMs have the advantages of calling the genotypes of several SNPs simultaneously and allow mapping error which, respectively, utilize the dependence among SNPs and correct the bias due to mapping error. In addition, HMM-ASE exploits ASE information to further improve genotype accuracy when the ASE is likely to be present.Simulation results indicate that the HMMs proposed demonstrate a very good prediction accuracy in terms of controlling both the false discovery rate (FDR) and the false negative rate (FNR). When ASE is present, the HMM-ASE had a lower FNR than HMM-NASE, while both can control the false discovery rate (FDR) at a similar level. By exploiting linkage disequilibrium (LD), a real data application demonstrate that the proposed methods have better sensitivity and similar FDR in calling heterozygous SNPs than the VarScan method. Sensitivity and FDR are similar to that of the BCFtools and Beagle methods. The resulting genotypes show good properties for the estimation of the genetic parameters and ASE ratios.ConclusionsWe introduce HMMs, which are able to exploit LD and account for the ASE and mapping errors, to simultaneously call SNPs from the next generation RNA sequence data. The method introduced can reliably call for cSNP genotypes even in the presence of ASE and under low sequencing coverage. As a byproduct, the proposed method is able to provide predictions of ASE ratios for the heterozygous genotypes, which can then be used for ASE testing.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-015-0479-2) contains supplementary material, which is available to authorized users.

High-depth sequencing of over 750 genes supports linear progression of primary tumors and metastases in most patients with liver-limited metastatic colorectal cancer

BackgroundColorectal cancer with metastases limited to the liver (liver-limited mCRC) is a distinct clinical subset characterized by possible cure with surgery. We performed high-depth sequencing of over 750 cancer-associated genes and copy number profiling in matched primary, metastasis and normal tissues to characterize genomic progression in 18 patients with liver-limited mCRC.ResultsHigh depth Illumina sequencing and use of three different variant callers enable comprehensive and accurate identification of somatic variants down to 2.5% variant allele frequency. We identify a median of 11 somatic single nucleotide variants (SNVs) per tumor. Across patients, a median of 79.3% of somatic SNVs present in the primary are present in the metastasis and 81.7% of all alterations present in the metastasis are present in the primary. Private alterations are found at lower allele frequencies; a different mutational signature characterized shared and private variants, suggesting distinct mutational processes. Using B-allele frequencies of heterozygous germline SNPs and copy number profiling, we find that broad regions of allelic imbalance and focal copy number changes, respectively, are generally shared between the primary tumor and metastasis.ConclusionsOur analyses point to high genomic concordance of primary tumor and metastasis, with a thick common trunk and smaller genomic branches in general support of the linear progression model in most patients with liver-limited mCRC. More extensive studies are warranted to further characterize genomic progression in this important clinical population.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-015-0589-1) contains supplementary material, which is available to authorized users.

Paternal X inactivation does not correlate with X chromosome evolutionary strata in marsupials.

BackgroundX chromosome inactivation is the transcriptional silencing of one X chromosome in the somatic cells of female mammals. In eutherian mammals (e.g. humans) one of the two X chromosomes is randomly chosen for silencing, with about 15% (usually in younger evolutionary strata of the X chromosome) of genes escaping this silencing. In contrast, in the distantly related marsupial mammals the paternally derived X is silenced, although not as completely as the eutherian X. A chromosome wide examination of X inactivation, using RNA-seq, was recently undertaken in grey short-tailed opossum (Monodelphis domestica) brain and extraembryonic tissues. However, no such study has been conduced in Australian marsupials, which diverged from their American cousins ~80 million years ago, leaving a large gap in our understanding of marsupial X inactivation.ResultsWe used RNA-seq data from blood or liver of a family (mother, father and daughter) of tammar wallabies (Macropus eugenii), which in conjunction with available genome sequence from the mother and father, permitted genotyping of 42 expressed heterozygous SNPs on the daughter’s X. These 42 SNPs represented 34 X loci, of which 68% (23 of the 34) were confirmed as inactivated on the paternally derived X in the daughter’s liver; the remaining 11 X loci escaped inactivation. Seven of the wallaby loci sampled were part of the old X evolutionary stratum, of which three escaped inactivation. Three loci were classified as part of the newer X stratum, of which two escaped inactivation. A meta-analysis of previously published opossum X inactivation data revealed that 5 of 52 genes in the old X stratum escaped inactivation.ConclusionsWe demonstrate that chromosome wide inactivation of the paternal X is common to an Australian marsupial representative, but that there is more escape from inactivation than reported for opossum (32% v 14%). We also provide evidence that, unlike the human X chromosome, the location of loci within the oldest evolutionary stratum on the marsupial X does not correlate with their probability of escape from inactivation.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-014-0267-z) contains supplementary material, which is available to authorized users.

Discovering chromatin motifs using FAIRE sequencing and the human diploid genome

BackgroundSpecific chromatin structures are associated with active or inactive gene transcription. The gene regulatory elements are intrinsically dynamic and alternate between inactive and active states through the recruitment of DNA binding proteins, such as chromatin-remodeling proteins.ResultsWe developed a unique genome-wide method to discover DNA motifs associated with chromatin accessibility using formaldehyde-assisted isolation of regulatory elements with high-throughput sequencing (FAIRE-seq). We aligned the FAIRE-seq reads to the GM12878 diploid genome and subsequently identified differential chromatin-state regions (DCSRs) using heterozygous SNPs. The DCSR pairs represent the locations of imbalances of chromatin accessibility between alleles and are ideal to reveal chromatin motifs that may directly modulate chromatin accessibility. In this study, we used DNA 6-10mer sequences to interrogate all DCSRs, and subsequently discovered conserved chromatin motifs with significant changes in the occurrence frequency. To investigate their likely roles in biology, we studied the annotated protein associated with each of the top ten chromatin motifs genome-wide, in the intergenic regions and in genes, respectively. As a result, we found that most of these annotated motifs are associated with chromatin remodeling, reflecting their significance in biology.ConclusionsOur method is the first one using fully phased diploid genome and FAIRE-seq to discover motifs associated with chromatin accessibility. Our results were collected to construct the first chromatin motif database (CMD), providing the potential DNA motifs recognized by chromatin-remodeling proteins and is freely available at http://syslab.nchu.edu.tw/chromatin.

No mutations in the serotonin related TPH1 and HTR1B genes in patients with monogenic sclerosing bone disorders

Since the identification of LRP5 as the causative gene for the osteoporosis pseudoglioma syndrome (OPPG) as well as the high bone mass (HBM) phenotype, LRP5 and the Wnt/β-catenin signaling have been extensively studied for their role in the differentiation and proliferation of osteoblasts, in the apoptosis of osteoblasts and osteocytes and in the response of bone to mechanical loading. However, more recently the direct effect of LRP5 on osteoblasts and bone formation has been questioned. Gene expression studies showed that mice lacking lrp5 have increased expression of tph1, the rate limiting enzyme for the production of serotonin in the gut. Furthermore mice lacking either tph1 or htr1B, the receptor for serotonin on the osteoblasts, were reported to have an increased bone mass due to increased bone formation. This led to the still controversial hypothesis that LRP5 influences bone formation indirectly by regulating the expression of thp1 and as a consequence influencing the production of serotonin in the gut. Based on these data we decided to evaluate the role of TPH1 and HTR1B in the development of craniotubular hyperostoses, a group of monogenic sclerosing bone dysplasias. We screened the coding regions of both genes in 53 patients lacking a mutation in the known causative genes LRP5, LRP4 and SOST. We could not find disease-causing coding variants in neither of the tested genes and therefore, we cannot provide support for an important function of TPH1 and HTR1B in the pathogenesis of sclerosing bone dysplasias in our tested patient cohort.

Whole-genome haplotyping by dilution, amplification, and sequencing

Standard whole-genome genotyping technologies are unable to determine haplotypes. Here we describe a method for rapid and cost-effective long-range haplotyping. Genomic DNA is diluted and distributed into multiple aliquots such that each aliquot receives a fraction of a haploid copy. The DNA template in each aliquot is amplified by multiple displacement amplification, converted into barcoded sequencing libraries using Nextera technology, and sequenced in multiplexed pools. To assess the performance of our method, we combined two male genomic DNA samples at equal ratios, resulting in a sample with diploid X chromosomes with known haplotypes. Pools of the multiplexed sequencing libraries were subjected to targeted pull-down of a 1-Mb contiguous region of the X-chromosome Duchenne muscular dystrophy gene. We were able to phase the Duchenne muscular dystrophy region into two contiguous haplotype blocks with a mean length of 494 kb. The haplotypes showed 99% agreement with the consensus base calls made by sequencing the individual DNAs. We subsequently used the strategy to haplotype two human genomes. Standard genomic sequencing to identify all heterozygous SNPs in the sample was combined with dilution-amplification-based sequencing data to resolve the phase of identified heterozygous SNPs. Using this procedure, we were able to phase >95% of the heterozygous SNPs from the diploid sequence data. The N50 for a Yoruba male DNA was 702 kb whereas the N50 for a European female DNA was 358 kb. Therefore, the strategy described here is suitable for haplotyping of a set of targeted regions as well as of the entire genome.

Do LQTS Gene Single Nucleotide Polymorphisms Alter QTc Intervals at Rest and during Exercise Stress Testing?

The impact of harboring, genetic variants or single nucleotide polymorphisms (LQT-PM) on the repolarization response during exercise and recovery is unknown. To assess the QTc interval adaptation during exercise stress testing (EST) in children with LQT polymorphisms compared to a group of age and gender matched normal controls. One hundred forty-eight patients were age and gender matched into two groups: LQT-PM and control. Each patient underwent a uniform exercise protocol employing a cycle ergometer followed by a 9 minute recovery phase with continuous 12-lead electrocardiogram (ECG) monitoring. Intervals (RR, QT and QTc) at rest (supine), peak exercise and in recovery (1, 3, 5, 7, and 9 minutes) were measured. Forty-three patients were positive for LQT-PM and the control group consisted of 105 patients. A total of 83 SNPs were identified: SCN5A n = 31 (37%), KCNE1 n = 29 (35%), KCNH2 n = 20 (24%), KCNQ1 n = 2 (2%) and KCNE2 n = 1 (1%). The QTc interval measurements of the LQT-PM were longer at rest, peak exercise and all phases of recovery when compared to the control group. Neither group demonstrated abnormal QTc interval adaptation in response to exercise. Patients with homozygous SNPs had longer resting QTc intervals when compared to patients with only heterozygous SNPs (435 ± 23 ms vs. 415 ± 20 ms, respectively, P value <0.006). Individuals with LQT-PM may have longer QTc intervals at rest as well as at peak exercise and all phases of the recovery period compared to normal controls. Additionally, subjects with homozygous SNPs had longer resting QTc intervals when compared to those with only heterozygous SNPs.

Exome Sequencing and Linkage Analysis Implicates Two Candidate Genes On Chromosome 3p in Familial Hodgkin Lymphoma

Abstract 53 BackgroundIt is well recognised that a genetic component exists in classical Hodgkin lymphoma. Higher rates are seen in first-degree relatives of affected individuals, with high concordance observed in monozygotic compared to dizygotic twins. Numerous associations with HLA alleles have been demonstrated in both Epstein-Barr virus-positive and negative forms of the disease and several non-HLA genes have now been implicated. However, these genetic factors are insufficient to account for all the observed inherited risk. MethodsWe report a family from North-East Scotland containing five related individuals with classical Hodgkin lymphoma. DNA was extracted from either whole peripheral blood or saliva. Affected individuals were genotyped using the Affymetrix 500Kb SNP array. Genome-wide linkage analysis was performed using easyLINKAGE plus. Targeted exome capture was performed on four affected individuals using the Agilent SureSelect 50Mb kit. Exome sequencing was performed on the Illumina HiSeq2000 to produce 101bp paired-end reads. Following quality control, raw fastq reads were aligned to the human genome (GRCh37) using bwa. SNP/indel calls were performed using Samtools. Variants identified were filtered by quality to include only those with a read coverage of ≥ 20 and phred score of ≥ 30. Given the high penetrance in this family, we hypothesised that a causative genetic variant would not have been described before so eliminated all SNPs present in dbSNP (version 132) or the 1000 genomes project. Finally, only SNPs affecting protein structure were considered further. SIFT and Polyphen-2 were used to predict impact of genetic variation upon protein structure/function. ResultsKaryotypic analysis was performed on one individual at diagnosis. This was normal, suggesting that a consitutional chromosomal disorder is not responsible for classical Hodgkin lymphoma in this family. The two regions exhibiting strongest linkage with disease were observed on chromosome 3p (LOD scores 1.5 & 1.4). Exome sequencing revealed seven novel, non-synonymous, heterozygous SNPs present in all four analysed individuals. Three were considered biologically plausible (FAM107A:A89S, ALS2CL:L440V, IGSF3:E414G) on the basis of known function and pattern of expression. Two of these three genes (FAM107A and ALS2CL) are on chromosome 3p, in regions similar to those identified in the linkage analysis, and have been implicated as candidate tumor suppressor genes in other malignancies. All three variants were located in evolutionary conserved regions. However, only the mutations in FAM107A and IGSF3 were predicted to disrupt protein function. DiscussionTwo genes on chromosome 3p represent candidate loci for causing familial classical Hodgkin lymphoma. FAM107A protein is downregulated in several malignant cell lines and primary tumor cells. Overexpression can result in suppression of tumor growth. A study of head and neck carcinoma identified frequent missense mutations and/or loss of heterozygosity at ALS2CL, suggesting that this gene can contribute to tumorigenesis. We are evaluating these genes further by assessing gene expression and protein levels in both Hodgkin cell lines and primary Reed-Sternberg cells. Finally, we intend to evaluate the presence of these candidate genetic variants in other, unaffected members of this family and sporadic cases. Disclosures:No relevant conflicts of interest to declare.