The accurate mapping of reads that span splice junctions is a critical component of all analytic techniques that work with RNA-seq data. We introduce a second generation splice detection algorithm, ...MapSplice, whose focus is high sensitivity and specificity in the detection of splices as well as CPU and memory efficiency. MapSplice can be applied to both short (<75 bp) and long reads (≥75 bp). MapSplice is not dependent on splice site features or intron length, consequently it can detect novel canonical as well as non-canonical splices. MapSplice leverages the quality and diversity of read alignments of a given splice to increase accuracy. We demonstrate that MapSplice achieves higher sensitivity and specificity than TopHat and SpliceMap on a set of simulated RNA-seq data. Experimental studies also support the accuracy of the algorithm. Splice junctions derived from eight breast cancer RNA-seq datasets recapitulated the extensiveness of alternative splicing on a global level as well as the differences between molecular subtypes of breast cancer. These combined results indicate that MapSplice is a highly accurate algorithm for the alignment of RNA-seq reads to splice junctions. Software download URL: http://www.netlab.uky.edu/p/bioinfo/MapSplice.
Genomic variants in both coding and non-coding sequences can have functionally important and sometimes deleterious effects on exon splicing of gene transcripts. For transcriptome profiling using ...RNA-seq, the accurate alignment of reads across exon junctions is a critical step. Existing algorithms that utilize a standard reference genome as a template sometimes have difficulty in mapping reads that carry genomic variants. These problems can lead to allelic ratio biases and the failure to detect splice variants created by splice site polymorphisms. To improve RNA-seq read alignment, we have developed a novel approach called iMapSplice that enables personalized mRNA transcriptome profiling. The algorithm makes use of personal genomic information and performs an unbiased alignment towards genome indices carrying both reference and alternative bases. Importantly, this breaks the dependency on reference genome splice site dinucleotide motifs and enables iMapSplice to discover personal splice junctions created through splice site polymorphisms. We report comparative analyses using a number of simulated and real datasets. Besides general improvements in read alignment and splice junction discovery, iMapSplice greatly alleviates allelic ratio biases and unravels many previously uncharacterized splice junctions created by splice site polymorphisms, with minimal overhead in computation time and storage. Software download URL: https://github.com/LiuBioinfo/iMapSplice.
Abstract
Background
Bone marrow (BM)- and adipose tissue (AT)-derived mesenchymal stromal cells (MSCs) have shown potential as cell-based therapies for cartilage and bone injuries and are used ...increasingly in human and veterinary practice to facilitate the treatment of orthopedic conditions. However, human and rodent studies have documented a sharp decline in chondrogenic and osteogenic differentiation potential with increasing donor age, which may be problematic for the important demographic of older orthopedic patients. The aim of this study was to identify the effect of donor age on the chondrogenic and osteogenic differentiation performance of equine BM- and AT-MSCs in vitro.
BM- and AT-MSCs and dermal fibroblasts (biological negative control) were harvested from horses in five different age groups (
n
= 4,
N
= 60); newborn (0 days), yearling (15–17 months), adult (5–8 years), middle-aged (12–18 years), and geriatric (≥ 22 years). Chondrogenic differentiation performance was assessed quantitatively by measuring pellet size, matrix proteoglycan levels, and gene expression of articular cartilage biomarkers. Osteogenic differentiation performance was assessed quantitatively by measuring alkaline phosphatase activity, calcium deposition, and gene expression of bone biomarkers.
Results
Chondrogenic and osteogenic differentiation performance of equine BM- and AT-MSCs declined with increasing donor age. BM-MSCs had a higher chondrogenic differentiation performance. AT-MSCs showed minimal chondrogenic differentiation performance in all age groups. For osteogenesis, alkaline phosphatase activity was also higher in BM-MSCs, but BM-MSCs calcium deposition was affected by donor age earlier than AT-MSCs. Chondrogenic and osteogenic differentiation performance of BM-MSCs exhibited a decline as early as between the newborn and yearling samples. Steady state levels of mRNA encoding growth factors, chondrogenic, and osteogenic biomarkers were lower with increasing donor age in both MSC types.
Conclusions
The data showed that chondrogenic and osteogenic differentiation performance of equine BM-MSCs declined already in yearlings, and that AT-MSCs showed minimal chondrogenic potential, but were affected later by donor age with regards to osteogenesis (calcium deposition). The results highlight the importance of donor age considerations and MSC selection for cell-based treatment of orthopedic injuries and will help inform clinicians on when to implement or potentially cryopreserve cells. Moreover, the study provides molecular targets affected by donor age.
Gene-annotation enrichment is a common method for utilizing ontology-based annotations in gene and gene-product centric knowledgebases. Effective utilization of these annotations requires inferring ...semantic linkages by tracing paths through edges in the ontological graph, referred to as relations. However, some relations are semantically problematic with respect to scope, necessitating their omission or modification lest erroneous term mappings occur. To address these issues, we created the Gene Ontology Categorization Suite, or GOcats-a novel tool that organizes the Gene Ontology into subgraphs representing user-defined concepts, while ensuring that all appropriate relations are congruent with respect to scoping semantics. Here, we demonstrate the improvements in annotation enrichment by re-interpreting edges that would otherwise be omitted by traditional ancestor path-tracing methods. Specifically, we show that GOcats' unique handling of relations improves enrichment over conventional methods in the analysis of two different gene-expression datasets: a breast cancer microarray dataset and several horse cartilage development RNAseq datasets. With the breast cancer microarray dataset, we observed significant improvement (one-sided binomial test p-value = 1.86E-25) in 182 of 217 significantly enriched GO terms identified from the conventional path traversal method when GOcats' path traversal was used. We also found new significantly enriched terms using GOcats, whose biological relevancy has been experimentally demonstrated elsewhere. Likewise, on the horse RNAseq datasets, we observed a significant improvement in GO term enrichment when using GOcat's path traversal: one-sided binomial test p-values range from 1.32E-03 to 2.58E-44.
Significance The domestication of the horse revolutionized warfare, trade, and the exchange of people and ideas. This at least 5,500-y-long process, which ultimately transformed wild horses into the ...hundreds of breeds living today, is difficult to reconstruct from archeological data and modern genetics alone. We therefore sequenced two complete horse genomes, predating domestication by thousands of years, to characterize the genetic footprint of domestication. These ancient genomes reveal predomestic population structure and a significant fraction of genetic variation shared with the domestic breeds but absent from Przewalski’s horses. We find positive selection on genes involved in various aspects of locomotion, physiology, and cognition. Finally, we show that modern horse genomes contain an excess of deleterious mutations, likely representing the genetic cost of domestication.
The domestication of the horse ∼5.5 kya and the emergence of mounted riding, chariotry, and cavalry dramatically transformed human civilization. However, the genetics underlying horse domestication are difficult to reconstruct, given the near extinction of wild horses. We therefore sequenced two ancient horse genomes from Taymyr, Russia (at 7.4- and 24.3-fold coverage), both predating the earliest archeological evidence of domestication. We compared these genomes with genomes of domesticated horses and the wild Przewalski’s horse and found genetic structure within Eurasia in the Late Pleistocene, with the ancient population contributing significantly to the genetic variation of domesticated breeds. We furthermore identified a conservative set of 125 potential domestication targets using four complementary scans for genes that have undergone positive selection. One group of genes is involved in muscular and limb development, articular junctions, and the cardiac system, and may represent physiological adaptations to human utilization. A second group consists of genes with cognitive functions, including social behavior, learning capabilities, fear response, and agreeableness, which may have been key for taming horses. We also found that domestication is associated with inbreeding and an excess of deleterious mutations. This genetic load is in line with the “cost of domestication” hypothesis also reported for rice, tomatoes, and dogs, and it is generally attributed to the relaxation of purifying selection resulting from the strong demographic bottlenecks accompanying domestication. Our work demonstrates the power of ancient genomes to reconstruct the complex genetic changes that transformed wild animals into their domesticated forms, and the population context in which this process took place.
Purkinje cell (PC) loss occurs at an early age in patients and animal models of Niemann-Pick Type C (NPC), a lysosomal storage disease caused by mutations in the Npc1 or Npc2 genes. Although ...degeneration of PCs occurs early in NPC, little is known about how NPC1 deficiency affects the postnatal development of PCs. Using the Npc1
mouse model, we found that NPC1 deficiency significantly affected the postnatal development of PC dendrites and synapses. The developing dendrites of Npc1
PCs were significantly deficient in mitochondria and lysosomes. Furthermore, anabolic (mTORC1) and catabolic (TFEB) signaling pathways were not only perturbed but simultaneously activated in NPC1-deficient PCs, suggesting a loss of metabolic balance. We also found that mice with conditional heterozygous deletion of the Phosphatase and Tensin Homolog Deleted on Chromosome 10 gene (Pten-cHet), an inhibitor of mTORC1, showed similar early dendritic alterations in PCs to those found in Npc1-deficient mice. However, in contrast to Npc1
mice, Pten-cHet mice exhibited the overactivation of the mTORC1 pathway but with a strong inhibition of TFEB signaling, along with no dendritic mitochondrial reductions by the end of their postnatal development. Our data suggest that disruption of the lysosomal-metabolic signaling in PCs causes dendritic and synaptic developmental deficits that precede and promote their early degeneration in NPC.
Recent advances in genomic sequencing technology and computational assembly methods have allowed scientists to improve reference genome assemblies in terms of contiguity and composition. EquCab2, a ...reference genome for the domestic horse, was released in 2007. Although of equal or better quality compared to other first-generation Sanger assemblies, it had many of the shortcomings common to them. In 2014, the equine genomics research community began a project to improve the reference sequence for the horse, building upon the solid foundation of EquCab2 and incorporating new short-read data, long-read data, and proximity ligation data. Here, we present EquCab3. The count of non-N bases in the incorporated chromosomes is improved from 2.33 Gb in EquCab2 to 2.41 Gb in EquCab3. Contiguity has also been improved nearly 40-fold with a contig N50 of 4.5 Mb and scaffold contiguity enhanced to where all but one of the 32 chromosomes is comprised of a single scaffold.
Current gene annotation of the horse genome is largely derived from in silico predictions and cross-species alignments. Only a small number of genes are annotated based on equine EST and mRNA ...sequences. To expand the number of equine genes annotated from equine experimental evidence, we sequenced mRNA from a pool of forty-three different tissues. From these, we derived the structures of 68,594 transcripts. In addition, we identified 301,829 positions with SNPs or small indels within these transcripts relative to EquCab2. Interestingly, 780 variants extend the open reading frame of the transcript and appear to be small errors in the equine reference genome, since they are also identified as homozygous variants by genomic DNA resequencing of the reference horse. Taken together, we provide a resource of equine mRNA structures and protein coding variants that will enhance equine and cross-species transcriptional and genomic comparisons.
Within developing synovial joints, interzone and anlagen cells progress through divergent chondrogenic pathways to generate stable articular cartilage and transient hypertrophic anlagen cartilage, ...respectively. Understanding the comparative cell biology between interzone and anlagen cells may provide novel insights into emergent cell-based therapies to support articular cartilage regeneration. The aim of this study was to assess the kinetics of gene expression profiles in these skeletal cell lines after inducing chondrogenesis in culture. Interzone and anlagen cells from seven equine fetuses were isolated and grown in a TGF-β1 chondrogenic inductive medium. Total RNA was isolated at ten time points (0, 1.5, 3, 6, 12, 24, 48, 96, 168, and 336 h), and gene expression for 93 targeted gene loci was measured in a microfluidic RT-qPCR system. Differential transcriptional responses were observed as early as 1.5 h after the initiation of chondrogenesis. Genes with functional annotations that include transcription regulation responded to the chondrogenic stimulation earlier (1.5–96 h) than genes involved in signal transduction (1.5–336 h) and the extracellular matrix biology (3–336 h). Between interzone and anlagen cell cultures, expression levels of 73 out of the 93 targeted genes were not initially different at 0 h, but 47 out of the 73 genes became differentially expressed under the chondrogenic stimulation. While interzone and anlagen cells are both chondrogenic, they display clear differences in response to the same TGF-β1 chondrogenic stimulation. This study provides new molecular insight into a timed sequence of the divergent developmental fates of interzone and anlagen cells in culture over 14 days.
Osteoarthritis (OA) may result from impaired ability of synovial macrophages to resolve joint inflammation. Increasing macrophage counts in inflamed joints through injection with bone marrow ...mononuclear cells (BMNC) induces lasting resolution of synovial inflammation. To uncover mechanisms by which BMNC may affect resolution, in this study, differential transcriptional signatures of BMNC in response to normal (SF) and inflamed synovial fluid (ISF) were analyzed. We demonstrate the temporal behavior of co-expressed gene networks associated with traits from related
and
studies. We also identified activated and inhibited signaling pathways and upstream regulators, further determining their protein expression in the synovium of inflamed joints treated with BMNC or DPBS controls. BMNC responded to ISF with an early pro-inflammatory response characterized by a short spike in the expression of a NF-ƙB- and mitogen-related gene network. This response was associated with sustained increased expression of two gene networks comprising known drivers of resolution (
, isoprenoid biosynthesis). These networks were common to SF and ISF, but more highly expressed in ISF. Most highly activated pathways in ISF included the mevalonate pathway and PPAR-γ signaling, with pro-resolving functional annotations that improve mitochondrial metabolism and deactivate NF-ƙB signaling. Lower expression of mevalonate kinase and phospho-PPARγ in synovium from inflamed joints treated with BMNC, and equivalent IL-1β staining between BMNC- and DPBS-treated joints, associates with accomplished resolution in BMNC-treated joints and emphasize the intricate balance of pro- and anti-inflammatory mechanisms required for resolution. Combined, our data suggest that BMNC-mediated resolution is characterized by constitutively expressed homeostatic mechanisms, whose expression are enhanced following inflammatory stimulus. These mechanisms translate into macrophage proliferation optimizing their capacity to counteract inflammatory damage and improving their general and mitochondrial metabolism to endure oxidative stress while driving tissue repair. Such effect is largely achieved through the synthesis of several lipids that mediate recovery of homeostasis. Our study reveals candidate mechanisms by which BMNC provide lasting improvement in patients with OA and suggests further investigation on the effects of PPAR-γ signaling enhancement for the treatment of arthritic conditions.
PDF
