Phylogenetic position and codon usage of two centrin genes from the rumen ciliate protozoan, Entodinium caudatum

Codon Usage Bias Analysis in Macronuclear Genomes of Ciliated Protozoa

Ciliated protozoa (ciliates) are unicellular eukaryotes, several of which are important model organisms for molecular biology research. Analyses of codon usage bias (CUB) of the macronuclear (MAC) genome of ciliates can promote a better understanding of the genetic mode and evolutionary history of these organisms and help optimize codons to improve gene editing efficiency in model ciliates. In this study, the following indices were calculated: the guanine-cytosine (GC) content, the frequency of the nucleotides at the third position of codons (T3, C3, A3, G3), the effective number of codons (ENc), GC content at the 3rd position of synonymous codons (GC3s), and the relative synonymous codon usage (RSCU). Parity rule 2 plot analysis, Neutrality plot analysis, ENc plot analysis, and correlation analysis were employed to explore the main influencing factors of CUB. The results showed that the GC content in the MAC genomes of each of 21 ciliate species, the genomes of which were relatively complete, was lower than 50%, and the base compositions of GC and GC3s were markedly distinct. Synonymous codon analysis revealed that the codons in most of the 21 ciliates ended with A or T and four codons were the general putative optimal codons. Collectively, our results indicated that most of the ciliates investigated preferred using the codons with anof AT-ending and that codon usage bias was affected by gene mutation and natural selection.

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.

Rumen Protozoa Play a Significant Role in Fungal Predation and Plant Carbohydrate Breakdown

The rumen protozoa, alongside fungi, comprise the eukaryotic portion of the rumen microbiome. Rumen protozoa may account for up to 50% of biomass, yet their role in this ecosystem remains unclear. Early experiments inferred a role in carbohydrate and protein metabolism, but due to their close association with bacteria, definitively attributing these functions to the protozoa was challenging. The advent of ‘omic technologies has created opportunities to broaden our understanding of the rumen protozoa. This study aimed to utilize these methods to further our understanding of the role that protozoa play in the rumen in terms of their metabolic capacities, and in doing so, contribute valuable sequence data to reduce the chance of mis or under-representation of the rumen protozoa in meta’omic datasets. Rumen protozoa were isolated and purified using glucose-based sedimentation and differential centrifugation, extracted RNA was Poly(A) fraction enriched and DNase treated before use in a phage-based, cDNA metatranscriptomic library. Biochemical activity testing of the phage library showed 6 putatively positive plaques in response to carboxymethyl cellulose agar (indicative of cellulose activity), and no positive results for tributyrin (indicative of esterase/lipase activity) or egg yolk agar (indicative of proteolysis). Direct sequencing of the cDNA was also conducted using the Illumina HiSeq 2500. The metatranscriptome identified a wealth of carbohydrate-active enzymes which accounted for 8% of total reads. The most highly expressed carbohydrate-active enzymes were glycosyl hydrolases 5 and 11, polysaccharide lyases and deacetylases, xylanases and enzymes active against pectin, mannan and chitin; the latter likely used to digest rumen fungi which contain a chitin-rich cell membrane. Codon usage analysis of expressed genes also showed evidence of horizontal gene transfer, suggesting that many of these enzymes were acquired from the rumen bacteria in an evolutionary response to the carbohydrate-rich environment of the rumen. This study provides evidence of the significant contribution that the protozoa make to carbohydrate breakdown in the rumen, potentially using horizontally acquired genes, and highlights their predatory capacity.

Single-cell transcriptome sequencing of rumen ciliates provides insight into their molecular adaptations to the anaerobic and carbohydrate-rich rumen microenvironment

Rumen ciliates are a specialized group of ciliates exclusively found in the anaerobic, carbohydrate-rich rumen microenvironment. However, the molecular and mechanistic basis of the physiological and behavioral adaptation of ciliates to the rumen microenvironment is undefined. We used single-cell transcriptome sequencing to explore the adaptive evolution of three rumen ciliates: two entodiniomorphids, Entodinium furca and Diplodinium dentatum; and one vestibuliferid, Isotricha intestinalis. We found that all three species are members of monophyletic orders within the class Litostomatea, with E. furca and D. dentatum in Entodiniomorphida and I. intestinalis in Vestibuliferida. The two entodiniomorphids might use H₂-producing mitochondria and the vestibuliferid might use anaerobic mitochondria to survive under strictly anaerobic conditions. Moreover, carbohydrate-active enzyme (CAZyme) genes were identified in all three species, including cellulases, hemicellulases, and pectinases. The evidence that all three species have acquired prokaryote-derived genes by horizontal gene transfer (HGT) to digest plant biomass includes a significant enrichment of gene ontology categories such as cell wall macromolecule catabolic process and carbohydrate catabolic process and the identification of genes in common between CAZyme and HGT groups. These findings suggest that HGT might be an important mechanism in the adaptive evolution of ciliates to the rumen microenvironment.

Characterisation of the effect of day length, and associated differences in dietary intake, on the gut microbiota of Soay sheep

Differences in the rumen bacterial community have been previously reported for Soay sheep housed under different day length conditions. This study extends this previous investigation to other organs of the digestive tract, as well as the analysis of ciliated protozoa and anaerobic fungi. The detectable concentrations of ciliated protozoa and anaerobic fungi decreased with increased day length in both the rumen and large colon, unlike those of bacteria where no effect was observed. Conversely, bacterial community composition was affected by day length in both the rumen and large colon, but the community composition of the detectable ciliated protozoa and anaerobic fungi was not affected. Day length-associated differences in the bacterial community composition extended to all of the organs examined, with the exception of the duodenum and the jejunum. It is proposed that differences in rumen fill and ruminal ‘by-pass’ nutrients together with endocrinological changes cause the observed effects of day length on the different gut microbial communities.

The Role of Ciliate Protozoa in the Rumen

First described in 1843, Rumen protozoa with their striking appearance were assumed to be important for the welfare of their host. However, despite contributing up to 50% of the bio-mass in the rumen, the role of protozoa in rumen microbial ecosystem remains unclear. Phylogenetic analysis of 18S rDNA libraries generated from the rumen of cattle, sheep, and goats has revealed an unexpected diversity of ciliated protozoa although variation in gene copy number between species makes it difficult to obtain absolute quantification. Despite repeated attempts it has proven impossible to maintain rumen protozoa in axenic culture. Thus it has been difficult to establish conclusively a role of ciliate protozoa in rumen fiber degradation. The development of techniques to clone and express ciliate genes in λ phage, together with bioinformatic indices to confirm the ciliate origin of the genes has allowed the isolation and characterization of fibrolytic genes from rumen protozoa. Elimination of the ciliate protozoa increases microbial protein supply by up to 30% and reduces methane production by up to 11%. Our recent findings suggest that holotrich protozoa play a disproportionate role in supporting methanogenesis whilst the small Entodinium are responsible for much of the bacterial protein turnover. As yet no method to control protozoa in the rumen that is safe and practically applicable has been developed, however a range of plant extract capable of controlling if not completely eliminating rumen protozoa have been described.

Activity-based metagenomic screening and biochemical characterization of bovine ruminal protozoan glycoside hydrolases

The rumen, the foregut of herbivorous ruminant animals such as cattle, functions as a bioreactor to process complex plant material. Among the numerous and diverse microbes involved in ruminal digestion are the ruminal protozoans, which are single-celled, ciliated eukaryotic organisms. An activity-based screen was executed to identify genes encoding fibrolytic enzymes present in the metatranscriptome of a bovine ruminal protozoan-enriched cDNA expression library. Of the four novel genes identified, two were characterized in biochemical assays. Our results provide evidence for the effective use of functional metagenomics to retrieve novel enzymes from microbial populations that cannot be maintained in axenic cultures.

The host species affects the microbial community in the goat rumen

AIMS

This study was carried out to determine whether bacterial and ciliate populations in goat rumen vary significantly between different goat species living in the same environment.

METHODS AND RESULTS

Bacterial and ciliate communities in the rumen of three goat species were analysed at the molecular level using denaturing gradient gel electrophoresis. The microbial community varied considerably among goats living in the same environment. Interspecies variation in the bacterial population was noticeably greater than intraspecies variation. In contrast, there was considerable variation in the ciliate population among goats within the same species, and intraspecies similarities were no greater than those observed across species.

CONCLUSIONS

Because environmental factors and diets were identical for all goats, differences in bacterial populations reflect species-specific differences in rumen microbes.

SIGNIFICANCE AND IMPACT OF THE STUDY

Factors related to the host species have an important effect on determining the bacterial composition in the goat rumen.

Horizontal gene transfer from Bacteria to rumen Ciliates indicates adaptation to their anaerobic, carbohydrates-rich environment

Background

The horizontal transfer of expressed genes from Bacteria into Ciliates which live in close contact with each other in the rumen (the foregut of ruminants) was studied using ciliate Expressed Sequence Tags (ESTs). More than 4000 ESTs were sequenced from representatives of the two major groups of rumen Cilates: the order Entodiniomorphida (Entodinium simplex, Entodinium caudatum, Eudiplodinium maggii, Metadinium medium, Diploplastron affine, Polyplastron multivesiculatum and Epidinium ecaudatum) and the order Vestibuliferida, previously called Holotricha (Isotricha prostoma, Isotricha intestinalis and Dasytricha ruminantium).

Results

A comparison of the sequences with the completely sequenced genomes of Eukaryotes and Prokaryotes, followed by large-scale construction and analysis of phylogenies, identified 148 ciliate genes that specifically cluster with genes from the Bacteria and Archaea. The phylogenetic clustering with bacterial genes, coupled with the absence of close relatives of these genes in the Ciliate Tetrahymena thermophila, indicates that they have been acquired via Horizontal Gene Transfer (HGT) after the colonization of the gut by the rumen Ciliates.

Conclusion

Among the HGT candidates, we found an over-representation (>75%) of genes involved in metabolism, specifically in the catabolism of complex carbohydrates, a rich food source in the rumen. We propose that the acquisition of these genes has greatly facilitated the Ciliates' colonization of the rumen providing evidence for the role of HGT in the adaptation to new niches.

An NAD+-dependent glutamate dehydrogenase cloned from the ruminal ciliate protozoan,Entodinium caudatum

An NAD(+)-dependent glutamate dehydrogenase (GDH; EC 1.4.1.24) was cloned from the ruminal ciliate protozoan, Entodinium caudatum. The gene had high sequence similarity to GDH genes from the Bacteroides (class)--a class of bacteria which is highly represented in the rumen. When expressed in Escherichia coli the enzyme had a high affinity for ammonia and alpha-ketoglutarate (apparent K(m) of 2.33 and 0.71 mM, respectively) and a low affinity for glutamate (apparent K(m) of 98 mM). GDH activity and GDH mRNA concentration were increased by incubating washed E. caudatum cells with ammonia and antibiotics. These results suggest that the GDH is an anabolic enzyme catalysing the assimilation of ammonia by E. caudatum in the rumen and that the gene was probably acquired by lateral gene transfer from a ruminal bacterium.

Xylanases and carboxymethylcellulases of the rumen protozoaPolyplastron multivesiculatum,Eudiplodinium maggiiandEntodiniumsp

Endoglucanase and xylanase activities of three rumen protozoa, Polyplastron multivesiculatum, Eudiplodinium maggii, and Entodinium sp. were compared qualitatively by zymograms and quantitatively by measuring specific activities against different polysaccharides. A set of carboxymethylcellulases and xylanases was produced by the large ciliates whereas no band of activity was observed for Entodinium sp. in zymograms. Specific activity of endoglucanases from P. multivesiculatum (1.3 micromol mg prot(-1) min(-1)) was twice that of E. maggii, whereas xylanase specific activity (4.5 micromol mg prot(-1) min(-1)) was only half. Very weak activities were observed for Entodinium sp. A new xylanase gene, xyn11D, from P. multivesiculatum was reported and its gene product compared to 33 other family 11 xylanases. Phylogenetic analysis showed that xylanase sequences from rumen protozoa are closely related to those of bacteria.

Characterization of XYN10B, a modular xylanase from the ruminal protozoan Polyplastron multivesiculatum, with a family 22 carbohydrate-binding module that binds to cellulose

A new xylanase gene, xyn10B, was isolated from the ruminal protozoan Polyplastron multivesiculatum and the gene product was characterized. XYN10B is the first protozoan family 10 glycoside hydrolase characterized so far and is a modular enzyme comprising a family 22 carbohydrate-binding module (CBM) preceding the catalytic domain. The CBM22 was shown to be a true CBM. It showed high affinity for soluble arabinoxylan and is the first example of a CBM22 that binds strongly to celluloses of various crystallinities. The enzymic properties of XYN10B were also analysed. Its optimal temperature and pH for activity were 39 degrees C and 7.0 respectively; these values being close to those of the ruminal ecosystem. The phylogenetic relationships between the XYN10B CBM22 or catalytic domain and related sequences from ruminal and non-ruminal bacteria and eukaryotes are reported. The xyn10B gene is shown to lack introns.

The 3' untranslated region of messages in the rumen protozoan Entodinium caudatum

The 3' untranslated regions of a number of cDNAs from the rumen protozoal species Entodinium caudatum were studied with a view to characterising their preference for stop codons, general length, nucleotide composition and polyadenylation signals. Unlike a number of ciliates, Entodinium caudatum uses UAA as a stop codon, rather than as a codon for glutamine. In addition, the 3' untranslated region of the message is generally less than 100 nucleotides in length, extremely A+T rich, and does not appear to utilise any of the conventional polyadenylation signals described in other organisms.

A xylanase produced by the rumen anaerobic protozoan Polyplastron multivesiculatum shows close sequence similarity to family 11 xylanases from gram-positive bacteria

We report for the first time the cloning and characterisation of a protozoal enzyme involved in plant cell wall polysaccharide degradation. A cDNA library was constructed from the ruminal protozoan Polyplastron multivesiculatum and a stable clone expressing xylanase activity was isolated. The encoded enzyme belongs to the glycoside hydrolase family 11, and phylogenetic analysis indicates a closer relationship with catalytic domains from Gram-positive bacteria than the other fibrolytic eukaryotes from the rumen, the anaerobic fungi.