Genome Sequence Data Research Articles

A repertoire of single nucleotide polymorphisms (SNPs) of major fecundity BMPR1B gene among 75 sheep breeds worldwide

The BMPR1B gene is a major determinant of sheep reproductive capacity. Previous studies revealed that Q249R (FecB) is a profound variant of BMPR1B that influences the ovulation rate and litter size in sheep. However, unlike Q249R locus, the full spectrum of single nucleotide polymorphisms (SNPs) within BMPR1B has not been extensively studied. A systematic screen of SNPs in BMPR1B would facilitate the discovery of novel variants that are associated with litter size. This study aimed to investigate SNPs in the BMPR1B gene via whole genome sequence (WGS) data from 2409 individuals of 75 sheep breeds worldwide. Herein, a total of 9688 variants were screened, among which 15 were coding variants and 8 were novel changes. Specifically, we presented the most comprehensive frequency distribution map of the well-known FecB mutation to date. Besides, among the above-mentioned SNPs, one synonymous mutation (g.30050773C > T) was found to be likely under selection and is potentially associated with fecundity in Duolang sheep. Thus, our study greatly expands the variation repertoire of the ovine BMPR1B gene and provides a valuable resource for exploring causative mutations and genetic markers associated with litter size.

Population epigenetics: DNA methylation in the plant omics era.

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

A comprehensive review of deep learning-based variant calling methods.

Genome sequencing data have become increasingly important in the field of personalized medicine and diagnosis. However, accurately detecting genomic variations remains a challenging task. Traditional variation detection methods rely on manual inspection or predefined rules, which can be time-consuming and prone to errors. Consequently, deep learning-based approaches for variation detection have gained attention due to their ability to automatically learn genomic features that distinguish between variants. In our review, we discuss the recent advancements in deep learning-based algorithms for detecting small variations and structural variations in genomic data, as well as their advantages and limitations.

Robust expansion of phylogeny for fast-growing genome sequence data.

Massive sequencing of SARS-CoV-2 genomes has urged novel methods that employ existing phylogenies to add new samples efficiently instead of de novo inference. 'TIPars' was developed for such challenge integrating parsimony analysis with pre-computed ancestral sequences. It took about 21 seconds to insert 100 SARS-CoV-2 genomes into a 100k-taxa reference tree using 1.4 gigabytes. Benchmarking on four datasets, TIPars achieved the highest accuracy for phylogenies of moderately similar sequences. For highly similar and divergent scenarios, fully parsimony-based and likelihood-based phylogenetic placement methods performed the best respectively while TIPars was the second best. TIPars accomplished efficient and accurate expansion of phylogenies of both similar and divergent sequences, which would have broad biological applications beyond SARS-CoV-2. TIPars is accessible from https://tipars.hku.hk/ and source codes are available at https://github.com/id-bioinfo/TIPars.

Detection and genetic diversity of subgroup K avian leukosis virus in local chicken breeds in Jiangxi from 2021 to 2023.

Avian leukosis virus subgroup K (ALV-K) is a new subgroup of avian leukosis virus (ALV) that was first identified in Chinese native chickens in recent years. To further understand the molecular epidemiology and evolutionary diversity of ALV-K, this study investigated the molecular epidemiology of 73,664 chicken plasma samples collected from Jiangxi native chicken flocks. The results showed that ALV-J was the most predominant ALV subtype in Jiangxi native chickens, with a high positivity rate of 4.34%. From 2021 to 2023, there was a gradual upward trend in the proportion of positive numbers of ALV-K among ALV-positive samples, and there was a trend of outbreaks. ALV-J and ALV-K were the main co-infection patterns. Genetic evolutionary analysis based on ALV-K gp85 gene showed that the isolated ALV-K in this study were distributed in various branches of the evolutionary tree with genetic diversity. The homology results showed that the amino acid homology of the isolated ALV-K gp85 gene ranged from 33.9 to 88.1% with the reference strains of subtypes A, B, C, D, E, and J, and from 91.9 to 100% with the other ALV-K reference strains. Multiple mutations were present in the ALV-K gp85, and especially significant mutations were found in the highly variable region hr2. The results of ALV-K replication efficiency showed that the replication efficiency of ALV-K was significantly lower than that of ALV-J. These results enriched the genome sequence data of ALV-K in Chinese geoducks, and laid the foundation for further research on the pathogenesis and prevention of ALV-K.

Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly.

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Abstract PR005: Tracing the clonal dynamic of metastatic castration-resistant prostate cancer over immunotherapy using circulating tumor DNA

Abstract Human cancer tissues are known to release their own DNA into body fluids such as plasma. These DNA molecules derived from tumors, collectively known as circulating tumor DNA (ctDNA), contain genetic information about disease progression at the time of sampling. Importantly, ctDNA has recently been utilized for capturing the whole-body disease profile and when sampled at various time points over the treatment, can inform the dynamic of clonal heterogeneity and its contribution to therapy failures. Metastatic castration-resistant prostate cancer (mCRPC) is a class of cancerous diseases with a high degree of clonal heterogeneity and high mortality which is largely due to the frequent emergence of therapy resistance. Therefore, reconstructing the clonal history for mCRPC diseases using ctDNA has great potential in identification of novel and targetable molecular mechanisms that drive mCRPC progression. To assess this possibility, we hypothesize that tumor clones that are responsible for therapy resistance carry distinct genomic profiles, which are captured in plasma ctDNA and could be computationally resolved through genomic sequencing and clonal reconstruction. To address this question, we obtained plasma cell-free DNA (cfDNA), a mixture of ctDNA and other tissue-derived DNA, from 38 individuals treated with combination PD-L1 and PARP inhibition involved in a recent clinical trial (NCT02484404). Whole-genome sequencing was performed using cfDNA and fragmented buffy coat DNA as germline control. cfDNA samples with a tumor fraction of lower than 10% were excluded. Various computational strategies were used to model the clonal structures of diseases and estimate the temporal order in which small and large structural variants emerged. Through this approach, we observed a negative association between cfDNA tumor fraction and therapy response. Multiple clonal evolutionary patterns leading to therapy failures were observed. For example, both gain and loss of activating androgen receptor mutations were detected in different cases following the treatment, which represents two distinct mechanisms of evading immunotherapy targeting. As a result, changes in the relevant transcriptional activities were also detected through ctDNA fragment-based computational inference. Additionally, clonal persistence was detected in multiple cases that exhibited stable diseases upon the treatment, which involves loss-of-function mutations in p53, ATM, and CDK12. To estimate the potential contribution to plasma ctDNA by metastatic diseases, multi-region whole genome and exome sequencing data of the tumors that are available for some cases prior to the treatment were also analyzed. Clones resolved from analyzing ctDNA were computationally mapped to the primary disease or metastatic tissue with a shared set of mutations in a subset of cases and may explain the patterns of treatment response. Our novel findings from this study will fill in a critical gap of knowledge about the complex genetic mechanisms driving mCRPC immunotherapy resistance. Citation Format: Chennan Li, Anna Baj, Clara C. Y. Seo, Nicholas T. Terrigino, John R. Bright, S. Thomas Hennigan, Isaiah M. King, Scott Wilkinson, Shana Y. Trostel, William D. Figg, William L. Dahut, Jung Min Lee, David Y. Takeda, Fatima Karzai, Adam G. Sowalsky. Tracing the clonal dynamic of metastatic castration-resistant prostate cancer over immunotherapy using circulating tumor DNA [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr PR005.

Clinical and public health implications of increasing notifications of LEE-negative Shiga toxin-producing Escherichia coli in England, 2014-2022.

Introduction. Shiga toxin-producing Escherichia coli (STEC) belong to a diverse group of gastrointestinal pathogens. The pathogenic potential of STEC is enhanced by the presence of the pathogenicity island called the Locus of Enterocyte Effacement (LEE), including the intimin encoding gene eae.Gap statement. STEC serotypes O128:H2 (Clonal Complex [CC]25), O91:H14 (CC33), and O146:H21 (CC442) are consistently in the top five STEC serotypes isolated from patients reporting gastrointestinal symptoms in England. However, they are eae/LEE-negative and perceived to be a low risk to public health, and we know little about their microbiology and epidemiology.Aim. We analysed clinical outcomes and genome sequencing data linked to patients infected with LEE-negative STEC belonging to CC25 (O128:H2, O21:H2), CC33 (O91:H14) and, and CC442 (O146:H21, O174:H21) in England to assess the risk to public health.Results. There was an almost ten-fold increase between 2014 and 2022 in the detection of all STEC belonging to CC25, CC33 and CC442 (2014 n=38, 2022 n=336), and a total of 1417 cases. There was a higher proportion of female cases (55-70 %) and more adults than children, with patients aged between 20-40 and >70 most at risk across the different serotypes. Symptoms were consistent across the three dominant serotypes O91:H14 (CC33), O146:H21 (CC442) and O128:H2 (CC25) (diarrhoea >75 %; bloody diarrhoea 25-32 %; abdominal pain 64-72 %; nausea 37-45 %; vomiting 10-24 %; and fever 27-30 %). Phylogenetic analyses revealed multiple events of acquisition and loss of different stx-encoding prophage. Additional putative virulence genes were identified including iha, agn43 and subA.Conclusions. Continued monitoring and surveillance of LEE-negative STEC infections is essential due to the increasing burden of infectious intestinal disease, and the risk that highly pathogenic strains may emerge following acquisition of the Shiga toxin subtypes associated with the most severe clinical outcomes.

OTSUCNV: an adaptive segmentation and OTSU-based anomaly classification method for CNV detection using NGS data

Copy-number variations (CNVs), which refer to deletions and duplications of chromosomal segments, represent a significant source of variation among individuals, contributing to human evolution and being implicated in various diseases ranging from mental illness and developmental disorders to cancer. Despite the development of several methods for detecting copy number variations based on next-generation sequencing (NGS) data, achieving robust detection performance for CNVs with arbitrary coverage and amplitude remains challenging due to the inherent complexity of sequencing samples. In this paper, we propose an alternative method called OTSUCNV for CNV detection on whole genome sequencing (WGS) data. This method utilizes a newly designed adaptive sequence segmentation algorithm and an OTSU-based CNV prediction algorithm, which does not rely on any distribution assumptions or involve complex outlier factor calculations. As a result, the effective detection of CNVs is achieved with lower computational complexity. The experimental results indicate that the proposed method demonstrates outstanding performance, and hence it may be used as an effective tool for CNV detection.

Molecular determinants of multidrug-resistant tuberculosis in Sierra Leone.

A substantial proportion of MDR-TB strains in Sierra Leone were resistant against all first line drugs; however this makes the all-oral-six-month BPaLM regimen or other 6-9 months all oral regimens still viable, mainly because there was no FQ resistance.Resistance to BDQ was detected, as well as RR, due to mutations outside of the hotspot region. While the prevalence of those resistances was low, it is still cause for concern and needs to be closely monitored.

Exploring the genetic factors behind the discrepancy in resistance to bovine tuberculosis between African zebu cattle and European taurine cattle

Caused by the pathogenic agent Mycobacterium bovis, bovine tuberculosis (bTB) is a major concern in cattle breeding due to both its zoonotic potential and economic impact. Greater resistance to this disease has been reported in certain African zebu breeds compared to European taurine breeds. However the genetic basis for the lower susceptibility to bTB infection observed in zebu cattle remains poorly explored. This study was conducted on whole genome sequencing data of three bTB infection-resistant African zebu breeds and two bTB infection-susceptible taurine breeds to decipher the genetic background. A set of four selection signature statistics based on linkage disequilibrium, site frequency spectrum, and population differentiation were used on SNPs whereas between population variance based VST and t-test were used on CNVs. As a complement, genes from previous literature reported as candidate genes for bTB resistance were also inspected to identify genetic variations. Interestingly, the resulting nine candidate genes had deleterious missense variants (SHC3, IFNGR1, TLR2, TLR6, IL1A, LRRK2, EP300 and IRAK4) or a CNV difference (CD48) segregating between the groups. The genes found in the study play a role in immune pathways activated during Mycobacterium infection, contributing to the proliferation of immune cells and the granuloma formation, ultimately modulating the outcome of the infectious event. In particular, a deleterious variant in the LRRK2 gene, whose deficiency has been linked to improved prognosis upon tuberculosis infection, was found in the bTB infection-resistant zebu breeds. Therefore, these genes constitute credible candidates in explaining the discrepancy in Mycobacterium bovis infection susceptibility among different breed.

Star allele search: a pharmacogenetic annotation database and user-friendly search tool of publicly available 1000 Genomes Project biospecimens

Here we describe a new public pharmacogenetic (PGx) annotation database of a large (n = 3,202) and diverse biospecimen collection of 1000 Genomes Project cell lines and DNAs. The database is searchable with a user friendly, web-based tool (www.coriell.org/StarAllele/Search). This resource leverages existing whole genome sequencing data and PharmVar annotations to characterize *alleles for each biospecimen in the collection. This new tool is designed to facilitate in vitro functional characterization of *allele haplotypes and diplotypes as well as support clinical PGx assay development, validation, and implementation.

Whole-genome sequencing of Fusarium spp. causing sugarcane root rot on both chewing cane and sugar-making cane

Previously we isolated three Fusarium strains (a F. sacchari strain namely GXUF-1, and another two F. commune strains namely GXUF-2 and GXUF-3), and we verified that GXUF-3 was able to cause sugarcane root rot to the chewing cane cultivar Badila. Considering that Fusarium spp. are a group of widely distributed fungal pathogens, we tested whether these three Fusarium isolates were able to cause root rot to Badila as well as sugar-making cane cultivar (Guitang42), using a suitable inoculation method established based on infection assays using Badila. We found that the three Fusarium strains were able to cause root rot symptoms to both Badila and Guitang42, to different extents. To better investigate the potential pathogenicity mechanisms, we performed Illumina high-throughput sequencing and analyzed the whole genomic sequence data of these three Fusarium strains. The results reveal that the assembly sizes of the three Fusarium strains were in a range of 44.7–48.2 Mb, with G + C contents of 48.0–48.5%, and 14,154–15,175 coding genes. The coding genes were annotated by multiple public databases, and potential pathogenic genes were predicted using proprietary databases (such as PHI, DFVF, CAZy, etc.). Furthermore, based on evolutionary analysis of the coding sequence, we found that contraction and expansion of gene families occurred in the three Fusarium strains. Overall, our results suggest a potential risk that the root rot disease may occur to the sugar-making canes although it was initially spotted from fruit cane, and provide clues to understand the pathogenic mechanisms of Fusarium spp. causing sugarcane root rot.

MARS an improved de novo peptide candidate selection method for non-canonical antigen target discovery in cancer

Understanding the nature and extent of non-canonical human leukocyte antigen (HLA) presentation in tumour cells is a priority for target antigen discovery for the development of next generation immunotherapies in cancer. We here employ a de novo mass spectrometric sequencing approach with a refined, MHC-centric analysis strategy to detect non-canonical MHC-associated peptides specific to cancer without any prior knowledge of the target sequence from genomic or RNA sequencing data. Our strategy integrates MHC binding rank, Average local confidence scores, and peptide Retention time prediction for improved de novo candidate Selection; culminating in the machine learning model MARS. We benchmark our model on a large synthetic peptide library dataset and reanalysis of a published dataset of high-quality non-canonical MHC-associated peptide identifications in human cancer. We achieve almost 2-fold improvement for high quality spectral assignments in comparison to de novo sequencing alone with an estimated accuracy of above 85.7% when integrated with a stepwise peptide sequence mapping strategy. Finally, we utilize MARS to detect and validate lncRNA-derived peptides in human cervical tumour resections, demonstrating its suitability to discover novel, immunogenic, non-canonical peptide sequences in primary tumour tissue.

Evaluating the efficacy of MEANGS for mitochondrial genome assembly of cartilaginous and ray-finned fish species.

The assembly of complete and circularized mitochondrial genomes (mitogenomes) is essential for population genetics, phylogenetics and evolution studies. Recently, Song etal. developed a seed-free tool called MEANGS for de novo mitochondrial assembly from whole genome sequencing (WGS) data in animals, achieving highly accurate and intact assemblies. However, the suitability of this tool for marine fish remains unexplored. Additionally, we have concerns regarding the overlap sequences in their original results, which may impact downstream analyses. In this Letter to the Editor, the effectiveness of MEANGS in assembling mitogenomes of cartilaginous and ray-finned fish species was assessed. Moreover, we also discussed the appropriate utilization of MEANGS in mitogenome assembly, including the implementation of the data-cut function and circular detection module. Our observations indicated that with the utilization of these modules, MEANGS efficiently assembled complete and circularized mitogenomes, even when handling large WGS datasets. Therefore, we strongly recommend users employ the data-cut function and circular detection module when using MEANGS, as the former significantly reduces runtime and the latter aids in the removal of overlapped sequences for improved circularization. Furthermore, our findings suggested that approximately 2× coverage of clean WGS data was sufficient for MEANGS to assemble mitogenomes in marine fish species. Moreover, due to its seed-free nature, MEANGS can be deemed one of the most efficient software tools for assembling mitogenomes from animal WGS data, particularly in studies with limited species or genetic background information.

Molecular characterization of lugdunin inactivation mechanisms and their association with Staphylococcus lugdunensis genetic types

Background and purposeOur previous studies showed that lugdunin activities are associated with Staphylococcus lugdunensis genotypes, and most isolates do not exhibit lugdunin activity. As a continuation of our previous analysis, we focused on the reasons for defects in lugdunin production in S. lugdunensis clinical isolates. MethodsA comparative analysis of 36 S. lugdunensis whole genome sequencing data revealed three major mutation types, unknown deletion mechanism that caused most of lug operon genes lost, mobile genetic element (MGE) insertion, and nonsense mutations, which potentially damaged lugdunin production. A total of 152 S. lugdunensis clinical isolates belonging to lugdunin nonproducers were further examined for the above three mutation types. PCR products were sequenced to examine these variations. ResultsForty-six of the 152 isolates were CRISPR-Cas IIC isolates, including 26 ST27, 14 ST4, and 6 ST29 isolates; further investigation confirmed that all of their lug operons had lost almost all lug operon genes except lugM. An IS256 insertion in lugA was identified in 16 isolates, and most isolates (15 over 16) belonged to ST3. In addition, three nonsense mutations caused by single nucleotide substitutions (an adenine deletion in lugB at the 361th and 1219th nucleotides and an adenine deletion in lugC at the 1612nd nucleotide) that were frequently observed among 36 S. lugdunensis whole genome sequencing data were further observed in our clinical isolates. These three nonsense mutations were frequently found in most of CRISPR-Cas IIIA strains, especially in ST6 isolates. ConclusionOur findings suggest that the mechanisms affecting lugdunin production are associated with S. lugdunensis molecular types.

Comprehensive landscape of non-CODIS STRs in global populations provides new insights into challenging DNA profiles

The worldwide implementation of short tandem repeats (STR) profiles in forensic genetics necessitated establishing and expanding the CODIS core loci set to facilitated efficient data management and exchange. Currently, the mainstay CODIS STRs are adopted in most general-purpose forensic kits. However, relying solely on these loci failed to yield satisfactory results for challenging tasks, such as bio-geographical ancestry inference, complex DNA mixture profile interpretation, and distant kinship analysis. In this context, non-CODIS STRs are potent supplements to enhance the systematic discriminating power, particularly when combined with the high-throughput next-generation sequencing (NGS) technique. Nevertheless, comprehensive evaluation on non-CODIS STRs in diverse populations was scarce, hindering their further application in routine caseworks. To address this gap, we investigated genetic variations of 178 historically available non-CODIS STRs from ethnolinguistically different worldwide populations and studied their characteristics and forensic potentials via high-coverage whole genome sequencing (WGS) data. Initially, we delineated the genomic properties of these non-CODIS markers through sequence searching, repeat structure scanning, and manual inspection. Subsequent population genetics analysis suggested that these non-CODIS STRs had comparable polymorphism levels and forensic utility to CODIS STRs. Furthermore, we constructed a theoretical next-generation sequencing (NGS) panel comprising 108 STRs (20 CODIS STRs and 88 non-CODIS STRs), and evaluated its performance in inferring bio-geographical ancestry origins, deconvoluting complex DNA mixtures, and differentiating distant kinships using real and simulated datasets. Our findings demonstrated that incorporating supplementary non-CODIS STRs enabled the extrapolation of multidimensional information from a single STR profile, thereby facilitating the analysis of challenging forensic tasks. In conclusion, this study presents an extensive genomic landscape of forensic non-CODIS STRs among global populations, and emphasized the imperative inclusion of additional polymorphic non-CODIS STRs in future NGS-based forensic systems.

Non-canonical antigens are the largest fraction of peptides presented by MHC class I in mismatch repair deficient murine colorectal cancer

BackgroundImmunotherapy based on checkpoint inhibitors is highly effective in mismatch repair deficient (MMRd) colorectal cancer (CRC). These tumors carry a high number of mutations, which are predicted to translate into a wide array of neoepitopes; however, a systematic classification of the neoantigen repertoire in MMRd CRC is lacking. Mass spectrometry peptidomics has demonstrated the existence of MHC class I associated peptides (MAPs) originating from non-coding DNA regions. Based on these premises we investigated DNA genomic regions responsible for generating MMRd-induced peptides.MethodsWe exploited mouse CRC models in which the MMR gene Mlh1 was genetically inactivated. Isogenic cell lines CT26 Mlh1+/+ and Mlh1-/- were inoculated in immunocompromised and immunocompetent mice. Whole genome and RNA sequencing data were generated from samples obtained before and after injection in murine hosts. First, peptide databases were built from transcriptomes of isogenic cell lines. We then compiled a database of peptides lost after tumor cells injection in immunocompetent mice, likely due to immune editing. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) and matched next-generation sequencing databases were employed to identify the DNA regions from which the immune-targeted MAPs originated. Finally, we adopted in vitro T cell assays to verify whether MAP-specific T cells were part of the in vivo immune response against Mlh1-/- cells.ResultsWhole genome sequencing analyses revealed an unbalanced distribution of immune edited alterations across the genome in Mlh1-/- cells grown in immunocompetent mice. Specifically, untranslated (UTR) and coding regions exhibited the largest fraction of mutations leading to highly immunogenic peptides. Moreover, the integrated computational and LC-MS/MS analyses revealed that MAPs originate mainly from atypical translational events in both Mlh1+/+ and Mlh1-/- tumor cells. In addition, mutated MAPs—derived from UTRs and out-of-frame translation of coding regions—were highly enriched in Mlh1-/- cells. The MAPs trigger T-cell activation in mice primed with Mlh1-/- cells.ConclusionsOur results suggest that—in comparison to MMR proficient CRC—MMRd tumors generate a significantly higher number of non-canonical mutated peptides able to elicit T cell responses. These results reveal the importance of evaluating the diversity of neoepitope repertoire in MMRd tumors.

Kuura-An automated workflow for analyzing WES and WGS data.

The advent of high-throughput sequencing technologies has revolutionized the field of genomic sciences by cutting down the cost and time associated with standard sequencing methods. This advancement has not only provided the research community with an abundance of data but has also presented the challenge of analyzing it. The paramount challenge in analyzing the copious amount of data is in using the optimal resources in terms of available tools. To address this research gap, we propose "Kuura-An automated workflow for analyzing WES and WGS data", which is optimized for both whole exome and whole genome sequencing data. This workflow is based on the nextflow pipeline scripting language and uses docker to manage and deploy the workflow. The workflow consists of four analysis stages-quality control, mapping to reference genome & quality score recalibration, variant calling & variant recalibration and variant consensus & annotation. An important feature of the DNA-seq workflow is that it uses the combination of multiple variant callers (GATK Haplotypecaller, DeepVariant, VarScan2, Freebayes and Strelka2), generating a list of high-confidence variants in a consensus call file. The workflow is flexible as it integrates the fragmented tools and can be easily extended by adding or updating tools or amending the parameters list. The use of a single parameters file enhances reproducibility of the results. The ease of deployment and usage of the workflow further increases computational reproducibility providing researchers with a standardized tool for the variant calling step in different projects. The source code, instructions for installation and use of the tool are publicly available at our github repository https://github.com/dhanaprakashj/kuura_pipeline.

Abstract B114: Mutational "lineage traces" of an acinar cell-of-origin in human pancreatic ductal adenocarcinoma

Abstract A clinical challenge of pancreatic ductal adenocarcinoma (PDAC) is that most cases are detected late, after intractable metastasis is well underway. This also poses a challenge to our basic understanding of PDAC development, as the early steps of neoplasia are separated by a gulf of years or decades from accessible biospecimens and cell lines. Although mouse genetic models of PDAC have helped to fill in this gap, key questions remain incompletely answered. For example, what is the cell of origin of human PDAC – and does the answer to this question actually matter at a clinical level? In mouse, PDAC can be induced from both acinar and duct cells, with the former apparently more susceptible to transformation while the latter appears harder to transform, but capable of generating more aggressive tumors. Meanwhile, human PDAC has been divided into two major subtypes based on genomic data, classical and basal-like, the latter of which appears to be more aggressive. Is there a relationship between cell-of-origin, tumor subtype and clinical outcomes? This question requires the ability to directly infer cell-of-origin in human tumors, an approach to which we report here. Studies in human non-pancreatic cancers have identified somatic “passenger” (non-selected) mutations in several highly-expressed genes of the normal organ – for example, mutations in the ALB (albumin) gene in hepatocellular carcinoma, or SFTP genes (surfactant proteins) in lung adenocarcinoma. These might indicate a mutational process occurring in normal cells, prior to transformation, and constituting a “lineage mark” for the cell of origin – however, as these genes are also highly-transcribed in the cancer cells themselves, their mutations may arise after tumorigenesis is underway. This problem does not apply to acinar genes and PDAC: we and others have established that an initiating event in Kras-induced transformation of acinar cells is complete downregulation of acinar-specific gene products. Therefore, if a mechanism of transcription-induced mutation applies in the pancreas, we would expect it to leave a trace in acinar-derived PDAC. We have completed a large-scale analysis (n>2,000) of human whole genome sequencing data, across several carcinoma types including pancreas and liver, focused on unbiased detection of tissue-specific mutations similar to those previously found in the ALB and SFTP genes. Remarkably, we find robust enrichment of mutations in acinar-specific genes in PDAC, affecting approximately 25% of tumors. This constitutes the first in vivo genetic evidence for an acinar cell-of-origin in human pancreatic cancer, and opens the door to comparative genomic, transcriptomic and outcomes analyses of acinar-derived tumors. These are underway, and the results will be described in our presentation. Citation Format: Charles Murtaugh. Mutational "lineage traces" of an acinar cell-of-origin in human pancreatic ductal adenocarcinoma [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Pancreatic Cancer; 2023 Sep 27-30; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(2 Suppl):Abstract nr B114.