Read Depth Analysis Research Articles

High incidence of AZF duplications in clan-structured Iranian populations detected through Y chromosome sequencing read depth analysis

The ampliconic region of the human Y chromosome consists of large duplicated sequences that can undergo non-allelic homologous recombination (NAHR), resulting in structural rearrangements that may cause infertility, especially when they occur in the azoospermia factor b/c (AZFb/c) region. Although AZF duplications have long been neglected due to the technical limitations of STS-based studies that focused mainly on deletions, recent next generation sequencing (NGS) technologies provided evidence for their importance in fertility. In this study, a NGS read depth approach was used to detect AZFb/c rearrangements in 87 Iranians from different ethnic groups. The duplication frequency in Iran proved to be twice as high as in the "1000 Genomes" dataset. Interestingly, most duplications were found in patrilineal ethnic groups, possibly as a consequence of their lower male effective population size which can counteract negative selection. Moreover, we found a large 8.0 Mb duplication, resulting in a fourfold increase in the copy number of AZFc genes, which to our knowledge is the largest duplication ever reported in this region. Overall, our results suggest that it is important to consider not only AZF deletions but also duplications to investigate the causes of male infertility, especially in patrilineal clan-based populations.

Copy Number Variations in Hereditary Spastic Paraplegia-Related Genes: Evaluation of an Iranian Hereditary Spastic Paraplegia Cohort and Literature Review

Introduction: In human genetic disorders, copy number variations (CNVs) are considered a considerable underlying cause. CNVs are generally detected by array-based methods but can also be discovered by read-depth analysis of whole-exome sequencing (WES) data. We performed WES-based CNV identification in a cohort of 35 Iranian families with hereditary spastic paraplegia (HSP) patients. Methods: Thirty-five patients whose routine single-nucleotide variants (SNVs) and insertion/deletion analyses from exome data were unrevealing underwent a pipeline of CNV analysis using the read-depth detection method. Subsequently, a comprehensive search about the existence of CNVs in all 84 known HSP-causing genes was carried out in all reported HSP cases, so far. Results and Discussion: CNV analysis of exome data indicated that 1 patient harbored a heterozygous deletion in exon 17 of the SPAST gene. Multiplex ligation-dependent probe amplification analysis confirmed this deletion in the proband and his affected father. Literature review demonstrated that, to date, pathogenic CNVs have been identified in 30 out of 84 HSP-causing genes (∼36%). However, CNVs in only 17 of these genes were specifically associated with the HSP phenotype. Among them, CNVs were more common in L1CAM, PLP1, SPAST, SPG7, SPG11, and REEP1 genes. The identification of the CNV in 1 of our patients suggests that WES allows the detection of both SNVs and CNVs from a single method without additional costs and execution time. However, because of intrinsic issues of WES in the detection of large rearrangements, it may not yet be exploited to replace the CNV detection methods in standard clinical practice.

Targeted nanopore sequencing enables complete characterisation of structural deletions initially identified using exon-based short-read sequencing strategies.

The widespread adoption of exome sequencing has greatly increased the rate of genetic diagnosis for inherited conditions. However, the detection and validation of large deletions remains challenging. While numerous bioinformatics approaches have been developed to detect deletions from whole - exome sequencing and targeted panels, further work is typically required to define the physical breakpoints or integration sites. Accurate characterisation requires either expensive follow - up whole - genome sequencing or the time - consuming, laborious process of PCR walking, both of which are challenging when dealing with the repeat sequences which frequently intersect deletion breakpoints. The aim of this study was to develop a cost-effective, long-range sequencing method to characterise deletions. Genomic DNA was amplified with primers spanning the deletion using long-range PCR and the products purified. Sequencing was performed on MinION flongle flowcells. The resulting fast5 files were basecalled using Guppy, trimmed using Porechop and aligned using Minimap2. Filtering was performed using NanoFilt. Nanopore sequencing results were verified by Sanger sequencing. Four cases with deletions detected following comparative read-depth analysis of targeted short-read sequencing were analysed. Nanopore sequencing defined breakpoints at the molecular level in all cases including homozygous breakpoints in EYS, CNGA1 and CNGB1 and a heterozygous deletion in PRPF31. All breakpoints were verified by Sanger sequencing. In this study, a quick, accurate and cost - effective method is described to characterise deletions identified from exome, and similar data, using nanopore sequencing.

Chromosome-level de novo genome assembly of Telopea speciosissima (New South Wales waratah) using long-reads, linked-reads and Hi-C.

Telopea speciosissima, the New South Wales waratah, is an Australian endemic woody shrub in the family Proteaceae. Waratahs have great potential as a model clade to better understand processes of speciation, introgression and adaptation, and are significant from a horticultural perspective. Here, we report the first chromosome-level genome for T. speciosissima. Combining Oxford Nanopore long-reads, 10x Genomics Chromium linked-reads and Hi-C data, the assembly spans 823Mb (scaffold N50 of 69.0Mb) with 97.8% of Embryophyta BUSCOs "Complete". We present a new method in Diploidocus (https://github.com/slimsuite/diploidocus) for classifying, curating and QC-filtering scaffolds, which combines read depths, k-mer frequencies and BUSCO predictions. We also present a new tool, DepthSizer (https://github.com/slimsuite/depthsizer), for genome size estimation from the read depth of single-copy orthologues and estimate the genome size to be approximately 900Mb. The largest 11scaffolds contained 94.1% of the assembly, conforming to the expected number of chromosomes (2n = 22). Genome annotation predicted 40,158 protein-coding genes, 351 rRNAs and 728 tRNAs. We investigated CYCLOIDEA (CYC) genes, which have a role in determination of floral symmetry, and confirm the presence of two copies in the genome. Read depth analysis of 180 "Duplicated" BUSCO genes using a new tool, DepthKopy (https://github.com/slimsuite/depthkopy), suggests almost all are real duplications, increasing confidence in the annotation and highlighting a possible need to revise the BUSCO set for this lineage. The chromosome-level T. speciosissima reference genome (Tspe_v1) provides an important new genomic resource of Proteaceae to support the conservation of flora in Australia and further afield.

Sequence Read Depth Analysis of a Monophyletic Cluster of Y Chromosomes Characterized by Structural Rearrangements in the AZFc Region Resulting in DYS448 Deletion and DYF387S1 Duplication.

The azoospermia factor c region (AZFc), located in the long arm of the human Y chromosome, is frequently involved in chromosome rearrangements, mainly due to non-allelic homologous recombination events that occur between the nearly identical sequences (amplicon) that comprises it. These rearrangements may have major phenotypic effects like spermatogenic failure or other pathologies linked to male infertility. Moreover, they may also be relevant in forensic genetics, since some of the Y chromosome short tandem repeats (Y-STRs) commonly used in forensic analysis are located in amplicons or in inter-amplicon sequences of the AZFc. In a previous study, we identified four phylogenetically related samples with a null allele at DYS448 and a tetrallelic pattern at DYF387S1, two Y-STRs located in the AZFc. Through NGS read depth analysis, we found that the unusual Y-STR pattern may be due to a 1.6 Mb deletion arising concurrently or after a 3.5 Mb duplication event. The observed large genomic rearrangement results in copy number reduction for the RBMY gene family as well as duplication of other AZFc genes. Based on the diversity of 16 additional Y-STRs, we estimated that the duplication/deletion event occurred at least twenty generations ago, suggesting that it has not been affected by negative selection.

Clinically relevant copy-number variants in exome sequencing data of patients with dystonia

IntroductionNext-generation sequencing is now used on a routine basis for molecular testing but studies on copy-number variant (CNV) detection from next-generation sequencing data are underrepresented. Utilizing an existing whole-exome sequencing (WES) dataset, we sought to investigate the contribution of rare CNVs to the genetic causality of dystonia. MethodsThe CNV read-depth analysis tool ExomeDepth was applied to the exome sequences of 953 unrelated patients with dystonia (600 with isolated dystonia and 353 with combined dystonia; 33% with additional neurological involvement). We prioritized rare CNVs that affected known disease genes and/or were known to be associated with defined microdeletion/microduplication syndromes. Pathogenicity assessment of CNVs was based on recently published standards of the American College of Medical Genetics and Genomics and the Clinical Genome Resource. ResultsWe identified pathogenic or likely pathogenic CNVs in 14 of 953 patients (1.5%). Of the 14 different CNVs, 12 were deletions and 2 were duplications, ranging in predicted size from 124bp to 17 Mb. Within the deletion intervals, BRPF1, CHD8, DJ1, EFTUD2, FGF14, GCH1, PANK2, SGCE, UBE3A, VPS16, WARS2, and WDR45 were determined as the most clinically relevant genes. The duplications involved chromosomal regions 6q21-q22 and 15q11-q13. CNV analysis increased the diagnostic yield in the total cohort from 18.4% to 19.8%, as compared to the assessment of single-nucleotide variants and small insertions and deletions alone. ConclusionsWES-based CNV analysis in dystonia is feasible, increases the diagnostic yield, and should be combined with the assessment of single-nucleotide variants and small insertions and deletions.

Nuclear and mitochondrial genome sequencing of North-African Leishmania infantum isolates from cured and relapsed visceral leishmaniasis patients reveals variations correlating with geography and phenotype.

Although several studies have investigated genetic diversity of Leishmania infantum in North Africa, genome-wide analyses are lacking. Here, we conducted comparative analyses of nuclear and mitochondrial genomes of seven L. infantum isolates from Tunisia with the aim to gain insight into factors that drive genomic and phenotypic adaptation. Isolates were from cured (n=4) and recurrent (n=3) visceral leishmaniasis (VL) cases, originating from northern (n=2) and central (n=5) Tunisia, where respectively stable and emerging VL foci are observed. All isolates from relapsed patients were from Kairouan governorate (Centre); one showing resistance to the anti-leishmanial drug Meglumine antimoniate. Nuclear genome diversity of the isolates was analysed by comparison to the L. infantum JPCM5 reference genome. Kinetoplast maxi and minicircle sequences (1 and 59, respectively) were extracted from unmapped reads and identified by blast analysis against public data sets. The genome variation analysis grouped together isolates from the same geographical origins. Strains from the North were very different from the reference showing more than 34 587 specific single nucleotide variants, with one isolate representing a full genetic hybrid as judged by variant frequency. Composition of minicircle classes within isolates corroborated this geographical population structure. Read depth analysis revealed several significant gene copy number variations correlating with either geographical origin (amastin and Hsp33 genes) or relapse (CLN3 gene). However, no specific gene copy number variation was found in the drug-resistant isolate. In contrast, resistance was associated with a specific minicircle pattern suggesting Leishmania mitochondrial DNA as a potential novel source for biomarker discovery.

Deletion of porcine BOLL is associated with defective acrosomes and subfertility in Yorkshire boars.

A recessive sperm defect of Yorkshire boars was detected more than a decade ago. Affected boars produce ejaculates that contain spermatozoa with defective acrosomes, resulting in low fertility. The acrosome defect was mapped to porcine chromosome 15 but the causal mutation has not been identified. We re-analyzed microarray-derived genotypes of affected boars and confirmed that the acrosome defect maps to a 12.24Mb segment of porcine chromosome 15. To detect the mutation causing defective acrosomes, we sequenced the genomes of two affected and three unaffected boars to an average coverage of 11-fold. Read depth analysis revealed a 55kb deletion that is associated with the acrosome defect. The deletion encompasses the BOLL gene encoding the boule homolog, an RNA binding protein which is an evolutionarily conserved member of the DAZ (Deleted in AZoospermia) gene family. Lack of BOLL expression causes spermatogenic arrest and sperm maturation failure in many species. Boars that carry the deletion in the homozygous state produce sperm but their acrosomes are defective, suggesting that lack of porcine BOLL compromises acrosome formation. Our findings warrant further research to investigate the role of BOLL during spermatogenesis and sperm maturation in pigs.

Long-read nanopore sequencing resolves a TMEM231 gene conversion event causing Meckel-Gruber syndrome.

The diagnostic deployment of massively parallel short-read next-generation sequencing (NGS) has greatly improved genetic test availability, speed, and diagnostic yield, particularly for rare inherited disorders. Nonetheless, diagnostic approaches based on short-read sequencing have a poor ability to accurately detect gene conversion events. We report on the genetic analysis of a family in which 3 fetuses had clinical features consistent with the autosomal recessive disorder Meckel-Gruber syndrome (MKS). Targeted NGS of 29 known MKS-associated genes revealed a heterozygous TMEM231 splice donor variant c.929+1A>G. Comparative read-depth analysis, performed to identify a second pathogenic allele, revealed an apparent heterozygous deletion of TMEM231 exon 4. To verify this result we performed single-molecule long-read sequencing of a long-range polymerase chain reaction product spanning this locus. We identified four missense variants that were absent from the short-read dataset due to the preferential mapping of variant-containing reads to a downstream TMEM231 pseudogene. Consistent with the parental segregation analysis, we demonstrate that the single-molecule long reads could be used to show that the variants are arranged in trans. Our experience shows that robust validation of apparent dosage variants remains essential to avoid the pitfalls of short-read sequencing and that new third-generation long-read sequencing technologies can already aid routine clinical care.

A clinically validated whole genome pipeline for structural variant detection and analysis

BackgroundWith the continuing decrease in cost of whole genome sequencing (WGS), we have already approached the point of inflection where WGS testing has become economically feasible, facilitating broader access to the benefits that are helping to define WGS as the new diagnostic standard. WGS provides unique opportunities for detection of structural variants; however, such analyses, despite being recognized by the research community, have not previously made their way into routine clinical practice.ResultsWe have developed a clinically validated pipeline for highly specific and sensitive detection of structural variants basing on 30X PCR-free WGS. Using a combination of breakpoint analysis of split and discordant reads, and read depth analysis, the pipeline identifies structural variants down to single base pair resolution. False positives are minimized using calculations for loss of heterozygosity and bi-modal heterozygous variant allele frequencies to enhance heterozygous deletion and duplication detection respectively. Compound and potential compound combinations of structural variants and small sequence changes are automatically detected. To facilitate clinical interpretation, identified variants are annotated with phenotype information derived from HGMD Professional and population allele frequencies derived from public and Variantyx allele frequency databases. Single base pair resolution enables easy visual inspection of potentially causal variants using the IGV genome browser as well as easy biochemical validation via PCR. Analytical and clinical sensitivity and specificity of the pipeline has been validated using analysis of Genome in a Bottle reference genomes and known positive samples confirmed by orthogonal sequencing technologies.ConclusionConsistent read depth of PCR-free WGS enables reliable detection of structural variants of any size. Annotation both on gene and variant level allows clinicians to match reported patient phenotype with detected variants and confidently report causative finding in all clinical cases used for validation.

Copy number determination of the gene for the human pancreatic polypeptide receptor NPY4R using read depth analysis and droplet digital PCR

BackgroundCopy number variation (CNV) plays an important role in human genetic diversity and has been associated with multiple complex disorders. Here we investigate a CNV on chromosome 10q11.22 that spans NPY4R, the gene for the appetite-regulating pancreatic polypeptide receptor Y4. This genomic region has been challenging to map due to multiple repeated elements and its precise organization has not yet been resolved. Previous studies using microarrays were interpreted to show that the most common copy number was 2 per genome.ResultsWe have investigated 18 individuals from the 1000 Genomes project using the well-established method of read depth analysis and the new droplet digital PCR (ddPCR) method. We find that the most common copy number for NPY4R is 4. The estimated number of copies ranged from three to seven based on read depth analyses with Control-FREEC and CNVnator, and from four to seven based on ddPCR. We suggest that the difference between our results and those published previously can be explained by methodological differences such as reference gene choice, data normalization and method reliability. Three high-quality archaic human genomes (two Neanderthal and one Denisova) display four copies of the NPY4R gene indicating that a duplication occurred prior to the human-Neanderthal/Denisova split.ConclusionsWe conclude that ddPCR is a sensitive and reliable method for CNV determination, that it can be used for read depth calibration in CNV studies based on already available whole-genome sequencing data, and that further investigation of NPY4R copy number variation and its consequences are necessary due to the role of Y4 receptor in food intake regulation.

A Review of 1000 Molecular Investigations of Rare Inherited Anaemia and Related Conditions with the Red Cell Gene Panel

The genetic diagnosis of inherited anaemias is an important aspect of the diagnostic pathway for patients with haematological disorders, allowing discrimination between conditions of overlapping phenotypes therefore enabling more effective clinical treatment. Next Generation sequencing platforms are now in widespread use in diagnostic settings and are facilitating more rapid, accurate and cost-effective molecular diagnosis. The Red Cell Gene Panel developed by the Viapath Molecular Pathology laboratory based at King's College Hospital, London has harnessed this technology with the aim of identifying genetic diagnoses of rare inherited causes of anaemia. Although originally setup to diagnose inherited red cell disorders, clinical demand has led to the inclusion of inherited bone marrow failure syndromes and other related conditions such that the panel now consists of 194 genes, divided into 16 subpanels (see table 1).Here we present the data from the first 1000 diagnostic cases reported under the following referral groups: 462 cases of unexplained anaemia (including haemolytic anaemia, sideroblastic anaemia, congenital dyserythropoietic anaemia, Diamond-Blackfan Anaemia), 232 cases of inherited bone marrow failure syndromes (including thrombocytopenia and neutropenia), 163 cases of congenital erythrocytosis and 143 other cases (including but not limited to iron regulation, haemophagocytic lymphohistiocytosis (HLH) and Criggler-Najjar ). Of these 1000 cases, we have achieved an overall diagnostic yield of approximately 25%. A diagnosed case is defined here as one in which a clear pathogenic or likely pathogenic variant that explains the phenotype has been detected. The unexplained anaemia cases have achieved the highest percentage of cases diagnosed with 47% diagnostic yield and data will be presented outlining the gene-by-gene breakdown of diagnoses made.Our bespoke bioinformatics pipeline has also allowed the detection of novel disease-causing structural variants in 20 cases, contributing 2% of our diagnostic yield. These are detected using three different methods; read-depth analysis, split-read mapping and discordant insert-size analysis. All reported structural variants have been confirmed with a second method, either breakpoint mapping or dosage-sensitive PCR.A significant proportion of cases (28%) have been reported with variants of uncertain clinical significance, highlighting the need for family studies and functional characterisation to be able to accurately ascertain the significance of these variants. Future developments of the service include functional characterisation of membrane disorders using next generation ektacytometry and preliminary data from this work will be presented here. [Display omitted] DisclosuresKulasekararaj:Alexion Pharmaceuticals, Inc.: Consultancy, Honoraria, Other: Travel Support . Pagliuca:Gentium: Honoraria, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau; Jazz Pharmaceuticals: Honoraria, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau.

Comparison of village dog and wolf genomes highlights the role of the neural crest in dog domestication

BackgroundDomesticated from gray wolves between 10 and 40 kya in Eurasia, dogs display a vast array of phenotypes that differ from their ancestors, yet mirror other domesticated animal species, a phenomenon known as the domestication syndrome. Here, we use signatures persisting in dog genomes to identify genes and pathways possibly altered by the selective pressures of domestication.ResultsWhole-genome SNP analyses of 43 globally distributed village dogs and 10 wolves differentiated signatures resulting from domestication rather than breed formation. We identified 246 candidate domestication regions containing 10.8 Mb of genome sequence and 429 genes. The regions share haplotypes with ancient dogs, suggesting that the detected signals are not the result of recent selection. Gene enrichments highlight numerous genes linked to neural crest and central nervous system development as well as neurological function. Read depth analysis suggests that copy number variation played a minor role in dog domestication.ConclusionsOur results identify genes that act early in embryogenesis and can confer phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development as the targets of selection during domestication. These differences reflect the phenotypes of the domestication syndrome, which can be explained by alterations in the migration or activity of neural crest cells during development. We propose that initial selection during early dog domestication was for behavior, a trait influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.

Read Depth Analysis to Identify CNV in Bacteria Using CNOGpro.

Whole-genome sequencing with short-read technologies is well suited for calling single nucleotide polymorphisms, but has major problems with the detection of structural variants larger than the read length. One such type of variation is copy number variation (CNV), which entails deletion or duplication of genomic regions, and the expansion or contraction of repeated elements. Duplicated and deleted regions will typically be collapsed during de novo assembly of sequence data, or ignored when mapping reads toward a reference. However, signatures of the copy number variation can be detected in the resultant read depth at each position in the genome. We here provide instructions on how to analyze this read depth signal with the R package CNOGpro, allowing for estimation of copy numbers with uncertainty for each feature in a genome.

DNA Mismatch Repair Deficiency Promotes Genomic Instability in a Subset of Papillary Thyroid Cancers.

Efficient DNA damage repair by MutL-homolog DNA mismatch repair (MMR) enzymes, MLH1, MLH3, PMS1 and PMS2, are required to maintain thyrocyte genomic integrity. We hypothesized that persistent oxidative stress and consequent transcriptional dysregulation observed in thyroid follicles will lead to MMR deficiency and potentiate papillary thyroid tumorigenesis. MMR gene expression was analyzed by targeted microarray in 18 papillary thyroid cancer (PTC), 9 paracarcinoma normal thyroid (PCNT) and 10 normal thyroid (NT) samples. The findings were validated by qRT-PCR, and in follicular thyroid cancers (FTC) and follicular thyroid adenomas (FTA) for comparison. FOXO transcription factor expression was also analyzed. Protein expression was assessed by immunohistochemistry. Genomic integrity was evaluated by whole-exome sequencing-derived read-depth analysis and Mann-Whitney U test. Clinical correlations were assessed using Fisher's exact and t tests. Microarray and qRT-PCR revealed reduced expression of all four MMR genes in PTC compared with PCNT and of PMS2 compared with NT. FTC and FTA showed upregulation in MLH1, MLH3 and PMS2. PMS2 protein expression correlated with the mRNA expression pattern. FOXO1 showed lower expression in PMS2-deficient PTCs (log2-fold change -1.72 vs. -0.55, U=11, p<0.05 two-tailed). Rate of LOH, a measure of genomic instability, was higher in PMS2-deficient PTCs (median 3 and 1, respectively; U=26, p<0.05 two-tailed). No correlation was noted between MMR deficiency and clinical characteristics. MMR deficiency, potentially promoted by FOXO1 suppression, may explain the etiology for PTC development in some patients. FTC and FTA retain MMR activity and are likely caused by a different tumorigenic pathway.

An enhanced method for targeted next generation sequencing copy number variant detection using ExomeDepth

Copy number variants (CNV) are a major cause of disease, with over 30,000 reported in the DECIPHER database. To use read depth data from targeted Next Generation Sequencing (NGS) panels to identify CNVs with the highest degree of sensitivity, it is necessary to account for biases inherent in the data. GC content and ambiguous mapping due to repetitive sequence elements and pseudogenes are the principal components of technical variability. In addition, the algorithms used favour the detection of multi-exon CNVs, and rely on suitably matched normal dosage samples for comparison. We developed a calling strategy that subdivides target intervals, and uses pools of historical control samples to overcome these limitations in a clinical diagnostic laboratory. We compared our enhanced strategy with an unmodified pipeline using the R software package ExomeDepth, using a cohort of 109 heterozygous CNVs (91 deletions, 18 duplications in 26 genes), including 25 single exon CNVs. The unmodified pipeline detected 104/109 CNVs, giving a sensitivity of 89.62% to 98.49% at the 95% confidence interval. The detection of all 109 CNVs by our enhanced method demonstrates 95% confidence the sensitivity is ≥96.67%, allowing NGS read depth analysis to be used for CNV detection in a clinical diagnostic setting.

Abstract 397: Whole genome copy number variation analysis using a SNP-focused targeted sequencing panel for tumor analysis

Abstract Accurate genome-wide copy number variation (CNV) analysis is critical for disease and cancer research. Current approaches for CNV analysis include fluorescence in situ hybridization (FISH), array comparative genomic hybridization (array CGH), and SNP arrays. Unfortunately, these methods are not sensitive enough for real world cancer samples because of tumor ploidy, purity and heterogeneity. NGS-based targeted sequencing is increasingly being used for CNV analysis due to throughput, coverage, cost, and sample input requirements. For CNV analysis, detection power is improved by combining both read depth and SNP allele frequency analysis, particularly for copy-neutral events such as loss of heterozygosity. A custom xGen Lockdown CNV backbone panel was developed for broad, uniform genome coverage and to enrich for population-based SNPs. We demonstrate use of the panel as an addition to the xGen Exome Research Panel and a custom cancer focused panel. Downstream analysis incorporates both read depth and observed minor allele frequencies to determine CNVs with enhanced sensitivity. To increase the resolution for large-scale alterations of chromosome 7, a hot-spot for disease-associated CNVs, probe density was increased 6 fold. A known standard, NA12878, was used to validate the panel’s ability to detect heterozygous SNPs with high confidence. In addition, mixtures of cancer cell lines from the Cancer Cell Line Encyclopedia (CCLE) were tested with varying levels of background copy-neutral genomic DNA. The sensitivity and specificity of the panel to detect CNV and LOH events with was assessed using deep exome and Affymetrix SNP array data. The ability to detect copy number alterations with high resolution and accuracy would be a valuable resource for disease and cancer research. Citation Format: Jiashi Wang, Kristina Giorda, Zhongwu Lai, Daniel Stetson, Mirna Jarosz. Whole genome copy number variation analysis using a SNP-focused targeted sequencing panel for tumor analysis [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 397. doi:10.1158/1538-7445.AM2017-397

Machine-learning-aided precise prediction of deletions with next-generation sequencing

When detecting deletions in complex human genomes, split-read approaches using short reads generated with next-generation sequencing still face the challenge that either false discovery rate is high, or sensitivity is low. To address the problem, an integrated strategy is proposed. It organically combines the fundamental theories of the three mainstream methods (read-pair approaches, split-read technologies and read-depth analysis) with modern machine learning algorithms, using the recipe of feature extraction as a bridge. Compared with the state-of-art split-read methods for deletion detection in both low and high sequence coverage, the machine-learning-aided strategy shows great ability in intelligently balancing sensitivity and false discovery rate and getting a both more sensitive and more precise call set at single-base-pair resolution. Thus, users do not need to rely on former experience to make an unnecessary trade-off beforehand and adjust parameters over and over again any more. It should be noted that modern machine learning models can play an important role in the field of structural variation prediction.

Statistical models for DNA copy number variation detection using read-depth data from next generation sequencing experiments

Summary In this ‘Big Data’ era, statisticians inevitably encounter data generated from various disciplines. In particular, advances in bio-technology have enabled scientists to produce enormous datasets in various biological experiments. In the last two decades, we have seen high-throughput microarray data resulting from various genomic studies. Recently, next generation sequencing (NGS) technology has been playing an important role in the study of genomic features, resulting in vast amount of NGS data. One frequent application of NGS technology is in the study of DNA copy number variants (CNVs). The resulting NGS read count data are then used by researchers to formulate their various scientific approaches to accurately detect CNVs. Computational and statistical approaches to the detection of CNVs using NGS data are, however, very limited at present. In this review paper, we will focus on read-depth analysis in CNV detection and give a brief summary of currently used statistical analysis methods in searching for CNVs using NGS data. In addition, based on the review, we discuss the challenges we face and future research directions. The ultimate goal of this review paper is to give a timely exposition of the surveyed statistical methods to researchers in related fields.

Large Deletions at the SHOX Locus in the Pseudoautosomal Region Are Associated with Skeletal Atavism in Shetland Ponies.

Skeletal atavism in Shetland ponies is a heritable disorder characterized by abnormal growth of the ulna and fibula that extend the carpal and tarsal joints, respectively. This causes abnormal skeletal structure and impaired movements, and affected foals are usually killed. In order to identify the causal mutation we subjected six confirmed Swedish cases and a DNA pool consisting of 21 control individuals to whole genome resequencing. We screened for polymorphisms where the cases and the control pool were fixed for opposite alleles and observed this signature for only 25 SNPs, most of which were scattered on genome assembly unassigned scaffolds. Read depth analysis at these loci revealed homozygosity or compound heterozygosity for two partially overlapping large deletions in the pseudoautosomal region (PAR) of chromosome X/Y in cases but not in the control pool. One of these deletions removes the entire coding region of the SHOX gene and both deletions remove parts of the CRLF2 gene located downstream of SHOX. The horse reference assembly of the PAR is highly fragmented, and in order to characterize this region we sequenced bacterial artificial chromosome (BAC) clones by single-molecule real-time (SMRT) sequencing technology. This considerably improved the assembly and enabled size estimations of the two deletions to 160−180 kb and 60−80 kb, respectively. Complete association between the presence of these deletions and disease status was verified in eight other affected horses. The result of the present study is consistent with previous studies in humans showing crucial importance of SHOX for normal skeletal development.