Discoba (Excavata) is an ancient group of eukaryotes with great morphological and ecological diversity. Unlike the other major divisions of Discoba (Jakobida and Euglenozoa), little is known about the mitochondrial DNAs (mtDNAs) of Heterolobosea. We have assembled a complete mtDNA genome from the aggregating heterolobosean amoeba, Acrasis kona, which consists of a single circular highly AT-rich (83.3%) molecule of 51.5 kb. Unexpectedly, A. kona mtDNA is missing roughly 40% of the protein-coding genes and nearly half of the transfer RNAs found in the only other sequenced heterolobosean mtDNAs, those of Naegleria spp. Instead, over a quarter of A. kona mtDNA consists of novel open reading frames. Eleven of the 16 protein-coding genes missing from A. kona mtDNA were identified in its nuclear DNA and polyA RNA, and phylogenetic analyses indicate that at least 10 of these 11 putative nuclear-encoded mitochondrial (NcMt) proteins arose by direct transfer from the mitochondrion. Acrasis kona mtDNA also employs C-to-U type RNA editing, and 12 homologs of DYW-type pentatricopeptide repeat (PPR) proteins implicated in plant organellar RNA editing are found in A. kona nuclear DNA. A mapping of mitochondrial gene content onto a consensus phylogeny reveals a sporadic pattern of relative stasis and rampant gene loss in Discoba. Rampant loss occurred independently in the unique common lineage leading to Heterolobosea + Tsukubamonadida and later in the unique lineage leading to Acrasis. Meanwhile, mtDNA gene content appears to be remarkably stable in the Acrasis sister lineage leading to Naegleria and in their distant relatives Jakobida.

Programmed -1 ribosomal frameshifting (-1 PRF) has been identified as a mechanism to regulate the expression of many viral genes and some cellular genes. The slippery site of -1 PRF has been well characterized, whereas the +1 PRF signal and the mechanism involved in +1 PRF remain poorly understood. Previous study confirmed that +1 PRF is required for the synthesis of protein products in several genes of ciliates from the genus Euplotes. To accurately assess the frequency of genes requiring frameshift in Euplotes, the macronuclear genome and transcriptome of Euplotes octocarinatus were analyzed in this study. A total of 3,700 +1 PRF candidate genes were identified from 32,353 transcripts, and the gene products of these putative +1 PRFs were mainly identified as protein kinases. Furthermore, we reported a putative suppressor tRNA of UAA which may provide new insights into the mechanism of +1 PRF in euplotids. For the first time, our transcriptome-wide survey of +1 PRF in E. octocarinatus provided a dataset which serves as a valuable resource for the future understanding of the mechanism underlying +1 PRF.

The class Myxosporea encompasses about 2,400 species, most of which are parasites of fish and cause serious damage in aquaculture. Due to the concerns about food safety issues and limited knowledge of Myxozoa life cycle and fish immune system, no chemicals, antibiotics or immune modulators are available to control myxozoa infection. Therefore, little can be done once Myxozoa establishment has occurred.

Scuticociliatosis is an invasive external or systemic infection caused by ciliated protozoa, mainly those within the subclass Scuticociliatia (scuticociliates). Many scuticociliates are fish pathogens, including Miamiensis avidus, Philasterides dicentrarchi, Pseudocohnilembus persalinus, and Uronema marinum. Our previous study showed that hemolysis-related genes derived from bacteria through horizontal gene transfer (HGT) may contribute to virulence in P. persalinus. Hemorrhagic lesions are a common feature of scuticociliatosis, but it is not known whether other scuticociliates also have bacteria-derived hemolysis-related genes. In this study, we constructed a high-quality macronuclear genome of another typical pathogenic scuticociliate, U. marinum. A total of 105 HGT genes were identified in this species, of which 35 were homologs of hemolysis-related genes (including hemolysin-like genes) that had previously been identified in P. persalinus. Sequencing of an additional five species from four scuticociliate families showed that bacteria-derived hemolysis-related genes (especially hemolysin-like genes) are widely distributed in scuticociliates. Based on these findings, we suggest that hemolysin-like genes may have originated before the divergence of scuticociliates.

The giant ciliate Stentor coeruleus (S. coeruleus) is a suitable model organism for studying morphogenesis and regeneration at the single-cell level. It contains a prominent structure on the anterior end of the cell, known as the oral apparatus (OA). OA can be induced to shed by urea treatment and then new OA regenerates via a series of defined morphological events and the cell resumes normal feeding activity. We identified OA constituents in S. coeruleus by mass spectrometry. A total of 882 OA-associated proteins were identified; the homologs of 181 of these are known OA constituents in other organisms. The expression pattern of OA-associated genes during regeneration was investigated using single-cell transcriptome sequencing. The expression of most OA-associated genes was high during regeneration, indicating their stable expression after OA shedding. We also identified OA-associated differentially expressed genes that may be involved in regulating OA reconstruction. In summary, this study gives preliminary insight into the molecular basis of OA in S. coeruleus.

To investigate the mechanisms underlying initiation of the sexual cell cycle in eukaryotes, we have focused on cyclins and cyclin-dependent kinases (CDKs) in the well-studied model ciliate, Tetrahymena thermophila We identified two genes, CDK19 and CYC9, which are highly co-expressed with the mating-associated factors MTA, MTB and HAP2. Both CDK19 and CYC9 were found to be essential for mating in T. thermophila Subcellular localization experiments suggested that these proteins are located at the oral area, including the conjugation junction area, and that CDK19 or CYC9 knockout prevents mating. We found that CDK19 and CYC9 form a complex, and also identified several additional subunits, which may have regulatory or constitutive functions. RNA sequencing analyses and cytological experiments showed that mating is abnormal in both ΔCDK19 and ΔCYC9, mainly at the entry to the co-stimulation stage. These results indicate that the CDK19-CYC9 complex initiates the sexual cell cycle in T. thermophila.

The peptidyl-prolyl cis-trans isomerase Pin1 is a pivotal therapeutic target in cancers, but the regulation of Pin1 protein stability is largely unknown. High Pin1 expression is associated with SUMO1-modified protein hypersumoylation in glioma stem cells (GSCs), but the underlying mechanisms remain elusive. Here we demonstrate that Pin1 is deubiquitinated and stabilized by USP34, which promotes isomerization of the sole SUMO E2 enzyme Ubc9, leading to SUMO1-modified hypersumoylation to support GSC maintenance. Pin1 interacts with USP34, a deubiquitinase with preferential expression and oncogenic function in GSCs. Such interaction is facilitated by Plk1-mediated phosphorylation of Pin1. Disruption of USP34 or inhibition of Plk1 promotes poly-ubiquitination and degradation of Pin1. Furthermore, Pin1 isomerizes Ubc9 to upregulate Ubc9 thioester formation with SUMO1, which requires CDK1-mediated phosphorylation of Ubc9. Combined inhibition of Pin1 and CDK1 with sulfopin and RO3306 most effectively suppresses orthotopic tumor growth. Our findings provide multiple molecular targets to induce Pin1 degradation and suppress hypersumoylation for cancer treatment.

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei-the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])-within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp-exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

It has become increasingly clear that the current taxonomy of clinical phenotypes is mixed with molecular heterogeneity, which potentially affects the treatment effect for involved patients. Defining the hidden molecular-distinct diseases using modern large-scale genomic approaches is therefore useful for refining clinical practice and improving intervention strategies. Given that microRNA expression profiling has provided a powerful way to dissect hidden genetic heterogeneity for complex diseases, the aim of the study was to develop a bioinformatics approach that identifies microRNA features leading to the hidden subtyping of complex clinical phenotypes. The basic strategy of the proposed method was to identify optimal miRNA clusters by iteratively partitioning the sample and feature space using the two-ways super-paramagnetic clustering technique. We evaluated the obtained optimal miRNA cluster by determining the consistency of co-expression and the chromosome location among the within-cluster microRNAs, and concluded that the optimal miRNA cluster could lead to a natural partition of disease samples. We applied the proposed method to a publicly available microarray dataset of breast cancer patients that have notoriously heterogeneous phenotypes. We obtained a feature subset of 13 microRNAs that could classify the 71 breast cancer patients into five subtypes with significantly different five-year overall survival rates (45%, 82.4%, 70.6%, 100% and 60% respectively; p = 0.008). By building a multivariate Cox proportional-hazards prediction model for the feature subset, we identified has-miR-146b as one of the most significant predictor (p = 0.045; hazard ratios = 0.39). The proposed algorithm is a promising computational strategy for dissecting hidden genetic heterogeneity for complex diseases, and will be of value for improving cancer diagnosis and treatment.

Tetrahymena thermophila is a widely used unicellular eukaryotic model organism in biological research and contains more than 1000 protein kinases and phosphatases with specificity for Ser/Thr/Tyr residues. However, only a few dozen phosphorylation sites in T. thermophila are known, presenting a major obstacle to further understanding of the regulatory roles of reversible phosphorylation in this organism. In this study, we used high-accuracy mass-spectrometry-based proteomics to conduct global and site-specific phosphoproteome profiling of T. thermophila. In total, 1384 phosphopeptides and 2238 phosphorylation sites from 1008 T. thermophila proteins were identified through the combined use of peptide prefractionation, TiO2 enrichment, and two-dimensional LC-MS/MS analysis. The identified phosphoproteins are implicated in the regulation of various biological processes such as transport, gene expression, and mRNA metabolic process. Moreover, integrated analysis of the T. thermophila phosphoproteome and gene network revealed the potential biological functions of many previously unannotated proteins and predicted some putative kinase-substrate pairs. Our data provide the first global survey of phosphorylation in T. thermophila using a phosphoproteomic approach and suggest a wide-ranging regulatory scope of this modification. The provided dataset is a valuable resource for the future understanding of signaling pathways in this important model organism.

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

Developmentally programmed genome rearrangement accompanies differentiation of the silent germline micronucleus into the transcriptionally active somatic macronucleus in the ciliated protozoan Tetrahymena thermophila. Internal eliminated sequences (IES) are excised, followed by rejoining of MAC-destined sequences, while fragmentation occurs at conserved chromosome breakage sequences, generating macronuclear chromosomes. Some macronuclear chromosomes, referred to as non-maintained chromosomes (NMC), are lost soon after differentiation. Large NMC contain genes implicated in development-specific roles. One such gene encodes the domesticated piggyBac transposase TPB6, required for heterochromatin-dependent precise excision of IES residing within exons of functionally important genes. These conserved exonic IES determine alternative transcription products in the developing macronucleus; some even contain free-standing genes. Examples of precise loss of some exonic IES in the micronucleus and retention of others in the macronucleus of related species suggest an evolutionary analogy to introns. Our results reveal that germline-limited sequences can encode genes with specific expression patterns and development-related functions, which may be a recurring theme in eukaryotic organisms experiencing programmed genome rearrangement during germline to soma differentiation.

Nonsense-mediated mRNA decay (NMD) is essential for removing premature termination codon-containing transcripts from cells. Studying the NMD pathway in model organisms can help to elucidate the NMD mechanism in humans and improve our understanding of how this biologically important process has evolved. Ciliates are among the earliest branching eukaryotes; their NMD mechanism is poorly understood and may be primordial. We demonstrate that highly conserved Upf proteins (Upf1a, Upf2 and Upf3) are involved in the NMD pathway of the ciliate, Tetrahymena thermophila. We further show that a novel protozoa-specific nuclease, Smg6L, is responsible for destroying many NMD-targeted transcripts. Transcriptome-wide identification and characterization of NMD-targeted transcripts in vegetative Tetrahymena cells showed that many have exon-exon junctions downstream of the termination codon. However, Tetrahymena may lack a functional exon junction complex (EJC), and the Tetrahymena ortholog of an EJC core component, Mago nashi (Mag1), is dispensable for NMD. Therefore, NMD is EJC independent in this early branching eukaryote.

Within the past decade, genomic studies have emerged as essential and highly productive tools to explore the biology of Tetrahymena thermophila. The current major resources, which have been extensively mined by the research community, are the annotated macronuclear genome assembly, transcriptomic data and the databases that house this information. Efforts in progress will soon improve these data sources and expand their scope, including providing annotated micronuclear and comparative genomic sequences. Future studies of Tetrahymena cell and molecular biology, development, physiology, evolution and ecology will benefit greatly from these resources and the advanced genomic technologies they enable.

The model eukaryote, Tetrahymena thermophila, is the first ciliated protozoan whose genome has been sequenced, enabling genome-wide analysis of gene expression.

Cytochrome P450 monooxygenases play key roles in the metabolism of a wide variety of substrates and they are closely associated with endocellular physiological processes or detoxification metabolism under environmental exposure. To date, however, none has been systematically characterized in the phylum Ciliophora. T. thermophila possess many advantages as a eukaryotic model organism and it exhibits rapid and sensitive responses to xenobiotics, making it an ideal model system to study the evolutionary and functional diversity of the P450 monooxygenase gene family.

The cell wall and chloroplast are two fundamental structures determining plant mechanical strength and grain yield. Therefore, understanding mechanisms that improve plants' ability to develop a robust cell wall and well-developed chloroplast is of utmost importance for agricultural activities.

The present study aimed to determine the inhibitory effects of angiotensin II (AngII) on angiopoietin‑like protein 2 (Angptl2) in rat primary cardiomyocytes, and to investigate the potential association between angiotensin II type 1 receptor (AT1R) and these effects. Cardiomyocytes were isolated from 3-day-old Wistar rats, and were cultured and identified. Subsequently, the expression levels of Angptl2 were detected following incubation with various concentrations of AngII for various durations using western blotting, reverse transcription‑quantitative polymerase chain reaction, enzyme-linked immunosorbent assay and immunofluorescence. Finally, under the most appropriate conditions (100 nmol/l AngII, 24 h), the cardiomyocytes were divided into six groups: Normal, AngII, AngII + losartan, normal + losartan, AngII + PD123319 and normal + PD123319 groups, in order to investigate the possible function of AT1R in Angptl2 suppression. Losartan and PD123319 are antagonists of AT1R and angiotensin II type 2 receptor, respectively. The statistical significance of the results was analyzed using Student's t‑test or one‑way analysis of variance. The results demonstrated that Angptl2 expression was evidently suppressed (P<0.05) following incubation with 100 nmol/l AngII for 24 h. Conversely, the expression levels of Angptl2 were significantly increased in the AngII + losartan group compared with the AngII group (P<0.01). However, no significant difference was detected between the AngII + PD123319, normal + losartan or normal + PD123319 groups and the normal group. The present in vitro study indicated that AngII was able to suppress Angptl2 expression, whereas losartan was able to significantly reverse this decrease by inhibiting AT1R.

Co-culture systems of rice and aquatic animals can contribute to the ecological intensification of agriculture by reducing nutrient loss and the need for N fertilizer application and by enhancing nutrient-use efficiency. However, the input of high-protein diets into paddy fields, to facilitate the growth of aquatic animals, has been found to increase N pollution and acidification of the soil. Although soil amendments have been widely used to ameliorate acidic soils, reduce N2O emissions, and improve agronomic production, the relationship between soil amendments and aquatic animal remains unclear. Therefore, this study investigated the effects of calcined dolomite (hereafter referred to as dolomite) as an acidic soil amendment and Ca-Mg supplement in rice-crab co-culture using Eriocheir sinensis crabs (Chinese mitten crabs). High-throughput sequencing was used to examine crab bacterial community composition and crab hepatopancreas biology. Although the water pH was significantly increased in the dolomite group, the number, composition, and diversity of bacteria identified in crab gut microbiome did not vary significantly between the dolomite and control groups. In the dolomite group, the probiotic agents Candidatus Hepatoplasma and Lactobacillus were highly abundant in the crab gut, and immune- and retinol metabolism-related genes were significantly upregulated in the crab hepatopancreas. Overall, dolomite application increased crab health and water pH. Dolomite is a low-cost amendment, with better stability, compared to other soil amendments, thus making it ideal for sustainable and clean rice-aquatic animal co-culture.

Ciliates are a large and diverse group of unicellular organisms characterized by having the following two distinct type of nuclei within a single cell: micronucleus (MIC) and macronucleus (MAC). Although the genomes of several ciliates in different groups have been sequenced, comparative genomics data for multiple species within a ciliate genus are not yet available. Here we collected the genome information and comparative genomics analysis results for 10 species in the Tetrahymena genus, including the previously sequenced model organism Tetrahymena thermophila and 9 newly sequenced species, and constructed a genus-level comparative analysis platform, the Tetrahymena Comparative Genomics Database (TCGD). Genome sequences, transcriptomic data, gene models, functional annotation, ortholog groups and synteny maps were built into this database and a user-friendly interface was developed for searching, visualizing and analyzing these data. In summary, the TCGD (http://ciliate.ihb.ac.cn) will be an important and useful resource for the ciliate research community.