Validated WGS and WES protocols proved saliva-derived gDNA as an equivalent to blood-derived gDNA for clinical and population genomic analyses.

BACKGROUND: Whole exome sequencing (WES) and whole genome sequencing (WGS) have become standard methods in human clinical diagnostics as well as in population genomics (POPGEN). Blood-derived genomic DNA (gDNA) is routinely used in the clinical environment. Conversely, many POPGEN studies and commercial tests benefit from easy saliva sampling. Here, we evaluated the quality of variant call sets and the level of genotype concordance of single nucleotide variants (SNVs) and small insertions and deletions (indels) for WES and WGS using paired blood- and saliva-derived gDNA isolates employing genomic reference-based validated protocols. METHODS: The genomic reference standard Coriell NA12878 was repeatedly analyzed using optimized WES and WGS protocols, and data calls were compared with the truth dataset published by the Genome in a Bottle Consortium. gDNA was extracted from the paired blood and saliva samples of 10 participants and processed using the same protocols. A comparison of paired blood-saliva call sets was performed in the context of WGS and WES genomic reference-based technical validation results. RESULTS: The quality pattern of called variants obtained from genomic-reference-based technical replicates correlates with data calls of paired blood-saliva-derived samples in all levels of tested examinations despite a higher rate of non-human contamination found in the saliva samples. The F1 score of 10 blood-to-saliva-derived comparisons ranged between 0.8030-0.9998 for SNVs and between 0.8883-0.9991 for small-indels in the case of the WGS protocol, and between 0.8643-0.999 for SNVs and between 0.7781-1.000 for small-indels in the case of the WES protocol. CONCLUSION: Saliva may be considered an equivalent material to blood for genetic analysis for both WGS and WES under strict protocol conditions. The accuracy of sequencing metrics and variant-detection accuracy is not affected by choosing saliva as the gDNA source instead of blood but much more significantly by the genomic context, variant types, and the sequencing technology used.

WarpSTR: determining tandem repeat lengths using raw nanopore signals.

MOTIVATION: Short tandem repeats (STRs) are regions of a genome containing many consecutive copies of the same short motif, possibly with small variations. Analysis of STRs has many clinical uses but is limited by technology mainly due to STRs surpassing the used read length. Nanopore sequencing, as one of long-read sequencing technologies, produces very long reads, thus offering more possibilities to study and analyze STRs. Basecalling of nanopore reads is however particularly unreliable in repeating regions, and therefore direct analysis from raw nanopore data is required. RESULTS: Here, we present WarpSTR, a novel method for characterizing both simple and complex tandem repeats directly from raw nanopore signals using a finite-state automaton and a search algorithm analogous to dynamic time warping. By applying this approach to determine the lengths of 241 STRs, we demonstrate that our approach decreases the mean absolute error of the STR length estimate compared to basecalling and STRique. AVAILABILITY AND IMPLEMENTATION: WarpSTR is freely available at https://github.com/fmfi-compbio/warpstr.

Evaluation and limitations of different approaches among COVID-19 fatal cases using whole-exome sequencing data.

BACKGROUND: COVID-19 caused by the SARS-CoV-2 infection may result in various disease symptoms and severity, ranging from asymptomatic, through mildly symptomatic, up to very severe and even fatal cases. Although environmental, clinical, and social factors play important roles in both susceptibility to the SARS-CoV-2 infection and progress of COVID-19 disease, it is becoming evident that both pathogen and host genetic factors are important too. In this study, we report findings from whole-exome sequencing (WES) of 27 individuals who died due to COVID-19, especially focusing on frequencies of DNA variants in genes previously associated with the SARS-CoV-2 infection and the severity of COVID-19. RESULTS: We selected the risk DNA variants/alleles or target genes using four different approaches: 1) aggregated GWAS results from the GWAS Catalog; 2) selected publications from PubMed; 3) the aggregated results of the Host Genetics Initiative database; and 4) a commercial DNA variant annotation/interpretation tool providing its own knowledgebase. We divided these variants/genes into those reported to influence the susceptibility to the SARS-CoV-2 infection and those influencing the severity of COVID-19. Based on the above, we compared the frequencies of alleles found in the fatal COVID-19 cases to the frequencies identified in two population control datasets (non-Finnish European population from the gnomAD database and genomic frequencies specific for the Slovak population from our own database). When compared to both control population datasets, our analyses indicated a trend of higher frequencies of severe COVID-19 associated risk alleles among fatal COVID-19 cases. This trend reached statistical significance specifically when using the HGI-derived variant list. We also analysed other approaches to WES data evaluation, demonstrating its utility as well as limitations. CONCLUSIONS: Although our results proved the likely involvement of host genetic factors pointed out by previous studies looking into severity of COVID-19 disease, careful considerations of the molecular-testing strategies and the evaluated genomic positions may have a strong impact on the utility of genomic testing.

APOC3 and ABCA1 variants in unusual combined hypolipidaemia showing premature peripheral vascular disease.

BACKGROUND: Familial combined hypolipidaemia is a condition characterised by very low concentrations of circulating very-low-density lipoprotein (VLDL), low-density lipoprotein cholesterol (LDL), and high-density lipoprotein cholesterol (HDL). It is thought that low LDL/combined hypolipidaemia can protect from cardiovascular disease (CVD), but this is not what we found in a case we present. OBJECTIVE: We report on a 57-years-old male patient with combined hypolipidaemia who presented with premature peripheral vascular disease. We investigated also his two sons, 32- and 27-years-old, who manifested a tendency to low lipid levels. METHODS AND RESULTS: We used Illumina exome analysis in all three individuals and in all of them we could exclude the major effect of the variants within the genes most frequently mutated in hypolipidaemia, including recently reported LIPC gene variant. Instead, in all three individuals we identified a novel ABCA1 variant, possibly responsible for the decreased HDL levels. The proband and one of his sons also share the splicing APOC3 variant rs138326449, known to be associated with decreased TG levels. CONCLUSION: The heterogeneous nature and the risk of atherosclerosis in combined hypolipidaemia seems to be variable, based on an interplay between low HDL and LDL levels, and it depends on the combination of variants that cause it (Tab. 2, Ref. 38).

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.

Privacy-preserving storage of sequenced genomic data.

BACKGROUND: The current and future applications of genomic data may raise ethical and privacy concerns. Processing and storing of this data introduce a risk of abuse by potential offenders since the human genome contains sensitive personal information. For this reason, we have developed a privacy-preserving method, named Varlock providing secure storage of sequenced genomic data. We used a public set of population allele frequencies to mask the personal alleles detected in genomic reads. Each personal allele described by the public set is masked by a randomly selected population allele with respect to its frequency. Masked alleles are preserved in an encrypted confidential file that can be shared in whole or in part using public-key cryptography. RESULTS: Our method masked the personal variants and introduced new variants detected in a personal masked genome. Alternative alleles with lower population frequency were masked and introduced more often. We performed a joint PCA analysis of personal and masked VCFs, showing that the VCFs between the two groups cannot be trivially mapped. Moreover, the method is reversible and personal alleles in specific genomic regions can be unmasked on demand. CONCLUSION: Our method masks personal alleles within genomic reads while preserving valuable non-sensitive properties of sequenced DNA fragments for further research. Personal alleles in the desired genomic regions may be restored and shared with patients, clinics, and researchers. We suggest that the method can provide an additional security layer for storing and sharing of the raw aligned reads.

Dante: genotyping of known complex and expanded short tandem repeats.

MOTIVATION: Short tandem repeats (STRs) are stretches of repetitive DNA in which short sequences, typically made of 2-6 nucleotides, are repeated several times. Since STRs have many important biological roles and also belong to the most polymorphic parts of the human genome, they became utilized in several molecular-genetic applications. Precise genotyping of STR alleles, therefore, was of high relevance during the last decades. Despite this, massively parallel sequencing (MPS) still lacks the analysis methods to fully utilize the information value of STRs in genome scale assays. RESULTS: We propose an alignment-free algorithm, called Dante, for genotyping and characterization of STR alleles at user-specified known loci based on sequence reads originating from STR loci of interest. The method accounts for natural deviations from the expected sequence, such as variation in the repeat count, sequencing errors, ambiguous bases and complex loci containing several different motifs. In addition, we implemented a correction for copy number defects caused by the polymerase induced stutter effect as well as a prediction of STR expansions that, according to the conventional view, cannot be fully captured by inherently short MPS reads. We tested Dante on simulated datasets and on datasets obtained by targeted sequencing of protein coding parts of thousands of selected clinically relevant genes. In both these datasets, Dante outperformed HipSTR and GATK genotyping tools. Furthermore, Dante was able to predict allele expansions in all tested clinical cases. AVAILABILITY AND IMPLEMENTATION: Dante is open source software, freely available for download at https://github.com/jbudis/dante. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Combining count- and length-based z-scores leads to improved predictions in non-invasive prenatal testing.

MOTIVATION: Non-invasive prenatal testing or NIPT is currently among the top researched topic in obstetric care. While the performance of the current state-of-the-art NIPT solutions achieve high sensitivity and specificity, they still struggle with a considerable number of samples that cannot be concluded with certainty. Such uninformative results are often subject to repeated blood sampling and re-analysis, usually after two weeks, and this period may cause a stress to the future mothers as well as increase the overall cost of the test. RESULTS: We propose a supplementary method to traditional z-scores to reduce the number of such uninformative calls. The method is based on a novel analysis of the length profile of circulating cell free DNA which compares the change in such profiles when random-based and length-based elimination of some fragments is performed. The proposed method is not as accurate as the standard z-score; however, our results suggest that combination of these two independent methods correctly resolves a substantial portion of healthy samples with an uninformative result. Additionally, we discuss how the proposed method can be used to identify maternal aberrations, thus reducing the risk of false positive and false negative calls. AVAILABILITY AND IMPLEMENTATION: The open-source code of the proposed methods, together with test data, is freely available for non-commercial users at github web page https://github.com/jbudis/lambda. SUPPLEMENTARY INFORMATION: Supplementary materials are available at Bioinformatics online.

Technical and Methodological Aspects of Cell-Free Nucleic Acids Analyzes.

Analyzes of cell-free nucleic acids (cfNAs) have shown huge potential in many biomedical applications, gradually entering several fields of research and everyday clinical care. Many biological properties of cfNAs can be informative to gain deeper insights into the function of the organism, such as their different types (DNA, RNAs) and subtypes (gDNA, mtDNA, bacterial DNA, miRNAs, etc.), forms (naked or vesicle bound NAs), fragmentation profiles, sequence composition, epigenetic modifications, and many others. On the other hand, the workflows of their analyzes comprise many important steps, from sample collection, storage and transportation, through extraction and laboratory analysis, up to bioinformatic analyzes and statistical evaluations, where each of these steps has the potential to affect the outcome and informational value of the performed analyzes. There are, however, no universal or standard protocols on how to exactly proceed when analyzing different cfNAs for different applications, at least according to our best knowledge. We decided therefore to prepare an overview of the available literature and products commercialized for cfNAs processing, in an attempt to summarize the benefits and limitations of the currently available approaches, devices, consumables, and protocols, together with various factors influencing the workflow, its processes, and outcomes.

Cell-Free Nucleic Acids and their Emerging Role in the Pathogenesis and Clinical Management of Inflammatory Bowel Disease.

Cell-free nucleic acids (cfNAs) are defined as any nucleic acids that are present outside the cell. They represent valuable biomarkers in various diagnostic protocols such as prenatal diagnostics, the detection of cancer, and cardiovascular or autoimmune diseases. However, in the current literature, little is known about their implication in inflammatory bowel disease (IBD). IBD is a group of multifactorial, autoimmune, and debilitating diseases with increasing incidence worldwide. Despite extensive research, their etiology and exact pathogenesis is still unclear. Since cfNAs were observed in other autoimmune diseases and appear to be relevant in inflammatory processes, their role in the pathogenesis of IBD has also been suggested. This review provides a summary of knowledge from the available literature about cfDNA and cfRNA and the structures involving them such as exosomes and neutrophil extracellular traps and their association with IBD. Current studies showed the promise of cfNAs in the management of IBD not only as biomarkers distinguishing patients from healthy people and differentiating active from inactive disease state, but also as a potential therapeutic target. However, the detailed biological characteristics of cfNAs need to be fully elucidated in future experimental and clinical studies.

Identification of Structural Variation from NGS-Based Non-Invasive Prenatal Testing.

Copy number variants (CNVs) are an important type of human genome variation, which play a significant role in evolution contribute to population diversity and human genetic diseases. In recent years, next generation sequencing has become a valuable tool for clinical diagnostics and to provide sensitive and accurate approaches for detecting CNVs. In our previous work, we described a non-invasive prenatal test (NIPT) based on low-coverage massively parallel whole-genome sequencing of total plasma DNA for detection of CNV aberrations ≥600 kbp. We reanalyzed NIPT genomic data from 5018 patients to evaluate CNV aberrations in the Slovak population. Our analysis of autosomal chromosomes identified 225 maternal CNVs (47 deletions; 178 duplications) ranging from 600 to 7820 kbp. According to the ClinVar database, 137 CNVs (60.89%) were fully overlapping with previously annotated variants, 66 CNVs (29.33%) were in partial overlap, and 22 CNVs (9.78%) did not overlap with any previously described variant. Identified variants were further classified with the AnnotSV method. In summary, we identified 129 likely benign variants, 13 variants of uncertain significance, and 83 likely pathogenic variants. In this study, we use NIPT as a valuable source of population specific data. Our results suggest the utility of genomic data from commercial CNV analysis test as background for a population study.

Combination of Fetal Fraction Estimators Based on Fragment Lengths and Fragment Counts in Non-Invasive Prenatal Testing.

The reliability of non-invasive prenatal testing is highly dependent on accurate estimation of fetal fraction. Several methods have been proposed up to date, utilizing different attributes of analyzed genomic material, for example length and genomic location of sequenced DNA fragments. These two sources of information are relatively unrelated, but so far, there have been no published attempts to combine them to get an improved predictor. We collected 2454 single euploid male fetus samples from women undergoing NIPT testing. Fetal fractions were calculated using several proposed predictors and the state-of-the-art SeqFF method. Predictions were compared with the reference Y-based method. We demonstrate that prediction based on length of sequenced DNA fragments may achieve nearly the same precision as the state-of-the-art methods based on their genomic locations. We also show that combination of several sample attributes leads to a predictor that has superior prediction accuracy over any single approach. Finally, appropriate weighting of samples in the training process may achieve higher accuracy for samples with low fetal fraction and so allow more reliability for subsequent testing for genomic aberrations. We propose several improvements in fetal fraction estimation with a special focus on the samples most prone to wrong conclusion.

Novel SCN1A variants in Dravet syndrome and evaluating a wide approach of patient selection.

Voltage-gated sodium channels are essential for generation and propagation of the action potential mainly in nerve and muscle cells. Causative variants in SCN1A gene which codes the main, pore-forming subunit of the channel expressed in central nervous system are associated predominantly with Dravet syndrome (DS), as well as with generalized epilepsy with febrile seizures plus (GEFS+) making it one of the most significant epilepsy gene. Our goal was to determine whether SCN1A screening is relevant in patients with a broad range of epileptic syndromes. 52 patients diagnosed with DS, generalized epilepsy with GEFS+ or similar types of epileptic syndromes were included. Sequencing of the protein coding parts of the gene complemented with MLPA analysis was carried out. One already described nonsense variant, four novel protein truncating variants and a deletion encompassing the whole SCN1A gene were revealed, all in heterozygous state. All identified variants were found in DS patients with 85.7% sensitivity, thus supporting the role of profound SCN1A gene variants in etiology of DS phenotype. No causative variants were identified in any of non-DS epileptic patients in our cohort, suggesting a minor, but not irrelevant role for SCN1A in patients with other types of childhood epilepsy.

In silico and experimental evaluation of DNA-based detection methods for the ability to discriminate almond from other Prunus spp.

Ten published DNA-based analytical methods aiming at detecting material of almond (Prunus dulcis) were in silico evaluated for potential cross-reactivity with other stone fruits (Prunus spp.), including peach, apricot, plum, cherry, sour cherry and Sargent cherry. For most assays, the analysis of nucleotide databases suggested none or insufficient discrimination of at least some stone fruits. On the other hand, the assay targeting non-specific lipid transfer protein (Röder et al., 2011, Anal Chim Acta 685:74-83) was sufficiently discriminative, judging from nucleotide alignments. Empirical evaluation was performed for three of the published methods, one modification of a commercial kit (SureFood allergen almond) and one attempted novel method targeting thaumatin-like protein gene. Samples of leaves and kernels were used in the experiments. The empirical results were favourable for the method from Röder et al. (2011) and a modification of SureFood allergen almond kit, both showing cross-reactivity <10(-3) compared to the model almond.

The origin of the giant liver fluke, Fascioloides magna (Trematoda: Fasciolidae) from Croatia determined by high-resolution melting screening of mitochondrial cox1 haplotypes.

The high-resolution melting (HRM) method, recently optimized as a reliable technique for population study of the European Fascioloides magna populations, was applied to determine an origin of F. magna individuals from Croatia. The structure and frequency of mitochondrial cytochrome c oxidase subunit I (439 bp; cox1) haplotypes of 200 Croatian flukes coming from 19 red deer (Cervus elaphus elaphus) livers were screened and compared with recently determined reference samples of F. magna from all European foci-Italy, Czech Republic, and Danube floodplain forests. While the reference haplotypes Ha1 and Ha2 were specific for flukes from the first European focus of fascioloidosis, the Natural Park La Mandria in Italy, the remaining three haplotypes (Ha3, Ha4, and Ha5) represented parasites from the second focus, Czech Republic. Besides, Ha3 and Ha4 were found also in the third, latest, and still expanding European focus, the Danube floodplain forests. The HRM screening of cox1 haplotypes of Croatian F. magna individuals resulted in classification of samples into the two mitochondrial haplogroups characterized by well-distinguished melting curves. They corresponded to Ha3 and Ha4 reference haplotypes that confirmed the Danube origin of F. magna from Croatia. The results support the theory that the Danube floodplain forests population of F. magna represents uniform genetic pool of the parasite. The spread of F. magna alongside the Danube River down to Croatia was possible due to suitable ecological conditions for definitive and intermediate hosts present in this unique biotope.

Microsatellite instability assessment is instrumental for Predictive, Preventive and Personalised Medicine: status quo and outlook.

A form of genomic alteration called microsatellite instability (MSI) occurs in a class of tandem repeats (TRs) called microsatellites (MSs) or short tandem repeats (STRs) due to the failure of a post-replicative DNA mismatch repair (MMR) system. Traditionally, the strategies for determining MSI events have been low-throughput procedures that typically require assessment of tumours as well as healthy samples. On the other hand, recent large-scale pan-tumour studies have consistently highlighted the potential of massively parallel sequencing (MPS) on the MSI scale. As a result of recent innovations, minimally invasive methods show a high potential to be integrated into the clinical routine and delivery of adapted medical care to all patients. Along with advances in sequencing technologies and their ever-increasing cost-effectiveness, they may bring about a new era of Predictive, Preventive and Personalised Medicine (3PM). In this paper, we offered a comprehensive analysis of high-throughput strategies and computational tools for the calling and assessment of MSI events, including whole-genome, whole-exome and targeted sequencing approaches. We also discussed in detail the detection of MSI status by current MPS blood-based methods and we hypothesised how they may contribute to the shift from conventional medicine to predictive diagnosis, targeted prevention and personalised medical services. Increasing the efficacy of patient stratification based on MSI status is crucial for tailored decision-making. Contextually, this paper highlights drawbacks both at the technical level and those embedded deeper in cellular/molecular processes and future applications in routine clinical testing.

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/ .

SnakeLines: integrated set of computational pipelines for sequencing reads.

With the rapid growth of massively parallel sequencing technologies, still more laboratories are utilising sequenced DNA fragments for genomic analyses. Interpretation of sequencing data is, however, strongly dependent on bioinformatics processing, which is often too demanding for clinicians and researchers without a computational background. Another problem represents the reproducibility of computational analyses across separated computational centres with inconsistent versions of installed libraries and bioinformatics tools. We propose an easily extensible set of computational pipelines, called SnakeLines, for processing sequencing reads; including mapping, assembly, variant calling, viral identification, transcriptomics, and metagenomics analysis. Individual steps of an analysis, along with methods and their parameters can be readily modified in a single configuration file. Provided pipelines are embedded in virtual environments that ensure isolation of required resources from the host operating system, rapid deployment, and reproducibility of analysis across different Unix-based platforms. SnakeLines is a powerful framework for the automation of bioinformatics analyses, with emphasis on a simple set-up, modifications, extensibility, and reproducibility. The framework is already routinely used in various research projects and their applications, especially in the Slovak national surveillance of SARS-CoV-2.

Repurposing non-invasive prenatal testing data: Population study of single nucleotide variants associated with colorectal cancer and Lynch syndrome.

In our previous work, genomic data generated through non-invasive prenatal testing (NIPT) based on low-coverage massively parallel whole-genome sequencing of total plasma DNA of pregnant women in Slovakia was described as a valuable source of population specific data. In the present study, these data were used to determine the population allele frequency of common risk variants located in genes associated with colorectal cancer (CRC) and Lynch syndrome (LS). Allele frequencies of identified variants were compared with six world populations to detect significant differences between populations. Finally, variants were interpreted, functional consequences were searched for and clinical significance of variants was investigated using publicly available databases. Although the present study did not identify any pathogenic variants associated with CRC or LS in the Slovak population using NIPT data, significant differences were observed in the allelic frequency of risk CRC variants previously reported in genome-wide association studies and common variants located in genes associated with LS. As Slovakia is one of the leading countries with the highest incidence of CRC among male patients in the world, there is a need for studies dedicated to investigating the cause of such a high incidence of CRC in Slovakia. The present study also assumed that extensive cross-country data aggregation of NIPT results would represent an unprecedented source of information concerning human genome variation in cancer research.

DNA copy number variation: Main characteristics, evolutionary significance, and pathological aspects.

Copy number variants (CNVs) were the subject of extensive research in the past years. They are common features of the human genome that play an important role in evolution, contribute to population diversity, development of certain diseases, and influence host-microbiome interactions. CNVs have found application in the molecular diagnosis of many diseases and in non-invasive prenatal care, but their full potential is only emerging. CNVs are expected to have a tremendous impact on screening, diagnosis, prognosis, and monitoring of several disorders, including cancer and cardiovascular disease. Here, we comprehensively review basic definitions of the term CNV, outline mechanisms and factors involved in CNV formation, and discuss their evolutionary and pathological aspects. We suggest a need for better defined distinguishing criteria and boundaries between known types of CNVs.