Weak correlation between sequence conservation in promoter regions and in protein-coding regions of human-mouse orthologous gene pairs

Endosomal disentanglement of a transducible artificial transcription factor targeting endothelin receptor A

Cell-penetrating peptides (CPPs) hold great promise for intracellular delivery of therapeutic proteins. However, endosomal entrapment of transduced cargo is a major bottleneck hampering their successful application. While developing a transducible zinc finger protein-based artificial transcription factor targeting the expression of endothelin receptor A, we identified interaction between the CPP and the endosomal membrane or endosomal entanglement as a main culprit for endosomal entrapment. To achieve endosomal disentanglement, we utilized endosome-resident proteases to sever the artificial transcription factor from its CPP upon arrival inside the endosome. Using this approach, we greatly enhanced the correct subcellular localization of the disentangled artificial transcription factor, significantly increasing its biological activity and distribution in vivo. With rational engineering of proteolytic sensitivity, we propose a new design principle for transducible therapeutic proteins, helping CPPs attain their full potential as delivery vectors for therapeutic proteins.

Genome-wide analysis of European sea bass provides insights into the evolution and functions of single-exon genes

Several studies have attempted to understand the origin and evolution of single-exon genes (SEGs) in eukaryotic organisms, including fishes, but few have examined the functional and evolutionary relationships between SEGs and multiple-exon gene (MEG) paralogs, in particular the conservation of promoter regions. Given that SEGs originate via the reverse transcription of mRNA from a "parental" MEGs, such comparisons may enable identifying evolutionarily-related SEG/MEG paralogs, which might fulfill equivalent physiological functions. Here, the relationship of SEG proportion with MEG count, gene density, intron count, and chromosome size was assessed for the genome of the European sea bass, Dicentrarchus labrax. Then, SEGs with an MEG parent were identified, and promoter sequences of SEG/MEG paralogs were compared, to identify highly conserved functional motifs. The results revealed a total count of 1,585 (8.3% of total genes) SEGs in the European sea bass genome, which was correlated with MEG count but not with gene density. The significant correlation of SEG content with the number of MEGs suggests that SEGs were continuously and independently generated over evolutionary time following species divergence through retrotranscription events, followed by tandem duplications. Functional annotation showed that the majority of SEGs are functional, as is evident from their expression in RNA-seq data used to support homology-based genome annotation. Differences in 5'UTR and 3'UTR lengths between SEG/MEG paralogs observed in this study may contribute to gene expression divergence between them and therefore lead to the emergence of new SEG functions. The comparison of nonsynonymous to synonymous changes (Ka/Ks) between SEG/MEG parents showed that 74 of them are under positive selection (Ka/Ks > 1; p = .0447). An additional fifteen SEGs with an MEG parent have a common promoter, which implies that they are under the influence of common regulatory networks.

Detection of long non-coding RNA homology, a comparative study on alignment and alignment-free metrics

BACKGROUND

Long non-coding RNAs (lncRNAs) represent a novel class of non-coding RNAs having a crucial role in many biological processes. The identification of long non-coding homologs among different species is essential to investigate such roles in model organisms as homologous genes tend to retain similar molecular and biological functions. Alignment-based metrics are able to effectively capture the conservation of transcribed coding sequences and then the homology of protein coding genes. However, unlike protein coding genes the poor sequence conservation of long non-coding genes makes the identification of their homologs a challenging task.

RESULTS

In this study we compare alignment-based and alignment-free string similarity metrics and look at promoter regions as a possible source of conserved information. We show that promoter regions encode relevant information for the conservation of long non-coding genes across species and that such information is better captured by alignment-free metrics. We perform a genome wide test of this hypothesis in human, mouse, and zebrafish.

CONCLUSIONS

The obtained results persuaded us to postulate the new hypothesis that, unlike protein coding genes, long non-coding genes tend to preserve their regulatory machinery rather than their transcribed sequence. All datasets, scripts, and the prediction tools adopted in this study are available at https://github.com/bioinformatics-sannio/lncrna-homologs .

The core promoter controls basal and inducible expression of duck retinoic acid inducible gene-I (RIG-I)

Retinoic acid inducible gene-I (RIG-I) is a cytoplasmic RNA sensor for detecting a variety of RNA viruses including influenza A viruses. Detection ultimately produces Type I interferon (IFN), which stimulates expression of interferon stimulated genes (ISGs), including RIG-I itself in a positive feedback loop. The structure and function of RIG-I is conserved across phylogeny, despite significant protein sequence divergence, however, the promoter sequences do not show the expected phylogenetic relationships and it is not known whether they are similarly regulated. We previously cloned duck RIG-I and showed it is highly induced during influenza A infection consistent with induction by the interferon produced. Here, we identified the Pekin duck RIG-I promoter and constructed promoter reporter vectors, which we transfected into duck embryonic fibroblasts or chicken DF-1 cells and tested in dual luciferase assays. We showed that activation of the Mitochondrial Antiviral Signalling (MAVS) pathway using the constitutively active N-terminal region of RIG-I or polyinosinic-polycytidylic acid (poly I:C) led to stimulation of duck RIG-I promoter activity. Using deletion constructs we showed the core promoter lies in the proximal 250 basepairs, and we identified essential cis-regulatory elements, a GC-box and an interferon-sensitive response element (ISRE), responsible for basal and inducible expression, respectively. Using mCherry-tagged interferon regulatory factors (IRFs) cloned from chickens and ducks, we show overexpression of chIRF7 induced the duck RIG-I promoter, and this required the ISRE site. Finally, we also demonstrated that overexpressed chIRF7 translocated to the nucleus, which was augmented by MAVS activation using RIG-I 2CARD. Our findings demonstrate that RIG-I expression is induced by chIRF7, in a positive regulatory loop. These studies show that the duck RIG-I promoter is appropriately regulated in chicken cells, necessary for the potential generation of transgenic chickens expressing RIG-I.

Maximal Expression of the Evolutionarily Conserved Slit2 Gene Promoter Requires Sp1

Slit2 is a neural axon guidance and chemorepellent protein that stimulates motility in a variety of cell types. The role of Slit2 in neural development and neoplastic growth and migration has been well established, while the genetic mechanisms underlying regulation of the Slit2 gene have not. We identified the core and proximal promoter of Slit2 by mapping multiple transcriptional start sites, analyzing transcriptional activity, and confirming sequence homology for the Slit2 proximal promoter among a number of species. Deletion series and transient transfection identified the Slit2 proximal promoter as within 399 base pairs upstream of the start of transcription. A crucial region for full expression of the Slit2 proximal promoter lies between 399 base pairs and 296 base pairs upstream of the start of transcription. Computer modeling identified three transcription factor-binding consensus sites within this region, of which only site-directed mutagenesis of one of the two identified Sp1 consensus sites inhibited transcriptional activity of the Slit2 proximal promoter (-399 to +253). Bioinformatics analysis of the Slit2 proximal promoter -399 base pair to -296 base pair region shows high sequence conservation over twenty-two species, and that this region follows an expected pattern of sequence divergence through evolution.

The relationship between gene isoform multiplicity, number of exons and protein divergence

At present we know that phenotypic differences between organisms arise from a variety of sources, like protein sequence divergence, regulatory sequence divergence, alternative splicing, etc. However, we do not have yet a complete view of how these sources are related. Here we address this problem, studying the relationship between protein divergence and the ability of genes to express multiple isoforms. We used three genome-wide datasets of human-mouse orthologs to study the relationship between isoform multiplicity co-occurrence between orthologs (the fact that two orthologs have more than one isoform) and protein divergence. In all cases our results showed that there was a monotonic dependence between these two properties. We could explain this relationship in terms of a more fundamental one, between exon number of the largest isoform and protein divergence. We found that this last relationship was present, although with variations, in other species (chimpanzee, cow, rat, chicken, zebrafish and fruit fly). In summary, we have identified a relationship between protein divergence and isoform multiplicity co-occurrence and explained its origin in terms of a simple gene-level property. Finally, we discuss the biological implications of these findings for our understanding of inter-species phenotypic differences.

Conservation of transcription factor binding events predicts gene expression across species

Recent technological advances have made it possible to determine the genome-wide binding sites of transcription factors (TFs). Comparisons across species have suggested a relatively low degree of evolutionary conservation of experimentally defined TF binding events (TFBEs). Using binding data for six different TFs in hepatocytes and embryonic stem cells from human and mouse, we demonstrate that evolutionary conservation of TFBEs within orthologous proximal promoters is closely linked to function, defined as expression of the target genes. We show that (i) there is a significantly higher degree of conservation of TFBEs when the target gene is expressed in both species; (ii) there is increased conservation of binding events for groups of TFs compared to individual TFs; and (iii) conserved TFBEs have a greater impact on the expression of their target genes than non-conserved ones. These results link conservation of structural elements (TFBEs) to conservation of function (gene expression) and suggest a higher degree of functional conservation than implied by previous studies.

Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.

Identification of cis-regulatory elements in the mammalian genome: the cREMaG database

Background

A growing number of gene expression-profiling datasets provides a reliable source of information about gene co-expression. In silico analyses of the properties shared among the promoters of co-expressed genes facilitates the identification of transcription factors (TFs) involved in the co-regulation of those genes. Our previous experience with microarray data led to the development of a database suitable for the examination of regulatory motifs in the promoters of co-expressed genes.

Methodology

We introduce the cREMaG (cis-Regulatory Elements in the Mammalian Genome) system designed for in silico studies of the promoter properties of co-regulated mammalian genes. The cREMaG system offers an analysis of data obtained from human, mouse, rat, bovine and canine gene expression-profiling studies. More than eight analysis parameters can be utilized in user-defined combinations. The selection of alternative transcription start sites and information about CpG islands are also available.

Conclusions

Using the cREMaG system, we successfully identified TFs mediating transcriptional responses in reference gene sets. The cREMaG system facilitates in silico studies of mammalian transcriptional gene regulation. The resource is freely available at http://www.cremag.org.

Divergence of exonic splicing elements after gene duplication and the impact on gene structures

An analysis of human exonic splicing elements in duplicated genes reveals their important role in the generation of new gene structures.

Background

The origin of new genes and their contribution to functional novelty has been the subject of considerable interest. There has been much progress in understanding the mechanisms by which new genes originate. Here we examine a novel way that new gene structures could originate, namely through the evolution of new alternative splicing isoforms after gene duplication.

Results

We studied the divergence of exonic splicing enhancers and silencers after gene duplication and the contributions of such divergence to the generation of new splicing isoforms. We found that exonic splicing enhancers and exonic splicing silencers diverge especially fast shortly after gene duplication. About 10% and 5% of paralogous exons undergo significantly asymmetric evolution of exonic splicing enhancers and silencers, respectively. When compared to pre-duplication ancestors, we found that there is a significant overall loss of exonic splicing enhancers and the magnitude increases with duplication age. Detailed examination reveals net gains and losses of exonic splicing enhancers and silencers in different copies and paralog clusters after gene duplication. Furthermore, we found that exonic splicing enhancer and silencer changes are mainly caused by synonymous mutations, though nonsynonymous changes also contribute. Finally, we found that exonic splicing enhancer and silencer divergence results in exon splicing state transitions (from constitutive to alternative or vice versa), and that the proportion of paralogous exon pairs with different splicing states also increases over time, consistent with previous predictions.

Conclusions

Our results suggest that exonic splicing enhancer and silencer changes after gene duplication have important roles in alternative splicing divergence and that these changes contribute to the generation of new gene structures.