Intron size in mammals: complexity comes to terms with economy

The evolutionary dynamics that retain long neutral genomic sequences in face of indel deletion bias: a model and its application to human introns

Insertions and deletions (indels) of short DNA segments are common evolutionary events. Numerous studies showed that deletions occur more often than insertions in both prokaryotes and eukaryotes. It raises the question why neutral sequences are not eradicated from the genome. We suggest that this is due to a phenomenon we term border-induced selection. Accordingly, a neutral sequence is bordered between conserved regions. Deletions occurring near the borders occasionally protrude to the conserved region and are thereby subject to strong purifying selection. Thus, for short neutral sequences, an insertion bias is expected. Here, we develop a set of increasingly complex models of indel dynamics that incorporate border-induced selection. Furthermore, we show that short conserved sequences within the neutrally evolving sequence help explain: (i) the presence of very long sequences; (ii) the high variance of sequence lengths; and (iii) the possible emergence of multimodality in sequence length distributions. Finally, we fitted our models to the human intron length distribution, as introns are thought to be mostly neutral and bordered by conserved exons. We show that when accounting for the occurrence of short conserved sequences within introns, we reproduce the main features, including the presence of long introns and the multimodality of intron distribution.

Intron size minimisation in teleosts

Background

Spliceosomal introns are parts of primary transcripts that are removed by RNA splicing. Although introns apparently do not contribute to the function of the mature transcript, in vertebrates they comprise the majority of the transcribed region increasing the metabolic cost of transcription. The persistence of long introns across evolutionary time suggests functional roles that can offset this metabolic cost. The teleosts comprise one of the largest vertebrate clades. They have unusually compact and variable genome sizes and provide a suitable system for analysing intron evolution.

Results

We have analysed intron lengths in 172 vertebrate genomes and show that teleost intron lengths are relatively short, highly variable and bimodally distributed. Introns that were long in teleosts were also found to be long in mammals and were more likely to be found in regulatory genes and to contain conserved sequences. Our results argue that intron length has decreased in parallel in a non-random manner throughout teleost evolution and represent a deviation from the ancestral state.

Conclusion

Our observations indicate an accelerated rate of intron size evolution in the teleosts and that teleost introns can be divided into two classes by their length. Teleost intron sizes have evolved primarily as a side-effect of genome size evolution and small genomes are dominated by short introns (<256 base pairs). However, a non-random subset of introns has resisted this process across the teleosts and these are more likely have functional roles in all vertebrate clades.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08760-w).

The transcriptional trajectories of pluripotency and differentiation comprise genes with antithetical architecture and repetitive-element content

Background

Extensive molecular differences exist between proliferative and differentiated cells. Here, we conduct a meta-analysis of publicly available transcriptomic datasets from preimplantation and differentiation stages examining the architectural properties and content of genes whose abundance changes significantly across developmental time points.

Results

Analysis of preimplantation embryos from human and mouse showed that short genes whose introns are enriched in Alu (human) and B (mouse) elements, respectively, have higher abundance in the blastocyst compared to the zygote. These highly expressed genes encode ribosomal proteins or metabolic enzymes. On the other hand, long genes whose introns are depleted in repetitive elements have lower abundance in the blastocyst and include genes from signaling pathways. Additionally, the sequences of the genes that are differentially expressed between the blastocyst and the zygote contain distinct collections of pyknon motifs that differ between up- and down-regulated genes. Further examination of the genes that participate in the stem cell-specific protein interaction network shows that their introns are short and enriched in Alu (human) and B (mouse) elements. As organogenesis progresses, in both human and mouse, we find that the primarily short and repeat-rich expressed genes make way for primarily longer, repeat-poor genes. With that in mind, we used a machine learning-based approach to identify gene signatures able to classify human adult tissues: we find that the most discriminatory genes comprising these signatures have long introns that are repeat-poor and include transcription factors and signaling-cascade genes. The introns of widely expressed genes across human tissues, on the other hand, are short and repeat-rich, and coincide with those with the highest expression at the blastocyst stage.

Conclusions

Protein-coding genes that are characteristic of each trajectory, i.e., proliferation/pluripotency or differentiation, exhibit antithetical biases in their intronic and exonic lengths and in their repetitive-element content. While the respective human and mouse gene signatures are functionally and evolutionarily conserved, their introns and exons are enriched or depleted in organism-specific repetitive elements. We posit that these organism-specific repetitive sequences found in exons and introns are used to effect the corresponding genes’ regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-020-00928-8.

Spliceosomal Introns: Features, Functions, and Evolution

Spliceosomal introns, which have been found in most eukaryotic genes, are non-coding sequences excised from pre-mRNAs by a special complex called spliceosome during mRNA splicing. Introns occur in both protein- and RNA-coding genes and can be found in coding and untranslated gene regions. Because intron sequences vary greatly due to a high rate of polymorphism, the functions of intron had been for a long time associated only with alternative splicing, while intron evolution had been viewed not as an evolution of an individual genomic element, but rather considered within a framework of the evolution of the gene intron-exon structure. Here, we review the theories of intron origin, evolutionary events in the exon-intron structure, such as intron gain, loss, and sliding, intron functions known to date, and mechanisms by which changes in the intron features (length and phase) can affect the regulation of gene-mediated processes.

Functionally Significant Features in the 5' Untranslated Region of the ABCA1 Gene and Their Comparison in Vertebrates

Single nucleotide polymorphisms located in 5′ untranslated regions (5′UTRs) can regulate gene expression and have clinical impact. Recognition of functionally significant sequences within 5′UTRs is crucial in next-generation sequencing applications. Furthermore, information about the behavior of 5′UTRs during gene evolution is scarce. Using the example of the ATP-binding cassette transporter A1 (ABCA1) gene (Tangier disease), we describe our algorithm for functionally significant sequence finding. 5′UTR features (upstream start and stop codons, open reading frames (ORFs), GC content, motifs, and secondary structures) were studied using freely available bioinformatics tools in 55 vertebrate orthologous genes obtained from Ensembl and UCSC. The most conserved sequences were suggested as hot spots. Exon and intron enhancers and silencers (sc35, ighg2 cgamma2, ctnt, gh-1, and fibronectin eda exon), transcription factors (TFIIA, TATA, NFAT1, NFAT4, and HOXA13), some of them cancer related, and microRNA (hsa-miR-4474-3p) were localized to these regions. An upstream ORF, overlapping with the main ORF in primates and possibly coding for a small bioactive peptide, was also detected. Moreover, we showed several features of 5′UTRs, such as GC content variation, hairpin structure conservation or 5′UTR segmentation, which are interesting from a phylogenetic point of view and can stimulate further evolutionary oriented research.

Dichotomy in redundant enhancers points to presence of initiators of gene regulation

Background

The regulatory landscape of a gene locus often consists of several functionally redundant enhancers establishing phenotypic robustness and evolutionary stability of its regulatory program. However, it is unclear what mechanisms are employed by redundant enhancers to cooperatively orchestrate gene expression.

Results

By comparing redundant enhancers to single enhancers (enhancers present in a single copy in a gene locus), we observed that the DNA sequence encryption differs between these two classes of enhancers, suggesting a difference in their regulatory mechanisms. Initiator enhancers, which are a subset of redundant enhancers and show similar sequence encryption to single enhancers, differ from the rest of redundant enhancers in their sequence encryption, evolutionary conservation and proximity to target genes. Genes hosting initiator enhancers in their loci feature elevated levels of expression. Initiator enhancers show a high level of 3D chromatin contacts with both transcription start sites and regular enhancers, suggesting their roles as primary activators and intermediate catalysts of gene expression, through which the regulatory signals of redundant enhancers are propagated to the target genes. In addition, GWAS and eQTLs variants are significantly enriched in initiator enhancers compared to redundant enhancers, suggesting a key functional role these sequences play in gene regulation.

Conclusions

The specific characteristics and widespread abundance of initiator enhancers advocate for a possible universal hierarchical mechanism of tissue-specific gene regulation involving multiple redundant enhancers acting through initiator enhancers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5335-0) contains supplementary material, which is available to authorized users.

Conserved Nonexonic Elements: A Novel Class of Marker for Phylogenomics

Noncoding markers have a particular appeal as tools for phylogenomic analysis because, at least in vertebrates, they appear less subject to strong variation in GC content among lineages. Thus far, ultraconserved elements (UCEs) and introns have been the most widely used noncoding markers. Here we analyze and study the evolutionary properties of a new type of noncoding marker, conserved nonexonic elements (CNEEs), which consists of noncoding elements that are estimated to evolve slower than the neutral rate across a set of species. Although they often include UCEs, CNEEs are distinct from UCEs because they are not ultraconserved, and, most importantly, the core region alone is analyzed, rather than both the core and its flanking regions. Using a data set of 16 birds plus an alligator outgroup, and ∼3600-∼3800 loci per marker type, we found that although CNEEs were less variable than bioinformatically derived UCEs or introns and in some cases exhibited a slower approach to branch resolution as determined by phylogenomic subsampling, the quality of CNEE alignments was superior to those of the other markers, with fewer gaps and missing species. Phylogenetic resolution using coalescent approaches was comparable among the three marker types, with most nodes being fully and congruently resolved. Comparison of phylogenetic results across the three marker types indicated that one branch, the sister group to the passerine + falcon clade, was resolved differently and with moderate (>70%) bootstrap support between CNEEs and UCEs or introns. Overall, CNEEs appear to be promising as phylogenomic markers, yielding phylogenetic resolution as high as for UCEs and introns but with fewer gaps, less ambiguity in alignments and with patterns of nucleotide substitution more consistent with the assumptions of commonly used methods of phylogenetic analysis.

Identification of minimal eukaryotic introns through GeneBase, a user-friendly tool for parsing the NCBI Gene databank

We have developed GeneBase, a full parser of the National Center for Biotechnology Information (NCBI) Gene database, which generates a fully structured local database with an intuitive user-friendly graphic interface for personal computers. Features of all the annotated eukaryotic genes are accessible through three main software tables, including for each entry details such as the gene summary, the gene exon/intron structure and the specific Gene Ontology attributions. The structuring of the data, the creation of additional calculation fields and the integration with nucleotide sequences allow users to make many types of comparisons and calculations that are useful for data retrieval and analysis. We provide an original example analysis of the existing introns across all the available species, through which the classic biological problem of the ‘minimal intron’ may find a solution using available data. Based on all currently available data, we can define the shortest known eukaryotic GT-AG intron length, setting the physical limit at the 30 base pair intron belonging to the human MST1L gene. This ‘model intron’ will shed light on the minimal requirement elements of recognition used for conventional splicing functioning. Remarkably, this size is indeed consistent with the sum of the splicing consensus sequence lengths.

Molecular characterization, phylogeny and expression of a hepcidin gene in the blotched snakehead Channa maculata

A hepcidin-like gene (cmHep) was cloned and characterized from the liver of the blotched snakehead Channa maculata. The complete cmHep cDNA was 756 bp in length, containing an open reading frame of 270 bp (encoding 89 amino acids), flanked by 210 bp and 276 bp of 5' and 3' untranslated regions, respectively. The deduced peptide of 89 amino acids consisted of 24 aa, 40 aa and 25 aa for signal peptide, prodomain and mature peptide, respectively. The mature peptide had eight cysteines at the identical conserved positions in common with most of other known hepcidins in vertebrates. cmHepc gene displayed a tripartite structure (three exons interrupted by two introns), which organisation was conserved between the blotched snakehead and other fish species. Phylogenetic analysis of hepcidins from C. maculata and other vertebrates showed that major phylogenetic grouping of fish hepcidin coincided with the current euteleosts classification, indicating the multiphyletic evolution of hepcidin in the teleosts. In the Acanthopterygii subclade, there were two distinct additional subclades named as HAMP-Ac1 and HAMP-Ac2. The blotched snakehead hepcidin was in the group HAMP-Ac1, which has the hypothetical iron regulatory sequence [Q-S/I-H-L/I-S/A] motif in N-terminal of mature peptide. The RT-PCR showed cmHep mRNA transcripts were widely distributed in all tissues tested in the blotched snakehead including the liver, gill, intestine, spleen, head kidney and peripheral white blood cell. The most abundant of cmHep mRNA was detected in liver. A significant up-regulation of cmHep expression was detected only in head kidney at 24h post-challenge with Vibrio parahaemolyticus in blotched snakehead adults, no significant differences found in liver, gill, intestine and spleen. The cmHep expression was up-regulated in spleen, head kidney and intestine at 24h post-injection with LPS in blotched snakehead juveniles, liver cmHep expression was not altered. Iron overloading and poly I:C stimulation down-regulated cmHep expression in liver, but did not significantly change cmHep expression in spleen, head kidney and intestine in blotched snakehead juveniles.

Genome-wide analysis of functional and evolutionary features of tele-enhancers

We investigated sequence features of enhancers separated from their target gene by at least one intermediate gene/exon (named tele-enhancers in this study) and enhancers residing inside their target gene locus. In this study, we used whole genome enhancer maps and gene expression profiles to establish a large panel of tele-enhancers. By contrasting tele-enhancers to proximal enhancers targeting heart genes, we observed that heart tele-enhancers use unique regulatory mechanisms based on the cardiac transcription factors SRF, TEAD, and NKX-2.5, whereas proximal heart enhancers rely on GATA4 instead. A functional analysis showed that tele-enhancers preferentially regulate house-keeping genes and genes with a metabolic role during heart development. In addition, tele-enhancers are significantly more conserved than their proximal counterparts. Similar trends have been observed for non-heart tissues and cell types, suggesting that our findings represent general characteristics of tele-enhancers.

DNase I-hypersensitive exons colocalize with promoters and distal regulatory elements

The precise splicing of genes confers an enormous transcriptional complexity to the human genome. The majority of gene splicing occurs cotranscriptionally, permitting epigenetic modifications to affect splicing outcomes. Here we show that select exonic regions are demarcated within the three-dimensional structure of the human genome. We identify a subset of exons that exhibit DNase I hypersensitivity and are accompanied by 'phantom' signals in chromatin immunoprecipitation and sequencing (ChIP-seq) that result from cross-linking with proximal promoter- or enhancer-bound factors. The capture of structural features by ChIP-seq is confirmed by chromatin interaction analysis that resolves local intragenic loops that fold exons close to cognate promoters while excluding intervening intronic sequences. These interactions of exons with promoters and enhancers are enriched for alternative splicing events, an effect reflected in cell type-specific periexonic DNase I hypersensitivity patterns. Collectively, our results connect local genome topography, chromatin structure and cis-regulatory landscapes with the generation of human transcriptional complexity by cotranscriptional splicing.

The evolution of intron size in amniotes: a role for powered flight?

Intronic DNA is a major component of eukaryotic genes and genomes and can be subject to selective constraint and have functions in gene regulation. Intron size is of particular interest given that it is thought to be the target of a variety of evolutionary forces and has been suggested to be linked ultimately to various phenotypic traits, such as powered flight. Using whole-genome analyses and comparative approaches that account for phylogenetic nonindependence, we examined interspecific variation in intron size variation in three data sets encompassing from 12 to 30 amniotes genomes and allowing for different levels of genome coverage. In addition to confirming that intron size is negatively associated with intron position and correlates with genome size, we found that on average mammals have longer introns than birds and nonavian reptiles, a trend that is correlated with the proliferation of repetitive elements in mammals. Two independent comparisons between flying and nonflying sister groups both showed a reduction of intron size in volant species, supporting an association between powered flight, or possibly the high metabolic rates associated with flight, and reduced intron/genome size. Small intron size in volant lineages is less easily explained as a neutral consequence of large effective population size. In conclusion, we found that the evolution of intron size in amniotes appears to be non-neutral, is correlated with genome size, and is likely influenced by powered flight and associated high metabolic rates.

Transposon-derived and satellite-derived repetitive sequences play distinct functional roles in Mammalian intron size expansion

Background

Repetitive sequences (RSs) are redundant, complex at times, and often lineage-specific, representing significant “building” materials for genes and genomes. According to their origins, sequence characteristics, and ways of propagation, repetitive sequences are divided into transposable elements (TEs) and satellite sequences (SSs) as well as related subfamilies and subgroups hierarchically. The combined changes attributable to the repetitive sequences alter gene and genome architectures, such as the expansion of exonic, intronic, and intergenic sequences, and most of them propagate in a seemingly random fashion and contribute very significantly to the entire mutation spectrum of mammalian genomes.

Principal findings

Our analysis is focused on evolutional features of TEs and SSs in the intronic sequence of twelve selected mammalian genomes. We divided them into four groups—primates, large mammals, rodents, and primary mammals—and used four non-mammalian vertebrate species as the out-group. After classifying intron size variation in an intron-centric way based on RS-dominance (TE-dominant or SS-dominant intron expansions), we observed several distinct profiles in intron length and positioning in different vertebrate lineages, such as retrotransposon-dominance in mammals and DNA transposon-dominance in the lower vertebrates, amphibians and fishes. The RS patterns of mouse and rat genes are most striking, which are not only distinct from those of other mammals but also different from that of the third rodent species analyzed in this study—guinea pig. Looking into the biological functions of relevant genes, we observed a two-dimensional divergence; in particular, genes that possess SS-dominant and/or RS-free introns are enriched in tissue-specific development and transcription regulation in all mammalian lineages. In addition, we found that the tendency of transposons in increasing intron size is much stronger than that of satellites, and the combined effect of both RSs is greater than either one of them alone in a simple arithmetic sum among the mammals and the opposite is found among the four non-mammalian vertebrates.

Conclusions

TE- and SS-derived RSs represent major mutational forces shaping the size and composition of vertebrate genes and genomes, and through natural selection they either fine-tune or facilitate changes in size expansion, position variation, and duplication, and thus in functions and evolutionary paths for better survival and fitness. When analyzed globally, not only are such changes significantly diversified but also comprehensible in lineages and biological implications.

Locating protein-coding sequences under selection for additional, overlapping functions in 29 mammalian genomes

The degeneracy of the genetic code allows protein-coding DNA and RNA sequences to simultaneously encode additional, overlapping functional elements. A sequence in which both protein-coding and additional overlapping functions have evolved under purifying selection should show increased evolutionary conservation compared to typical protein-coding genes--especially at synonymous sites. In this study, we use genome alignments of 29 placental mammals to systematically locate short regions within human ORFs that show conspicuously low estimated rates of synonymous substitution across these species. The 29-species alignment provides statistical power to locate more than 10,000 such regions with resolution down to nine-codon windows, which are found within more than a quarter of all human protein-coding genes and contain ∼2% of their synonymous sites. We collect numerous lines of evidence that the observed synonymous constraint in these regions reflects selection on overlapping functional elements including splicing regulatory elements, dual-coding genes, RNA secondary structures, microRNA target sites, and developmental enhancers. Our results show that overlapping functional elements are common in mammalian genes, despite the vast genomic landscape.

Evolutionary analysis of three gibberellin oxidase genes in rice, Arabidopsis, and soybean

GAs are plant hormones that play fundamental roles in plant growth and development. GA2ox, GA3ox, and GA20ox are three key enzymes in GA biosynthesis. These enzymes belong to the 2OG-Fe (II) oxygenase superfamily and are independently encoded by different gene families. To date, genome-wide comparative analyses of GA oxidases in plant species have not been thoroughly carried out. In the present work, 61 GA oxidase family genes from rice (Oryza sativa), Arabidopsis, and soybean (Glycine max) were identified and a full study of these genes including phylogenetic tree construction, gene structure, gene family expansion and analysis of functional motifs was performed. Based on phylogeny, most of the GA oxidases were divided into four subgroups that reflected functional classifications. Intron/intron average length of GA oxidase genes in rice analysis revealed that GA oxidase genes in rice experienced substantial evolutionary divergence. Segmental duplication events were mainly found in soybean genome. However, in rice and Arabidopsis, no single expansion pattern exhibited dominance, indicating that GA oxidase genes from these species might have been subjected to a more complex evolutionary mechanism. In addition, special functional motifs were discovered in GA20ox, GA3ox, and GA2ox, which suggested that different functional motifs are associated with differences in protein function. Taken together our results suggest that GA oxidase family genes have undergone divergent evolutionary routes, especially at the monocot-dicot split, with dynamic evolution occurring in Arabidopsis thaliana and soybean.

Gene expression patterns are correlated with genomic and genic structure in soybean

Studies have indicated that exon and intron size and intergenic distance are correlated with gene expression levels and expression breadth. Previous reports on these correlations in plants and animals have been conflicting. In this study, next-generation sequence data, which has been shown to be more sensitive than previous expression profiling technologies, were generated and analyzed from 14 tissues. Our results revealed a novel dichotomy. At the low expression level, an increase in expression breadth correlated with an increase in transcript size because of an increase in the number of exons and introns. No significant changes in intron or exon sizes were noted. Conversely, genes expressed at the intermediate to high expression levels displayed a decrease in transcript size as their expression breadth increased. This was due to smaller exons, with no significant change in the number of exons. Taking advantage of the known gene space of soybean, we evaluated the positioning of genes and found significant clustering of similarly expressed genes. Identifying the correlations between the physical parameters of individual genes could lead to uncovering the role of regulation owing to nucleotide composition, which might have potential impacts in discerning the role of the noncoding regions.

Temperature and length-dependent modulation of the MH class II beta gene expression in brook charr (Salvelinus fontinalis) by a cis-acting minisatellite

It is widely recognized that the variation in gene regulation is an important factor from which evolutionary changes in diverse aspects of phenotype can be observed in all organisms. Distinctive elements with functional roles on gene regulation have been identified within the non-coding part of the genome, including repeated elements. Major histocompatibility complex (MHC) genes have been the subject of an abundant literature which made them unique candidates for studies of adaptation in natural populations. Yet, the vast majority of studies on MHC genes have dealt with patterns of polymorphism in sequence variation while very few paid attention to the possible implication of differential expression in adaptive responses. In this paper, we report the identification of a polymorphic minisatellite formed of a 32 nucleotides motif (38% G+C) involved in regulation of the major histocompatibility class II beta gene (MHII beta) of brook charr (Salvelinus fontinalis). Our main objectives were: to analyze the variability of this minisatellite found in the second intron of the MHII beta gene and to document its effect to the variation of expression level of this gene under different environmental conditions. Distinctive number of the minisatellite repeats were associated with each different MHII beta alleles identified from exon 2 sequences. Relative expression levels of specific alleles in heterozygous individuals were determined from fish lymphocytes in different genotypes. We found that alleles carrying the longest minisatellite showed a significant 1.67-2.56-fold reduction in the transcript expression relatively to the shortest one. Results obtained in three different genotypes also indicated that the repressive activity associated to the longest minisatellite was more effective at 18 degrees C compared to 6 degrees C. In contrast, no significant difference was observed in transcript levels between alleles with comparable minisatellite length at both temperatures. We also depicted a significant up-regulation of the total MHII beta transcript at 6 degrees C relative to 18 degrees C. These results reveal for the first time that a temperature-sensitive minisatellite could potentially play an important role in the gene regulation of the adaptive immune response in fishes.

Insight into microRNA regulation by analyzing the characteristics of their targets in humans

Background

microRNAs (miRNAs) are believed to regulate their targets through posttranscriptional gene regulation and have the potential to silence gene expression via multiple mechanisms. Despite previous advances on miRNA regulation of gene expression, little has been investigated from a genome scale.

Results

To gain new insight into miRNA regulation in humans, we used large scale data and carried out a series of studies to compare various features of miRNA target genes to that of non-miRNA target genes. We observed significant differences between miRNA and non-miRNA target genes for a number of characteristics, including higher and broader mRNA expression, faster mRNA decay rate, longer protein half-life, and longer gene structures. Based on these features and by analyzing their relationships we found that miRNA target genes, other than having miRNA repression, were most likely under more complex regulation than non-miRNA target genes, which was evidenced by their higher and broader gene expression but longer gene structures. Our results of higher and broader gene expression but fast mRNA decay rates also provide evidence that miRNA dampening of the output of preexisting transcripts facilitates a more rapid and robust transition to new expression programs. This could be achieved by enhancing mRNA degradation through an additive effect from multiple miRNA targeting.

Conclusion

Genome-scale analysis on the nature of miRNA target genes has revealed a general mechanism for miRNA regulation of human gene expression. The results of this study also indicate that miRNA target genes, other than having miRNA repression, are under more complex gene regulation than non-miRNA target genes. These findings provide novel insight into miRNA regulation of human gene expression.

Entropic contributions to the splicing process

It has been recently argued that depletion attraction may play an important role in different aspects of cellular organization, ranging from the organization of transcriptional activity in transcription factories to the formation of nuclear bodies. In this paper, we suggest a new application of these ideas in the context of the splicing process, a crucial step of messenger RNA maturation in eukaryotes. We shall show that entropy effects and the resulting depletion attraction may explain the relevance of the aspecific intron length variable in the choice of splice-site recognition modality. On top of that, some qualitative features of the genome architecture of higher eukaryotes can find evolutionary realistic motivation in the light of our model.

In plants, expression breadth and expression level distinctly and non-linearly correlate with gene structure

BACKGROUND

Compactness of highly/broadly expressed genes in human has been explained as selection for efficiency, regional mutation biases or genomic design. However, highly expressed genes in flowering plants were shown to be less compact than lowly expressed ones. On the other hand, opposite facts have also been documented that pollen-expressed Arabidopsis genes tend to contain shorter introns and highly expressed moss genes are compact. This issue is important because it provides a chance to compare the selectionism and the neutralism views about genome evolution. Furthermore, this issue also helps to understand the fates of introns, from the angle of gene expression.

RESULTS

In this study, I used expression data covering more tissues and employ new analytical methods to reexamine the correlations between gene expression and gene structure for two flowering plants, Arabidopsis thaliana and Oryza sativa. It is shown that, different aspects of expression pattern correlate with different parts of gene sequences in distinct ways. In detail, expression level is significantly negatively correlated with gene size, especially the size of non-coding regions, whereas expression breadth correlates with non-coding structural parameters positively and with coding region parameters negatively. Furthermore, the relationships between expression level and structural parameters seem to be non-linear, with the extremes of structural parameters possibly scale as power-laws or logrithmic functions of expression levels.

CONCLUSION

In plants, highly expressed genes are compact, especially in the non-coding regions. Broadly expressed genes tend to contain longer non-coding sequences, which may be necessary for complex regulations. In combination with previous studies about other plants and about animals, some common scenarios about the correlation between gene expression and gene structure begin to emerge. Based on the functional relationships between extreme values of structural characteristics and expression level, an effort was made to evaluate the relative effectiveness of the energy-cost hypothesis and the time-cost hypothesis.

[Recent progress in plant genome size evolution]

It has been known that eukaryotic genomes span a wide range of sizes regardless of organism complexity. The observed differences in genome size are primarily due to polyploidy level and abundance of non-coding DNA, especially the contribution of transposable elements (TEs). Here we reviewed the current progress in genome size variation of plant species and the underlying evolutionary forces that contribute to genome expansion or contraction. Polyploidization and the accumulation of transposable element are the primary contributors to genome expansion. As to the mechanisms of DNA loss, unequal homologous recombination and illegitimate recombination are thought to be the counterbalances to the unlimited expansion of a genome. The evolutionary direction of plant genome size is also discussed, which tends to favor larger genomes with deletion mechanisms acting to only attenuate genome expansion but not reverse.

TopHat: discovering splice junctions with RNA-Seq

Motivation: A new protocol for sequencing the messenger RNA in a cell, known as RNA-Seq, generates millions of short sequence fragments in a single run. These fragments, or ‘reads’, can be used to measure levels of gene expression and to identify novel splice variants of genes. However, current software for aligning RNA-Seq data to a genome relies on known splice junctions and cannot identify novel ones. TopHat is an efficient read-mapping algorithm designed to align reads from an RNA-Seq experiment to a reference genome without relying on known splice sites.

Results: We mapped the RNA-Seq reads from a recent mammalian RNA-Seq experiment and recovered more than 72% of the splice junctions reported by the annotation-based software from that study, along with nearly 20 000 previously unreported junctions. The TopHat pipeline is much faster than previous systems, mapping nearly 2.2 million reads per CPU hour, which is sufficient to process an entire RNA-Seq experiment in less than a day on a standard desktop computer. We describe several challenges unique to ab initio splice site discovery from RNA-Seq reads that will require further algorithm development.

Availability: TopHat is free, open-source software available from http://tophat.cbcb.umd.edu

Contact:cole@cs.umd.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Evidence against the energetic cost hypothesis for the short introns in highly expressed genes

BACKGROUND

In animals, the moss Physcomitrella patens and the pollen of Arabidopsis thaliana, highly expressed genes have shorter introns than weakly expressed genes. A popular explanation for this is selection for transcription efficiency, which includes two sub-hypotheses: to minimize the energetic cost or to minimize the time cost.

RESULTS

In an individual human, different organs may differ up to hundreds of times in cell number (for example, a liver versus a hypothalamus). Considered at the individual level, a gene specifically expressed in a large organ is actually transcribed tens or hundreds of times more than a gene with a similar expression level (a measure of mRNA abundance per cell) specifically expressed in a small organ. According to the energetic cost hypothesis, the former should have shorter introns than the latter. However, in humans and mice we have not found significant differences in intron length between large-tissue/organ-specific genes and small-tissue/organ-specific genes with similar expression levels. Qualitative estimation shows that the deleterious effect (that is, the energetic burden) of long introns in highly expressed genes is too negligible to be efficiently selected against in mammals.

CONCLUSION

The short introns in highly expressed genes should not be attributed to energy constraint. We evaluated evidence for the time cost hypothesis and other alternatives.

The calcium channel beta2 (CACNB2) subunit repertoire in teleosts

BACKGROUND

Cardiomyocyte contraction is initiated by influx of extracellular calcium through voltage-gated calcium channels. These oligomeric channels utilize auxiliary beta subunits to chaperone the pore-forming alpha subunit to the plasma membrane, and to modulate channel electrophysiology 1. Several beta subunit family members are detected by RT-PCR in the embryonic heart. Null mutations in mouse beta2, but not in the other three beta family members, are embryonic lethal at E10.5 due to defects in cardiac contractility 2. However, a drawback of the mouse model is that embryonic heart rhythm is difficult to study in live embryos due to their intra-uterine development. Moreover, phenotypes may be obscured by secondary effects of hypoxia. As a first step towards developing a model for contributions of beta subunits to the onset of embryonic heart rhythm, we characterized the structure and expression of beta2 subunits in zebrafish and other teleosts.

RESULTS

Cloning of two zebrafish beta2 subunit genes (beta2.1 and beta2.2) indicated they are membrane-associated guanylate kinase (MAGUK)-family genes. Zebrafish beta2 genes show high conservation with mammals within the SH3 and guanylate kinase domains that comprise the "core" of MAGUK proteins, but beta2.2 is much more divergent in sequence than beta2.1. Alternative splicing occurs at the N-terminus and within the internal HOOK domain. In both beta2 genes, alternative short ATG-containing first exons are separated by some of the largest introns in the genome, suggesting that individual transcript variants could be subject to independent cis-regulatory control. In the Tetraodon nigrovidis and Fugu rubripes genomes, we identified single beta2 subunit gene loci. Comparative analysis of the teleost and human beta2 loci indicates that the short 5' exon sequences are highly conserved. A subset of 5' exons appear to be unique to teleost genomes, while others are shared with mammals. Alternative splicing is temporally and spatially regulated in embryo and adult. Moreover, a different subset of spliced beta2 transcript variants is detected in the embryonic heart compared to the adult.

CONCLUSION

These studies refine our understanding of beta2 subunit diversity arising from alternative splicing, and provide the groundwork for functional analysis of beta2 subunit diversity in the embryonic heart.

Both selective and neutral processes drive GC content evolution in the human genome

Background

Mammalian genomes consist of regions differing in GC content, referred to as isochores or GC-content domains. The scientific debate is still open as to whether such compositional heterogeneity is a selected or neutral trait.

Results

Here we analyze SNP allele frequencies, retrotransposon insertion polymorphisms (RIPs), as well as fixed substitutions accumulated in the human lineage since its divergence from chimpanzee to indicate that biased gene conversion (BGC) has been playing a role in within-genome GC content variation. Yet, a distinct contribution to GC content evolution is accounted for by a selective process. Accordingly, we searched for independent evidences that GC content distribution does not conform to neutral expectations. Indeed, after correcting for possible biases, we show that intron GC content and size display isochore-specific correlations.

Conclusion

We consider that the more parsimonious explanation for our results is that GC content is subjected to the action of both weak selection and BGC in the human genome with features such as nucleosome positioning or chromatin conformation possibly representing the final target of selective processes. This view might reconcile previous contrasting findings and add some theoretical background to recent evidences suggesting that GC content domains display different behaviors with respect to highly regulated biological processes such as developmentally-stage related gene expression and programmed replication timing during neural stem cell differentiation.

Excess of microRNAs in large and very 5' biased introns

Many of microRNAs (miRNAs) and small nucleolar RNAs (snoRNAs) are located within the introns of genes in eukaryotes. Contrary to intronic snoRNAs, intronic miRNAs are processed from unspliced intronic regions before the catalysis of splicing in vertebrates. By analyzing the distribution patterns of the length and position of the introns hosting these two groups of small RNA genes, we observed that both human and mouse intronic miRNAs tended to be present in large introns, and miRNA host introns have a more 5'-biased position distribution compared with all other introns among the two genomes. These observations indicate that the negative selection of functional constraints might affect the intron size in both genomes. Interestingly, the very 5'-biased positions of miRNA host introns may be necessary for the transcription and regulation of intronic miRNAs to utilize the regulatory signals within the 5'-UTRs of their host genes.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

Selection for the miniaturization of highly expressed genes

Most widely expressed genes are also highly expressed. Based on high or wide expression, different models were proposed to explain the small sizes of highly/widely expressed genes. We found that housekeeping genes are not more compact than narrowly expressed genes with similar expression levels, but compactness and expression level are correlated in housekeeping genes (except that highly expressed Arabidopsis HK genes have longer intron length). Meanwhile, we found evidence that genes with high functional/regulatory complexity do not have longer introns and longer proteins. The genome design hypothesis is thus not supported. Furthermore, we found that housekeeping genes are not more compact than the narrowly expressed somatic genes with similar average expression levels. Because housekeeping genes are expected to have much higher germline expression levels than narrowly expressed somatic genes, transcription-associated deletion bias is not supported. Selection of the compactness of highly expressed genes for economy is supported.

Gene function and expression level influence the insertion/fixation dynamics of distinct transposon families in mammalian introns

BACKGROUND

Transposable elements (TEs) represent more than 45% of the human and mouse genomes. Both parasitic and mutualistic features have been shown to apply to the host-TE relationship but a comprehensive scenario of the forces driving TE fixation within mammalian genes is still missing.

RESULTS

We show that intronic multispecies conserved sequences (MCSs) have been affecting TE integration frequency over time. We verify that a selective economizing pressure has been acting on TEs to decrease their frequency in highly expressed genes. After correcting for GC content, MCS density and intron size, we identified TE-enriched and TE-depleted gene categories. In addition to developmental regulators and transcription factors, TE-depleted regions encompass loci that might require subtle regulation of transcript levels or precise activation timing, such as growth factors, cytokines, hormones, and genes involved in the immune response. The latter, despite having reduced frequencies of most TE types, are significantly enriched in mammalian-wide interspersed repeats (MIRs). Analysis of orthologous genes indicated that MIR over-representation also occurs in dog and opossum immune response genes, suggesting, given the partially independent origin of MIR sequences in eutheria and metatheria, the evolutionary conservation of a specific function for MIRs located in these loci. Consistently, the core MIR sequence is over-represented in defense response genes compared to the background intronic frequency.

CONCLUSION

Our data indicate that gene function, expression level, and sequence conservation influence TE insertion/fixation in mammalian introns. Moreover, we provide the first report showing that a specific TE family is evolutionarily associated with a gene function category.