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.

Determination of the Prevalence of Microsatellite Instability, BRAF and KRAS/NRAS Mutation Status in Patients with Colorectal Cancer in Slovakia.

Slovakia has one of the highest rates of colorectal cancer among the developed countries, ranking as the second highest in the incidence of this disease for men worldwide. Despite the significant burden on both quality of life and the healthcare system this disease imposes, data on molecular analysis of biomarkers in CRC-diagnosed patients is scarce. In our study, we analyzed confirmed CRC patients from the database of the National Cancer Institute (NCI) and evaluated the presence of 4 biomarkers in tumor tissues. Altogether, 83 FFPE tumor tissues from CRC patients listed in the NCI database were analyzed for microsatellite instability status, presence of BRAF and KRAS/NRAS mutations, and neoplastic cell percentage in tissue samples. We identified 4 MSI-high samples, 39 KRAS/NRAS mutations, and 5 BRAF p.V600E mutations, with one case of coexistence of all three markers in a single tumor sample. We also evaluated possible relationships between biomarkers, their coexistence, and the age and sex of the studied population.

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.

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.

Carbapenem-Resistant Klebsiella pneumoniae in COVID-19 Era-Challenges and Solutions.

The COVID-19 era brought about new medical challenges, which, together with nosocomial bacterial infections, resulted in an enormous burden for the healthcare system. One of the most alarming nosocomial threats was carbapenem-resistant Klebsiella pneumoniae (CRKP). Monitoring CRKP incidence and antimicrobial resistance globally and locally is vitally important. In a retrospective study, the incidence of CRKP in the pre-COVID-19 period (2017-2019) and the COVID-19 pandemic (2020-2022) was investigated in the Central Military Hospital in Ruzomberok, Slovak Republic. The relative incidence of CRKP significantly increased during the COVID-19 period-by 4.8 times, from 0.18 to 0.76%. At the same time, 47% of CRKP-positive patients also had COVID-19. Twenty-six KPC and sixty-nine NDM-producing isolates were identified. CRKPs isolated in the year 2022 were submitted to whole genome sequencing, and their susceptibility was tested to cefiderocol, ceftazidime-avibactam, imipenem-relebactam and meropenem-vaborbactam, with excellent results. KPC-producing isolates were also highly susceptible to colistin (92%). The NDM isolates revealed lower susceptibility rates, including only 57% colistin susceptibility. ST-307 prevailed in KPC and ST-11 in NDM isolates. Despite the excellent activity of new antimicrobials, rational antibiotic policy must be thoroughly followed, supported by complementary treatments and strict anti-epidemic precautions.

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.

New Possible Ways to Use Exosomes in Diagnostics and Therapy via JAK/STAT Pathways.

Exosomes have the potential to be the future of personalized diagnostics and therapy. They are nano-sized particles between 30 and 100 nm flowing in the extracellular milieu, where they mediate cell-cell communication and participate in immune system regulation. Tumor-derived exosomes (TDEs) secreted from different types of cancer cells are the key regulators of the tumor microenvironment. With their immune suppressive cargo, TDEs prevent the antitumor immune response, leading to reduced effectiveness of cancer treatment by promoting a pro-tumorigenic microenvironment. Involved signaling pathways take part in the regulation of tumor proliferation, differentiation, apoptosis, and angiogenesis. Signal transducers and activators of transcription factors (STATs) and Janus kinase (JAK) signaling pathways are crucial in malignancies and autoimmune diseases alike, and their potential to be manipulated is currently the focus of interest. In this review, we aim to discuss exosomes, TDEs, and the JAK/STAT pathways, along with mediators like interleukins, tripartite motif proteins, and interferons.

Cross-Kingdom Interaction of miRNAs and Gut Microbiota with Non-Invasive Diagnostic and Therapeutic Implications in Colorectal Cancer.

Colorectal cancer (CRC) has one of the highest incidences among all types of malignant diseases, affecting millions of people worldwide. It shows slow progression, making it preventable. However, this is not the case due to shortcomings in its diagnostic and management procedure and a lack of effective non-invasive biomarkers for screening. Here, we discuss CRC-associated microRNAs (miRNAs) and gut microbial species with potential as CRC diagnostic and therapy biomarkers. We provide rich evidence of cross-kingdom miRNA-mediated interactions between the host and gut microbiome. miRNAs have emerged with the ability to shape the composition and dynamics of gut microbiota. Intestinal microbes can uptake miRNAs, which in turn influence microbial growth and provide the ability to regulate the abundance of various microbial species. In the context of CRC, targeting miRNAs could aid in manipulating the balance of the microbiota. Our findings suggest the need for correlation analysis between the composition of the gut microbiome and the miRNA expression profile.

PHERI-Phage Host ExploRation Pipeline.

Antibiotic resistance is becoming a common problem in medicine, food, and industry, with multidrug-resistant bacterial strains occurring in all regions. One of the possible future solutions is the use of bacteriophages. Phages are the most abundant form of life in the biosphere, so we can highly likely purify a specific phage against each target bacterium. The identification and consistent characterization of individual phages was a common form of phage work and included determining bacteriophages' host-specificity. With the advent of new modern sequencing methods, there was a problem with the detailed characterization of phages in the environment identified by metagenome analysis. The solution to this problem may be to use a bioinformatic approach in the form of prediction software capable of determining a bacterial host based on the phage whole-genome sequence. The result of our research is the machine learning algorithm-based tool called PHERI. PHERI predicts the suitable bacterial host genus for the purification of individual viruses from different samples. In addition, it can identify and highlight protein sequences that are important for host selection.

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

Complete Genome Sequence of Cronobacter dublinensis Bacteriophage vB_Cdu_VP8.

Cronobacter dublinensis is a Gram-negative pathogen that is capable of causing infection in humans. In this announcement, we describe the characterization of bacteriophage vB_Cdu_VP8, which is able to lyse a Cronobacter dublinensis strain. Related to phages belonging to the genus Muldoonvirus, such as Muldoon and SP1, vB_Cdu_VP8 contains 264 predicted protein-coding genes and 3 tRNAs.

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

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.

Complete Genome Sequence of New Cronobacter-Specific Bacteriophage Dev_CS701.

Cronobacter is a ubiquitous Gram-negative bacterium linked with serious infections. In this report, we describe the characterization of Cronobacter phage Dev_CS701, which was isolated from wastewater. Related to phages belonging to the Straboviridae family and Pseudotevenvirus genus, such as vB_CsaM_IeB, Dev_CS701 contains 257 predicted protein-coding genes and a tRNA gene.

Telomere Length Changes in Cancer: Insights on Carcinogenesis and Potential for Non-Invasive Diagnostic Strategies.

Telomere dynamics play a crucial role in the maintenance of chromosome integrity; changes in telomere length may thus contribute to the development of various diseases including cancer. Understanding the role of telomeric DNA in carcinogenesis and detecting the presence of cell-free telomeric DNA (cf-telDNA) in body fluids offer a potential biomarker for novel cancer screening and diagnostic strategies. Liquid biopsy is becoming increasingly popular due to its undeniable benefits over conventional invasive methods. However, the organization and function of cf-telDNA in the extracellular milieu are understudied. This paper provides a review based on 3,398,017 cancer patients, patients with other conditions, and control individuals with the aim to shed more light on the inconsistent nature of telomere lengthening/shortening in oncological contexts. To gain a better understanding of biological factors (e.g., telomerase activation, alternative lengthening of telomeres) affecting telomere homeostasis across different types of cancer, we summarize mechanisms responsible for telomere length maintenance. In conclusion, we compare tissue- and liquid biopsy-based approaches in cancer assessment and provide a brief outlook on the methodology used for telomere length evaluation, highlighting the advances of state-of-the-art approaches in the field.

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.

Extracellular vesicles derived from dental mesenchymal stem/stromal cells with gemcitabine as a cargo have an inhibitory effect on the growth of pancreatic carcinoma cell lines in vitro.

Extracellular vesicles (EVs) are nowadays a target of interest in cancer therapy as a successful drug delivering tool. Based on their many beneficial biocompatible properties are designed to transport nucleic acids, proteins, various nanomaterials or chemotherapeutics. Extracellular vesicles derived from mesenchymal stem/stromal cells (MSCs) possess their tumor-homing abilities. This inspired us to engineer the MSC's EVs to be packed with chemotherapeutic agents and deliver it as a Trojan horse directly into tumor cells. In our study, human dental pulp MSCs (DP-MSCs) were cultivated with gemcitabine (GCB), which led to its absorption by the cells and subsequent secretion of the drug out into conditioned media in EVs. Concentrated conditioned media containing small EVs (potentially exosomes) significantly inhibited the cell growth of pancreatic carcinoma cell lines in vitro. DP-MSCs were simultaneously engineered to express a suicide gene fused yeast cytosinedeaminase:uracilphosphoribosyltransferase (yCD::UPRT). The product of the suicide gene converts non-toxic prodrug 5-fluorocytosine (5-FC) to highly cytotoxic chemotherapeutic drug 5-fluorouracil (5-FU) in the recipient cancer cells. Conversion of 5-FC to 5-FU had an additional effect on cancer cell's growth inhibition. Our results showed a therapeutic potential for DP-MSC-EVs to be designed for successful delivering of chemotherapeutic drugs, together with prodrug suicide gene therapy system.

Mitochondrial DNA copy number changes, heteroplasmy, and mutations in plasma-derived exosomes and brain tissue of glioblastoma patients.

Glioblastoma is the most common malignant tumor of the central nervous system (CNS) in adults. Glioblastoma cells show increased glucose consumption associated with poor prognosis. Since mitochondria play a crucial role in energy metabolism, mutations and copy number changes of mitochondrial DNA may serve as biomarkers. As the brain is difficult to access, analysis of mitochondria directly from the brain tissue represents a challenge. Exosome analysis is an alternative (still poorly explored) approach to investigate molecular changes in CNS tumors. We analyzed brain tissue DNA and plasma-derived exosomal DNA (exoDNA) of 44 glioblastoma patients and 40 control individuals. Quantitative real-time PCR was performed to determine mtDNA copy numbers and the Kruskal-Wallis and Mann-Whitney U test were used for statistical analysis of data. Subsequently, sequencing libraries were prepared and sequenced on the MiSeq platform to identify mtDNA point mutations. Tissue mtDNA copy number was different among controls and patients in multiple comparisons. A similar tendency was detected in exosomes. Based on NGS analysis, several mtDNA point mutations showed slightly different frequencies between cases and controls, but the clinical relevance of these observations is difficult to assess and likely less than that of overall mtDNA copy number changes. Allele frequencies of variants were used to determine the level of heteroplasmy (found to be higher in exo-mtDNA of control individuals). Despite the suggested potential, the use of such biomarkers for the screening and/or diagnosis of glioblastomas is still limited, thus further studies are needed.

Systematic Genomic Surveillance of SARS-CoV-2 Virus on Illumina Sequencing Platforms in the Slovak Republic-One Year Experience.

To explore a genomic pool of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during the pandemic, the Ministry of Health of the Slovak Republic formed a genomics surveillance workgroup, and the Public Health Authority of the Slovak Republic launched a systematic national epidemiological surveillance using whole-genome sequencing (WGS). Six out of seven genomic centers implementing Illumina sequencing technology were involved in the national SARS-CoV-2 virus sequencing program. Here we analyze a total of 33,024 SARS-CoV-2 isolates collected from the Slovak population from 1 March 2021, to 31 March 2022, that were sequenced and analyzed in a consistent manner. Overall, 28,005 out of 30,793 successfully sequenced samples met the criteria to be deposited in the global GISAID database. During this period, we identified four variants of concern (VOC)-Alpha (B.1.1.7), Beta (B.1.351), Delta (B.1.617.2) and Omicron (B.1.1.529). In detail, we observed 165 lineages in our dataset, with dominating Alpha, Delta and Omicron in three major consecutive incidence waves. This study aims to describe the results of a routine but high-level SARS-CoV-2 genomic surveillance program. Our study of SARS-CoV-2 genomes in collaboration with the Public Health Authority of the Slovak Republic also helped to inform the public about the epidemiological situation during the pandemic.

Non-intuitive trends of fetal fraction development related to gestational age and fetal gender, and their practical implications for non-invasive prenatal testing.

Discovery of fetal cell-free DNA fragments in maternal blood revolutionized prenatal diagnostics. Although non-invasive prenatal testing (NIPT) is already a matured screening test with high specificity and sensitivity, the accurate estimation of the proportion of fetal fragments, called fetal fraction, is crucial to avoid false-negative results. In this study, we collected 6999 samples from women undergoing NIPT testing with a single male fetus to demonstrate the influence of fetal fraction by the maternal and fetal characteristics. We show several fetal fraction discrepancies that contradict the generally presented conventional view. At first, the fetal fraction is not consistently rising with the maturity of the fetus due to a drop in 15 weeks of maturation. Secondly, the male samples have a lower fetal fraction than female fetuses, arguably due to the smaller gonosomal chromosomes. Finally, we discuss not only the possible reasons why this inconsistency exists but we also outline why these differences have not yet been identified and published. We demonstrate two non-intuitive trends to better comprehend the fetal fraction development and more precise selection of patients with sufficient fetal fraction for accurate testing.