Functional characterization of enhancer activity during a long terminal repeat's evolution.

Many transposable elements (TEs) contain transcription factor binding sites and are implicated as potential regulatory elements. However, TEs are rarely functionally tested for regulatory activity, which in turn limits our understanding of how TE regulatory activity has evolved. We systematically tested the human LTR18A subfamily for regulatory activity using massively parallel reporter assay (MPRA) and found AP-1- and CEBP-related binding motifs as drivers of enhancer activity. Functional analysis of evolutionarily reconstructed ancestral sequences revealed that LTR18A elements have generally lost regulatory activity over time through sequence changes, with the largest effects occurring owing to mutations in the AP-1 and CEBP motifs. We observed that the two motifs are conserved at higher rates than expected based on neutral evolution. Finally, we identified LTR18A elements as potential enhancers in the human genome, primarily in epithelial cells. Together, our results provide a model for the origin, evolution, and co-option of TE-derived regulatory elements.

The OGCleaner: filtering false-positive homology clusters.

Detecting homologous sequences in organisms is an essential step in protein structure and function prediction, gene annotation and phylogenetic tree construction. Heuristic methods are often employed for quality control of putative homology clusters. These heuristics, however, usually only apply to pairwise sequence comparison and do not examine clusters as a whole. We present the Orthology Group Cleaner (the OGCleaner), a tool designed for filtering putative orthology groups as homology or non-homology clusters by considering all sequences in a cluster. The OGCleaner relies on high-quality orthologous groups identified in OrthoDB to train machine learning algorithms that are able to distinguish between true-positive and false-positive homology groups. This package aims to improve the quality of phylogenetic tree construction especially in instances of lower-quality transcriptome assemblies. AVAILABILITY AND IMPLEMENTATION: https://github.com/byucsl/ogcleaner CONTACT: sfujimoto@gmail.comSupplementary information: Supplementary data are available at Bioinformatics online.

Opsins have evolved under the permanent heterozygote model: insights from phylotranscriptomics of Odonata.

Gene duplication plays a central role in adaptation to novel environments by providing new genetic material for functional divergence and evolution of biological complexity. Several evolutionary models have been proposed for gene duplication to explain how new gene copies are preserved by natural selection, but these models have rarely been tested using empirical data. Opsin proteins, when combined with a chromophore, form a photopigment that is responsible for the absorption of light, the first step in the phototransduction cascade. Adaptive gene duplications have occurred many times within the animal opsins' gene family, leading to novel wavelength sensitivities. Consequently, opsins are an attractive choice for the study of gene duplication evolutionary models. Odonata (dragonflies and damselflies) have the largest opsin repertoire of any insect currently known. Additionally, there is tremendous variation in opsin copy number between species, particularly in the long-wavelength-sensitive (LWS) class. Using comprehensive phylotranscriptomic and statistical approaches, we tested various evolutionary models of gene duplication. Our results suggest that both the blue-sensitive (BS) and LWS opsin classes were subjected to strong positive selection that greatly weakens after multiple duplication events, a pattern that is consistent with the permanent heterozygote model. Due to the immense interspecific variation and duplicability potential of opsin genes among odonates, they represent a unique model system to test hypotheses regarding opsin gene duplication and diversification at the molecular level.

Detecting false positive sequence homology: a machine learning approach.

BACKGROUND: Accurate detection of homologous relationships of biological sequences (DNA or amino acid) amongst organisms is an important and often difficult task that is essential to various evolutionary studies, ranging from building phylogenies to predicting functional gene annotations. There are many existing heuristic tools, most commonly based on bidirectional BLAST searches that are used to identify homologous genes and combine them into two fundamentally distinct classes: orthologs and paralogs. Due to only using heuristic filtering based on significance score cutoffs and having no cluster post-processing tools available, these methods can often produce multiple clusters constituting unrelated (non-homologous) sequences. Therefore sequencing data extracted from incomplete genome/transcriptome assemblies originated from low coverage sequencing or produced by de novo processes without a reference genome are susceptible to high false positive rates of homology detection. RESULTS: In this paper we develop biologically informative features that can be extracted from multiple sequence alignments of putative homologous genes (orthologs and paralogs) and further utilized in context of guided experimentation to verify false positive outcomes. We demonstrate that our machine learning method trained on both known homology clusters obtained from OrthoDB and randomly generated sequence alignments (non-homologs), successfully determines apparent false positives inferred by heuristic algorithms especially among proteomes recovered from low-coverage RNA-seq data. Almost ~42 % and ~25 % of predicted putative homologies by InParanoid and HaMStR respectively were classified as false positives on experimental data set. CONCLUSIONS: Our process increases the quality of output from other clustering algorithms by providing a novel post-processing method that is both fast and efficient at removing low quality clusters of putative homologous genes recovered by heuristic-based approaches.

Functional trajectories during innate spinal cord repair.

Adult zebrafish are capable of anatomical and functional recovery following severe spinal cord injury. Axon growth, glial bridging and adult neurogenesis are hallmarks of cellular regeneration during spinal cord repair. However, the correlation between these cellular regenerative processes and functional recovery remains to be elucidated. Whereas the majority of established functional regeneration metrics measure swim capacity, we hypothesize that gait quality is more directly related to neurological health. Here, we performed a longitudinal swim tracking study for 60 individual zebrafish spanning 8 weeks of spinal cord regeneration. Multiple swim parameters as well as axonal and glial bridging were integrated. We established rostral compensation as a new gait quality metric that highly correlates with functional recovery. Tensor component analysis of longitudinal data supports a correspondence between functional recovery trajectories and neurological outcomes. Moreover, our studies predicted and validated that a subset of functional regeneration parameters measured 1 to 2 weeks post-injury is sufficient to predict the regenerative outcomes of individual animals at 8 weeks post-injury. Our findings established new functional regeneration parameters and generated a comprehensive correlative database between various functional and cellular regeneration outputs.

Functional Trajectories during innate spinal cord repair.

Adult zebrafish are capable of anatomical and functional recovery following severe spinal cord injury. Axon growth, glial bridging and adult neurogenesis are hallmarks of cellular regeneration during spinal cord repair. However, the correlation between these cellular regenerative processes and functional recovery remains to be elucidated. Whereas the majority of established functional regeneration metrics measure swim capacity, we hypothesize that gait quality is more directly related to neurological health. Here, we performed a longitudinal swim tracking study for sixty individual zebrafish spanning eight weeks of spinal cord regeneration. Multiple swim parameters as well as axonal and glial bridging were integrated. We established rostral compensation as a new gait quality metric that highly correlates with functional recovery. Tensor component analysis of longitudinal data supports a correspondence between functional recovery trajectories and neurological outcomes. Moreover, our studies predicted and validated that a subset of functional regeneration parameters measured 1 to 2 weeks post-injury is sufficient to predict the regenerative outcomes of individual animals at 8 weeks post-injury. Our findings established new functional regeneration parameters and generated a comprehensive correlative database between various functional and cellular regeneration outputs.

Fractality of tics as a quantitative assessment tool for Tourette syndrome.

Tics manifest as brief, purposeless and unintentional movements or noises that, for many individuals, can be suppressed temporarily with effort. Previous work has hypothesized that the chaotic temporal nature of tics could possess an inherent fractality, that is, have neighbour-to-neighbour correlation at all levels of timescale. However, demonstrating this phenomenon has eluded researchers for more than two decades, primarily because of the challenges associated with estimating the scale-invariant, power law exponent-called the fractal dimension Df-from fractional Brownian noise. Here, we confirm this hypothesis and establish the fractality of tics by examining two tic time series datasets collected 6-12 months apart in children with tics, using random walk models and directional statistics. We find that Df is correlated with tic severity as measured by the YGTTS total tic score, and that Df is a sensitive parameter in examining the effect of several tic suppression conditions on the tic time series. Our findings pave the way for using the fractal nature of tics as a robust quantitative tool for estimating tic severity and treatment effectiveness, as well as a possible marker for differentiating typical from functional tics.