Chromosome-level genome assembly defines female-biased genes associated with sex determination and differentiation in the human blood fluke Schistosoma japonicum.

Schistosomiasis is a neglected tropical disease of humans caused by blood flukes of the genus Schistosoma, the only dioecious parasitic flatworm. Although aspects of sex determination, differentiation and reproduction have been studied in some Schistosoma species, almost nothing is known for Schistosoma japonicum, the causative agent of schistosomiasis japonica. This mainly reflects the lack of high-quality genomic and transcriptomic resources for this species. As current genomes for S. japonicum are highly fragmented, we assembled and report a chromosome-level reference genome (seven autosomes, the Z-chromosome and partial W-chromosome), achieving a substantially enhanced gene annotation. Utilizing this genome, we discovered that the sex chromosomes of S. japonicum and its congener S. mansoni independently suppressed recombination during evolution, forming five and two evolutionary strata, respectively. By exploring the W-chromosome and sex-specific transcriptomes, we identified 35 W-linked genes and 257 female-preferentially transcribed genes (FTGs) from our chromosomal assembly and uncovered a signature for sex determination and differentiation in S. japonicum. These FTGs clustering within autosomes or the Z-chromosome exhibit a highly dynamic transcription profile during the pairing of female and male schistosomula, thereby representing a critical phase for the maturation of the female worms and suggesting distinct layers of regulatory control of gene transcription at this development stage. Collectively, these data provide a valuable resource for further functional genomic characterization of S. japonicum, shed light on the evolution of sex chromosomes in this highly virulent human blood fluke, and provide a pathway to identify novel targets for development of intervention tools against schistosomiasis.

Full-Length Transcriptome Analysis of Plasmodium falciparum by Single-Molecule Long-Read Sequencing.

Malaria, an infectious disease caused by Plasmodium parasites, still accounts for amounts of deaths annually in last decades. Despite the significance of Plasmodium falciparum as a model organism of malaria parasites, our understanding of gene expression of this parasite remains largely elusive since lots of progress on its genome and transcriptome are based on assembly with short sequencing reads. Herein, we report the new version of transcriptome dataset containing all full-length transcripts over the whole asexual blood stages by adopting a full-length sequencing approach with optimized experimental conditions of cDNA library preparation. We have identified a total of 393 alternative splicing (AS) events, 3,623 long non-coding RNAs (lncRNAs), 1,555 alternative polyadenylation (APA) events, 57 transcription factors (TF), 1,721 fusion transcripts in P. falciparum. Furthermore, the shotgun proteome was performed to validate the full-length transcriptome of P. falciparum. More importantly, integration of full-length transcriptomic and proteomic data identified 160 novel small proteins in lncRNA regions. Collectively, this full-length transcriptome dataset with high quality and accuracy and the shotgun proteome analyses shed light on the complex gene expression in malaria parasites and provide a valuable resource for related functional and mechanistic researches on P. falciparum genes.

Rrp6 Regulates Heterochromatic Gene Silencing via ncRNA RUF6 Decay in Malaria Parasites.

The heterochromatin environment plays a central role in silencing genes associated with the malaria parasite's development, survival in the host, and transmission to the mosquito vector. However, the underlying mechanism regulating the dynamic chromatin structure is not understood yet. Here, we have uncovered that Plasmodium falciparum Rrp6, an orthologue of eukaryotic RNA exosome-associated RNase, controls the silencing of heterochromatic genes. PfRrp6 knockdown disrupted the singular expression of the GC-rich ncRNA RUF6 family, a known critical regulator of virulence gene expression, through the stabilization of the nascent transcripts. Mechanistic investigation showed that the accumulation of the multiple RUF6 ncRNAs triggered local chromatin remodeling in situ, which activated their adjacent var genes. Strikingly, chromatin isolation by RNA purification analysis (ChIRP-seq) revealed that a remarkable RUF6 ncRNA had interacted with distal heterochromatin regions directly and stimulated a global derepression effect on heterochromatic genes, including all variant gene families and the sexual commitment-associated regulator ap2-g gene. Collectively, Rrp6 appears to conduct the epigenetic surveillance of heterochromatic gene expression through controlling RUF6 levels, thereby securing antigenic variation and sexual commitment of malaria parasites during the infection of the host.IMPORTANCE Malaria remains a major public health and economic burden. The heterochromatin environment controls the silencing of genes associated with the fate of malaria parasites. Previous studies have demonstrated that a group of GC-rich ncRNAs (RUF6) is associated with the mutually exclusive expression of var genes, but the underlying mechanisms remain elusive. Here, through a series of genetic manipulation and genome-wide multiomics analysis, we have identified the plasmodial orthologue of RNA exosome-associated Rrp6 as an upstream regulator of RUF6 expression and revealed that the dysregulation of RUF6 upon Rrp6 knockdown triggered local chromatin alteration, thereby activating most heterochromatic genes via direct interaction of RUF6 and distal gene loci. This finding not only uncovered the in-depth mechanism of RUF6-mediated regulation of heterochromatic genes but also identified Rrp6 as a novel regulator of gene expression in human malaria parasites, which provides a new target for developing intervention strategies against malaria.

The cryptic unstable transcripts are associated with developmentally regulated gene expression in blood-stage Plasmodium falciparum.

The tight gene expression regulation controls the development and pathogenesis of human malaria parasite Plasmodium falciparum throughout the complex life cycle. Recent studies have revealed the pervasive nascent transcripts in the genome of P. falciparum, suggesting the existence of a hidden transcriptome involved in the dynamic gene expression. However, the landscape and related biological functions of nascent non-coding RNAs (ns-ncRNAs) are still poorly explored. Here we profiled the transcription dynamics of nascent RNAs by rRNA-depleted and stranded RNA sequencing over the course of 48-h intraerythrocytic developmental cycle (IDC). We identified the genome-wide sources of a total of 2252 ns-ncRNAs, mostly originating from intergenic and untranslated regions of annotated genes. By integrating the nascent RNA abundances with ATAC-seq and ChIP-seq analysis, we uncovered the euchromatic microenvironment surrounding the ns-ncRNA loci, and revealed a positive correlation between ns-ncRNAs and corresponding mRNA abundances. Finally, by gene knock-down strategy, we showed that the cooperation of RNA exosome catalytic subunit PfDis3 and PfMtr4 cofactor played a major role in ns-ncRNAs degradation. Collectively, this study contributes to understanding of the potential roles of short-lived nascent ncRNAs in regulating gene expression in malaria parasites.

Disruption of the RNA exosome reveals the hidden face of the malaria parasite transcriptome.

Antisense transcription emerges as a key regulator of important biological processes in the human malaria parasite Plasmodium falciparum. RNA-processing factors, however, remain poorly characterized in this pathogen. Here, we purified the multiprotein RNA exosome complex of malaria parasites by affinity chromatography, using HA-tagged PfRrp4 and PfDis3 as the ligands. Seven distinct core exosome subunits (PfRrp41, PfMtr3, PfRrp42, PfRrp45, PfRrp4, PfRrp40, PfCsl4) and two exoribonuclease proteins PfRrp6 and PfDis3 are identified by mass spectrometry. Western blot analysis detects Dis3 and Rrp4 predominantly in the cytoplasmic fraction during asexual blood stage development. An inducible gene knock out of the PfDis3 subunit reveals the upregulation of structural and coding RNA, but the vast majority belongs to antisense RNA. Furthermore, we detect numerous types of cryptic unstable transcripts (CUTs) linked to virulence gene families including antisense RNA in the rif gene family. Our work highlights the limitations of steady-state RNA analysis to predict transcriptional activity and link the RNA surveillance machinery directly with post-transcriptional control and gene expression in malaria parasites.

Dual regulatory effects of non-coding GC-rich elements on the expression of virulence genes in malaria parasites.

As the primary virulence factor of falciparum malaria, var genes harboring mutually exclusive expression pattern lead to antigenic variation and immune evasion of this pathogen in human host. Although various mechanisms contribute to silence of var genes, little is known of transcriptional activation pathways of a single var gene and maintenance of its active state with other silent var loci. Here, we report a monoallelic expression pattern of the non-coding GC-elements flanking chromosomal internal var genes, and transcript from the active one was required for activation of the var gene in the same array. Meanwhile, GFP reporter assays revealed a repressive effect on the adjacent gene induced by DNA motifs of the insulator-like GC-element, which was linked to heterochromatin subnuclear localization. Taken together, these data for the first time provide experimental evidence of the dual cis- and trans-acting regulatory functions of the GC-elements in both silence and activation of var genes, which would advance our understanding of the complex regulatory network of the virulence gene family in P. falciparum.

Evaluation of the population structure and genetic diversity of Plasmodium falciparum in southern China.

BACKGROUND: Yunnan and Hainan provinces are the two major endemic regions for Plasmodium falciparum malaria in China. However, few studies have investigated the characteristics of this parasite. Therefore, this study aimed to evaluate the genetic diversity and population structure of P. falciparum to predict the geographic origin of falciparum malaria. METHODS: Thirteen highly polymorphic microsatellite loci were studied to estimate the genetic diversity and population structure of 425 P. falciparum isolates obtained from blood samples collected from Yunnan and Hainan provinces of South China. The isolates were analysed for genetic diversity, linkage disequilibrium, and population structure. The parasite populations were clustered into two subgroups (i.e., Yunnan and Hainan) and a classification algorithm was used to identify molecular markers for classifying the P. falciparum populations. RESULTS: All 13 microsatellite loci were highly polymorphic, with the number of alleles per locus varying from 5 to 20. The mean expected heterozygosity (He) in Yunnan and Hainan was 0.766 ± 0.036 and 0.677 ± 0.039, respectively, revealing a moderate high level of genetic diversity. Significant linkage disequilibrium was found for some regions of Yunnan (Lazan county and Xishuangbanna region) and Hainan (Dongfang city and Sanya city) province. According to the classification algorithm, a combination of three microsatellites could be used as a discriminatory marker to identify the origin of P. falciparum isolates. CONCLUSIONS: The results on the genetic structure of P. falciparum populations from South China provide a basis for developing a genetic marker-based tool to trace the source of the parasite infections and consequently improve malaria control and elimination strategies.

Distinctive origin and spread route of pyrimethamine-resistant Plasmodium falciparum in southern China.

Southeast Asia (the Thailand-Cambodia border) has been considered the primal epicenter for most antimalarial drug resistance; however, numerous molecular epidemiological studies have successively reported multiple independent origins of sulfadoxine-pyrimethamine (SP) resistance-associated Plasmodium falciparum dhfr (pfdhfr) and pfdhps alleles in other areas. To better understand the origin and evolutionary pathway of the SP resistance in Southeast Asia, a total of 374 P. falciparum field isolates from the Yunnan-Burma border and Hainan Island in southern China have been collected for comprehensive investigations on the mutation patterns of the pfdhfr/pfdhps genes as well as their microsatellite haplotypes. By comparative analysis of single-nucleotide polymorphism (SNP) genotyping and flanking microsatellite haplotypes, we reveal a unique origin of pyrimethamine-resistant mutations in Pfdhfr gene in Hainan Island and an oriented spread route of the pyrimethamine resistance from the Thailand-Cambodia border into the Hainan area, which reflects the geographical traits and SP administration histories in the two geographically independent areas. Moreover, genetic linkages between the high-level SP resistance-conferring pfdhfr/pfdhps alleles have been established in the isolates from the Yunnan-Burma border, raising the concern of a genetic basis in adopting combination chemotherapies against falciparum malaria.

PfSETvs methylation of histone H3K36 represses virulence genes in Plasmodium falciparum.

The variant antigen Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1), which is expressed on the surface of P. falciparum-infected red blood cells, is a critical virulence factor for malaria. Each parasite has 60 antigenically distinct var genes that each code for a different PfEMP1 protein. During infection the clonal parasite population expresses only one gene at a time before switching to the expression of a new variant antigen as an immune-evasion mechanism to avoid the host antibody response. The mechanism by which 59 of the 60 var genes are silenced remains largely unknown. Here we show that knocking out the P. falciparum variant-silencing SET gene (here termed PfSETvs), which encodes an orthologue of Drosophila melanogaster ASH1 and controls histone H3 lysine 36 trimethylation (H3K36me3) on var genes, results in the transcription of virtually all var genes in the single parasite nuclei and their expression as proteins on the surface of individual infected red blood cells. PfSETvs-dependent H3K36me3 is present along the entire gene body, including the transcription start site, to silence var genes. With low occupancy of PfSETvs at both the transcription start site of var genes and the intronic promoter, expression of var genes coincides with transcription of their corresponding antisense long noncoding RNA. These results uncover a previously unknown role of PfSETvs-dependent H3K36me3 in silencing var genes in P. falciparum that might provide a general mechanism by which orthologues of PfSETvs repress gene expression in other eukaryotes. PfSETvs knockout parasites expressing all PfEMP1 proteins may also be applied to the development of a malaria vaccine.