Nucleotide Substitutions Research Articles

Vanderwaltozyma japonica sp. nov.: a novel yeast isolated from subtropical regions in Japan.

Seven yeast strains, NBRC 11072, NBRC 11076, NBRC 11078T, NBRC 11080, NBRC 116645, NBRC 116646 and NBRC 116647, were isolated from the subtropical regions of Japan. The seven strains share identical nucleotide sequences in the D1/D2 domains of the large subunit (LSU) rRNA gene and internal transcribed spacer (ITS) regions. Six strains (NBRC 11072, NBRC 11076, NBRC 11078T, NBRC 11080, NBRC 116645 and NBRC 116646) share identical mitochondrial cytochrome c oxidase subunit II (COX II) sequences as well, while NBRC 116647 differs from them by one substitution in this region. The yeast strains were closely related to Vanderwaltozyma verrucispora NBRC 1884T in the ITS regions. However, in the D1/D2 domains and COX II, they were closely related to Torulaspora globosa CBS 764T and Naumovozyma castellii NRRL Y-12630T, respectively. They differed from V. verrucispora by 49 nucleotide substitutions with 26 gaps in the ITS region, from T. globosa by 19 nucleotide substitutions with 3 gaps in the D1/D2 domains and from N. castellii by 49 nucleotide substitutions with 6 gaps in the COX II. The phylogenetic analysis of the seven strains based on the sequences of the D1/D2 domains of LSU and ITS regions among related species showed that these strains represent a single novel species within the genus Vanderwaltozyma. The seven strains could be differentiated from other known species of the genus Vanderwaltozyma because of their inability to assimilate glycerol as a carbon source. Therefore, we suggest that these seven strains represent a novel species, designated Vanderwaltozyma japonica sp. nov. The holotype is NBRC 11078T (isotype: TBRC 19102T). The MycoBank number is MB 855476.

Abstract B057: Base-error model (BEM) for improved detection of ctDNA in lymph samples of HPV-negative head and neck cancer patients

Abstract Introduction: Previously we demonstrated that ctDNA present in lymphatic exudate collected via surgical drains (“lymph”) outperformed plasma for detecting minimal residual disease (MRD) in head and neck squamous cell carcinoma (HNSCC) patients through a targeted sequencing (TS) approach. However, artifacts introduced during library preparation and sequencing remain challenges to sensitive low allelic fraction mutation detection. To enhance the performance of the TS approach, we characterized the background noise in lymph samples and developed a base-error model (BEM) that reduces sequencing artifacts and maximizes ctDNA signal. Methods: Lymph and blood were collected from 44 HPV-negative HNSCC patients postoperatively at 24 hours along with resected tumor. Cell-free DNA was extracted from lymph and sequenced using a custom TS panel with unique molecular identifiers (UMI). Somatic mutations were identified by tumor and blood exome sequencing (200x). Five patients had <2 mutations in tumor and were excluded. Two patients were censored due to lack of clinical data, yielding 17 patients with disease recurrence (REC) and 20 with no evidence of disease (NED) with >1 year of follow-up. Tumor-specific variants were force-called in lymph using a custom pipeline. A BEM of each tumor variant was built using a series of high-quality lymph reference samples to quantify the background noise. BEM was fit by Weibull distribution if more than 5 non-zero non-reference alleles were observed at a tumor variant position across the reference samples. Otherwise, Gaussian distribution was used to fit BEM. Force-called variants were considered artifacts if the variant allele fraction (VAF) was not greater than BEM cutoff controlled by false discovery rate. Patients were considered MRD positive if the mean VAF was greater than the estimated limit of detection. The Kaplan-Meier (KM) estimator with log-rank test and Cox proportional-hazards model were used for survival analyses. Results: Background noise was observed in reference lymph samples across all 12 nucleotide substitution classes. Out of 299 single nucleotide variants (SNVs), 94 were determined as artifacts by BEM. Among the artifacts, 77 were C/G to A/T mutations, indicating that the majority of recurrent background artifact in TS is likely the result of oxidative damage. The artifacts had VAFs ranging from 0.004% to 0.12% with median of 0.02%. Incorporation of BEM in lymph TS showed significantly improved specificity (sensitivity (SN) = 71%, specificity (SP) = 65%; p = 0.027, Hazard ratio (HR) = 3.07) compared to the same cohort without BEM (SN = 76%, SP = 30%, p = 0.60, HR = 1.34). Conclusions: BEM is a sensitive and robust method to improve ctDNA detection accuracy in TS. We demonstrated its utility in improving recurrence prediction in a HNSCC cohort. Given its advantages, we envision that BEM will prove useful as a general strategy for deep sequencing applications requiring precise digital quantification of low frequency alleles in TS. Citation Format: Zhuosheng Gu, Maciej Pacula, Seka Lazare, Noah Earland, Marra S. Francis, Adam Harmon, Megan Long, Aadel A. Chaudhuri, Jose P. Zevallos, Wendy Winckler. Base-error model (BEM) for improved detection of ctDNA in lymph samples of HPV-negative head and neck cancer patients [abstract]. In: Proceedings of the AACR Special Conference: Liquid Biopsy: From Discovery to Clinical Implementation; 2024 Nov 13-16; San Diego, CA. Philadelphia (PA): AACR; Clin Cancer Res 2024;30(21_Suppl):Abstract nr B057.

PATH-54. HIGH-THROUGHPUT SINGLE NUCLEAR DNA SEQUENCING OF HUMAN SPORADICNF1 MUTANT IDH-WILDTYPE GLIOBLASTOMAS REVEALS PATTERNS OF TUMOR EVOLUTION AND RECURRENCE

Abstract The advent of single cell techniques has advanced our understanding of glioblastoma (GBM) evolution and cell lineage relationships. However, study of mutational co-occurrence and clonal phylogeny has been hampered by limitations in single cell genotyping. Here, we perform semi-automated, rapid single nuclear dissociation followed by targeted single nucleus DNA sequencing (snDNA-seq) of eleven retrospectively identified somatic NF1 mutant IDH-wildtype glioblastomas. Bulk DNA sequencing was performed as part of routine clinical care using a CLIA-certified targeted sequencing panel. For snDNA-seq, tissue cores with greater than 30% tumor were dissociated on a S2 genomics S100 Singulator followed by snDNA-seq using a targeted 361 amplicon panel (Tapestri, Mission Bio, USA) to sequence for recurrently observed single nucleotide substitutions and copy number alterations. A mean of 3,661 cells were recovered per sample with a mean of 96 reads per cell per amplicon and 87% mean panel uniformity. SnDNA-seq validated point mutations and copy number alterations observed in bulk DNA-sequencing and revealed novel alterations and subclonal copy number alterations not detected on bulk analysis. Within individual samples, snDNA-seq based phylogenetic reconstruction identified mutually exclusive patterns of mutation between NF1 and PI3K signaling alterations, suggesting these events may occur in distinct tumor subclones. With regard to copy number alterations (CNAs), snDNA-seq demonstrated intra-tumoral heterogeneity for chromosome 9p loss, chromosome 7 gain, and chromosome 10 monosomy. Subclonal CNA clones not detected by bulk sequencing were identified in 5 samples. Analysis of a matched primary-recurrent tumor pair revealed expansion of a specific mutational clone (loss of 9p, 10p, gains of 1p, 7, 15, and 19q) at recurrence. SnDNA-seq recapitulated bulk DNA sequencing alterations and identified distinct patterns of mutational co-occurrence, CNAs, and tumor evolution over time. Reconstructing GBM phylogeny at single cell resolution has potential translational implications for understanding tumor evolution and treatment resistance.

Genetic diversity and spread of recombinant coxsackievirus A4 in hand, foot, and mouth disease cases in Bangkok, Thailand: 2017–2023

Coxsackievirus A4 (CVA4) has recently become one of the most common causative agents of hand, foot, and mouth disease. The current study investigated the genetic diversity and spread of recombinant CVA4 by analyzing circulating genotypes and recombinant strains in Bangkok, Thailand, from 2017 to 2023. Partial VP1, 3Dpol, and whole genome sequencing of CVA4 samples collected from collaborating hospitals were conducted. Phylogenetic analysis of CVA4 VP1 and 3Dpol genome regions revealed discordance, indicating recombination. The predominant CVA4 genotype was C3, primarily observed in 2019. The predominant genotype in 2017 was C1. D2, commonly found in China, was occasionally observed. In nucleotide similarity analysis, intertypic recombination between CVA4 and EV-A during the evolutionary history of the virus was evident, particularly in the nonstructural region. The estimated emergence of genotypes C1 and C3 in Thailand occurred around 2014, with an evolutionary rate of 5.8 × 10− 3 nucleotide substitutions per site per year. Genotype D2 exhibited notable variability across both the entire genome and the structural protein region compared to genotype C. Monitoring the genetic diversity and circulation of recombinant CVA4 is crucial for identifying newly emerging virus strains, enabling prompt public health responses and containment efforts, and enhancing surveillance in Thailand.

Novel HLA-DQB1, DQA1, DPA1 and DPB1 alleles identified in the Australian New South Wales tissue typing laboratory

Introduction of NGS into clinical histocompatibility laboratories has greatly increased the frequency of novel HLA allele discovery. We identified 9 DQA1, 7 DQB1, 8 DPA1 and 2 DPB1 novel alleles over the last 2 years. 92% (24 of 26) differed from their closest allele by single nucleotide substitutions detected in coding regions and 84% (22 of 26) contained polymorphism outside of exon 2. Identification of novel alleles in non-coding intronic regions were restricted to splice sites. Two novel alleles arising from nucleotide substitution affecting splicing were identified and named with questionable expression status. Novel alleles with polymorphism detected only in non-coding regions (introns and UTRs) were not reported. 21 additional DQA1, DQB1 and DPA1 alleles detected in multiple individuals as novel at the time of testing, were also detected in other laboratories. The number of novel alleles detected further emphasises the role of NGS in novel allele discovery.

Characterisation of the Novel HLA-A*02:1189 Allele by Sequencing-Based Typing.

HLA-A2*02:1189 differs from HLA-A*02:01:01:01 by one nucleotide substitution in codon 275 in exon 5.

Discovery of the Novel HLA-DRB1*16:02:13 Allele, a Variant of HLA-DRB1*16, in a Taiwanese Individual.

One nucleotide substitution in codon 206 of HLA-DRB1*16:02:01:01 results in a novel allele, HLA-DRB1*16:02:13.

Characterisation of the Novel HLA-B*40:587 Allele by Sequencing-Based Typing.

HLA-B*40:587 differs from HLA-B*40:02:01:01 by one nucleotide substitution in codon 275 in exon 5.

The Novel HLA-DRB1*04:388 Allele, Identified by Sanger Dideoxy Nucleotide Sequencing in a Chinese Individual.

The HLA-DRB1*04:388 allele differs from HLA-DRB1*04:01:01:01 by one nucleotide substitution in codon 44 in exon 2.

Genotypic Antimicrobial Resistance Profiles of Diarrheagenic Escherichia coli and Nontyphoidal Salmonella Strains Isolated from Children with Diarrhea and Their Exposure Environments in Ethiopia.

Antimicrobial resistance (AMR) poses a significant global threat, particularly in low- and middle-income countries, such as Ethiopia, where surveillance is limited. This study aimed to predict and characterize the AMR profiles of diarrheagenic Escherichia coli (DEC) and nontyphoidal Salmonella (NTS) strains isolated from human, animal, food, and environmental samples using whole genome sequencing. A total of 57 NTS and 50 DEC isolates were sequenced on an Illumina NextSeq 550. The ResFinder and PointFinder tools were employed to identify antimicrobial resistance genes (ARGs) and point mutations. Salmonella serotypes were determined using SeqSero. The analysis identified at least one ARG in every NTS sample and 78% of the DEC isolates, with 22 distinct ARGs in the NTS samples and 40 in the DEC samples. The most prevalent ARGs were aac(6')-Iaa and aph(3')-Ib, which predict aminoglycoside resistance in 100% of NTS and 54% of DEC isolates, respectively. Other commonly identified ARGs include sul2, aph(6)-Id, blaTEM-1B , and tet(A), which confer resistance to folate inhibitors, aminoglycosides, β-lactams, and tetracycline. Some ARGs predicted phenotypic multidrug resistance in both DEC and NTS isolates. All identified β-lactam ARGs, except for blaTEM -1D, conferred resistance to more than three antibiotics. Interestingly, blaCTX- M-15 was found to confer resistance to nine antibiotics, including third-generation cephalosporins, in 18% of DEC and 3.5% of NTS isolates. DEC isolates from children exhibited the highest ARG diversity. Notably, genes such as aph(3″)-Ib, aph(6)-Id, sul2, and tet(A) were detected across all sample types, including water sources, although some ARGs were exclusive to specific sample types. Point mutations mediating AMR were detected in several genes, with mutations associated with nucleotide substitution being the most frequent. This genotypic AMR profiling revealed the presence of widespread drug-resistant NTS and DEC strains in Ethiopia. Robust and sustained AMR surveillance is essential for monitoring the emergence and spread of these resistant pathogens.

The Novel HLA-B*40:585 Allele, Identified by Sanger Dideoxy Nucleotide Sequencing in a Chinese Individual.

The HLA-B*40:585 allele differs from HLA-B*40:01:02:01 by one nucleotide substitution in codon 130 in exon 3.

Field-evolved resistance to neonicotinoids in the mosquito, Anopheles gambiae, is associated with mutations of nicotinic acetylcholine receptor subunits combined with cytochrome P450-mediated detoxification

New insecticides prequalified for malaria control interventions include modulators of nicotinic acetylcholine receptors that act selectively on different subunits leading to variable sensitivity among arthropods. This study aimed to investigate the molecular mechanisms underlying contrasting susceptibility to neonicotinoids observed in wild populations of two mosquito sibling species. Bioassays and a synergist test with piperonyl butoxide revealed that the sister taxa, Anopheles gambiae and An. coluzzii, from Yaounde, Cameroon, both have the potential to develop resistance to acetamiprid through cytochrome P450-mediated detoxification. However, contrary to An. coluzzii, An. gambiae populations are evolving cross-resistance to several active ingredients facilitated by mutations of nicotinic acetylcholine receptors (nAChRs). We sequenced coding regions on the β1 and α6 nAChR subunits where variants associated with resistance to neonicotinoids or to spinosyns have been found in agricultural pests and detected no mutation in An. coluzzii. By contrast, six nucleotide substitutions including an amino acid change in one of the loops that modulate ligand binding and affect sensitivity were present in the resistant species, An. gambiae. Allele frequency distributions were consistent with the spread of beneficial mutations that likely reduce the affinity of An. gambiae nAChRs for synthetic modulators. Our findings provide critical information for the application and resistance management of nAChR modulators in malaria prevention.

Number of human protein interactions correlates with structural, but not regulatory conservation of the respective genes.

The differential ratio of nonsynonymous to synonymous nucleotide substitutions (dN/dS) is a common measure of the rate of structural evolution in proteincoding genes. In addition, we recently suggested that the proportion of transposable elements in gene promoters that host functional genomic sites serves as a marker of the rate of regulatory evolution of genes. Such functional genomic regions may include transcription factor binding sites and modified histone binding loci. Here, we constructed a model of the human interactome based on 600,136 documented molecular interactions and investigated the overall relationship between the number of interactions of each protein and the rate of structural and regulatory evolution of the corresponding genes. By evaluating a total of 4,505 human genes and 1,936 molecular pathways we found a general correlation between structural and regulatory evolution rate metrics (Spearman 0.08-0.16 and 0.25-0.37 for gene and pathway levels, respectively, p < 0.01). Further exploration revealed in the established human interactome model lack of correlation between the rate of gene regulatory evolution and the number of protein interactions on gene level, and weak negative correlation (∼0.15) on pathway level. We also found a statistically significant negative correlation between the rate of gene structural evolution and the number of protein interactions (Spearman -0.11 and -0.3 for gene and pathway levels, respectively, p < 0.01). Our result suggests stronger structural rather than regulatory conservation of genes whose protein products have multiple interaction partners.

Genomic epidemiology reveals the variation and transmission properties of SARS-CoV-2 in a single-source community outbreak.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the coronavirus disease 2019 (COVID-19) pandemic, which is still a global public health concern. During March 2022, a rapid and confined single-source outbreak of SARS-CoV-2 was identified in a community in Nanjing municipal city. Overall, 95 individuals had laboratory-confirmed SARS-CoV-2 infection. The whole genomes of 61 viral samples were obtained, which were all members of the BA.2.2 lineage and clearly demonstrated the presence of one large clade, and all the infections could be traced back to the original index case. The most distant sequence from the index case presented a difference of 4 SNPs, and 118 intrahost single-nucleotide variants (iSNVs) at 74 genomic sites were identified. Some minor iSNVs can be transmitted and subsequently rapidly fixed in the viral population. The minor iSNVs transmission resulted in at least two nucleotide substitutions among all seven SNPs identified in the outbreak, generating genetically diverse populations. We estimated the overall transmission bottleneck size to be 3 using 11 convincing donor-recipient transmission pairs. Our study provides new insights into genomic epidemiology and viral transmission, revealing how iSNVs become fixed in local clusters, followed by viral transmission across the community, which contributes to population diversity.

Effect of C-to-T transition at CpG sites on tumor suppressor genes in tumor development in cattle evaluated by somatic mutation analysis in enzootic bovine leukosis.

Oncogenic transformation of normal cells is caused by mutations and chromosomal abnormalities in cancer-related genes. Enzootic bovine leukosis (EBL) is a malignant B-cell lymphoma caused by bovine leukemia virus (BLV) infection in cattle. Although a small fraction of BLV-infected cattle develops EBL after a long latent period, the mechanisms for oncogenesis in EBL cattle remain largely unknown. In this study, we analyzed the types and patterns of somatic mutations in cancer cells from 36 EBL cases, targeting 21 cancer-related genes. Various somatic mutations were identified in eight genes, TP53, KMT2D, CREBBP, KRAS, PTEN, NOTCH1, MYD88, and CARD11. In addition, TP53 gene was found to be mutated in 69.4% of EBL cases, with most being biallelic mutations. In some cases, associations were observed between the ages at which cattle had developed EBL and somatic mutation patterns; young onset of EBL possibly occurs due to high impact mutations affecting protein translation and biallelic mutations. Furthermore, nucleotide substitution patterns indicated that cytosine at CpG sites tended to be converted to thymine in many EBL cases, which was considered to be the result of spontaneous deamination of 5-methylcytosine. These results demonstrate how somatic mutations have occurred in cancer cells leading to EBL development, thereby explaining its pathogenic mechanism. These findings will contribute to a better understanding and future elucidation of disease progression in BLV infection.IMPORTANCEEnzootic bovine leukosis (EBL) is a malignant and lethal disease in cattle. Currently, there are no effective vaccines or therapeutic methods against bovine leukemia virus (BLV) infection, resulting in severe economic losses in livestock industry. This study provides a renewed hypothesis to explain the general mechanisms of EBL onset by combining the previous finding that several integration sites of BLV provirus can affect the increase in survival and proliferation of infected cells. We demonstrate that two additional random events are necessary for oncogenic transformation in infected cell clones, elucidating the reason why only few infected cattle develop EBL. Further exploration of somatic mutation and BLV integration sites could support this hypothesis more firmly, potentially contributing to the development of novel control methods for EBL onset.

Copy Number Analysis of Whole Exome Sequencing Data

Abstract Background Genetic alterations implicated in constitutional disorders range in size from single nucleotide substitutions to losses and gains of millions of bases. Whole exome sequencing (WES) captures many clinically relevant variants, but an analytic pipeline optimized for detection of small substitutions, insertions, and deletions may miss pathogenic copy number variants (CNVs). Historically, our laboratory performed separate cytogenetic tests to evaluate for larger genomic imbalances. Addition of copy number analysis to our WES pipeline could improve the assay’s diagnostic yield. This study describes part of our strategy to validate copy number analysis of our constitutional whole exome sequencing data. Methods We identified 18 CNV specimens that underwent both chromosomal microarray analysis (CMA) and WES between 2017 and 2023. Variants included 13 deletions (size range 14Kb-2Mb) and 5 duplications (size range 325Kb-10Mb). The standard method of constitutional copy number analysis in our clinical laboratory is CMA using the Affymetrix Cytoscan HD Array and analysis software. CNVs identified on CMA were reviewed with reference to the Database of Genomic Variants and multiple academic hospital databases. WES was performed on an Illumina sequencing system following capture-based library enrichment with the Integrated DNA Technologies xGen Exome Research, CNV Backbone, Human mtDNA Research, and Human ID Research Panel probes. CNVs identified on WES were called using the Broad Institute’s Genomic Analysis Toolkit (GATK) with a quality score cutoff of 120. Whole exome sequencing data was visualized on the Integrative Genomics Viewer. Results 14 out of 18 CNVs reported on CMA were also detected on WES with a quality score at or above 120 (78%). This includes 10 out of 13 deletions (77%) and 3 out of 5 duplications (60%). For both deletions and duplications, the size ranges of detected variants at or above the cutoff score and undetected variants overlapped. One duplication reported on CMA was detected on WES with a quality score below 120. Visualization of whole exome sequencing data showed variation in probe density that may have affected CNV calling. Conclusions The majority of CNVs reported after CMA (14/18, 78%) were also detected on WES data using a GATK quality score cutoff of 120. The data does not suggest an association between CNV size and detection by WES, although interpretation is limited by the small number of cases. Data visualization may be required to identify some CNVs. Future steps to characterize the limitations of copy number analysis of WES data include assessments of variance in probe density and read depth across the exome. Copy number analysis of WES data from cases with normal CMA results will also be performed for assay validation.

Inference of Phylogenetic Networks from Sequence Data using Composite Likelihood.

While phylogenies have been essential in understanding how species evolve, they do not adequately describe some evolutionary processes. For instance, hybridization, a common phenomenon where interbreeding between two species leads to formation of a new species, must be depicted by a phylogenetic network, a structure that modifies a phylogenetic tree by allowing two branches to merge into one, resulting in reticulation. However, existing methods for estimating networks become computationally expensive as the dataset size and/or topological complexity increase. The lack of methods for scalable inference hampers phylogenetic networks from being widely used in practice, despite accumulating evidence that hybridization occurs frequently in nature. Here, we propose a novel method, PhyNEST (Phylogenetic Network Estimation using SiTe patterns), that estimates binary, level-1 phylogenetic networks with a fixed, user-specified number of reticulations directly from sequence data. By using the composite likelihood as the basis for inference, PhyNEST is able to use the full genomic data in a computationally tractable manner, eliminating the need to summarize the data as a set of gene trees prior to network estimation. To search network space, PhyNEST implements both hill climbing and simulated annealing algorithms. PhyNEST assumes that the data are composed of coalescent independent sites that evolve according to the Jukes-Cantor substitution model and that the network has a constant effective population size. Simulation studies demonstrate that PhyNEST is often more accurate than two existing composite likelihood summary methods (SNaQ and PhyloNet) and that it is robust to at least one form of model misspecification (assuming a less complex nucleotide substitution model than the true generating model). We applied PhyNEST to reconstruct the evolutionary relationships among Heliconius butterflies and Papionini primates, characterized by hybrid speciation and widespread introgression, respectively. PhyNEST is implemented in an open-source Julia package and is publicly available at https://github.com/sungsik-kong/PhyNEST.jl.

Changes in the genome of tick-borne encephalitis virus during cultivation

The tick-borne encephalitis virus (TBEV) of the Siberian genotype was previously isolated from the brain of a deceased person. TBEV variants obtained at passages 3 and 8 on SPEV cells were inoculated into the brains of white mice for subsequent passages. Full-genome sequences of all virus variants were analyzed by high-throughput sequencing. The analysis showed the presence of 41 point nucleotide substitutions, which were mainly localized in the genes of non-structural proteins NS3 and NS5 of TBEV. In the deduced virus protein sequences, 12 amino acid substitutions were identified. After three passages through mouse brains, reverse nucleotide and amino acid substitutions were detected. Most of them were mapped in the NS5 protein gene, where 5 new nucleotide substitutions also appeared. At the same time, there was an increase in the length of the 3′-untranslated region (3′-UTR) of the viral genome by 306 nucleotides. The Y3 and Y2 3’-UTR elements were found to contain imperfect L and R repeats, which probably associated with inhibition of the activity of cellular XRN1 RNase and thus involved in the formation of sfRNA1 and sfRNA2. For all TBEV variants, high levels of infectious virus were detected both in cell culture and in the brain of white mice. The revealed changes in the genome that occur during successive passages of TBEV are most likely due to the significant genetic variability of the virus, which ensures its efficient reproduction in different hosts and active circulation in nature.

Natural selection on apical membrane antigen 1 (AMA1) of an emerging zoonotic malaria parasite Plasmodium inui

Apical membrane antigen 1 (AMA1) of malaria parasites plays an important role in host cell invasion. Antibodies to AMA1 can inhibit malaria merozoite invasion of erythrocytes while vaccine-induced specific cytotoxic T cell responses to this protein are associated with clinical protection. Polymorphisms in AMA1 of Plasmodium falciparum (PfAMA1) and P. vivax (PvAMA1) are of concern for vaccine development. To date, little is known about sequence diversity in ama1 of P. inui (Piama1), an emerging zoonotic malaria parasite. In this study, 80 complete Piama1 coding sequences were obtained from 57 macaques in Thailand that defined 60 haplotypes clustering in two phylogenetic lineages. In total, 74 nucleotide substitutions were identified and distributed unevenly across the gene. Blockwise analysis of the rates of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions did not show a significant deviation from neutrality among Thai isolates. However, significantly negative Tajima’s D values were detected in domain I and the loop region of domain II, implying purifying selection. Codon-based analysis of dN/dS has identified 12 and 14 codons under positive and negative selections, respectively. Meanwhile, 85 amino acid substitutions were identified among 80 Thai and 11 non-Thai PiAMA1 sequences. Of these, 48 substituted residues had a significant alteration in physicochemical properties, suggesting positive selection. More than half of these positively selected amino acids (32 of 48) corresponded to the predicted B-cell or T-cell epitopes, suggesting that selective pressure could be mediated by host immunity. Importantly, 14 amino acid substitutions were singletons and predicted to be deleterious that could be subject to ongoing purifying selection or elimination. Besides genetic drift and natural selection, intragenic recombination identified in domain II could generate sequence variation in Piama1. It is likely that malarial ama1 exhibits interspecies differences in evolutionary histories. Knowledge of the sequence diversity of the Piama1 locus further provides an evolutionary perspective of this important malaria vaccine candidate.

Biology and genetic diversity of Candida krusei isolates from fermented vegetables and clinical samples in China

ABSTRACT Candida krusei, also known as Pichia kudriavzevii, is an emerging non-albicans Candida (NAC) species causing both superficial and deep-seated infections in humans. This fungal pathogen is inherently resistant to the first-line antifungal drug, fluconazole, and is widely distributed in natural environments such as soil, foods, vegetables, and fruits. In this study, we collected 86 C. krusei strains from clinical settings and traditional fermented vegetables from different areas of China. Compared to C. krusei strains from fermented vegetables, clinical isolates exhibited a higher ability to undergo filamentation and biofilm development, which could facilitate its host colonization and infections. Isolates from fermented vegetables showed higher resistance to several antifungal drugs including fluconazole, voriconazole, itraconazole, amphotericin B, and caspofungin, than clinical strains, while they were more susceptible to posaconazole than clinical strains. Although C. krusei has been thought to be a diploid organism, we found that one-fourth of clinical strains and the majority of isolates from fermented vegetables (87.5%) are triploid. Whole-genome sequencing and population genetic analyses demonstrated that isolates from clinical settings and fermented food are genetically associated, and distributed across a wide range of genetic clusters. Additionally, we found that six nucleotide substitutions at the promoter region of the ABC11 gene, encoding a multidrug efflux pump, could play a critical role in antifungal resistance in this species. Given the ubiquitous distribution of C. krusei strains in fermented vegetables and their genetic association with clinical strains, a One Health approach will be necessary to control the prevalence of this pathogen.