Low Mutation Rate Research Articles

Complete Plastid Genome Sequences of Three Tropical African Beilschmiediineae Trees (Lauraceae: Crytocaryeae)

Millions of years of isolation have given Madagascar a unique flora that still reflects some of its relationship with the continents of Africa and India. Here, the complete chloroplast sequence of Beilschmiedia moratii, a tropical tree in Madagascar, was determined. The plastome, with a length of 158,410 bp, was 143 bp and 187 bp smaller than those of two closely related species, B. pierreana and Potameia microphylla, in sub-Saharan Africa and Madagascar with published sequences, respectively. A total of 124 repeats and 114 simple sequence repeats (SSRs) were detected in the plastome of B. moratii. Six highly variable regions, including ndhF, ndhF-rpl32, trnC-petN, pebE-petL, rpl32-trnL, and ycf1, among the three African species were identified and 1151 mutation events, including 14 SVs, 351 indels, and 786 substitutions, were accurately located. There were 634 mutation events between B. moratii and P. microphylla with a mean nucleotide variability (π) value of 0.00279, while there were 827 mutation events between B. moratii and B. pierreana with a mean π value of 0.00385. The Ka/Ks ratios of 86 protein-coding genes in the three African species were less than 1, and the mean value between B. moratii and P. microphylla was 0.184, while the mean value between B. moratii and B. pierreana was 0.286. In this study, the plastid genomes of the three African Beilschmiediineae species were compared for the first time and revealed that B. moratii and P. microphylla from Madagascar were relatively conserved, with low mutation rates and slower evolutionary rates.

Non-Poissonian bursts in the arrival of phenotypic variation can strongly affect the dynamics of adaptation.

Modelling the rate at which adaptive phenotypes appear in a population is a key to predicting evolutionary processes. Given random mutations, should this rate be modelled by a simple Poisson process, or is a more complex dynamics needed? Here we use analytic calculations and simulations of evolving populations on explicit genotype-phenotype maps to show that the introduction of novel phenotypes can be 'bursty' or overdispersed. In other words, a novel phenotype either appears multiple times in quick succession, or not at all for many generations. These bursts are fundamentally caused by statistical fluctuations and other structure in the map from genotypes to phenotypes. Their strength depends on population parameters, being highest for 'monomorphic' populations with low mutation rates. They can also be enhanced by additional inhomogeneities in the mapping from genotypes to phenotypes. We mainly investigate the effect of bursts using the well-studied genotype-phenotype map for RNA secondary structure, but find similar behaviour in a lattice protein model and in Richard Dawkins's biomorphs model of morphological development. Bursts can profoundly affect adaptive dynamics. Most notably, they imply that fitness differences play a smaller role in determining which phenotype fixes than would be the case for a Poisson process without bursts.

Clonal expansion of shared T cell receptors reveals the existence of immune commonality among different lesions of synchronous multiple primary lung cancer

The increase in the detection rate of synchronous multiple primary lung cancer (MPLC) has posed remarkable clinical challenges due to the limited understanding of its pathogenesis and molecular features. Here, comprehensive comparisons of genomic and immunologic features between MPLC and solitary lung cancer nodule (SN), as well as different lesions of the same patient, were performed. Compared with SN, MPLC displayed a lower rate of EGFR mutation but higher rates of BRAF, MAP2K1, and MTOR mutation, which function exactly in the upstream and downstream of the same signaling pathway. Considerable heterogeneity in T cell receptor (TCR) repertoire exists among not only different patients but also among different lesions of the same patient. Invasive lesions of MPLC exhibited significantly higher TCR diversity and lower TCR expansion than those of SN. Intriguingly, different lesions of the same patient always shared a certain proportion of TCR clonotypes. Significant clonal expansion could be observed in shared TCR clonotypes, particularly in those existing in all lesions of the same patient. In conclusion, this study provided evidences of the distinctive mutational landscape, activation of oncogenic signaling pathways, and TCR repertoire in MPLC as compared with SN. The significant clonal expansion of shared TCR clonotypes demonstrated the existence of immune commonality among different lesions of the same patient and shed new light on the individually tailored precision therapy for MPLC.

Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae).

Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.

Endoplasmic reticulum stress-related genes as prognostic and immunogenic biomarkers in prostate cancer

BackgroundThe metastasis and aggressive nature of prostate cancer (PCa) has become a major malignancy related threat that concerns men’s health. The efficacy of immune monotherapy against PCa is questionable due to its lymphocyte-suppressive nature.MethodEndoplasmic reticulum stress- (ERS-) and PCa-prognosis-related genes were obtained from the Molecular Signatures Database and the Cancer Genome Atlas database. The expression, prognosis and immune infiltration values of key genes were explored by “survival R package”, “rms”, “xCELL algorithm”, and univariate–multivariate Cox and LASSO regression analyses. The “consensus cluster plus R package” was used for cluster analysis.ResultAs ERS-related genes, ERLIN2 and CDK5RAP3 showed significant expressional, prognostic and clinic-pathologic values. They were defined as the key genes significantly correlated with immune infiltration and response. The nomogram was constructed with T-stage and primary treatment outcome, and the risk-prognostic model was constructed in the following way: Riskscore = (− 0.1918) * ERLIN2 + (0.5254) * CDK5RAP3. Subsequently, prognostic subgroups based on key genes classified the high-risk group as a pro-cancer subgroup that had lower mutation rates of critical genes (SPOP and MUC16), multiple low-expression immune-relevant molecules, and differences in macrophages (M1 and M2) expressions. Finally, ERLIN2 as an anti-oncogene and CDK5RAP3 as a pro-oncogene were further confirmed by cell phenotype assays and immunohistochemistry.ConclusionWe identified ERLIN2 and CDK5RAP3 as ERS-related genes with important prognostic and immunologic values, and classified patients between high- and low-risk subgroups, which provided new prognostic markers, immunotherapeutic targets, and basis for prognostic assessments.

Spontaneous mutations and mutational responses to penicillin treatment in the bacterial pathogen Streptococcus pneumoniae D39

Bacteria with functional DNA repair systems are expected to have low mutation rates due to strong natural selection for genomic stability. However, our study of the wild-type Streptococcus pneumoniae D39, a pathogen responsible for many common diseases, revealed a high spontaneous mutation rate of 0.02 per genome per cell division in mutation-accumulation (MA) lines. This rate is orders of magnitude higher than that of other non-mutator bacteria and is characterized by a high mutation bias in the A/T direction. The high mutation rate may have resulted from a reduction in the overall efficiency of selection, conferred by the tiny effective population size in nature. In line with this, S. pneumoniae D39 also exhibited the lowest DNA mismatch-repair (MMR) efficiency among bacteria. Treatment with the antibiotic penicillin did not elevate the mutation rate, as penicillin did not induce DNA damage and S. pneumoniae lacks a stress response pathway. Our findings suggested that the MA results are applicable to within-host scenarios and provide insights into pathogen evolution.

How mutation accumulation depends on the structure of the cell lineage tree.

All the cells of a multicellular organism are the product of cell divisions that trace out a single binary tree, the so-called cell lineage tree. Because cell divisions are accompanied by replication errors, the shape of the cell lineage tree is a key determinant of how somatic evolution, which can potentially lead to cancer, proceeds. Carcinogenesis requires the accumulation of a certain number of driver mutations. By mapping the accumulation of mutations into a graph theoretical problem, we present an exact numerical method to calculate the probability of collecting a given number of mutations and show that for low mutation rates it can be approximated with a simple analytical formula, which depends only on the distribution of the lineage lengths, and is dominated by the longest lineages. Our results are crucial in understanding how natural selection can shape the cell lineage trees of multicellular organisms and curtail somatic evolution.

Low mutation rate of spontaneous mutants enables detection of causative genes by comparing whole genome sequences.

In the early 1900s, mutation breeding to select varieties with desirable traits using spontaneous mutation was actively conducted around the world, including Japan. In rice, the number of fixed mutations per generation was estimated to be 1.38-2.25. Although this low mutation rate was a major problem for breeding in those days, in the modern era with the development of next-generation sequencing (NGS) technology, it was conversely considered to be an advantage for efficient gene identification. In this paper, we proposed an in silico approach using NGS to compare the whole genome sequence of a spontaneous mutant with that of a closely related strain with a nearly identical genome, to find polymorphisms that differ between them, and to identify the causal gene by predicting the functional variation of the gene caused by the polymorphism. Using this approach, we found four causal genes for the dwarf mutation, the round shape grain mutation and the awnless mutation. Three of these genes were the same as those previously reported, but one was a novel gene involved in awn formation. The novel gene was isolated from Bozu-Aikoku, a mutant of Aikoku with the awnless trait, in which nine polymorphisms were predicted to alter gene function by their whole-genome comparison. Based on the information on gene function and tissue-specific expression patterns of these candidate genes, Os03g0115700/LOC_Os03g02460, annotated as a short-chain dehydrogenase/reductase SDR family protein, is most likely to be involved in the awnless mutation. Indeed, complementation tests by transformation showed that it is involved in awn formation. Thus, this method is an effective way to accelerate genome breeding of various crop species by enabling the identification of useful genes that can be used for crop breeding with minimal effort for NGS analysis.

InDEL instability in two different tumoral tissues and its forensic significance.

There may be cases where malignant tumor samples can be used for forensic DNA profiling studies. STRs are the first systems preferred in forensic science laboratories for identification purposes. However, genetic instability in tumoral tissues causes STR polymorphism to change, leading to erroneous results. On the other hand, insertion/deletion polymorphism (InDels) are used as genetic markers in forensic science, as they have features that make both STR and SNPs preferable. Although previous studies approved that STR instability is observed in many different tumors, there are only a few studies that have displayed the instability of InDels in tumoral tissues before. In this study, it was aimed to determine whether instability is observed in formalin-fixed paraffin-embedded breast and thyroid tumoral tissues at 36plex InDel Panel. A total of 47 cases, 26 of which were diagnosed as breast cancer and 21 as thyroid cancer, were included in the study. In 21 of 26 (80.76%) breast cancers mutational changes were observed, however only 6 of 21 (28.57%) thyroid carcinoma cases displayed instability.Moreover, in these six cases, mutations were detected at only 1 or 2 loci. The most common change in both tissues was loss of heterozygosity. These findings suggest that paraffin embedded tissues of thyroid tumor can be used in cases of forensic genetic identification, however paraffin embedded breast cancer tissues should be examined with care. In conclusion, low InDel mutation rates compared to STR instability, make InDel analysis from paraffin blocks suitable for forensic genetic identification. However, researchers should keep in mind that there may be differences between the profiles of the tumoral tissues taken as reference and the actual case. In addition, by incorporating additional markers such as SNPs and microhaplotypes with low mutation rates into the study alongside Indels, researchers can significantly enhance the discrimination power in identification processes.

Major Role of S-Glycoprotein in Providing Immunogenicity and Protective Immunity in mRNA Lipid Nanoparticle Vaccines Based on SARS-CoV-2 Structural Proteins.

SARS-CoV-2 variants have evolved over time in recent years, demonstrating immune evasion of vaccine-induced neutralizing antibodies directed against the original S protein. Updated S-targeted vaccines provide a high level of protection against circulating variants of SARS-CoV-2, but this protection declines over time due to ongoing virus evolution. To achieve a broader protection, novel vaccine candidates involving additional antigens with low mutation rates are currently needed. Based on our recently studied mRNA lipid nanoparticle (mRNA-LNP) platform, we have generated mRNA-LNP encoding SARS-CoV-2 structural proteins M, N, S from different virus variants and studied their immunogenicity separately or in combination in vivo. As a result, all mRNA-LNP vaccine compositions encoding the S and N proteins induced excellent titers of RBD- and N-specific binding antibodies. The T cell responses were mainly specific CD4+ T cell lymphocytes producing IL-2 and TNF-alpha. mRNA-LNP encoding the M protein did not show a high immunogenicity. High neutralizing activity was detected in the sera of mice vaccinated with mRNA-LNP encoding S protein (alone or in combinations) against closely related strains, but was undetectable or significantly lower against an evolutionarily distant variant. Our data showed that the addition of mRNAs encoding S and M antigens to mRNA-N in the vaccine composition enhanced the immunogenicity of mRNA-N and induced a more robust immune response to the N protein. Based on our results, we suggested that the S protein plays a key role in enhancing the immune response to the N protein when they are both encoded in the mRNA-LNP vaccine.

Detection of circulating tumor DNA in plasma of patients with primary CNS lymphoma by digital droplet PCR

BackgroundPrimary central nervous system lymphoma (PCNSL) are rare mature B-cell lymphoproliferative diseases characterized by a high incidence of MYD88 L265P and CD79B Y196 hotspot mutations. Diagnosis of PCNSL can be challenging. The aim of the study was to analyze the detection rate of the MYD88 L265P and CD79B Y196 mutation in cell free DNA (cfDNA) in plasma of patients with PCNSL.MethodsWe analyzed by digital droplet PCR (ddPCR) to determine presence of the MYD88 L265P and CD79B Y196 hotspot mutations in cfDNA isolated from plasma of 24 PCNSL patients with active disease. Corresponding tumor samples were available for 14 cases. Based on the false positive rate observed in 8 healthy control samples, a stringent cut-off for the MYD88 L265P and CD79B Y196 mutation were set at 0.3% and 0.5%, respectively.ResultsMYD88 L265P and CD79B Y196 mutations were detected in 9/14 (64%) and 2/13 (15%) tumor biopsies, respectively. In cfDNA samples, the MYD88 L265P mutation was detected in 3/24 (12.5%), while the CD79B Y196 mutation was not detected in any of the 23 tested cfDNA samples. Overall, MYD88 L265P and/or CD79B Y196 were detected in cfDNA in 3/24 cases (12.5%). The detection rate of the combined analysis did not improve the single detection rate for either MYD88 L265P or CD79B Y196.ConclusionThe low detection rate of MYD88 L265P and CD79B Y196 mutations in cfDNA in the plasma of PCNSL patients argues against its use in routine diagnostics. However, detection of MYD88 L265P by ddPCR in cfDNA in the plasma could be considered in challenging cases.

Gene Profile and Clinical Significance of Concomitant Mutations in CN-AML Patients with CEBPA Mutation

To analyze the occurrence of concomitant gene mutations in cytogenetically normal acute myeloid leukemia (CN-AML) patients with CEBPA mutation and its impact on the clinical characteristics and prognosis of the patients. 151 newly diagnosed patients with CN-AML in the Second Hospital of Shanxi Medical University from June 2013 to June 2020 were analyzed retrospectively. 34 common genetic mutations associated with hematologic malignancies were detected by next-generation sequencing technology. The occurrence of concomitant gene mutations in patients with CEBPA positive and negative groups was compared, and the correlation between concomitant mutations in different functional groups and the clinical characteristics and prognosis of CN-AML patients with CEBPA mutation was analyzed. In 151 patients with CN-AML, 55 (36.42%) were positive for CEBPA mutation (including 36 cases of CEBPAdm and 19 cases of CEBPAsm), of which 41 (74.55%) had co-mutations with other genes. The main mutated genes were GATA2 (25.45%, 14/55), TET2 (21.82%, 12/55), FLT3 (20.00%, 11/55), NRAS (12.73%, 7/55) and WT1 (9.09%, 9/55), etc. Some cases had two or more concomitant gene mutations. Grouping the mutant genes according to their functions showed that CEBPA+ group had lower mutation rates of histone methylation (P =0.002) and chromatin modification genes (P =0.002, P =0.033), and higher mutation rates of transcription factors (P =0.037) than CEBPA- group. In 55 patients with CEBPA+ CN-AML, the platelet count at diagnosis in signaling pathway gene mutation-positive group was lower than that in the mutation-negative group (P =0.005), the proportion of bone marrow blasts in transcription factor mutation-positive group was higher than that in the mutation-negative group (P =0.003), and the onset age in DNA methylation gene mutation-positive group and chromatin modifier mutation-positive group was older than that in the mutation-negative group, respectively (P =0.002, P =0.008). DFS of CEBPA+ CN-AML patients in signaling pathway gene mutation group was shorter than that in signaling pathway gene mutation-negative group (median DFS: 12 months vs not reached) (P =0.034). Compared with DNA methylation gene mutation-negative group, CEBPA+ CN-AML patients with DNA methylation gene mutation had lower CR rate (P =0.025) significantly shorter OS and DFS (median OS: 20 months vs not reached, P =0.006; median DFS: 15 months vs not reached, P =0.049). OS in patients with histone methylation gene mutation was significantly shorter than that in the histone methylation gene mutation-negative group (median OS: 12 months vs 40 months) (P =0.008). Multivariate analysis of prognostic factors showed that the proportion of bone marrow blasts (P =0.046), concomitant DNA methylation gene mutation (P =0.006) and histone methylation gene mutation (P =0.036) were independent risk factors affecting the prognosis. CN-AML patients with CEBPA mutation have specific concomitant gene profile, and the concomitant mutations of different functional genes have a certain impact on the clinical characteristics and prognosis of the patients.

Bacillus subtilis stress-associated mutagenesis and developmental DNA repair.

SUMMARYThe metabolic conditions that prevail during bacterial growth have evolved with the faithful operation of repair systems that recognize and eliminate DNA lesions caused by intracellular and exogenous agents. This idea is supported by the low rate of spontaneous mutations (10-9) that occur in replicating cells, maintaining genome integrity. In contrast, when growth and/or replication cease, bacteria frequently process DNA lesions in an error-prone manner. DNA repairs provide cells with the tools needed for maintaining homeostasis during stressful conditions and depend on the developmental context in which repair events occur. Thus, different physiological scenarios can be anticipated. In nutritionally stressed bacteria, different components of the base excision repair pathway may process damaged DNA in an error-prone approach, promoting genetic variability. Interestingly, suppressing the mismatch repair machinery and activating specific DNA glycosylases promote stationary-phase mutations. Current evidence also suggests that in resting cells, coupling repair processes to actively transcribed genes may promote multiple genetic transactions that are advantageous for stressed cells. DNA repair during sporulation is of interest as a model to understand how transcriptional processes influence the formation of mutations in conditions where replication is halted. Current reports indicate that transcriptional coupling repair-dependent and -independent processes operate in differentiating cells to process spontaneous and induced DNA damage and that error-prone synthesis of DNA is involved in these events. These and other noncanonical ways of DNA repair that contribute to mutagenesis, survival, and evolution are reviewed in this manuscript.

KRASG12C Inhibitor as a Treatment Option for Non-Small-Cell Lung Cancer with Comorbid Interstitial Pneumonia.

Non-small-cell lung cancer (NSCLC) with comorbid interstitial pneumonia (IP) is a population with limited treatment options and a poor prognosis. Patients with comorbid IP are at high risk of developing fatal drug-induced pneumonitis, and data on the safety and efficacy of molecularly targeted therapies are lacking. KRAS mutations have been frequently detected in patients with NSCLC with comorbid IP. However, the low detection rate of common driver gene mutations, such as epidermal growth factor receptor and anaplastic lymphoma kinase, in patients with comorbid IP frequently results in inadequate screening for driver mutations, and KRAS mutations may be overlooked. Recently, sotorasib and adagrasib were approved as treatment options for advanced NSCLC with KRASG12C mutations. Although patients with comorbid IP were not excluded from clinical trials of these KRASG12C inhibitors, the incidence of drug-induced pneumonitis was low. Therefore, KRASG12C inhibitors may be a safe and effective treatment option for NSCLC with comorbid IP. This review article discusses the promise and prospects of molecular-targeted therapies, especially KRASG12C inhibitors, for NSCLC with comorbid IP, along with our own clinical experience.

Abstract 2151: Additional impact of genetic ancestry over self-reported race/ethnicity to prevalence of KRAS mutations and allele-specific subtypes in non-small cell lung cancer

Abstract Background: KRAS mutation is the most common oncogenic driver in patients with non-small cell lung cancer (NSCLC). However, the interplay between self-reported race and/or ethnicity (SIRE) and genetically inferred ancestry (GIA) on KRAS mutations and allele-specific subtypes is largely unknown. Methods: A total of 3918 NSCLC patients from the Chinese OrigiMed (OM) cohort (n=2039) and Boston Lung Cancer Survival (BLCS) cohort (n=1879) were included, all of whom had baseline covariates including age, sex, clinical stage, SIRE, histology and smoking status. KRAS mutation and subtypes were determined by targeted NGS panels in each cohort, and KRAS mutation outcomes, including KRAS mutation, KRAS subtypes, and KRAS transversion and transition were analyzed. GIA was inferred from paired germline data for patients in the BLCS cohort. Multivariable and multinomial logistic regressions were employed to assess the associations between SIRE, GIA and KRAS outcomes. Results: In the entire cohort, 804 (20.5%) patients harbored KRAS mutation, with a considerably higher prevalence of 31.5% observed in BLCS compared to 10.4% in OM cohort. Patients with self-identified Asian was associated with a significantly lower rate of KRAS mutation (OR = 0.44, 95% CI:0.22 - 0.81, P = 0.01), transversion substitutions (OR = 0.24, 95% CI:0.10 - 0.62, P = 0.003) and KRASG12C (OR = 0.17, 95% CI:0.04 - 0.72, P = 0.02) compared to SIRE-White patients after adjusting for confounders including smoking status, sex and histology. The association of SIRE and KRAS driver mutation was further modified by sex, with Asian women demonstrating lower risk compared to Asian men, while the opposite trend was observed in SIRE-White patients. Increased European ancestry was associated with a higher risk of KRAS mutation (OR per 10% increase = 1.05, 95% CI:1.00 - 1.09, P = 0.03), particularly with more transition substitutions (OR per 10% increase = 1.10, 95% CI:1.01 - 1.21, P = 0.03) and KRASG12D (OR per 10% increase = 1.18, 95% CI:1.04 - 1.34, P = 0.01). This association was even more pronounced in patients of increased Northern European ancestry. Furthermore, even among SIRE-White patients, a 10% increase in European ancestry was associated with an additional 6% increase in the likelihood of KRAS mutations (OR = 1.06, 95% CI: 1.01 - 1.12, P = 0.03), while a 10% increase in Admixed American ancestry was associated with a 6% decreased likelihood of KRAS mutations (OR = 0.94, 95% CI: 0.88 - 0.99, P = 0.04). Results were validated in an independent cohort of 1450 patients. Conclusion: Self-reported race/ethnicity and genetic ancestry are both associated with KRAS driver mutations and allele-specific subtypes of NSCLC. This work highlights the significance of integrating both into population-based cancer research, aiming to refine clinical decision-making processes and mitigate health disparities. Citation Format: Xinan Wang, Kangcheng Hou, Biagio Ricciuti, Joao V. Alessi, Xihao Li, Federica Pecci, Rounak Dey, Jia Luo, Mark M. Awad, Alexander Gusev, Xihong Lin, Bruce E. Johnson, David C. Christiani. Additional impact of genetic ancestry over self-reported race/ethnicity to prevalence of KRAS mutations and allele-specific subtypes in non-small cell lung cancer [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 2151.

Abstract 3861: Associations of intragenic rearrangement burden with immune cell infiltration and response to immune checkpoint blockade in cancer

Abstract While immune checkpoint blockade (ICB) induces durable cancer remission, only a small subset of patients can receive benefits, whereas a substantial portion of patients encounter severe immune-related adverse events. Further, about 4-29% of patients experience dramatic acceleration of disease progression after ICB treatment. Predictive biomarkers of ICB responses are in high demand to ensure patient benefits and avoid tragic clinical outcomes. In addition to PD-L1 tests, several genetic biomarkers are associated with ICB effectiveness, such as tumor mutation burden (TMB), somatic copy number aberrations, indels, and microsatellite instability. However, none of these neoantigen markers are sufficient to discriminate ICB responders in diverse cancer entities of distinct genetic profiles. In particular, TMB fails to predict in the cancer types with low simple mutation rates, such as triple-negative breast cancer (TNBC), ovarian cancer, and esophageal cancers. In TNBC, high TMB only accounts for 5-10% of patients, whereas a much larger proportion of tumors are inflamed (26-50%), which suggests a possible uncharacterized major source of neoantigens. Similarly, TMB is not predictive of ICB response in ovarian cancers according to the IMagyn050 trial. In esophageal cancer, TMB is not associated with T cell infiltration and does not predict ICB response. Our analyses of Pan-cancer whole genomes revealed that intragenic rearrangement (IGR) burden is a pivotal contributor to immune infiltration in breast, ovarian, esophageal, and endometrial cancers, which correlates with increased Macrophage M1 and CD8+ T cell signatures. IGRs result in exon duplications or deletions through rearrangements of a genomic distance of tens of kilobases. Multivariate regression against spatially counted tumor infiltrated lymphocytes in breast, endometrial, and ovarian cancers suggest that IGR burden is a more influential covariate than other genetic aberrations in these cancers. In the MEDI4736 trial evaluating durvalumab in esophageal adenocarcinoma, IGR burden correlates with patient benefits. In the IMVigor210 trial evaluating atezolizumab in urothelial carcinoma, IGR burden increases with platinum exposure and predicts patient benefits for TMB-low, platinum-exposed tumors. Altogether, we demonstrate that IGR burden contributes to T-cell inflammation and predicts ICB benefits in TMB-low, IGR-dominant tumors, or in platinum-exposed tumors. IGRs as a cryptic class of structural variants have been heretofore overlooked in cancer genomics studies. The association of IGRs with ICB response has not been previously reported in cancer in general. The development of IGR burden-based predictive marker could increase effectiveness of ICB therapy, reduce severe side-effects and unnecessary costs. Citation Format: Han Zhang, Sanghoon Lee, Renee R. Muthakana, Binfeng Lu, David N Boone, Daniel Lee, Xiaosong Wang. Associations of intragenic rearrangement burden with immune cell infiltration and response to immune checkpoint blockade in cancer [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 3861.

A Study of Alternative TrkA Splicing Identifies TrkAIII as a Novel Potentially Targetable Participant in PitNET Progression.

Pituitary neuroendocrine tumors (PitNETs) are generally benign but comprise an aggressive, invasive, therapy-resistant, metastatic subset, underpinning a need for novel therapeutic targets. PitNETs exhibit low mutation rates but are associated with conditions linked to alternative splicing, an alternative oncogene pathway activation mechanism. PitNETs express the neurotrophin receptor TrkA, which exhibits oncogenic alternative TrkAIII splicing in other neuroendocrine tumors. We, therefore, assessed whether TrkAIII splicing represents a potential oncogenic participant in PitNETs. TrkAIII splicing was RT-PCR assessed in 53 PitNETs and TrkA isoform(s) expression and activation were assessed by confocal immunofluorescence. TrkAIII splicing was also compared to HIF1α, HIF2α, SF3B1, SRSF2, U2AF1, and JCPyV large T antigen mRNA expression, Xbp1 splicing, and SF3B1 mutation. TrkAIII splicing was detected in all invasive and most non-invasive PitNETs and was significantly elevated in invasive cases. In PitNET lineages, TrkAIII splicing was significantly elevated in invasive PIT1 PitNETs and high in invasive and non-invasive SF1 and TPIT lineages. Immunoreactivity consistent with TrkAIII activation characterized PitNET expressing TrkAIII mRNA, and invasive Pit1 PitNETs exhibited elevated HIF2α expression. TrkAIII splicing did not associate with SF3B1 mutations, altered SF3B1, SRSF2, and U2AF1 or JCPyV large T antigen expression, or Xbp1 splicing. Therefore, TrkAIII splicing is common in PitNETs, is elevated in invasive, especially PIT1 tumors, can result in intracellular TrkAIII activation, and may involve hypoxia. The data support a role for TrkAIII splicing in PitNET pathogenesis and progression and identify TrkAIII as a novel potential target in refractory PitNETs.

Comparative analysis of 326 chloroplast genomes in Chinese jujube (Ziziphus jujuba): Structural variations, horizontal gene transfer events, and evolutionary patterns impacting its domestication from wild jujube

AbstractJujube (Ziziphus jujuba Mill.), renowned for its nutritional value and health benefits, is believed to have originated in the middle and lower reaches of the Yellow River in China, where it underwent domestication from wild jujube. Nonetheless, the evolutionary trajectory and species differentiation between wild jujube and cultivated jujube still require further elucidation. The chloroplast genome (plastome), characterized by its relatively lower mutation rate compared to the nuclear genome, serves as an excellent model for evolutionary and comparative genomic research. In this study, we analyzed 326 nonredundant plastomes, encompassing 133 jujube cultivars and 193 wild jujube genotypes distributed throughout China. Noteworthy variations in the large single copy region primarily account for the size differences among these plastomes, impacting the evolution from wild jujube to cultivated varieties. Horizontal gene transfer (HGT) unveiled a unique chloroplast‐to‐nucleus transfer event, with transferred fragments predominantly influencing the evolution of the nuclear genome while leaving the plastome relatively unaffected. Population genetics analysis revealed two distinct evolutionary pathways from wild jujube to cultivated jujube: one driven by natural selection with minimal human interference, and the other resulting from human domestication and cultivation. Molecular dating, based on phylogenetic analysis, supported the likelihood that wild jujube and cultivated jujube fall within the same taxonomic category, Z. jujuba. In summary, our study comprehensively examined jujube plastome structures and HGT events, simultaneously contributing novel insights into the intricate processes that govern the evolution and domestication of jujube species.

Within-host evolution of SARS-CoV-2: how often are de novo mutations transmitted from symptomatic infections?

Despite a relatively low mutation rate, the large number of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections has allowed for substantial genetic change, leading to a multitude of emerging variants. Using a recently determined mutation rate (per site replication), as well as within-host parameter estimates for symptomatic SARS-CoV-2 infection, we apply a stochastic transmission-bottleneck model to describe the survival probability of de novo SARS-CoV-2 mutations as a function of bottleneck size and selection coefficient. For narrow bottlenecks, we find that mutations affecting per-target-cell attachment rate (with phenotypes associated with fusogenicity and ACE2 binding) have similar transmission probabilities to mutations affecting viral load clearance (with phenotypes associated with humoral evasion). We further find that mutations affecting the eclipse rate (with phenotypes associated with reorganization of cellular metabolic processes and synthesis of viral budding precursor material) are highly favoured relative to all other traits examined. We find that mutations leading to reduced removal rates of infected cells (with phenotypes associated with innate immune evasion) have limited transmission advantage relative to mutations leading to humoral evasion. Predicted transmission probabilities, however, for mutations affecting innate immune evasion are more consistent with the range of clinically estimated household transmission probabilities for de novo mutations. This result suggests that although mutations affecting humoral evasion are more easily transmitted when they occur, mutations affecting innate immune evasion may occur more readily. We examine our predictions in the context of a number of previously characterized mutations in circulating strains of SARS-CoV-2. Our work offers both a null model for SARS-CoV-2 mutation rates and predicts which aspects of viral life history are most likely to successfully evolve, despite low mutation rates and repeated transmission bottlenecks.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution.

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.