Quantitative profiling of the UGT transcriptome in human drug-metabolizing tissues.

Alternative splicing as a mean to control gene expression and diversify function is suspected to considerably influence drug response and clearance. We report the quantitative expression profiles of the human UGT genes including alternatively spliced variants not previously annotated established by deep RNA-sequencing in tissues of pharmacological importance. We reveal a comprehensive quantification of the alternative UGT transcriptome that differ across tissues and among individuals. Alternative transcripts that comprise novel in-frame sequences associated or not with truncations of the 5'- and/or 3'- termini, significantly contribute to the total expression levels of each UGT1 and UGT2 gene averaging 21% in normal tissues, with expression of UGT2 variants surpassing those of UGT1. Quantitative data expose preferential tissue expression patterns and remodeling in favor of alternative variants upon tumorigenesis. These complex alternative splicing programs have the strong potential to contribute to interindividual variability in drug metabolism in addition to diversify the UGT proteome.

Unravelling the transcriptomic landscape of the major phase II UDP-glucuronosyltransferase drug metabolizing pathway using targeted RNA sequencing.

A comprehensive view of the human UDP-glucuronosyltransferase (UGT) transcriptome is a prerequisite to the establishment of an individual's UGT metabolic glucuronidation signature. Here, we uncover the transcriptome landscape of the 10 human UGT gene loci in normal and tumoral metabolic tissues by targeted RNA next-generation sequencing. Alignment on the human hg19 reference genome identifies 234 novel exon-exon junctions. We recover all previously known UGT1 and UGT2 enzyme-coding transcripts and identify over 130 structurally and functionally diverse novel UGT variants. We further expose a revised genomic structure of UGT loci and provide a comprehensive repertoire of transcripts for each UGT gene. Data also uncover a remodelling of the UGT transcriptome occurring in a tissue- and tumor-specific manner. The complex alternative splicing program regulating UGT expression and protein functions is likely critical in determining detoxification capacity of an organ and stress-related responses, with significant impact on drug responses and diseases.

A novel UGT1 marker associated with better tolerance against irinotecan-induced severe neutropenia in metastatic colorectal cancer patients.

The risk of severe irinotecan-induced neutropenia has been shown to be related to the UGT1 variant UGT1A1*28, which increases exposure to the potent metabolite SN-38. Our goal was to identify a novel UGT1 marker(s) using 28 haplotype-tagged single nucleotide polymorphisms genotyped by mass spectrometry. By characterizing the UGT1 sequence from a cohort of 167 Canadian metastatic colorectal cancer (mCRC) patients and a validation cohort of 250 Italian mCRC patients, we found rs11563250G, located in the intergenic region downstream of UGT1, to be significantly associated with reduced risk of severe neutropenia (odds ratio (OR)=0.21; P=0.043 and OR=0.27; P=0.036, respectively, and OR=0.31 when combined; P=0.001), which remained significant upon correction for multiple testing in the combined cohort (P=0.041). For the two-marker haplotype rs11563250G and UGT1A1*1 (rs8175347 TA6), the OR was of 0.17 (P=0.0004). Genetic testing of this marker may identify patients who might benefit from increased irinotecan dosing.

Keys for microRNA expression profiling of single rat hypothalamic nuclei and multiplex sequencing strategies.

Integrative research has taken on the challenge of addressing questions in physiology by using novel knowledge and novel techniques. Recently, small and long non-coding RNAs have emerged as key regulators of gene expression, while next-generation sequencing technologies have revolutionized the characterization of genomes and gene expression. For a decade, it has been known that microRNAs (miRNAs) are RNAs of 18-24 bases that regulate gene expression in mammals. Here, we first describe the nature of miRNAs and the advantages of high-throughput sequencing technologies for establishing miRNA expression profiles. The hypothalamus harbours a dozen specialized areas or nuclei, the sampling of which is required to establish physiologically relevant miRNA expression profiles. MicroRNA expression profiling from single animals is also important for investigating potential genetic or epigenetic differences between individuals. Establishing a large number of miRNA expression profiles of individual hypothalamic nuclei of single rats at a cost compatible with laboratory finance can be achieved by using tagged cDNA libraries constructed from purified small RNAs and a multiplex sequencing strategy. We continue this report by surveying specificities of the different strategies that are used at present for constructing tagged cDNA libraries and provide a comparative analysis of miRNA expression profiles from hypothalamic arcuate nuclei of seven male Wistar rats.

MicroRNA expression profiling of hypothalamic arcuate and paraventricular nuclei from single rats using Illumina sequencing technology.

MicroRNAs (miRNAs) finely tune messenger RNA (mRNA) expression. As the brain is a highly heterogeneous tissue, physiologically relevant miRNA expression profiling greatly benefits from sampling brain regions or nuclei. MiRNA expression profiling from individual samples is also important for investigating potential differences between animals according to their physiological and pathophysiological status. We have punched the arcuate (ARC) and paraventricular (PVN) nuclei from the hypothalamus of seven male Wistar rats and used them to establish a novel method for the characterization of the miRNA expression profile of individual rat brain nuclei. The identity of the ARC and PVN samples was checked for proopiomelanocortin and arginine vasopressin mRNA expression, respectively. Individual cDNA libraries were constructed from purified RNAs between 16 and 26 bases, using barcoded adapters. Libraries were multiplexed and sequenced using Illumina technology to a read depth >10(5). The ARC and PVN profiles displayed similar expression from a set of more than 210 miRNA genes. Expression was high or moderate for about twenty miRNAs that may be used to define a common ARC/PVN prototype profile of male Wistar rats. These miRNAs included seven of the eight genes of the let-7 family, the two miR-7 genes, miR-9 gene and 5' copy of the three miR-30 loci. Our method shows that the ARC and PVN from a single rat are accessible for miRNA digital characterization. This method will allow miRNA transcriptome characterization for any rat brain substructure or nuclei that can be microdissected.

Searching for excystment-regulated genes in Sterkiella histriomuscorum (Ciliophora, Oxytrichidae): a mRNA differential display analysis of gene expression in excysting cells.

In the absence of food, the oxytrichid Sterkiella histriomuscorum transforms like many ciliates into resting cysts. When transferred back into feeding medium, the cyst re-transforms into a vegetative cell. The entry into and exit from the dormant cyst stage are complex developmental processes still poorly investigated at the molecular level. Assuming that these changes in state could involve changes in gene expression, we have used the technique of mRNA differential display to detect differentially expressed genes in cysts and two different stages of excysting cell. Variation in the temporal expression pattern of transcripts could be detected and, in using an inverse-PCR strategy on circularized macronuclear DNA, we have sequenced the macronuclear genes of three of the isolated cDNAs. which correspond to 1) a nucleotide-binding domain-encoding gene, 2) a DHHC-domain-carrying gene, and 3) a phosphatase type 2C-encoding gene. For the first two genes, Northern blot analyses supported an excystment-associated regulated gene expression. We discuss their possible role during excystment and we show that the combination of differential display and inverse PCR constitutes a powerful approach to isolate excystment-regulated genes in hypotrichs.

Messenger RNA in dormant cells of Sterkiella histriomuscorum (Oxytrichiade): indentification of putative regulatory gene transcripts.

In the absence of food, the oxytrichid Sterkiella histriomuscorum, like many ciliates, enters into dormancy and transforms into a round and walled encysted cell. When transferred back into a feeding medium, the cyst re-transforms into a vegetative cell in a few hours. This encystment-excystment pathway, which is common to many free-living and parasitic protists, is still poorly understood at the molecular level. In order to identify potential dormant transcripts in the cysts of Sterkiella, we have constructed cDNA libraries from mature cysts. Transcripts have been isolated confirming the presence of a mRNA pool in the dormant cells. The sequence analysis of two cDNA indicates open reading frames which show significant similarities to known proteins involved in mechanisms of regulation: 1) nifR3, an element of the nitrogen regulatory system in bacteria and 2) CROC-1, a newly identified human transcription factor. The two corresponding macronuclear genes represent the first putative regulatory genes isolated in ciliates. From a differential screening of the cDNA library against vegetative cDNA, one cyst-specific (and very abundant) transcript has been isolated but the product has not yet been identified. The possible involvment of these new ciliate genes in the excystment process is discussed.

Protein coding gene trees in ciliates: comparison with rRNA-based phylogenies.

We have reexamined the phylogeny of the ciliates using alpha-tubulin and phosphoglycerate kinase gene sequences. For alpha-tubulin, we have compared the amino acid and nucleotide sequences of 20 species representing seven of the nine classes of the phylum (Karyorelictea, Heterotrichea, Hypotrichea, Oligohymenophorea, Colpodea, Nassophorea, and Litostomatea). The phylogenetic tree resembles a bush from which three monophyletic lineages can be distinguished which correspond to the three classes Hypotrichea, Oligohymenophorea, and Litostomatea. For phosphoglycerate kinase, we have compared the amino acid sequences from 7 species representing three classes (Heterotrichea, Hypotrichea, and Oligohymenophorea). The branching pattern is resolved in three deeply separated branches with an early emergence of the heterotrich. Our comparative analysis shows that if alpha-tubulin phylogeny is not informative at the interclass level, the preliminary data from the phosphoglycerate kinase molecule appear more promising. Nevertheless, at low taxonomic level and at the class level, the resolved phylogenetic relationships inferred from both protein and rRNA sequence data are congruent.

Genetic code deviations in the ciliates: evidence for multiple and independent events.

In several species of ciliates, the universal stop codons UAA and UAG are translated into glutamine, while in the euplotids, the glutamine codon usage is normal, but UGA appears to be translated as cysteine. Because the emerging position of this monophyletic group in the eukaryotic lineage is relatively late, this deviant genetic code represents a derived state of the universal code. The question is therefore raised as to how these changes arose within the evolutionary pathways of the phylum. Here, we have investigated the presence of stop codons in alpha tubulin and/or phosphoglycerate kinase gene coding sequences from diverse species of ciliates scattered over the phylogenetic tree constructed from 28S rRNA sequences. In our data set, when deviations occur they correspond to in frame UAA and UAG coding for glutamine. By combining these new data with those previously reported, we show that (i) utilization of UAA and UAG codons occurs to different extents between, but also within, the different classes of ciliates and (ii) the resulting phylogenetic pattern of deviations from the universal code cannot be accounted for by a scenario involving a single transition to the unusual code. Thus, contrary to expectations, deviations from the universal genetic code have arisen independently several times within the phylum.

Isolation and characterization of libraries of monoclonal antibodies directed against various forms of tubulin in Paramecium.

Ciliates are very good models for studying post-translationally generated tubulin heterogeneity because they exhibit highly differentiated microtubular networks in combination with reduced genetic diversity. We have approached the analysis of tubulin heterogeneity in Paramecium through extensive isolation and characterization of monoclonal antibodies using various antigens and several immunization protocols. Eight monoclonal antibodies and 10 hybridoma supernatants were characterized by: i) immunoblotting on ciliate and pig brain tubulins as well as on peptide maps of Paramecium axonemal tubulin; ii) immunoblotting on ciliate tubulin fusion peptides generated in E coli, a procedure which allows in principle to discriminate antibodies that are directed against tubulin sequence (reactive on fusion peptides) from those directed against a post-translational epitope (non-reactive); and iii) immunofluorescence on Paramecium, 3T3 and PtK2 cells. Twelve antibodies labeled all microtubules in Paramecium cells and were found to be directed against tubulin primary sequences (nine of them being located in the alpha N-terminal domain, one in the beta C-terminal one, and two in alpha and beta central stretches). The remaining ones decorated only a specific subset of microtubules within the cell and were presumably directed against post-translational modifications. Among these, three antibodies are directed against an N-terminal acetylated epitope of alpha-tubulin whereas the epitopes of three other ones (TAP 952 degrees, AXO 58 and AXO 49 degrees) apparently correspond to still unidentified post-translational modifications, located in the C-terminal domain of both alpha- and beta-tubulins. The AXO 49 degrees specificity is similar to that of a previously described polyclonal serum raised against Paramecium axonemal tubulin [2]. The results are discussed in terms of identification and accessibility of the epitopes and immunogenicity of ciliate tubulin with reference to mammalian and ciliate tubulin sequences.

A broad molecular phylogeny of ciliates: identification of major evolutionary trends and radiations within the phylum.

The cellular architecture of ciliates is one of the most complex known within eukaryotes. Detailed systematic schemes have thus been constructed through extensive comparative morphological and ultrastructural analysis of the ciliature and of its internal cytoskeletal derivatives (the infraciliature), as well as of the architecture of the oral apparatus. In recent years, a consensus was reached in which the phylum was divided in eight classes as defined by Lynn and Corliss [Lynn, D. H. & Corliss, J. O. (1991) in Microscopic Anatomy of Invertebrates: Protozoa (Wiley-Liss, New York), Vol. 1, pp. 333-467]. By comparing partial sequences of the large subunit rRNA molecule, and by using both distance-matrix and maximum-parsimony-tree construction methods (checked by boot-strapping), we examine the phylogenetic relationships of 22 species belonging to seven of these eight classes. At low taxonomic levels, the traditional grouping of the species is generally confirmed. At higher taxonomic levels, the branching pattern of these seven classes is resolved in several deeply separated major branches. Surprisingly, the first emerging one contains the heterotrichs and is strongly associated with a karyorelictid but deeply separated from hypotrichs. The litostomes, the oligohymenophorans, and the hypotrichs separate later in a bush-like topology hindering the resolution of their order of diversification. These results show a much more ancient origin of heterotrichs than was classically assumed, indicating that asymmetric, abundantly ciliated oral apparatuses do not correspond to "highly evolved" traits as previously thought. They also suggest the occurrence of a major radiative explosion in the evolutionary history of the ciliates, yielding five of the eight classes of the phylum. These classes appear to differ essentially according to the cytoskeletal architecture used to shape and sustain the cellular cortex (a process of essential adaptative and morphogenetic importance in ciliates).

[Eukaryogenesis: a model derivated from ribosomal RNA molecular phylogenise]. / L'Eucaryogenèse: un modèle dérivé des phylogénies moléculaires basées sur les ARN ribosomiques.

We have undertaken the construction of a broad molecular phylogeny of protists through the comparison of 28S rRNA molecules. The sequences from several major protistan phyla were aligned and combined with a broad database of metazoans, metaphytes, fungi and bacteria and we have derived dendrograms from both distance matrix and parsimony methods. In agreement with classical systematics, a number of monophyletic groups separated by large evolutionary distances were observed (those of the ciliates, the chlorophytes, etc.). From this analysis, several inferences on the eukaryogenesis can be made among which the ancient origin of the cytoskeleton, the late occurrence of the chloroplastic endosymbiosis and the simultaneous emergence of the triploblastic and diploblastic metazoan patterns.