Molecular identification of vivax malaria relapse patients in the Yunnan Province based on homology analysis of the Plasmodium vivax circumsporozoite protein gene.

More than 85% of the malaria burden in the Yunnan Province is caused by imported vivax malaria, and Yunnan is also where the majority of vivax malaria patients are diagnosed in China. Timely removal of the infection sources of Plasmodium vivax and its breeding environment remains the key to eliminating the secondary transmission of imported malaria. To that end, blood samples were collected from cases diagnosed and revalidated as single species infection with P. vivax in the Yunnan Province from 2013 to 2020. Specifically, samples from vivax malaria patients with suspected relapses episodes were subjected to PCR amplification, product sequencing, and analysis of the P. vivax circumsporozoite protein (pvcsp) gene. In total, 77 suspected relapse patients were identified out of 2484 cases infected with P. vivax, with a total of 81 recurrent episodes. A total of 156 CDS (coding DNA sequence) chains were obtained through PCR amplification and sequencing of the pvcsp gene from 159 blood samples, 121 of which can be matched to the paired sequences of 59 vivax malaria patients with both primary attack and recurrent experience. Of the 59 pairs of pvcsp gene sequences, every one of 31 pairs showed only one haplotype and no variant sites (VS), meaning every two paired sequence was completely homologous. Every one of the remaining 28 paired sequences had two haplotypes but no length polymorphism, indicating that the paired sequences was "weakly heterologous" with no fragment insertions (or deletions). All 59 vivax malaria patients with recurrences were caused by the activation of P. vivax hypnozoites originated from the same population as the primary infection. The paired analysis of the similarity between high variant genes allowed the identification of relapse episodes caused by P. vivax homologous hypnozoites and also demonstrated pvcsp gene as one of the candidate molecular markers for tracing infection origin.

Prediction of the CYP2D6 enzymatic activity based on investigating of the CYP2D6 genotypes around the vivax malaria patients in Yunnan Province, China.

BACKGROUND: In recent years, the incidence rate of vivax malaria recurrence still had 3.1% in Yunnan Province population after eradication therapy using primaquine (PQ). In order to understand the specific failure reasons for preventing vivax malaria relapses, a preliminary exploration on the CYP2D6 enzyme activity was carried out in the vivax malaria patients in Yunnan Province population by analysing mutational polymorphism in the coding region of CYP2D6 gene. METHODS: Blood samples were collected from vivax malaria patients with suspected relapse (SR) and non-relapsed (NR) malaria in Yunnan Province. The DNA fragments containing 9 exons regions of human CYP2D6 gene were amplified by performing PCR and sequenced. The sequencing results were aligned by using DNAStar 11.0 to obtain the coding DNA sequence (CDS) of CYP2D6 gene. DnaSP 6.11.01 software was used to identify mutant polymorphisms and haplotypes of the CDS chain. The waterfall function of GenVisR package in R was utilized to visualize the mutational landscape. The alleles of CYP2D6 gene were identified according to the criteria prescribed by Human Cytochrome P450 (CYP) Allele Nomenclature Committee Database and the CYP2D6 enzyme activity was predicted based on diploid genotype. RESULTS: A total of 320 maternal CDS chains, including 63 from SR group and 257 from NR group, were obtained. Twelve mutant loci, including c.31 (rs769259), c.100 (rs1065852), c.271 (rs28371703), c.281 (rs28371704), c.294 (rs28371705), c.297 (rs200269944), c.336 (rs1081003), c.408 (rs1058164), c.505 (rs5030865), c.801 (rs28371718), c.886 (rs16947), and c.1,457 (rs1135840) were observed on the 640 CDS chains (including 320 maternal and 320 paternal chains). The high-frequency mutation at rs1135840 (0.703) and low-frequency mutation, such as rs28371703, were detected only in the SR group. The frequency of mutant rs1058164 and rs1135840 were significantly increased in the SR group ([Formula: see text]= 4.468, 5.889, P < 0.05), as opposed to the NR group. Of the 23 haplotypes (from Hap_1 to Hap_23), the nomenclatures of 11 allelic forms could be found: Hap_3 was non-mutant, Hap_2 accounted for the highest frequency (36.9%, 236/640), and Hap_9 had the most complex sequence structure, containing 7 loci mutations. Allele *10 was the most frequent among these genotypes (0.423). Among the allele *10 standard named genotypes, *1/*10, *1/*1 and *2/*10 were significantly more frequent in the NR group ([Formula: see text]= 3.911, P < 0.05) and all showed uncompromised enzyme activity; the impaired genotype *10/*39 was more frequent in the SR group ([Formula: see text]= 10.050, P < 0.05), and genotype *4/*4was detected only in the SR group. CONCLUSION: In the patients receiving PQ dosage in Yunnan Province population, both rs1135840 single nucleotide polymorphism and *10 allele form was common in the CYP2D6 gene. Low-frequency mutation sites, such as rs28371703, were only presented in patients with vivax malaria relapse.

The association of CYP2D6 gene polymorphisms in the full-length coding region with higher recurrence rate of vivax malaria in Yunnan Province, China.

BACKGROUND: Accumulating evidence suggest that compromised CYP2D6 enzyme activity caused by gene mutation could contribute to primaquine failure for the radical cure of vivax malaria. The current study aims to preliminarily reveal the association between the recurrence of vivax malaria in Yunnan Province and CYP2D6 gene mutation by analysing polymorphisms in the entire coding region of human CYP2D6 gene. METHODS: Blood samples were collected from patients with vivax malaria, who received "chloroquine and 8-day course of primaquine therapy" in Yunnan Province. The suspected relapsed cases were determined by epidemiological approaches and gene sequence alignment. PCR was conducted to amplify the CYP2D6 gene in the human genome, and the amplified products were then sequenced to compare with the non-mutation "reference" sequence, so as to ensure correct sequencing results and to determine 9 exon regions. Subsequently, the DNA sequences of 9 exons were spliced into the coding DNA sequence (CDS), which, by default, is known as maternal CDS. The paternal CDS was obtained by adjusting the bases according to the sequencing peaks. The mutation loci, haplotypes (star alleles), genotypes and odds ratios (OR) of all the CDSs were analysed. RESULTS: Of the119 maternal CDS chains in total with 1491 bp in length, 12 mutation sites in the 238 maternal and paternal CDS chains were detected. The c.408G > C mutation was most frequently detected in both suspected relapsed group (SR) and non-relapsed group (NR), reaching 85.2% (75/88) and 76.0% (114/150), respectively. The c.886C > T mutation was most closely related to the recurrence of vivax malaria (OR = 2.167, 95% CI 1.104-4.252, P < 0.05). Among the 23 haplotypes (Hap_1 ~ Hap_23), Hap_3 was non-mutant, and the sequence structure of Hap_9 was the most complicated one. Five star alleles, including *1, *2, *4, *10 and *39, were confirmed by comparison, and CYP2D6*10 allele accounted for the largest percentage (45.4%, 108/238). The frequency of CYP2D6*2 allele in the SR group was significantly higher than that in the NR group (Χ2 = 16.177, P < 0.05). Of the defined 24 genotypes, 8 genotypes, including *4/*4, *4/*o, *2/*39, *39/*m, *39/*x, *1/*r, *1/*n, and *v/*10, were detected only in the SR group. CONCLUSION: Mutation of CYP2D6*10 allele accounts for the highest proportion of vivax malaria cases in Yunnan Province. The mutations of c. 886C > T and CYP2D6*2 allele, which correspond to impaired PQ metabolizer phenotype, are most closely related to the relapse of vivax malaria. In addition, the genotype *4/*4 with null CYP2D6 enzyme function was only detected in the SR group. These results reveal the risk of defected CYP2D6 enzyme activity that diminishes the therapeutic effect of primaquine on vivax malaria.

Polymorphism analysis of propeller domain of k13 gene in Plasmodium ovale curtisi and Plasmodium ovale wallikeri isolates original infection from Myanmar and Africa in Yunnan Province, China.

BACKGROUND: Eighteen imported ovale malaria cases imported from Myanmar and various African countries have been reported in Yunnan Province, China from 2013 to 2018. All of them have been confirmed by morphological examination and 18S small subunit ribosomal RNA gene (18S rRNA) based PCR in YNRL. Nevertheless, the subtypes of Plasmodium ovale could not be identified based on 18S rRNA gene test, thus posing challenges on its accurate diagnosis. To help establish a more sensitive and specific method for the detection of P. ovale genes, this study performs sequence analysis on k13-propeller polymorphisms in P. ovale. METHODS: Dried blood spots (DBS) from ovale malaria cases were collected from January 2013 to December 2018, and the infection sources were confirmed according to epidemiological investigation. DNA was extracted, and the coding region (from 206th aa to 725th aa) in k13 gene propeller domain was amplified using nested PCR. Subsequently, the amplified products were sequenced and compared with reference sequence to obtain CDS. The haplotypes and mutation loci of the CDS were analysed, and the spatial structure of the amino acid peptide chain of k13 gene propeller domain was predicted by SWISS-MODEL. RESULTS: The coding region from 224th aa to 725th aa of k13 gene from P. ovale in 83.3% of collected samples (15/18) were amplified. Three haplotypes were observed in 15 samples, and the values of Ka/Ks, nucleic acid diversity index (π) and expected heterozygosity (He) were 3.784, 0.0095, and 0.4250. Curtisi haplotype, Wallikeri haplotype, and mutant type accounted for 73.3% (11/15), 20.0% (3/15), and 6.7% (1/15). The predominant haplotypes of P. ovale curtisi were determined in all five Myanmar isolates. Of the ten African isolates, six were identified as P. o. curtisi, three were P. o. wallikeri and one was mutant type. Base substitutions between the sequences of P. o. curtisi and P. o. wallikeri were determined at 38 loci, such as c.711. Moreover, the A > T base substitution at c.1428 was a nonsynonymous mutation, resulting in amino acid variation of T476S in the 476th position. Compared with sequence of P. o. wallikeri, the double nonsynonymous mutations of G > A and A > T at the sites of c.1186 and c.1428 leads to the variations of D396N and T476S for the 396th and 476th amino acids positions. For P. o. curtisi and P. o. wallikeri, the peptide chains in the coding region from 224th aa to 725th aa of k13 gene merely formed a monomeric spatial model, whereas the double-variant peptide chains of D396N and T476S formed homodimeric spatial model. CONCLUSION: The propeller domain of k13 gene in the P. ovale isolates imported into Yunnan Province from Myanmar and Africa showed high differentiation. The sequences of Myanmar-imported isolates belong to P. o. curtisi, while the sequences of African isolates showed the sympatric distribution from P. o. curtisi, P. o. wallikeri and mutant isolates. The CDS with a double base substitution formed a dimeric spatial model to encode the peptide chain, which is completely different from the monomeric spatial structure to encode the peptide chain from P. o. curtisi and P. o. wallikeri.

Identification and Comprehensive Analysis of circRNA-miRNA-mRNA Regulatory Networks in A2780 Cells Treated with Resveratrol.

Ovarian cancer (OC) is one of the most commonplace gynecological malignancies. This study explored the effects of resveratrol (RES) on OC cell proliferation and apoptosis. Proliferation activity was measured for A2780 cells treated with RES for 24 h and 48 h at concentrations of 0, 10, 25, 50, 75, 100, 150, 200, and 300 µM. RNA sequencing (RNA-seq) was performed to analyze the circular RNA (circRNA), microRNA (miRNA), and messenger RNA (mRNA) expression spectrum. The differentially expressed genes included 460 circRNAs, 1988 miRNAs, and 1671 mRNAs, and they were subjected to analyses including Gene Ontology, the Kyoto Encyclopedia of Genes and Genomes (KEGG), and Reactome enrichment. We selected signaling pathways enriched in the cell processes by mRNA KEGG, comprehensively analyzed the circRNA-miRNA-mRNA regulatory network, and verified several miRNAs expressed in the regulatory network diagram using the quantitative real-time polymerase chain reaction. The data showed that the cell proliferation of A2780 cells treated with RES for 24 h or 48 h decreased with increasing concentrations of RES. The circRNA-miRNA-mRNA regulatory network that we constructed provides new insights into the ability of RES to inhibit cell proliferation and promote apoptosis in A2780 cells.