Pseudogene-derived small interference RNAs regulate gene expression in African Trypanosoma brucei

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

Pseudogenes in Cancer: State of the Art

Pseudogenes are duplicates of protein-coding genes that have accumulated multiple detrimental alterations, rendering them unable to produce the protein they encode. Initially disregarded as "junk DNA" due to their perceived lack of functionality, research on their biological roles has been hindered by this assumption. Nevertheless, recent focus has shifted towards these molecules due to their abnormal expression in cancer phenotypes. In this review, our objective is to provide a thorough overview of the current understanding of pseudogene formation, the mechanisms governing their expression, and the roles they may play in promoting tumorigenesis.

Enhanced genome-wide association reveals the role of YABBY11-NGATHA-LIKE1 in leaf serration development of Populus

Leaf margins are complex plant morphological features that contribute to leaf shape diversity, which affects plant structure, yield, and adaptation. Although several leaf margin regulators have been identified to date, the genetic basis of their natural variation has not been fully elucidated. In this study, we profiled two distinct leaf morphology types (serrated and smooth) using the persistent homology mathematical framework (PHMF) in two poplar species (Populus tomentosa and Populus simonii, respectively). A combined genome-wide association study (GWAS) and expression quantitative trait nucleotide (eQTN) mapping were applied to create a leaf morphology control module using data from P. tomentosa and P. simonii populations. Natural variation in leaf margins was associated with YABBY11 (YAB11) transcript abundance in poplar. In P. tomentosa, PtoYAB11 carries a premature stop codon (PtoYAB11PSC), resulting in the loss of its positive regulation of NGATHA-LIKE1 (PtoNGAL-1) and RIBULOSE BISPHOSPHATE CARBOXYLASE LARGE SUBUNIT (PtoRBCL). Overexpression of PtoYAB11PSC promoted serrated leaf margins, enlarged leaves, enhanced photosynthesis, and increased biomass. Overexpression of PsiYAB11 in P. tomentosa promoted smooth leaf margins, higher stomatal density, and greater light damage repair ability. In poplar, YAB11-NGAL1 is sensitive to environmental conditions, acts as a positive regulator of leaf margin serration, and may also link environmental signaling to leaf morphological plasticity.

PAG3 promotes the differentiation of bloodstream forms in Trypanosoma brucei and reveals the evolutionary relationship among the Trypanozoon trypanosomes

Introduction

Trypanosoma brucei, T. evansi and T. equiperdum are members of the subgenus Trypanozoon and are highly similar morphologically and genetically. The main differences between these three species are their differentiation patterns in the hosts and the role of vectors in their life cycles. However, the mechanisms causing these differences are still controversial.

Methods

PAG3 gene was accessed by PCR amplification in 26 strains of Trypanozoon and sequences were then analyzed by BLAST accompanied with T. evansitype B group. RNA interference and CRISPR/Cas9 were used for revealing possible role of PAG3 in slender to stumpy transformation.

Results

The procyclin associated gene 3 (PAG3) can be found in the pleomorphicspecies, T.brucei, which undergoes differentiation of slender forms to the stumpy form. This differentiation process is crucial for transmission to the tsetse fly vector. However, a homologue of PAG3 was not detected in either T. evansi or in the majority of T. equiperdum strains which are allmonomorphic. Furthere xperiments in T. brucei demonstrated that, when PAG3 was down-regulated or absent, there was a significant reduction in the differentiation from slender to stumpy forms.

Conclusion

Therefore, we conclude that PAG3 is a key nuclear gene involved in the slender to stumpy differentiation pathway of T.brucei in the mammalian host. Loss of this gene might also offer a simple evolutionary mechanism explaining why T. evansi and some T. equiperdum have lost the ability to differentiate and have been driven to adapt to transmission cycles that by pass the tsetse vector or mechanical contact.

Uncovering Pseudogenes and Intergenic Protein-coding Sequences in TriTryps' Genomes

Trypanosomatids belong to a remarkable group of unicellular, parasitic organisms of the order Kinetoplastida, an early diverging branch of the phylogenetic tree of eukaryotes, exhibiting intriguing biological characteristics affecting gene expression (intronless polycistronic transcription, trans-splicing, and RNA editing), metabolism, surface molecules, and organelles (compartmentalization of glycolysis, variation of the surface molecules, and unique mitochondrial DNA), cell biology and life cycle (phagocytic vacuoles evasion and intricate patterns of cell morphogenesis). With numerous genomic-scale data of several trypanosomatids becoming available since 2005 (genomes, transcriptomes, and proteomes), the scientific community can further investigate the mechanisms underlying these unusual features and address other unexplored phenomena possibly revealing biological aspects of the early evolution of eukaryotes. One fundamental aspect comprises the processes and mechanisms involved in the acquisition and loss of genes throughout the evolutionary history of these primitive microorganisms. Here, we present a comprehensive in silico analysis of pseudogenes in three major representatives of this group: Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Pseudogenes, DNA segments originating from altered genes that lost their original function, are genomic relics that can offer an essential record of the evolutionary history of functional genes, as well as clues about the dynamics and evolution of hosting genomes. Scanning these genomes with functional proteins as proxies to reveal intergenic regions with protein-coding features, relying on a customized threshold to distinguish statistically and biologically significant sequence similarities, and reassembling remnant sequences from their debris, we found thousands of pseudogenes and hundreds of open reading frames, with particular characteristics in each trypanosomatid: mutation profile, number, content, density, codon bias, average size, single- or multi-copy gene origin, number and type of mutations, putative primitive function, and transcriptional activity. These features suggest a common process of pseudogene formation, different patterns of pseudogene evolution and extant biological functions, and/or distinct genome organization undertaken by those parasites during evolution, as well as different evolutionary and/or selective pressures acting on distinct lineages.

Current Status of Regulatory Non-Coding RNAs Research in the Tritryp

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

Pseudogenes: Four Decades of Discovery

A pseudogene is defined as a genomic DNA sequence that looks like a mutated or truncated version of a known functional gene. Nearly four decades since their first discovery it has been estimated that between ~12,000 and ~20,000 pseudogenes exist in the human genome. Early efforts to characterize functions for pseudogenes were unsuccessful, thus they were considered functionless relics of evolutionary selection, junk DNA or genetic fossils. Remarkably, an increasing number of pseudogenes have been reported to be expressed as RNA transcripts above and beyond levels considered accidental or spurious transcription. There is emerging evidence that some expressed pseudogene transcripts have biological functions and should be defined as a subclass of functional long noncoding RNAs (lncRNA). In this introductory chapter, I briefly summarize the history and the current knowledge of pseudogenes, and highlight the emerging functions of some pseudogenes in human biology and disease. This second iteration of Pseudogenes in Methods in Molecular Biology highlights new methodological approaches to investigate this intriguing family of lncRNAs and the extent of their biological function.

Cancer, Retrogenes, and Evolution

This review summarizes the knowledge about retrogenes in the context of cancer and evolution. The retroposition, in which the processed mRNA from parental genes undergoes reverse transcription and the resulting cDNA is integrated back into the genome, results in additional copies of existing genes. Despite the initial misconception, retroposition-derived copies can become functional, and due to their role in the molecular evolution of genomes, they have been named the “seeds of evolution”. It is convincing that retrogenes, as important elements involved in the evolution of species, also take part in the evolution of neoplastic tumors at the cell and species levels. The occurrence of specific “resistance mechanisms” to neoplastic transformation in some species has been noted. This phenomenon has been related to additional gene copies, including retrogenes. In addition, the role of retrogenes in the evolution of tumors has been described. Retrogene expression correlates with the occurrence of specific cancer subtypes, their stages, and their response to therapy. Phylogenetic insights into retrogenes show that most cancer-related retrocopies arose in the lineage of primates, and the number of identified cancer-related retrogenes demonstrates that these duplicates are quite important players in human carcinogenesis.

Methods to Study Protein-Binding to Pseudogene Transcripts

One of the most commonly described biological feature of processed pseudogenes is the ability to influence the expression of their parental coding genes. As evidenced in several studies, the high sequence similarity between these RNA pairs sets up a certain level of competition for posttranscriptional regulators, including, among others, RNA-binding proteins (RBPs). RBPs may affect, positively or negatively, the stability of bound mRNAs, so that, if an overexpressed pseudogene competes with its homologous coding gene, the downstream protein synthesis would change, with potential pathological consequences. Given these premises, a rigorous and comprehensive understanding of interactions between pseudogene-parental gene RNA pairs and RBPs could provide further insights into the biological bases of complex diseases, such as cancer, cardiovascular disease, and type 2 diabetes, identifying novel predictive and/or prognostic biomarkers.Herein, we detail easily adaptable protocols of plasmid-based molecular cloning and RNA-electrophoretic mobility shift assay (EMSA) used in our laboratory for determining the interaction between a cytoplasmatic stabilizing protein (αCP1) and the pseudogene-parental gene RNA pair HMGA1-p /HMGA1. We also offer a general overview of RNA immunoprecipitation procedures and present novel bioinformatic tools for predicting RBPs binding sites on pseudogene transcripts.

Genomic Organization and Generation of Genetic Variability in the RHS (Retrotransposon Hot Spot) Protein Multigene Family in Trypanosoma cruzi

Retrotransposon Hot Spot (RHS) is the most abundant gene family in Trypanosoma cruzi, with unknown function in this parasite. The aim of this work was to shed light on the organization and expression of RHS in T. cruzi. The diversity of the RHS protein family in T. cruzi was demonstrated by phylogenetic and recombination analyses. Transcribed sequences carrying the RHS domain were classified into ten distinct groups of monophyletic origin. We identified numerous recombination events among the RHS and traced the origins of the donors and target sequences. The transcribed RHS genes have a mosaic structure that may contain fragments of different RHS inserted in the target sequence. About 30% of RHS sequences are located in the subtelomere, a region very susceptible to recombination. The evolution of the RHS family has been marked by many events, including gene duplication by unequal mitotic crossing-over, homologous, as well as ectopic recombination, and gene conversion. The expression of RHS was analyzed by immunofluorescence and immunoblotting using anti-RHS antibodies. RHS proteins are evenly distributed in the nuclear region of T. cruzi replicative forms (amastigote and epimastigote), suggesting that they could be involved in the control of the chromatin structure and gene expression, as has been proposed for T. brucei.

Identification of Potential Diagnostic and Prognostic Pseudogenes in Hepatocellular Carcinoma Based on Pseudogene-miRNA-mRNA Competitive Network

Background

It is widely known that hepatocellular carcinoma (HCC) has high rates of morbidity and mortality. A large number of studies have indicated that pseudogenes have an important effect on the carcinogenesis of HCC. Pseudogenes can play a role through the ceRNA network. There have been numerous studies on lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA networks. However, the pseudogene-miRNA-mRNA network in HCC has rarely been researched or reported on.

Material/Methods

The Cancer Genome Atlas (TCGA) database was researched and differences between selected genes were studied. A pseudogene-miRNA-mRNA network was then constructed and clustering of pseudogenes was studied. The diagnostic value of the selected pseudogenes, their functions, and pathways were investigated using available databases to understand their possible pathogenic mechanism in HCC. The protein-protein interaction network of target genes was found and the top 10 hub genes were identified. Expression of hub genes in HCC tissues was then detected by RT-qPCR.

Results

By analyzing the gene difference and clinical data of HCC, we constructed a ceRNA network composed of 4 pseudogenes, 8 miRNAs, and 30 mRNAs. The pseudogenes AP000769.1, KRT16P1, KRT16P3, and RPLP0P2 were all correlated with the diagnosis and prognosis of HCC. Functional analyses through the Kyoto Encyclopedia of Genes and Genomes and the Gene Ontology databases indicated that pseudogenes can affect the physiological process of HCC through the p53 pathway. The top 10 hub genes identified were all highly expressed in HCC tissues and affected the patient survival rate.

Conclusions

In this study, 4 pseudogenes related to the diagnosis and prognosis of liver cancer were found through the construction of a ceRNA network. These 4 pseudogenes might constitute new therapeutic targets for liver cancer patients.

Re-recognition of pseudogenes: From molecular to clinical applications

Pseudogenes were initially regarded as "nonfunctional" genomic elements that did not have protein-coding abilities due to several endogenous inactivating mutations. Although pseudogenes are widely expressed in prokaryotes and eukaryotes, for decades, they have been largely ignored and classified as gene "junk" or "relics". With the widespread availability of high-throughput sequencing analysis, especially omics technologies, knowledge concerning pseudogenes has substantially increased. Pseudogenes are evolutionarily conserved and derive primarily from a mutation or retrotransposon, conferring the pseudogene with a "gene repository" role to store and expand genetic information. In contrast to previous notions, pseudogenes have a variety of functions at the DNA, RNA and protein levels for broadly participating in gene regulation to influence the development and progression of certain diseases, especially cancer. Indeed, some pseudogenes have been proven to encode proteins, strongly contradicting their "trash" identification, and have been confirmed to have tissue-specific and disease subtype-specific expression, indicating their own value in disease diagnosis. Moreover, pseudogenes have been correlated with the life expectancy of patients and exhibit great potential for future use in disease treatment, suggesting that they are promising biomarkers and therapeutic targets for clinical applications. In this review, we summarize the natural properties, functions, disease involvement and clinical value of pseudogenes. Although our knowledge of pseudogenes remains nascent, this field deserves more attention and deeper exploration.

Extensive profiling in Arabidopsis reveals abundant polysome-associated 24-nt small RNAs including AGO5-dependent pseudogene-derived siRNAs

In a reductionist perspective, plant silencing small (s)RNAs are often classified as mediating nuclear transcriptional gene silencing (TGS) or cytosolic posttranscriptional gene silencing (PTGS). Among the PTGS diagnostics is the association of AGOs and their sRNA cargos with the translation apparatus. In Arabidopsis, this is observed for AGO1 loaded with micro(mi)RNAs and, accordingly, translational-repression (TR) is one layer of plant miRNA action. Using AGO1:miRNA-mediated TR as a paradigm, we explored, with two unrelated polysome-isolation methods, which, among the ten Arabidopsis AGOs and numerous sRNA classes, interact with translation. We found that representatives of all three AGO-clades associate with polysomes, including the TGS-effector AGO4 and stereotypical 24-nt sRNAs that normally mediate TGS of transposons/repeats. Strikingly, approximately half of these annotated 24-nt siRNAs displayed unique matches in coding regions/introns of genes, and in pseudogenes, but not in transposons/repeats commonly found in their vicinity. Protein-coding gene-derived 24-nt sRNAs correlate with gene-body methylation. Those derived from pseudogenes belong to two main clusters defined by their parental-gene expression patterns, and are vastly enriched in AGO5, itself found on polysomes. Based on their tight expression pattern in developing and mature siliques, their biogenesis, and genomic/epigenomic features of their loci-of-origin, we discuss potential roles for these hitherto unknown polysome-enriched, pseudogene-derived siRNAs.

lncRNA-PDPK2P promotes hepatocellular carcinoma progression through the PDK1/AKT/Caspase 3 pathway

Hepatocellular carcinoma (HCC) is a malignancy with one of the worst prognoses. Long noncoding RNA (lncRNA) are emerging as an important regulator of gene expression and function, leading to the development of cancer. The aim of this study was to determine the relationship between lncRNA and HCC and to further guide clinical therapy. lncRNA in HCC and adjacent tissues were screened, and the correlation between lncRNA-PDPK2P expression in liver tissues and the pathological characteristics and severity of HCC was assessed. The effects of PDPK2P on HCC proliferation, apoptosis, metastasis, and invasion were also systematically investigated via CCK-8 assay, flow cytometry, scratch wound healing, and transwell assay, respectively. The relationship between PDPK2P and PDK1 was verified by RNA pull-down, rescue experiments and western blot. lncRNA-PDPK2P was highly expressed in HCC tissues with a distinct positive correlation between PDPK2P and PDK1, and the upregulation was clinically associated with a larger tumor embolus, low differentiation, and poor survival. Mechanistically, lncRNA-PDPK2P interacted with PDK1 and promoted HCC progression through the PDK1/AKT/caspase 3 signaling pathway. lncRNA-PDPK2P can promote HCC progression, suggesting it may be a clinically valuable biomarker and serve as a molecular target for the diagnosis, prognosis, and therapy of hepatocellular carcinoma.

dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease

Although thousands of pseudogenes have been annotated in the human genome, their transcriptional regulation, expression profiles and functional mechanisms are largely unknown. In this study, we developed dreamBase (http://rna.sysu.edu.cn/dreamBase) to facilitate the investigation of DNA modification, RNA regulation and protein binding of potential expressed pseudogenes from multidimensional high-throughput sequencing data. Based on ∼5500 ChIP-seq and DNase-seq datasets, we identified genome-wide binding profiles of various transcription-associated factors around pseudogene loci. By integrating ∼18 000 RNA-seq data, we analysed the expression profiles of pseudogenes and explored their co-expression patterns with their parent genes in 32 cancers and 31 normal tissues. By combining microRNA binding sites, we demonstrated complex post-transcriptional regulation networks involving 275 microRNAs and 1201 pseudogenes. We generated ceRNA networks to illustrate the crosstalk between pseudogenes and their parent genes through competitive binding of microRNAs. In addition, we studied transcriptome-wide interactions between RNA binding proteins (RBPs) and pseudogenes based on 458 CLIP-seq datasets. In conjunction with epitranscriptome sequencing data, we also mapped 1039 RNA modification sites onto 635 pseudogenes. This database will provide insights into the transcriptional regulation, expression, functions and mechanisms of pseudogenes as well as their roles in biological processes and diseases.

Suppressors of RNA silencing encoded by geminiviruses and associated DNA satellites

In plants, RNA silencing provides a major line of defence against viruses. This antiviral immunity involves production of virus-derived small interfering RNAs (vsiRNAs) and results in specific silencing of viruses by vsiRNAs-guided effector complexes. As a counterattack against RNA silencing, many plant viruses encode suppressors of RNA silencing called viral suppressors of RNA silencing (VSRs), which interfere with the silencing pathway by various mechanisms. This review describes various methods that are being used to characterize viral proteins for suppressor function, VSRs found in geminiviruses and associated DNA satellites and their mechanisms of action.

Evolutionary Origins of Pseudogenes and Their Association with Regulatory Sequences in Plants

Pseudogenes (Ψs), nonfunctional relatives of functional genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Evolutionary analysis provides clues to Ψ origins and effects on gene regulation. However, few systematic studies of plant Ψs have been conducted, hampering comparative analyses. Here, we examined the origin, evolution, and expression patterns of Ψs and their relationships with noncoding sequences in seven angiosperm plants. We identified ∼250,000 Ψs, most of which are more lineage specific than protein-coding genes. The distribution of Ψs on the chromosome indicates that genome recombination may contribute to Ψ elimination. Most Ψs evolve rapidly in terms of sequence and expression levels, showing tissue- or stage-specific expression patterns. We found that a surprisingly large fraction of nontransposable element regulatory noncoding RNAs (microRNAs and long noncoding RNAs) originate from transcription of Ψ proximal upstream regions. We also found that transcription factor binding sites preferentially occur in putative Ψ proximal upstream regions compared with random intergenic regions, suggesting that Ψs have conditioned genome evolution by providing transcription factor binding sites that serve as promoters and enhancers. We therefore propose that rapid rewiring of Ψ transcriptional regulatory regions is a major mechanism driving the origin of novel regulatory modules.

Pseudogene CHIAP2 inhibits proliferation and invasion of lung adenocarcinoma cells by means of the WNT pathway

Recently, pseudogenes have emerged as critical regulators in the onset of human neoplasia. Here, we present a comprehensive analysis of pseudogene alterations at transcriptional levels in lung adenocarcinoma (LUAD) from The Cancer Genome Atlas. By combinations of differential expression analysis, survival analysis, and univariate and multivariable Cox proportional hazards regression models, we identified four dysregulated pseudogenes, whose expression level was closely related to LUAD patients' prognosis and the four pseudogene signature could act as an independent prognostic indicator for LUAD patients. We further characterized CHIAP2, one of those four pseudogenes, whose expression level was the most closely linked to LUAD patients' prognosis. Consistent with our analysis, the expression of CHIAP2 was abnormally downregulated in LUAD tissues compared with that in normal tissues in our 50 pairs of clinical samples. Functional assays demonstrated that upregulation of CHIAP2 significantly impaired cell proliferation and invasion. After performing RNA sequencing (RNA-seq) and small RNA-seq between CHIAP2 overexpression and negative control LUAD cell lines, we identified differentially expressed messenger RNAs and microRNAs (miRNAs), among which six miRNAs were downregulated. Target genes of six downregulated miRNAs were predicted with online miRNA target prediction tools and significant pathways including the WNT signal pathway were identified with Gene Set Enrichment Analysis. By combining predictor genes of six downregulated miRNAs and dysregulated genes of the WNT pathway, we inferred that overexpression of CHAP2 may inhibit LUAD cell proliferation and invasion via modulation of NFATC2 or GSK3B (WNT signal pathway) targeted by miR-3614-5p or miR-873-3p.

Origination and Function of Plant Pseudogenes

Pseudogenes, nonfunctional genomic sequences derived from functional protein-coding genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Studies on the evolution and functional aspects of plant pseudogenes are limited, despite their abundance in the plant genome. To date, most researches on pseudogenes focus on mammals. Here, we summarized current knowledge on pseudogenes including the historical and recent progress, analyzes their essential roles in gene regulation in hope of further stimulating researches in plant species for understanding gene regulation and evolution.

Role of Pseudogenes in Tumorigenesis

Functional genomics has provided evidence that the human genome transcribes a large number of non-coding genes in addition to protein-coding genes, including microRNAs and long non-coding RNAs (lncRNAs). Among the group of lncRNAs are pseudogenes that have not been paid attention in the past, compared to other members of lncRNAs. However, increasing evidence points the important role of pseudogenes in diverse cellular functions, and dysregulation of pseudogenes are often associated with various human diseases including cancer. Like other types of lncRNAs, pseudogenes can also function as master regulators for gene expression and thus, they can play a critical role in various aspects of tumorigenesis. In this review we discuss the latest developments in pseudogene research, focusing on how pseudogenes impact tumorigenesis through different gene regulation mechanisms. Given the high sequence homology with the corresponding parent genes, we also discuss challenges for pseudogene research.

Genome Plasticity in Cultured Leishmania donovani: Comparison of Early and Late Passages

Leishmania donovani possesses a complex heteroxenic life cycle where infective metacyclic promastigotes are pre-adapted to infect their host and cope up with intracellular stress. Exploiting the similarities between cultured and sandfly derived promastigotes, we used early and late passage cultured promastigotes to show specific changes at genome level which compromise pathogen fitness reflected in gene expression and infection studies. The pathogen loses virulence mostly via transcriptional and translational regulations and long-time cultivation makes them struggle to convert to virulent metacyclics. At the genomic level very subtle plasticity was observed between the early and the late passages mostly in defense-related, nutrient acquisition and signal transduction genes. Chromosome Copy number variation is seen in the early and late passages involving several genes that may be playing a role in pathogenicity. Our study highlights the importance of ABC transporters and calpain like cysteine proteases in parasite virulence in cultured promastigotes. Interestingly, these proteins are emerging as important patho-adaptive factors in clinical isolates of Leishmania. We found that the currently available genome of Leishmania in the NCBI database are from late passages. Our early passage genome can act as a reference for future studies on virulent isolates of Leishmania. The annotated leads from this study can be used for virulence surveillance and therapeutic studies in the Indian subcontinent.

HMGA1-pseudogenes and cancer

Pseudogenes are DNA sequences with high homology to the corresponding functional gene, but, because of the accumulation of various mutations, they have lost their initial functions to code for proteins. Consequently, pseudogenes have been considered until few years ago dysfunctional relatives of the corresponding ancestral genes, and then useless in the course of genome evolution. However, several studies have recently established that pseudogenes are owners of key biological functions. Indeed, some pseudogenes control the expression of functional genes by competitively binding to the miRNAs, some of them generate small interference RNAs to negatively modulate the expression of functional genes, and some of them even encode functional mutated proteins. Here, we concentrate our attention on the pseudogenes of the HMGA1 gene, that codes for the HMGA1a and HMGA1b proteins having a critical role in development and cancer progression. In this review, we analyze the family of HMGA1 pseudogenes through three aspects: classification, characterization, and their possible function and involvement in cancer.

Landscape and variation of novel retroduplications in 26 human populations

Retroduplications come from reverse transcription of mRNAs and their insertion back into the genome. Here, we performed comprehensive discovery and analysis of retroduplications in a large cohort of 2,535 individuals from 26 human populations, as part of 1000 Genomes Phase 3. We developed an integrated approach to discover novel retroduplications combining high-coverage exome and low-coverage whole-genome sequencing data, utilizing information from both exon-exon junctions and discordant paired-end reads. We found 503 parent genes having novel retroduplications absent from the reference genome. Based solely on retroduplication variation, we built phylogenetic trees of human populations; these represent superpopulation structure well and indicate that variable retroduplications are effective population markers. We further identified 43 retroduplication parent genes differentiating superpopulations. This group contains several interesting insertion events, including a SLMO2 retroduplication and insertion into CAV3, which has a potential disease association. We also found retroduplications to be associated with a variety of genomic features: (1) Insertion sites were correlated with regular nucleosome positioning. (2) They, predictably, tend to avoid conserved functional regions, such as exons, but, somewhat surprisingly, also avoid introns. (3) Retroduplications tend to be co-inserted with young L1 elements, indicating recent retrotranspositional activity, and (4) they have a weak tendency to originate from highly expressed parent genes. Our investigation provides insight into the functional impact and association with genomic elements of retroduplications. We anticipate our approach and analytical methodology to have application in a more clinical context, where exome sequencing data is abundant and the discovery of retroduplications can potentially improve the accuracy of SNP calling.

Protein-Coding Genes' Retrocopies and Their Functions

Transposable elements, often considered to be not important for survival, significantly contribute to the evolution of transcriptomes, promoters, and proteomes. Reverse transcriptase, encoded by some transposable elements, can be used in trans to produce a DNA copy of any RNA molecule in the cell. The retrotransposition of protein-coding genes requires the presence of reverse transcriptase, which could be delivered by either non-long terminal repeat (non-LTR) or LTR transposons. The majority of these copies are in a state of “relaxed” selection and remain “dormant” because they are lacking regulatory regions; however, many become functional. In the course of evolution, they may undergo subfunctionalization, neofunctionalization, or replace their progenitors. Functional retrocopies (retrogenes) can encode proteins, novel or similar to those encoded by their progenitors, can be used as alternative exons or create chimeric transcripts, and can also be involved in transcriptional interference and participate in the epigenetic regulation of parental gene expression. They can also act in trans as natural antisense transcripts, microRNA (miRNA) sponges, or a source of various small RNAs. Moreover, many retrocopies of protein-coding genes are linked to human diseases, especially various types of cancer.

Exo-miRExplorer: A Comprehensive Resource for Exploring and Comparatively Analyzing Exogenous MicroRNAs

MicroRNAs (miRNAs) are small regulatory RNAs that play important roles in animals, plants, and viruses. Deep-sequencing technology has been widely adopted in miRNA investigations. However, it is still a big mysterious why nearly all sequencing data contain miRNA sequences from exogenous species, called exo-miRNAs. In this study, we developed a novel platform, exo-miRExplorer, for mining and identifying exo-miRNAs from high-throughput small RNA sequencing experiments which originated from tissues and cell lines of multiple organisms. Thousands of exo-miRNAs are characterized with their expression abundance, the RNA families, original organisms and the sequencing platforms presented in exo-miRExplorer. Subsequently, we used exo-miRExplorer to perform further analysis. Comparative analysis of the exo-miRNAs between different sequencing datasets revealed significant correlation of exo-miRNAs between experiments in the same study. The plant-derived exo-miRNAs analysis provided robust evidence for non-diet source of exo-miRNAs. Virus-derived exo-miRNA analysis showed that pathogen RNAs could transfer to host cells and exist in deep-sequencing result at abundance level. In conclusion, exo-miRExplorer provides users with an integrative resource to facilitate detection and analysis of exo-miRNAs. exo-miRExplorer is available at the following URL: http://rna.sysu.edu.cn/exomiRDB/.

Emergence and evolution of inter-specific segregating retrocopies in cynomolgus monkey (Macaca fascicularis) and rhesus macaque (Macaca mulatta)

Retroposition is an RNA-mediated mechanism to generate gene duplication, and is believed to play an important role in genome evolution and phenotypic adaptation in various species including primates. Previous studies suggested an elevated rate of recent retroposition in the rhesus macaque genome. To better understand the impact of retroposition on macaque species which have undergone an adaptive radiation approximately 3–6 million years ago, we developed a bioinformatics pipeline to identify recently derived retrocopies in cynomolgus monkeys. As a result, we identified seven experimentally validated young retrocopies, all of which are polymorphic in cynomolgus monkeys. Unexpectedly, five of them are also present in rhesus monkeys and are still segregating. Molecular evolutionary analysis indicates that the observed inter-specific polymorphism is attribute to ancestral polymorphism. Further population genetics analysis provided strong evidence of balancing selection on at least one case (Crab-eating monkey retrocopy 6, or CER6) in both species. CER6 is in adjacent with an immunoglobulin related gene and may be involved in host-pathogen interaction, a well-known target of balancing selection. Altogether, our data support that retroposition is an important force to shape genome evolution and species adaptation.

Evolutionary direction of processed pseudogenes

While some pseudogenes have been reported to play important roles in gene regulation, little is known about the possible relationship between pseudogene functions and evolutionary process of pseudogenes, or about the forces responsible for the pseudogene evolution. In this study, we characterized human processed pseudogenes in terms of evolutionary dynamics. Our results show that pseudogenes tend to evolve toward: lower GC content, strong dinucleotide bias, reduced abundance of transcription factor binding motifs and short palindromes, and decreased ability to form nucleosomes. We explored possible evolutionary forces that shaped the evolution pattern of pseudogenes, and concluded that mutations in pseudogenes are likely determined, at least partially, by neighbor-dependent mutational bias and recombination-associated selection.

Two sets of RNAi components are required for heterochromatin formation in trans triggered by truncated transgenes

Across kingdoms, RNA interference (RNAi) has been shown to control gene expression at the transcriptional- or the post-transcriptional level. Here, we describe a mechanism which involves both aspects: truncated transgenes, which fail to produce intact mRNA, induce siRNA accumulation and silencing of homologous loci in trans in the ciliate Paramecium. We show that silencing is achieved by co-transcriptional silencing, associated with repressive histone marks at the endogenous gene. This is accompanied by secondary siRNA accumulation, strictly limited to the open reading frame of the remote locus. Our data shows that in this mechanism, heterochromatic marks depend on a variety of RNAi components. These include RDR3 and PTIWI14 as well as a second set of components, which are also involved in post-transcriptional silencing: RDR2, PTIWI13, DCR1 and CID2. Our data indicates differential processing of nascent un-spliced and long, spliced transcripts thus suggesting a hitherto-unrecognized functional interaction between post-transcriptional and co-transcriptional RNAi. Both sets of RNAi components are required for efficient trans-acting RNAi at the chromatin level and our data indicates similar mechanisms contributing to genome wide regulation of gene expression by epigenetic mechanisms.

Histone H3 Variant Regulates RNA Polymerase II Transcription Termination and Dual Strand Transcription of siRNA Loci in Trypanosoma brucei

Base J, β-D-glucosyl-hydroxymethyluracil, is a chromatin modification of thymine in the nuclear DNA of flagellated protozoa of the order Kinetoplastida. In Trypanosoma brucei, J is enriched, along with histone H3 variant (H3.V), at sites involved in RNA Polymerase (RNAP) II termination and telomeric sites involved in regulating variant surface glycoprotein gene (VSG) transcription by RNAP I. Reduction of J in T. brucei indicated a role of J in the regulation of RNAP II termination, where the loss of J at specific sites within polycistronic gene clusters led to read-through transcription and increased expression of downstream genes. We now demonstrate that the loss of H3.V leads to similar defects in RNAP II termination within gene clusters and increased expression of downstream genes. Gene derepression is intensified upon the subsequent loss of J in the H3.V knockout. mRNA-seq indicates gene derepression includes VSG genes within the silent RNAP I transcribed telomeric gene clusters, suggesting an important role for H3.V in telomeric gene repression and antigenic variation. Furthermore, the loss of H3.V at regions of overlapping transcription at the end of convergent gene clusters leads to increased nascent RNA and siRNA production. Our results suggest base J and H3.V can act independently as well as synergistically to regulate transcription termination and expression of coding and non-coding RNAs in T. brucei, depending on chromatin context (and transcribing polymerase). As such these studies provide the first direct evidence for histone H3.V negatively influencing transcription elongation to promote termination.

Nucleosome Organization around Pseudogenes in the Human Genome

Pseudogene, disabled copy of functional gene, plays a subtle role in gene expression and genome evolution. The first step in deciphering RNA-level regulation of pseudogenes is to understand their transcriptional activity. So far, there has been no report on possible roles of nucleosome organization in pseudogene transcription. In this paper, we investigated the effect of nucleosome positioning on pseudogene transcription. For transcribed pseudogenes, the experimental nucleosome occupancy shows a prominent depletion at the regions both upstream of pseudogene start positions and downstream of pseudogene end positions. Intriguingly, the same depletion is also observed for nontranscribed pseudogenes, which is unexpected since nucleosome depletion in those regions is thought to be unnecessary in light of the nontranscriptional property of those pseudogenes. The sequence-dependent prediction of nucleosome occupancy shows a consistent pattern with the experimental data-based analysis. Our results indicate that nucleosome positioning may play important roles in both the transcription initiation and termination of pseudogenes.

The genome of Leishmania panamensis: insights into genomics of the L. (Viannia) subgenus

Kinetoplastid parasites of the Leishmania genus cause several forms of leishmaniasis. Leishmania species pathogenic to human are separated into two subgenera, Leishmania (Leishmania) and L. (Viannia). Species from the Viannia subgenus cause predominantly cutaneous leishmaniasis in Central and South America, occasionally leading to more severe clinical presentations. Although the genomes of several species of Leishmania have been sequenced to date, only one belongs to this rather different subgenus. Here we explore the unique features of the Viannia subgenus by sequencing and analyzing the genome of L. (Viannia) panamensis. Against a background of conservation in gene content and synteny, we found key differences at the genomic level that may explain the occurrence of molecular processes involving nucleic acid manipulation and differential modification of surface glycoconjugates. These differences may in part explain some phenotypic characteristics of the Viannia parasites, including their increased adaptive capacity and enhanced metastatic ability.

GENCODE pseudogenes

Historically pseudogenes were believed to represent nonfunctional genomic fossils; however, there is emerging evidence that many of them could be biologically active. This possibility has ignited interest in pseudogene loci and made the need for their high-quality annotation more pressing as an accurate knowledge of all pseudogenes in the human reference genome sequence facilitates confident functional analysis. GENCODE have undertaken the first genome-wide pseudogene assignment for protein-coding genes combining both large-scale manual annotation and computational pseudogene prediction pipelines. Multiple computational predictions provide an unbiased set of hints for manual annotators to investigate, both during first-pass annotation and as part of QC to identify any potential missing pseudogene loci. Where a pseudogene is identified, the extent of its homology to the parent locus is fully investigated by a manual annotator; a pseudogene model is built and assigned to one of eight pseudogene biotypes depending on the mechanism of creation and on the presence of locus-specific transcriptional or proteomic data. The high-quality, information-rich set of pseudogenes created has been integrated with ENCODE functional genomics data, specifically expression level, transcription factor and RNA polymerase II binding, and chromatin marks. In this way we have been able to identify some pseudogenes that possess conventional characteristics of functionality as well as others with interesting patterns of partial activity, which might suggest that putatively inactive loci could be gaining a novel function, for example as long noncoding RNAs. The activity data associated with every pseudogene is stored in the psiDR resource.

Pseudogene-derived endogenous siRNAs and their function

Pseudogenes were once considered genomic fossils, but recent studies indicate that they may function as gene regulators through the generation of endogenous small interfering RNAs (esiRNAs), antisense RNAs, and decoys for microRNAs. In this review, we summarize pseudogene study methods, emphasizing relevant publicly available resources, and we describe a systematic pipeline to identify pseudogene-derived esiRNAs and their targets, which can lead to a deeper understanding of pseudogene function.

Making novel proteins from pseudogenes

MOTIVATION

Recently, we made synthetic proteins from non-coding DNA of Escherichia coli. Encouraged by this, we asked: can we artificially express pseudogenes into novel and functional proteins? What kind of structures would be generated? Would these proteins be stable? How would the organism respond to the artificial reactivation of pseudogenes?

RESULTS

To answer these questions, we studied 16 full-length protein equivalents of pseudogenes. The sequence-based predictions indicated interesting molecular and cellular functional roles for pseudogene-derived proteins. Most of the proteins were predicted to be involved in the amino acid biosynthesis, energy metabolism, purines and pyrimidine biosynthesis, central intermediary metabolism, transport and binding. Interestingly, many of the pseudogene-derived proteins were predicted to be enzymes. Furthermore, proteins showed strong evidence of stable tertiary structures. The prediction scores for structure, function and stability were found to be favorable in most of the cases.

IMPACT

To our best knowledge, this is the first such report that predicts the possibility of making functional and stable proteins from pseudogenes. In future, it would be interesting to experimentally synthesize and validate these predictions.

Pseudogenes as regulators of biological function

A pseudogene arises when a gene loses the ability to produce a protein, which can be due to mutation or inaccurate duplication. Previous dogma has dictated that because the pseudogene no longer produces a protein it becomes functionless and evolutionarily inert, being neither conserved nor removed. However, recent evidence has forced a re-evaluation of this view. Some pseudogenes, although not translated into protein, are at least transcribed into RNA. In some cases, these pseudogene transcripts are capable of influencing the activity of other genes that code for proteins, thereby altering expression and in turn affecting the phenotype of the organism. In the present chapter, we will define pseudogenes, describe the evidence that they are transcribed into non-coding RNAs and outline the mechanisms by which they are able to influence the machinery of the eukaryotic cell.

Characterization of human pseudogene-derived non-coding RNAs for functional potential

Thousands of pseudogenes exist in the human genome and many are transcribed, but their functional potential remains elusive and understudied. To explore these issues systematically, we first developed a computational pipeline to identify transcribed pseudogenes from RNA-Seq data. Applying the pipeline to datasets from 16 distinct normal human tissues identified ∼ 3,000 pseudogenes that could produce non-coding RNAs in a manner of low abundance but high tissue specificity under normal physiological conditions. Cross-tissue comparison revealed that the transcriptional profiles of pseudogenes and their parent genes showed mostly positive correlations, suggesting that pseudogene transcription could have a positive effect on the expression of their parent genes, perhaps by functioning as competing endogenous RNAs (ceRNAs), as previously suggested and demonstrated with the PTEN pseudogene, PTENP1. Our analysis of the ENCODE project data also found many transcriptionally active pseudogenes in the GM12878 and K562 cell lines; moreover, it showed that many human pseudogenes produced small RNAs (sRNAs) and some pseudogene-derived sRNAs, especially those from antisense strands, exhibited evidence of interfering with gene expression. Further integrated analysis of transcriptomics and epigenomics data, however, demonstrated that trimethylation of histone 3 at lysine 9 (H3K9me3), a posttranslational modification typically associated with gene repression and heterochromatin, was enriched at many transcribed pseudogenes in a transcription-level dependent manner in the two cell lines. The H3K9me3 enrichment was more prominent in pseudogenes that produced sRNAs at pseudogene loci and their adjacent regions, an observation further supported by the co-enrichment of SETDB1 (a H3K9 methyltransferase), suggesting that pseudogene sRNAs may have a role in regional chromatin repression. Taken together, our comprehensive and systematic characterization of pseudogene transcription uncovers a complex picture of how pseudogene ncRNAs could influence gene and pseudogene expression, at both epigenetic and post-transcriptional levels.

Characterization of the small RNA content of Trypanosoma cruzi extracellular vesicles

A growing body of evidence in mammalian cells indicates that secreted vesicles can be used to mediate intercellular communication processes by transferring various bioactive molecules, including mRNAs and microRNAs. Based on these findings, we decided to analyze whether Trypanosoma cruzi-derived extracellular vesicles contain RNA molecules and performed a deep sequencing and genome-wide analysis of a size-fractioned cDNA library (16-40nt) from extracellular vesicles secreted by noninfective epimastigote and infective metacyclic trypomastigote forms. Our data show that the small RNAs contained in these extracellular vesicles originate from multiple sources, including tRNAs. In addition, our results reveal that the variety and expression of small RNAs are different between parasite stages, suggesting diverse functions. Taken together, these observations call attention to the potential regulatory functions that these RNAs might play once transferred between parasites and/or to mammalian host cells.

A contig-based strategy for the genome-wide discovery of microRNAs without complete genome resources

MicroRNAs (miRNAs) are important regulators of many cellular processes and exist in a wide range of eukaryotes. High-throughput sequencing is a mainstream method of miRNA identification through which it is possible to obtain the complete small RNA profile of an organism. Currently, most approaches to miRNA identification rely on a reference genome for the prediction of hairpin structures. However, many species of economic and phylogenetic importance are non-model organisms without complete genome sequences, and this limits miRNA discovery. Here, to overcome this limitation, we have developed a contig-based miRNA identification strategy. We applied this method to a triploid species of edible banana (GCTCV-119, Musa spp. AAA group) and identified 180 pre-miRNAs and 314 mature miRNAs, which is three times more than those were predicted by the available dataset-based methods (represented by EST+GSS). Based on the recently published miRNA data set of Musa acuminate, the recall rate and precision of our strategy are estimated to be 70.6% and 92.2%, respectively, significantly better than those of EST+GSS-based strategy (10.2% and 50.0%, respectively). Our novel, efficient and cost-effective strategy facilitates the study of the functional and evolutionary role of miRNAs, as well as miRNA-based molecular breeding, in non-model species of economic or evolutionary interest.

Analysis of variable retroduplications in human populations suggests coupling of retrotransposition to cell division

In primates and other animals, reverse transcription of mRNA followed by genomic integration creates retroduplications. Expressed retroduplications are either “retrogenes” coding for functioning proteins, or expressed “processed pseudogenes,” which can function as noncoding RNAs. To date, little is known about the variation in retroduplications in terms of their presence or absence across individuals in the human population. We have developed new methodologies that allow us to identify “novel” retroduplications (i.e., those not present in the reference genome), to find their insertion points, and to genotype them. Using these methods, we catalogued and analyzed 174 retroduplication variants in almost one thousand humans, which were sequenced as part of Phase 1 of The 1000 Genomes Project Consortium. The accuracy of our data set was corroborated by (1) multiple lines of sequencing evidence for retroduplication (e.g., depth of coverage in exons vs. introns), (2) experimental validation, and (3) the fact that we can reconstruct a correct phylogenetic tree of human subpopulations based solely on retroduplications. We also show that parent genes of retroduplication variants tend to be expressed at the M-to-G1 transition in the cell cycle and that M-to-G1 expressed genes have more copies of fixed retroduplications than genes expressed at other times. These findings suggest that cell division is coupled to retrotransposition and, perhaps, is even a requirement for it.

Comparative transcriptome analysis of small noncoding RNAs in different stages of Trypanosoma brucei

Trypanosoma brucei, a pathogen of human and domestic animals, is an early evolved parasitic protozoan with a complex life cycle. Most genes of this parasite are post-transcriptionally regulated. However, the mechanisms and the molecules involved remain largely unknown. We have deep-sequenced the small RNAs of two life stages of this parasite--the bloodstream form and the procyclic form. Our results show that the small RNAs of T. brucei could derive from multiple sources, including NATs (natural antisense transcripts), tRNAs, and rRNAs. Most of these small RNAs in the two stages were found to share uniform characteristics. However, our results demonstrate that their variety and expression show significant differences between different stages, indicating possible functional differentiation. Dicer-knockdown evidence further proved that some of the small interfering RNAs (siRNAs) could regulate the expression of genes. Based on the genome-wide analysis of the small RNAs in the two stages of T. brucei, our results not only provide evidence to study their differentiation but also shed light on questions regarding the origins and evolution of small RNA-based mechanisms in early eukaryotes.

Genome-wide association identifies genetic variants associated with lentiform nucleus volume in N = 1345 young and elderly subjects

Deficits in lentiform nucleus volume and morphometry are implicated in a number of genetically influenced disorders, including Parkinson's disease, schizophrenia, and ADHD. Here we performed genome-wide searches to discover common genetic variants associated with differences in lentiform nucleus volume in human populations. We assessed structural MRI scans of the brain in two large genotyped samples: the Alzheimer's Disease Neuroimaging Initiative (ADNI; N = 706) and the Queensland Twin Imaging Study (QTIM; N = 639). Statistics of association from each cohort were combined meta-analytically using a fixed-effects model to boost power and to reduce the prevalence of false positive findings. We identified a number of associations in and around the flavin-containing monooxygenase (FMO) gene cluster. The most highly associated SNP, rs1795240, was located in the FMO3 gene; after meta-analysis, it showed genome-wide significant evidence of association with lentiform nucleus volume (P MA  = 4.79 × 10(-8)). This commonly-carried genetic variant accounted for 2.68 % and 0.84 % of the trait variability in the ADNI and QTIM samples, respectively, even though the QTIM sample was on average 50 years younger. Pathway enrichment analysis revealed significant contributions of this gene to the cytochrome P450 pathway, which is involved in metabolizing numerous therapeutic drugs for pain, seizures, mania, depression, anxiety, and psychosis. The genetic variants we identified provide replicated, genome-wide significant evidence for the FMO gene cluster's involvement in lentiform nucleus volume differences in human populations.

Sex bias in copy number variation of olfactory receptor gene family depends on ethnicity

Gender plays a pivotal role in the human genetic identity and is also manifested in many genetic disorders particularly mental retardation. In this study its effect on copy number variation (CNV), known to cause genetic disorders was explored. As the olfactory receptor (OR) repertoire comprises the largest human gene family, it was selected for this study, which was carried out within and between three populations, derived from 150 individuals from the 1000 Genome Project. Analysis of 3872 CNVs detected among 791 OR loci, in which 307 loci showed CNV, revealed the following novel findings: Sex bias in CNV was significantly more prevalent in uncommon than common CNV variants of OR pseudogenes, in which the male genome showed more CNVs; and in one-copy number loss compared to complete deletion of OR pseudogenes; both findings implying a more recent evolutionary role for gender. Sex bias in copy number gain was also detected. Another novel finding was that the observed sex bias was largely dependent on ethnicity and was in general absent in East Asians. Using a CNV public database for sick children (International Standard Cytogenomic Array Consortium) the application of these findings for improving clinical molecular diagnostics is discussed by showing an example of sex bias in CNV among kids with autism. Additional clinical relevance is discussed, as the most polymorphic CNV-enriched OR cluster in the human genome, located on chr 15q11.2, is found near the Prader–Willi syndrome/Angelman syndrome bi-directionally imprinted region associated with two well-known mental retardation syndromes. As olfaction represents the primitive cognition in most mammals, arguably in competition with the development of a larger brain, the extensive retention of OR pseudogenes in females of this study, might point to a parent-of-origin indirect regulatory role for OR pseudogenes in the embryonic development of human brain. Thus any perturbation in the temporal regulation of olfactory system could lead to developmental delay disorders including mental retardation.

Retrotransposition of gene transcripts leads to structural variation in mammalian genomes

BACKGROUND

Retroposed processed gene transcripts are an important source of material for new gene formation on evolutionary timescales. Most prior work on gene retrocopy discovery compared copies in reference genome assemblies to their source genes. Here, we explore gene retrocopy insertion polymorphisms (GRIPs) that are present in the germlines of individual humans, mice, and chimpanzees, and we identify novel gene retrocopy insertions in cancerous somatic tissues that are absent from patient-matched non-cancer genomes.

RESULTS

Through analysis of whole-genome sequence data, we found evidence for 48 GRIPs in the genomes of one or more humans sequenced as part of the 1,000 Genomes Project and The Cancer Genome Atlas, but which were not in the human reference assembly. Similarly, we found evidence for 755 GRIPs at distinct locations in one or more of 17 inbred mouse strains but which were not in the mouse reference assembly, and 19 GRIPs across a cohort of 10 chimpanzee genomes, which were not in the chimpanzee reference genome assembly. Many of these insertions are new members of existing gene families whose source genes are highly and widely expressed, and the majority have detectable hallmarks of processed gene retrocopy formation. We estimate the rate of novel gene retrocopy insertions in humans and chimps at roughly one new gene retrocopy insertion for every 6,000 individuals.

CONCLUSIONS

We find that gene retrocopy polymorphisms are a widespread phenomenon, present a multi-species analysis of these events, and provide a method for their ascertainment.

The structure and repertoire of small interfering RNAs in Leishmania (Viannia) braziliensis reveal diversification in the trypanosomatid RNAi pathway

Among trypanosomatid protozoa the mechanism of RNA interference (RNAi) has been investigated in Trypanosoma brucei and to a lesser extent in Leishmania braziliensis. Although these two parasitic organisms belong to the same family, they are evolutionarily distantly related raising questions about the conservation of the RNAi pathway. Here we carried out an in-depth analysis of small interfering RNAs (siRNAs) associated with L. braziliensis Argonaute1 (LbrAGO1). In contrast to T. brucei, Leishmania siRNAs are sensitive to 3' end oxidation, indicating the absence of blocking groups, and the Leishmania genome does not code for a HEN1 RNA 2'-O-methyltransferase, which modifies small RNA 3' ends. Consistent with this observation, ~20% of siRNA 3' ends carry non-templated uridines. Thus siRNA biogenesis, and most likely their metabolism, is different in these organisms. Similarly to T. brucei, putative mobile elements and repeats constitute the major Leishmania siRNA-producing loci and AGO1 ablation leads to accumulation of long transcripts derived from putative mobile elements. However, contrary to T. brucei, no siRNAs were detected from other genomic regions with the potential to form double-stranded RNA, namely sites of convergent transcription and inverted repeats. Thus, our results indicate that organism-specific diversification has occurred in the RNAi pathway during evolution of the trypanosomatid lineage.

Small interfering RNA-producing loci in the ancient parasitic eukaryote Trypanosoma brucei

Background

At the core of the RNA interference (RNAi) pathway in Trypanosoma brucei is a single Argonaute protein, TbAGO1, with an established role in controlling retroposon and repeat transcripts. Recent evidence from higher eukaryotes suggests that a variety of genomic sequences with the potential to produce double-stranded RNA are sources for small interfering RNAs (siRNAs).

Results

To test whether such endogenous siRNAs are present in T. brucei and to probe the individual role of the two Dicer-like enzymes, we affinity purified TbAGO1 from wild-type procyclic trypanosomes, as well as from cells deficient in the cytoplasmic (TbDCL1) or nuclear (TbDCL2) Dicer, and subjected the bound RNAs to Illumina high-throughput sequencing. In wild-type cells the majority of reads originated from two classes of retroposons. We also considerably expanded the repertoire of trypanosome siRNAs to encompass a family of 147-bp satellite-like repeats, many of the regions where RNA polymerase II transcription converges, large inverted repeats and two pseudogenes. Production of these newly described siRNAs is strictly dependent on the nuclear DCL2. Notably, our data indicate that putative centromeric regions, excluding the CIR147 repeats, are not a significant source for endogenous siRNAs.

Conclusions

Our data suggest that endogenous RNAi targets may be as evolutionarily old as the mechanism itself.

Specific inactivation of two immunomodulatory SIGLEC genes during human evolution

Sialic acid-recognizing Ig-like lectins (Siglecs) are signaling receptors that modulate immune responses, and are targeted for interactions by certain pathogens. We describe two primate Siglecs that were rendered nonfunctional by single genetic events during hominin evolution after our common ancestor with the chimpanzee. SIGLEC13 was deleted by an Alu-mediated recombination event, and a single base pair deletion disrupted the ORF of SIGLEC17. Siglec-13 is expressed on chimpanzee monocytes, innate immune cells that react to bacteria. The human SIGLEC17P pseudogene mRNA is still expressed at high levels in human natural killer cells, which bridge innate and adaptive immune responses. As both resulting pseudogenes are homozygous in all human populations, we resurrected the originally encoded proteins and examined their functions. Chimpanzee Siglec-13 and the resurrected human Siglec-17 recruit a signaling adapter and bind sialic acids. Expression of either Siglec in innate immune cells alters inflammatory cytokine secretion in response to Toll-like receptor-4 stimulation. Both Siglecs can also be engaged by two potentially lethal sialylated bacterial pathogens of newborns and infants, agents with a potential impact on reproductive fitness. Neanderthal and Denisovan genomes show human-like sequences at both loci, corroborating estimates that the initial pseudogenization events occurred in the common ancestral population of these hominins. Both loci also show limited polymorphic diversity, suggesting selection forces predating the origin of modern humans. Taken together, these data suggest that genetic elimination of Siglec-13 and/or Siglec-17 represents signatures of infectious and/or other inflammatory selective processes contributing to population restrictions during hominin origins.

Pseudogenes are not pseudo any more

Recent significant progress toward understanding the function of pseudogenes in protozoa (Trypanosoma brucei), metazoa (mouse) and plants, make it pertinent to provide a brief overview on what has been learned about this fascinating subject. We discuss the regulatory mechanisms of pseudogenes at the post-transcriptional level and advance new ideas toward understanding the evolution of these, sometimes called "garbage genes" or "junk DNA," seeking to stimulate the interest of scientists and additional research on the subject. We hope this point-of-view can be helpful to scientists working or seeking to work on these and related issues.

The short non-coding transcriptome of the protozoan parasite Trypanosoma cruzi

The pathway for RNA interference is widespread in metazoans and participates in numerous cellular tasks, from gene silencing to chromatin remodeling and protection against retrotransposition. The unicellular eukaryote Trypanosoma cruzi is missing the canonical RNAi pathway and is unable to induce RNAi-related processes. To further understand alternative RNA pathways operating in this organism, we have performed deep sequencing and genome-wide analyses of a size-fractioned cDNA library (16–61 nt) from the epimastigote life stage. Deep sequencing generated 582,243 short sequences of which 91% could be aligned with the genome sequence. About 95–98% of the aligned data (depending on the haplotype) corresponded to small RNAs derived from tRNAs, rRNAs, snRNAs and snoRNAs. The largest class consisted of tRNA-derived small RNAs which primarily originated from the 3′ end of tRNAs, followed by small RNAs derived from rRNA. The remaining sequences revealed the presence of 92 novel transcribed loci, of which 79 did not show homology to known RNA classes.

Chagas' disease is a major health problem in Latin America and is caused by the protozoan parasite Trypanosoma cruzi. T. cruzi lacks the pathway for RNA interference, which is widespread among eukaryotes, and is therefore unable to induce RNAi-related processes. In many organisms, small RNAs play an important role in regulating gene expression and other cellular processes. In order to understand if other small RNA pathways are operating in this organism, we performed high throughput sequencing and genome-wide analyses of the short transcriptome. We identified an abundance of small RNAs derived from non-coding RNA genes, including transfer RNAs, ribosomal RNAs as well as small nucleolar RNAs and small nuclear RNAs. Certain tRNA types were overrepresented as precursors for small RNAs. Further, we identified 79 novel small non-coding RNAs, not previously reported. We did not identify canonical small RNAs, like microRNAs and small interfering RNAs, and concluded that these do not exist in T. cruzi. This study has provided insights into the short transcriptome of a major human pathogen and provided starting points for further functional investigation of small RNAs and their biological roles.

RNA interference in protozoan parasites: achievements and challenges

Protozoan parasites that profoundly affect mankind represent an exceptionally diverse group of organisms, including Plasmodium, Toxoplasma, Entamoeba, Giardia, trypanosomes, and Leishmania. Despite the overwhelming impact of these parasites, there remain many aspects to be discovered about mechanisms of pathogenesis and how these organisms survive in the host. Combined with the ever-increasing availability of sequenced genomes, RNA interference (RNAi), discovered a mere 13 years ago, has enormously facilitated the analysis of gene function, especially in organisms that are not amenable to classical genetic approaches. Here we review the current status of RNAi in studies of parasitic protozoa, with special emphasis on its use as a postgenomic tool.