Understanding genetic variability: exploring large-scale copy number variants through non-invasive prenatal testing in European populations.

Large-scale copy number variants (CNVs) are structural alterations in the genome that involve the duplication or deletion of DNA segments, contributing to genetic diversity and playing a crucial role in the evolution and development of various diseases and disorders, as they can lead to the dosage imbalance of one or more genes. Massively parallel sequencing (MPS) has revolutionized the field of genetic analysis and contributed significantly to routine clinical diagnosis and screening. It offers a precise method for detecting CNVs with exceptional accuracy. In this context, a non-invasive prenatal test (NIPT) based on the sequencing of cell-free DNA (cfDNA) from pregnant women's plasma using a low-coverage whole genome MPS (WGS) approach represents a valuable source for population studies. Here, we analyzed genomic data of 12,732 pregnant women from the Slovak (9,230), Czech (1,583), and Hungarian (1,919) populations. We identified 5,062 CNVs ranging from 200 kbp and described their basic characteristics and differences between the subject populations. Our results suggest that re-analysis of sequencing data from routine WGS assays has the potential to obtain large-scale CNV population frequencies, which are not well known and may provide valuable information to support the classification and interpretation of this type of genetic variation. Furthermore, this could contribute to expanding knowledge about the central European genome without investing in additional laboratory work, as NIPTs are a relatively widely used screening method.

Association between glutathione S-transferases M1 expression and treatment outcome in germ cell tumor patients.

Cisplatin-based chemotherapy is the mainstay in the treatment of germ cell tumors (GCTs). Glutathione S-transferases (GSTs) are polymorphic enzymes that catalyze the glutathione conjugation of alkylating agents, platinum compounds, and free radicals formed by chemotherapy and are thus implicated in developing treatment resistance. This study aimed to assess the expression level of GST mu 1 (GSTM1) and its association with treatment outcomes in patients with GCT. This translational study included tumor specimens from 207 patients with newly diagnosed GCTs, as well as cisplatin-sensitive GCT cell line xenografts and their resistant variants for all histological variants of GCTs. GSTM1 expression was detected by reverse transcription-quantitative PCR and immunohistochemistry using monoclonal antibodies, scored by the multiplicative quickscore (QS) method. GSTM1 expression was correlated with patient/tumor characteristics and treatment outcomes. The highest GSTM1 expression was observed in seminoma, followed by choriocarcinoma, embryonal carcinoma, and yolk sac tumor, while the lowest was observed in teratoma (p<0.0001). There was no association between GSTM1 expression in tumor tissue and patient/tumor characteristics. The low GSTM1 expression was associated with significantly better relapse-free survival compared with high GSTM1 (HR=0.50, 95% CI 0.23-1.09, p=0.03) but not overall survival (HR=0.61, 95% CI 0.24-1.54, p=0.22). Multivariate analysis showed that the prognostic value of GSTM1 was independent of the International Germ Cell Cancer Collaborative Group (IGCCCG) score. These data revealed the prognostic value of GSTM1 in GCTs, with a high GSTM1 expression associated with worse outcomes, suggesting that GSTM1 could be responsible, in part, for treatment resistance in GCTs.

The Potential of Liquid Biopsy in Detection of Endometrial Cancer Biomarkers: A Pilot Study.

Endometrial cancer belongs to the most common gynecologic cancer types globally, with increasing incidence. There are numerous ways of classifying different cases. The most recent decade has brought advances in molecular classification, which show more accurate prognostic factors and the possibility of personalised adjuvant treatment. In addition, diagnostic approaches lag behind these advances, with methods causing patients discomfort while lacking the reproducibility of tissue sampling for biopsy. Minimally invasive liquid biopsies could therefore represent an alternative screening and diagnostic approach in patients with endometrial cancer. The method could potentially detect molecular changes in this cancer type and identify patients at early stages. In this pilot study, we tested such a detection method based on circulating tumour DNA isolated from the peripheral blood plasma of 21 Slovak endometrial cancer patients. We successfully detected oncomutations in the circulating DNA of every single patient, although the prognostic value of the detected mutations failed to offer certainty. Furthermore, we detected changes associated with clonal hematopoiesis, including DNMT3A mutations, which were present in the majority of circulating tumour DNA samples.

Molecularly confirmed pontocerebellar hypoplasia in a large family from Slovakia with four severely affected children.

BACKGROUND: Pontocerebellar hypoplasia type 1 (PCH1) is characterized by a central and peripheral motor dysfunction associated with anterior horn cell degeneration, similar to spinal muscular atrophy (SMA). OBJECTIVES: We analysed three probands (later discovered to be siblings) suspected to have severe SMA, however, not confirmed by genetic test. METHODS: Clinical-exome analysis (Illumina) was performed to identify causative variants, followed by Sanger sequencing confirmation in probands and other 10 family members. RESULTS: Homozygous pathogenic variant c.92G>C (p.(Gly31Ala)) in the Exosome Component 3 (EXOSC3) gene was found in all 3 probands, thus confirming the diagnosis of a severe form of PCH1B. The parents and six siblings were carriers, while one sibling was homozygous for a reference allele. This variant is frequent in the Czech Roma population, where it is considered a founder mutation. Haplotype analysis in this largest reported PCH1B family showed that our patients inherited from their father (of Roma origin) a haplotype identical to that found in the Czech Roma population, thus indicating these alleles have a common origin. CONCLUSION: This EXOSC3 variant is rare among the general population but most likely frequent also among Roma people in Slovakia. PCH1B should be considered for a differential diagnosis in infants manifesting SMA-like phenotype, especially if of Roma origin (Tab. 1, Fig. 1, Ref. 22). Text in PDF www.elis.sk Keywords: pontocerebellar hypoplasia, PCH1B, EXOSC3, SMA plus syndromes, pathogenic sequence variant.

Analysis of transcriptomics data from COVID-19 patients: a pilot research.

During SARS-CoV-2 infection, the virus transforms the infected host cell into factories that produce new viral particles. As infection progresses, the infected cells undergo numerous changes in various pathways. One of these changes is the occurrence of a cytokine storm, which leads to severe symptoms. In this study, we examined the transcriptomic changes caused by COVID-19 by analyzing RNA-seq data obtained from COVID-19-positive patients as well as COVID-19-negative donors. RNA-seq data were collected for the purpose of identification of potential biomarkers associated with a different course of the disease. We analyzed the first datasets, consisting of 96 samples to validate our methods. The objective of this publication is to report the pilot results. To explore potential biomarkers related to disease severity, we conducted a differential expression analysis of human transcriptome, focusing on COVID-19 positivity and symptom severity. Given the large number of potential biomarkers we identified, we further performed pathway enrichment analysis with terms from Kyoto Encyclopedia of Genes and Genomics (KEGG) to obtain a more profound understanding of altered pathways. Our results indicate that pathways related to immune processes, response to infection, and multiple signaling pathways were affected. These findings align with several previous studies that also reported the influence of SARS-CoV-2 infection on these pathways.

Insights into non-informative results from non-invasive prenatal screening through gestational age, maternal BMI, and age analyses.

The discovery of cell-free fetal DNA fragments in the maternal plasma initiated a novel testing method in prenatal care, called non-invasive prenatal screening (NIPS). One of the limitations of NIPS is the necessity for a sufficient proportion of fetal fragments in the analyzed circulating DNA mixture (fetal fraction), otherwise, the sample is uninterpretable. We present the effect of gestational age, maternal body mass index (BMI), and maternal age on the fetal fraction (FF) of the sample. We retrospectively analyzed data from 5543 pregnant women with a single male fetus who underwent NIPS from which 189 samples received a repeat testing due to an insufficient FF. We showed the relationship between the failure rate of the samples after the repeated analysis, the FF, and the gestational age at the first sampling. Next, we found that different maternal BMI categories affect the FF and thus the chance of an informative redraw. A better understanding of the factors affecting the FF will reduce the number of non-informative calls from repeated analyzes. In this study, we provide helpful information to clinicians on how to approach non-informative analyses.

Association Between Cytometric Biomarkers, Clinical Phenotype, and Complications of Common Variable Immunodeficiency.

BACKGROUND: Common variable immunodeficiency (CVID) is a heterogeneous group of immune disorders. The patients are classified according to the clinical manifestation with the infection-only phenotype (CVIDinf) and CVID with immune dysregulation (CVIDid). METHODS: We performed a retrospective clinical analysis of 64 CVID patients (34 males, 53.13%; mean age: 41.4 years; SD: ±21.4 years). We divided the patients into subgroups according to the clinical manifestation (CVIDinf and CVIDid) and according to B cell phenotypic profiling after performing flow cytometry with the use of the EUROclass classification. We compared clinical manifestations, selected laboratory parameters, and therapy in these groups. All CVIDid patients were tested after the manifestation of complications associated with immune dysregulation and in eight patients during the immunosuppressive treatment (systemic corticosteroids and hydroxychloroquine). RESULTS: Two-thirds of patients in our cohort had symptoms resulting from immune dysregulation. Almost half of the patients had autoimmune complications. A higher proportion of marginal zone B cells was associated with autoimmune complications. A lower percentage of naïve B cells was connected to autoimmunity, whereas a lower proportion of transitional B cells was associated with rheumatic diseases and splenomegaly. Patients with lymphadenopathy had a higher percentage of double-negative T cells and a lower percentage of switched memory B cells. We performed molecular-genetic testing in 28% (n = 17) of patients and found a causal pathogenic variant in 23.5% (n = 4) of this group. CONCLUSION: Based on our results, there is an association between specific cytometric parameters, clinical phenotype, and complications of CVID. The use of the subpopulations of B cells can be helpful in the diagnosis of these specific clinical complications in CVID patients and could help to personalise the therapeutic approach.

Comparison of the ABC and ACMG systems for variant classification.

The ABC and ACMG variant classification systems were compared by asking mainly European clinical laboratories to classify variants in 10 challenging cases using both systems, and to state if the variant in question would be reported as a relevant result or not as a measure of clinical utility. In contrast to the ABC system, the ACMG system was not made to guide variant reporting but to determine the likelihood of pathogenicity. Nevertheless, this comparison is justified since the ACMG class determines variant reporting in many laboratories. Forty-three laboratories participated in the survey. In seven cases, the classification system used did not influence the reporting likelihood when variants labeled as "maybe report" after ACMG-based classification were included. In three cases of population frequent but disease-associated variants, there was a difference in favor of reporting after ABC classification. A possible reason is that ABC step C (standard variant comments) allows a variant to be reported in one clinical setting but not another, e.g., based on Bayesian-based likelihood calculation of clinical relevance. Finally, the selection of ACMG criteria was compared between 36 laboratories. When excluding criteria used by less than four laboratories (<10%), the average concordance rate was 46%. Taken together, ABC-based classification is more clear-cut than ACMG-based classification since molecular and clinical information is handled separately, and variant reporting can be adapted to the clinical question and phenotype. Furthermore, variants do not get a clinically inappropriate label, like pathogenic when not pathogenic in a clinical context, or variant of unknown significance when the significance is known.

Metatranscriptome Analysis of Nasopharyngeal Swabs across the Varying Severity of COVID-19 Disease Demonstrated Unprecedented Species Diversity.

The recent global emergence of the SARS-CoV-2 pandemic has accelerated research in several areas of science whose valuable outputs and findings can help to address future health challenges in the event of emerging infectious agents. We conducted a comprehensive shotgun analysis targeting multiple aspects to compare differences in bacterial spectrum and viral presence through culture-independent RNA sequencing. We conducted a comparative analysis of the microbiome between healthy individuals and those with varying degrees of COVID-19 severity, including a total of 151 participants. Our findings revealed a noteworthy increase in microbial species diversity among patients with COVID-19, irrespective of disease severity. Specifically, our analysis revealed a significant difference in the abundance of bacterial phyla between healthy individuals and those infected with COVID-19. We found that Actinobacteria, among other bacterial phyla, showed a notably higher abundance in healthy individuals compared to infected individuals. Conversely, Bacteroides showed a lower abundance in the latter group. Infected people, regardless of severity and symptoms, have the same proportional representation of Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, and Fusobacteriales. In addition to SARS-CoV-2 and numerous phage groups, we identified sequences of clinically significant viruses such as Human Herpes Virus 1, Human Mastadenovirus D, and Rhinovirus A in several samples. Analyses were performed retrospectively, therefore, in the case of SARS-CoV-2 various WHO variants such as Alpha (B.1.1.7), Delta (B.1.617.2), Omicron (B.1.1.529), and 20C strains are represented. Additionally, the presence of specific virus strains has a certain effect on the distribution of individual microbial taxa.

Long COVID Complicated by Fatal Cytomegalovirus and Aspergillus Infection of the Lungs: An Autopsy Case Report.

After the acute phase of COVID-19, some patients develop long COVID. This term is used for a variety of conditions with a complex, yet not fully elucidated etiology, likely including the prolonged persistence of the virus in the organism and progression to lung fibrosis. We present a unique autopsy case of a patient with severe COVID-19 with prolonged viral persistence who developed interstitial lung fibrosis complicated by a fatal combination of cytomegalovirus and Aspergillus infection. SARS-CoV-2 virus was detected at autopsy in the lungs more than two months after the acute infection, although tests from the nasopharynx were negative. Immune dysregulation after COVID-19 and the administration of corticoid therapy created favorable conditions for the cytomegalovirus and Aspergillus infection that were uncovered at autopsy. These pathogens may represent a risk for opportunistic infections, complicating not only the acute coronavirus infection but also long COVID, as was documented in the presented case.

Combination of expert guidelines-based and machine learning-based approaches leads to superior accuracy of automated prediction of clinical effect of copy number variations.

Clinical interpretation of copy number variants (CNVs) is a complex process that requires skilled clinical professionals. General recommendations have been recently released to guide the CNV interpretation based on predefined criteria to uniform the decision process. Several semiautomatic computational methods have been proposed to recommend appropriate choices, relieving clinicians of tedious searching in vast genomic databases. We have developed and evaluated such a tool called MarCNV and tested it on CNV records collected from the ClinVar database. Alternatively, the emerging machine learning-based tools, such as the recently published ISV (Interpretation of Structural Variants), showed promising ways of even fully automated predictions using broader characterization of affected genomic elements. Such tools utilize features additional to ACMG criteria, thus providing supporting evidence and the potential to improve CNV classification. Since both approaches contribute to evaluation of CNVs clinical impact, we propose a combined solution in the form of a decision support tool based on automated ACMG guidelines (MarCNV) supplemented by a machine learning-based pathogenicity prediction (ISV) for the classification of CNVs. We provide evidence that such a combined approach is able to reduce the number of uncertain classifications and reveal potentially incorrect classifications using automated guidelines. CNV interpretation using MarCNV, ISV, and combined approach is available for non-commercial use at https://predict.genovisio.com/ .

Maternal Copy Number Imbalances in Non-Invasive Prenatal Testing: Do They Matter?

Non-invasive prenatal testing (NIPT) has become a routine practice in screening for common aneuploidies of chromosomes 21, 18, and 13 and gonosomes X and Y in fetuses worldwide since 2015 and has even expanded to include smaller subchromosomal events. In fact, the fetal fraction represents only a small proportion of cell-free DNA on a predominant background of maternal DNA. Unlike fetal findings that have to be confirmed using invasive testing, it has been well documented that NIPT provides information on maternal mosaicism, occult malignancies, and hidden health conditions due to copy number variations (CNVs) with diagnostic resolution. Although large duplications or deletions associated with certain medical conditions or syndromes are usually well recognized and easy to interpret, very little is known about small, relatively common copy number variations on the order of a few hundred kilobases and their potential impact on human health. We analyzed data from 6422 NIPT patient samples with a CNV detection resolution of 200 kb for the maternal genome and identified 942 distinct CNVs; 328 occurred repeatedly. We defined them as multiple occurring variants (MOVs). We scrutinized the most common ones, compared them with frequencies in the gnomAD SVs v2.1, dbVar, and DGV population databases, and analyzed them with an emphasis on genomic content and potential association with specific phenotypes.

Dilated cardiomyopathy is a part of the ARV1-associated phenotype: a case report.

BACKGROUND: ACAT-related enzyme 2 required for viability 1 (ARV1) encodes a transmembrane lipid transporter of the endoplasmic reticulum, which is presented in all eukaryotes and in plants. Deficiency of ARV1 is clinically presented as autosomal recessive developmental and epileptic encephalopathy 38 (DEE38) in humans and in mice. So far, three different homozygous and two compound heterozygous ARV1 mutations in humans have been reported in 15 children. CASE PRESENTATION: In this case report we present a novel homozygous in-frame ARV1-deletion (c.554_556delTAT, p.L185del) in a 21-year old Caucasian man with developmental delay, intellectual disability, seizures, walking and speech impairments, as well as with a dilated cardiomyopathy (DCM), which has not yet been firmly related to the ARV1-associated phenotype. Interestingly, this novel variant lies in the proximity of the p.G189R mutation, which was previously described in two brothers with DEE38 and dilated cardiomyopathy. CONCLUSION: The finding of dilated cardiomyopathy in the presented as well as in three previously reported patients from two different families indicates that dilated cardiomyopathy is a part of the ARV1-induced DEE38 phenotype. However, more data are needed to make this conclusion definitive.

Copy Number Variant Detection with Low-Coverage Whole-Genome Sequencing Represents a Viable Alternative to the Conventional Array-CGH.

Copy number variations (CNVs) represent a type of structural variant involving alterations in the number of copies of specific regions of DNA that can either be deleted or duplicated. CNVs contribute substantially to normal population variability, however, abnormal CNVs cause numerous genetic disorders. At present, several methods for CNV detection are applied, ranging from the conventional cytogenetic analysis, through microarray-based methods (aCGH), to next-generation sequencing (NGS). In this paper, we present GenomeScreen, an NGS-based CNV detection method for low-coverage, whole-genome sequencing. We determined the theoretical limits of its accuracy and obtained confirmation in an extensive in silico study and in real patient samples with known genotypes. In theory, at least 6 M uniquely mapped reads are required to detect a CNV with the length of 100 kilobases (kb) or more with high confidence (Z-score > 7). In practice, the in silico analysis required at least 8 M to obtain >99% accuracy (for 100 kb deviations). We compared GenomeScreen with one of the currently used aCGH methods in diagnostic laboratories, which has mean resolution of 200 kb. GenomeScreen and aCGH both detected 59 deviations, while GenomeScreen furthermore detected 134 other (usually) smaller variations. When compared to aCGH, overall performance of the proposed GenemoScreen tool is comparable or superior in terms of accuracy, turn-around time, and cost-effectiveness, thus providing reasonable benefits, particularly in a prenatal diagnosis setting.

Characterisation of Non-Pathogenic Premutation-Range Myotonic Dystrophy Type 2 Alleles.

Myotonic dystrophy type 2 (DM2) is caused by expansion of a (CCTG)n repeat in the cellular retroviral nucleic acid-binding protein (CNBP) gene. The sequence of the repeat is most commonly interrupted and is stably inherited in the general population. Although expanded alleles, premutation range and, in rare cases, also non-disease associated alleles containing uninterrupted CCTG tracts have been described, the threshold between these categories is poorly characterised. Here, we describe four families with members reporting neuromuscular complaints, in whom we identified altogether nine ambiguous CNBP alleles containing uninterrupted CCTG repeats in the range between 32 and 42 repeats. While these grey-zone alleles are most likely not pathogenic themselves, since other pathogenic mutations were identified and particular family structures did not support their pathogenic role, they were found to be unstable during intergenerational transmission. On the other hand, there was no observable general microsatellite instability in the genome of the carriers of these alleles. Our results further refine the division of CNBP CCTG repeat alleles into two major groups, i.e., interrupted and uninterrupted alleles. Both interrupted and uninterrupted alleles with up to approximately 30 CCTG repeats were shown to be generally stable during intergenerational transmission, while intergenerational as well as somatic instability seems to gradually increase in uninterrupted alleles with tract length growing above this threshold.

Validation of Copy Number Variants Detection from Pregnant Plasma Using Low-Pass Whole-Genome Sequencing in Noninvasive Prenatal Testing-Like Settings.

Detection of copy number variants as an integral part of noninvasive prenatal testing is increasingly used in clinical practice worldwide. We performed validation on plasma samples from 34 pregnant women with known aberrations using cell-free DNA sequencing to evaluate the sensitivity for copy number variants (CNV) detection using an in-house CNV fraction-based detection algorithm. The sensitivity for CNVs smaller than 3 megabases (Mb), larger than 3Mb, and overall was 78.57%, 100%, and 90.6%, respectively. Regarding the fetal fraction, detection sensitivity in the group with a fetal fraction of less than 10% was 57.14%, whereas there was 100% sensitivity in the group with fetal fraction exceeding 10%. The assay is also capable of indicating whether the origin of an aberration is exclusively fetal or fetomaternal/maternal. This validation demonstrated that a CNV fraction-based algorithm was applicable and feasible in clinical settings as a supplement to testing for common trisomies 21, 18, and 13.

Expanding the phenotype of thrombocytopenia absent radius syndrome with hypospadias.

Rare genetic diseases and syndromes may appear with unique features in some patients. In genetically-solved cases, this situation indicates a phenotypic expansion of the syndrome with additional features (i.e. the disease-associated gene gives rise to unusual clinical presentation). However, this situation can also hide a multilocus pathogenic variation that cannot be solved genetically except by a massive sequencing approach, such as exome sequencing. Here we describe the case of a child with bilateral radial aplasia, transient thrombocytopenia and anemia, cow's milk intolerance, hypospadias, facial dysmorphism, mild hypothyroidism and umbilical and inguinal hernia. Bilaterally absent radius, presence of thumbs and low platelet count are pathognomonic of thrombocytopenia absent radius (TAR) syndrome, but the child also showed other features beyond those reported in the literature. Since various diseases resembling the proband's phenotype required differential diagnosis, clinical exome sequencing was performed. The results showed compound heterozygous mutations in the RBM8A gene, confirming the suspicion of TAR syndrome. A truncating heterozygous variant in the DUOX2 gene, known to be associated with transient thyroid dyshormonogenesis type 6 (TDH6), was also detected and may explain the proband's mild hypothyroidism.

Non-invasive prenatal testing (NIPT) by low coverage genomic sequencing: Detection limits of screened chromosomal microdeletions.

To study the detection limits of chromosomal microaberrations in non-invasive prenatal testing with aim for five target microdeletion syndromes, including DiGeorge, Prader-Willi/Angelman, 1p36, Cri-Du-Chat, and Wolf-Hirschhorn syndromes. We used known cases of pathogenic deletions from ISCA database to specifically define regions critical for the target syndromes. Our approach to detect microdeletions, from whole genome sequencing data, is based on sample normalization and read counting for individual bins. We performed both an in-silico study using artificially created data sets and a laboratory test on mixed DNA samples, with known microdeletions, to assess the sensitivity of prediction for varying fetal fractions, deletion lengths, and sequencing read counts. The in-silico study showed sensitivity of 79.3% for 10% fetal fraction with 20M read count, which further increased to 98.4% if we searched only for deletions longer than 3Mb. The test on laboratory-prepared mixed samples was in agreement with in-silico results, while we were able to correctly detect 24 out of 29 control samples. Our results suggest that it is possible to incorporate microaberration detection into basic NIPT as part of the offered screening/diagnostics procedure, however, accuracy and reliability depends on several specific factors.

Spinal muscular atrophy caused by a novel Alu-mediated deletion of exons 2a-5 in SMN1 undetectable with routine genetic testing.

BACKGROUND: Spinal muscular atrophy (SMA) is an inherited neuromuscular disease affecting 1 in 8,000 newborns. The majority of patients carry bi-allelic variants in the survival of motor neuron 1 gene (SMN1). SMN1 is located in a duplicated region on chromosome 5q13 that contains Alu elements and is predisposed to genomic rearrangements. Due to the genomic complexity of the SMN region and genetic heterogeneity, approximately 50% of SMA patients remain without genetic diagnosis that is a prerequisite for genetic treatments. In this work we describe the diagnostic odyssey of one SMA patient in whom routine diagnostics identified only a maternal heterozygous SMN1Δ(7-8) deletion. METHODS: We characterized SMN transcripts, assessed SMN protein content in peripheral blood mononuclear cells (PBMC), estimated SMN genes dosage, and mapped genomic rearrangement in the SMN region. RESULTS: We identified an Alu-mediated deletion encompassing exons 2a-5 of SMN1 on the paternal allele and a complete deletion of SMN1 on the maternal allele as the cause of SMA in this patient. CONCLUSION: Alu-mediated rearrangements in SMN1 can escape routine diagnostic testing. Parallel analysis of SMN gene dosage, SMN transcripts, and total SMN protein levels in PBMC can identify genomic rearrangements and should be considered in genetically undefined SMA cases.

On the critical evaluation and confirmation of germline sequence variants identified using massively parallel sequencing.

Although massively parallel sequencing (MPS) is becoming common practice in both research and routine clinical care, confirmation requirements of identified DNA variants using alternative methods are still topics of debate. When evaluating variants directly from MPS data, different read depth statistics, together with specialized genotype quality scores are, therefore, of high relevance. Here we report results of our validation study performed in two different ways: 1) confirmation of MPS identified variants using Sanger sequencing; and 2) simultaneous Sanger and MPS analysis of exons of selected genes. Detailed examination of false-positive and false-negative findings revealed typical error sources connected to low read depth/coverage, incomplete reference genome, indel realignment problems, as well as microsatellite associated amplification errors leading to base miss-calling. However, all these error types were identifiable with thorough manual revision of aligned reads according to specific patterns of distributions of variants and their corresponding reads. Moreover, our results point to dependence of both basic quantitative metrics (such as total read counts, alternative allele read counts and allelic balance) together with specific genotype quality scores on the used bioinformatics pipeline, stressing thus the need for establishing of specific thresholds for these metrics in each laboratory and for each involved pipeline independently.