Molecular characterization of QTL for grain zinc and iron concentrations in wheat landrace Chinese Spring.

KEY MESSAGE: Three stable QTL for grain zinc concentration were identified in wheat landrace Chinese Spring. Favorable alleles were more frequent in landraces than in modern wheat cultivars. Wheat is a major source of dietary energy for the growing world population. Developing cultivars with enriched zinc and iron can potentially alleviate human micronutrient deficiency. In this study, a recombinant inbred line (RIL) population with 245 lines derived from cross Zhou 8425B/Chinese Spring was used to detect quantitative trait loci (QTL) for grain zinc concentration (GZnC) and grain iron concentration (GFeC) across four environments. Three stable QTL for GZnC with all favorable alleles from Chinese Spring were identified on chromosomes 3BL, 5AL, and 5BL. These QTL explaining maxima of 8.7%, 5.8%, and 7.1% of phenotypic variances were validated in 125 resequenced wheat accessions encompassing both landraces and modern cultivars using six kompetitive allele specific PCR (KASP) assays. The frequencies of favorable alleles for QGZnCzc.caas-3BL, QGZnCzc.caas-5AL and QGZnCzc.caas-5BL were higher in landraces (90.4%, 68.0%, and 100.0%, respectively) compared to modern cultivars (45.9%, 35.4%, and 40.9%), suggesting they were not selected in breeding programs. Candidate gene association studies on GZnC in the cultivar panel further delimited the QTL into 8.5 Mb, 4.1 Mb, and 47.8 Mb regions containing 46, 4, and 199 candidate genes, respectively. The 5BL QTL located in a region where recombination was suppressed. Two stable and three less stable QTL for GFeC with favorable alleles also from Chinese Spring were identified on chromosomes 4BS (Rht-B1a), 4DS (Rht-D1a), 1DS, 3AS, and 6DS. This study sheds light on the genetic basis of GZnC and GFeC in Chinese Spring and provides useful molecular markers for wheat biofortification.

Construction of high-density genetic map based on SLAF-seq and QTL analysis of major traits in sweetpotato [Ipomoea batatas (L.) Lam.].

Sweetpotato, Ipomoea batatas (L.) Lam., is an important worldwide crop used as feed, food, and fuel. However, its polyploidy, high heterozygosity and self-incompatibility makes it difficult to study its genetics and genomics. Longest vine length (LVL), yield per plant (YPP), dry matter content (DMC), starch content (SC), soluble sugar content (SSC), and carotenoid content (CC) are some of the major agronomic traits being used to evaluate sweetpotato. However limited research has actually examined how these traits are inherited. Therefore, after selecting 212 F1 from a Xin24 × Yushu10 crossing as the mapping population, this study applied specific-locus amplified fragment sequencing (SLAF-seq), at an average sequencing depth of 26.73 × (parents) and 52.25 × (progeny), to detect single nucleotide polymorphisms (SNPs). This approach generated an integrated genetic map of length 2441.56 cM and a mean distance of 0.51 cM between adjacent markers, encompassing 15 linkage groups (LGs). Based on the linkage map, 26 quantitative trait loci (QTLs), comprising six QTLs for LVL, six QTLs for YPP, ten QTLs for DMC, one QTL for SC, one QTL for SSC, and two QTLs for CC, were identified. Each of these QTLs explained 6.3-10% of the phenotypic variation. It is expected that the findings will be of benefit for marker-assisted breeding and gene cloning of sweetpotato.

Optical Genome Mapping as a New Tool to Overcome Conventional Cytogenetics Limitations in Patients with Bone Marrow Failure.

Cytogenetic studies are essential in the diagnosis and follow up of patients with bone marrow failure syndromes (BMFSs), but obtaining good quality results is often challenging due to hypocellularity. Optical Genome Mapping (OGM), a novel technology capable of detecting most types chromosomal structural variants (SVs) at high resolution, is being increasingly used in many settings, including hematologic malignancies. Herein, we compared conventional cytogenetic techniques to OGM in 20 patients with diverse BMFSs. Twenty metaphases for the karyotype were only obtained in three subjects (15%), and no SVs were found in any of the samples. One patient with culture failure showed a gain in chromosome 1q by fluorescence in situ hybridization, which was confirmed by OGM. In contrast, OGM provided good quality results in all subjects, and SVs were detected in 14 of them (70%), mostly corresponding to cryptic submicroscopic alterations not observed by standard techniques. Therefore, OGM emerges as a powerful tool that provides complete and evaluable results in hypocellular BMFSs, reducing multiple tests into a single assay and overcoming some of the main limitations of conventional techniques. Furthermore, in addition to confirming the abnormalities detected by conventional techniques, OGM found new alterations beyond their detection limits.

[Detection of cytogenetic abnormalities in multiple myeloma by using optical genome mapping].

Multiple myeloma (MM) is a plasma cell neoplasm characterized by numerous chromosomal number and structural abnormalities, which are of great significance for risk stratification and response evaluation of MM patients. Optical genome mapping (OGM) is a novel technology that has the potential to resolve many of the issues and limitations associated with traditional cytogenetic methods. To date, the clinical utility of OGM has been validated in the fields of cancer, reproduction, and embryonic dysplasia, et al. In this study, we compared OGM to traditional techniques for the first time in five newly diagnosed MM patients, and evaluated the potential of OGM for detecting cytogenetic aberrations and its clinical application value in MM.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions.

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Characterization of a new greenbug resistance gene Gb9 in a synthetic hexaploid wheat.

Greenbug [Schizaphis graminum (Rondani)] is a serious insect pest that not only damages cereal crops, but also transmits several destructive viruses. The emergence of new greenbug biotypes in the field makes it urgent to identify novel greenbug resistance genes in wheat. CWI 76364 (PI 703397), a synthetic hexaploid wheat (SHW) line, exhibits greenbug resistance. Evaluation of an F2:3 population from cross OK 14319 × CWI 76364 indicated that a dominant gene, designated Gb9, conditions greenbug resistance in CWI 76364. Selective genotyping of a subset of F2 plants with contrasting phenotypes by genotyping-by-sequencing identified 25 SNPs closely linked to Gb9 on chromosome arm 7DL. Ten of these SNPs were converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers for genotyping the entire F2 population. Genetic analysis delimited Gb9 to a 0.6-Mb interval flanked by KASP markers located at 599,835,668 bp (Stars-KASP872) and 600,471,081 bp (Stars-KASP881) on 7DL. Gb9 was 0.5 cM distal to Stars-KASP872 and 0.5 cM proximal to Stars-KASP881. Allelism tests indicated that Gb9 is a new greenbug resistance gene which confers resistance to greenbug biotypes C, E, H, I, and TX1. TX1 is one of the most widely virulent biotypes and has overcome most known wheat greenbug resistance genes. The introgression of Gb9 into locally adapted wheat cultivars is of economic importance, and the KASP markers developed in this study can be used to tag Gb9 in cultivar development.

Combined deletion of MEN1, ATRX and PTEN triggers development of high-grade pancreatic neuroendocrine tumors in mice.

Pancreatic neuroendocrine tumors (PanNETs) are a heterogeneous group of tumors that exhibit an unpredictable and broad spectrum of clinical presentations and biological aggressiveness. Surgical resection is still the only curative therapeutic option for localized PanNET, but the majority of patients are diagnosed at an advanced and metastatic stage with limited therapeutic options. Key factors limiting the development of new therapeutics are the extensive heterogeneity of PanNETs and the lack of appropriate clinically relevant models. In that context, genomic sequencing of human PanNETs revealed recurrent mutations and structural alterations in several tumor suppressors. Here, we demonstrated that combined loss of MEN1, ATRX, and PTEN, tumor suppressors commonly mutated in human PanNETs, triggers the development of high-grade pancreatic neuroendocrine tumors in mice. Histopathological evaluation and gene expression analyses of the developed tumors confirm the presence of PanNET hallmarks and significant overlap in gene expression patterns found in human disease. Thus, we postulate that the presented novel genetically defined mouse model is the first clinically relevant immunocompetent high-grade PanNET mouse model.

High-quality chromosome-level genome assembly of Nicotiana benthamiana.

Nicotiana benthamiana is a fundamental model organism in plant research. Recent advancements in genomic sequencing have revealed significant intraspecific genetic variations. This study addresses the pressing need for a precise genome sequence specific to its geographic origin by presenting a comprehensive genome assembly of the N. benthamiana LAB strain from the Republic of Korea (NbKLAB). We compare this assembly with the widely used NbLAB360 strain, shedding light on essential genomic differences between them. The outcome is a high-quality, chromosome-level genome assembly comprising 19 chromosomes, spanning 2,762 Mb, with an N50 of 142.6 Mb. Comparative analyses revealed notable variations, including 46,215 protein-coding genes, with an impressive 99.5% BUSCO completeness score. Furthermore, the NbKLAB assembly substantially improved the QV from 33% for NbLAB360 to 49%. This refined chromosomal genome assembly for N. benthamiana, in conjunction with comparative insights, provides a valuable resource for genomics research and molecular biology. This accomplishment forms a strong foundation for in-depth exploration into the intricacies of plant genetics and genomics, improved precision, and a comparative framework.

The clinical value of optical genome mapping in the rapid characterization of RB1 duplication and 15q23q24.2 triplication, for more appropriate prenatal genetic counselling.

BACKGROUND: Despite recent advances in prenatal genetic diagnosis, medical geneticists still face considerable difficulty in interpreting the clinical outcome of copy-number-variant duplications and defining the mechanisms underlying the formation of certain chromosomal rearrangements. Optical genome mapping (OGM) is an emerging cytogenomic tool with proved ability to identify the full spectrum of cytogenetic aberrations. METHODS: Here, we report on the use of OGM in a prenatal diagnosis setting. Detailed breakpoint mapping was used to determine the relative orientations of triplicated and duplicated segments in two unrelated foetuses harbouring chromosomal aberrations: a de novo 15q23q24.2 triplication and a paternally inherited 13q14.2 duplication that overlapped partially with the RB1 gene. RESULTS: OGM enabled us to suggest a plausible mechanism for the triplication and confirmed that the RB1 duplication was direct oriented and in tandem. This enabled us to predict the pathogenic consequences, refine the prognosis and adapt the follow-up and familial screening appropriately. CONCLUSION: Along with an increase in diagnostic rates, OGM can rapidly highlight genotype-phenotype correlations, improve genetic counselling and significantly influence prenatal management.

Human telomere length is chromosome end-specific and conserved across individuals.

Short telomeres cause age-related disease, and long telomeres contribute to cancer; however, the mechanisms regulating telomere length are unclear. We developed a nanopore-based method, which we call Telomere Profiling, to determine telomere length at nearly single-nucleotide resolution. Mapping telomere reads to chromosome ends showed chromosome end-specific length distributions that could differ by more than six kilobases. Examination of telomere lengths in 147 individuals revealed that certain chromosome ends were consistently longer or shorter. The same rank order was found in newborn cord blood, suggesting that telomere length is determined at birth and that chromosome end-specific telomere length differences are maintained as telomeres shorten with age. Telomere Profiling makes precision investigation of telomere length widely accessible for laboratory, clinical, and drug discovery efforts and will allow deeper insights into telomere biology.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

A Genome-Wide Analysis of the CEP Gene Family in Cotton and a Functional Study of GhCEP46-D05 in Plant Development.

C-TERMINALLY ENCODED PEPTIDEs (CEPs) are a class of peptide hormones that have been shown in previous studies to play an important role in regulating the development and response to abiotic stress in model plants. However, their role in cotton is not well understood. In this study, we identified 54, 59, 34, and 35 CEP genes from Gossypium hirsutum (2n = 4x = 52, AD1), G. barbadense (AD2), G. arboreum (2n = 2X = 26, A2), and G. raimondii (2n = 2X = 26, D5), respectively. Sequence alignment and phylogenetic analyses indicate that cotton CEP proteins can be categorized into two subgroups based on the differentiation of their CEP domain. Chromosomal distribution and collinearity analyses show that most of the cotton CEP genes are situated in gene clusters, suggesting that segmental duplication may be a critical factor in CEP gene expansion. Expression pattern analyses showed that cotton CEP genes are widely expressed throughout the plant, with some genes exhibiting specific expression patterns. Ectopic expression of GhCEP46-D05 in Arabidopsis led to a significant reduction in both root length and seed size, resulting in a dwarf phenotype. Similarly, overexpression of GhCEP46-D05 in cotton resulted in reduced internode length and plant height. These findings provide a foundation for further investigation into the function of cotton CEP genes and their potential role in cotton breeding.

An R2R3-MYB Transcriptional Factor LuMYB314 Associated with the Loss of Petal Pigmentation in Flax (Linum usitatissimum L.).

The loss of anthocyanin pigments is one of the most common evolutionary transitions in petal color, yet the genetic basis for these changes in flax remains largely unknown. In this study, we used crossing studies, a bulk segregant analysis, genome-wide association studies, a phylogenetic analysis, and transgenic testing to identify genes responsible for the transition from blue to white petals in flax. This study found no correspondence between the petal color and seed color, refuting the conclusion that a locus controlling the seed coat color is associated with the petal color, as reported in previous studies. The locus controlling the petal color was mapped using a BSA-seq analysis based on the F2 population. However, no significantly associated genomic regions were detected. Our genome-wide association study identified a highly significant QTL (BP4.1) on chromosome 4 associated with flax petal color in the natural population. The combination of a local Manhattan plot and an LD heat map identified LuMYB314, an R2R3-MYB transcription factor, as a potential gene responsible for the natural variations in petal color in flax. The overexpression of LuMYB314 in both Arabidopsis thaliana and Nicotiana tabacum resulted in anthocyanin deposition, indicating that LuMYB314 is a credible candidate gene for controlling the petal color in flax. Additionally, our study highlights the limitations of the BSA-seq method in low-linkage genomic regions, while also demonstrating the powerful detection capabilities of GWAS based on high-density genomic variation mapping. This study enhances our genetic insight into petal color variations and has potential breeding value for engineering LuMYB314 to develop colored petals, bast fibers, and seeds for multifunctional use in flax.

Cas9-directed long-read sequencing to resolve optical genome mapping findings in leukemia diagnostics.

Leukemias are genetically heterogeneous and diagnostics therefore includes various standard-of-care (SOC) techniques, including karyotyping, SNP-array and FISH. Optical genome mapping (OGM) may replace these as it detects different types of structural aberrations simultaneously and additionally detects much smaller aberrations (500 bp vs 5-10 Mb with karyotyping). However, its resolution may still be too low to define clinical relevance of aberrations when they are located between two OGM labels or when labels are not distinct enough. Here, we test the potential of Cas9-directed long-read sequencing (LRS) as an additional technique to resolve such potentially relevant new findings. From an internal Bionano implementation study we selected ten OGM calls that could not be validated with SOC methods. Per variant we designed crRNAs for Cas9 enrichment, prepared libraries and sequenced them on a MinION/GridION device. We could confirm all aberrations and, importantly, the actual breakpoints of the OGM calls were located between 0.2 and 5.5 kb of the OGM-estimated breakpoints, confirming the high reliability of OGM. Furthermore, we show examples of redefinition of aberrations between labels that enable judgment of clinical relevance. Our results suggest that Cas9-directed LRS can be a relevant and flexible secondary technique in diagnostic workflows including OGM.

Dynamics of accessible chromatin regions and subgenome dominance in octoploid strawberry.

Subgenome dominance has been reported in diverse allopolyploid species, where genes from one subgenome are preferentially retained and are more highly expressed than those from other subgenome(s). However, the molecular mechanisms responsible for subgenome dominance remain poorly understood. Here, we develop genome-wide map of accessible chromatin regions (ACRs) in cultivated strawberry (2n = 8x = 56, with A, B, C, D subgenomes). Each ACR is identified as an MNase hypersensitive site (MHS). We discover that the dominant subgenome A contains a greater number of total MHSs and MHS per gene than the submissive B/C/D subgenomes. Subgenome A suffers fewer losses of MHS-related DNA sequences and fewer MHS fragmentations caused by insertions of transposable elements. We also discover that genes and MHSs related to stress response have been preferentially retained in subgenome A. We conclude that preservation of genes and their cognate ACRs, especially those related to stress responses, play a major role in the establishment of subgenome dominance in octoploid strawberry.

Genome-Wide and Expression Pattern Analysis of the HIT4 Gene Family Uncovers the Involvement of GHHIT4_4 in Response to Verticillium Wilt in Gossypium hirsutum.

Chromatin remodelers are essential for regulating plant growth, development, and responses to environmental stresses. HIT4 (HEAT-INTOLERANT 4) is a novel stress-induced chromatin remodeling factor that has been less studied in abiotic stress and stress resistance, particularly in cotton. In this study, we conducted a comprehensive analysis of the members of the HIT4 gene family in Gossypium hirsutum using bioinformatics methods, including phylogenetic relationships, gene organization, transcription profiles, phylogenetic connections, selection pressure, and stress response. A total of 18 HIT4 genes were identified in four cotton species, with six HIT4 gene members in upland cotton. Based on the evolutionary relationships shown in the phylogenetic tree, the 18 HIT4 protein sequences were classified into four distinct subgroups. Furthermore, we conducted chromosome mapping to determine the genomic locations of these genes and visually represented the structural characteristics of HIT4 in G. hirsutum. In addition, we predicted the regulatory elements in HIT4 in G. hirsutum and conducted an analysis of repetitive sequences and gene collinearity among HIT4 in four cotton species. Moreover, we calculated the Ka/Ks ratio for homologous genes to assess the selection pressure acting on HIT4. Using RNA-seq, we explored the expression patterns of HIT4 genes in G. hirsutum and Gossypium barbadense. Through weighted gene co-expression network analysis (WGCNA), we found that GHHIT4_4 belonged to the MEblue module, which was mainly enriched in pathways such as DNA replication, phagosome, pentose and glucuronate interconversions, steroid biosynthesis, and starch and sucrose metabolism. This module may regulate the mechanism of upland cotton resistance to Verticillium wilt through DNA replication, phagosome, and various metabolic pathways. In addition, we performed heterologous overexpression of GH_D11G0591 (GHHIT4_4) in tobacco, and the results showed a significant reduction in disease index compared to the wild type, with higher expression levels of disease resistance genes in the transgenic tobacco. After conducting a VIGS (virus-induced gene silencing) experiment in cotton, the results indicated that silencing GHHIT4_4 had a significant impact, the resistance to Verticillium wilt weakened, and the internode length of the plants significantly decreased by 30.7% while the number of true leaves increased by 41.5%. qRT-PCR analysis indicated that GHHIT4_4 mainly enhanced cotton resistance to Verticillium wilt by indirectly regulating the PAL, 4CL, and CHI genes. The subcellular localization results revealed that GHHIT4_4 was predominantly distributed in the mitochondria and nucleus. This study offers preliminary evidence for the involvement of the GHHIT4_4 in cotton resistance to Verticillium wilt and lays the foundation for further research on the disease resistance mechanism of this gene in cotton.

Decryption of the survival "black box": gene family expansion promotes the encystment in ciliated protists.

BACKGROUND: Encystment is an important survival strategy extensively employed by microbial organisms to survive unfavorable conditions. Single-celled ciliated protists (ciliates) are popular model eukaryotes for studying encystment, whereby these cells degenerate their ciliary structures and develop cyst walls, then reverse the process under more favorable conditions. However, to date, the evolutionary basis and mechanism for encystment in ciliates is largely unknown. With the rapid development of high-throughput sequencing technologies, genome sequencing and comparative genomics of ciliates have become effective methods to provide insights into above questions. RESULTS: Here, we profiled the MAC genome of Pseudourostyla cristata, a model hypotrich ciliate for encystment studies. Like other hypotrich MAC genomes, the P. cristata MAC genome is extremely fragmented with a single gene on most chromosomes, and encodes introns that are generally small and lack a conserved branch point for pre-mRNA splicing. Gene family expansion analyses indicate that multiple gene families involved in the encystment are expanded during the evolution of P. cristata. Furthermore, genomic comparisons with other five representative hypotrichs indicate that gene families of phosphorelay sensor kinase, which play a role in the two-component signal transduction system that is related to encystment, show significant expansion among all six hypotrichs. Additionally, cyst wall-related chitin synthase genes have experienced structural changes that increase them from single-exon to multi-exon genes during evolution. These genomic features potentially promote the encystment in hypotrichs and enhance their ability to survive in adverse environments during evolution. CONCLUSIONS: We systematically investigated the genomic structure of hypotrichs and key evolutionary phenomenon, gene family expansion, for encystment promotion in ciliates. In summary, our results provided insights into the evolutionary mechanism of encystment in ciliates.

Whole-Genome Mapping of Epigenetic Modification of 5-Formylcytosine at Single-Base Resolution by Chemical Labeling Enrichment and Deamination Sequencing.

DNA cytosine methylation (5-methylcytosine, 5mC) is a predominant epigenetic modification that plays a critical role in a variety of biological and pathological processes in mammals. In active DNA demethylation, the 10-11 translocation (TET) dioxygenases can sequentially oxidize 5mC to generate three modified forms of cytosine, 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC). Beyond being a demethylation intermediate, recent studies have shown that 5fC has regulatory functions in gene expression and chromatin organization. While some methods have been developed to detect 5fC, genome-wide mapping of 5fC at base resolution is still highly desirable. Herein, we propose a chemical labeling enrichment and deamination sequencing (CLED-seq) method for detecting 5fC in genomic DNA at single-base resolution. The CLED-seq method utilizes selective labeling and enrichment of 5fC-containing DNA fragments, followed by deamination mediated by apolipoprotein B mRNA-editing catalytic polypeptide-like 3A (APOBEC3A or A3A) and sequencing. In the CLED-seq process, while all C, 5mC, and 5hmC are interpreted as T during sequencing, 5fC is still read as C, enabling the precise detection of 5fC in DNA. Using the proposed CLED-seq method, we accomplished genome-wide mapping of 5fC in mouse embryonic stem cells. The mapping study revealed that promoter regions enriched with 5fC overlapped with H3K4me1, H3K4me3, and H3K27ac marks. These findings suggest a correlation between 5fC marks and active gene expression in mESCs. In conclusion, CLED-seq is a straightforward, bisulfite-free method that offers a valuable tool for detecting 5fC in genomes at a single-base resolution.

Whole-genome mapping of APOBEC mutagenesis in metastatic urothelial carcinoma identifies driver hotspot mutations and a novel mutational signature.

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) enzymes mutate specific DNA sequences and hairpin-loop structures, challenging the distinction between passenger and driver hotspot mutations. Here, we characterized 115 whole genomes of metastatic urothelial carcinoma (mUC) to identify APOBEC mutagenic hotspot drivers. APOBEC-associated mutations were detected in 92% of mUCs and were equally distributed across the genome, while APOBEC hotspot mutations (ApoHMs) were enriched in open chromatin. Hairpin loops were frequent targets of didymi (twins in Greek), two hotspot mutations characterized by the APOBEC SBS2 signature, in conjunction with an uncharacterized mutational context (Ap[C>T]). Next, we developed a statistical framework that identified ApoHMs as drivers in coding and non-coding genomic regions of mUCs. Our results and statistical framework were validated in independent cohorts of 23 non-metastatic UCs and 3,744 samples of 17 metastatic cancers, identifying cancer-type-specific drivers. Our study highlights the role of APOBEC in cancer development and may contribute to developing novel targeted therapy options for APOBEC-driven cancers.