Synonymous Nucleotide Substitutions Research Articles

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.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity.

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Synonymous variant at the terminal nucleotide in exon 3 of F7 causes abnormal splicing: A case report.

Synonymous variants are non-pathogenic due to non-substitution of amino acids. However, synonymous exonic terminal nucleotide substitutions may affect splicing. Splicing variants are easily analyzed at RNA level for genes expressed in blood cells. Minigene analysis provides another method for splicing variant analysis of genes that are poorly or not expressed in peripheral blood. Whole exome sequencing was performed to screen for potential pathogenic mutations in the proband, which were validated within the family by Sanger sequencing. The pathogenicity of the synonymous mutation was analyzed using the minigene technology. The proband harbored the compound heterogeneous variants c. [291G >A; 572-50C >T] and c.681 + 1G >T in F7, of which the synonymous variant c.291G >A was located at the terminal position of exon 3. Minigene analysis revealed exon3 skipping due to this mutation, which may have subsequently affected protein sequence, structure, and function. Our finding confirmed the pathogenicity of c.291G >A, thus extending the pathogenic mutation spectrum of F7, and providing insights for effective reproductive counseling.

Molecular characterization, function, tissue differential expression, and single-nucleotide polymorphism of buffalo TP53 gene

Abstract. TP53 has been shown to be involved in lactation in cattle. However, the role of TP53 in buffalo lactation remains unknown. To this end, we isolated and identified the complete coding sequence (CDS) of the TP53 gene from the buffalo mammary gland and further analyzed its molecular characteristics, function, tissue differential expression, and single-nucleotide polymorphism (SNP). A transcript of this gene was cloned with a CDS length of 1161 bp, encoding a protein consisting of 386 amino acid residues. Bioinformatics analysis showed that buffalo TP53 CDS and the physicochemical characteristics, conserved domains, structure, and function of its encoded protein are highly similar to those of other species in Bovidae. The buffalo TP53 protein contains an N-terminal activation domain, a DNA-binding domain, and a tetrameric domain, and it plays a functional role in the nucleus. TP53 was found to express in all 11 detected buffalo tissues, and its expression in the heart, kidney, brain, muscle, and rumen during lactation was significantly higher than that during non-lactation (p<0.05), while in the liver, lung, and mammary gland, its expression was the opposite (p<0.05). Interference experiments in buffalo mammary epithelial cells (BuMECs) showed that TP53 inhibits the expression of genes related to milk protein and milk fat synthesis through the PI3K–AKT–mTOR pathway. A synonymous nucleotide substitution (c.204C > T) was found in the TP53 CDS of river buffalo, which is the CC homozygote in swamp buffalo. The results indicate that the TP53 gene is involved in buffalo lactation by negatively regulating the synthesis of milk protein and milk fat.

The novel HLA-C*07:02:147 allele characterised by two different sequencing-based typing techniques.

The novel allele HLA-C*07:02:147 differs from HLA-C*07:02:01:01 by one synonymous nucleotide substitution in exon 2.

The novel HLA-A*30:01:23 allele characterised by two different sequencing-based typing techniques.

The novel allele HLA-A*30:01:23 differs from HLA-A*30:01:01:01 by one synonymous nucleotide substitution in exon 2.

The novel HLA-B*44:48:02 allele characterised by two different sequencing-based typing techniques.

The novel allele HLA-B*44:48:02 differs from HLA-B*44:48:01 by one synonymous nucleotide substitution in exon 3.

Genomic full-length sequence of the novel HLA-B*56:94 allele.

HLA-B*56:94 differs from HLA-B*56:05:01 by two non-synonymous nucleotide substitutions in exon 1 and synonymous nucleotide substitutions in exon 1 and exon 2.

The HLA-B*58:01:43 allele identified in a Korean individual awaiting hematopoietic stem cell transplant.

HLA-B*58:01:43 differs from HLA-B*58:01:01:01 by one synonymous nucleotide substitution in codon 197.

The novel HLA-C allele, HLA-C*01:02:84 allele, identified in a solid organ transplantation donor.

HLA-C*01:02:84 differs from HLA-C*01:02:01:01 by one synonymous nucleotide substitution in codon 48.

Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.

The novel HLA-B*53:01:30 allele characterised by two different sequencing-based typing techniques.

The novel allele HLA-B*53:01:30 differs from HLA-B*53:01:01:01 by one synonymous nucleotide substitution in exon 3.

The novel HLA-B*14:118 allele characterized by two different sequencing-based typing techniques.

HLA-B*14:118 differs from HLA-B*14:93 by two non-synonymous and one synonymous nucleotide substitution in exon 2.

Distinct Traits of Structural and Regulatory Evolutional Conservation of Human Genes with Specific Focus on Major Cancer Molecular Pathways.

The evolution of protein-coding genes has both structural and regulatory components. The first can be assessed by measuring the ratio of non-synonymous to synonymous nucleotide substitutions. The second component can be measured as the normalized proportion of transposable elements that are used as regulatory elements. For the first time, we characterized in parallel the regulatory and structural evolutionary profiles for 10,890 human genes and 2972 molecular pathways. We observed a ~0.1 correlation between the structural and regulatory metrics at the gene level, which appeared much higher (~0.4) at the pathway level. We deposited the data in the publicly available database RetroSpect. We also analyzed the evolutionary dynamics of six cancer pathways of two major axes: Notch/WNT/Hedgehog and AKT/mTOR/EGFR. The Hedgehog pathway had both components slower, whereas the Akt pathway had clearly accelerated structural evolution. In particular, the major hub nodes Akt and beta-catenin showed both components strongly decreased, whereas two major regulators of Akt TCL1 and CTMP had outstandingly high evolutionary rates. We also noticed structural conservation of serine/threonine kinases and the genes related to guanosine metabolism in cancer signaling: GPCRs, G proteins, and small regulatory GTPases (Src, Rac, Ras); however, this was compensated by the accelerated regulatory evolution.

Characterization of the novel HLA-DQB1*06:03:47 allele by next generation sequencing.

HLA-DQB1*06:03:47 differs from HLA-DQB1*06:03:01 by one synonymous nucleotide substitution in codon 158 in exon 3.

Geographical origin of Plasmodium vivax in the Hainan Island, China: insights from mitochondrial genome

BackgroundHainan Province, China, has been an endemic region with high transmission of Plasmodium falciparum and Plasmodium vivax. Indigenous malaria caused by P. vivax was eliminated in Hainan in 2011, while imported vivax malaria remains. However, the geographical origin of P. vivax cases in Hainan remains unclear.MethodsIndigenous and imported P. vivax isolates (n = 45) were collected from Hainan Province, and the 6 kb mitochondrial genome was obtained. Nucleotide (π) and haplotype (h) diversity were estimated using DnaSP. The numbers of synonymous nucleotide substitutions per synonymous site (dS) and nonsynonymous nucleotide substitutions per nonsynonymous site (dN) were calculated using the SNAP program. Arlequin software was used to estimate the genetic diversity index and assess population differentiation. Bayesian phylogenetic analysis of P. vivax was performed using MrBayes. A haplotype network was generated using the NETWORK program.ResultsIn total, 983 complete mitochondrial genome sequences were collected, including 45 from this study and 938 publicly available from the NCBI. Thirty-three SNPs were identified, and 18 haplotypes were defined. The haplotype (0.834) and nucleotide (0.00061) diversity in the Hainan populations were higher than China’s Anhui and Guizhou population, and the majority of pairwise FST values in Hainan exceeded 0.25, suggesting strong differentiation among most populations except in Southeast Asia. Most Hainan haplotypes were connected to South/East Asian and China’s others haplotypes, but less connected with populations from China's Anhui and Guizhou provinces. Mitochondrial lineages of Hainan P. vivax belonged to clade 1 of four well-supported clades in a phylogenetic tree, most haplotypes of indigenous cases formed a subclade of clade 1, and the origin of seven imported cases (50%) could be inferred from the phylogenetic tree, but five imported cases (42.8%) could not be traced using the phylogenetic tree alone, necessitating epidemiological investigation.ConclusionsIndigenous cases in Hainan display high genetic (haplotype and nucleotide) diversity. Haplotype network analysis also revealed most haplotypes in Hainan were connected to the Southeast Asian populations and divergence to a cluster of China’s other populations. According to the mtDNA phylogenetic tree, some haplotypes were shared between geographic populations, and some haplotypes have formed lineages. Multiple tests are needed to further explore the origin and expansion of P. vivax populations.

Parallel Evolution of Sex-Linked Genes across XX/XY and ZZ/ZW Sex Chromosome Systems in the Frog Glandirana rugosa.

Genetic sex-determination features male (XX/XY) or female heterogamety (ZZ/ZW). To identify similarities and differences in the molecular evolution of sex-linked genes between these systems, we directly compared the sex chromosome systems existing in the frog Glandirana rugosa. The heteromorphic X/Y and Z/W sex chromosomes were derived from chromosomes 7 (2n = 26). RNA-Seq, de novo assembly, and BLASTP analyses identified 766 sex-linked genes. These genes were classified into three different clusters (XW/YZ, XY/ZW, and XZ/YW) based on sequence identities between the chromosomes, probably reflecting each step of the sex chromosome evolutionary history. The nucleotide substitution per site was significantly higher in the Y- and Z-genes than in the X- and W- genes, indicating male-driven mutation. The ratio of nonsynonymous to synonymous nucleotide substitution rates was higher in the X- and W-genes than in the Y- and Z-genes, with a female bias. Allelic expression in gonad, brain, and muscle was significantly higher in the Y- and W-genes than in the X- and Z-genes, favoring heterogametic sex. The same set of sex-linked genes showed parallel evolution across the two distinct systems. In contrast, the unique genomic region of the sex chromosomes demonstrated a difference between the two systems, with even and extremely high expression ratios of W/Z and Y/X, respectively.

Characterization of the novel HLA-DQB1*02:02:23 allele in a Russian hematopoietic stem cell donor.

HLA-DQB1*02:02:23 differs from HLA-DQB1*02:02:01:01 by one synonymous nucleotide substitution at codon -6 in exon 1.

Allotetraploidization event of Coptis chinensis shared by all Ranunculales

Coptis chinensis Franch., also named Chinese goldthread is a member of Ranunculaceae in the order Ranunculales and represents an important lineage of early eudicots with traditional medicinal value. In our study, by using syntenic analysis combined with phylogenomic analysis of C. chinensis and four other representative genomes from basal and core eudicots, we confirmed that the WGD event in C. chinensis was shared by Aquilegia coerulea and Papaver somniferum L. and quickly occurred after Ranunculales diverged from other eudicots, likely a Ranunculales common tetraploidization (RCT). The synonymous nucleotide substitutions at synonymous sites distribution of syntenic blocks across these genomes showed that the evolutionary rate of the P. somniferum genome is faster than that of the C. chinensis genome by approximately 13.7%, possibly due to Papaveraceaes having an additional special tetraploidization event (PST). After Ks correction, the RCT dated to 115–130 million years ago (MYA), which was close to the divergence of Ranunculaceaes and Papaveraceaes approximately 115.45–130.51 MYA. Moreover, we identified homologous genes related to polyploidization and speciation and constructed multiple sequence alignments with different reference genomes. Notably, the event-related subgenomes in the basal genomes all showed genomic fractionation bias, suggesting a likely allopolyploid nature of the RCT, PST and T-Alpha and T-Beta events in Tetracentron sinense. In addition, we detected that the sixteen P450 subfamilies were markedly expanded in the genomes of Ranunculales, and most of them were related to the RCT and PST events. We constructed a new platform for Early Eudicot Comparative Genomic Research (http://www.cgrpoee.top/index.html) to store more information. In summary, our findings support the WGD of C. chinensis shared by Ranunculales, which is likely an allotetraploidization event. This present effort offered new insights into the evolution of key polyploidization events and the genes related to secondary metabolites during the diversification of early eudicots.

Cyanophages as an important factor in the early evolution of oxygenic photosynthesis

Cyanophages are viruses that infect cyanobacteria. An interesting feature of many of them is the presence of psbA and psbD, genes that encode D1 and D2 proteins, respectively. The D1 and D2 are core proteins of the photosystem II (PSII) in cyanobacteria, algae and plants and influence the proper function of oxygenic photosynthesis (OP) in all oxyphototrophs on Earth. The frequent occurrence of psbA and psbD in cyanophages raises the question whether these genes coevolved with hosts during the early stages of cyanophage and cyanobacterial evolution, or whether they are direct descendants of genes adopted from the genomes of cyanobacterial hosts. The phylogeny of D1/D2 proteins encoded in the genomes of selected cyanophages and oxyphototrophs was reconstructed. In addition, common ancestral sequences of D1 and D2 proteins were predicted for cyanophages and oxyphototrophs. Based on this, the reconstruction of the 3D structures of D1 and D2 proteins was performed. In addition, the ratio of non-synonymous to synonymous (dN/dS) nucleotide substitutions in the coding sequences (CDSs) of psbA and psbD was determined. The results of the predicted spatial structures of the D1 and D2 proteins and purifying selection for the CDSs of psbA and psbD suggest that they belong to the ancient proteins, which may have formed the primordial PSII. It cannot be ruled out that they involved in water oxidation in cyanobacteria-like organisms at early stages of the evolution of life on Earth and coevolved with ancient cyanophages. The data are also discussed in the context of the origin of viruses.