Nonsynonymous Sites Research Articles

Enhanced Detection and Molecular Modeling of Adaptive Mutations in SARS-CoV-2 Coding and Non-Coding Regions Using the c/µ Test

Abstract Accurately identifying mutations under beneficial selection in viral genomes is crucial for understanding their molecular evolution and pathogenicity. Traditional methods like the Ka/Ks test, which assesses non-synonymous (Ka) versus synonymous (Ks) substitution rates, assume that synonymous substitutions at synonymous sites are neutral and thus is equal to the mutation rate (µ). Yet, evidence suggests that synonymous sites in translated (TR) regions and untranslated regions (UTR) can be under strong beneficial selection (Ks>µ) and strongly conserved (Ks ≈ 0), leading to false predictions of adaptive mutations from codon-by-codon Ka/Ks analysis. Our previous work used a relative substitution rate test (c/µ, c: substitution rate in UTR/TR, and µ: mutation rate) to identify adaptive mutations in SARS-CoV-2 genome without the neutrality assumption of the synonymous sites. This study refines the c/µ test by optimizing µ value, leading to a smaller set of nucleotide and amino acid sites under beneficial selection in both UTR (11 sites with c/µ>3) and TR (69 nonsynonymous sites: c/µ>3 and Ka/Ks>2.5; 107 synonymous sites: Ks/µ>3). Encouragingly, the top 2 mutations in UTR and 70% of the top nonsynonymous mutations in TR had reported or predicted effects in the literature. Molecular modeling of top adaptive mutations for some critical proteins (S, NSP11 and NSP5) was carried out to elucidate the possible molecular mechanism of their adaptivity.

Purifying Selection Influences the Comparison of Heterozygosities between Populations.

Heterozygosity is a fundamental measure routinely used to compare between populations to infer the level of genetic variation and their relative effective population sizes. However, such comparison is highly influenced by the magnitude of selection pressure on the genomic regions used. Using over 2 million Single Nucleotide Variants (SNVs) from chimpanzee and mouse populations, this study shows that the heterozygosities estimated using neutrally evolving sites of large populations were two times higher than those of small populations. However, this difference was only ~1.6 times for the heterozygosities estimated using nonsynonymous sites. This suggests an excess in the nonsynonymous heterozygosities due to the segregation of deleterious variants in small populations. This excess in the nonsynonymous heterozygosities of the small populations was estimated to be 23-31%. Further analysis revealed that the magnitude of the excess is modulated by effective population size (Ne) and selection intensity (s). Using chimpanzee populations, this investigation found that the excess in nonsynonymous diversity in the small population was little (6%) when the difference between the Ne values of large and small populations was small (2.4 times). Conversely, this was high (23%) when the difference in Ne was large (5.9 times). Analysis using mouse populations showed that the excess in the nonsynonymous diversity of highly constrained genes of the small population was much higher (38%) than that observed for the genes under relaxed selective constraints (21%). Similar results were observed when the expression levels of genes were used as a proxy for selection intensity. These results emphasize the use of neutral regions, less constrained genes, or lowly expressed genes when comparing the heterozygosities between populations.

MhGeneS: An Analytical Pipeline to Allow for Robust Microhaplotype Genotyping.

Microhaplotypes are small linked genomic regions comprising two or more single-nucleotide polymorphisms (SNPs) that are being applied in forensics and are emerging in wildlife monitoring studies and genomic epidemiology. Typically, targeted in non-coding regions, microhaplotypes in exonic regions can be designed with larger amplicons to capture functional non-synonymous sites and minimise insertion/deletion (indel) polymorphisms. Quality control is an important first step for high-confidence genotyping to counteract such false-positive variants. As genetic markers with higher polymorphism compared to biallelic SNPs, it is critical to ensure sequencing errors across the microhaplotype amplicon are filtered out to avoid introducing false-haplotypes. We developed the MhGeneS pipeline which works in tandem with Seq2Sat to help validate microhaplotype genotyping of the coding region of genes, with broader applicability to any microhaplotype profiling. We genotyped microhaplotype regions of the Zfx (≅ 160 bp) and Zfy (≅ 140 bp) genes, as well as an exon of the prion protein (Prnp) gene (≅ 370 bp) in caribou (Rangifer tarandus) using paired-end Illumina technology. As important quality metrics affecting microhaplotype calling, we identified the sequencing error rate profile related to the overlap or non-overlap of paired-end reads as well as the read depth as significant. In the case of Prnp, we achieved confident microhaplotype calling through MhGeneS by removing small sections of the 5' and 3' amplicons and using a minimum read depth of 20. Read depth and sequence trimming may be locus-specific, and validation of these parameters is recommended before the high-throughput profiling of samples.

Mitogenomics of a declining species, the Rio Grande silvery minnow (Hybognathus amarus), with boom–bust population dynamics

Abstract Large mitochondrial DNA fragments and nearly complete mitogenomes were obtained from nextera-tagmented reductively amplified DNA (nextRAD) data from the Rio Grande silvery minnow (Hybognathus amarus). We used these data to characterize variation across the mitogenome and to evaluate changes in mitochondrial diversity across a 20 year time series punctuated by marked species decline. We observed different evolutionary rates among protein-coding genes and inflated values of the ratio of non-synonymous change to non-synonymous sites to synonymous changes per synonymous site for several protein-coding genes, violating assumptions of selective neutrality. Analyses suggested that non-synonymous substitutions are effectively neutral, and our results were not consistent with genetic draft. Mitogenomes revealed large long-term effective population size for Rio Grande silvery minnow, and time-series data indicated recent reductions in mitochondrial diversity via genetic drift associated with population bottlenecks. The results presented in this study show that in some cases, nextRAD sequences contain informative mitochondrial DNA data.

Discovering novel germline genetic variants linked to severe fluoropyrimidine-related toxicity in- and outside DPYD

BackgroundThe Alpe-DPD study (NCT02324452) demonstrated that prospective genotyping and dose-individualization using four alleles in DPYD (DPYD*2A/rs3918290, c.1236G > A/rs75017182, c.2846A > T/rs67376798 and c.1679 T > G/rs56038477) can mitigate the risk of severe fluoropyrimidine toxicity. However, this could not prevent all toxicities. The goal of this study was to identify additional genetic variants, both inside and outside DPYD, that may contribute to fluoropyrimidine toxicity.MethodsBiospecimens and data from the Alpe-DPD study were used. Exon sequencing was performed to identify risk variants inside DPYD. In silico and in vitro analyses were used to classify DPYD variants. A genome-wide association study (GWAS) with severe fluoropyrimidine-related toxicity was performed to identify variants outside DPYD. Association with severe toxicity was assessed using matched-pair analyses for the exon sequencing and logistic, Cox, and ordinal regression analyses for GWAS.ResultsTwenty-four non-synonymous, frameshift, and splice site DPYD variants were detected in ten of 986 patients. Seven of these variants (c.1670C > T, c.1913 T > C, c.1925 T > C, c.506delC, c.731A > C, c.1740 + 1G > T, c.763 − 2A > G) were predicted to be deleterious. The carriers of either of these variants showed a trend towards a 2.14-fold (95% CI, 0.41–11.3, P = 0.388) increased risk of severe toxicity compared to matched controls (N = 30). After GWAS of 942 patients, no individual single nucleotide polymorphisms achieved genome-wide significance (P ≤ 5 × 10−8), however, five variants were suggestive of association (P < 5 × 10−6) with severe toxicity.ConclusionsResults from DPYD exon sequencing and GWAS analysis did not identify additional genetic variants associated with severe toxicity, which suggests that testing for single markers at a population level currently has limited clinical value. Identifying additional variants on an individual level is still promising to explain fluoropyrimidine-related severe toxicity. In addition, studies with larger samples sizes, in more diverse cohorts are needed to identify potential clinically relevant genetic variants related to severe fluoropyrimidine toxicity.

Mining for QTL controlling maize low-phosphorus response genes combined with deep resequencing of RIL parental genomes and in silico GWAS analysis.

Extensive and comprehensive phenotypic data from a maize RIL population under both low- and normal-Pi treatments were used to conduct QTL mapping. Additionally, we integrated parental resequencing data from the RIL population, GWAS results, and transcriptome data to identify candidate genes associated with low-Pi stress in maize. Phosphorus (Pi) is one of the essential nutrients that greatly affect the maize yield. However, the genes underlying the QTL controlling maize low-Pi response remain largely unknown. In this study, a total of 38 traits at both seedling and maturity stages were evaluated under low- and normal-Pi conditions using a RIL population constructed from X178 (tolerant) and 9782 (sensitive), and most traits varied significantly between low- and normal-Pi treatments. Twenty-nine QTLs specific to low-Pi conditions were identified after excluding those with common intervals under both low- and normal-Pi conditions. Furthermore, 45 additional QTLs were identified based on the index value ((Trait_under_LowPi-Trait_under_NormalPi)/Trait_under_NormalPi) of each trait. These 74 QTLs collectively were classified as Pi-dependent QTLs. Additionally, 39 Pi-dependent QTLs were clustered in nine HotspotQTLs. The Pi-dependent QTL interval contained 19,613 unique genes, 6,999 of which exhibited sequence differences with non-synonymous mutation sites between X178 and 9782. Combined with in silico GWAS results, 277 consistent candidate genes were identified, with 124 genes located within the HotspotQTL intervals. The transcriptome analysis revealed that 21 genes, including the Pi transporter ZmPT7 and the strigolactones pathway-related gene ZmPDR1, exhibited consistent low-Pi stress response patterns across various maize inbred lines or tissues. It is noteworthy that ZmPDR1 in maize roots can be sharply up-regulated by low-Pi stress, suggesting its potential importance as a candidate gene for responding to low-Pi stress through the strigolactones pathway.

New comparative genomic evidence supporting the proteomic diversification role of A-to-I RNA editing in insects.

Adenosine-to-inosine (A-to-I) RNA editing, resembling A-to-G mutation, confers adaptiveness by increasing proteomic diversity in a temporal-spatial manner. This evolutionary theory named "proteomic diversifying hypothesis" has only partially been tested in very few organisms like Drosophila melanogaster, mainly by observing the positive selection on nonsynonymous editing events. To find additional genome-wide evidences supporting this interesting assumption, we retrieved the genomes of four Drosophila species and collected 20 deep-sequenced transcriptomes of different developmental stages and neuron populations of D. melanogaster. We systematically profiled the RNA editomes in these samples and performed meticulous comparative genomic analyses. Further evidences were found to support the diversifying hypothesis. (1) None of the nonsynonymous editing sites in D. melanogaster had ancestral G-alleles, while the silent editing sites had an unignorable fraction of ancestral G-alleles; (2) Only very few nonsynonymous editing sites in D. melanogaster had corresponding G-alleles derived in the genomes of sibling species, and the fraction of such situation was significantly lower than that of silent editing sites; (3) The few nonsynonymous editing with corresponding G-alleles had significantly more variable editing levels (across samples) than other nonsynonymous editing sites in D. melanogaster. The proteomic diversifying nature of RNA editing in Drosophila excludes the restorative role which favors an ancestral G-allele. The few fixed G-alleles in sibling species might facilitate the adaptation to particular environment and the corresponding nonsynonymous editing in D. melanogaster would introduce stronger advantage of flexible proteomic diversification. With multi-Omics data, our study consolidates the nature of evolutionary significance of A-to-I RNA editing sites in model insects.

Quantitative trait loci identification reveals zinc finger protein CONSTANS-LIKE 4 as the key candidate gene of stigma color in watermelon (Citrullus lanatus)

Stigma color is a critical agronomic trait in watermelon that plays an important role in pollination. However, there are few reports on the regulation of stigma color in watermelon. In this study, a genetic analysis of the F2 population derived from ZXG1553 (P1, with orange stigma) and W1-17 (P2, with yellow stigma) indicated that stigma color is a quantitative trait and the orange stigma is recessive compared with the yellow stigma. Bulk segregant analysis sequencing (BSA-seq) revealed a 3.75 Mb segment on chromosome 6 that is related to stigma color. Also, a major stable effective QTL Clqsc6.1 (QTL stigma color) was detected in two years between cleaved amplified polymorphic sequencing (CAPS) markers Chr06_8338913 and Chr06_9344593 spanning a ~1.01 Mb interval that harbors 51 annotated genes. Cla97C06G117020 (annotated as zinc finger protein CONSTANS-LIKE 4) was identified as the best candidate gene for the stigma color trait through RNA-seq, quantitative real-time PCR (qRT-PCR), and gene structure alignment analysis among the natural watermelon panel. The expression level of Cla97C06G117020 in the orange stigma accession was lower than in the yellow stigma accessions with a significant difference. A nonsynonymous SNP site of the Cla97C06G117020 coding region that causes amino acid variation was related to the stigma color variation among nine watermelon accessions according to their re-sequencing data. Stigma color formation is often related to carotenoids, and we also found that the expression trend of ClCHYB (annotated as β-carotene hydroxylase) in the carotenoid metabolic pathway was consistent with Cla97C06G117020, and it was expressed in low amounts in the orange stigma accession. These data indicated that Cla97C06G117020 and ClCHYB may interact to form the stigma color. This study provides a theoretical basis for gene fine mapping and mechanisms for the regulation of stigma color in watermelon.

Identification of Quantitative Trait Loci and Candidate Genes Controlling Seed Dormancy in Eggplant (Solanum melongena L.).

Seed dormancy is a life adaptation trait exhibited by plants in response to environmental changes during their growth and development. The dormancy of commercial seeds is the key factor affecting seed quality. Eggplant seed dormancy is controlled by quantitative trait loci (QTLs), but reliable QTLs related to eggplant dormancy are still lacking. In this study, F2 populations obtained through the hybridization of paternally inbred lines with significant differences in dormancy were used to detect regulatory sites of dormancy in eggplant seeds. Three QTLs (dr1.1, dr2.1, and dr6.1) related to seed dormancy were detected on three chromosomes of eggplant using the QTL-Seq technique. By combining nonsynonymous sites within the candidate regions and gene functional annotation analysis, nine candidate genes were selected from three QTL candidate regions. According to the germination results on the eighth day, the male parent was not dormant, but the female parent was dormant. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the expression of nine candidate genes, and the Smechr0201082 gene showed roughly the same trend as that in the phenotypic data. We proposed Smechr0201082 as the potential key gene involved in regulating the dormancy of eggplant seeds. The results of seed experiments with different concentrations of gibberellin A3 (GA3) showed that, within a certain range, the higher the gibberellin concentration, the earlier the emergence and the higher the germination rate. However, higher concentrations of GA3 may have potential effects on eggplant seedlings. We suggest the use of GA3 at a concentration of 200-250 mg·L-1 to treat dormant seeds. This study provides a foundation for the further exploration of genes related to the regulation of seed dormancy and the elucidation of the molecular mechanism of eggplant seed dormancy and germination.

Purifying selection drives distinctive arsenic metabolism pathways in prokaryotic and eukaryotic microbes.

Microbes play a crucial role in the arsenic biogeochemical cycle through specific metabolic pathways to adapt to arsenic toxicity. However, the different arsenic-detoxification strategies between prokaryotic and eukaryotic microbes are poorly understood. This hampers our comprehension of how microbe-arsenic interactions drive the arsenic cycle and the development of microbial methods for remediation. In this study, we utilized conserved protein domains from 16 arsenic biotransformation genes (ABGs) to search for homologous proteins in 670 microbial genomes. Prokaryotes exhibited a wider species distribution of arsenic reduction- and arsenic efflux-related genes than fungi, whereas arsenic oxidation-related genes were more prevalent in fungi than in prokaryotes. This was supported by significantly higher acr3 (arsenite efflux permease) expression in bacteria (upregulated 3.72-fold) than in fungi (upregulated 1.54-fold) and higher aoxA (arsenite oxidase) expression in fungi (upregulated 5.11-fold) than in bacteria (upregulated 2.05-fold) under arsenite stress. The average values of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site (dN/dS) of homologous ABGs were higher in archaea (0.098) and bacteria (0.124) than in fungi (0.051). Significant negative correlations between the dN/dS of ABGs and species distribution breadth and gene expression levels in archaea, bacteria, and fungi indicated that microbes establish the distinct strength of purifying selection for homologous ABGs. These differences contribute to the distinct arsenic metabolism pathways in prokaryotic and eukaryotic microbes. These observations facilitate a significant shift from studying individual or several ABGs to characterizing the comprehensive microbial strategies of arsenic detoxification.

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes.

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Detection and editing of the updated Arabidopsis plastid- and mitochondrial-encoded proteomes through PeptideAtlas.

Arabidopsis (Arabidopsis thaliana) ecotype Col-0 has plastid and mitochondrial genomes encoding over 100 proteins. Public databases (e.g. Araport11) have redundancy and discrepancies in gene identifiers for these organelle-encoded proteins. RNA editing results in changes to specific amino acid residues or creation of start and stop codons for many of these proteins, but the impact of RNA editing at the protein level is largely unexplored due to the complexities of detection. Here, we assembled the nonredundant set of identifiers, their correct protein sequences, and 452 predicted nonsynonymous editing sites of which 56 are edited at lower frequency. We then determined accumulation of edited and/or unedited proteoforms by searching ∼259 million raw tandem MS spectra from ProteomeXchange, which is part of PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/). We identified all mitochondrial proteins and all except 3 plastid-encoded proteins (NdhG/Ndh6, PsbM, and Rps16), but no proteins predicted from the 4 ORFs were identified. We suggest that Rps16 and 3 of the ORFs are pseudogenes. Detection frequencies for each edit site and type of edit (e.g. S to L/F) were determined at the protein level, cross-referenced against the metadata (e.g. tissue), and evaluated for technical detection challenges. We detected 167 predicted edit sites at the proteome level. Minor frequency sites were edited at low frequency at the protein level except for cytochrome C biogenesis 382 at residue 124 (Ccb382-124). Major frequency sites (>50% editing of RNA) only accumulated in edited form (>98% to 100% edited) at the protein level, with the exception of Rpl5-22. We conclude that RNA editing for major editing sites is required for stable protein accumulation.

The Many Roles of A-to-I RNA Editing in Animals: Functional or Adaptive?

Metazoan adenosine-to-inosine (A-to-I) RNA editing is a highly conserved mechanism that diversifies the transcriptome by post-transcriptionally converting adenosine to inosine. Millions of editing sites have been identified in different species and, based on abnormal editing observed in various disorders, it is intuitive to conclude that RNA editing is both functional and adaptive. In this review, we propose the following major points: (1) "Function/functional" only represents a molecular/phenotypic consequence and is not necessarily connected to "adaptation/adaptive"; (2) Adaptive editing should be judged in the light of evolution and emphasize advantages of temporal-spatial flexibility; (3) Adaptive editing could, in theory, be extended from nonsynonymous sites to all potentially functional sites. This review seeks to conceptually bridge the gap between molecular biology and evolutionary biology and provide a more objective understanding on the biological functions and evolutionary significance of RNA editing.

Differential adaptive RNA editing signals between insects and plants revealed by a new measurement termed haplotype diversity

BackgroundC-to-U RNA editing in plants is believed to confer its evolutionary adaptiveness by reversing unfavorable DNA mutations. This “restorative hypothesis” has not yet been tested genome-wide. In contrast, A-to-I RNA editing in insects like Drosophila and honeybee is already known to benefit the host by increasing proteomic diversity in a spatial-temporal manner (namely “diversifying hypothesis”).MethodsWe profiled the RNA editomes of multiple tissues of Arabidopsis thaliana, Drosophila melanogaster, and Apis melifera. We unprecedentedly defined the haplotype diversity (HD) of RNA molecules based on nonsynonymous editing events (recoding sites).ResultsSignals of adaptation is confirmed in Arabidopsis by observing higher frequencies and levels at nonsynonymous editing sites over synonymous sites. Compared to A-to-I recoding sites in Drosophila, the C-to-U recoding sites in Arabidopsis show significantly lower HD, presumably due to the stronger linkage between C-to-U events.ConclusionsC-to-U RNA editing in Arabidopsis is adaptive but it is not designed for diversifying the proteome like A-to-I editing in Drosophila. Instead, C-to-U recoding sites resemble DNA mutations. Our observation supports the restorative hypothesis of plant C-to-U editing which claims that editing is used for fixing unfavorable genomic sequences.

The Ka/Ks and πa/πs ratios under different models of gametophytic and sporophytic selection.

Alternation of generations in plant life cycle provides a biological basis for natural selection occurring in either the gametophyte or the sporophyte phase, or in both. Divergent biphasic selection could yield distinct evolutionary rates for phase-specific or pleiotropic genes. Here we analyze models that deal with antagonistic and synergistic selection between alternative generations in terms of the ratio of nonsynonymous to synonymous divergence (Ka/Ks). Effects of biphasic selection are opposite under antagonistic selection, but cumulative under synergistic selection for pleiotropic genes. Under the additive and comparable strengths of biphasic allelic selection, the absolute Ka/Ks for the gametophyte-gene is equal to in outcrossing but smaller than in a mixed mating system that for the sporophyte-gene under antagonistic selection. The same pattern is predicted for Ka/Ks under synergistic selection. Selfing reduces efficacy of gametophytic selection. Other processes, including pollen and seed flow and genetic drift, reduce selection efficacy. The polymorphism (πa) at a nonsynonymous site is affected by the joint effects of selfing with gametophytic or sporophytic selection. Likewise, the ratio of nonsynonymous to synonymous polymorphism (πa/πs) is also affected by the same joint effects. Gene flow and genetic drift have opposite effects on πa or πa/πs in interacting with gametophytic and sporophytic selection. We discuss implications of this theory for detecting natural selection in terms of Ka/Ks and for interpreting the evolutionary divergence among gametophyte-specific, sporophyte-specific, and pleiotropic genes.

Relationship between Sugarcane eIF4E Gene and Resistance against Sugarcane Streak Mosaic Virus.

Sugarcane mosaic disease, mainly caused by Sugarcane streak mosaic virus (SCSMV), has serious adverse effects on the yield and quality of sugarcane. Eukaryotic translation initiation factor 4E (eIF4E) is a natural resistance gene in plants. The eIF4E-mediated natural recessive resistance results from non-synonymous mutations of the eIF4E protein. In this study, two sugarcane varieties, CP94-1100 and ROC22, were selected for analysis of their differences in resistance to SCSMV. Four-base missense mutations in the ORF region of eIF4E resulted in different conserved domains. Therefore, the differences in resistance to SCSMV are due to the inherent differences in eIF4E of the sugarcane varieties. The coding regions of eIF4E included 28 SNP loci and no InDel loci, which were affected by negative selection and were relatively conserved. A total of 11 haploids encoded 11 protein sequences. Prediction of the protein spatial structure revealed three non-synonymous mutation sites for amino acids located in the cap pocket of eIF4E; one of these sites existed only in a resistant material (Yuetang 55), whereas the other site existed only in a susceptible material (ROC22), suggesting that these two sites might be related to the resistance to SCSMV. The results provide a strong basis for further analysis of the functional role of eIF4E in regulating mosaic resistance in sugarcane.

Identification and fine mapping of gummy stem blight resistance gene Gsb-7(t) in Melon.

Gummy stem blight (GSB), caused by the Didymella bryoniae (Auersw.) Rehm, is a devastating fungal disease of melon worldwide. Breeding GSB-resistant cultivars with host resistance genes is considered to be the most economic and effective strategy to control this disease. In this study, 260 melon germplasm resources were screened for resistance to GSB, and an inbred line H55R exhibited immunity to GSB was identified. To further understand the resistance mechanism of H55R against GSB, an F2 population was obtained from a cross between the GSB-susceptible line A15 and H55R, and genetic analysis indicated that the GSB resistance in H55R was controlled by a single dominant gene, tentatively named Gsb-7(t). The Gsb-7(t) gene was finally delimited to a 140 kb interval on chromosome 7 using bulked-segregant analysis and chromosome walking strategies. Ten putative genes were annotated in this region that contains a wall-associated receptor kinase (WAK) gene MELO3C010403. The MELO3C010403 gene contains two alternative transcripts, MELO3C010403-T1 and MELO3C010403-T2, with five and seven non-synonymous mutation sites, respectively. Gene expression analysis showed that expression of MELO3C010403-T2 but not the MELO3C010403-T1 was significantly induced by the D. bryoniae at 24 hours post-inoculation (hpi), indicating that the MELO3C010403-T2 transcript of MELO3C010403 was the most likely candidate gene of Gsb-7(t). Our results offer new genetic resources and will be helpful for the development of GSB-resistant melon cultivars in the future.

Weaker selection on genes with treatment-specific expression consistent with a limit on plasticity evolution in Arabidopsis thaliana.

Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene's expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.

Abstract 954: PIK3CA mutation induces immunogenicity and increases the immune checkpoint inhibitor response in urothelial carcinoma

Abstract Background: Immune checkpoint inhibitor (ICI) has been widely used for advanced urothelial carcinoma (mUC) globally. However, little is known about the potential predictive factors of the ICI response. The study aims to identify a potential biomarker to predict ICI response in UC. Material and methods: We performed targeted next generation sequencing (NGS) of pre-treatment tumor Formalin-fixed paraffin-embedded (FFPE) samples by using ACTOnco® comprehensive genomic panel (ACT Genomics, Co. Ltd., Taiwan) to target all coding exons of 440 cancer-associated genes. We retrospective analyzed a total of 70 patients with mUC underwent at least one cycle of ICI with or without chemotherapy at Kaohsiung Chang Gung Memorial Hospital between April 2016 and November 2019. Patients were defined as responders (CR/PR) and non-responders (SD/PD) according to the RECIST evaluation. TMB was estimated by calculating all coding exons variants, including somatic non-synonymous, frameshift and splice site variants. Results: A total of 55 UC tumor samples passed QC process were analyzed. The median age was 67, and 74% of all patients were in ICI monotherapy cohort. The most common genetic alterations were TP53 (66%), KMT2D (42%) and MUC16 (32%). In ICI monotherapy group, responders had a higher TMB than non-responders (median 19.5 vs 6.5 mt/Mb; p = 0.004). Patients with high TMB on ICI monotherapy had a significantly longer PFS (14.6 vs 2.9 months; p = 0.003) and OS (NR vs 11.0 months; p = 0.01) than TMB low patients. We found that PIK3CA (32% vs 3%), UBR5 (23% vs 0%), PRKDC (32% vs 10%) were the top 3 enriched mutations between responders and non-responders. All PIK3CA mutation were hot-spot mutation, and tumor harbored PIK3CA mutation had a median high TMB than PIK3CA wild type in UC, NSCLC, breast and colorectal cancer. By analyzing tumor and immune microenvironment by TISIDB platform, UC with PIK3CA mutation was highly associated with increasing MHC molecules expression (HLA-A, HLA-B, HLA-C, B2M) in UC and gastric cancer. In vitro study, we confirmed anti-PDL1 with PBMC had a synergistic cytotoxicity in PIK3CA mutated cell lines (BFTC-909, AGS and DLD-1), but not in PIK3CA wild type cell lines (T24, RT4, UMUC-14, N87 and SW837). Knocking-down PIK3CA expression diminished the cytotoxic activity of anti-PDL1 with PBMC in BFTC-909 and AGS cells. PIK3CA down-regulation inhibits antigen presentation machinery process. We confirmed that knockout PIK3CA gene function by CRISPR/Cas9 decreased MHC1 and B2M molecules expression, and showing modest resistance to ICI than mice with preserved PIK3CA function. Conclusion: UC with PIK3CA missense mutation exhibits a higher immunogenicity through enhancing MHC class I and B2M expression. PIK3CA missense mutation is a potential predictive marker for ICI response in UC. Citation Format: Harvey Yu-Li Su, Ling-Yi Xiao, Shih-Yu Huang, Li-Chung Chang, Yi-Hua Chen, Hao-Lun Luo, Chia-Ling Wu. PIK3CA mutation induces immunogenicity and increases the immune checkpoint inhibitor response in urothelial carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 954.

Deleterious mutation load in the admixed mice population

Deleterious mutation loads are known to correlate negatively with effective population size (Ne). Due to this reason, previous studies observed a higher proportion of harmful mutations in small populations than that in large populations. However, the mutational load in an admixed population that derived from introgression between individuals from two populations with vastly different Ne is not known. We investigated this using the whole genome data from two subspecies of the mouse (Mus musculus castaneus and Mus musculus musculus) with significantly different Ne. We used the ratio of diversities at nonsynonymous and synonymous sites (dN/dS) to measure the harmful mutation load. Our results showed that this ratio observed for the admixed population was intermediate between those of the parental populations. The dN/dS ratio of the hybrid population was significantly higher than that of M. m. castaneus but lower than that of M. m. musculus. Our analysis revealed a significant positive correlation between the proportion of M. m. musculus ancestry in admixed individuals and their dN/dS ratio. This suggests that the admixed individuals with high proportions of M. m. musculus ancestry have large dN/dS ratios. We also used the proportion of deleterious nonsynonymous SNVs as a proxy for deleterious mutation load, which also produced similar results. The observed results were in concordance with those expected by theory. We also show a shift in the distribution of fitness effects of nonsynonymous SNVs in the admixed genomes compared to the parental populations. These findings suggest that the deleterious mutation load of the admixed population is determined by the proportion of the ancestries of the subspecies. Therefore, it is important to consider the status and the level of genetic admixture of the populations whilst estimating the mutation loads.