Sequencing Depth Coverage Research Articles

P-528 Establishing sequencing coverage depth for preimplantation genetic testing for monogenic disorders (PGT-M) using a specific targeted next generation sequencing (tNGS)

Abstract Study question What is the sufficient sequencing coverage read depth (RD) for direct variant analysis for PGT-M using a specific tNGS assay? Summary answer Downsampling modelling experiments identified a minimum sequencing depth of 50X that is recommended for PGT-M with this specific tNGS assay. What is known already Properly developed and validated multiplex PCR amplification and NGS-based analysis reliably produces PGT-M diagnoses with low allele and amplification failure at RD of hundreds of reads when analyzing a 3-10 cell trophectoderm (TE) biopsy. However, the RD can be affected by amplification efficiency in a specific genomic region and therefore affect the reproducibility of variant detection: the higher the number of aligned reads, the higher the confidence of base calling at a particular position. Downsampling is a powerful bioinformatics tool for modelling comparisons between samples or between replicates, identifying analytical sequencing depth limits for accurate clinical diagnosis. Study design, size, duration This is a retrospective comparative diagnostic test accuracy (DTA) and verification study. Twenty-three embryos from 8 patients that had undergone PGT-M with clinically established Taqman qPCR genotyping were used as positive controls. The primary outcome was to establish the lowest RD point of down sampling where allelic frequency is reliably detected at heterozygosity between 20-80% in known positive controls. A further 30 PGT-M clinical cases were used to verify analytical performances under clinical conditions. Participants/materials, setting, methods Embryo biopsies were pre-amplified with a multiplex primer pool for PGT-M and PGT-A, using a two-step PCR strategy. Sequencing was performed on Illumina NextSeq550. Reads were aligned to a human reference genome (GRCh37/hg19) and processed with an in house bioinformatic pipeline. The RD and alternative allele frequencies (AAF) for PGT-M was calculated. The reads for the heterozygous variants were downsampled to 10X, 20X, 30X, 50X, 100X, 200X, 300X, 500X, and 1000X to establish the cutoff. Main results and the role of chance One or two TE biopsies (N = 44) of 23 embryos that were 100% concordant with previous PGT-M results were tested for 1 to 3 of the 12 variants (MLH1 c.199G>A, RB1 c.1589A>G, SLC16A2 c.604G>A, GNPTAB c.1123C>T, GNPTAB c.3314A>G, WFS1 c.1391C>G, USH2A c.10073G>A, USH2A del exons 5-11, F5 c.1601G>A, GAA c.693-1G>C, GAA c.875A>G, and HEXA c.1274_1277dupTATC). PGT-A analysis of multiple biopsies from the same embryo demonstrated consistent karyotypes (N = 44). The 19 heterozygous variants from 17 biopsies exhibited average RD of 763X (238X -1980X) with average AAF of 50% (44%- 56%). The reads for the heterozygous variants were downsampled to different levels (100,000 replicates). The distribution of the AAF were compared among the different depth levels. At a RD of 50X, the AAF displayed an interquartile range (IQR) of 8% (46%-54%) between the maximum and minimum values from 20% to 78%. Compared to 30X (IQR: 43%-57%, with minimum 13% and maximum 90%), and 100X (IQR: 47%-54%, with minimum 30% and maximum 71%), the RD at 50X generated reliable heterozygous calls with minimum depth requirement. 245 biopsies from 30 clinic cases were analyzed by tNGS-based PGT-M. 90 samples with heterozygous variant calls showed average RD of 979X (57X-5803X) with average AAF of 47% (29%-585). Limitations, reasons for caution The 50X cutoff is based on the heterozygous variants with theoretical 50% AFF and is specific for this PGT-A assay. This evaluation should be corroborated for complicated variants with homologous genes, which might require different thresholds. The aneuploidy at specific chromosome can cause the deviation of AFF from 50%. Wider implications of the findings This specific tNGS-based PGT-M can provide direct variant analysis at a 50X read depth with a high level of consistency compared to established methodologies while allowing for simultaneous aneuploidy screening using the same biopsy and procedural optimization. Trial registration number not applicable

Placental somatic mutation in human stillbirth and live birth: A pilot case-control study of paired placental, fetal, and maternal whole genomes

IntroductionA high frequency of single nucleotide somatic mutations in the placenta has been recently described, but its relationship to placental dysfunction is unknown. MethodsWe performed a pilot case-control study using paired fetal, maternal, and placental samples collected from healthy live birth controls (n = 10), live births with fetal growth restriction (FGR) due to placental insufficiency (n = 7), and stillbirths with FGR and placental insufficiency (n = 11). We quantified single nucleotide and structural somatic variants using bulk whole genome sequencing (30-60X coverage) in four biopsies from each placenta. We also assessed their association with clinical and histological evidence of placental dysfunction. ResultsSeventeen pregnancies had sufficiently high-quality placental, fetal, and maternal DNA for analysis. Each placenta had a median of 473 variants (range 111–870), with 95 % arising in just one biopsy within each placenta. In controls, live births with FGR, and stillbirths, the median variant counts per placenta were 514 (IQR 381–779), 582 (450–735), and 338 (245–441), respectively. After adjusting for depth of sequencing coverage and gestational age at birth, the somatic mutation burden was similar between groups (FGR live births vs. controls, adjusted diff. 59, 95 % CI -218 to +336; stillbirths vs controls, adjusted diff. −34, −351 to +419), and with no association with placental dysfunction (p = 0.7). DiscussionWe confirmed the high prevalence of somatic mutation in the human placenta and conclude that the placenta is highly clonal. We were not able to identify any relationship between somatic mutation burden and clinical or histologic placental insufficiency.

Abstract 405: Translocation detection in cancer using low-pass pore-c sequencing

Abstract Complex chromosomal rearrangements, including translocations, play a critical role in oncogenesis and are often identified as recurrent genetic aberrations in hematologic malignancies and solid tumors[1]. Translocation detection by karyotyping of an individual’s metaphase chromosomes remains challenging for detection of balanced rearrangements, which have no gain or loss of genetic material. Rearrangement detection using whole-genome sequencing is a viable alternative, but it relies on breakpoint-spanning reads and requires high sequencing depth (30-60x) with long reads to achieve high sensitivity, especially in repetitive regions of the genome[2]. Here, we describe a workflow for translocation detection at low sequencing depth with Oxford Nanopore’s long read chromatin conformation capture technique, Pore-C[3]. Importantly, for translocation calling, Pore-C does not rely on breakpoint-spanning reads but rather the high intrachromosomal interaction frequency of genomic regions around the breakpoint, thus largely lowering sequencing depth requirements and reducing mapping issues in low complexity regions. We prepared, barcoded, and sequenced Pore-C libraries of 3 cancer cell lines including lung cancer (A549), Acute Monocytic Leukemia (THP-1), melanoma (COLO 829) and its matching normal (COLO 829 BL) on a single Q20+ MinION flow cell (FLO-MIN114). Genome-wide contact maps of low-pass data (<1.5 Gbps per sample) were compared against higher depth (>30 Gbps per sample) Hi-C maps [4], revealing that low-pass Pore-C successfully captured large-scale genomic rearrangement. For example, chromosomal pairs chr8-chr11 and chr15-chr19 in A549, chr9-chr11 and chr1-chr20 in THP1, and chr7-chr15 in Colo829 were detected at <0.5X depth of sequencing coverage. These translocation events were further validated by breakpoint analysis using adaptive sampling, where targeted regions were previously confirmed by independent studies[4],[5]. Low-pass Pore-C detects translocations in an unbiased manner and does not require prior knowledge of the translocation structure. Combined with its simple sample-preparation workflow and the capability to provide genome-wide copy-number information in a single experiment, it can serve as a cost-effective and comprehensive tool for cancer genomic studies. [1] Chromosomal translocations in human cancer, Nature, 372, 143 (1994) [2] Hi-C as a tool for precise detection and characterisation of chromosomal rearrangements and copy number variation in human tumors, Genome Biology, 18, 125 (2017) [3] Identifying synergistic high-order 3D chromatin conformations from genome-scale nanopore concatemer sequencing, Nature Biotechnology, 40, 1488(2022) [4] Chromosomal translocations detection in cancer cells using chromosomal conformation capture data, Genes (Basel), 13, 7 (2022) [5] A multi-platform reference for somatic structural variation detection, Cell Genomics, 2, 6 (2022) Citation Format: Scott Hickey, Xiaoguang Dai, Sergey Aganezov, John Beaulaurier, Eoghan Harrington, Sissel Juul. Translocation detection in cancer using low-pass pore-c sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 405.

Generation of High-Quality African Swine Fever Virus Complete Genome from Field Samples by Next-Generation Sequencing.

African swine fever (ASF) is a lethal contagious viral disease of domestic pigs and wild boars caused by the African swine fever virus (ASFV). The pandemic spread of ASF has caused severe effects on the global pig industry. Whole-genome sequencing provides crucial information for virus strain characterization, epidemiology analysis and vaccine development. Here, we evaluated the performance of next-generation sequencing (NGS) in generating ASFV genome sequences from clinical samples. Thirty-four ASFV-positive field samples including spleen, lymph node, lung, liver and blood with a range of Ct values from 14.73 to 25.95 were sequenced. For different tissue samples collected from the same sick pigs, the proportion of ASFV reads obtained from the spleen samples was 3.69-9.86 times higher than other tissues. For the high-viral-load spleen samples (Ct < 20), a minimum of a 99.8% breadth of ≥10× coverage was revealed for all the samples. For the spleen samples with Ct ≥ 20, 6/12 samples had a minimum of a 99.8% breadth of ≥10× coverage. A high average depth of sequencing coverage was also achieved from the blood samples. According to our results, high-quality ASFV whole-genome sequences could be obtained from the spleen or blood samples with Ct < 20. The high-quality ASFV genome sequence generated in this study was further used for the high-resolution phylogenetic analysis of the ASFV genomes in the early stage of the ASF epidemic in China. Our study demonstrates that NGS may act as a useful tool for efficient ASFV genome characterization, providing valuable information for disease control.

Next Generation Sequencing of Red Blood Cell Antigens in Transfusion Medicine: Systematic Review and Meta-Analysis

Molecular analysis of blood groups is important in transfusion medicine, allowing the prediction of red blood cell (RBC) antigens. Many blood banks use single nucleotide variant (SNV) based methods for blood group analysis. While this is a well-established approach, it is limited to the polymorphisms included in genotyping panels. Thus, variants that alter antigenic expression may be ignored, resulting in incorrect prediction of phenotypes. The popularization of next-generation sequencing (NGS) has led to its application in transfusion medicine, including for RBC antigens determination. The present review/meta-analysis aimed to evaluate the applicability of the NGS for the prediction of RBC antigens. A systematic review was conducted following a comprehensive literature search in accordance with the Preferred Reporting Items for Systematic Review and Meta-Analysis guidelines. Studies were selected based on predefined criteria and evaluated using Strengthening the Reporting of Observational studies in Epidemiology guidelines. The characteristics and results of the studies were extracted and meta-analysis was performed to verify the agreement between results from standard molecular methods and NGS. Kell (rs8176058), Duffy (rs2814778, rs12078), or Kidd (rs1085396) alleles were selected as a model for comparisons. Additionally, results are presented for other blood group systems. Of the 864 eligible studies identified, 10 met the inclusion criteria and were selected for meta-analysis. The pooled concordance proportion for NGS compared to other methods ranged from 0.982 to 0.994. The sequencing depth coverage was identified as crucial parameters for the reliability of the results. Some studies reported difficulty in analyzing more complex systems, such as Rh and MNS, requiring the adoption of specific strategies. NGS is a technology capable of predicting blood group phenotypes and has many strengths such as the possibility of simultaneously analyzing hundred individuals and gene regions, and the ability to provide comprehensive genetic analysis, which is useful in the description of new alleles and a better understanding of the genetic basis of blood groups. The implementation of NGS in the routine of blood banks depends on several factors such as cost reduction, the availability of widely validated panels, the establishment of clear quality parameters and access to bioinformatics analysis tools that are easy to access and operate.

Integrating Internal Fragments in the Interpretation of Top-Down Sequencing Data of Larger Oligonucleotides.

In the context of direct top-down analysis or concerted bottom-up characterization of nucleic acid samples, the waning yield of terminal fragments as a function of precursor ion size poses a significant challenge to the gas-phase sequencing of progressively larger oligonucleotides. In this report, we examined the behavior of oligoribonucleotide samples ranging from 20 to 364 nt upon collision-induced dissociation (CID). The experimental data showed a progressive shift from terminal to internal fragments as a function of size. The systematic evaluation of experimental factors, such as collision energy, precursor charge, sample temperature, and the presence of chaotropic agents, showed that this trend could be modestly alleviated but not suppressed. This inexorable effect, which has been reported also for other activation techniques, prompted a re-examination of the features that have traditionally discouraged the utilization of internal fragments as a source of sequence information in data interpretation procedures. Our simulations highlighted the ability of internal fragments to produce self-consistent ladders with either end corresponding to each nucleotide in the sequence, which enables both proper alignment and correct recognition of intervening nucleotides. In turn, contiguous ladders display extensive overlaps with one another and with the ladders formed by terminal fragments, which unambiguously constrain their mutual placement within the analyte sequence. The experimental data borne out the predictions by showing ladders with extensive overlaps, which translated into uninterrupted "walks" covering the entire sequence with no gaps from end to end. More significantly, the results showed that combining the information afforded by internal and terminal ladders resulted in much a greater sequence coverage and nucleotide coverage depth than those achievable when either type of information was considered separately. The examination of a series of 58-mer oligonucleotides with high sequence homology showed that the assignment ambiguities engendered by internal fragments did not significantly exceed those afforded by the terminal ones. Therefore, the balance between potential benefits and perils of including the former makes a compelling argument for the development of integrated data interpretation strategies, which are better equipped for dealing with the changing fragmentation patterns obtained from progressively larger oligonucleotides.

Oxford nanopore long-read sequencing enables the generation of complete bacterial and plasmid genomes without short-read sequencing.

Genome-based analysis is crucial in monitoring antibiotic-resistant bacteria (ARB)and antibiotic-resistance genes (ARGs). Short-read sequencing is typically used to obtain incomplete draft genomes, while long-read sequencing can obtain genomes of multidrug resistance (MDR) plasmids and track the transmission of plasmid-borne antimicrobial resistance genes in bacteria. However, long-read sequencing suffers from low-accuracy base calling, and short-read sequencing is often required to improve genome accuracy. This increases costs and turnaround time. In this study, a novel ONT sequencing method is described, which uses the latest ONT chemistry with improved accuracy to assemble genomes of MDR strains and plasmids from long-read sequencing data only. Three strains of Salmonella carrying MDR plasmids were sequenced using the ONT SQK-LSK114 kit with flow cell R10.4.1, and de novo genome assembly was performed with average read accuracy (Q > 10) of 98.9%. For a 5-Mb-long bacterial genome, finished genome sequences with accuracy of >99.99% could be obtained at 75× sequencing coverage depth using Flye and Medaka software. Thus, this new ONT method greatly improves base-calling accuracy, allowing for the de novo assembly of high-quality finished bacterial or plasmid genomes without the need for short-read sequencing. This saves both money and time and supports the application of ONT data in critical genome-based epidemiological analyses. The novel ONT approach described in this study can take the place of traditional combination genome assembly based on short- and long-read sequencing, enabling pangenomic analyses based on high-quality complete bacterial and plasmid genomes to monitor the spread of antibiotic-resistant bacteria and antibiotic resistance genes.

Metagenomic sequencing, molecular characterization, and Bayesian phylogenetics of imported type 2 vaccine-derived poliovirus, Spain, 2021.

In 2021, a type 2 vaccine-derived poliovirus (VDPV2) was isolated from the stool of a patient with acute flaccid paralysis (AFP) admitted to Spain from Senegal. A virological investigation was conducted to characterize and trace the origin of VDPV2. We used an unbiased metagenomic approach for the whole-genome sequencing of VDPV2 from the stool (pre-treated with chloroform) and from the poliovirus-positive supernatant. Phylogenetic analyses and molecular epidemiological analyses relying on the Bayesian Markov Chain Monte Carlo methodology were used to determine the geographical origin and estimate the date of the initiating dose of the oral poliovirus vaccine for the imported VDPV2. We obtained a high percentage of viral reads per total reads mapped to the poliovirus genome (69.5% for pre-treated stool and 75.8% for isolate) with a great depth of sequencing coverage (5,931 and 11,581, respectively) and complete genome coverage (100%). The two key attenuating mutations in the Sabin 2 strain had reverted (A481G in the 5'UTR and Ile143Thr in VP1). In addition, the genome had a recombinant structure between type-2 poliovirus and an unidentified non-polio enterovirus-C (NPEV-C) strain with a crossover point in the protease-2A genomic region. VP1 phylogenetic analysis revealed that this strain is closely related to VDPV2 strains circulating in Senegal in 2021. According to Bayesian phylogenetics, the most recent common ancestor of the imported VDPV2 could date back 2.6 years (95% HPD: 1.7-3.7) in Senegal. We suggest that all VDPV2s circulating in 2020-21 in Senegal, Guinea, Gambia, and Mauritania have an ancestral origin in Senegal estimated around 2015. All 50 stool samples from healthy case contacts collected in Spain (n = 25) and Senegal (n = 25) and four wastewater samples collected in Spain were poliovirus negative. By using a whole-genome sequencing protocol with unbiased metagenomics from the clinical sample and viral isolate with high sequence coverage, efficiency, and throughput, we confirmed the classification of VDPV as a circulating type. The close genomic linkage with strains from Senegal was consistent with their classification as imported. Given the scarce number of complete genome sequences for NPEV-C in public databases, this protocol could help expand poliovirus and NPEV-C sequencing capacity worldwide.

Noninvasive genomic profiling of somatic mutations in oral cavity cancers

ObjectivesSomatic mutations may predict prognosis, therapeutic response, or cancer progression. We evaluated targeted sequencing of oral rinse samples (ORS) for non-invasive mutational profiling of oral squamous cell carcinomas (OSCC). Materials and methodsA custom hybrid capture panel targeting 42 frequently mutated genes in OSCC was used to identify DNA sequence variants in matched ORS and fresh-frozen tumors from 120 newly-diagnosed patients. Receiver operating characteristic (ROC) curves determined the optimal variant allele fraction (VAF) cutoff for variant discrimination in ORS. Behavioral, clinical, and analytical factors were evaluated for impacts on assay performance. ResultsHalf of tumors involved oral tongue (50 %), and a majority were T1-T2 tumor stage (55 %). Median depth of sequencing coverage was 260X for OSCC and 1,563X for ORS. Frequencies of single nucleotide variants (SNVs) at highly mutated genes (including TP53, FAT1, HRAS, NOTCH1, CDKN2A, CASP8, NFE2L2, and PIK3CA) in OSCC were highly correlated with TCGA data (R = 0.96, p = 2.5E−22). An ROC curve with area-under-the-curve (AUC) of 0.80 showed that, at an optimal VAF cutoff of 0.10 %, ORS provided 76 % sensitivity, 96 % specificity, but precision of only 2.6E−4. At this VAF cutoff, 206 of 270 SNVs in OSCC were detected in matched ORS. Sensitivity varied by patient, T stage and target gene. Neither downsampled ORS as matched control nor a naïve Bayesian classifier adjusting for sequencing bias appreciably improved assay performance. ConclusionTargeted sequencing of ORS provides moderate assay performance for noninvasive detection of SNVs in OSCC. Our findings strongly rationalize further clinical and laboratory optimization of this assay, including strategies to improve precision.

XY sex determination in a cnidarian

BackgroundSex determination occurs across animal species, but most of our knowledge about its mechanisms comes from only a handful of bilaterian taxa. This limits our ability to infer the evolutionary history of sex determination within animals.ResultsIn this study, we generated a linkage map of the genome of the colonial cnidarian Hydractinia symbiolongicarpus and used it to demonstrate that this species has an XX/XY sex determination system. We demonstrate that the X and Y chromosomes have pseudoautosomal and non-recombining regions. We then use the linkage map and a method based on the depth of sequencing coverage to identify genes encoded in the non-recombining region and show that many of them have male gonad-specific expression. In addition, we demonstrate that recombination rates are enhanced in the female genome and that the haploid chromosome number in Hydractinia is n = 15.ConclusionsThese findings establish Hydractinia as a tractable non-bilaterian model system for the study of sex determination and the evolution of sex chromosomes.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans.

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks.

139. Feasibility of comprehensive whole genome profiling in hematological malignancies

Whole genome sequencing (WGS) has been recently proposed as a replacement for routine karyotype and FISH analyses of myeloid and lymphoid neoplasms. WGS allows for the unbiased and comprehensive detection of multiple types of genomic aberrations, from single nucleotide to structural variants, which could streamline current genomic testing algorithms and clinical laboratory workflows for hematologic malignancies. In this study, we explored the feasibility of whole genome profiling of myeloid and lymphoid neoplasms by sequencing eleven AML, MPN, MDS, B-ALL, and T-ALL genomes with an average depth of sequencing coverage of ∼40X, and comparing genome analyses data to their available cytogenetic and molecular panel results. 81% of clinically relevant SNVs were detected by WGS, with false negatives being mostly variants at low frequencies (VAF <10%). Nine new variants were detected by WGS in an AML, B-ALL and T-ALL samples, eight of which were classified as VUS, and one (FLT3 p.Tyr572Cys) in B-ALL as likely pathogenic. 95% of clinically relevant copy-number variants (CNVs) were detected by WGS, with false negatives being mostly small (<150Kb) or mosaic CNVs present in <15% of the sample. The identification of the clinically relevant translocations in these samples (such as ETV6/RUNX1 and CRLF2 rearrangements) using WGS data is currently being investigated. This feasibility study demonstrated the potential of using WGS as the replacement of molecular panels and chromosome microarray in hematologic malignancies, particularly at the time of initial diagnosis.

The population genomic legacy of the second plague pandemic

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

A high scale SARS-CoV-2 profiling by its whole-genome sequencing using Oxford Nanopore Technology in Kazakhstan.

Severe acute respiratory syndrome (SARS-CoV-2) is responsible for the worldwide pandemic, COVID-19. The original viral whole-genome was sequenced by a high-throughput sequencing approach from the samples obtained from Wuhan, China. Real-time gene sequencing is the main parameter to manage viral outbreaks because it expands our understanding of virus proliferation, spread, and evolution. Whole-genome sequencing is critical for SARS-CoV-2 variant surveillance, the development of new vaccines and boosters, and the representation of epidemiological situations in the country. A significant increase in the number of COVID-19 cases confirmed in August 2021 in Kazakhstan facilitated a need to establish an effective and proficient system for further study of SARS-CoV-2 genetic variants and the development of future Kazakhstan’s genomic surveillance program. The SARS-CoV-2 whole-genome was sequenced according to SARS-CoV-2 ARTIC protocol (EXP-MRT001) by Oxford Nanopore Technologies at the National Laboratory Astana, Kazakhstan to track viral variants circulating in the country. The 500 samples kindly provided by the Republican Diagnostic Center (UMC-NU) and private laboratory KDL “Olymp” were collected from individuals in Nur-Sultan city diagnosed with COVID-19 from August 2021 to May 2022 using real-time reverse transcription-quantitative polymerase chain reaction (RT-qPCR). All samples had a cycle threshold (Ct) value below 20 with an average Ct value of 17.03. The overall average value of sequencing depth coverage for samples is 244X. 341 whole-genome sequences that passed quality control were deposited in the Global initiative on sharing all influenza data (GISAID). The BA.1.1 (n = 189), BA.1 (n = 15), BA.2 (n = 3), BA.1.15 (n = 1), BA.1.17.2 (n = 1) omicron lineages, AY.122 (n = 119), B.1.617.2 (n = 8), AY.111 (n = 2), AY.126 (n = 1), AY.4 (n = 1) delta lineages, one sample B.1.1.7 (n = 1) belongs to alpha lineage, and one sample B.1.637 (n = 1) belongs to small sublineage were detected in this study. This is the first study of SARS-CoV-2 whole-genome sequencing by the ONT approach in Kazakhstan, which can be expanded for the investigation of other emerging viral or bacterial infections on the country level.

Discovery of Novel Variants on the CHD7 Gene: A Case Series of CHARGE Syndrome.

Background: CHARGE syndrome (CS) is a single-gene genetic disorder with multiple organ malformations caused by a variant of the chromodomain helicase DNA-binding protein 7 (CHD7) gene on chromosome 8q12.1. In this study, we aimed to investigate new variants that have emerged in these cases compared with typical CS and the relationship between the genes and phenotypes. Methods: Patients with suspected genetic diseases were subjected to Whole Exome Sequencing (WES) at a genetics laboratory in Guangzhou. The average sequencing coverage depth was >200 ×, and 96% was >20 ×. The variant interpretation was manipulated according to the American College of Medical Genetics (ACMG) guidelines. Molecular data on databases for ClinVar and CHD7 were also collected and collated. We reviewed the currently described CHD7 variants and analyzed the genetic variation and phenotypic heterogeneity. Results: Data of 12 patients with CS from four hospitals in China were collected. According to gestational age, most of them (8/12) were near-term babies with a lower birth weight than their peers, averaging 2.62 kg. In this study, the most common phenotypes were respiratory tract malformations (11/12), heart malformations (10/12), and central nervous system malformations (9/12). Two fetuses were confirmed to have brain or heart abnormalities during prenatal testing, while 10/12 were found to have abnormalities during prenatal testing. The maximum Acute Physiology and Chronic Health Evaluation (APACHE II) score at admission was 19, and the average was 11.58. Five variants in the CHD7 gene c.7012C > T (p.Q2338*), c.7868delC (p.P2623Rfs*16), c.5405-3C > G, c.6936 + 2T > C, and c.8077-2A > G) were novel and were located in exons 33, 36, and introns 25, 32, and 37, respectively. There may be a positive correlation between exon location and phenotype. Conclusion: Five novel variants were discovered. These expanded the mutational spectrum of the CHD7 gene and the phenotype of CS. There may be a correlation between the new mutation sites and the phenotype, which has some reference value for the evaluation of mutation sites.

Abstract LB130: High-throughput RNA sequencing directly from cell lysates enables reproducible phenotypic profiling for assessing novel compounds and treatments in cancer research

Abstract New oncology therapies are driven by multiple approaches, with a heavy emphasis on small-molecule therapeutics to target cancer cells and provide a therapeutic benefit to cancer patients. These small-molecule compounds can be used as a stand-alone treatment approach or may be combined with other compounds to bring enhanced effects. Compound therapeutic discovery relies heavily on high-throughput screening assays for measurement of phenotypic responses following compound treatment of oncology cell cultures which can range from very targeted assays to wide-ranging unbiased approaches. Standard whole transcriptome RNA-seq is a preferred method for unbiased measurement of phenotypic responses but is often overlooked due to limited scalability and high sample screening costs. To address these challenges, we have developed a high-throughput gene expression (HT-GEx) assay that combines a variety of strategies to enable a low-cost and high-throughput alternative to standard RNA-seq. First, a simplified workflow removes upstream RNA isolation steps and tags transcripts directly from cell lysate. This is achieved by incorporating both a sample barcode and a unique molecular index (UMI) during the reverse transcription reaction. Next, we leverage a 3’ end counting approach to enable a reduction in sequencing depth coverage (i.e. sample cost) without compromising gene detection sensitivity. With these key advances, a high-plex and high-throughput screening assay is achievable at a reduced cost by removing the need to purify RNA, tagging transcripts early in the workflow to allow pooling of samples and reducing sequencing depth. We have combined these high-throughput and cost-saving strategies and present results which confirm HT-GEx offers results on par with standard RNA-seq. Gene detection sensitivity based on number of genes detected per millions of reads is highly similar, with gene detection sensitivity saturated at approximately 2 million reads. Additionally, the number of genes detected between replicates for RNA samples and lysate samples demonstrate strong linear correlation, implying highly reproducible results. Thus, high-throughput gene expression is an ideal method for phenotypic screening in oncology cell lines following stand-alone or combined compound treatments. Citation Format: Andrea O'Hara, Michael Stephens, Ilaria DeVito, Laure Turner, Christopher Mozdzierz, Haythem Latif, Ginger Zhou. High-throughput RNA sequencing directly from cell lysates enables reproducible phenotypic profiling for assessing novel compounds and treatments in cancer research [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr LB130.

Identification of SARS-CoV-2 variants using viral sequencing for the Centers for Disease Control and Prevention genomic surveillance program

BackgroundThe Centers for Disease Control and Prevention contracted with laboratories to sequence the SARS-CoV-2 genome from positive samples across the United States to enable public health officials to investigate the impact of variants on disease severity as well as the effectiveness of vaccines and treatment. Herein we present the initial results correlating RT-PCR quality control metrics with sample collection and sequencing methods from full SARS-CoV-2 viral genomic sequencing of 24,441 positive patient samples between April and June 2021.MethodsRT-PCR confirmed (N Gene Ct value < 30) positive patient samples, with nucleic acid extracted from saliva, nasopharyngeal and oropharyngeal swabs were selected for viral whole genome SARS-CoV-2 sequencing. Sequencing was performed using Illumina COVIDSeq™ protocol on either the NextSeq550 or NovaSeq6000 systems. Informatic variant calling, and lineage analysis were performed using DRAGEN COVID Lineage applications on Illumina’s Basespace cloud analytical system. All sequence data and variant calls were uploaded to NCBI and GISAID.ResultsAn association was observed between higher sequencing coverage, quality, and samples with a lower Ct value, with < 27 being optimal, across both sequencing platforms and sample collection methods. Both nasopharyngeal swabs and saliva samples were found to be optimal samples of choice for SARS-CoV-2 surveillance sequencing studies, both in terms of strain identification and sequencing depth of coverage, with NovaSeq 6000 providing higher coverage than the NextSeq 550. The most frequent variants identified were the B.1.617.2 Delta (India) and P.1 Gamma (Brazil) variants in the samples sequenced between April 2021 and June 2021. At the time of submission, the most common variant > 99% of positives sequenced was Omicron.ConclusionThese initial analyses highlight the importance of sequencing platform, sample collection methods, and RT-PCR Ct values in guiding surveillance efforts. These surveillance studies evaluating genetic changes of SARS-CoV-2 have been identified as critical by the CDC that can affect many aspects of public health including transmission, disease severity, diagnostics, therapeutics, and vaccines.

Ultra-deep sequencing data from a liquid biopsy proficiency study demonstrating analytic validity

Recently we reported the accuracy and reproducibility of circulating tumor DNA (ctDNA) assays using a unique set of reference materials, associated analytical framework, and suggested best practices. With the rapid adoption of ctDNA sequencing in precision oncology, it is critical to understand the analytical validity and technical limitations of this cutting-edge and medical-practice-changing technology. The SEQC2 Oncopanel Sequencing Working Group has developed a multi-site, cross-platform study design for evaluating the analytical performance of five industry-leading ctDNA assays. The study used tailor-made reference samples at various levels of input material to assess ctDNA sequencing across 12 participating clinical and research facilities. The generated dataset encompasses multiple key variables, including a broad range of mutation frequencies, sequencing coverage depth, DNA input quantity, etc. It is the most comprehensive public-facing dataset of its kind and provides valuable insights into ultra-deep ctDNA sequencing technology. Eventually the clinical utility of ctDNA assays is required and our proficiency study and corresponding dataset are needed steps towards this goal.

Chromosome-level genome of the globe skimmer dragonfly (Pantala flavescens).

ABSTRACTBackgroundThe globe skimmer dragonfly (Pantala flavescens) is a notable Odonata insect distributed in nature fields and farmlands worldwide, and it is commonly recognized as a natural enemy because it preys on agricultural pests and health pests. As one of the sister groups of winged insects, odonatan species are key to understanding the evolution of insect wings.FindingsWe present a high-quality reference genome of P. flavescens, which is the first chromosome-level genome in the Palaeoptera (Odonata and Ephemeroptera). The assembled genome size was 662 Mb, with a contig N50 of 16.2 Mb. Via Hi-C scaffolding, 648 Mb (97.9%) of contig sequences were clustered, ordered, and assembled into 12 large scaffolds, each corresponding to a natural chromosome. The X chromosome was identified by sequence coverage depth. The repetitive sequences and gene density of the X chromosome are similar to those of autosomal sequences, but the X chromosome shows a much lower degree of heterozygosity. Our analysis shows that the effective population size experienced 3 declining events, which may have been caused by climate change and environmental pollution.ConclusionsThe genome of P. flavescens provides more information on the biology and evolution of insects and will help for the use of this species in pest control.

Evaluation of Germline Structural Variant Calling Methods for Nanopore Sequencing Data.

Structural variants (SVs) are genomic rearrangements that involve at least 50 nucleotides and are known to have a serious impact on human health. While prior short-read sequencing technologies have often proved inadequate for a comprehensive assessment of structural variation, more recent long reads from Oxford Nanopore Technologies have already been proven invaluable for the discovery of large SVs and hold the potential to facilitate the resolution of the full SV spectrum. With many long-read sequencing studies to follow, it is crucial to assess factors affecting current SV calling pipelines for nanopore sequencing data. In this brief research report, we evaluate and compare the performances of five long-read SV callers across four long-read aligners using both real and synthetic nanopore datasets. In particular, we focus on the effects of read alignment, sequencing coverage, and variant allele depth on the detection and genotyping of SVs of different types and size ranges and provide insights into precision and recall of SV callsets generated by integrating the various long-read aligners and SV callers. The computational pipeline we propose is publicly available at https://github.com/davidebolo1993/EViNCe and can be adjusted to further evaluate future nanopore sequencing datasets.