Allelic dropout in PAH affecting the results of genetic diagnosis in phenylketonuria

Objectives Phenylketonuria (PKU) is an inherited autosomal recessive disorder of phenylalanine metabolism. It is mainly caused by a deficiency in phenylalanine hydroxylase (PAH) and frequently diagnosed with Sanger sequencing. To some extent, allelic dropout can explain the inconsistency in genotype and phenotype. Methods Three families were evaluated through DNA sequence analysis, multiplex ligation-dependent probe amplification (MLPA) and prenatal diagnosis technologies. The possibility of inconsistency in phenotype and genotype with c.331C>T variant was analysed. Results Through pedigree analysis, three mothers carried a homozygous c.331C>T variant, which was a false-positive result. New primers were used, and this error was caused by allelic dropout. In this case, c.158G>A was likely a benign variant. Conclusions Sequence variants in primer-binding regions could cause allelic dropout, creating unpredictable errors in genotyping. Our results emphasised the need for careful measures to treat genotype–phenotype inconsistencies.

Body Remains Left by Bird Predators as a Reliable Source for Population Genetic Studies in the Great Capricorn Beetle Cerambyx cerdo, a Veteran Oak Specialist

Obtaining biological material for DNA extraction is often challenging in organisms of conservation interest. Non-invasive sampling (i.e., sampling without injuring or disturbing an animal) is preferred as it carries no risk to the population’s survival. Here, we tested the possibility of using the body remains left by bird predators for microsatellite genotyping in Cerambyx cerdo, a veteran oak specialist. We compared results obtained from such potentially degraded samples with samples of fresh beetle tarsi (i.e., invasive and destructive but non-lethal samples). Using 10 SSR loci, we genotyped 28 fresh, and 28 remains samples. The analysis indicated that PCR amplification efficiency was not influenced by sample type but allele length and individual heterozygosity. Allele frequencies were perfectly correlated for both sample types (R2 = 0.94). Additionally, null allele frequencies and genotyping failure rates were not significantly different from zero. Although the point estimates of individual inbreeding rates (fi) were higher in remains than fresh samples (medians 0.08 vs. 0.02, respectively), both groups were not significantly different from each other and zero. Our study demonstrated that non-invasive remains samples could provide satisfactory data for population–genetic studies. However, we highlight the problem of biased inbreeding estimates, which may result from samples affected by allelic dropout.

Assessing genotyping errors in mammalian museum study skins using high-throughput genotyping-by-sequencing

The use of museum specimens held in natural history repositories for population and conservation genetic research is increasing in tandem with the use of massively parallel sequencing technologies. Short Tandem Repeats (STRs), or microsatellite loci, are commonly used genetic markers in wildlife and population genetic studies. However, they traditionally suffered from a host of issues including length homoplasy, high costs, low throughput, and difficulties in reproducibility across laboratories. Massively parallel sequencing technologies can address these problems, but the incorporation of museum specimen derived DNA suffers from significant fragmentation and exogenous DNA contamination. Combatting these issues requires extra measures of stringency in the lab and during data analysis, yet there have not been any high-throughput sequencing studies evaluating microsatellite allelic dropout from museum specimen extracted DNA. In this study, we evaluate genotyping errors derived from mammalian museum skin DNA extracts for previously characterized microsatellites across PCR replicates utilizing high-throughput sequencing. We found it useful to classify samples based on DNA concentration, which determined the rate by which genotypes were accurately recovered. Longer microsatellites performed worse in all museum specimens. Allelic dropout rates across loci were dependent on sample quantity, with high concentration museum specimens performing as well and recovering quality metrics nearly as high as the frozen tissue sample. Based on our results, we provide a set of best practices for quality assurance and incorporation of reliable genotypes from museum specimens.

Comment on “Individual heterozygosity predicts translocation success in threatened desert tortoises”

Scott et al. (Reports, 27 November 2020, p. 1086) suggest, on the basis of conclusions obtained from a desert tortoise reintroduction program, that higher genomic heterozygosity should be used to identify individuals for successful translocation. I contend that this recommendation is questionable given these relocated tortoises’ unknown origin, their high mortality, insufficient data on resident tortoises and other components of fitness, and potential allelic dropout.

MIPP-Seq: ultra-sensitive rapid detection and validation of low-frequency mosaic mutations

Background Mosaic mutations contribute to numerous human disorders. As such, the identification and precise quantification of mosaic mutations is essential for a wide range of research applications, clinical diagnoses, and early detection of cancers. Currently, the low-throughput nature of single allele assays (e.g., allele-specific ddPCR) commonly used for genotyping known mutations at very low alternate allelic fractions (AAFs) have limited the integration of low-level mosaic analyses into clinical and research applications. The growing importance of mosaic mutations requires a more rapid, low-cost solution for mutation detection and validation. Methods To overcome these limitations, we developed Multiple Independent Primer PCR Sequencing (MIPP-Seq) which combines the power of ultra-deep sequencing and truly independent assays. The accuracy of MIPP-seq to quantifiable detect and measure extremely low allelic fractions was assessed using a combination of SNVs, insertions, and deletions at known allelic fractions in blood and brain derived DNA samples. Results The Independent amplicon analyses of MIPP-Seq markedly reduce the impact of allelic dropout, amplification bias, PCR-induced, and sequencing artifacts. Using low DNA inputs of either 25 ng or 50 ng of DNA, MIPP-Seq provides sensitive and quantitative assessments of AAFs as low as 0.025% for SNVs, insertion, and deletions. Conclusions MIPP-Seq provides an ultra-sensitive, low-cost approach for detecting and validating known and novel mutations in a highly scalable system with broad utility spanning both research and clinical diagnostic testing applications. The scalability of MIPP-Seq allows for multiplexing mutations and samples, which dramatically reduce costs of variant validation when compared to methods like ddPCR. By leveraging the power of individual analyses of multiple unique and independent reactions, MIPP-Seq can validate and precisely quantitate extremely low AAFs across multiple tissues and mutational categories including both indels and SNVs. Furthermore, using Illumina sequencing technology, MIPP-seq provides a robust method for accurate detection of novel mutations at an extremely low AAF.

Allelic Dropout Is a Common Phenomenon That Reduces the Diagnostic Yield of PCR-Based Sequencing of Targeted Gene Panels

Primary cardiomyopathies (CMPs) are monogenic but multi-allelic disorders with dozens of genes involved in pathogenesis. The implementation of next-generation sequencing (NGS) approaches has resulted in more time- and cost-efficient DNA diagnostics of cardiomyopathies. However, the diagnostic yield of genetic testing for each subtype of CMP fails to exceed 60%. The aim of this study was to demonstrate that allelic dropout (ADO) is a common phenomenon that reduces the diagnostic yield in primary cardiomyopathy genetic testing based on targeted gene panels assayed on the Ion Torrent platform. We performed mutational screening with three custom targeted gene panels based on sets of oligoprimers designed automatically using AmpliSeq Designer® containing 1049 primer pairs for 37 genes with a total length of 153 kb. DNA samples from 232 patients were screened with at least one of these targeted gene panels. We detected six ADO events in both IonTorrent PGM (three cases) and capillary Sanger sequencing (three cases) data, identifying ADO-causing variants in all cases. All ADO events occurred due to common or rare single nucleotide variants (SNVs) in the oligoprimer binding sites and were detected because of the presence of “marker” SNVs in the target DNA fragment. We ultimately identified that PCR-based NGS involves a risk of ADO that necessitates the use of Sanger sequencing to validate NGS results. We assume that oligoprimer design without ADO data affects the amplification efficiency of up to 0.77% of amplicons.

Factors influencing genotyping success and genotyping error rate of Eurasian otter (Lutra lutra) faeces collected in temperate Central Europe

The use of non-invasively collected DNA source material for genetic and genomic applications is usually characterized by low target DNA concentration and quality, genotyping errors and cost-intensive lab procedures. However, for otters (Lutrinae) as elusive species of conservation concern, genetic non-invasive sampling has become an important tool to study their ecology and demography. To increase cost-efficiency of monitoring programmes and to promote the expansion of genomic approaches to non-invasive samples, we aimed to refine sample collection and preparation. Therefore, we examined the effects of intrinsic sample characteristics (including diet), environmental conditions in the field and sample treatment in the molecular laboratory on the success of genotyping and allelic dropout (ADO) rates using microsatellite markers in 1970 fresh Eurasian otter (Lutra lutra) scats. Using fresh samples only, we probably eliminated one of the most important impediments of genotyping DNA from otter faecal samples beforehand. But, we observed higher genotyping success and lower ADO rates for anal glad secretions and faecal samples containing high proportions of mucus. Moist conditions during sample collection may promote DNA degradation and PCR inhibition, leading to decreased genotyping success rates. ADO was further affected by the type of extraction kit. However, a high proportion of variance remaining unexplained by our models implied that additional parameters were acting (amount of PCR inhibitors, non-uniform distribution of intestinal cells, efficiency of PCRs, specific microclimate at marking sites). We summarized influential factors maximizing genotyping quality of otter scats and give recommendations for sample collection, storage and DNA extraction based on our results and current literature.

Sanger Validation of High-Throughput Sequencing in Genetic Diagnosis: Still the Best Practice?

Next-generation sequencing (NGS)’s crucial role in supporting genetic diagnosis and personalized medicine leads to the definition of Guidelines for Diagnostic NGS by the European Society of Human Genetics. Factors of different nature producing false-positive/negative NGS data together with the paucity of internationally accepted guidelines providing specified NGS quality metrics to be followed for diagnostics purpose made the Sanger validation of NGS variants still mandatory. We reported the analysis of three cases of discrepancy between NGS and Sanger sequencing in a cohort of 218 patients. NGS was performed by Illumina MiSeq® and Haloplex/SureSelect protocols targeting 97 or 57 or 10 gene panels usually applied for diagnostics. Variants called following guidelines suggested by the Broad Institute and identified according to MAF <0.01 and allele balance >0.2 were Sanger validated. Three out of 945 validated variants showed a discrepancy between NGS and Sanger. In all three cases, a deep evaluation of the discrepant gene variant results and methodological approach allowed to confirm the NGS datum. Allelic dropout (ADO) occurrence during polymerase chain or sequencing reaction was observed, mainly related to incorrect variant zygosity. Our study extends literature data in which almost 100% “high quality” NGS variants are confirmed by Sanger; moreover, it demonstrates that in case of discrepancy between a high-quality NGS variant and Sanger validation, NGS call should not be a priori assumed to represent the source of the error. Actually, difficulties (i.e., ADO, unpredictable presence of private variants on primer-binding regions) of the so-called gold standard direct sequencing should be considered especially in light of the constantly implemented and accurate high-throughput technologies. Our data along with literature raise a discussion on the opportunity to establish a standardized quality threshold by International Guidelines for clinical NGS in order to limit Sanger confirmation to borderline conditions of variant quality parameters and verification of correct gene variant call/patient coupling on a different blood sample aliquot.

MHC Genotyping by SSCP and Amplicon-Based NGS Approach in Chamois

Genes of the major histocompatibility complex (MHC) code for cell surface proteins essential for adaptive immunity. They show the most outstanding genetic diversity in vertebrates, which has been connected with various fitness traits and thus with the long-term persistence of populations. In this study, polymorphism of the MHC class II DRB locus was investigated in chamois with Single-Strand Conformation Polymorphism (SSCP)/Sanger genotyping and Ion Torrent S5 next-generation sequencing (NGS). From eight identified DRB variants in 28 individuals, five had already been described, and three were new, undescribed alleles. With conventional SSCP/Sanger sequencing, we were able to detect seven alleles, all of which were also detected with NGS. We found inconsistencies in the individual genotypes between the two methods, which were mainly caused by allelic dropout in the SSCP/Sanger method. Six out of 28 individuals were falsely classified as homozygous with SSCP/Sanger analysis. Overall, 25% of the individuals were identified as genotyping discrepancies between the two methods. Our results show that NGS technologies are better performing in sequencing highly variable regions such as the MHC, and they also have a higher detection capacity, thus allowing a more accurate description of the genetic composition, which is crucial for evolutionary and population genetic studies.

CellPhy: accurate and fast probabilistic inference of single-cell phylogenies from scDNA-seq data

We have developed a maximum likelihood framework called CellPhy for inferring phylogenetic trees from single-cell DNA sequencing (scDNA-seq) data, that can be directly applied to somatic cells and clones. CellPhy is based on a finite-site Markov nucleotide substitution model with 10 diploid states, akin to those typically used in statistical phylogenetics. It includes a dedicated error function for single cells that explicitly incorporates amplification/sequencing error and allelic dropout (ADO). Moreover, it can explicitly consider the uncertainty of the variant calling process by using genotype likelihoods as input. We implemented CellPhy in a widely used open-source phylogenetic inference package (RAxML-NG) that provides statistical confidence measurements on the estimated tree and scales particularly well on large phylogenies with hundreds or even thousands of cells. To benchmark CellPhy, we carried out 19,400 coalescent simulations of cell samples from exponentially-growing tumors for which the true phylogeny was known. We evolved single-cell diploid DNA genotypes along the simulated genealogies under different scenarios including infinite- and finite-sites nucleotide mutation models, trinucleotide mutational signatures, sequencing and amplification errors, allele dropouts, and doublet cells. Our simulations suggest that CellPhy is robust to amplification/sequencing errors and to ADO and that it outperforms the state-of-the-art methods under realistic scDNA-seq scenarios both in terms of accuracy and speed. In addition, we sequenced 24 single-cell whole genomes from a colorectal cancer, and together with three published scDNA-seq data sets, analyzed them to illustrate how CellPhy can provide more reliable biological insights than competing methods. CellPhy is freely available at https://github.com/amkozlov/cellphy.