3'-UTR SIRF: a database for identifying clusters of whort interspersed repeats in 3' untranslated regions

Expression Localization of the KRT32 Gene and Its Association of Genetic Variation with Wool Traits

Changes in keratin gene expression and spatiotemporal regulation determine the compositional content and cellular localization of wool keratin, thereby affecting wool traits. Therefore, keratin gene family member 32 (KRT32) was selected for a study using RT-qPCR, immunofluorescence, and penta-primer amplification refractory mutation system (PARMS) techniques. The results showed that KRT32 mRNA was highly expressed in the skin and localized to the inner root sheath (IRS), outer root sheath (ORS) and dermal papilla (DP). Sequencing results identified eight SNPs in KRT32, and association analyses revealed that the variations were significantly associated with multiple traits in wool (p < 0.05), including MFD, CF and MFC. The constructed haplotype combination H2H3 has higher CF and smaller MFD than other haplotype combination (p < 0.05). In conclusion, KRT32 can be used as a candidate gene for molecular genetic improvement of wool in Gansu Alpine Fine-wool sheep.

Identification of the Keratin-Associated Protein 22-2 Gene in the Capra hircus and Association of Its Variation with Cashmere Traits

The Cashmere goat is an excellent local goat breed in Gansu Province of China, and it is expected to improve cashmere production and cashmere quality through selection and breeding to enhance its commercial value. Keratin-associated proteins (KAPs) play an important role in maintaining wool structure. The gene encoding the keratin-associated protein 22-2 (KAP22-2) gene has been identified in selected species other than goats, such as humans, mice, and sheep. In this study, the sequence of the sheep KAP22-2 gene (KRTAP22-2) was aligned into the goat genome, and the sequence with the highest homology was assumed to be the goat KRTAP22-2 sequence and used to design primers to amplify the goat gene sequence. A total of 356 Longdong Cashmere goats (Gansu Province, China) were used for screening of genetic variants. Four specific bands were detected by polymerase chain reaction-single-stranded conformational polymorphism (PCR-SSCP) analysis, and they formed a total of six band types individually or in combination. Four alleles were identified by DNA sequencing of PCR amplification products. A total of four single nucleotide polymorphic sites (SNPs) were detected in the four sequenced KRTAP22-2 alleles. Two of them are in the 5'UTR region and the other two are in the coding region, and the variants in the coding region are all non-synonymous mutations. In addition, there was a 6 bp length variation in allele C. The gene was expressed in the cortical layer of primary and secondary hair follicles, the inner root sheath, as well as hair papillae and hair maternal cells in goats. The results of the correlation analysis between genotypes and cashmere traits showed that after excluding genotypes with a gene frequency of less than 5%, the mean fiber diameter (MFD) of cashmere was significantly higher in the AB genotype than in the AA and AC genotypes. That is, the KRTAP22-2 gene variants are associated with mean fiber diameter in cashmere. The above results suggest that the goat KRTAP22-2 variant can be utilized as a molecular marker candidate gene for cashmere traits.

Native functions of short tandem repeats

Over a third of the human genome is comprised of repetitive sequences, including more than a million short tandem repeats (STRs). While studies of the pathologic consequences of repeat expansions that cause syndromic human diseases are extensive, the potential native functions of STRs are often ignored. Here, we summarize a growing body of research into the normal biological functions for repetitive elements across the genome, with a particular focus on the roles of STRs in regulating gene expression. We propose reconceptualizing the pathogenic consequences of repeat expansions as aberrancies in normal gene regulation. From this altered viewpoint, we predict that future work will reveal broader roles for STRs in neuronal function and as risk alleles for more common human neurological diseases.

Deep-BSC: Predicting Raw DNA Binding Pattern in Arabidopsis Thaliana

Background:

With the rapid development of the sequencing methods in recent years, binding sites have been systematically identified in such projects as Nested-MICA and MEME. Prediction of DNA motifs with higher accuracy and precision has been a very important task for bioinformaticians. Nevertheless, experimental approaches are still time-consuming for big data set, making computational identification of binding sites indispensable.

Objective:

To facilitate the identification of the binding site, we proposed a deep learning architecture, named Deep-BSC (Deep-Learning Binary Search Classification), to predict binding sites in a raw DNA sequence with more precision and accuracy.

Methods:

Our proposed architecture purely relies on the raw DNA sequence to predict the binding sites for protein by using a convolutional neural network (CNN). We trained our deep learning model on binding sites at the nucleotide level. DNA sequence of A. thaliana is used in this study because it is a model plant.

Results:

The results demonstrate the effectiveness and efficiency of our method in the classification of binding sites against random sequences, using deep learning. We construct a CNN with different layers and filters to show the usefulness of max-pooling technique in the proposed method. To gain the interpretability of our approach, we further visualized binding sites in the saliency map and successfully identified similar motifs in the raw sequence. The proposed computational framework is time and resource efficient.

Conclusion:

Deep-BSC enables the identification of binding sites in the DNA sequences via a highly accurate CNN. The proposed computational framework can also be applied to problems such as operator, repeats in the genome, DNA markers, and recognition sites for enzymes, thereby promoting the use of Deep-BSC method in life sciences.

Drosophila 3' UTRs are more complex than protein-coding sequences

The 3′ UTRs of eukaryotic genes participate in a variety of post-transcriptional (and some transcriptional) regulatory interactions. Some of these interactions are well characterised, but an undetermined number remain to be discovered. While some regulatory sequences in 3′ UTRs may be conserved over long evolutionary time scales, others may have only ephemeral functional significance as regulatory profiles respond to changing selective pressures. Here we propose a sensitive segmentation methodology for investigating patterns of composition and conservation in 3′ UTRs based on comparison of closely related species. We describe encodings of pairwise and three-way alignments integrating information about conservation, GC content and transition/transversion ratios and apply the method to three closely related Drosophila species: D. melanogaster, D. simulans and D. yakuba. Incorporating multiple data types greatly increased the number of segment classes identified compared to similar methods based on conservation or GC content alone. We propose that the number of segments and number of types of segment identified by the method can be used as proxies for functional complexity. Our main finding is that the number of segments and segment classes identified in 3′ UTRs is greater than in the same length of protein-coding sequence, suggesting greater functional complexity in 3′ UTRs. There is thus a need for sustained and extensive efforts by bioinformaticians to delineate functional elements in this important genomic fraction. C code, data and results are available upon request.

RNA transport and long-term memory storage

Several studies have shown that synthesis of new proteins at the synapse is a prerequisite for the storage of long-term memories. Relatively little is known about the availability of distinct mRNA populations for translation at specific synapses, the process that determines mRNA localization, and the temporal designations of localized mRNA translation during memory storage. Techniques such as synaptosome preparation and microdissection of distal neuronal processes of cultured neurons and dendritic layers in brain slices are general approaches used to identify localized RNAs. Exploration of the association of RNA-binding proteins to the axonal transport machinery has led to the development of a strategy to identify RNAs that are transported from the cell body to synapses by molecular motor kinesin. In this article, RNA localization at the synapse, as well as its mechanisms and significance in understanding long-term memory storage, are discussed.

Tuning the engine: an introduction to resources on post-transcriptional regulation of gene expression

In the last years post-transcriptional regulation (PTR) of gene expression has been increasingly recognized to be a powerful and general determinant of the quantitative changes in proteomes, and therefore a driving force for cell phenotypes. By means of networks of trans-factors on one hand, and cis-elements found primarily in untranslated regions (UTRs) of mRNA on the other hand, mRNA availability to translation and translation rates are tightly controlled and can be rapidly tuned according to the changing state of the cell. A number of dedicated resources and tools, including databases and predictive algorithms, have been proposed as bioinformatics aids for the study of this fundamental layer of gene expression regulation. Their use, however, is rendered difficult by heterogeneity and fragmentation.   This review aims to locate these resources in their proper space, classifying them according to their goals, limitations and integration capabilities and, in the end, to provide the user with an initial toolbox for the bioinformatic analysis of post-transcriptional regulation of gene expression. The accompanying website, available at www.ptrguide.org, lists all resources, provides summary and features for each one and will be regularly updated in the future.

Full-length 3'-untranslated region reporter construction with recombineering

Complexity in higher animals derives in part from various modalities of protein-coding gene expression regulation, including microRNA repression by binding to 3'-untranslated regions (UTRs) of specific genes. Reporter constructs containing candidate microRNA target sites are a popular approach of functional studies, and full-length 3'-UTR sequences are preferred because they contain all regulatory elements and preserve higher order structure as much as possible. However, this approach is often handicapped by the extreme length of the 3'-UTR. Here, we present a rapid and accurate cloning procedure to generate full-length 3'-UTR reporter constructs by recombinogenic engineering (recombineering) in vivo cloning. The approach includes making retrieval constructs by sequence- and ligation-independent cloning (SLIC) and retrieving the full-length 3'-UTR in one exon to the retrieval construct from a bacterial artificial chromosome (BAC) by recombineering to generate the final full-length 3'-UTR reporter construct for the gene of interest. This method is successfully implemented with mouse full-length 3'-UTRs of Igf1 (6.5 kb), Igf1r (7.5 kb), and Sp1 (5.5 kb). Expansion of this method is adaptable to retrieve 3'-UTRs encoded in more than one exon by removing the introns from the BAC first with recombineering. This method will advance functional studies of regulation of gene expression at the post-transcriptional level through microRNA suppression.

Isoform diversity and its importance for axon regeneration

Axon regeneration is a fundamental problem facing neuroscientists and clinicians. Failure of axon regeneration is caused by both extrinsic and intrinsic mechanisms. New techniques to examine gene expression such as Next Generation Sequencing of the Transcriptome (RNA-Seq) drastically increase our knowledge of both gene expression complexity (RNA isoforms) and gene expression regulation. By utilizing RNA-Seq, gene expression can now be defined at the level of isoforms, an essential step for understanding the mechanisms governing cell identity, growth and ultimately cellular responses to injury and disease.

Identification of germ plasm-associated transcripts by microarray analysis of Xenopus vegetal cortex RNA

RNA localization is a common mechanism for regulating cell structure and function. Localized RNAs in Xenopus oocytes are critical for early development, including germline specification by the germ plasm. Despite the importance of these localized RNAs, only approximately 25 have been identified and fewer are functionally characterized. Using microarrays, we identified a large set of localized RNAs from the vegetal cortex. Overall, our results indicate a minimum of 275 localized RNAs in oocytes, or 2-3% of maternal transcripts, which are in general agreement with previous findings. We further validated vegetal localization for 24 candidates and further characterized three genes expressed in the germ plasm. We identified novel germ plasm expression for reticulon 3.1, exd2 (a novel exonuclease-domain encoding gene), and a putative noncoding RNA. Further analysis of these and other localized RNAs will likely identify new functions of germ plasm and facilitate the identification of cis-acting RNA localization elements.

Genomic features of the human dopamine transporter gene and its potential epigenetic States: implications for phenotypic diversity

Human dopamine transporter gene (DAT1 or SLC6A3) has been associated with various brain-related diseases and behavioral traits and, as such, has been investigated intensely in experimental- and clinical-settings. However, the abundance of research data has not clarified the biological mechanism of DAT regulation; similarly, studies of DAT genotype-phenotype associations yielded inconsistent results. Hence, our understanding of the control of the DAT protein product is incomplete; having this knowledge is critical, since DAT plays the major role in the brain's dopaminergic circuitry. Accordingly, we reevaluated the genomic attributes of the SLC6A3 gene that might confer sensitivity to regulation, hypothesizing that its unique genomic characteristics might facilitate highly dynamic, region-specific DAT expression, so enabling multiple regulatory modes. Our comprehensive bioinformatic analyzes revealed very distinctive genomic characteristics of the SLC6A3, including high inter-individual variability of its sequence (897 SNPs, about 90 repeats and several CNVs spell out all abbreviations in abstract) and pronounced sensitivity to regulation by epigenetic mechanisms, as evident from the GC-bias composition (0.55) of the SLC6A3, and numerous intragenic CpG islands (27 CGIs). We propose that this unique combination of the genomic features and the regulatory attributes enables the differential expression of the DAT1 gene and fulfills seemingly contradictory demands to its regulation; that is, robustness of region-specific expression and functional dynamics.

Microtubule-dependent mRNA transport in fungi

The localization and local translation of mRNAs constitute an important mechanism to promote the correct subcellular targeting of proteins. mRNA localization is mediated by the active transport of mRNPs, large assemblies consisting of mRNAs and associated factors such as RNA-binding proteins. Molecular motors move mRNPs along the actin or microtubule cytoskeleton for short-distance or long-distance trafficking, respectively. In filamentous fungi, microtubule-based long-distance transport of vesicles, which are involved in membrane and cell wall expansion, supports efficient hyphal growth. Recently, we discovered that the microtubule-mediated transport of mRNAs is essential for the fast polar growth of infectious filaments in the corn pathogen Ustilago maydis. Combining in vivo UV cross-linking and RNA live imaging revealed that the RNA-binding protein Rrm4, which constitutes an integral part of the mRNP transport machinery, mediates the transport of distinct mRNAs encoding polarity factors, protein synthesis factors, and mitochondrial proteins. Moreover, our results indicate that microtubule-dependent mRNA transport is evolutionarily conserved from fungi to higher eukaryotes. This raises the exciting possibility of U. maydis as a model system to uncover basic concepts of long-distance mRNA transport.

Evolutionary dynamics of a conserved sequence motif in the ribosomal genes of the ciliate Paramecium

Background

In protozoa, the identification of preserved motifs by comparative genomics is often impeded by difficulties to generate reliable alignments for non-coding sequences. Moreover, the evolutionary dynamics of regulatory elements in 3' untranslated regions (both in protozoa and metazoa) remains a virtually unexplored issue.

Results

By screening Paramecium tetraurelia's 3' untranslated regions for 8-mers that were previously found to be preserved in mammalian 3' UTRs, we detect and characterize a motif that is distinctly conserved in the ribosomal genes of this ciliate. The motif appears to be conserved across Paramecium aurelia species but is absent from the ribosomal genes of four additional non-Paramecium species surveyed, including another ciliate, Tetrahymena thermophila. Motif-free ribosomal genes retain fewer paralogs in the genome and appear to be lost more rapidly relative to motif-containing genes. Features associated with the discovered preserved motif are consistent with this 8-mer playing a role in post-transcriptional regulation.

Conclusions

Our observations 1) shed light on the evolution of a putative regulatory motif across large phylogenetic distances; 2) are expected to facilitate the understanding of the modulation of ribosomal genes expression in Paramecium; and 3) reveal a largely unexplored--and presumably not restricted to Paramecium--association between the presence/absence of a DNA motif and the evolutionary fate of its host genes.

Triplet repeat length bias and variation in the human transcriptome

Length variation in short tandem repeats (STRs) is an important family of DNA polymorphisms with numerous applications in genetics, medicine, forensics, and evolutionary analysis. Several major diseases have been associated with length variation of trinucleotide (triplet) repeats including Huntington's disease, hereditary ataxias and spinobulbar muscular atrophy. Using the reference human genome, we have catalogued all triplet repeats in genic regions. This data revealed a bias in noncoding DNA repeat lengths. It also enabled a survey of repeat-length polymorphisms (RLPs) in human genomes and a comparison of the rate of polymorphism in humans versus divergence from chimpanzee. For short repeats, this analysis of three human genomes reveals a relatively low RLP rate in exons and, somewhat surprisingly, in introns. All short RLPs observed in multiple genomes are biallelic (at least in this small sample). In contrast, long repeats are highly polymorphic and some long RLPs are multiallelic. For long repeats, the chimpanzee sequence frequently differs from all observed human alleles. This suggests a high expansion/contraction rate in all long repeats. Expansions and contractions are not, however, affected by natural selection discernable from our comparison of human-chimpanzee divergence with human RLPs. Our catalog of human triplet repeats and their surrounding flanking regions can be used to produce a cost-effective whole-genome assay to test individuals. This repeat assay could someday complement SNP arrays for producing tests that assess the risk of an individual to develop a disease, or become part of personalized genomic strategy that provides therapeutic guidance with respect to drug response.

The fungal RNA-binding protein Rrm4 mediates long-distance transport of ubi1 and rho3 mRNAs

Cytoskeletal transport promotes polar growth in filamentous fungi. In Ustilago maydis, the RNA-binding protein Rrm4 shuttles along microtubules and is crucial for polarity in infectious filaments. Mutations in the RNA-binding domain cause loss of function. However, it was unclear which RNAs are bound and transported. Here, we applied in vivo RNA binding studies and live imaging to determine the molecular function of Rrm4. This new combination revealed that Rrm4 mediates microtubule-dependent transport of distinct mRNAs encoding, for example, the ubiquitin fusion protein Ubi1 and the small G protein Rho3. These transcripts accumulate in ribonucleoprotein particles (mRNPs) that move bidirectionally along microtubules and co-localise with Rrm4. Importantly, the 3' untranslated region of ubi1 containing a CA-rich binding site functions as zipcode during mRNA transport. Furthermore, motile mRNPs are not formed when the RNA-binding domain of Rrm4 is deleted, although the protein is still shuttling. Thus, Rrm4 constitutes an integral component of the transport machinery. We propose that microtubule-dependent mRNP trafficking is crucial for hyphal growth introducing U. maydis as attractive model for studying mRNA transport in higher eukaryotes.

RNA localization to the Balbiani body in Xenopus oocytes is regulated by the energy state of the cell and is facilitated by kinesin II

Xenopus oocytes provide an excellent model system for understanding the cis-elements and protein factors that carry out mRNA localization in vertebrate cells. More than 20 mRNAs have been identified that localize to the vegetal cortex during stages II-IV of oogenesis. The earliest localizing RNAs are presorted to a subcellular structure, the Balbiani body (also called the mitochondrial cloud in Xenopus), of stage I oocytes prior to entering the vegetal cortex. While some evidence has suggested that diffusion drives RNA localization to the Balbiani body, a role for temperature and metabolic energy in this process has not been explored. To address this issue, we developed a quantitative assay to monitor RNA localization in stage I oocytes. Here we show that the rate of RNA accumulation to the Balbiani body is highly dependent on temperature and the intracellular concentration of ATP. In fact, while ATP depletion severely impairs RNA localization, increasing the intracellular concentration of ATP by a factor of two doubles the localization rate, indicating that ATP is limiting under normal conditions. We also show that RNA localization in stage I oocytes is reduced by inhibition of kinesin II, and that the Xcat-2 RNA localization element recruits kinesin II to the Balbiani body. We conclude from these studies that the energy state of the cell regulates the rate of RNA localization to the Balbiani body and that this process, at least to some extent, involves kinesin II.