PD-1 Targeted Antibody Discovery Using AI Protein Diffusion.

The programmed cell death protein 1 (PD-1, CD279) is an important therapeutic target in many oncological diseases. This checkpoint protein inhibits T lymphocytes from attacking other cells in the body and thus blocking it improves the clearance of tumor cells by the immune system. While there are already multiple FDA-approved anti-PD-1 antibodies, including nivolumab (Opdivo® from Bristol-Myers Squibb) and pembrolizumab (Keytruda® from Merck), there are ongoing efforts to discover new and improved checkpoint inhibitor therapeutics. In this study, we present multiple anti-PD-1 antibody fragments that were derived computationally using protein diffusion and evaluated through our scalable, in silico pipeline. Here we present nine synthetic Fv structures that are suitable for further empirical testing of their anti-PD-1 activity due to desirable predicted binding performance.

Human cytokine and coronavirus nucleocapsid protein interactivity using large-scale virtual screens.

Understanding the interactions between SARS-CoV-2 and the human immune system is paramount to the characterization of novel variants as the virus co-evolves with the human host. In this study, we employed state-of-the-art molecular docking tools to conduct large-scale virtual screens, predicting the binding affinities between 64 human cytokines against 17 nucleocapsid proteins from six betacoronaviruses. Our comprehensive in silico analyses reveal specific changes in cytokine-nucleocapsid protein interactions, shedding light on potential modulators of the host immune response during infection. These findings offer valuable insights into the molecular mechanisms underlying viral pathogenesis and may guide the future development of targeted interventions. This manuscript serves as insight into the comparison of deep learning based AlphaFold2-Multimer and the semi-physicochemical based HADDOCK for protein-protein docking. We show the two methods are complementary in their predictive capabilities. We also introduce a novel algorithm for rapidly assessing the binding interface of protein-protein docks using graph edit distance: graph-based interface residue assessment function (GIRAF). The high-performance computational framework presented here will not only aid in accelerating the discovery of effective interventions against emerging viral threats, but extend to other applications of high throughput protein-protein screens.

Vertical and temporal flight patterns of coffee berry borer (Coleoptera: Curculionidae) in Hawaii.

Coffee berry borer (Hypothenemus hampei Ferrari) (Coleoptera: Curculionidae) is the most damaging insect pest of coffee worldwide, causing significant losses in coffee yields and quality. Knowledge of vertical and temporal flight patterns in coffee berry borer could be used to optimize spray timing and precision targeting of areas within the coffee tree, which may be more susceptible. In the present study, we estimated the vertical distribution of coffee berry borer females using traps set at 1-m intervals up to 5 m in height. We also quantified coffee berry borer infestation in the low, mid, and high canopy and documented fruit availability. Temporal flight patterns were estimated using timer traps, and correlation analyses were conducted to determine the relationship between the timing of daily flight and weather variables. Across the 4 study sites, we observed that 77%-84% of the trap catch was at 1 m, 11%-20% was at 2 m, and 1%-4% was at 3-5 m in height. Fruit infestation was significantly higher in the low branches (35%) relative to the high branches (17%). Flight height remained the same year-round, regardless of fruit availability. Coffee berry borer flew in low numbers during the day and night but peaked from 12 to 4 PM. Daily flight was positively correlated with an increase in air temperature and wind speed and negatively correlated with relative humidity. Findings from this study suggest that pesticide sprays should target low- to mid-level branches at 1-2 m in height and aim to be conducted in the early afternoon when coffee berry borer are actively flying and most vulnerable to chemical controls.

The role of lung biopsy for diagnosis and prognosis of interstitial lung disease in systemic sclerosis: a systematic literature review.

BACKGROUND: The prognostic and theragnostic role of histopathological subsets in systemic sclerosis interstitial lung disease (SSc-ILD) have been largely neglected due to the paucity of treatment options and the risks associated with surgical lung biopsy. The novel drugs for the treatment of ILDs and the availability of transbronchial cryobiopsy provide a new clinical scenario making lung biopsy more feasible and a pivotal guide for treatment. The aim of our study was to investigate the usefulness of lung biopsy in SSc ILD with a systematic literature review (SLR). METHODS: PubMed, Embase and Cochrane databases were searched up to June 30, 2023. Search terms included both database-specific controlled vocabulary terms and free-text terms relating to lung biopsy and SSc-ILD diagnostic and prognosis. The SLR was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA). Studies were selected according to the PEO (population, exposure, and outcomes) framework and Quality assessment of diagnostic accuracy studies (QUADAS) were reported. RESULTS: We selected 14 articles (comprising 364 SSc-ILD patients). The paucity and heterogeneity of the studies prevented a systematic analysis. Diffuse cutaneous SSc was present in 30-100% of cases. Female predominance was observed in all studies (ranging from 64 to 100%). Mean age ranged from 42 to 64 years. Mean FVC was 73.98 (+/-17.3), mean DLCO was 59.49 (+/-16.1). Anti-Scl70 antibodies positivity was detected in 33% of cases (range: 0-69.6). All patients underwent surgical lung biopsies, and multiple lobes were biopsied in a minority of studies (4/14). Poor HRCT-pathologic correlation was reported with HRCT-NSIP showing histopathologic UIP in up to 1/3 of cases. Limited data suggest that SSc-UIP patients may have a worse prognosis and response to immunosuppressive treatment compared to other histopathologic patterns. CONCLUSIONS: The data from this SLR clearly show the paucity and heterogeneity of the studies reporting lung biopsy in SSc ILD. Moreover, they highlight the need for further research to address whether the lung biopsy can be helpful to refine prognostic prediction and guide therapeutic choices.

Bronchioalveolar organoids: A preclinical tool to screen toxicity associated with antibody-drug conjugates.

Despite extensive preclinical testing, cancer therapeutics can result in unanticipated toxicity to non-tumor tissue in patients. These toxicities may pass undetected in preclinical experiments due to modeling limitations involving poor biomimicry of 2-dimensional in vitro cell cultures and due to lack of interspecies translatability in in vivo studies. Instead, primary cells can be grown into miniature 3-dimensional structures that recapitulate morphological and functional aspects of native tissue, termed "organoids." Here, human bronchioalveolar organoids grown from primary alveolar epithelial cells were employed to model lung epithelium and investigate off-target toxicities associated with antibody-drug conjugates (ADCs). ADCs with three different linker-payload combinations (mafodotin, vedotin, and deruxtecan) were tested in bronchioalveolar organoids generated from human, rat, and nonhuman primate lung cells. Organoids demonstrated antibody uptake and changes in viability in response to ADC exposure that model in vivo drug sensitivity. RNA sequencing identified inflammatory activation in bronchioalveolar cells in response to deruxtecan. Future studies will explore specific cell populations involved in interstitial lung disease and incorporate immune cells to the culture.

Behavioral and transcriptomic analyses of mecp2 function in zebrafish.

Rett syndrome (RTT), a human neurodevelopmental disorder characterized by severe cognitive and motor impairments, is caused by dysfunction of the conserved transcriptional regulator Methyl-CpG-binding protein 2 (MECP2). Genetic analyses in mouse Mecp2 mutants, which exhibit key features of human RTT, have been essential for deciphering the mechanisms of MeCP2 function; nonetheless, our understanding of these complex mechanisms is incomplete. Zebrafish mecp2 mutants exhibit mild behavioral deficits but have not been analyzed in depth. Here, we combine transcriptomic and behavioral assays to assess baseline and stimulus-evoked motor responses and sensory filtering in zebrafish mecp2 mutants from 5 to 7 days post-fertilization (dpf). We show that zebrafish mecp2 function is required for normal thigmotaxis but is dispensable for gross movement, acoustic startle response, and sensory filtering (habituation and sensorimotor gating), and reveal a previously unknown role for mecp2 in behavioral responses to visual stimuli. RNA-seq analysis identified a large gene set that requires mecp2 function for correct transcription at 4 dpf, and pathway analysis revealed several pathways that require MeCP2 function in both zebrafish and mammals. These findings show that MeCP2's function as a transcriptional regulator is conserved across vertebrates and supports using zebrafish to complement mouse modeling in elucidating these conserved mechanisms.

Comparative transcriptomics reveal differential gene expression among Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms.

The documentation of Plasmodium vivax malaria across Africa especially in regions where Duffy negatives are dominant suggests possibly alternative erythrocyte invasion mechanisms. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa is unclear. In this study, we examined the expression of 4,404 gene transcripts belong to 12 functional groups and 43 erythrocyte binding gene candidates in Ethiopian isolates and compared them with the Cambodian and Brazilian P. vivax transcriptomes. Overall, there were 10-26% differences in the gene expression profile amongst geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian isolates. Members of the reticulocyte binding protein PvRBP2a and PvRBP3 expressed six-fold higher than Duffy binding protein PvDBP1 and 60-fold higher than PvEBP/DBP2 in the Ethiopian isolates. Other genes including PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression patterns were observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, gametocyte genes including PvAP2-G, PvGAP (female gametocytes), and Pvs47 (male gametocytes) were highly expressed across geographical samples.

Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects.

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7, TPM1, and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes such as MYH2, TPM2, and TNNI2 that encode parts of the skeletal muscle sarcomere cause muscle diseases affecting skeletal muscle, such as distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7) encoding sarcomeric proteins in which the same pathogenic variant affects skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain pathogenic variants that also cause cardiac abnormalities. We report five families with DA because of heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 (ACTC1). ACTC1 encodes a highly conserved actin that binds to myosin in cardiac and skeletal muscle. Pathogenic variants in ACTC1 have been found previously to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition because of variants in ACTC1 and suggests that some functions of ACTC1 are shared in cardiac and skeletal muscle.

Phylogeny, ancestral ranges and reclassification of sand dollars.

Classification of the Class Echinoidea is under significant revision in light of emerging molecular phylogenetic evidence. In particular, the sister-group relationships within the superorder Luminacea (Echinoidea: Irregularia) have been considerably updated. However, the placement of many families remains largely unresolved due to a series of incongruent evidence obtained from morphological, paleontological, and genetic data for the majority of extant representatives. In this study, we investigated the phylogenetic relationships of 25 taxa, belonging to eleven luminacean families. We proposed three new superfamilies: Astriclypeoidea, Mellitoidea, and Taiwanasteroidea (including Dendrasteridae, Taiwanasteridae, Scutellidae, and Echinarachniidae), instead of the currently recognized superfamily Scutelloidea Gray, 1825. In light of the new data obtained from ten additional species, the historical biogeography reconstructed shows that the tropical western Pacific and eastern Indian Oceans are the cradle for early sand dollar diversification. Hothouse conditions during the late Cretaceous and early Paleogene were coupled with diversification events of major clades of sand dollars. We also demonstrate that Taiwan fauna can play a key role in terms of understanding the major Cenozoic migration and dispersal events in the evolutionary history of Luminacea.

Variants in ACTC1 underlie distal arthrogryposis accompanied by congenital heart defects.

Contraction of the human sarcomere is the result of interactions between myosin cross-bridges and actin filaments. Pathogenic variants in genes such as MYH7 , TPM1 , and TNNI3 that encode parts of the cardiac sarcomere cause muscle diseases that affect the heart, such as dilated cardiomyopathy and hypertrophic cardiomyopathy. In contrast, pathogenic variants in homologous genes MYH2 , TPM2 , and TNNI2 , that encode parts of the skeletal muscle sarcomere, cause muscle diseases affecting skeletal muscle, such as the distal arthrogryposis (DA) syndromes and skeletal myopathies. To date, there have been few reports of genes (e.g., MYH7 ) encoding sarcomeric proteins in which the same pathogenic variant affects both skeletal and cardiac muscle. Moreover, none of the known genes underlying DA have been found to contain mutations that also cause cardiac abnormalities. We report five families with DA due to heterozygous missense variants in the gene actin, alpha, cardiac muscle 1 ( ACTC1 ). ACTC1 encodes a highly conserved actin that binds to myosin in both cardiac and skeletal muscle. Mutations in ACTC1 have previously been found to underlie atrial septal defect, dilated cardiomyopathy, hypertrophic cardiomyopathy, and left ventricular noncompaction. Our discovery delineates a new DA condition due to mutations in ACTC1 and suggests that some functions of actin, alpha, cardiac muscle 1 are shared in cardiac and skeletal muscle.

Comparative transcriptomics reveal differential gene expression in Plasmodium vivax geographical isolates and implications on erythrocyte invasion mechanisms.

Plasmodium vivax uses Duffy binding protein (PvDBP1) to bind to the Duffy Antigen-Chemokine Receptor (DARC) to invade human erythrocytes. Individuals who lack DARC expression (Duffy-negative) are thought to be resistance to P. vivax. In recent years, P. vivax malaria is becoming more prevalent in Africa with a portion of these cases detected in Duffy-negatives. Apart from DBP1, members of the reticulocyte binding protein (RBP) and tryptophan-rich antigen (TRAg) families may also play a role in erythrocyte invasion. While the transcriptomes of the Southeast Asian and South American P. vivax are well documented, the gene expression profile of P. vivax in Africa and more specifically the expression level of several erythrocyte binding gene candidates as compared to DBP1 are largely unknown. This paper characterized the first P. vivax transcriptome in Africa and compared with those from the Southeast Asian and South American isolates. The expression of 4,404 gene transcripts belong to 12 functional groups including 43 specific erythrocyte binding gene candidates were examined. Overall, there were 10-26% differences in the gene expression profile amongst the geographical isolates, with the Ethiopian and Cambodian P. vivax being most similar. Majority of the gene transcripts involved in protein transportation, housekeeping, and host interaction were highly transcribed in the Ethiopian P. vivax. Erythrocyte binding genes including PvRBP2a and PvRBP3 expressed six-fold higher than PvDBP1and 60-fold higher than PvEBP/DBP2. Other genes including PvRBP1a, PvMSP3.8, PvMSP3.9, PvTRAG2, PvTRAG14, and PvTRAG22 also showed relatively high expression. Differential expression was observed among geographical isolates, e.g., PvDBP1 and PvEBP/DBP2 were highly expressed in the Cambodian but not the Brazilian and Ethiopian isolates, whereas PvRBP2a and PvRBP2b showed higher expression in the Ethiopian and Cambodian than the Brazilian isolates. Compared to Pvs25, the standard biomarker for detecting female gametocytes, PvAP2-G (PVP01_1440800), GAP (PVP01_1403000), and Pvs47 (PVP01_1208000) were highly expressed across geographical samples. These findings provide an important baseline for future comparisons of P. vivax transcriptomes from Duffy-negative infections and highlight potential biomarkers for improved gametocyte detection.

Genetic variations of Plasmodium falciparum circumsporozoite protein and the impact on interactions with human immunoproteins and malaria vaccine efficacy.

In October 2021, the world's first malaria vaccine RTS,S was endorsed by WHO for broad use in children, despite its low efficacy. This study examined polyclonal infections and the associations of parasite genetic variations with binding affinity to human leukocyte antigen (HLA). Multiplicity of infection was determined by amplicon deep sequencing of PfMSP1. Genetic variations in PfCSP were examined across 88 samples from Ghana and analyzed together with 1655 PfCSP sequences from other African and non-African isolates. Binding interactions of PfCSP peptide variants and HLA were predicted using NetChop and HADDOCK. High polyclonality was detected among infections, with each infection harboring multiple non-3D7 PfCSP variants. Twenty-seven PfCSP haplotypes were detected in the Ghanaian samples, and they broadly represented PfCSP diversity across Africa. The number of genetic differences between 3D7 and non-3D7 PfCSP variants does not influence binding to HLA. However, CSP peptide length after proteolytic degradation significantly affects its molecular weight and binding affinity to HLA. Despite the high diversity of HLA, the majority of the HLAI and II alleles interacted/bound with all Ghana CSP peptides. Multiple non-3D7 strains among P. falciparum infections could impact the effectiveness of RTS,S. Longer peptides of the Th2R/Th3R CSP regions should be considered in future versions of RTS,S.

Nationwide Surveillance of Pfhrp2 Exon 2 Diversity in Plasmodium falciparum Circulating in Symptomatic Malaria Patients Living in Ghana.

Reports of increasing false-negative HRP2-based rapid diagnostic test results across Africa require constant monitoring of factors associated with these false-negative outcomes, as failure of this diagnostic tool will have severe consequences on malaria treatment and control programs. This study characterized the extent of genetic diversity in the Plasmodium falciparum histidine-rich protein 2 (Pfhrp2) gene in P. falciparum isolates from symptomatic malaria patients across the regions of Ghana. Exon 2 of Pfhrp2 was amplified from gDNA using polymerase chain reaction. All Pfhrp2-negative samples were subjected to Pf18S rRNA and Pfmsp2 gene amplifications. The amplified Pfhrp2 exon 2 fragments from clonal samples were sent for commercial Sanger sequencing. The type and number of PfHRP2 repeats, classified based on repeat types previously reported, were estimated from the sequence data and compared among geographical regions. About 81% (2,333/2,890) of the original microscopy positive DBS were available and used in this study. The Pfhrp2 exon 2 amplification was successful in 98.5% (2,297/2,333) of the tested samples, with band size ranging from 400 bp to 1,050 bp. A total of 13 out of the 24 previously reported repeat types were identified among the samples, with three samples lacking both type 2 and type 7 repeat motifs. This study suggested that the genetic diversity of Pfhrp2 exon 2 identified in P. falciparum circulating in symptomatic malaria patients in Ghana is unlikely to influence the sensitivity and specificity of HRP2 RDT-based diagnosis.

Okra (Abelmoschus esculentus) in a refugee context in East Africa: Kitchen gardening helps with mineral provision.

Kitchen gardening is considered a way to reconnect with agriculture and complement the cereal-based relief food offered to refugees in East Africa. This work aimed at profiling mineral content of okra in four refugee camps and settlements located in Ethiopia and Uganda and its contribution to adequate intake (AIs) or recommended dietary allowances (RDAs) for young children and pregnant and lactating women (PLW). The study also evaluated the applicability of portable X-ray fluorescence (PXRF) as compared with inductively coupled plasma mass spectrometry (ICP-MS) for mineral profiling of okra powder samples. The contents of minerals (mg kg-1) from the ICP-MS readings were in the following ranges: K (14,385-33,294), Ca (2610-14,090), P (3178-13,248), Mg (3896-7986), Cu (3.81-19.3), Fe (75.7-1243), Zn (33-141) and Mn (23.1-261). Regardless of geographic origin, at low-end consumption probability (17 g day-1 for young children and 68 g day-1 for PLW), okra could contribute ˂ 15% (2.7-12.9%) AI for macro-minerals (K and Ca). In addition, the contributions to RDA values for Fe and Zn, elements of known public health interest, ranged from 4.5 to 34.7% for young children. Interestingly, regression lines revealed strong agreement between ICP-MS and PXRF readings for Mn and Zn, with R2 values > 0.91. This information is useful in support of nutrition-sensitive kitchen gardening programs through scaling culturally important crops in refugee settings. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s42452-021-04898-6.

Further delineation of van den Ende-Gupta syndrome: Genetic heterogeneity and overlap with congenital heart defects and skeletal malformations syndrome.

Van den Ende-Gupta syndrome (VDEGS) is a rare autosomal recessive condition characterized by distinctive facial and skeletal features, and in most affected persons, by biallelic pathogenic variants in SCARF2. We review the type and frequency of the clinical features in 36 reported individuals with features of VDEGS, 15 (42%) of whom had known pathogenic variants in SCARF2, 6 (16%) with negative SCARF2 testing, and 15 (42%) not tested. We also report three new individuals with pathogenic variants in SCARF2 and clinical features of VDEGS. Of the six persons without known pathogenic variants in SCARF2, three remain unsolved despite extensive genetic testing. Three were found to have pathogenic ABL1 variants using whole exome sequencing (WES) or whole genome sequencing (WGS). Their phenotype was consistent with the congenital heart disease and skeletal malformations syndrome (CHDSKM), which has been associated with ABL1 variants. Of the three unsolved cases, two were brothers who underwent WGS and targeted long-range sequencing of both SCARF2 and ABL1, and the third person who underwent WES and RNA sequencing for SCARF2. Because these affected individuals with classical features of VDEGS lacked a detectable pathogenic SCARF2 variant, genetic heterogeneity is likely. Our study shows the importance of performing genetic testing on individuals with the VDEGS "phenotype," either as a targeted gene analysis (SCARF2, ABL1) or WES/WGS. Additionally, individuals with the combination of arachnodactyly and blepharophimosis should undergo echocardiography while awaiting results of molecular testing due to the overlapping physical features of VDEGS and CHDSKM.

Gene Polymorphisms Among Plasmodium vivax Geographical Isolates and the Potential as New Biomarkers for Gametocyte Detection.

The unique biological features of Plasmodium vivax not only make it difficult to control but also to eliminate. For the transmission of the malaria parasite from infected human to the vector, gametocytes play a major role. The transmission potential of a malarial infection is inferred based on microscopic detection of gametocytes and molecular screening of genes in the female gametocytes. Microscopy-based detection methods could grossly underestimate the reservoirs of infection as gametocytes may occur as submicroscopic or as micro- or macro-gametocytes. The identification of genes that are highly expressed and polymorphic in male and female gametocytes is critical for monitoring changes not only in their relative proportions but also the composition of gametocyte clones contributing to transmission over time. Recent transcriptomic study revealed two distinct clusters of highly correlated genes expressed in the P. vivax gametocytes, indicating that the male and female terminal gametocytogeneses are independently regulated. However, the detective power of these genes is unclear. In this study, we compared genetic variations of 15 and 11 genes expressed, respectively, in the female and male gametocytes among P. vivax isolates from Southeast Asia, Africa, and South America. Further, we constructed phylogenetic trees to determine the resolution power and clustering patterns of gametocyte clones. As expected, Pvs25 (PVP01_0616100) and Pvs16 (PVP01_0305600) expressed in the female gametocytes were highly conserved in all geographical isolates. In contrast, genes including 6-cysteine protein Pvs230 (PVP01_0415800) and upregulated in late gametocytes ULG8 (PVP01_1452800) expressed in the female gametocytes, as well as two CPW-WPC family proteins (PVP01_1215900 and PVP01_1320100) expressed in the male gametocytes indicated considerably high nucleotide and haplotype diversity among isolates. Parasite samples expressed in male and female gametocyte genes were observed in separate phylogenetic clusters and likely represented distinct gametocyte clones. Compared to Pvs25, Pvs230 (PVP01_0415800) and a CPW-WPC family protein (PVP01_0904300) showed higher expression in a subset of Ethiopian P. vivax samples. Thus, Pvs230, ULG8, and CPW-WPC family proteins including PVP01_0904300, PVP01_1215900, and PVP01_1320100 could potentially be used as novel biomarkers for detecting both sexes of P. vivax gametocytes in low-density infections and estimating transmission reservoirs.