Genome-wide Bioinformatics Approaches Research Articles

Regulation potential of transcribed simple repeated sequences in developing neurons

Simple repeated sequences (SRSs), defined as tandem iterations of microsatellite- to satellite-sized DNA units, occupy a substantial part of the human genome. Some of these elements are known to be transcribed in the context of repeat expansion disorders. Mounting evidence suggests that the transcription of SRSs may also contribute to normal cellular functions. Here, we used genome-wide bioinformatics approaches to systematically examine SRS transcriptional activity in cells undergoing neuronal differentiation. We identified thousands of long noncoding RNAs containing >200-nucleotide-long SRSs (SRS-lncRNAs), with hundreds of these transcripts significantly upregulated in the neural lineage. We show that SRS-lncRNAs often originate from telomere-proximal regions and that they have a strong potential to form multivalent contacts with a wide range of RNA-binding proteins. Our analyses also uncovered a cluster of neurally upregulated SRS-lncRNAs encoded in a centromere-proximal part of chromosome 9, which underwent an evolutionarily recent segmental duplication. Using a newly established in vitro system for rapid neuronal differentiation of induced pluripotent stem cells, we demonstrate that at least some of the bioinformatically predicted SRS-lncRNAs, including those encoded in the segmentally duplicated part of chromosome 9, indeed increase their expression in developing neurons to readily detectable levels. These and other lines of evidence suggest that many SRSs may be expressed in a cell type and developmental stage-specific manner, providing a valuable resource for further studies focused on the functional consequences of SRS-lncRNAs in the normal development of the human brain, as well as in the context of neurodevelopmental disorders.

The C. elegans 3' UTRome v2 resource for studying mRNA cleavage and polyadenylation, 3'-UTR biology, and miRNA targeting.

3′ Untranslated regions (3′ UTRs) of mRNAs emerged as central regulators of cellular function because they contain important but poorly characterized cis-regulatory elements targeted by a multitude of regulatory factors. The model nematode Caenorhabditis elegans is ideal to study these interactions because it possesses a well-defined 3′ UTRome. To improve its annotation, we have used a genome-wide bioinformatics approach to download raw transcriptome data for 1088 transcriptome data sets corresponding to the entire collection of C. elegans trancriptomes from 2015 to 2018 from the Sequence Read Archive at the NCBI. We then extracted and mapped high-quality 3′-UTR data at ultradeep coverage. Here, we describe and release to the community the updated version of the worm 3′ UTRome, which we named 3′ UTRome v2. This resource contains high-quality 3′-UTR data mapped at single-base ultraresolution for 23,084 3′-UTR isoform variants corresponding to 14,788 protein-coding genes and is updated to the latest release of WormBase. We used this data set to study and probe principles of mRNA cleavage and polyadenylation in C. elegans. The worm 3′ UTRome v2 represents the most comprehensive and high-resolution 3′-UTR data set available in C. elegans and provides a novel resource to investigate the mRNA cleavage and polyadenylation reaction, 3′-UTR biology, and miRNA targeting in a living organism.

Genome-wide characterization of microsatellites and marker development in the carcinogenic liver fluke Clonorchis sinensis

Clonorchis sinensis is an important carcinogenic human liver fluke endemic in East and Southeast Asia. There are several conventional molecular markers that have been used for identification and genetic diversity; however, no information about microsatellites of this liver fluke is published so far. We here report microsatellite characterization and marker development for a genetic diversity study in C. sinensis, using a genome-wide bioinformatics approach. Based on our search criteria, a total of 256,990 microsatellites (≥12 base pairs) were identified from a genome database of C. sinensis, with hexanucleotide motif being the most abundant (51%) followed by pentanucleotide (18.3%) and trinucleotide (12.7%). The tetranucleotide, dinucleotide, and mononucleotide motifs accounted for 9.75, 7.63, and 0.14%, respectively. The total length of all microsatellites accounts for 0. 72% of 547 Mb of the whole genome size, and the frequency of microsatellites was found to be one microsatellite in every 2.13 kb of DNA. For the di-, tri-, and tetranucleotide, the repeat numbers redundant are six (28%), four (45%), and three (76%), respectively. The ATC repeat is the most abundant microsatellites followed by AT, AAT, and AC, respectively. Within 40 microsatellite loci developed, 24 microsatellite markers showed potential to differentiate between C. sinensis and Opisthorchis viverrini. Seven out of 24 loci showed to be heterozygous with observed heterozygosity that ranged from 0.467 to 1. Four primer sets could amplify both C. sinensis and O. viverrini DNA with different sizes. This study provides basic information of C. sinensis microsatellites, and the genome-wide markers developed may be a useful tool for the genetic study of C. sinensis.

In silico analysis and developmental expression of ubiquitin-conjugating enzymes in Schistosoma mansoni.

Ubiquitin-conjugating enzymes (Ub-E2) perform the second step of ubiquitination and, consequently, are essential for regulating proteolysis and for modulating protein function, interactions and trafficking. Previously, our group demonstrated the crucial role of ubiquitination and the Ub-proteasome pathway during the Schistosoma mansoni life cycle. In the present investigation, we used a homology-based genome-wide bioinformatics approach to identify and molecularly characterise the Ub-E2 enzymes in S. mansoni. The putative functions were further investigated through molecular phylogenetic and expression profile analyses using cercariae, adult worms, eggs and mechanically transformed schistosomula (MTS) cultured in vitro for 3.5h or 1 or 3days. We identified, via in silico analysis, 17 Ub-E2 enzymes with conserved structural characteristics: the beta-sheet and the helix-2 form a central core bordered by helix-1 at one side and helix-3 and helix-4 at the other. The observed quantitative differences in the steady-state transcript levels between the cercariae and adult worms may contribute to the differential protein ubiquitination observed during the parasite's life cycle. This study is the first to identify and characterise the E2 ubiquitin conjugation family in S. mansoni and provides fundamental information regarding their molecular phylogenetics and developmental expression during intra-mammalian stages.

In silico analysis of the fucosylation-associated genome of the human blood fluke Schistosoma mansoni: cloning and characterization of the enzymes involved in GDP-L-fucose synthesis and Golgi import.

BackgroundCarbohydrate structures of surface-expressed and secreted/excreted glycoconjugates of the human blood fluke Schistosoma mansoni are key determinants that mediate host-parasite interactions in both snail and mammalian hosts. Fucose is a major constituent of these immunologically important glycans, and recent studies have sought to characterize fucosylation-associated enzymes, including the Golgi-localized fucosyltransferases that catalyze the transfer of L-fucose from a GDP-L-fucose donor to an oligosaccharide acceptor. Importantly, GDP-L-fucose is the only nucleotide-sugar donor used by fucosyltransferases and its availability represents a bottleneck in fucosyl-glycotope expression.MethodsA homology-based genome-wide bioinformatics approach was used to identify and molecularly characterize the enzymes that contribute to GDP-L-fucose synthesis and Golgi import in S. mansoni. Putative functions were further investigated through molecular phylogenetic and immunocytochemical analyses.ResultsWe identified homologs of GDP-D-mannose-4,6-dehydratase (GMD) and GDP-4-keto-6-deoxy-D-mannose-3,5-epimerase-4-reductase (GMER), which constitute a de novo pathway for GDP-L-fucose synthesis, in addition to a GDP-L-fucose transporter (GFT) that putatively imports cytosolic GDP-L-fucose into the Golgi. In silico primary sequence analyses identified characteristic Rossman loop and short-chain dehydrogenase/reductase motifs in GMD and GMER as well as 10 transmembrane domains in GFT. All genes are alternatively spliced, generating variants of unknown function. Observed quantitative differences in steady-state transcript levels between miracidia and primary sporocysts may contribute to differential glycotope expression in early larval development. Additionally, analyses of protein expression suggest the occurrence of cytosolic GMD and GMER in the ciliated epidermal plates and tegument of miracidia and primary sporocysts, respectively, which is consistent with previous localization of highly fucosylated glycotopes.ConclusionsThis study is the first to identify and characterize three key genes that are putatively involved in the synthesis and Golgi import of GDP-L-fucose in S. mansoni and provides fundamental information regarding their genomic organization, genetic variation, molecular phylogenetics, and developmental expression in intramolluscan larval stages.

A database for poplar gene co-expression analysis for systematic understanding of biological processes, including stress responses

Reforestation in the humid tropics and arid zones, where trees are often subject to stresses, is an effective approach for mitigating global warming. Forestation with Populus species that are tolerant to the stresses in such regions has been conducted. The selection of poplar trees with higher stress tolerance leads to more efficient reforestation. The genome-wide bioinformatics approaches of gene function have been used for revealing the mechanisms of biological processes, including such stress tolerance. The decoding of the poplar genome has been followed by the genome-wide identification of genes and then the inference of gene function for systematic understanding of biological processes. To predict gene function in poplar, we analyzed poplar gene expression data using DNA microarray datasets obtained from the Gene Expression Omnibus database and developed a database for poplar gene co-expression analysis. Using the database, we illustrate the steps to retrieve two groups of co-expressed genes that are specifically expressed in experiments of hypoxic stress response in gray poplar, a flooding-tolerant tree species. Our database allows users to extract genes involved in biological processes, such as stress reaction, and then is useful for understanding such mechanisms in tree species.

Activity, Splice Variants, Conserved Peptide Motifs, and Phylogeny of Two New α1,3-Fucosyltransferase Families (FUT10 and FUT11)

We report the cloning of three splice variants of the FUT10 gene, encoding for active alpha-l-fucosyltransferase-isoforms of 391, 419, and 479 amino acids, and two splice variants of the FUT11 gene, encoding for two related alpha-l-fucosyltransferases of 476 and 492 amino acids. The FUT10 and FUT11 appeared 830 million years ago, whereas the other alpha1,3-fucosyltransferases emerged 450 million years ago. FUT10-391 and FUT10-419 were expressed in human embryos, whereas FUT10-479 was cloned from adult brain and was not found in embryos. Recombinant FUT10-419 and FUT10-479 have a type II trans-membrane topology and are retained in the endoplasmic reticulum (ER) by a membrane retention signal at their NH(2) termini. The FUT10-479 has, in addition, a COOH-ER membrane retention signal. The FUT10-391 is a soluble protein without a trans-membrane domain or ER retention signal that transiently localizes to the Golgi and then is routed to the lysosome. After transfection in COS7 cells, the three FUT10s and at least one FUT11, link alpha-l-fucose onto conalbumin glycopeptides and biantennary N-glycan acceptors but not onto short lactosaminyl acceptor substrates as do classical monoexonic alpha1,3-fucosyltransferases. Modifications of the innermost core GlcNAc of the N-glycan, by substitution with ManNAc or with an opened GlcNAc ring or by the addition of an alpha1,6-fucose, suggest that the FUT10 transfer is performed on the innermost GlcNAc of the core chitobiose. We can exclude alpha1,3-fucosylation of the two peripheral GlcNAcs linked to the trimannosyl core of the acceptor, because the FUT10 fucosylated biantennary N-glycan product loses both terminal GlcNAc residues after digestion with human placenta alpha-N-acetylglucosaminidase.

A gene co-expression database for understanding biological processes in soybean

The sobean genome has been decoded independently in the United States, China, South Korea, and Japan using varieties indigenous to those regions. Using an EST library of the variety indigenous to the United States, a DNA microarray technique was developed and gene expression data in microarray experiments have been accumulated in public databases such as the Gene Expression Omnibus database (GEO). Such data can be used for genome-wide bioinformatics approaches to co-expression analysis for predicting gene function. To extract co-expressed genes in soybean, we analyzed soybean gene expression data based on DNA microarray datasets obtained from GEO and developed a database for soybean gene co-expression analysis. Using the database, we illustrate the steps to retrieve co-expressed genes that may be involved in isoflavonoid biosynthesis and secondary cell wall biosynthesis. Our database is useful for extracting co-expressed genes, which may be involved in a particular biological process, and then for understanding the mechanisms of such processes in legumes.