Short-read Whole-genome Sequencing Research Articles

Cas9-targeted-based long-read sequencing for genetic screening of RPE65 locus.

Long-read sequencing (LRS) enables accurate structural variant detection and variant phasing. When a molecular diagnosis is suspected, target enrichment can reduce the cost and duration of sequencing. LRS was conducted in five inherited retinal dystrophy (IRD) patients harboring a monoallelic variant in RPE65 that remained uncharacterized after clinical exome sequencing (CES). CRISPR-Cas9 guide RNA probes were designed to target a 31kb region, including the entire RPE65 locus. The DNA was sequenced on a MinION platform. Short-read ×30 whole-genome sequencing (WGS) was performed for five patients to validate nanopore results. The nanopore sequencing process yielded a median of 271 reads within the targeted region, with a mean depth of 109 and a median read size of 8kb. All variants identified by CES have been detected using this approach, and no additional RPE65 gene causative variants were found. Nanopore variant detection demonstrated performance akin to short-read WGS at similar coverage levels, although exhibiting increased false positive calls at lower coverage. In this study, we explore the advantages of using a targeted approach together with long-read sequencing to identify variants associated with IRD. The results underscore the utility of targeted long reads for characterizing patients affected by rare diseases when first-tier diagnostic tests are non-conclusive.

Yield of genetic association signals from genomes, exomes and imputation in the UK Biobank

Whole-genome sequencing (WGS), whole-exome sequencing (WES) and array genotyping with imputation (IMP) are common strategies for assessing genetic variation and its association with medically relevant phenotypes. To date, there has been no systematic empirical assessment of the yield of these approaches when applied to hundreds of thousands of samples to enable the discovery of complex trait genetic signals. Using data for 100 complex traits from 149,195 individuals in the UK Biobank, we systematically compare the relative yield of these strategies in genetic association studies. We find that WGS and WES combined with arrays and imputation (WES + IMP) have the largest association yield. Although WGS results in an approximately fivefold increase in the total number of assayed variants over WES + IMP, the number of detected signals differed by only 1% for both single-variant and gene-based association analyses. Given that WES + IMP typically results in savings of lab and computational time and resources expended per sample, we evaluate the potential benefits of applying WES + IMP to larger samples. When we extend our WES + IMP analyses to 468,169 UK Biobank individuals, we observe an approximately fourfold increase in association signals with the threefold increase in sample size. We conclude that prioritizing WES + IMP and large sample sizes rather than contemporary short-read WGS alternatives will maximize the number of discoveries in genetic association studies.

Easy analysis of bacterial whole-genome sequencing data for clinical microbiologists using open-source Galaxy platform: Characterization of ESBL-producing Enterobacterales from bloodstream infections

ObjectivesClinical microbiologists require easy-to-use open access tools with graphical interfaces to perform bacterial whole-genome sequencing (WGS) in routine practice. This study aimed to build a bioinformatics pipeline on the open-source Galaxy platform, facilitating comprehensive and reproducible analysis of bacterial WGS data in a few steps. We then used it to characterize our local epidemiology of ESBL-producing Enterobacterales isolated from patients with bacteremia. MethodsWe built a bioinformatics pipeline consisting of the following sequential tools: Fastp (input data trimming); FastQC (read quality control); SPAdes (genome assembly); Quast (quality control of genome assembly); Prokka (gene annotation); Staramr (ResFinder database) and ABRicate (CARD database) for antimicrobial resistance (AMR) gene screening and molecular strain typing. Paired-end short read WGS data from all ESBL-producing Enterobacterales strains isolated from patients with bacteremia over one year were analysed. ResultsThe Galaxy platform does not require command line tools. The bioinformatics pipeline was constructed within one hour. It only required uploading fastq files and facilitated systematization of de novo assembly of genomes, MLST typing, and AMR gene screening in one step. Among the 66 ESBL-producing strains analysed, the two most frequent ESBL genes were blaCTX−M-15 (62.1%) and blaCTX−M-27 (13.6%). ConclusionsThe open-access Galaxy platform provides a graphical interface and easy-to-use tools suitable for routine use in clinical microbiology laboratories without bioinformatics specialists. We believe that this platform will facilitate fast and low-cost analysis of bacterial WGS data, especially in resource-limited settings.

Prevalence of mutations associated with artemisinin partial resistance and sulfadoxine–pyrimethamine resistance in 13 regions in Tanzania in 2021: a cross-sectional survey

The emergence of the artemisinin partial resistance (ART-R) mutation in the Plasmodium falciparum kelch13 gene (k13), Arg561His, in Rwanda and the regional presence of polymorphisms affecting sulfadoxine-pyrimethamine have raised concern in neighbouring Tanzania. The goal of this study was to assess the status of antimalarial resistance in Tanzania, with a focus on the border with Rwanda, to understand the distribution of the Arg561His mutation, partner drug resistance, and resistance to chemoprevention drugs. In this cross-sectional survey, capillary dried blood spots were collected from malaria positive asymptomatic individuals in the community and symptomatic individuals in health facilities aged 6months and older, in 13regions of mainland Tanzania from Jan 31to June 26, 2021. Exclusion criteria included residence of the areas other than the target sites, presenting to the health facility for care and treatment of conditions other than malaria, and not providing informed consent. Samples were assessed for antimalarial resistance polymorphisms and genetic relatedness using molecular inversion probes targeting Pfalciparum and short-read whole-genome sequencing. The primary outcome was the prevalence of molecular markers of antimalarial resistance at the region level, as well as at the district level in Kagera, a region in the northwest of the country at the border with Rwanda. 6855 (88·1%) of 7782capillary dried blood spot samples collected were successfully genotyped. The overall prevalence of k13 Arg561His in Kagera was 7·7% (90% CI 6·0-9·4; 50of 649), with the highest prevalence in the districts near the Rwandan border (22·8% [31 of 136] in Karagwe, 14·4% [17 of 118]) in Kyerwa, and 1·4% [two of 144] in Ngara). k13 Arg561His was uncommon in the other regions. Haplotype analysis suggested that some of these parasites are related to isolates collected in Rwanda in 2015, supporting regional spread of Arg561His. However, a novel k13 Arg561His haplotype was observed, potentially indicating a second origin in the region. Other validated k13 resistance polymorphisms (one Arg622Ile and two Ala675Val isolates) were also identified. Aregion of prevalent dihydrofolate reductase Ile164Leu mutation, associated with sulfadoxine-pyrimethamine resistance, was also identified in Kagera (15·2% [12·6-17·8%]; 80of 526). The mutant crt Lys76Thr mutation, associated with chloroquine and amodiaquine resistance, was uncommon, occurring only in 75of 2861genotyped isolates, whereases the wild-type mdr1 Asn86Tyr allele, associated with reduced sensitivity to lumefantrine, was found in 99·7% (3819 of 3830) of samples countrywide. These findings show that the k13 Arg561His mutation is common in northwest Tanzania and that multiple emergences of ART-R, similar as to what was seen in southeast Asia, have occurred. Mutations associated with high levels of sulfadoxine-pyrimethamine resistance are common. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region. Understanding how multiple emergences interact with drivers of regional spread is essential for combating ART-R in Africa. This study was funded by the Bill & Melinda Gates Foundation and the National Institutes of Health.

Prevalence of Transcription Factor 4 Gene Triplet Repeat Expansion Associated with Fuchs’ Endothelial Corneal Dystrophy in the United States and Global Populations

ObjectiveAn intronic CTG triplet repeat expansion in the transcription factor 4 gene (TCF4) gene (CTG18.1) confers significant risk for the development of Fuchs’ endothelial corneal dystrophy (FECD). The objective of this study was to conduct an unbiased survey of the CTG18.1 repeat expansion allele frequencies in a multi-ethnic, population-based cohort from the United States and in global populations. DesignCross-sectional study SubjectsDallas Heart Study (DHS) cohort including 1,599 African Americans (AAs), 1,028 European Americans (EAs), and 458 Latinos; 2,500 individuals from the 1000 Genomes Project (1KGP) sampled from 26 populations across 5 continents. MethodsWe genotyped the CTG18.1 short tandem repeat in DHS using targeted polymerase chain reaction amplification followed by fragment analysis. We also inferred the CTG18.1 repeat genotype based on short-read whole-genome sequencing in 1KGP using the computational tool ExpansionHunter. Main Outcome MeasuresThe prevalence of an expanded CTG18.1 allele with ≥40 repeats was determined in U.S. and global populations. ResultsThe carrier rates of the expanded allele were 3.1%, 8.1%, and 3.3% in AAs, EAs, and Latinos, respectively, in the DHS, and 2.7%, 9.5%, 5.2%, 7.2%, and 5.2% in the African (AFR), European (EUR), East Asian (EAS), South Asian (SAS), and admixed American (AMR) continental populations, respectively, in the 1KGP. The distributions of the CTG18.1 repeat in DHS and in 1KGP are similar. The median repeat length was ∼17 with the interquartile range (IQR) between 12 and 23 in the DHS populations. The median repeat length was ∼19 in all the 1KGP populations with the IQR between 13 and 26. The highest prevalence of the expanded allele carriers ranging from 12.1% to 12.5% was observed in some EUR and AMR subpopulations. The frequency of expanded alleles carriers was absent or low (0 - 1.9%) in subpopulations of East Africa but was present at 6.2% in a Kenyan subpopulation in West Africa. ConclusionsThe TCF4 repeat expansion is most prevalent in people of European ancestry and least in African ancestry, which is consistent with FECD prevalence. The expanded TCF4 CTG18.1 allele is the most common disease-causing short tandem repeat in humans with worldwide implications for corneal disease.

Analysis of genomic copy number variations through whole-genome scan in Yunling cattle.

Yunling cattle is a new breed of beef cattle bred in Yunnan Province, China, which has the advantages of fast growth, excellent meat quality, improved tolerance ability, and important landscape value. Copy number variation (CNV) is a significant source of gene structural variation and plays a crucial role in evolution and phenotypic diversity. Based on the latest reference genome ARS-UCD2.0, this study analyzed the genome-wide distribution of CNVs in Yunling cattle using short-read whole-genome sequencing data (n = 129) and single-molecule long-read sequencing data (n = 1), and a total of 16,507 CNVs were detected. After merging CNVs with overlapping genomic positions, 3,728 CNV regions (CNVRs) were obtained, accounting for 0.61% of the reference genome. The functional analysis indicated significant enrichment of CNVRs in 96 GO terms and 57 KEGG pathways, primarily related to cell adhesion, signal transduction, neuromodulation, and nutritional metabolism. Additionally, 111 CNVRs overlapped with 76 quantitative trait loci (QTLs), including Subcutaneous fat thickness QTL, Longissimus muscle area QTL, and Marbling score QTL. Several CNVR-overlapping genes, including BZW1, AOX1, and LOC100138449, overlap with regions associated with meat color and quality QTLs. Furthermore, Vst analysis showed that PSMB4, ERICH1, SMC2, and PPP4R3A were highly divergent between Yunling and Brahman cattle. In summary, we have constructed the genomic CNV map of Yunling cattle for the first time using whole-genome resequencing. This provides valuable genetic variation resources for the study of the Yunling cattle genome and contributes to the study of economic traits in Yunling cattle.

A Comparison of Structural Variant Calling from Short-Read and Nanopore-Based Whole-Genome Sequencing Using Optical Genome Mapping as a Benchmark.

The identification of structural variants (SVs) in genomic data represents an ongoing challenge because of difficulties in reliable SV calling leading to reduced sensitivity and specificity. We prepared high-quality DNA from 9 parent-child trios, who had previously undergone short-read whole-genome sequencing (Illumina platform) as part of the Genomics England 100,000 Genomes Project. We reanalysed the genomes using both Bionano optical genome mapping (OGM; 8 probands and one trio) and Nanopore long-read sequencing (Oxford Nanopore Technologies [ONT] platform; all samples). To establish a "truth" dataset, we asked whether rare proband SV calls (n = 234) made by the Bionano Access (version 1.6.1)/Solve software (version 3.6.1_11162020) could be verified by individual visualisation using the Integrative Genomics Viewer with either or both of the Illumina and ONT raw sequence. Of these, 222 calls were verified, indicating that Bionano OGM calls have high precision (positive predictive value 95%). We then asked what proportion of the 222 true Bionano SVs had been identified by SV callers in the other two datasets. In the Illumina dataset, sensitivity varied according to variant type, being high for deletions (115/134; 86%) but poor for insertions (13/58; 22%). In the ONT dataset, sensitivity was generally poor using the original Sniffles variant caller (48% overall) but improved substantially with use of Sniffles2 (36/40; 90% and 17/23; 74% for deletions and insertions, respectively). In summary, we show that the precision of OGM is very high. In addition, when applying the Sniffles2 caller, the sensitivity of SV calling using ONT long-read sequence data outperforms Illumina sequencing for most SV types.

RFC1 Repeat Distribution in the Cypriot Population: Study of a Large Cohort of Patients With Undiagnosed Ataxia and Non-Disease Controls.

The intronic biallelic AAGGG expansion in the replication factor C subunit 1 (RFC1) gene was recently associated with a phenotype combining cerebellar ataxia, neuropathy, and vestibular areflexia syndrome, as well as with late-onset ataxia. Following this discovery, studies in multiple populations extended the phenotypic and genotypic spectrum of this locus. Multiple benign and additional pathogenic configurations are currently known. Our main objectives were to study the prevalence of the pathogenic AAGGG expansion in the Cypriot population, to further characterize the RFC1 repeat locus allele distribution, and to search for possible novel repeat configurations. Cypriot undiagnosed patients, in the majority presenting at least with cerebellar ataxia and non-neurologic disease controls, were included in this study. A combination of conventional methods was used, including standard PCR flanking the repeat region, repeat-primed PCR, long-range PCR, and Sanger sequencing. Bioinformatics analysis of already available in-house short-read whole-genome sequencing data was also performed. A large group of undiagnosed patients (n = 194), mainly presenting with pure ataxia or with ataxia accompanied by neuropathy or additional symptoms, as well as a group of non-disease controls (n = 100), were investigated in the current study. Our findings include the diagnosis of 10 patients homozygous for the pathogenic AAGGG expansion and a high percentage of heterozygous AAGGG carriers in both groups. The benign AAAAGn, AAAGGn, and AAGAGn configurations were also identified in our cohorts. We also report and discuss the identification of 2 recently reported novel and possibly benign repeat configurations, AAAGGGn and AAGACn, thus confirming their existence in another distinct population, and we highlight an increased frequency of the AAAGGGn in the patient group, including a single case of homozygosity. Our findings indicate the existence of genetic heterogeneity regarding the RFC1 repeat configurations and that the AAGGG pathogenic expansion is a frequent cause of ataxia in the Cypriot population.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BackgroundThe short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.ResultsTo fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.ConclusionsAligning denovo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

Exploring the resistome, virulome, and mobilome of multidrug-resistant Klebsiella pneumoniae isolates: deciphering the molecular basis of carbapenem resistance

BackgroundKlebsiella pneumoniae, a notorious pathogen for causing nosocomial infections has become a major cause of neonatal septicemia, leading to high morbidity and mortality worldwide. This opportunistic bacterium has become highly resistant to antibiotics due to the widespread acquisition of genes encoding a variety of enzymes such as extended-spectrum beta-lactamases (ESBLs) and carbapenemases. We collected Klebsiella pneumoniae isolates from a local tertiary care hospital from February 2019–February 2021. To gain molecular insight into the resistome, virulome, and genetic environment of significant genes of multidrug-resistant K. pneumoniae isolates, we performed the short-read whole-genome sequencing of 10 K. pneumoniae isolates recovered from adult patients, neonates, and hospital tap water samples.ResultsThe draft genomes of the isolates varied in size, ranging from 5.48 to 5.96 Mbp suggesting the genome plasticity of this pathogen. Various genes conferring resistance to different classes of antibiotics e.g., aminoglycosides, quinolones, sulfonamides, tetracycline, and trimethoprim were identified in all sequenced isolates. The highest resistance was observed towards carbapenems, which has been putatively linked to the presence of both class B and class D carbapenemases, blaNDM, and blaOXA, respectively. Moreover, the biocide resistance gene qacEdelta1 was found in 6/10 of the sequenced strains. The sequenced isolates exhibited a broad range of sequence types and capsular types. The significant antibiotic resistance genes (ARGs) were bracketed by a variety of mobile genetic elements (MGEs). Various spontaneous mutations in genes other than the acquired antibiotic-resistance genes were observed, which play an indirect role in making these bugs resistant to antibiotics. Loss or deficiency of outer membrane porins, combined with ESBL production, played a significant role in carbapenem resistance in our sequenced isolates. Phylogenetic analysis revealed that the study isolates exhibited evolutionary relationships with strains from China, India, and the USA suggesting a shared evolutionary history and potential dissemination of similar genes amongst the isolates of different origins.ConclusionsThis study provides valuable insight into the presence of multiple mechanisms of carbapenem resistance in K. pneumoniae strains including the acquisition of multiple antibiotic-resistance genes through mobile genetic elements. Identification of rich mobilome yielded insightful information regarding the crucial role of insertion sequences, transposons, and integrons in shaping the genome of bacteria for the transmission of various resistance-associated genes. Multi-drug resistant isolates that had the fewest resistance genes exhibited a significant number of mutations. K. pneumoniae isolate from water source displayed comparable antibiotic resistance determinants to clinical isolates and the highest number of virulence-associated genes suggesting the possible interplay of ARGs amongst bacteria from different sources.

A novel framework for human leukocyte antigen (HLA) genotyping using probe capture-based targeted next-generation sequencing and computational analysis

Human leukocyte antigen (HLA) genes play pivotal roles in numerous immunological applications. Given the immense number of polymorphisms, achieving accurate high-throughput HLA typing remains challenging. This study aimed to harness the human pan-genome reference consortium (HPRC) resources as a potential benchmark for HLA reference materials. We meticulously annotated specific four field-resolution alleles for 11 HLA genes (HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3, -DRB4 and -DRB5) from 44 high-quality HPRC personal genome assemblies. For sequencing, we crafted HLA-specific probes and conducted capture-based targeted sequencing of the genomic DNA of the HPRC cohort, ensuring focused and comprehensive coverage of the HLA region of interest. We used publicly available short-read whole-genome sequencing (WGS) data from identical samples to offer a comparative perspective. To decipher the vast amount of sequencing data, we employed seven distinct software tools: OptiType, HLA-VBseq, HISAT genotype, SpecHLA, T1K, QzType, and DRAGEN. Each tool offers unique capabilities and algorithms for HLA genotyping, allowing comprehensive analysis and validation of the results. We then compared these results with benchmarks derived from personal genome assemblies. Our findings present a comprehensive four-field-resolution HLA allele annotation for 44 HPRC samples. Significantly, our innovative targeted next-generation sequencing (NGS) approach for HLA genes showed superior accuracy compared with conventional short-read WGS. An integrated analysis involving QzType, T1K, and DRAGEN was developed, achieving 100% accuracy for all 11 HLA genes. In conclusion, our study highlighted the combination of targeted short-read sequencing and astute computational analysis as a robust approach for HLA genotyping. Furthermore, the HPRC cohort has emerged as a valuable assembly-based reference in this realm.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Landscape of genomic structural variations in Indian population-based cohorts: Deeper insights into their prevalence and clinical relevance

Structural variations (SV) are large (> 50 bp) genomic rearrangements comprising deletions, duplications, insertions, inversions, and translocations. Studying SVs is important because they play active and critical roles in regulating gene expression, determining disease predispositions, and identifying population-specific differences among individuals of diverse ancestries. However, SV discoveries in the Indian population using whole-genome sequencing (WGS) has been limited. In this study, using short-read WGS having an average 42X depth of coverage, we identify and characterize 36,210 SVs from 529 individuals enrolled in population-based cohorts in India. These SVs include 24,574 deletions, 2913 duplications, 8710 insertions, and 13 inversions. 1.26% (456 out of 36,210) of the identified SVs can potentially impact the coding regions of genes. Furthermore, 56 of these SVs are highly intolerant to loss-of-function changes to the mapped genes, and five SVs impacting ADAMTS17, CCDC40, and RHCE are common in our study individuals. Seven rare SVs significantly impact dosage sensitivity of genes known to be associated with various clinical phenotypes. Most of the SVs in our study are rare and heterozygous. This fine-scale SV discovery in the under-represented Indian population provides valuable insights that extends beyond Euro-centric human genetic studies.

Abstract 6247: Comprehensive characterization of cancers with epithelial mesenchymal transition phenotype through long-read sequencing

Abstract Introduction: Pancreatic ductal adenocarcinoma (PDAC) and linitis-plastica type gastric cancer (STAD) are the most representative refractory cancers even in the era of molecularly targeted therapy and immunotherapy. Extensive characterization of these tumors, especially regarding the genomic and epigenomic characteristics of pure cancer cells, has been hampered by the high stromal content in the primary tumor samples. Here we took advantage of cancer cells purified from malignant ascites and established corresponding cell lines to overcome the limitations associated with the analysis and to discover their vulnerabilities. Methods: We developed cell lines from malignant ascitic fluid samples obtained from 60 patients with PDAC and 59 patients with linitis-plastica type STAD. We performed comprehensive multi-omic analyses of these cell lines using long-read whole-genome sequencing (WGS) for complex somatic genomic alterations and for methylation aberration. We also conducted long-read RNA sequencing (iso-seq) for full-length isoform exploration, and chromatin immunoprecipitation coupled with sequencing (ChIP-seq) for enhancer analysis. Complementary short-read whole-genome sequencing and RNA sequencing were also performed. Results: Comprehensive transcriptome analyses identified a cluster with high epithelial-mesenchymal transition (EMT) characteristics, consisting of the majority of the PDAC cell lines and a half of the STAD cell lines. ChIP-seq analysis against H3K27ac revealed the distinct transcriptional circuit that generated EMT cluster-specific super-enhancers governed by SMAD3, independent of cancer type. Most of the transcriptional activity was controlled by the gene methylation. We also identified the differential expression levels of the splicing factors between the EMT-activated and the non-activated groups. Comprehensive isoform analysis by iso-seq revealed that the isoform usage was completely different between the EMT group and the other, especially in genes of the actomyosin pathway which confer characteristic features of the EMT cluster in the cell morphology and cell motility. Using the CHM13 reference genome, long-read WGS analyses revealed the recurrent structural alterations specific to the EMT group such as those involving CDKN2A. We also identified PDAC-specific structural alterations. Conclusion: We identified the distinct molecular characteristics of the EMT-activated group in PDAC and linitis-plastica type STAD which may confer the therapeutic resistance. Our findings provide an interesting insight for the further innovation of the treatment strategy for these intractable cancers. Citation Format: Yosuke Tanaka, Hiroyuki Mano. Comprehensive characterization of cancers with epithelial mesenchymal transition phenotype through long-read 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 6247.

On the distribution and diversity of tissue-specific somatic mutations in honey bee (Apis mellifera) drones

Somatic mutations originate from both exogenous (e.g. UV radiation, chemical agents) and endogenous (e.g., DNA replication, defective DNA repair) sources and can have significant impacts on an animal’s reproductive success. This may be especially true for haploid organisms that are susceptible to any deleterious alleles inherited from their parent and any that arise over their lifetime. Unfortunately, little is known about the rate of somatic mutation accumulation across individuals and tissues of haplodiploid animal populations, the functional processes through which they arise, and their distribution across tissues and the genome. Here, we generated short-read whole-genome sequencing data for four tissues of haploid honey bee males. We paired this with estimates of telomere length and tissue-specific DNA content to address three major questions: is there variance in somatic mutational load across haploid individuals and specific tissues therein, does increased DNA content in a tissue contribute to somatic mutational load, and does telomere length correlate with mutational load? Our results suggest that variance in somatic mutational load is better captured across individuals than across tissues, that tissue-specific DNA content is not associated with somatic mutation load, and that variance in telomere length does not correlate with somatic mutation loads across tissues. To our knowledge, this is the first observational study on somatic mutational load in Apoidea and likely Hymenoptera. It serves as a useful advent for additional studies understanding the processes through which haploids tolerate or repair somatic mutations.

HAT: de novo variant calling for highly accurate short-read and long-read sequencing data.

de novo variants (DNVs) are variants that are present in offspring but not in their parents. DNVs are both important for examining mutation rates as well as in the identification of disease-related variation. While efforts have been made to call DNVs, calling of DNVs is still challenging from parent-child sequenced trio data. We developed Hare And Tortoise (HAT) as an automated DNV detection workflow for highly accurate short-read and long-read sequencing data. Reliable detection of DNVs is important for human genomics and HAT addresses this need. HAT is a computational workflow that begins with aligned read data (i.e. CRAM or BAM) from a parent-child sequenced trio and outputs DNVs. HAT detects high-quality DNVs from Illumina short-read whole-exome sequencing, Illumina short-read whole-genome sequencing, and highly accurate PacBio HiFi long-read whole-genome sequencing data. The quality of these DNVs is high based on a series of quality metrics including number of DNVs per individual, percent of DNVs at CpG sites, and percent of DNVs phased to the paternal chromosome of origin. https://github.com/TNTurnerLab/HAT.

Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus.

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Novel characterization of Alzheimer’s disease genetic loci using long‐read sequencing

AbstractBackgroundAlzheimer’s disease (AD) is up to 60‐80% heritable, but less than ∼20% is explained by studies analyzing single nucleotide variants (SNVs). One limitation of short‐read whole‐genome sequencing (SRS) is the standard read length of 150 base pairs, which does not enable the detection of longer structural variants (SVs). Here we sequenced individuals using long‐read sequencing (LRS), which can sequence reads with an average length of ∼20 kilobases, allowing us to identify SVs previously uncaptured.MethodAll participants underwent whole‐genome LRS (∼15x coverage) and a subgroup (84%) underwent SRS. Out of 576 participants (47% males, age = 70.6±7.9 y.o.), 115 were diagnosed with AD or mild cognitive impairment, 365 were healthy controls, and 96 were diagnosed with a synucleinopathy (either Parkinson or Lewy Body disease). Eighty‐three index SNVs from loci associated with AD risk through GWAS (Bellenguez et al., Nature Genetics 2022) were genotyped with 30x coverage SRS, in addition to APOE2‐4. A 1Mbp window was defined around these SNVs to construct the discovery range for SVs (Sniffles2 population mode). Linkage disequilibrium (LD) was assessed between SNVs and SVs (CubeX).ResultA total of 14854 SVs were found across the AD risk loci (Figure1). After LD calculation, N = 197 SVs had a R2>0.1 (Figure2). The SVs with the highest LD (R2>0.7) are reported in Table1, among which there is a 322 bp deletion in the 3’ UTR region of the TMEM106B in high LD (R2 = 0.918) with the intronic SNV rs13237518. This TMEM106B locus has been previously associated with the risk of frontotemporal lobar dementia with TDP‐43 pathology in addition to AD, but the causative SNV has not been identified yet. This large deletion may mediate TMEM106B’s risk‐modulating role in AD and FTLD‐TDP (Chemparathy et al. medRxiv 2023). At the complex MAPT locus, three SVs showed a high LD with rs199515.ConclusionInsights into the role of SVs in neurodegenerative disorders have been hampered due to limitations with SRS. Using LRS in a large AD‐related sample, we characterized for the first time the genetic variation of SVs in known AD risk loci and provide a roadmap to identify potential causal SVs driving the AD association signal.

Routine Whole Genome Sequencing Offers a More Comprehensive Assessment of IG/TCR Gene Rearrangements for MRD Targets in Acute Lymphoblastic Leukemia - Data from the First 62 Cases from the NHS England Genomic Medicine Service

Introduction Minimal residual disease (MRD) monitoring is fundamental for risk stratification and treatment in acute lymphoblastic leukemia (ALL). Sensitive molecular assay requires presentation samples to be screened for immunoglobulin and T-cell receptor (IG/TCR) gene rearrangements. Previously samples were screened by PCR, heteroduplex analysis and Sanger sequencing to identify clone defining IG/TCR gene rearrangements. Subsequently these sequences were used to develop qPCR assays for MRD follow-up monitoring. More recently, targeted IG/TCR high-throughput sequencing (HTS) has been used with similar panels to identify gene rearrangements. Both these methods have limitations, in that only rearrangements covered by the panels are detected, and the techniques are labour intensive and expensive. Whole genome sequencing (WGS) has potential to discover all IG/TCR rearrangements in a single assay giving greater potential to detect markers for follow up MRD monitoring. In England all patients with ALL are eligible for WGS at diagnosis via the National Health Service, Genomic Medicine Service (NHS GMS). Methods Sixty-two patients, presenting with ALL at either Great Ormond Street Hospital or Bristol Children's Hospital, had short-read WGS via the NHS GMS performed concurrently with standard of care (SOC) targeted HTS to identify IG/TCR rearrangements. A bespoke bioinformatic pipeline was developed to identify IG/TCR rearrangements in the WGS samples. SOC targeted HTS analysis was performed using the open source Vidjil bioinformatic software. A cut-off of 5% (clone proportion of total clonotypes) was used to define an acceptable disease-specific clone for monitoring. Results All 62 samples had a greater number of rearrangements detected by WGS than targeted HTS. The mean number of clones >5% detected was 2.7 by targeted HTS and 3.1 by WGS. For rearrangements identified by targeted HTS at a 5% cut-off, 97% (185 out of 190) were identified by WGS. All targets were identified by WGS using a 10% cut-off (Table 1).Despite the lower depth of coverage of the WGS assay, it retrieved 89% of sequences seen in SOC down to a level of 1%, as well as additional re-arrangements not seen by the SOC assay (Figure 1). As well as additional rearrangements being discovered by WGS, we observed novel rearrangements (dd3/ja rearrangements) not previously identified using a targeted approach. To establish if the additional sequences identified by WGS were clinically applicable as MRD targets, we designed allele specific oligonucleotides for a representative 32 rearrangements for testing by qPCR. Over half (n=16) were able to be used to track MRD in follow-up samples. New primer/probe sets were designed for the novel dd3/ja rearrangements identified by WGS, these were also tested by qPCR and performed well enough to be clinically applicable for MRD monitoring. Conclusions WGS approach outperformed targeted HTS, identifying more IG/TCR rearrangements suitable for molecular MRD monitoring in all patients. This is of particularly benefit in cases where targeted HTS only identifies 1 or 2 potential MRD targets. WGS also has the advantage of identifying other key gene rearrangements (BCR-ABL, KMT2A, SIL-TAL, IKZF1) that make excellent targets for MRD monitoring. Beyond this, WGS delivers all diagnostic information, germline and additional somatic findings that aid diagnosis. For ALL, in one test WGS can now replace all molecular diagnostics. In England all patients with leukemia are eligible for centralised WGS. With continued improvements in turnaround times, and reducing costs, we propose that WGS should be the new diagnostic standard for pediatric ALL

Cost-Effective Cas9-Mediated Targeted Sequencing of Spinocerebellar Ataxia Repeat Expansions

Hereditary repeat diseases are caused by an abnormal expansion of short tandem repeats in the genome. Among them, spinocerebellar ataxia (SCA) is a heterogeneous disease, and currently, 16 responsible repeats are known. Genetic diagnosis is obtained by analyzing the number of repeats through separate testing of each repeat. Although simultaneous detection of candidate repeats using current massively parallel sequencing technologies has been developed to avoid complicated multiple experiments, these methods are generally expensive. This study developed a cost-effective SCA repeat panel [Flongle SCA repeat panel sequencing (FLO-SCAp)] using Cas9-mediated targeted long-read sequencing and the smallest long-read sequencing apparatus, Flongle. This panel enabled the detection of repeat copy number changes, internal repeat sequences, and DNA methylation in seven patients with different repeat expansion diseases. The median (interquartile range) values of coverage and on-target rate were 39.5 (12 to 72) and 11.6% (7.5% to 16.5%), respectively. This approach was validated by comparing repeat copy number changes measured by FLO-SCAp and short-read whole-genome sequencing. A high correlation was observed between FLO-SCAp and short-read whole-genome sequencing when the repeat length was ≤250 bp (r=0.98; P<0.001). Thus, FLO-SCAp represents the most cost-effective method for conducting multiplex testing of repeats and can serve as the first-line diagnostic tool for SCA.