Genome Coverage Research Articles

Advancing equine genomics: the development of a high density Axiom_Ashwa SNP chip for Indian horses and ponies.

The unique horse and pony breeds of India are declining at an alarming rate. These horses have been integral to the Indian culture and customs for centuries and represent a valuable genetic resource. It is imperative to harness the potential of this equine genetic resource that urgently needs conservation. The study highlights the design and development of a high density SNP array, the Axiom_Ashwa to aid in the genetic analysis and conservation efforts for Indian horse and pony breeds. With 613,950 SNPs, this chip offers extensive genome coverage having an average inter-marker distance of 4kb. The Axiom_Ashwa has been validated on a larger set of diverse indigenous samples as well as Thoroughbreds, demonstrating a high call rate of 99.4% and robustness for genotyping indigenous breeds. Linkage disequilibrium (LD) analysis showed higher average LD in Indian breeds compared to exotic breeds, suggesting a limited effective population size and recent bottlenecks. Phylogenetic and population stratification analyses using PCA and DAPC clearly distinguished horses, ponies and Thoroughbreds, confirming the efficacy of the Axiom_Ashwa chip. These findings underscore the urgent need for conservation efforts for Indian horse breeds, which have experienced significant drop in population size. The Axiom_Ashwa SNP chip offers advantages such as cost-effectiveness and high throughput, providing a more accurate genetic representation of Indian horses.

Empirical assessment of the enrichment-based metagenomic methods in identifying diverse respiratory pathogens.

Probe-based nucleic acid enrichment represents an effective route to enhance the detection capacity of next-generation sequencing (NGS) in a set of clinically diverse and relevant microbial species. In this study, we assessed the effect of the enrichment-based sequencing on identifying respiratory infections using tiling RNA probes targeting 76 respiratory pathogens and sequenced using both Illumina and Oxford Nanopore platforms. Forty respiratory swab samples pre-tested for a panel of respiratory pathogens by qPCR were used to benchmark the sequencing data. We observed a general improvement in sensitivity after enrichment. The overall detection rate increased from 73 to 85% after probe capture detected by Illumina. Moreover, enrichment with probe sets boosted the frequency of unique pathogen reads by 34.6 and 37.8-fold for Illumina DNA and cDNA sequencing, respectively. This also resulted in significant improvements on genome coverage especially in viruses. Despite these advantages, we found that library pooling may cause reads mis-assignment, probably due to crosstalk issues arise from post-capture PCR and from pooled sequencing, thus increasing the risk of bleed-through signal. Taken together, an overall improvement in the breadth and depth of pathogen coverage is achieved using enrichment-based sequencing method. For future applications, automated library processing and pooling-free sequencing could enhance the precision and timeliness of probe enrichment-based clinical metagenomics.

Draft genome sequencing of multidrug-resistant Pseudomonas asiatica strains isolated from dairy cows with clinical mastitis and their farm environment.

We report, for the first time, the draft genomes of four multidrug-resistant Pseudomonas asiatica strains isolated from milk (2M1), feces (2F1 and 2F2), and farm soil (2S1) samples of dairy cows with clinical mastitis. The assembled genomes of these strains were 5.7 Mbp, 62.8% G + C content, and 50× genome coverage.

Sequence-matching adapter trimmers generate consistent quality and assembly metrics for Illumina sequencing of RNA viruses

Trimming adapters and low-quality bases from next-generation sequencing (NGS) data is crucial for optimal analysis. We evaluated six trimming programs, implementing five different algorithms, for their effectiveness in trimming adapters and improving quality, contig assembly, and single-nucleotide polymorphism (SNP) quality and concordance for poliovirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and norovirus paired data sequenced on Illumina iSeq and MiSeq platforms. Trimmomatic and BBDuk effectively removed adapters from all datasets, unlike FastP, AdapterRemoval, SeqPurge, and Skewer. All trimmers improved read quality (Q ≥ 30, 87.8 − 96.1%) compared to raw reads (83.6 − 93.2%). Trimmers implementing traditional sequence-matching (Trimmomatic and AdapterRemoval) and overlapping algorithm (FastP) retained the highest-quality reads. While all trimmers improved the maximum contig length and genome coverage for iSeq and MiSeq viral assemblies, BBDuk-trimmed reads assembled the shortest contigs. SNP concordance was consistently high (> 97.7 − 100%) across trimmers. However, BBDuk-trimmed reads had the lowest quality SNPs. Overall, the two adapter trimmers that utilized the traditional sequence-matching algorithm performed consistently across the viral datasets analyzed. Our findings guide software selection and inform future versatile trimmer development for viral genome analysis.

First Report of Paprika mild mottle virus Infecting pepper (Capsicum annuum L.) in China.

Paprika mild mottle virus (PaMMV), a tobamovirus from the Virgaviridae family, has been reported to infect pepper (Capsicum annuum L.) in the Netherlands (García-Luque et al., 1993), Japan (Hamada et al., 2003), Bulgaria (Ruíz del Pino et al., 2003), and Israel (Luria et al., 2018). Pepper has become the most widely planted vegetable and the most heavily consumed spicy condiment, with a large planting area (>21,000 km2) in China (Zou et al., 2021). However, the crop's productivity is limited by viral diseases. In August 2021, a survey of viral diseases was conducted across five pepper-growing regions in Linyi, Shandong Province, China, revealed that 5% to 10% of plants in most fields exhibited symptoms such as stunting, leaf narrowing, chlorosis, and crinkling of leaf tissue (Fig. S1-a). Subsequently, seven symptomatic samples (Fig. S1-b) were collected from a plot (118° 29' E, 34° 65' N) where thin-skinned peppers were grown. All samples tested positive in a Western blot assay using a polyclonal antibody for tobacco mosaic virus (TMV), indicating the presence of a tobamovirus (Fig. S1-c). To identify the virus, total RNAs were extracted from 10 different symptomatic pepper plant samples, using the RNeasy Plant Mini Kit (QIAGEN, Hilden, Germany) and combined into a single sample in equal amounts (100 ng/μL each). RNA-Seq was performed on this mixture using the Illumina NovaSeq 6000 (Illumina, San Diego, USA). Raw reads (6.92G) in fastq format were processed with in-house perl scripts to yield clean reads (6.68G), which were assembled into 780 contigs (mean length: 2036 bp) using Trinity software (v2.6.6). These contigs were then analysed using the NCBI BLASTX program (http://www.ncbi.nlm.nih.gov/blast) against the viral RefSeq database. The results demonstrated that the contigs exhibited an average genome coverage of 40.16%. Notably, one unique contig (6453 bp) was mapped to the PaMMV genome (reference sequence KX187305.1) with 98.34% identity. To verify the RNA-Seq result, a PaMMV-specific primer pair (forward: 5'-GAGTTCATAGAGGCAGTACC-3'; reverse: 5'-CTTCGATTTAAGTGGAGGGAT-3') was designed for amplifying an 800-nt fragment of the PaMMV coat protein (CP) gene by RT-PCR. The 800-nt fragment was successfully amplified from all the symptomatic samples (Fig. S1-d). Sequenced RT-PCR products (GenBank No. OR365081.1) showed 99.75% nucleotide identity with PaMMV isolates (OQ198318.1, China, goat) according to BLAST analysis. Mechanical transmission of PaMMV from three infected pepper plants collected in Linyi to healthy pepper plants confirmed pathogenicity. Three pepper cultivars (trade name: Qiemen, Tianyu, Haonong11) were tested and showed varying symptoms post-inoculation (Fig. S1-e). All cultivars were confirmed in the lab to be PMMoV-sensitive and exhibited significant viral disease symptoms after PMMoV inoculation (data not shown). Each variety was inoculated with at least 3 seedlings from different field-collected virulence sources. Qiemen (cultivar: sweet pepper) exhibited severe symptoms with necrotic spots, Tianyu (cultivar: Chaotian pepper) had mild leaf wrinkling, and Haonong11 (cultivar: Pickled pepper) showed yellow mottle and leaf shrinking. All tested seedlings were positive for PaMMV by RT-PCR (Fig. S1-f). To our knowledge, this is the first report of PaMMV infecting pepper in China, highlighting the need for vigilant monitoring to protect the chili industry and prevent the virus's spread to other crops.

Analysis of Genetic Diversity of Some Olive Cultivars Olea Europoea L. Using ISSR, SSR

In the current investigation, the genetic relations, and the genetic dimension of seven olive varieties (Olea europaea L.) were identified using, ISSR, and SSR markers. The study involved uses leaves from different samples for DNA extraction. The DNA amount ranged from 150 to 400 micrograms with purity ranging between 1.6 to 1.9. Agarose gel electrophoresis is applied to evaluate PCR reaction success. The ISSR marker, which used 4 primers, showed different results for the multiplication. Various bands were observed that differed from each other, resulting in a total of 93 bands - 2 general and 91 different. The primer also distinguished unique bands, with 4 unique bands and 7 absent bands. The primer with the highest molecular size (1500bp) was UBC-817, while the lowest molecular volume (200bp) was observed in the UBC-826 primer. The marker results indicated that the two types, Suranie (3) and Frantoio (4), had the lowest hereditary dimension of (0.030). The Frantoio and Santacatrina varieties showed the highest genetic difference, with a value of 0.622. According to the genetic relationship analysis, there were three main groups. The first group consisted of a single category (5), while There were two subgroups in the second group, B1 and B2. Subgroup B1 had only one category, while the rest of the categories were in Subgroup B2.. Even though these markers use different mechanisms to detect variance and genome coverage, they complement each other. Moreover, five SSR markers were used in this study to describe seven olive cultivars and analyze their genetic relationship. Four out of five primers showed positive results, while one was not present. The SSR markers were effective in identifying the similarity of collected species, as they are specialized indicators ISSR markers.

Genome resequencing and genome-wide polymorphisms in mosquito vectors Aedes aegypti and Aedes albopictus from south India

Aedes aegypti and Aedes albopictus mosquitoes spread major vector-borne viral diseases in tropical and sub-tropical regions of the globe. In this study, we sequenced the genome of Indian Ae. aegypti and Ae. albopictus and mapped to their reference genomes. Comparative genomics were performed between our strain and the reference strains. A total of 14,416,484 single nucleotide polymorphisms (SNPs) and 156,487 insertions and deletions (InDels) were found in Ae. aegypti, and 28,940,433 SNPs and 188,987 InDels in Ae. albopictus. Particular emphasis was given to gene families involved in mosquito digestion, development, and innate immunity, which could be putative candidates for vector control. Serine protease cascades and their inhibitors called serpins, play a central role in these processes. We extracted high-impact variants in genes associated with serine proteases and serpins. This study reports for the first time a high coverage genome sequence data of an Indian Ae. albopictus mosquito. The results from this study will provide insights into Indian Aedes specific polymorphisms and the evolution of immune related genes in mosquitoes, which can serve as a resource for future comparative genomics and those pursuing the development of targeted biopesticides for effective mosquito control strategies.

MVP: a modular viromics pipeline to identify, filter, cluster, annotate, and bin viruses from metagenomes.

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.IMPORTANCEThe significance of our work lies in the development of Modular Viromics Pipeline (MVP), an integrated and user-friendly pipeline tailored exclusively for viromics analyses. MVP stands out due to its modular design, which ensures easy installation, execution, and integration of new tools and databases. By combining state-of-the-art tools such as geNomad and CheckV, MVP provides high-quality viral genome recovery and taxonomy and host assignment, and functional annotation, addressing the limitations of existing pipelines. MVP's ability to handle diverse sample types, including environmental, human microbiome, and plant-associated samples, makes it a versatile tool for the broader microbiome research community. By standardizing the analysis process and providing easily interpretable results, MVP enables researchers to perform comprehensive studies of viral communities, significantly advancing our understanding of viral ecology and its impact on various ecosystems.

A novel isothermal whole genome sequencing approach for Monkeypox Virus

Monkeypox virus (MPXV) is the zoonotic agent responsible for mpox, an often-self-limiting pox-like disease. Since May 2022, an outbreak characterized by increased human-to-human transmission was detected outside the endemic regions. Whole genome sequencing (WGS) has been successfully used to keep track of viral evolution during outbreaks or for surveillance of multiple pathogens of public health significance. Current WGS protocols for MPXV are either based on metagenomic sequencing or tiled-PCR amplification. The latter allows multiplexing due to the efficient enrichment of the viral DNA, however, mutations or the presence of different clades can negatively influence genome coverage yield. Here, we present the establishment of a novel isothermal WGS method for MPXV based on Phi29 DNA polymerase-based multiple displacement amplification (MDA) properties making use of only 6 primers. This approach yielded from 88% up to 100% genome coverage using either alkaline denatured extracted DNA or clinical material as starting material, with the highest coverage generated by clinical material. We demonstrate that this novel isothermal WGS protocol is suitable for monitoring viral evolution during MPXV outbreaks and surveillance in any conventional laboratory setting.

Construction and characterization of DNA libraries from cultured phages and environmental viromes.

Despite many efforts to understand and leverage the functional potential of environmental viromes, most bacteriophage genes are largely uncharacterized. To explore novel biology from uncultivated microbes like phages, metagenomics has emerged as a powerful tool to directly mine new genes without the need to culture the diverse microbiota and the viruses within. When a pure computational approach cannot infer gene function, it may be necessary to create a DNA library from environmental genomic DNA, followed by the screening of that library for a particular function. However, these screens are often initiated without a metagenomic analysis of the completed DNA library being reported. Here, we describe the construction and characterization of DNA libraries from a single cultured phage (ΦT4), five cultured Escherichia coli phages, and three metagenomic viral sets built from freshwater, seawater, and wastewater samples. Through next-generation sequencing of five independent samplings of the libraries, we found a consistent number of recovered genes per replicate for each library, with many genes classifiable via the KEGG and Pharokka databases. By characterizing the size of the genes and inserts, we found that our libraries contain a median of one to two genes per contig with a median gene length of 303-381 bp for all libraries, reflective of the small genomes of viruses. The environmental libraries were genetically diverse compared to the single phage and multi-phage libraries. Additionally, we found reduced coverage of individual genomes when five phages were used as opposed to one. Taken together, this work provides a comprehensive analysis of the DNA libraries from phage genomes that can be used for metagenomic exploration and functional screens to infer and identify new biology.IMPORTANCEFunctional metagenomics is an approach that aims to characterize the putative biological function of genes in the microbial world. This includes an examination of the sequencing data collected from a pooled source of diverse microbes and inference of gene function by comparison to annotated and studied genes from public databases. At times, DNA libraries are made from these genes, and the library is screened for a specific function. Hits are validated using a combination of biological, computational, and structural analysis. Left unresolved is a detailed characterization of the library, both its diversity and content, for the purposes of imputing function entirely by computational means, a process that may yield findings that aid in designing useful screens to identify novel gene functions. In this study, we constructed libraries from cultured phages and uncultured viromes from the environment and characterized some important parameters, such as gene number, genes per contig, ratio of hypothetical to known proteins, total genomic coverage and recovery, and the effect of pooling genetic information from multiple sources, to provide a better understanding of the nature of these libraries. This work will aid the design and implementation of future screens of pooled DNA libraries to discover and isolate viral genes with novel biology across various biomes.

Optimized genome-wide CRISPR screening enables rapid engineering of growth-based phenotypes in Yarrowia lipolytica

CRISPR-Cas9 functional genomic screens uncover gene targets linked to various phenotypes for metabolic engineering with remarkable efficiency. However, these genome-wide screens face a number of design challenges, including variable guide RNA activity, ensuring sufficient genome coverage, and maintaining high transformation efficiencies to ensure full library representation. These challenges are prevalent in non-conventional yeast, many of which exhibit traits that are well suited to metabolic engineering and bioprocessing. To address these hurdles in the oleaginous yeast Yarrowia lipolytica, we designed a compact, high-activity genome-wide sgRNA library. The library was designed using DeepGuide, a sgRNA activity prediction algorithm and a large dataset of ∼50,000 sgRNAs with known activity. Three guides per gene enables redundant targeting of 98.8% of genes in the genome in a library of 23,900 sgRNAs. We deployed the optimized library to uncover genes essential to the tolerance of acetate, a promising alternative carbon source, and various hydrocarbons present in many waste streams. Our screens yielded several gene knockouts that improve acetate tolerance on their own and as double knockouts in media containing acetate as the sole carbon source. Analysis of the hydrocarbon screens revealed genes related to fatty acid and alkane metabolism in Y. lipolytica. The optimized CRISPR gRNA library and its successful use in Y. lipolytica led to the discovery of alternative carbon source-related genes and provides a workflow for creating high-activity, compact genome-wide libraries for strain engineering.

AMAR-seq: Automated Multimodal Sequencing of DNA Methylation, Chromatin Accessibility, and RNA Expression with Single-Cell Resolution.

Parallel single-cell multimodal sequencing is the most intuitive and precise tool for cellular status research. In this study, we propose AMAR-seq to automate methylation, chromatin accessibility, and RNA expression coanalysis with single-cell precision. We validated the accuracy and robustness of AMAR-seq in comparison with standard single-omics methods. The high gene detection rate and genome coverage of AMAR-seq enabled us to establish a genome-wide gene expression regulatory atlas and triple-omics landscape with single base resolution and implement single-cell copy number variation analysis. Applying AMAR-seq to investigate the process of mouse embryonic stem cell differentiation, we revealed the dynamic coupling of the epigenome and transcriptome, which may contribute to unraveling the molecular mechanisms of early embryonic development. Collectively, we propose AMAR-seq for the in-depth and accurate establishment of single-cell multiomics regulatory patterns in a cost-effective, efficient, and automated manner, paving the way for insightful dissection of complex life processes.

Reconstructing SARS-CoV-2 lineages from mixed wastewater sequencing data

Wastewater surveillance of SARS-CoV-2 has emerged as a critical tool for tracking the spread of COVID-19. In addition to estimating the relative case numbers using quantitative PCR, SARS-CoV-2 genomic RNA can be extracted from wastewater and sequenced. There are many existing techniques for using the sequenced RNA to determine the relative abundance of known lineages in a sample. However, it is very challenging to predict novel lineages from wastewater data due to its mixed composition and unreliable genomic coverage. In this work, we present a novel technique based on non-negative matrix factorization which is able to reconstruct lineage definitions by analyzing data from across different samples. We test the method both on synthetic and real wastewater sequencing data. We show that the technique is able to determine major lineages such as Omicron and Delta as well as sub-lineages such as BA.5.2.1. We provide a method for determining emerging lineages in wastewater without the need for genomic data from clinical samples. This could be used for routine monitoring of SARS-CoV-2 as well as other emerging viral pathogens in wastewater. Additionally, it may be used to determine more full-genome sequences for viruses with fewer available genomes.

Insights from homozygous signatures of cervus nippon revealed genetic architecture for components of fitness.

This study investigates the genomic landscape of Sika deer populations, emphasizing the detection and characterization of runs of homozygosity (ROH) and their contribution towards components of fitness. Using 85,001 high-confidence SNPs, the investigation into ROH distribution unveiled nuanced patterns of autozygosity across individuals especially in 2 out of the 8 farms, exhibiting elevated ROH levels and mean genome coverage under ROH segments. The prevalence of shorter ROH segments (0.5-4Mb) suggests historical relatedness and potential selective pressures within these populations. Intriguingly, despite observed variations in ROH profiles, the overall genomic inbreeding coefficient (FROH) remained relatively low across all farms, indicating a discernible degree of genetic exchange and effective mitigation of inbreeding within the studied Sika deer populations. Consensus ROH (cROH) were found to harbor genes for important functions viz., EGFLAM gene which is involved in the vision function of the eye, SKP2 gene which regulates cell cycle, CAPSL involved in adipogenesis, SPEF2 which is essential for sperm flagellar assembly, DCLK3 involved in the heat stress. This first ever study on ROH in Sika deer, to shed light on the adaptive role of genes in these homozygous regions. The insights garnered from this study have broader implications in the management of genetic diversity in this vulnerable species.

Telomere-to-Telomere Haplotype-Resolved Genomes of Agrocybe chaxingu Reveals Unique Genetic Features and Developmental Insights.

Agrocybe chaxingu is a widely cultivated edible fungus in China, which is rich in nutrients and medicinal compounds. However, the lack of a high-quality genome hinders further research. In this study, we assembled the telomere-to-telomere genomes of two sexually compatible monokaryons (CchA and CchB) derived from a primarily cultivated strain AS-5. The genomes of CchA and CchB were 50.60 Mb and 51.66 Mb with contig N50 values of 3.95 Mb and 3.97 Mb, respectively. Each contained 13 complete chromosomes with telomeres at both ends. The high mapping rate, uniform genome coverage, high LAI score, all BUSCOs with 98.5%, and all base accuracy exceeding 99.999% indicated the high level of integrity and quality of these two assembled genomes. Comparison of the two genomes revealed that approximately 30% of the nucleotide sequences between homologous chromosomes were non-syntenic, including 19 translocations, 36 inversions, and 15 duplications. An additional gene CchA_000467 was identified at the Mat A locus of CchA, which was observed exclusively in the Cyclocybe cylindracea species complex. A total of 613 (4.26%) and 483 (3.4%) unique genes were identified in CchA and CchB, respectively, with over 80% of these being hypothetical proteins. Transcriptomic analysis revealed that the expression levels of unique genes in CchB were significantly higher than those in CchA, and both CchA and CchB had unique genes specifically expressed at stages of mycelium and fruiting body. It was indicated that the growth and development of the A. chaxingu strain AS-5 required the coordinated action of two different nuclei, with CchB potentially playing a more significant role. These findings contributed to a more profound comprehension of the growth and developmental processes of basidiomycetes.

Evaluation of Hi-C Sequencing for Detection of Gene Fusions in Hematologic and Solid Tumor Pediatric Cancer Samples.

Hi-C sequencing is a DNA-based next-generation sequencing method that preserves the 3D genome conformation and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate Hi-C as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens. Archived viable and non-viable frozen leukemic cells and formalin-fixed paraffin-embedded (FFPE) tumor specimens were analyzed. Pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to Hi-C to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases without known genomic rearrangements based on prior clinical diagnostic testing was evaluated to determine whether Hi-C could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, we observed 100% concordance between Hi-C and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study provides an institutional proof of principle evaluation of Hi-C sequencing to medical diagnostic testing as it identified several clinically relevant rearrangements, including those that were missed by current clinical testing workflows.

Capturing clinically relevant Campylobacter attributes through direct whole genome sequencing of stool.

Campylobacter is the leading bacterial cause of infectious intestinal disease, but the pathogen typically accounts for a very small proportion of the overall stool microbiome in each patient. Diagnosis is even more difficult due to the fastidious nature of Campylobacter in the laboratory setting. This has, in part, driven a change in recent years, from culture-based to rapid PCR-based diagnostic assays which have improved diagnostic detection, whilst creating a knowledge gap in our clinical and epidemiological understanding of Campylobacter genotypes - no isolates to sequence. In this study, direct metagenomic sequencing approaches were used to assess the possibility of replacing genome sequences with metagenome sequences; metagenomic sequencing outputs were used to describe clinically relevant attributes of Campylobacter genotypes. A total of 37 diarrhoeal stool samples with Campylobacter and five samples with an unknown pathogen result were collected and processed with and without filtration, DNA was extracted, and metagenomes were sequenced by short-read sequencing. Culture-based methods were used to validate Campylobacter metagenome-derived genome (MDG) results. Sequence output metrics were assessed for Campylobacter genome quality and accuracy of characterization. Of the 42 samples passing quality checks for analysis, identification of Campylobacter to the genus and species level was dependent on Campylobacter genome read count, coverage and genome completeness. A total of 65% (24/37) of samples were reliably identified to the genus level through Campylobacter MDG, 73% (27/37) by culture and 97% (36/37) by qPCR. The Campylobacter genomes with a genome completeness of over 60% (n=21) were all accurately identified at the species level (100%). Of those, 72% (15/21) were identified to sequence types (STs), and 95% (20/21) accurately identified antimicrobial resistance (AMR) gene determinants. Filtration of stool samples enhanced Campylobacter MDG recovery and genome quality metrics compared to the corresponding unfiltered samples, which improved the identification of STs and AMR profiles. The phylogenetic analysis in this study demonstrated the clustering of the metagenome-derived with culture-derived genomes and revealed the reliability of genomes from direct stool sequencing. Furthermore, Campylobacter genome spiking percentages ranging from 0 to 2% total metagenome abundance in the ONT MinION sequencer, configured to adaptive sequencing, exhibited better assembly quality and accurate identification of STs, particularly in the analysis of metagenomes containing 2 and 1% of Campylobacter jejuni genomes. Direct sequencing of Campylobacter from stool samples provides clinically relevant and epidemiologically important genomic information without the reliance on cultured genomes.

Evidence of Grapevine enamovirus 1 Infecting Grapevine (Vitis vinifera L.) in Canada.

Grapevine enamovirus 1 (GEV1) belongs to the genus Enamovirus, in the family Solemoviridae. It has been reported from several countries infecting grapevines including Brazil (Silva et al. 2017), China (Ren et al. 2021) and France (Hily et al. 2022). To assess the prevalence and diversity of economically important grapevine viruses in nine Canadian vineyards, total RNA and double-stranded RNA (dsRNA) (Fall et al. 2020) were extracted from 30 and 100 composite samples respectively, with each consisting of five vines of the same cultivars. The cultivars included in this study are Frontenac noir (n=34), Vidal (n=32), Marquette (n=33), Riesling (n=31), and Pinot noir (n=31). The total RNA and dsRNA samples were subsequently multiplexed and diagnosed by high-throughput sequencing (HTS) on NovaSeq (600 S4 PE100) and MiSeq (2 × 250 cycle PE) respectively. From NovaSeq and MiSeq sequencing, an average of 410,000 to 1.3 million reads/sample were obtained, respectively, with mapped viral reads representing 10.92% to 12.48% of the total reads. After sequence quality was verified using Trimmomatic v.0.40 (Bolger et al. 2014), the clean sequences were screened against all possible viruses in the databases using the Virtool (Rott et al. 2017) and VirFind virus detection pipelines (Ho and Tzanetakis 2014). GEV1 was detected in clean sequences from two, three, and two leaf samples of cultivars 'Marquette' 'Riesling' and 'Frontenac noir' respectively. Six of the seven HTS-assembled GEV1 genomes were partial, ranging from 4,523 to 6,000 nucleotide (nt) with genome coverage varying from 71% to 89%. Only one 6,314 nt long assembled contig (Accession No. OR021829), represented a nearly complete genome, being only 53 and 3 nt shorter than Sd-CG (MT536978) at 5' and 3' untranslated regions (UTR), respectively. Isolate 3- Riesling-CAN (OR021829) shares 90.56 to 94.19% nt identities with several GEV1isolates at 96-99% of query coverage. Phylogenetically, OR021829 is closer to GEV1 isolates from France and China (Figure S1). To validate the HTS results, the developed primer pair SetF and Set1R (Silva et al., 2017) was used for RT-PCR detection. The amplicons from all seven HTS-positive samples were sequenced using Sanger sequencing, confirming the presence of GEV-1 in three studied grape cultivars in Canadian vineyards. Symptoms associated with the specific GEV1-infected vines could not be explained as composite samples were used. Each of the combined samples HTS library also tested positive for at least one of the known grape virus/viroids, namely grapevine leafroll associated-virus -3, grapevine pinot gris virus, grapevine rupestris stem pitting-associated virus, Marafivirus syrahense grapevine Syrah virus-1 and hop stunt viroid. To our knowledge, this is the first report of GEV1 being detected in grapevines in Canada, or in any North American vineyard. GEV1 is a relatively new virus, and its biology remains largely unknown. Based on this sequence new GEV1 primers can be developed to know the genetic variability among GEV-1 and improve the detection of this virus in vineyards.

VISTA: an integrated framework for structural variant discovery.

Structural variation (SV) refers to insertions, deletions, inversions, and duplications in human genomes. SVs are present in approximately 1.5% of the human genome. Still, this small subset of genetic variation has been implicated in the pathogenesis of psoriasis, Crohn's disease and other autoimmune disorders, autism spectrum and other neurodevelopmental disorders, and schizophrenia. Since identifying structural variants is an important problem in genetics, several specialized computational techniques have been developed to detect structural variants directly from sequencing data. With advances in whole-genome sequencing (WGS) technologies, a plethora of SV detection methods have been developed. However, dissecting SVs from WGS data remains a challenge, with the majority of SV detection methods prone to a high false-positive rate, and no existing method able to precisely detect a full range of SVs present in a sample. Previous studies have shown that none of the existing SV callers can maintain high accuracy across various SV lengths and genomic coverages. Here, we report an integrated structural variant calling framework, Variant Identification and Structural Variant Analysis (VISTA), that leverages the results of individual callers using a novel and robust filtering and merging algorithm. In contrast to existing consensus-based tools which ignore the length and coverage, VISTA overcomes this limitation by executing various combinations of top-performing callers based on variant length and genomic coverage to generate SV events with high accuracy. We evaluated the performance of VISTA on comprehensive gold-standard datasets across varying organisms and coverage. We benchmarked VISTA using the Genome-in-a-Bottle gold standard SV set, haplotype-resolved de novo assemblies from the Human Pangenome Reference Consortium, along with an in-house polymerase chain reaction (PCR)-validated mouse gold standard set. VISTA maintained the highest F1 score among top consensus-based tools measured using a comprehensive gold standard across both mouse and human genomes. VISTA also has an optimized mode, where the calls can be optimized for precision or recall. VISTA-optimized can attain 100% precision and the highest sensitivity among other variant callers. In conclusion, VISTA represents a significant advancement in structural variant calling, offering a robust and accurate framework that outperforms existing consensus-based tools and sets a new standard for SV detection in genomic research.

Beyond the base pairs: comparative genome-wide DNA methylation profiling across sequencing technologies.

Deoxyribonucleic acid (DNA) methylation plays a key role in gene regulation and is critical for development and human disease. Techniques such as whole-genome bisulfite sequencing (WGBS) and reduced representation bisulfite sequencing (RRBS) allow DNA methylation analysis at the genome scale, with Illumina NovaSeq 6000 and MGI Tech DNBSEQ-T7 being popular due to their efficiency and affordability. However, detailed comparative studies of their performance are not available. In this study, we constructed 60 WGBS and RRBS libraries for two platforms using different types of clinical samples and generated approximately 2.8 terabases of sequencing data. We systematically compared quality control metrics, genomic coverage, CpG methylation levels, intra- and interplatform correlations, and performance in detecting differentially methylated positions. Our results revealed that the DNBSEQ platform exhibited better raw read quality, although base quality recalibration indicated potential overestimation of base quality. The DNBSEQ platform also showed lower sequencing depth and less coverage uniformity in GC-rich regions than did the NovaSeq platform and tended to enrich methylated regions. Overall, both platforms demonstrated robust intra- and interplatform reproducibility for RRBS and WGBS, with NovaSeq performing better for WGBS, highlighting the importance of considering these factors when selecting a platform for bisulfite sequencing.