Evaluating and Correcting Inherent Bias of microRNA Expression in Illumina Sequencing Analysis.

microRNA (miRNA) expression profiles based on the highly powerful Illumina sequencing technology rely on the construction of cDNA libraries in which adaptor ligation is known to deeply favor some miRNAs over others. This introduces erroneous measurements of the miRNA abundances and relative miRNA quantities in biological samples. Here, by using the commercial miRXplore Universal Reference that contains an equimolar mixture of 963 animal miRNAs and TruSeq or bulged adaptors, we describe a method for correcting ligation biases in expression profiles obtained with standard protocols of cDNA library construction and provide data for quantifying the true miRNA abundances in biological samples. Ligation biases were evaluated at three ratios of miRNA to 3'-adaptor and four numbers of polymerase chain reaction amplification cycles by calculating efficiency captures/correcting factors for each miRNA. We show that ligation biases lead to over- or under-expression covering a 105 amplitude range. We also show that, at each miRNA:3'-adaptor ratio, coefficients of variation (CVs) of efficiency captures calculated over the four number of amplification cycles using sliding windows of 10 values ranged from 0.1 for the miRNAs of high expression to 0.6 for the miRNAs of low expression. Efficiency captures of miRNAs of high and low expression in profiles are therefore differently impacted by the number of amplification cycles. Importantly, we observed that at a given number of amplification cycles, CVs of efficiency captures calculated over the three miRNA:3'-adaptor ratios displayed a steady value of 0.3 +/- 0.05 STD for miRNAs of high and low expression. This allows, at a given number of amplification cycles, accurate comparison of miRNA expression between biological samples over a substantial expression range. Finally we provide tables of correcting factors that allow to measure the abundances of 963 miRNAs in biological samples from TruSeq-based expression profiles and, an example of their use by characterizing miRNAs of the let-7, miR-26, miR-29, and miR-30 families as the more abundant miRNAs of the rat adult cerebellum.

miRNA Long-Term Response to Early Metabolic Environmental Challenge in Hypothalamic Arcuate Nucleus.

Epidemiological reports and studies using rodent models indicate that early exposure to nutrient and/or hormonal challenges can reprogram metabolism at adulthood. Hypothalamic arcuate nucleus (ARC) integrates peripheral and central signals to adequately regulate energy homeostasis. microRNAs (miRNAs) participate in the control of gene expression of large regulatory networks including many signaling pathways involved in epigenetics regulations. Here, we have characterized and compared the miRNA population of ARC of adult male rats continuously exposed to a balanced metabolic environment to the one of adult male rats exposed to an unbalanced high-fat/high-carbohydrate/moderate-protein metabolic environment during the perinatal period and/or at adulthood that consequently displayed hyperinsulinemia and/or hyperleptinemia. We identified more than 400 miRNA species in ARC of adult male rats. By comparing the miRNA content of six biological replicates in each of the four perinatal/adult environments/rat groups, we identified the 10 miRNAs specified by clusters miR-96/182/183, miR-141/200c, and miR-200a/200b/429 as miRNAs of systematic and uncommonly high variation of expression. This uncommon variation of expression may underlie high individual differences in aging disease susceptibilities. By comparing the miRNA content of the adult ARC between the rat groups, we showed that the miRNA population was not affected by the unbalanced adult environment while, in contrast, the expression of 11 miRNAs was repeatedly impacted by the perinatal unbalanced environment. Our data revealed a miRNA response of adult ARC to early metabolic environmental challenge.

Holistic and Affordable Analyses of MicroRNA Expression Profiles Using Tagged cDNA Libraries and a Multiplex Sequencing Strategy.

Small and long noncoding RNAs (ncRNAs) are key regulators of gene expression. Variations in ncRNA expression patterns can consequently affect the control of many cellular processes. Not just since 2006, when Andrew Z Fire and Craig C Mello were jointly awarded The Nobel Prize in Physiology or Medicine for their discovery of RNA interference, great efforts were undertaken to unleash the biomedical applicability of small noncoding RNAs, in particular microRNAs. With the technological evolution of massive parallel sequencing technologies over the last years, which now are available for an increasing number of scientists, there is a demand for comprehensible and efficient workflows reliable even for unique and valuable clinical specimens. Here we describe a highly reproducible low-cost protocol for analyses of miRNA expression patterns using tagged cDNA libraries and a multiplex sequencing strategy following an Illumina-like protocol. This protocol easily allows the identification of expression differences from samples of tissues of 1-2 mm3 and fluids of 50-200 µL. We further provide entry points into useful computational biology applications, whose target groups explicitly involve non-bioinformaticians.

Substantial and robust changes in microRNA transcriptome support postnatal development of the hypothalamus in rat.

MicroRNAs (miRNAs) modulate gene expression in male germ cells and somatic tissues of mammals on a genome-wide scale. Hundreds of miRNAs are encoded by mammalian genomes, a large fraction of which is expressed in brain. Here we have investigated the complexity and dynamics of miRNA transcriptomes that associate with neuronal network maturation of hypothalamic arcuate nucleus and median eminence (ARC/ME) in rat by analysing more than 300 miRNAs from 3-7 biological replicates at 5 postnatal time-points. The network connecting ARC/ME to other hypothalamic and extra-hypothalamic regions maturates in an environment-dependent manner. We therefore analyzed miRNA transcriptomes of progeny of dams fed either a balanced or unbalanced diet during gestation and lactation. More than 30% of the miRNAs displayed significative changes of expression between stages P8 and P14, and P21 and P28; half of the changes were greater than 3-fold. Among those miRNAs were well-known and dozens of still poorly documented miRNAs. Progeny of dams fed an unbanced diet displayed a severe growth retardation phenotype, lower levels of plasma leptin but almost identical miRNA transcriptomes. Together these data demonstrate that two substantial and robust changes in miRNA transcriptome of ARC/ME occur at a period crucial for neuronal network functional organization.

Multiple uses of Lys63-polyubiquitination in the ciliate Sterkiella histriomuscorum.

Poly-ubiquitination through Lys63, a post-translational modification system, is emerging as a mechanism implicated in a variety of non-proteolytic functions in eukaryotic cells. This modification is mainly carried out via a unique complex formed between a canonical ubiquitin-conjugating enzyme (E2), Ubc13, and an E2-variant, Uev. To date, a wealth of functional information has been obtained mostly from animals and yeast. Here, we show that the ciliate Sterkiella histriomuscorum harbors a Ubc13 gene that encodes a putative protein exhibiting 60-64% identity with Ubc13 from organisms of diverse eukaryotic phyla. A yeast double-hybrid assay and 3D-structure modelling gave evidence that ShUbc13 is able to physically interact with ShUev; the existence of such a complex in a ciliate attests for its ancient origin. Western blotting reveals that Ubc13 is present in different stages of the vegetative cycle of S. histriomuscorum. Immunolocalization of Ubc13 and K63 polyubiquitinated proteins shows three main locations, namely the nucleus, cilia and plasma membrane. Possible functions of ShUbc13 and K63 polyubiquitination in this single-celled organism are considered with reference to current knowledge. In particular their likely involvement in DNA replication and DNA damage response is especially discussed.

Differentially expressed genes during the encystment-excystment cycle of the ciliate Sterkiella histriomuscorum.

Under unfavourable environmental conditions, many ciliates transform into resting cysts through a developmental process called encystment. Excystment is the reverse transformation of the resting cyst into a vegetative cell when favourable conditions are restored. In the oxytrichid Sterkiella histriomuscorum, the encystment - excystment (E-E) cycle involves extensive morphological changes since the whole cytoskeleton is disassembled during encystment. Assuming that these changes in cellular organization may be significantly reflected in the gene expression pattern, we used a "DNA macroarray" strategy to measure the transcript levels of 37 selected genes present at four distinct cellular stages (starved, encysted, excysting and vegetative cells). Differential expression was observed for 16 genes; four transcripts appeared to be markedly accumulated in a stage-specific manner. For most of the differentially expressed genes, the mRNA level was increased in cysts and excysting cells. When these mRNA are transcribed and when they are used, are still open questions. We showed that the copy number of the differentially expressed genes is the same in the macronuclei of cysts and vegetative cells ruling out a modulation of gene expression through a variation in the gene copy number upon encystment.

Gene structure of the ciliate Sterkiella histriomuscorum based on a combined analysis of DNA and cDNA sequences from 21 macronuclear chromosomes.

Macronuclear deoxyribonucleic acid (DNA) in hypotrichous ciliates consists of a set of linear molecules ranging in size from 0.5 to several tens of kilobases and typically carrying a single gene. Each minichromosome is present at a ploidy of >or=1,000 per macronucleus. These molecules are known as gene-sized molecules. Multigene molecules are also present, but are still poorly described. In analyzing the encystment-excystment cycle of Sterkiella histriomuscorum, we have characterized a set of 21 macronuclear molecules both at the DNA and complementary DNA (cDNA) levels. On a total of 23 validated coding sequences, we mapped the 5' and 3' untranslated regions for a subset of 10 and 18 transcripts, respectively. A combination of DNA and cDNA data allows us to precisely determine several structural features of macronuclear chromosomes, such as the organization of multigene molecules, an intron content higher than expected, and a conserved sequence surrounding the initiation transcription site. It also reveals one coding sequence containing a transcribed 10-bp element that displays the characteristic features of internal eliminated sequences (IES). Its presence in a fraction of the minichromosomes carrying this gene raises the possibility of an incomplete IES excision process during the development of the S. histriomuscorum macronucleus.

Proposed function of the accumulation of plasma membrane-type Ca2+-ATPase mRNA in resting cysts of the ciliate Sterkiella histriomuscorum.

From an mRNA differential-display analysis of the encystment-excystment cycle of the ciliate Sterkiella histriomuscorum, we have isolated an expressed sequence tag encoding a plasma membrane-type Ca2+-ATPase (PMCA). PMCAs are located either in the plasma membranes or in the membranes of intracellular organelles, and their function is to pump calcium either out of the cell or into the intracellular calcium stores, respectively. The S. histriomuscorum macronuclear PMCA gene (ShPMCA) and its corresponding cDNA were cloned; it is the first member of the Ca2+-ATPase family identified in Sterkiella. The predicted protein of 1,065 amino acids exhibits 37% identity with PMCAs of diverse organisms. A phylogenetic analysis showed its relatedness to homologs of two alveolates: the ciliate Paramecium tetraurelia and the apicomplexan Toxoplasma gondii. Overexpression of the protein ShPMCA failed to rescue the wild-type phenotype of three Ca2+-ATPase-defective mutant strains of Saccharomyces cerevisiae; this failure contrasts with the reported ability of the PMCAs of parasites to complement defects in yeast. ShPMCA mRNA is markedly accumulated during encystment and in resting cysts, suggesting a function during excystment. To address the possibility of a signaling role for calcium at excystment, the capacity of calcium to induce excystment was examined.

Myb genes in ciliates: a common origin with the myb protooncogene?

In animals, the protooncogene myb family is characterized by a DNA-binding domain (so-called MYB domain), which consists of 3 imperfect tandem repeats of a helix-turn-helix motif. Homologous genes have been characterized in plants and also in Dictyostelium discoideum. However, in plants, the myb family is more diverse and displays 2 types of MYB domains: the animal-like 3 repeats (MYB-3R) and the 2 repeats (MYB-2R) domains. The question is therefore raised as to the putative existence of genes with MYB-3R and/or MYB-2R domains in their last common unicellular ancestor. Here, we present evidence that in ciliates like in plants, both types of domain exist. A gene having a MYB-3R domain has been identified in the oxytrichid Sterkiella histriomuscorum and a gene having a MYB-2R domain has been identified in the euplotid Euplotes aediculatus. Both genes are expressed during the vegetative growth of the cells. A conserved intron exists in the gene of Sterkiella and phylogenetical analyses show that the 2 ciliate genes belong to the myb protooncogene family as deeply split lineages. This is the first report of a myb homolog in a ciliated protist, thus, confirming its origin in strict unicellular eukaryotes.

Plateins: a novel family of signal peptide-containing articulins in euplotid ciliates.

In euplotid ciliates, the cortex is reinforced by alveolar plates--proteinaceous scales located within the membranous alveolar sacs, forming a monolayer just below the plasma membrane. This system appears to play a cytoskeletal role analogous to that provided by the fibrous epiplasm found beneath the cortical alveoli in other ciliates. In Euplotes aediculatus, the major alveolar plate proteins (termed alpha-, beta-, and gamma-plateins) have been identified. Using anti-platein antibodies, an expression library of Euplotes genes was screened, and a platein gene identified, cloned, and completely sequenced. Comparison of its derived amino acid sequence with microsequences obtained directly from purified plateins identified this gene as encoding one of the closely related beta- or gamma-plateins. The derived protein, of 644 amino acids (74.9 kDa), is very acidic (pI = 4.88). Microsequences from authentic alpha-platein were then used to design oligonucleotide primers, which yielded, via a PCR-based approach, the sequences of two alpha-platein genes from E. aediculatus. Even more acidic proteins, the derived alpha1- and alpha2-plateins contain 536 and 501 residues, respectively. Analyses of their amino acid sequences revealed the plateins to be members of the articulin superfamily of cytoskeletal proteins, first described in Euglena and now identified in the ciliate Pseudomicrothorax and in Plasmodium. The hallmark articulin repetitive motifs (based on degenerate valine- and proline-rich 12-mers) are present in all three plateins. In beta/gamma-platein this primary motif domain (27 repeats) is central in the molecule, whereas the primary repeats in the alpha-plateins lie near their C-termini. A cluster of proline-rich pentameric secondary repeats is found in the C-terminus of beta/gamma-platein, but near the N-terminus of alpha-plateins. All three plateins contain canonical N-terminal signal sequences, unique among known cytoskeletal proteins. The presence of start-transfer sequences correlates well with the final intra-alveolar location of these proteins. This feature, and significant differences from known articulins in amino acid usage and arrangement within the repeat domains, lead us to propose that the plateins comprise a new family of articulin-related proteins. Efforts to follow microscopically the assembly of plateins into new alveolar plates during pre-fission morphogenesis are underway.

Cytoskeletal proteins with N-terminal signal peptides: plateins in the ciliate Euplotes define a new family of articulins.

Protistan cells employ a wide variety of strategies to reinforce and give pattern to their outermost cortical layers. Whereas some use common cytoskeletal elements such as microtubules, others are based on novel cytoskeletal proteins that are as-yet-unknown in higher eukaryotes. The hypotrich ciliate Euplotes possesses a continuous monolayer of scales or plates, located within flattened membranous sacs ('alveoli') just below the plasma membrane, and this provides rigidity and form to the cell. Using immunological techniques, the major proteins comprising these 'alveolar plates' have been identified and termed alpha-, beta-, and gamma-plateins. The present report describes work leading to the molecular characterization of three plateins, alpha 1 and alpha 2 (predicted M(r)s of 61 and 56 kDa) and a beta/gamma form (M(r)=73 kDa). All three proteins have features that are hallmarks of articulins, a class of cytoskeletal proteins that has been identified in the cortex of a wide variety of protistan cells, including certain flagellates, ciliates, dinoflagellates and PLASMODIUM: Chief among these common features are a prominent primary domain of tandem 12-amino acid repeats, rich in valine and proline, and a secondary domain of fewer, shorter repeating units. However, variations in amino acid use within both primary and secondary repetitive domains, and a much more acidic character (predicted pIs of 4.7-4.9), indicate that the plateins represent the first proteins in a new subclass or family of articulins. This conclusion is supported by another novel feature of the plateins, the presence of a canonical hydrophobic signal peptide at the N-terminus of each derived platein sequence. This correlates well with the final cellular location of the plateins, which are assembled into plates within the membrane-limited alveolar sacs. To our knowledge, this is the first report in any eukaryote of cytoskeletal proteins with such start-transfer sequences. Confocal immunofluorescence microscopy, using antibodies to the plateins as probes, reveals that new alveolar plates (enlarging in cortical zones undergoing morphogenesis) label more faintly than mature parental plates. During plate assembly (or polymerization), the plateins thus appear to exist in a more soluble form.

A homologue of CROC-1 in a ciliated protist (Sterkiella histriomuscorum) testifies to the ancient origin of the ubiquitin-conjugating enzyme variant family.

Resting cysts of Sterkiella histriomuscorum (Ciliophora, Oxytrichidae) have been shown to contain messenger RNA, one of which codes for a protein significantly similar to CROC-1. CROC-1 is a human regulatory protein capable of transactivating the promoter of c-fos and belongs to a newly characterized family of ubiquitin-conjugating enzyme (E2) variants (UEV). We have determined the corresponding macronuclear gene sequence, which is the first protistan UEV sequence available. The phylogenetic analysis indicates the deep separation and solid clustering of all the UEV sequences within the E2 tree showing the ancient origin of these regulatory genes and their high structural conservation during evolution. Furthermore, overexpression of the ciliate UEV is able to rescue the Saccharomyces cerevisiae mms2 null mutant from killing by DNA damaging agents, implying that the UEV family proteins are functionally conserved. In S. histriomuscorum, expression of UEV is correlated with the growth of the cells as transcripts are present in excysting and vegetative cells but are rapidly down-regulated during starvation. These data support the high conservation of the UEV family in eukaryotes, and a regulatory role of the gene is discussed in relation to known functions of UEVs. This analysis may promote the search for homologues of other regulatory genes (metazoan regulators of differentiation) in ciliates.