Properties of non-coding DNA and identification of putative cis-regulatory elements in Theileria parva

AraC transcriptional regulator, aspartate semialdehyde dehydrogenase and acyltransferase: Three putative genes in phenol catabolic pathway of Acinetobacter sp. Strain DF4

The objective of this study was to identify genes associated with the biodegradation of phenol by Acinetobacter sp. strain DF4 through the use of differential display (DD) methodology. The bacteria were grown in YEPG medium, and total RNA was extracted and analyzed using labeled primers to detect gene expression differences. Three distinctively expressed cDNA bands (ph1, ph2, and ph3) were identified, cloned, and sequenced. DNA analysis involved searching for open reading frames (ORFs), verifying results with the NCBI database, predicting promoter regions, and constructing phylogenetic trees using bioinformatics tools. The ph1 gene displayed a 97% identity with the AraC transcriptional regulator, suggesting its potential role in regulating the ortho-catabolic pathway of phenol. The ph2 gene showed a 98% identity with aspartate semialdehyde dehydrogenase, which is involved in phenol degradation. The ph3 gene had a 93% identity with acetyltransferase. Essential transcription factors, such as TATA, GTGTGT, CACA, and CTTTT, were detected, and the three genes promoter regions were predicted. This study successfully identified functional genes involved in the metabolism of cyclic chemicals, particularly phenol, using the DD technique. These findings provide insights into the biodegradation pathways of phenol by Acinetobacter sp. Strain DF4 and may contribute to the development of more efficient bioremediation strategies for phenol-contaminated environments.

Transcriptomics for Drought Stress Mediated by Biological Processes in-relation to Key Regulated Pathways in Gossypium darwinii

BACKGROUND

Wild cotton Gossypium darwinii, an allotetraploid harbours important traits useful for tolerating abiotic stress, i.e., drought, salt and good genetic stability, hence these characteristics can be transferred to cultivated cotton for genetic improvement.

MATERIALS AND METHODS

In this study, we analyzed the RNA-seq transcriptomes from leaves of G. darwinii seedlings with and without drought stress. A total of 86.7 million valid reads with an average length of 95.79 bp were generated from the two samples and 58,960 transcripts with a length of more than 500 bp were assembled. We searched the known proteins on the strength of sequence similarity; these transcripts were annotated with COG, KEGG and GO functional categories. According to gene expression abundance RPKM value, we carried out RT-qPCR analysis to determine the expression pattern of the obtained transcription factors.

RESULTS

A total of 58,960 genes was differentially expressed (DEG), with 32,693 and 25,919 genes found to be upregulated and downregulated, respectively. Through gene ontology and KEGG pathways, the upregulated genes were found to associate with all the GO terms, molecular functions (MF), biological process (BP) and cellular components (CC), which are highly linked to enhancing drought stress tolerance.

CONCLUSION

The study provides an in-depth knowledge of regulation of pathways and genes involved in photosynthesis during drought stress in G. darwinii. These pathways and genes were found to be significantly downregulated and this information could be further utilized by cotton breeders in developing a more drought tolerant cotton germplasm.

A Field Guide to Eukaryotic Transposable Elements

Transposable elements (TEs) are mobile DNA sequences that propagate within genomes. Through diverse invasion strategies, TEs have come to occupy a substantial fraction of nearly all eukaryotic genomes, and they represent a major source of genetic variation and novelty. Here we review the defining features of each major group of eukaryotic TEs and explore their evolutionary origins and relationships. We discuss how the unique biology of different TEs influences their propagation and distribution within and across genomes. Environmental and genetic factors acting at the level of the host species further modulate the activity, diversification, and fate of TEs, producing the dramatic variation in TE content observed across eukaryotes. We argue that cataloging TE diversity and dissecting the idiosyncratic behavior of individual elements are crucial to expanding our comprehension of their impact on the biology of genomes and the evolution of species.

Cis regulatory motifs and antisense transcriptional control in the apicomplexan Theileria parva

Background

Theileria parva is an intracellular parasite that causes a lymphoproliferative disease in cattle. It does so by inducing cancer-like phenotypes in the host cells it infects, although the molecular and regulatory mechanisms involved remain poorly understood. RNAseq data, and the resulting updated genome annotation now available for this parasite, offer an unprecedented opportunity to characterize the genomic features associated with gene regulation in this species. Our previous analyses revealed a T. parva genome even more gene-dense than previously thought, with many adjacent loci overlapping each other, not only at the level of untranslated sequences (UTRs) but even in coding sequences.

Results

Despite this compactness, Theileria intergenic regions show a pattern of size distribution indicative of monocistronic gene transcription. Three previously described motifs are conserved among Theileria species and highly prevalent in promoter regions near or at the transcription start sites. We found novel motifs at many transcription termination sites, as well as upstream of parasite genes thought to be critical for host transformation. Adjacent genes that could be regulated by antisense transcription from an overlapping transcriptional unit are syntenic between T. parva and P. falciparum at a frequency higher than expected by chance, suggesting the presence of common, and evolutionary old, regulatory mechanisms in the phylum Apicomplexa.

Conclusions

We propose a model of transcription with conserved sense and antisense transcription from a few taxonomically ubiquitous and several species-specific promoter motifs. Interestingly, the gene networks regulated by conserved promoters are themselves, in most cases, not conserved between species or genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2444-5) contains supplementary material, which is available to authorized users.

ApiAP2 Factors as Candidate Regulators of Stochastic Commitment to Merozoite Production in Theileria annulata

Background

Differentiation of one life-cycle stage to the next is critical for survival and transmission of apicomplexan parasites. A number of studies have shown that stage differentiation is a stochastic process and is associated with a point that commits the cell to a change over in the pattern of gene expression. Studies on differentiation to merozoite production (merogony) in T. annulata postulated that commitment involves a concentration threshold of DNA binding proteins and an auto-regulatory loop.

Principal Findings

In this study ApiAP2 DNA binding proteins that show changes in expression level during merogony of T. annulata have been identified. DNA motifs bound by orthologous domains in Plasmodium were found to be enriched in upstream regions of stage-regulated T. annulata genes and validated as targets for the T. annulata AP2 domains by electrophoretic mobility shift assay (EMSA). Two findings were of particular note: the gene in T. annulata encoding the orthologue of the ApiAP2 domain in the AP2-G factor that commits Plasmodium to gametocyte production, has an expression profile indicating involvement in transmission of T. annulata to the tick vector; genes encoding related domains that bind, or are predicted to bind, sequence motifs of the type 5'-(A)CACAC(A) are implicated in differential regulation of gene expression, with one gene (TA11145) likely to be preferentially up-regulated via auto-regulation as the cell progresses to merogony.

Conclusions

We postulate that the Theileria factor possessing the AP2 domain orthologous to that of Plasmodium AP2-G may regulate gametocytogenesis in a similar manner to AP2-G. In addition, paralogous ApiAP2 factors that recognise 5'-(A)CACAC(A) type motifs could operate in a competitive manner to promote reversible progression towards the point that commits the cell to undergo merogony. Factors possessing AP2 domains that bind (or are predicted to bind) this motif are present in the vector-borne genera Theileria, Babesia and Plasmodium, and other Apicomplexa; leading to the proposal that the mechanisms that control stage differentiation will show a degree of conservation.

The ability of vector-borne Apicomplexan parasites (Babesia, Plasmodium and Theileria) to change from one life-cycle stage to the next is critical for establishment of infection and transmission to new hosts. Stage differentiation steps of both Plasmodium and Theileria are known to involve stochastic transition through an intermediate form to a point that commits the cell to generate the next stage in the life-cycle. In this study we have identified genes encoding ApiAP2 DNA binding proteins in Theileria annulata that are differentially expressed during differentiation from the macroschizont stage, through merozoite production (merogony) to the piroplasm stage. The results provide evidence that the ApiAp2 factor in Theileria that possesses the orthologue of the Plasmodium AP2-G domain may also operate to regulate gametocytogenesis, and that progression to merogony is promoted by the ability of a merozoite DNA binding protein to preferentially up-regulate its own production. In addition, identification of multiple ApiAP2 DNA binding domains that bind related motifs within and across vector-borne Apicomplexan genera lead to the proposal that the mechanisms that promote the transition from asexual to sexual replication will show a degree of conservation.

The Babesia bovis gene and promoter model: an update from full-length EST analysis

BACKGROUND

Babesia bovis is an apicomplexan parasite that causes babesiosis in infected cattle. Genomes of pathogens contain promising information that can facilitate the development of methods for controlling infections. Although the genome of B. bovis is publically available, annotated gene models are not highly reliable prior to experimental validation. Therefore, we validated a preproposed gene model of B. bovis and extended the associated annotations on the basis of experimentally obtained full-length expressed sequence tags (ESTs).

RESULTS

From in vitro cultured merozoites, 12,286 clones harboring full-length cDNAs were sequenced from both ends using the Sanger method, and 6,787 full-length cDNAs were assembled. These were then clustered, and a nonredundant referential data set of 2,115 full-length cDNA sequences was constructed. The comparison of the preproposed gene model with our data set identified 310 identical genes, 342 almost identical genes, 1,054 genes with potential structural inconsistencies, and 409 novel genes. The median length of 5' untranslated regions (UTRs) was 152 nt. Subsequently, we identified 4,086 transcription start sites (TSSs) and 2,023 transcriptionally active regions (TARs) by examining 5' ESTs. We identified ATGGGG and CCCCAT sites as consensus motifs in TARs that were distributed around -50 bp from TSSs. In addition, we found ACACA, TGTGT, and TATAT sites, which were distributed periodically around TSSs in cycles of approximately 150 bp. Moreover, related periodical distributions were not observed in mammalian promoter regions.

CONCLUSIONS

The observations in this study indicate the utility of integrated bioinformatics and experimental data for improving genome annotations. In particular, full-length cDNAs with one-base resolution for TSSs enabled the identification of consensus motifs in promoter sequences and demonstrated clear distributions of identified motifs. These observations allowed the illustration of a model promoter composition, which supports the differences in transcriptional regulation frameworks between apicomplexan parasites and mammals.

Genome-wide upstream motif analysis of Cryptosporidium parvum genes clustered by expression profile

Background

There are very few molecular genetic tools available to study the apicomplexan parasite Cryptosporidium parvum. The organism is not amenable to continuous in vitro cultivation or transfection, and purification of intracellular developmental stages in sufficient numbers for most downstream molecular applications is difficult and expensive since animal hosts are required. As such, very little is known about gene regulation in C. parvum.

Results

We have clustered whole-genome gene expression profiles generated from a previous study of seven post-infection time points of 3,281 genes to identify genes that show similar expression patterns throughout the first 72 hours of in vitro epithelial cell culture. We used the algorithms MEME, AlignACE and FIRE to identify conserved, overrepresented DNA motifs in the upstream promoter region of genes with similar expression profiles. The most overrepresented motifs were E2F (5′-TGGCGCCA-3′); G-box (5′-G.GGGG-3′); a well-documented ApiAP2 binding motif (5′-TGCAT-3′), and an unknown motif (5′-[A/C] AACTA-3′). We generated a recombinant C. parvum DNA-binding protein domain from a putative ApiAP2 transcription factor [CryptoDB: cgd8_810] and determined its binding specificity using protein-binding microarrays. We demonstrate that cgd8_810 can putatively bind the overrepresented G-box motif, implicating this ApiAP2 in the regulation of many gene clusters.

Conclusion

Several DNA motifs were identified in the upstream sequences of gene clusters that might serve as potential cis-regulatory elements. These motifs, in concert with protein DNA binding site data, establish for the first time the beginnings of a global C. parvum gene regulatory map that will contribute to our understanding of the development of this zoonotic parasite.

Evolution of transcriptional regulation in closely related bacteria

BACKGROUND

The exponential growth of the number of fully sequenced genomes at varying taxonomic closeness allows one to characterize transcriptional regulation using comparative-genomics analysis instead of time-consuming experimental methods. A transcriptional regulatory unit consists of a transcription factor, its binding site and a regulated gene. These units constitute a graph which contains so-called "network motifs", subgraphs of a given structure. Here we consider genomes of closely related Enterobacteriales and estimate the fraction of conserved network motifs and sites as well as positions under selection in various types of non-coding regions.

RESULTS

Using a newly developed technique, we found that the highest fraction of positions under selection, approximately 50%, was observed in synvergon spacers (between consecutive genes from the same strand), followed by ~45% in divergon spacers (common 5'-regions), and ~10% in convergon spacers (common 3'-regions). The fraction of selected positions in functional regions was higher, 60% in transcription factor-binding sites and ~45% in terminators and promoters. Small, but significant differences were observed between Escherichia coli and Salmonella enterica. This fraction is similar to the one observed in eukaryotes.The conservation of binding sites demonstrated some differences between types of regulatory units. In E. coli, strains the interactions of the type "local transcriptional factor gene" turned out to be more conserved in feed-forward loops (FFLs) compared to non-motif interactions. The coherent FFLs tend to be less conserved than the incoherent FFLs. A natural explanation is that the former imply functional redundancy.

CONCLUSIONS

A naïve hypothesis that FFL would be highly conserved turned out to be not entirely true: its conservation depends on its status in the transcriptional network and also from its usage. The fraction of positions under selection in intergenic regions of bacterial genomes is roughly similar to that of eukaryotes. Known regulatory sites explain 20±5% of selected positions.

High-resolution genotyping and mapping of recombination and gene conversion in the protozoan Theileria parva using whole genome sequencing

BACKGROUND

Theileria parva is a tick-borne protozoan parasite, which causes East Coast Fever, a disease of cattle in sub-Saharan Africa. Like Plasmodium falciparum, the parasite undergoes a transient diploid life-cycle stage in the gut of the arthropod vector, which involves an obligate sexual cycle. As assessed using low-resolution VNTR markers, the crossover (CO) rate in T. parva is relatively high and has been reported to vary across different regions of the genome; non-crossovers (NCOs) and CO-associated gene conversions have not yet been characterised due to the lack of informative markers. To examine all recombination events at high marker resolution, we sequenced the haploid genomes of two parental strains, and two recombinant clones derived from ticks fed on cattle that had been simultaneously co-infected with two different parasite isolates.

RESULTS

By comparing the genome sequences, we were able to genotype over 64 thousand SNP markers with an average spacing of 127 bp in the two progeny clones. Previously unrecognized COs in sub-telomeric regions were detected. About 50% of CO breakpoints were accompanied by gene conversion events. Such a high fraction of COs accompanied by gene conversions demonstrated the contributions of meiotic recombination to the diversity and evolutionary success of T. parva, as the process not only redistributed existing genetic variations, but also altered allelic frequencies. Compared to COs, NCOs were more frequently observed and more uniformly distributed across the genome. In both progeny clones, genomic regions with more SNP markers had a reduced frequency of COs or NCOs, suggesting that the sequence divergence between the parental strains was high enough to adversely affect recombination frequencies. Intra-species polymorphism analysis identified 81 loci as likely to be under selection in the sequenced genomes.

CONCLUSIONS

Using whole genome sequencing of two recombinant clones and their parents, we generated maps of COs, NCOs, and CO-associated gene conversion events for T. parva. The data comprises one of the highest-resolution genome-wide analyses of the multiple outcomes of meiotic recombination for this pathogen. The study also demonstrates the usefulness of high throughput sequencing typing for detailed analysis of recombination in organisms in which conventional genetic analysis is technically difficult.

Conservation and divergence of known apicomplexan transcriptional regulons

BACKGROUND

The apicomplexans are a diverse phylum of parasites causing an assortment of diseases including malaria in a wide variety of animals and lymphoproliferation in cattle. Little is known about how these varied parasites regulate their transcriptional regulons. Even less is known about how regulon systems, consisting of transcription factors and target genes together with their associated biological process, evolve in these diverse parasites.

RESULTS

In order to obtain insights into the differences in transcriptional regulation between these parasites we compared the orthology profiles of putative malaria transcription factors across species and examined the enrichment patterns of four binding sites across eleven apicomplexans. About three-fifths of the factors are broadly conserved in several phylogenetic orders of sequenced apicomplexans. This observation suggests the existence of regulons whose regulation is conserved across this ancient phylum. Transcription factors not broadly conserved across the phylum are possibly involved in regulon systems that have diverged between species. Examining binding site enrichment patterns in light of transcription factor conservation patterns suggests a second mode via which regulon systems may diverge - rewiring of existing transcription factors and their associated binding sites in specific ways. Integrating binding sites with transcription factor conservation patterns also facilitated prediction of putative regulators for one of the binding sites.

CONCLUSIONS

Even though transcription factors are underrepresented in apicomplexans, the distribution of these factors and their associated regulons reflect common and family-specific transcriptional regulatory processes.

Role of cis-regulatory elements on the ring-specific hrp3 promoter in the human parasite Plasmodium falciparum

Identification of promoter elements responsible for regulation of gene expression has been hampered by the AT richness of P. falciparum intergenic regions. Nested deletions of histidine-rich protein 3 (hrp3) promoter suggested the presence of a multipartite ring-specific element. Linker scanning (LS) of this ring-specific promoter showed that the alteration of several promoter regions decreased the luciferase activity compared to the wild-type configuration, indicating that these regions played a role in gene expression. No homology was observed by comparison of putative regulatory elements of other genes identified by bioinformatic analysis with the hrp3 enhancer, implying a different mechanism of gene regulation by the hrp3 promoter. LS and deletion analysis of the 5' untranslated region (UTR) of the hrp3 suggested that this region contains elements which interact with promoter elements to regulate gene expression. Analysis of the intron in the UTR region suggested that this region does not play a role in stage specificity in the hrp3 promoter. Together, our results indicate the presence of multiple mechanisms of gene regulation in the parasite.