SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants.

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

An Explainable Deep Learning Classifier of Bovine Mastitis Based on Whole-Genome Sequence Data-Circumventing the p >> n Problem.

The serious drawback underlying the biological annotation of whole-genome sequence data is the p >> n problem, which means that the number of polymorphic variants (p) is much larger than the number of available phenotypic records (n). We propose a way to circumvent the problem by combining a LASSO logistic regression with deep learning to classify cows as susceptible or resistant to mastitis, based on single nucleotide polymorphism (SNP) genotypes. Among several architectures, the one with 204,642 SNPs was selected as the best. This architecture was composed of two layers with, respectively, 7 and 46 units per layer implementing respective drop-out rates of 0.210 and 0.358. The classification of the test data resulted in AUC = 0.750, accuracy = 0.650, sensitivity = 0.600, and specificity = 0.700. Significant SNPs were selected based on the SHapley Additive exPlanation (SHAP). As a final result, one GO term related to the biological process and thirteen GO terms related to molecular function were significantly enriched in the gene set that corresponded to the significant SNPs. Our findings revealed that the optimal approach can correctly predict susceptibility or resistance status for approximately 65% of cows. Genes marked by the most significant SNPs are related to the immune response and protein synthesis.

Whole-genome sequencing of Ganoderma boninense, the causal agent of basal stem rot disease in oil palm, via combined short- and long-read sequencing.

The hemibiotrophic Basidiomycete pathogen Ganoderma boninense (Gb) is the dominant causal agent of oil palm basal stem rot disease. Here, we report a complete chromosomal genome map of Gb using a combination of short-read Illumina and long-read Pacific Biosciences (PacBio) sequencing platforms combined with chromatin conformation capture data from the Chicago and Hi-C platforms. The genome was 55.87 Mb in length and assembled to a high contiguity (N50: 304.34 kb) of 12 chromosomes built from 112 scaffolds, with a total of only 4.34 Mb (~ 7.77%) remaining unplaced. The final assemblies were evaluated for completeness of the genome by using Benchmarking Universal Single Copy Orthologs (BUSCO) v4.1.4, and based on 4464 total BUSCO polyporales group searches, the assemblies yielded 4264 (95.52%) of the conserved orthologs as complete and only a few fragmented BUSCO of 42 (0.94%) as well as a missing BUSCO of 158 (3.53%). Genome annotation predicted a total of 21,074 coding genes, with a GC content ratio of 59.2%. The genome features were analyzed with different databases, which revealed 2471 Gene Ontology/GO (11.72%), 5418 KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthologous/KO (25.71%), 13,913 Cluster of Orthologous Groups of proteins/COG (66.02%), 60 ABC transporter (0.28%), 1049 Carbohydrate-Active Enzymes/CAZy (4.98%), 4005 pathogen-host interactions/PHI (19%), and 515 fungal transcription factor/FTFD (2.44%) genes. The results obtained in this study provide deep insight for further studies in the future.

First detection and characterization of mcr-1 colistin resistant E. coli from wild rat in Bangladesh.

Colistin resistance is a global concern warning for a one health approach to combat the challenge. Colistin resistant E. coli and their resistance determinants are widely distributed in the environment, and rats could be a potential source of these isolates and resistant determinants to a diverse environmental setting. This study was aimed to determine the presence of colistin resistant E. coli (CREC) in wild rats, their antimicrobial resistance (AMR) phenotypes, and genotypic analysis of mcr-1 CREC through whole genome sequencing (WGS). A total of 39 rats were examined and CREC was isolated from their fecal pellets onto MacConkey agar containing colistin sulfate (1 µg/ mL). AMR of the CREC was determined by disc diffusion and broth microdilution was employed to determine MIC to colistin sulfate. CREC were screened for mcr genes (mcr-1 to mcr-8) and phylogenetic grouping by PCR. Finally, WGS of one mcr-1 CREC was performed to explore its genetic characteristics especially resistomes and virulence determinants. 43.59% of the rats carried CREC with one (2.56%) of them carrying CREC with mcr-1 gene among the mcr genes examined. Examination of seventeen (17) isolates from the CREC positive rats (n = 17) revealed that majority of them belonging to the pathogenic phylogroup D (52.94%) and B2 (11.76%). 58.82% of the CREC were MDR on disc diffusion test. Shockingly, the mcr-1 CREC showed phenotypic resistance to 16 antimicrobials of 8 different classes and carried the ARGs in its genome. The mcr-1 gene was located on a 60 kb IncI2 plasmid. On the other hand, ARGs related to aminoglycosides, phenicols, sulfonamides, tetracyclines and trimethoprims were located on a 288 kb mega-plasmid separately. The mcr-1 CREC carried 58 virulence genes including genes related to adhesion, colonization, biofilm formation, hemolysis and immune-evasion. The isolate belonged to ST224 and closely related to E. coli from different sources including UPEC clinical isolates from human based on cgMLST analysis. The current research indicates that rats might be a possible source of CREC, and the presence of mcr-1 and other ARGs on plasmid increases the risk of ARGs spreading and endangering human health and other environmental components through this infamous pest.

Recurrent neural network for predicting absence of heterozygosity from low pass WGS with ultra-low depth.

BACKGROUND: The absence of heterozygosity (AOH) is a kind of genomic change characterized by a long contiguous region of homozygous alleles in a chromosome, which may cause human genetic disorders. However, no method of low-pass whole genome sequencing (LP-WGS) has been reported for the detection of AOH in a low-pass setting of less than onefold. We developed a method, termed CNVseq-AOH, for predicting the absence of heterozygosity using LP-WGS with ultra-low sequencing data, which overcomes the sparse nature of typical LP-WGS data by combing population-based haplotype information, adjustable sliding windows, and recurrent neural network (RNN). We tested the feasibility of CNVseq-AOH for the detection of AOH in 409 cases (11 AOH regions for model training and 863 AOH regions for validation) from the 1000 Genomes Project (1KGP). AOH detection using CNVseq-AOH was also performed on 6 clinical cases with previously ascertained AOHs by whole exome sequencing (WES). RESULTS: Using SNP-based microarray results as reference (AOHs detected by CNVseq-AOH with at least a 50% overlap with the AOHs detected by chromosomal microarray analysis), 409 samples (863 AOH regions) in the 1KGP were used for concordant analysis. For 784 AOHs on autosomes and 79 AOHs on the X chromosome, CNVseq-AOH can predict AOHs with a concordant rate of 96.23% and 59.49% respectively based on the analysis of 0.1-fold LP-WGS data, which is far lower than the current standard in the field. Using 0.1-fold LP-WGS data, CNVseq-AOH revealed 5 additional AOHs (larger than 10 Mb in size) in the 409 samples. We further analyzed AOHs larger than 10 Mb, which is recommended for reporting the possibility of UPD. For the 291 AOH regions larger than 10 Mb, CNVseq-AOH can predict AOHs with a concordant rate of 99.66% with only 0.1-fold LP-WGS data. In the 6 clinical cases, CNVseq-AOH revealed all 15 known AOH regions. CONCLUSIONS: Here we reported a method for analyzing LP-WGS data to accurately identify regions of AOH, which possesses great potential to improve genetic testing of AOH.

Evaluation of the accuracy of bacterial genome reconstruction with Oxford Nanopore R10.4.1 long-read-only sequencing.

Whole-genome reconstruction of bacterial pathogens has become an important tool for tracking transmission and antimicrobial resistance gene spread, but highly accurate and complete assemblies have largely only historically been achievable using hybrid long- and short-read sequencing. We previously found the Oxford Nanopore Technologies (ONT) R10.4/kit12 flowcell/chemistry produced improved assemblies over the R9.4.1/kit10 combination, however long-read only assemblies contained more errors compared to Illumina-ONT hybrid assemblies. ONT have since released an R10.4.1/kit14 flowcell/chemistry upgrade and recommended the use of Bovine Serum Albumin (BSA) during library preparation, both of which reportedly increase accuracy and yield. They have also released updated basecallers trained using native bacterial DNA containing methylation sites intended to fix systematic basecalling errors, including common adenosine (A) to guanine (G) and cytosine (C) to thymine (T) substitutions. To evaluate these improvements, we successfully sequenced four bacterial reference strains, namely Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, and nine genetically diverse E. coli bloodstream infection-associated isolates from different phylogroups and sequence types, both with and without BSA. These sequences were de novo assembled and compared against Illumina-corrected reference genomes. In this small evaluation of 13 isolates we found that nanopore long-read-only R10.4.1/kit 14 assemblies with updated basecallers trained using bacterial methylated DNA produce accurate assemblies with ≥40×depth, sufficient to be cost-effective compared with hybrid ONT/Illumina sequencing in our setting.

Rapid detection of the source of a Listeria monocytogenes outbreak in Switzerland through routine interviewing of patients and whole-genome sequencing.

AIMS OF THE STUDY: Listeriosis is a notifiable disease in Switzerland. In summer 2022, the Swiss Federal Office of Public Health noticed an increase in reports of listeriosis cases, indicating a possible ongoing outbreak. Here we present the approaches applied for rapidly confirming the outbreak, detecting the underlying source of infection and the measures put in place to eliminate it and contain the outbreak. METHODS: For close surveillance and early detection of outbreak situations with their possible sources, listeriosis patients in Switzerland are systematically interviewed about risk behaviours and foods consumed prior to the infection. Listeria monocytogenes isolates derived from patients in medical laboratories are sent to the National Reference Laboratory for Enteropathogenic Bacteria and Listeria, where they routinely undergo whole-genome sequencing. Interview and whole-genome sequencing data are continuously linked for comparison and analysis. RESULTS: In summer 2022, 20 patient-derived L. monocytogenes serotype 4b sequence type 388 strains were found to belong to an outbreak cluster (≤10 different alleles between neighbouring isolates) based on core genome multilocus sequence typing analysis. Geographically, 18 of 20 outbreak cases occurred in northeastern Switzerland. The median age of patients was 77.4 years (range: 58.1-89.7), with both sexes equally affected. Rolling analysis of the interview data revealed smoked trout from a local producer as a suspected infection source, triggering an on-site investigation of the production facility and sampling of the suspected products by the responsible cantonal food inspection team on 15 July 2022. Seven of ten samples tested positive for L. monocytogenes and the respective cantonal authority ordered a ban on production and distribution as well as a product recall. The Federal Food Safety and Veterinary Office released a nationwide public alert covering the smoked fish products concerned. Whole-genome sequencing analysis confirmed the interrelatedness of the L. monocytogenes smoked trout product isolates and the patient-derived isolates. Following the ban on production and distribution and the product recall, reporting of new outbreak-related cases rapidly dropped to zero. CONCLUSIONS: This listeriosis outbreak could be contained within a relatively short time thanks to identification of the source of contamination through the established combined approach of timely interviewing of every listeriosis patient or a representative and continuous molecular analysis of the patient- and food-derived L. monocytogenes isolates. These findings highlight the effectiveness of this well-established, joint approach involving the federal and cantonal authorities and the research institutions mandated to contain listeriosis outbreaks in Switzerland.

Multimodal analysis of cfDNA methylomes for early detecting esophageal squamous cell carcinoma and precancerous lesions.

Detecting early-stage esophageal squamous cell carcinoma (ESCC) and precancerous lesions is critical for improving survival. Here, we conduct whole-genome bisulfite sequencing (WGBS) on 460 cfDNA samples from patients with non-metastatic ESCC or precancerous lesions and matched healthy controls. We develop an expanded multimodal analysis (EMMA) framework to simultaneously identify cfDNA methylation, copy number variants (CNVs), and fragmentation markers in cfDNA WGBS data. cfDNA methylation markers are the earliest and most sensitive, detectable in 70% of ESCCs and 50% of precancerous lesions, and associated with molecular subtypes and tumor microenvironments. CNVs and fragmentation features show high specificity but are linked to late-stage disease. EMMA significantly improves detection rates, increasing AUCs from 0.90 to 0.99, and detects 87% of ESCCs and 62% of precancerous lesions with >95% specificity in validation cohorts. Our findings demonstrate the potential of multimodal analysis of cfDNA methylome for early detection and monitoring of molecular characteristics in ESCC.

The complete genome sequence of unculturable Mycoplasma faucium obtained through clinical metagenomic next-generation sequencing.

Introduction: Diagnosing Mycoplasma faucium poses challenges, and it's unclear if its rare isolation is due to infrequent occurrence or its fastidious nutritional requirements. Methods: This study analyzes the complete genome sequence of M. faucium, obtained directly from the pus of a sternum infection in a lung transplant patient using metagenomic sequencing. Results: Genome analysis revealed limited therapeutic options for the M. faucium infection, primarily susceptibility to tetracyclines. Three classes of mobile genetic elements were identified: two new insertion sequences, a new prophage (phiUMCG-1), and a species-specific variant of a mycoplasma integrative and conjugative element (MICE). Additionally, a Type I Restriction-Modification system was identified, featuring 5'-terminally truncated hsdS pseudogenes with overlapping repeats, indicating the potential for forming alternative hsdS variants through recombination. Conclusion: This study represents the first-ever acquisition of a complete circularized bacterial genome directly from a patient sample obtained from invasive infection of a primary sterile site using culture-independent, PCR-free clinical metagenomics.

Whole-genome sequences restore the original classification of dabbling ducks (genus Anas).

Anas, is a genus of dabbling ducks and encompasses a considerable number of species, among which some are the progenitors of domestic ducks. However, the taxonomic position of the Anas genus remains uncertain because several of its species, initially categorized as Anas based on morphological characteristics, were subsequently reclassified and grouped with the South American genus Tachyeres, primarily based on analysis of their mitochondrial gene sequences. Here, we constructed a phylogenetic tree using nine of our recently assembled Anas genomes, two Tachyeres genomes, and one Cairina genome that are publicly available. The results showed that the Northern shoveler (Anas clypeata) and Baikal teal (Anas formosa) clustered with the other Anas species at the whole-genome level rather than with the Steamer ducks (genus Tachyeres). Therefore, we propose to restore the original classification of the Anas genus, which includes the Northern shoveler and Baikal teal species, 47 species in total. Moreover, our study unveiled extensive incomplete lineage sorting and an ancient introgression event from Tachyeres to Anas, which has led to notable phylogenetic incongruence within the Anas genome. This ancient introgression event not only supports the theory that Anas originated in South America but also that it played a significant role in shaping the evolutionary trajectory of Anas, including the domestic duck.

DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus.

BACKGROUND: The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes. RESULTS: We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 10-100 RNA copies/µL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments. CONCLUSIONS: DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

[Carbapenemase Genes, Virulence Genes, and Molecular Epidemiology of Carbapenem-Resistant Klebsiella pneumoniae Derived From Bloodstream Infections]. / è¡€æµæ„ŸæŸ“ç¢³é’éœ‰çƒ¯è€è¯è‚ºç‚Žå ‹é›·ä¼¯èŒçš„è€è¯åŸºå› å’Œæ¯’åŠ›åŸºå› åŠåˆ†å­æµè¡Œç— å­¦ç ”ç©¶.

Objective: To investigate the clinical characteristics and molecular epidemiology of carbapenem-resistant Klebsiella pneumoniae (CRKP) isolated from patients with bloodstream infections in a large tertiary-care general hospital in Southwest China. Methods: A total of 131 strains of non-repeating CRKP were collected from the blood cultures of patients who had bloodstream infections in 2015-2019. The strains were identified by VITEK-2, a fully automated microbial analyzer, and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry. The minimum inhibitory concentration (MIC) was determined by microbroth dilution method. The common carbapenemase resistant genes and virulence factors were identified by PCR. Homology analysis was performed by multilocus sequencing typing. Whole genome sequencing was performed to analyze the genomic characteristics of CRKP without carbapenemase. Results: The 131 strains of CRKP showed resistance to common antibiotics, except for polymyxin B (1.6% resistance rate) and tigacycline (8.0% resistance rate). A total of 105 (80.2%) CRKP strains carried the Klebsiella pneumoniae carbapenemase (KPC) resistance gene, 15 (11.4%) strains carried the New Delhi Metallo-ß-lactamase (NDM) gene, and 4 (3.1%) isolates carried both KPC and NDM genes. Sequence typing (ST) 11 (74.0%) was the dominant sequence type. High detection rates for mrkD (96.2%), fimH (98.5%), entB (100%), and other virulence genes were reported. One hypervirulent CRKP strain was detected. The seven strains of CRKP that did not produce carbapenemase were shown to carry ESBL or AmpC genes and had anomalies in membrane porins OMPK35 and OMPK36, according to whole genome sequencing. Conclusion: In a large-scale tertiary-care general hospital, CRKP mainly carries the KPC gene, has a high drug resistance rate to a variety of antibiotics, and possesses multiple virulence genes. Attention should be paid to CRKP strains with high virulence.

Identification of Novel Induced Mutations in Seed and Vegetatively Propagated Plants from Reduced Representation or Whole Genome Sequencing Data.

Treatment of plants with chemical mutagens results primarily in the production of novel single nucleotide variants. Mutagenesis is a mostly random process and as such plants derived from mutagenesis of different seeds or in vitro material are expected to accumulate different mutations. An important step in the creation of a mutant population for forward or reverse genetics is the choice of treatment conditions (e.g., dosage) such that sufficient mutations accumulate while not adversely affecting propagation of the plant. DNA sequencing provides a quick method to evaluate the effect of different treatment conditions and their effect on the density and spectrum of accumulated mutations. Whole genome sequencing or reduced representation sequencing is carried out followed by mapping to a reference genome and production of a Variant Call Format (VCF) file. We provide here a method for generating a multi-sample VCF from mutagenized plants and describe a new tool to streamline the process of recovering unique induced mutations and determining their possible effect on gene function.

Exploring the tumor genomic landscape of aggressive prostate cancer by whole-genome sequencing of tissue or liquid biopsies.

Treatment resistance remains a major issue in aggressive prostate cancer (PC), and novel genomic biomarkers may guide better treatment selection. Circulating tumor DNA (ctDNA) can provide minimally invasive information about tumor genomes, but the genomic landscape of aggressive PC based on whole-genome sequencing (WGS) of ctDNA remains incompletely characterized. Thus, we here performed WGS of tumor tissue (n = 31) or plasma ctDNA (n = 10) from a total of 41 aggressive PC patients, including 11 hormone-naïve, 15 hormone-sensitive, and 15 castration-resistant patients. Across all variant types, we found progressively more altered tumor genomic profiles in later stages of aggressive PC. The potential driver genes most frequently affected by single-nucleotide variants or insertions/deletions included the known PC-related genes TP53, CDK12, and PTEN and the novel genes COL13A1, KCNH3, and SENP3. Etiologically, aggressive PC was associated with age-related and DNA repair-related mutational signatures. Copy number variants most frequently affected 14q11.2 and 8p21.2, where no well-recognized PC-related genes are located, and also frequently affected regions near the known PC-related genes MYC, AR, TP53, PTEN, and BRCA1. Structural variants most frequently involved not only the known PC-related genes TMPRSS2 and ERG but also the less extensively studied gene in this context, PTPRD. Finally, clinically actionable variants were detected throughout all stages of aggressive PC and in both plasma and tissue samples, emphasizing the potential clinical applicability of WGS of minimally invasive plasma samples. Overall, our study highlights the feasibility of using liquid biopsies for comprehensive genomic characterization as an alternative to tissue biopsies in advanced/aggressive PC.

Whole-genome assembly of a novel invertebrate herpesvirus from the gastropod Babylonia areolata.

Molluscan herpesviruses cause disease in species of major importance to aquaculture and are the only known herpesviruses to infect invertebrates, which lack an adaptive immune system. Understanding the evolution of malacoherpesviruses in relation to their hosts will likely require comparative genomic studies on multiple phylogenetic scales. Currently, only two malacoherpesvirus species have genomes that have been fully assembled, which limits the ability to perform comparative genomic studies on this family of viruses. In the present study, we fully assemble a herpesvirus from Illumina and Nanopore sequence data that were previously used to assemble the genome of the gastropod Babylonia areolata. We tentatively assign this novel herpesvirus to the genus Aurivirus within the family Malacoherpesviridae based on a phylogenetic analysis of DNA polymerase. While structurally similar to other malacoherpesvirus genomes, a synteny analysis of the novel herpesvirus with another Aurivirus species indicates that genomic rearrangements might be an important process in the evolution of this genus. We anticipate that future complete assemblies of malacoherpesviruses will be a valuable resource in comparative herpesvirus research.

DNA Methylation in Ctenophores.

Epigenomic regulation and dynamic DNA methylation, in particular, are widespread mechanisms orchestrating the genome operation across time and species. Whole-genome bisulfite sequencing (WGBS) is currently the only method for unbiasedly capturing the presence of 5-methylcytosine (5-mC) DNA methylation patterns across an entire genome with single-nucleotide resolution. Bisulfite treatment converts unmethylated cytosines to uracils but leaves methylated cytosines intact, thereby creating a map of all methylated cytosines across a genome also known as a methylome. These epigenomic patterns of DNA methylation have been found to regulate gene expression and influence gene evolution rates between species. While protocols have been optimized for vertebrate methylome production, little adaptation has been done for invertebrates. Creating a methylome reference allows comparisons to be made between rates of transcription and epigenomic patterning in animals. Here we present a method of library construction for bisulfite sequencing optimized for non-bilateral metazoans such as the ctenophore, Mnemiopsis leidyi. We have improved upon our previously published method by including spike-in genomic DNA controls to measure methylation conversion rates. By pooling two bisulfite conversion reactions from the same individual, we also produced sequencing libraries that yielded a higher percentage of sequenced reads uniquely mapping to the reference genome. We successfully detected 5-mC in whole-animal methylomes at CpG, CHG, and CHH sites and visualized datasets using circos diagrams. The proof-of-concept tests were performed both under control conditions and following injury tests with changes in methylation patterns of genes encoding innexins, toxins and neuropeptides. Our approach can be easily adapted to produce epigenomes from other fragile marine animals.

Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

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.

Marker Density and Models to Improve the Accuracy of Genomic Selection for Growth and Slaughter Traits in Meat Rabbits.

The selection and breeding of good meat rabbit breeds are fundamental to their industrial development, and genomic selection (GS) can employ genomic information to make up for the shortcomings of traditional phenotype-based breeding methods. For the practical implementation of GS in meat rabbit breeding, it is necessary to assess different marker densities and GS models. Here, we obtained low-coverage whole-genome sequencing (lcWGS) data from 1515 meat rabbits (including parent herd and half-sibling offspring). The specific objectives were (1) to derive a baseline for heritability estimates and genomic predictions based on randomly selected marker densities and (2) to assess the accuracy of genomic predictions for single- and multiple-trait linear mixed models. We found that a marker density of 50 K can be used as a baseline for heritability estimation and genomic prediction. For GS, the multi-trait genomic best linear unbiased prediction (GBLUP) model results in more accurate predictions for virtually all traits compared to the single-trait model, with improvements greater than 15% for all of them, which may be attributed to the use of information on genetically related traits. In addition, we discovered a positive correlation between the performance of the multi-trait GBLUP and the genetic correlation between the traits. We anticipate that this approach will provide solutions for GS, as well as optimize breeding programs, in meat rabbits.