Intron size and exon evolution in Drosophila

Wnt/β-catenin signalling is required for pole-specific chromatin remodeling during planarian regeneration

For successful regeneration, the identity of the missing tissue must be specified according to the pre-existing tissue. Planarians are ideal for the study of the mechanisms underlying this process; the same field of cells can regrow a head or a tail according to the missing body part. After amputation, the differential activation of the Wnt/β-catenin signal specifies anterior versus posterior identity. Initially, both wnt1 and notum (Wnt inhibitor) are expressed in all wounds, but 48 hours later they are restricted to posterior or anterior facing wounds, respectively, by an unknown mechanism. Here we show that 12 hours after amputation, the chromatin accessibility of cells in the wound region changes according to the polarity of the pre-existing tissue in a Wnt/β-catenin-dependent manner. Genomic analyses suggest that homeobox transcription factors and chromatin-remodeling proteins are direct Wnt/β-catenin targets, which trigger the expression of posterior effectors. Finally, we identify FoxG as a wnt1 up-stream regulator, probably via binding to its first intron enhancer region.

The position of the longest intron is related to biological functions in some human genes

The evidence that introns can influence different levels of transfer of genetic information between DNA and the final product is increasing. Longer first introns were found to be a general property of eukaryotic gene structure and shown to contain a higher fraction of conserved sequence and different functional elements. Our work brings more precise information about the position of the longest introns in human protein-coding genes and possible connection with biological function and gene expression. According to our results, the position of the longest intron can be localized to the first third of introns in 64%, the second third in 19%, and the third in 17%, with notable peaks at the middle and last introns of approximately 5% and 6%, respectively. The median lengths of the longest introns decrease with increasing distance from the start of the gene from approximately 15,000 to 5,000 bp. We have shown that the position of the longest intron is in some cases linked to the biological function of the given gene. For example, DNA repair genes have the longest intron more often in the second or third. In the distribution of gene expression according to the position of the longest intron, tissue-specific profiles can be traced with the highest expression usually at the absolute positions of intron 1 and 2. In this work, we present arguments supporting the hypothesis that the position of the longest intron in a gene is another biological factor modulating the transmission of genetic information. The position of the longest intron is related to biological functions in some human genes.

A novel approach to analyze the association characteristics between post-spliced introns and their corresponding mRNA

Studies have shown that post-spliced introns promote cell survival when nutrients are scarce, and intron loss/gain can influence many stages of mRNA metabolism. However, few approaches are currently available to study the correlation between intron sequences and their corresponding mature mRNA sequences. Here, based on the results of the improved Smith-Waterman local alignment-based algorithm method (SW method) and binding free energy weighted local alignment algorithm method (BFE method), the optimal matched segments between introns and their corresponding mature mRNAs in Caenorhabditis elegans (C.elegans) and their relative matching frequency (RF) distributions were obtained. The results showed that although the distributions of relative matching frequencies on mRNAs obtained by the BFE method were similar to the SW method, the interaction intensity in 5'and 3'untranslated regions (UTRs) regions was weaker than the SW method. The RF distributions in the exon-exon junction regions were comparable, the effects of long and short introns on mRNA and on the five functional sites with BFE method were similar to the SW method. However, the interaction intensity in 5'and 3'UTR regions with BFE method was weaker than with SW method. Although the matching rate and length distribution shape of the optimal matched fragment were consistent with the SW method, an increase in length was observed. The matching rates and the length of the optimal matched fragments were mainly in the range of 60%-80% and 20-30bp, respectively. Although we found that there were still matching preferences in the 5'and 3'UTR regions of the mRNAs with BFE, the matching intensities were significantly lower than the matching intensities between introns and their corresponding mRNAs with SW method. Overall, our findings suggest that the interaction between introns and mRNAs results from synergism among different types of sequences during the evolutionary process.

STREAMING-tag system reveals spatiotemporal relationships between transcriptional regulatory factors and transcriptional activity

Transcription is a dynamic process. To detect the dynamic relationship among protein clusters of RNA polymerase II and coactivators, gene loci, and transcriptional activity, we insert an MS2 repeat, a TetO repeat, and inteins with a selection marker just downstream of the transcription start site. By optimizing the individual elements, we develop the Spliced TetO REpeAt, MS2 repeat, and INtein sandwiched reporter Gene tag (STREAMING-tag) system. Clusters of RNA polymerase II and BRD4 are observed proximal to the transcription start site of Nanog when the gene is transcribed in mouse embryonic stem cells. In contrast, clusters of MED19 and MED22 tend to be located near the transcription start site, even without transcription activity. Thus, the STREAMING-tag system reveals the spatiotemporal relationships between transcriptional activity and protein clusters near the gene. This powerful tool is useful for quantitatively understanding transcriptional regulation in living cells.

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).

Recombinant Spider Silk: Promises and Bottlenecks

Spider silk threads have exceptional mechanical properties such as toughness, elasticity and low density, which reach maximum values compared to other fibre materials. They are superior even compared to Kevlar and steel. These extraordinary properties stem from long length and specific protein structures. Spider silk proteins can consist of more than 20,000 amino acids. Polypeptide stretches account for more than 90% of the whole protein, and these domains can be repeated more than a hundred times. Each repeat unit has a specific function resulting in the final properties of the silk. These properties make them attractive for innovative material development for medical or technical products as well as cosmetics. However, with livestock breeding of spiders it is not possible to reach high volumes of silk due to the cannibalistic behaviour of these animals. In order to obtain spider silk proteins (spidroins) on a large scale, recombinant production is attempted in various expression systems such as plants, bacteria, yeasts, insects, silkworms, mammalian cells and animals. For viable large-scale production, cost-effective and efficient production systems are needed. This review describes the different types of spider silk, their proteins and structures and discusses the production of these difficult-to-express proteins in different host organisms with an emphasis on plant systems.

Using GC Content to Compare Recombination Patterns on the Sex Chromosomes and Autosomes of the Guppy, Poecilia reticulata, and Its Close Outgroup Species

Genetic and physical mapping of the guppy (Poecilia reticulata) have shown that recombination patterns differ greatly between males and females. Crossover events occur evenly across the chromosomes in females, but in male meiosis they are restricted to the tip furthest from the centromere of each chromosome, creating very high recombination rates per megabase, as in pseudoautosomal regions of mammalian sex chromosomes. We used GC content to indirectly infer recombination patterns on guppy chromosomes, based on evidence that recombination is associated with GC-biased gene conversion, so that genome regions with high recombination rates should be detectable by high GC content. We used intron sequences and third positions of codons to make comparisons between sequences that are matched, as far as possible, and are all probably under weak selection. Almost all guppy chromosomes, including the sex chromosome (LG12), have very high GC values near their assembly ends, suggesting high recombination rates due to strong crossover localization in male meiosis. Our test does not suggest that the guppy XY pair has stronger crossover localization than the autosomes, or than the homologous chromosome in the close relative, the platyfish (Xiphophorus maculatus). We therefore conclude that the guppy XY pair has not recently undergone an evolutionary change to a different recombination pattern, or reduced its crossover rate, but that the guppy evolved Y-linkage due to acquiring a male-determining factor that also conferred the male crossover pattern. We also identify the centromere ends of guppy chromosomes, which were not determined in the genome assembly.

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.

Genome-wide identification and transcriptional modulation of histone variants and modification related genes in the low pH-exposed marine rotifer Brachionus koreanus

Histone modification is considered to be a major epigenetic control mechanism. These modifications (e.g. acetylation, phosphorylation, and methylation) may affect the interaction of histones with DNA and/or regulate DNA-based processes (e.g., recombination, repair, replication, and transcription) and chromatin remodeling complexes. Despite their significance in metazoan life and evolution, few studies have been conducted to identify genes undergoing epigenetic control modification in aquatic invertebrates. In this study, we identified whole core histones (70 total genes) and post-translational modification (PTM) histone genes (63 total genes) in the marine rotifer Brachionus koreanus through whole-genome analysis, and annotated them according to the human nomenclature. Notably, upon comparative analysis of cis-regulatory motif sequences, we found that B. koreanus core histone protein structures were similar to those of mammals. Furthermore, to examine the effect of parental low pH stress on the offspring's epigenetic regulation, we investigated the expression of PTM genes in two generations of B. koreanus exposed to low pH conditions. Given that the B. koreanus genome does not possess DNA methyltransferase 1 and 3 genes, we concluded that histone genes could be involved as an important epigenetic mechanism in B. koreanus. Therefore, the histone-associated genes identified in this study could be useful for ecotoxicological studies and facilitate the application of chromatin immunoprecipitation sequencing using high-throughput DNA sequencing based on the genome-wide identification of transcription factor binding sites in rotifers.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Genome-Wide Identification, Gene Structure and Expression Analysis of the MADS-Box Gene Family Indicate Their Function in the Development of Tobacco (Nicotiana tabacum L.)

MADS-box genes play a pivotal role in various processes, including floral and seed development, controlling flowering time, regulation of fruits ripening, and respond to abiotic and biotic stressors in planta. Tobacco (Nicotiana tabacum) has been widely used as a model plant for analyzing the gene function, however, there has been less information on the regulation of flowering, and the associated genes. In the present study, a total of 168 NtMADS-box genes were identified from tobacco, and their phylogenetic relationship, chromosome locations, and gene structures were further analyzed. NtMADS-box genes can be clustered into four sub-families of Mα, Mγ, MIKC*, and MIKCC. A total of 111 NtMADS-box genes were distributed on 20 chromosomes, and 57 NtMADS-box genes were located on the unanchored scaffolds due to the complex and incomplete assembly of the tobacco genome. Expression profiles of NtMADS-box genes by microarray from 23 different tissues indicated that members in different NtMADS-box gene subfamilies might play specific roles in the growth and flower development, and the transcript levels of 24 NtMADS-box genes were confirmed by quantitative real-time PCR. Importantly, overexpressed NtSOC1/NtMADS133 could promote early flowering and dwarfism in transgenic tobacco plants. Therefore, our findings provide insights on the characterization of NtMADS-box genes to further study their functions in plant development.

Adaptation and Conservation throughout the Drosophila melanogaster Life-Cycle

Previous studies of the evolution of genes expressed at different life-cycle stages of Drosophila melanogaster have not been able to disentangle adaptive from nonadaptive substitutions when using nonsynonymous sites. Here, we overcome this limitation by combining whole-genome polymorphism data from D. melanogaster and divergence data between D. melanogaster and Drosophila yakuba. For the set of genes expressed at different life-cycle stages of D. melanogaster, as reported in modENCODE, we estimate the ratio of substitutions relative to polymorphism between nonsynonymous and synonymous sites (α) and then α is discomposed into the ratio of adaptive (ω_a) and nonadaptive (ω_na) substitutions to synonymous substitutions.

We find that the genes expressed in mid- and late-embryonic development are the most conserved, whereas those expressed in early development and postembryonic stages are the least conserved. Importantly, we found that low conservation in early development is due to high rates of nonadaptive substitutions (high ω_na), whereas in postembryonic stages it is due, instead, to high rates of adaptive substitutions (high ω_a).

By using estimates of different genomic features (codon bias, average intron length, exon number, recombination rate, among others), we also find that genes expressed in mid- and late-embryonic development show the most complex architecture: they are larger, have more exons, more transcripts, and longer introns. In addition, these genes are broadly expressed among all stages. We suggest that all these genomic features are related to the conservation of mid- and late-embryonic development. Globally, our study supports the hourglass pattern of conservation and adaptation over the life-cycle.

Small Segmental Duplications in Drosophila-High Rate of Emergence and Elimination

Segmental duplications are an important class of mutations. Because a large proportion of segmental duplications may often be strongly deleterious, high frequency or fixed segmental duplications may represent only a tiny fraction of the mutational input. To understand the emergence and elimination of segmental duplications, we survey polymorphic duplications, including tandem and interspersed duplications, in natural populations of Drosophila by haploid embryo genomes. As haploid embryos are not expected to be heterozygous, the genome, sites of heterozygosity (referred to as pseudoheterozygous sites [PHS]), may likely represent recent duplications that have acquired new mutations. Among the 29 genomes of Drosophila melanogaster, we identify 2,282 polymorphic PHS duplications (linked PHS regions) in total or 154 PHS duplications per genome. Most PHS duplications are small (83.4% < 500 bp), Drosophila melanogaster lineage specific, and strain specific (72.6% singletons). The excess of the observed singleton PHS duplications deviates significantly from the neutral expectation, suggesting that most PHS duplications are strongly deleterious. In addition, these small segmental duplications are not evenly distributed in genomic regions and less common in noncoding functional element regions. The underrepresentation in RNA polymerase II binding sites and regions with active histone modifications is correlated with ages of duplications. In conclusion, small segmental duplications occur frequently in Drosophila but rapidly eliminated by natural selection.

Evolutionary and functional analyses of the 2-oxoglutarate-dependent dioxygenase genes involved in the flavonoid biosynthesis pathway in tobacco

MAIN CONCLUSION

This study illustrates the differences in the gene structure of 2-oxoglutarate-dependent oxygenase involved in flavonoid biosynthesis (2ODD-IFB), and their potential roles in regulating tobacco flavonoid biosynthesis and plant growth. Flavonol synthase (FLS), anthocyanidin synthase (ANS), and flavanone 3β-hydroxylase belong to the 2-oxoglutarate-dependent (2ODD) oxygenase family, and each performs crucial functions in the biosynthesis of flavonoids. We identified two NtFLS genes, two NtANS genes, and four NtF3H genes from Nicotiana tabacum genome, as well as their homologous genes in the N. sylvestris and N. tomentosiformis genomes. Our phylogenetic analysis indicated that these three types of genes split from each other before the divergence of gymnosperms and angiosperms. FLS evolved faster in the eudicot plants, whereas ANS evolved faster in the monocot plants. Gene structure analysis revealed two fragment insertions occurred at different times in the intron one position of tobacco FLS genes. Homologous protein modeling revealed distinct structures in the N terminus of the tobacco 2ODD oxygenases. We found that the expression patterns of genes encoding tobacco 2ODD oxygenases in flavonoids biosynthesis (2ODD-IFB) did not determine the accumulation patterns of flavonoids among various tobacco tissues, but strongly affected the concentration of flavonoids in the tissues, where they were biosynthesized. More carbon resource flowed to the flavonol biosynthesis when NtANS gene was silenced, otherwise more anthocyanidin accumulated when NtFLS gene was repressed. This study illustrates the 2ODD-IFB gene structure evolution, differences among their protein structures, and provides a foundation for regulating plant development and altering flavonoid content and/or composition through the manipulation of plant 2ODD-IFB genes.

SNP-based susceptibility-resistance association and mRNA expression regulation analyses of tlr7 to grass carp Ctenopharyngodon idella reovirus

Eleven single nucleotide polymorphisms (SNP) in Ctenopharyngodon idella toll-like receptor 7 (citlr7) gene, containing two in the 5'-flanking region, three within the single intron and six distributed in the coding sequence (CDS), were identified. A case-control study of 73 susceptible individuals and 67 resistant individuals was conducted to test the SNPs-based susceptibility-resistance association and mRNA expression of citlr7 to grass carp reovirus (GCRV), showing that both 820 A/G and 1726 A/G were significantly correlative sites in genotype (P < 0·05). Multifactor dimensionality reduction (MDR) analysis suggested the exertion of antiviral effects of 820 A/G might rely on SNPs interactions of citlr7 and C. idella toll-like receptor 8 (citlr8). Combining the mortality rate and citlr7 mRNA expression, it was suggested that 1726 GG-genotyped individuals might be more resistant than 1726 A/G genotyped individuals, indicating the selection on synonymous mutations in 1726 A/G might be susceptibility-resistance-type specific. In addition, haplotype analysis uncovered no significantly correlative haplotypes in citlr7. These findings may provide an in-depth insight for the further functional research of citlr7. The potential genetic markers identified may contribute to the molecular and transgenic breeding of C. idella.

Inferring the determinants of protein evolutionary rates in mammals

Understanding the determinants of protein evolutionary rates is one of the most fundamental evolutionary questions. Previous studies have revealed that many biological variables are tightly associated with protein evolutionary rates in mammals. However, the dominant role of these biological variables and their combinatorial effects to evolutionary rates of mammalian proteins are still less understood. In this work, we derived a quantitative model to correlate protein evolutionary rates with the levels of these variables. The result showed that only a small number of variables are necessary to accurately predict protein evolutionary rates, among which miRNA regulation plays the most important role. Our result suggested that biological variables are extensively interrelated and suffer from hidden redundancies in determining protein evolutionary rates. Various variables should be considered in a natural ensemble to comprehensively assess the determinants of protein evolutionary rate.

ERK1 and ERK2 Map Kinases: Specific Roles or Functional Redundancy?

The MAP kinase signaling cascade Ras/Raf/MEK/ERK has been involved in a large variety of cellular and physiological processes that are crucial for life. Many pathological situations have been associated to this pathway. More than one isoform has been described at each level of the cascade. In this review we devoted our attention to ERK1 and ERK2, which are the effector kinases of the pathway. Whether ERK1 and ERK2 specify functional differences or are in contrast functionally redundant, constitutes an ongoing debate despite the huge amount of studies performed to date. In this review we compiled data on ERK1 vs. ERK2 gene structures, protein sequences, expression levels, structural and molecular mechanisms of activation and substrate recognition. We have also attempted to perform a rigorous analysis of studies regarding the individual roles of ERK1 and ERK2 by the means of morpholinos, siRNA, and shRNA silencing as well as gene disruption or gene replacement in mice. Finally, we comment on a recent study of gene and protein evolution of ERK isoforms as a distinct approach to address the same question. Our review permits the evaluation of the relevance of published studies in the field especially when measurements of global ERK activation are taken into account. Our analysis favors the hypothesis of ERK1 and ERK2 exhibiting functional redundancy and points to the concept of the global ERK quantity, and not isoform specificity, as being the essential determinant to achieve ERK function.

Intron retention and rhythmic diel pattern regulation of carotenoid cleavage dioxygenase 2 during crocetin biosynthesis in saffron

The carotenoid cleavage dioxygenase 2, a new member of the CCD family, catalyzes the conversion of zeaxanthin into crocetin-dialdehyde in Crocus. CCD2 is expressed in flowers, being responsible for the yellow, orange and red colorations displayed by tepals and stigma. Three CsCCD2 genes were identified in Crocus sativus, the longest contains ten exons and the shorter is a truncated copy with no introns and which lacks one exon sequence. Analysis of RNA-seq datasets of three developmental stages of saffron stigma allowed the determination of alternative splicing in CsCCD2, being intron retention (IR) the prevalent form of alternative splicing in CsCCD2. Further, high IR was observed in tissues that do not accumulate crocetin. The analysis of one CsCCD2 promoter showed cis-regulatory motifs involved in the response to light, temperature, and circadian regulation. The light and circadian regulation are common elements shared with the previously characterized CsLycB2a promoter, and these shared common cis-acting elements may represent binding sites for transcription factors responsible for co-regulation of these genes during the development of the stigma in saffron. A daily coordinated rhythmic regulation for CsCCD2 and CsLycB2a was observed, with higher levels of mRNA occurring at low temperatures during darkness, confirming the results obtained in the in silico promoter analysis. In addition, to the light and temperature dependent regulation of CsCCD2 expression, the apocarotenoid β-cyclocitral up-regulated CsCCD2 expression and could acts as a mediator of chromoplast-to-nucleus signalling, coordinating the expression of CsCCD2 with the developmental state of the chromoplast in the developing stigma.

Comparative genomics of flatworms (platyhelminthes) reveals shared genomic features of ecto- and endoparastic neodermata

The ectoparasitic Monogenea comprise a major part of the obligate parasitic flatworm diversity. Although genomic adaptations to parasitism have been studied in the endoparasitic tapeworms (Cestoda) and flukes (Trematoda), no representative of the Monogenea has been investigated yet. We present the high-quality draft genome of Gyrodactylus salaris, an economically important monogenean ectoparasite of wild Atlantic salmon (Salmo salar). A total of 15,488 gene models were identified, of which 7,102 were functionally annotated. The controversial phylogenetic relationships within the obligate parasitic Neodermata were resolved in a phylogenomic analysis using 1,719 gene models (alignment length of >500,000 amino acids) for a set of 16 metazoan taxa. The Monogenea were found basal to the Cestoda and Trematoda, which implies ectoparasitism being plesiomorphic within the Neodermata and strongly supports a common origin of complex life cycles. Comparative analysis of seven parasitic flatworm genomes identified shared genomic features for the ecto- and endoparasitic lineages, such as a substantial reduction of the core bilaterian gene complement, including the homeodomain-containing genes, and a loss of the piwi and vasa genes, which are considered essential for animal development. Furthermore, the shared loss of functional fatty acid biosynthesis pathways and the absence of peroxisomes, the latter organelles presumed ubiquitous in eukaryotes except for parasitic protozoans, were inferred. The draft genome of G. salaris opens for future in-depth analyses of pathogenicity and host specificity of poorly characterized G. salaris strains, and will enhance studies addressing the genomics of host–parasite interactions and speciation in the highly diverse monogenean flatworms.

Introns: The Functional Benefits of Introns in Genomes

The intron has been a big biological mystery since it was first discovered in several aspects. First, all of the completely sequenced eukaryotes harbor introns in the genomic structure, whereas no prokaryotes identified so far carry introns. Second, the amount of total introns varies in different species. Third, the length and number of introns vary in different genes, even within the same species genome. Fourth, all introns are copied into RNAs by transcription and DNAs by replication processes, but intron sequences do not participate in protein-coding sequences. The existence of introns in the genome should be a burden to some cells, because cells have to consume a great deal of energy to copy and excise them exactly at the correct positions with the help of complicated spliceosomal machineries. The existence throughout the long evolutionary history is explained, only if selective advantages of carrying introns are assumed to be given to cells to overcome the negative effect of introns. In that regard, we summarize previous research about the functional roles or benefits of introns. Additionally, several other studies strongly suggesting that introns should not be junk will be introduced.

Identification and expression analysis of a novel stylicin antimicrobial peptide from Kuruma shrimp (Marsupenaeus japonicus)

Antimicrobial peptides (AMPs) are important components of the innate immune system and function as the first line of defense against invading pathogens. In current study we identified, cloned and characterized a novel stylicin AMP from Kuruma shrimp Marsupenaeus japonicus (Mj-sty). The full-length cDNA of Mj-sty was 428 bp with an open reading frame of 315 bp that encoded 104 amino acids. The theoretical molecular mass of mature Mj-sty was 8.693 kDa with an isoelectric point (pI) of 4.79. A proline-rich N-terminal region and a C-terminal region contained 13 cysteine residues were identified. Genomic sequence analysis with respect to its cDNA showed that Mj-sty was organized into two exons interrupted by one intron. Tissue-specific expression revealed that Mj-sty was mainly transcribed in gills and hemocytes. Expression of Mj-sty in early developmental stages demonstrated that Mj-sty mRNA were present from fertilized eggs to post-larvae of 17 days (PL17), and the expression levels showed a significant variation in different developmental stages. After challenge of white spot syndrome virus (WSSV), the time-dependent expression pattern of Mj-sty in both gills and hepatopancrease showed down-regulation at the early hours of infection, subsequently up-regulation and down-regulation, and then up-regulation at the end hours to almost the half of the controls. The results indicate that Mj-sty is potentially involved in the ontogenesis and immune responses against WSSV.

ERK1 and ERK2 present functional redundancy in tetrapods despite higher evolution rate of ERK1

Background

The Ras/Raf/MEK/ERK signaling pathway is involved in essential cell processes and it is abnormally activated in ~30 % of cancers and cognitive disorders. Two ERK isoforms have been described, ERK1 and ERK2; ERK2 being regarded by many as essential due to the embryonic lethality of ERK2 knock-out mice, whereas mice lacking ERK1 are viable and fertile. The controversial question of why we have two ERKs and whether they have differential functions or display functional redundancy has not yet been resolved.

Results

To investigate this question we used a novel approach based on comparing the evolution of ERK isoforms’ sequences and protein expression across vertebrates. We gathered and cloned erk1 and erk2 coding sequences and we examined protein expression of isoforms in brain extracts in all major clades of vertebrate evolution. For the first time, we measured each isoforms’ relative protein level in phylogenetically distant animals using anti-phospho antibodies targeting active ERKs. We demonstrate that squamates (lizards, snakes and geckos), despite having both genes, do not express ERK2 protein whereas other tetrapods either do not express ERK1 protein or have lost the erk1 gene. To demonstrate the unexpected squamates’ lack of ERK2 expression, we targeted each ERK isoform in lizard primary fibroblasts by specific siRNA-mediated knockdown. We also found that undetectable expression of ERK2 in lizard is compensated by a greater strength of lizard’s erk1 promoter. Finally, phylogenetic analysis revealed that ERK1 amino acids sequences evolve faster than ERK2’s likely due to genomic factors, including a large difference in gene size, rather than from functional differences since amino acids essential for function are kept invariant.

Conclusions

ERK isoforms appeared by a single gene duplication at the onset of vertebrate evolution at least 400 Mya. Our results demonstrate that tetrapods can live by expressing either one or both ERK isoforms, supporting the notion that ERK1/2 act interchangeably. Substrate recognition sites and catalytic cleft are nearly invariant in all vertebrate ERKs further suggesting functional redundancy. We suggest that future ERK research should shift towards understanding the role and regulation of total ERK quantity, especially in light of newly described erk2 gene amplification identified in tumors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0450-x) contains supplementary material, which is available to authorized users.

Highly variable recombinational landscape modulates efficacy of natural selection in birds

Determining the rate of protein evolution and identifying the causes of its variation across the genome are powerful ways to understand forces that are important for genome evolution. By using a multitissue transcriptome data set from great tit (Parus major), we analyzed patterns of molecular evolution between two passerine birds, great tit and zebra finch (Taeniopygia guttata), using the chicken genome (Gallus gallus) as an outgroup. We investigated whether a special feature of avian genomes, the highly variable recombinational landscape, modulates the efficacy of natural selection through the effects of Hill-Robertson interference, which predicts that selection should be more effective in removing deleterious mutations and incorporating beneficial mutations in high-recombination regions than in low-recombination regions. In agreement with these predictions, genes located in low-recombination regions tend to have a high proportion of neutrally evolving sites and relaxed selective constraint on sites subject to purifying selection, whereas genes that show strong support for past episodes of positive selection appear disproportionally in high-recombination regions. There is also evidence that genes located in high-recombination regions tend to have higher gene expression specificity than those located in low-recombination regions. Furthermore, more compact genes (i.e., those with fewer/shorter introns or shorter proteins) evolve faster than less compact ones. In sum, our results demonstrate that transcriptome sequencing is a powerful method to answer fundamental questions about genome evolution in nonmodel organisms.

Genome-wide analysis of alternative splicing in Volvox carteri

Background

Alternative splicing is an essential mechanism for increasing transcriptome and proteome diversity in eukaryotes. Particularly in multicellular eukaryotes, this mechanism is involved in the regulation of developmental and physiological processes like growth, differentiation and signal transduction.

Results

Here we report the genome-wide analysis of alternative splicing in the multicellular green alga Volvox carteri. The bioinformatic analysis of 132,038 expressed sequence tags (ESTs) identified 580 alternative splicing events in a total of 426 genes. The predominant type of alternative splicing in Volvox is intron retention (46.5%) followed by alternative 5′ (17.9%) and 3′ (21.9%) splice sites and exon skipping (9.5%). Our analysis shows that in Volvox at least ~2.9% of the intron-containing genes are subject to alternative splicing. Considering the total number of sequenced ESTs, the Volvox genome seems to provide more favorable conditions (e.g., regarding length and GC content of introns) for the occurrence of alternative splicing than the genome of its close unicellular relative Chlamydomonas. Moreover, many randomly chosen alternatively spliced genes of Volvox do not show alternative splicing in Chlamydomonas. Since the Volvox genome contains about the same number of protein-coding genes as the Chlamydomonas genome (~14,500 protein-coding genes), we assumed that alternative splicing may play a key role in generation of genomic diversity, which is required to evolve from a simple one-cell ancestor to a multicellular organism with differentiated cell types (Mol Biol Evol 31:1402-1413, 2014). To confirm the alternative splicing events identified by bioinformatic analysis, several genes with different types of alternatively splicing have been selected followed by experimental verification of the predicted splice variants by RT-PCR.

Conclusions

The results show that our approach for prediction of alternative splicing events in Volvox was accurate and reliable. Moreover, quantitative real-time RT-PCR appears to be useful in Volvox for analyses of relationships between the appearance of specific alternative splicing variants and different kinds of physiological, metabolic and developmental processes as well as responses to environmental changes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1117) contains supplementary material, which is available to authorized users.

Intron evolution in Neurospora: the role of mutational bias and selection

We used comparative and population genomics to study intron evolutionary dynamics in the fungal model genus Neurospora. For our investigation, we used well-annotated genomes of N. crassa, N. discreta, and N. tetrasperma, and 92 resequenced genomes of N. tetrasperma from natural populations. By analyzing the four well-annotated genomes, we identified 9495 intron sites in 7619 orthologous genes. Our data supports nonhomologous end joining (NHEJ) and tandem duplication as mechanisms for intron gains in the genus and the RT-mRNA process as a mechanism for intron loss. We found a moderate intron gain rate (5.78–6.89 × 10⁻¹³ intron gains per nucleotide site per year) and a high intron loss rate (7.53–13.76 × 10⁻¹⁰ intron losses per intron sites per year) as compared to other eukaryotes. The derived intron gains and losses are skewed to high frequencies, relative to neutral SNPs, in natural populations of N. tetrasperma, suggesting that selection is involved in maintaining a high intron turnover. Furthermore, our analyses of the association between intron population-level frequency and genomic features suggest that selection is involved in shaping a 5′ intron position bias and a low intron GC content. However, intron sequence analyses suggest that the gained introns were not exposed to recent selective sweeps. Taken together, this work contributes to our understanding of the importance of mutational bias and selection in shaping the intron distribution in eukaryotic genomes.

Gene expression variation in Drosophila melanogaster due to rare transposable element insertion alleles of large effect

Transposable elements are a common source of genetic variation that may play a substantial role in contributing to gene expression variation. However, the contribution of transposable elements to expression variation thus far consists of a handful of examples. We used previously published gene expression data from 37 inbred Drosophila melanogaster lines from the Drosophila Genetic Reference Panel to perform a genome-wide assessment of the effects of transposable elements on gene expression. We found thousands of transcripts with transposable element insertions in or near the transcript and that the presence of a transposable element in or near a transcript is significantly associated with reductions in expression. We estimate that within this example population, ∼2.2% of transcripts have a transposable element insertion, which significantly reduces expression in the line containing the transposable element. We also find that transcripts with insertions within 500 bp of the transcript show on average a 0.67 standard deviation decrease in expression level. These large decreases in expression level are most pronounced for transposable element insertions close to transcripts and the effect diminishes for more distant insertions. This work represents the first genome-wide analysis of gene expression variation due to transposable elements and suggests that transposable elements are an important class of mutation underlying expression variation in Drosophila and likely in other systems, given the ubiquity of these mobile elements in eukaryotic genomes.

Menzerath-Altmann law in mammalian exons reflects the dynamics of gene structure evolution

Genomic sequences exhibit self-organization properties at various hierarchical levels. One such is the gene structure of higher eukaryotes with its complex exon/intron arrangement. Exon sizes and exon numbers in genes have been shown to conform to a law derived from statistical linguistics and formulated by Menzerath and Altmann, according to which the mean size of the constituents of an entity is inversely related to the number of these constituents. We herein perform a detailed analysis of this property in the complete exon set of the mouse genome in correlation to the sequence conservation of each exon and the transcriptional complexity of each gene locus. We show that extensive linear fits, representative of accordance to Menzerath-Altmann law are restricted to a particular subset of genes that are formed by exons under low or intermediate sequence constraints and have a small number of alternative transcripts. Based on this observation we propose a hypothesis for the law of Menzerath-Altmann in mammalian genes being predominantly due to genes that are more versatile in function and thus, more prone to undergo changes in their structure. To this end we demonstrate one test case where gene categories of different functionality also show differences in the extent of conformity to Menzerath-Altmann law.

Conservation in first introns is positively associated with the number of exons within genes and the presence of regulatory epigenetic signals

Background

Genomes of higher eukaryotes have surprisingly long first introns and in some cases, the first introns have been shown to have higher conservation relative to other introns. However, the functional relevance of conserved regions in the first introns is poorly understood. Leveraging the recent ENCODE data, here we assess potential regulatory roles of conserved regions in the first intron of human genes.

Results

We first show that relative to other downstream introns, the first introns are enriched for blocks of highly conserved sequences. We also found that the first introns are enriched for several chromatin marks indicative of active regulatory regions and this enrichment of regulatory marks is correlated with enrichment of conserved blocks in the first intron; the enrichments of conservation and regulatory marks in first intron are not entirely explained by a general, albeit variable, bias for certain marks toward the 5’ end of introns. Interestingly, conservation as well as proportions of active regulatory chromatin marks in the first intron of a gene correlates positively with the numbers of exons in the gene but the correlation is significantly weakened in second introns and negligible beyond the second intron. The first intron conservation is also positively correlated with the gene’s expression level in several human tissues. Finally, a gene-wise analysis shows significant enrichments of active chromatin marks in conserved regions of first introns, relative to the conserved regions in other introns of the same gene.

Conclusions

Taken together, our analyses strongly suggest that first introns are enriched for active transcriptional regulatory signals under purifying selection.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-526) contains supplementary material, which is available to authorized users.

Impacts of pretranscriptional DNA methylation, transcriptional transcription factor, and posttranscriptional microRNA regulations on protein evolutionary rate

Gene expression is largely regulated by DNA methylation, transcription factor (TF), and microRNA (miRNA) before, during, and after transcription, respectively. Although the evolutionary effects of TF/miRNA regulations have been widely studied, evolutionary analysis of simultaneously accounting for DNA methylation, TF, and miRNA regulations and whether promoter methylation and gene body (coding regions) methylation have different effects on the rate of gene evolution remain uninvestigated. Here, we compared human–macaque and human–mouse protein evolutionary rates against experimentally determined single base-resolution DNA methylation data, revealing that promoter methylation level is positively correlated with protein evolutionary rates but negatively correlated with TF/miRNA regulations, whereas the opposite was observed for gene body methylation level. Our results showed that the relative importance of these regulatory factors in determining the rate of mammalian protein evolution is as follows: Promoter methylation ≈ miRNA regulation > gene body methylation > TF regulation, and further indicated that promoter methylation and miRNA regulation have a significant dependent effect on protein evolutionary rates. Although the mechanisms underlying cooperation between DNA methylation and TFs/miRNAs in gene regulation remain unclear, our study helps to not only illuminate the impact of these regulatory factors on mammalian protein evolution but also their intricate interaction within gene regulatory networks.

Hybridization reveals the evolving genomic architecture of speciation

The rate at which genomes diverge during speciation is unknown, as are the physical dynamics of the process. Here, we compare full genome sequences of 32 butterflies, representing five species from a hybridizing Heliconius butterfly community, to examine genome-wide patterns of introgression and infer how divergence evolves during the speciation process. Our analyses reveal that initial divergence is restricted to a small fraction of the genome, largely clustered around known wing-patterning genes. Over time, divergence evolves rapidly, due primarily to the origin of new divergent regions. Furthermore, divergent genomic regions display signatures of both selection and adaptive introgression, demonstrating the link between microevolutionary processes acting within species and the origin of species across macroevolutionary timescales. Our results provide a uniquely comprehensive portrait of the evolving species boundary due to the role that hybridization plays in reducing the background accumulation of divergence at neutral sites.

Aberrantly spliced HTT, a new player in Huntington's disease pathogenesis

Huntington's disease (HD) is an adult-onset neurodegenerative disorder caused by a mutated CAG repeat in the huntingtin gene that is translated into an expanded polyglutamine tract. The clinical manifestation of HD is a progressive physical, cognitive, and psychiatric deterioration that is eventually fatal. The mutant huntingtin protein is processed into several smaller fragments, which have been implicated as critical factors in HD pathogenesis. The search for proteases responsible for their production has led to the identification of several cleavage sites on the huntingtin protein. However, the origin of the small N-terminal fragments that are found in HD postmortem brains has remained elusive. Recent mapping of huntingtin fragments in a mouse model demonstrated that the smallest N-terminal fragment is an exon 1 protein. This discovery spurred our hypothesis that mis-splicing as opposed to proteolysis could be generating the smallest huntingtin fragment. We demonstrated that mis-splicing of mutant huntingtin intron 1 does indeed occur and results in a short polyadenylated mRNA, which is translated into an exon 1 protein. The exon 1 protein fragment is highly pathogenic. Transgenic mouse models containing just human huntingtin exon 1 develop a rapid onset of HD-like symptoms. Our finding that a small, mis-spliced HTT transcript and corresponding exon 1 protein are produced in the context of an expanded CAG repeat has unraveled a new molecular mechanism in HD pathogenesis. Here we present detailed models of how mis-splicing could be facilitated, what challenges remain in this model, and implications for therapeutic studies.

The impact of trans-regulation on the evolutionary rates of metazoan proteins

Transcription factor (TF) and microRNA (miRNA) are two crucial trans-regulatory factors that coordinately control gene expression. Understanding the impacts of these two factors on the rate of protein sequence evolution is of great importance in evolutionary biology. While many biological factors associated with evolutionary rate variations have been studied, evolutionary analysis of simultaneously accounting for TF and miRNA regulations across metazoans is still uninvestigated. Here, we provide a series of statistical analyses to assess the influences of TF and miRNA regulations on evolutionary rates across metazoans (human, mouse and fruit fly). Our results reveal that the negative correlations between trans-regulation and evolutionary rates hold well across metazoans, but the strength of TF regulation as a rate indicator becomes weak when the other confounding factors that may affect evolutionary rates are controlled. We show that miRNA regulation tends to be a more essential indicator of evolutionary rates than TF regulation, and the combination of TF and miRNA regulations has a significant dependent effect on protein evolutionary rates. We also show that trans-regulation (especially miRNA regulation) is much more important in human/mouse than in fruit fly in determining protein evolutionary rates, suggesting a considerable variation in rate determinants between vertebrates and invertebrates.

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.

Three independent determinants of protein evolutionary rate

One of the most widely accepted ideas related to the evolutionary rates of proteins is that functionally important residues or regions evolve slower than other regions, a reasonable outcome of which should be a slower evolutionary rate of the proteins with a higher density of functionally important sites. Oddly, the role of functional importance, mainly measured by essentiality, in determining evolutionary rate has been challenged in recent studies. Several variables other than protein essentiality, such as expression level, gene compactness, protein-protein interactions, etc., have been suggested to affect protein evolutionary rate. In the present review, we try to refine the concept of functional importance of a gene, and consider three factors-functional importance, expression level, and gene compactness, as independent determinants of evolutionary rate of a protein, based not only on their known correlation with evolutionary rate but also on a reasonable mechanistic model. We suggest a framework based on these mechanistic models to correctly interpret the correlations between evolutionary rates and the various variables as well as the interrelationships among the variables.

Full-length minor ampullate spidroin gene sequence

Spider silk includes seven protein based fibers and glue-like substances produced by glands in the spider's abdomen. Minor ampullate silk is used to make the auxiliary spiral of the orb-web and also for wrapping prey, has a high tensile strength and does not supercontract in water. So far, only partial cDNA sequences have been obtained for minor ampullate spidroins (MiSps). Here we describe the first MiSp full-length gene sequence from the spider species Araneus ventricosus, using a multidimensional PCR approach. Comparative analysis of the sequence reveals regulatory elements, as well as unique spidroin gene and protein architecture including the presence of an unusually large intron. The spliced full-length transcript of MiSp gene is 5440 bp in size and encodes 1766 amino acid residues organized into conserved nonrepetitive N- and C-terminal domains and a central predominantly repetitive region composed of four units that are iterated in a non regular manner. The repeats are more conserved within A. ventricosus MiSp than compared to repeats from homologous proteins, and are interrupted by two nonrepetitive spacer regions, which have 100% identity even at the nucleotide level.

Structural and functional characterization of ribosomal protein gene introns in sponges

Ribosomal protein genes (RPGs) are a powerful tool for studying intron evolution. They exist in all three domains of life and are much conserved. Accumulating genomic data suggest that RPG introns in many organisms abound with non-protein-coding-RNAs (ncRNAs). These ancient ncRNAs are small nucleolar RNAs (snoRNAs) essential for ribosome assembly. They are also mobile genetic elements and therefore probably important in diversification and enrichment of transcriptomes through various mechanisms such as intron/exon gain/loss. snoRNAs in basal metazoans are poorly characterized. We examined 449 RPG introns, in total, from four demosponges: Amphimedon queenslandica, Suberites domuncula, Suberites ficus and Suberites pagurorum and showed that RPG introns from A. queenslandica share position conservancy and some structural similarity with "higher" metazoans. Moreover, our study indicates that mobile element insertions play an important role in the evolution of their size. In four sponges 51 snoRNAs were identified. The analysis showed discrepancies between the snoRNA pools of orthologous RPG introns between S. domuncula and A. queenslandica. Furthermore, these two sponges show as much conservancy of RPG intron positions between each other as between themselves and human. Sponges from the Suberites genus show consistency in RPG intron position conservation. However, significant differences in some of the orthologous RPG introns of closely related sponges were observed. This indicates that RPG introns are dynamic even on these shorter evolutionary time scales.

Origin and evolution of spliceosomal introns

Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded ‘introns first’ held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers’ Reports section.

Shrimp invertebrate lysozyme i-lyz: gene structure, molecular model and response of c and i lysozymes to lipopolysaccharide (LPS)

The invertebrate lysozyme (i-lyz or destabilase) is present in shrimp. This protein may have a function as a peptidoglycan-breaking enzyme and as a peptidase. Shrimp is commonly infected with Vibrio sp., a Gram-negative bacteria, and it is known that the c-lyz (similar to chicken lysozyme) is active against these bacteria. To further understand the regulation of lysozymes, we determined the gene sequence and modeled the protein structure of i-lyz. In addition, the expression of i-lyz and c-lyz in response to lipopolysaccharide (LPS) was studied. The shrimp i-lyz gene is interrupted by two introns with canonical splice junctions. The expression of the shrimp i-lyz was transiently down-regulated after LPS injection followed by induction after 6 h in hepatopancreas. In contrast, c-lyz was up-regulated in hepatopancreas 4 h post-injection and slightly down-regulated in gills. The L. vannamei i-lyz does not contain the catalytic residues for muramidase (glycohydrolase) neither isopeptidase activities; however, it is known that the antibacterial activity does not solely rely on the enzymatic activity of the protein. The study of invertebrate lysozyme will increase our understanding of the regulatory process of the defense mechanisms.

Mutation bias is the driving force of codon usage in the Gallus gallus genome

Synonymous codons are used with different frequencies both among species and among genes within the same genome and are controlled by neutral processes (such as mutation and drift) as well as by selection. Up to now, a systematic examination of the codon usage for the chicken genome has not been performed. Here, we carried out a whole genome analysis of the chicken genome by the use of the relative synonymous codon usage (RSCU) method and identified 11 putative optimal codons, all of them ending with uracil (U), which is significantly departing from the pattern observed in other eukaryotes. Optimal codons in the chicken genome are most likely the ones corresponding to highly expressed transfer RNA (tRNAs) or tRNA gene copy numbers in the cell. Codon bias, measured as the frequency of optimal codons (Fop), is negatively correlated with the G + C content, recombination rate, but positively correlated with gene expression, protein length, gene length and intron length. The positive correlation between codon bias and protein, gene and intron length is quite different from other multi-cellular organism, as this trend has been only found in unicellular organisms. Our data displayed that regional G + C content explains a large proportion of the variance of codon bias in chicken. Stepwise selection model analyses indicate that G + C content of coding sequence is the most important factor for codon bias. It appears that variation in the G + C content of CDSs accounts for over 60% of the variation of codon bias. This study suggests that both mutation bias and selection contribute to codon bias. However, mutation bias is the driving force of the codon usage in the Gallus gallus genome. Our data also provide evidence that the negative correlation between codon bias and recombination rates in G. gallus is determined mostly by recombination-dependent mutational patterns.

Comparative analysis of teleost genome sequences reveals an ancient intron size expansion in the zebrafish lineage

We have developed a bioinformatics pipeline for the comparative evolutionary analysis of Ensembl genomes and have used it to analyze the introns of the five available teleost fish genomes. We show our pipeline to be a powerful tool for revealing variation between genomes that may otherwise be overlooked with simple summary statistics. We identify that the zebrafish, Danio rerio, has an unusual distribution of intron sizes, with a greater number of larger introns in general and a notable peak in the frequency of introns of approximately 500 to 2,000 bp compared with the monotonically decreasing frequency distributions of the other fish. We determine that 47% of D. rerio introns are composed of repetitive sequences, although the remainder, over 331 Mb, is not. Because repetitive elements may be the origin of the majority of all noncoding DNA, it is likely that the remaining D. rerio intronic sequence has an ancient repetitive origin and has since accumulated so many mutations that it can no longer be recognized as such. To study such an ancient expansion of repeats in the Danio, lineage will require further comparative analysis of fish genomes incorporating a broader distribution of teleost lineages.

Contrasting 5' and 3' evolutionary histories and frequent evolutionary convergence in Meis/hth gene structures

Organisms show striking differences in genome structure; however, the functional implications and fundamental forces that govern these differences remain obscure. The intron–exon organization of nuclear genes is involved in a particularly large variety of structures and functional roles. We performed a 22-species study of Meis/hth genes, intron-rich homeodomain-containing transcription factors involved in a wide range of developmental processes. Our study revealed three surprising results that suggest important and very different functions for Meis intron–exon structures. First, we find unexpected conservation across species of intron positions and lengths along most of the Meis locus. This contrasts with the high degree of structural divergence found in genome-wide studies and may attest to conserved regulatory elements residing within these conserved introns. Second, we find very different evolutionary histories for the 5′ and 3′ regions of the gene. The 5′-most 10 exons, which encode the highly conserved Meis domain and homeodomain, show striking conservation. By contrast, the 3′ of the gene, which encodes several domains implicated in transcriptional activation and response to cell signaling, shows a remarkably active evolutionary history, with diverse isoforms and frequent creation and loss of new exons and splice sites. This region-specific diversity suggests evolutionary “tinkering,” with alternative splicing allowing for more subtle regulation of protein function. Third, we find a large number of cases of convergent evolution in the 3′ region, including 1) parallel losses of ancestral coding sequence, 2) parallel gains of external and internal splice sites, and 3) recurrent truncation of C-terminal coding regions. These results attest to the importance of locus-specific splicing functions in differences in structural evolution across genes, as well as to commonalities of forces shaping the evolution of individual genes along different lineages.

Spliceosomal intron size expansion in domesticated grapevine (Vitis vinifera)

BACKGROUND

Spliceosomal introns are important components of eukaryotic genes as their structure, sizes and contents reflect the architecture of gene and genomes. Intron size, determined by both neutral evolution, repetitive elements activities and potential functional constraints, varies significantly in eukaryotes, suggesting unique dynamics and evolution in different lineages of eukaryotic organisms. However, the evolution of intron size, is rarely studied. To investigate intron size dynamics in flowering plants, in particular domesticated grapevines, a survey of intron size and content in wine grape (Vitis vinifera Pinot Noir) genes was conducted by assembling and mapping the transcriptome of V. vinifera genes from ESTs to characterize and analyze spliceosomal introns.

RESULTS

Uncommonly large size of spliceosomal intron was observed in V. vinifera genome, otherwise inconsistent with overall genome size dynamics when comparing Arabidopsis, Populus and Vitis. In domesticated grapevine, intron size is generally not related to gene function. The composition of enlarged introns in grapevines indicated extensive transposable element (TE) activity within intronic regions. TEs comprise about 80% of the expanded intron space and in particular, recent LTR retrotransposon insertions are enriched in these intronic regions, suggesting an intron size expansion in the lineage leading to domesticated grapevine, instead of size contractions in Arabidopsis and Populus. Comparative analysis of selected intronic regions in V. vinifera cultivars and wild grapevine species revealed that accelerated TE activity was associated with grapevine domestication, and in some cases with the development of specific cultivars.

CONCLUSIONS

In this study, we showed intron size expansion driven by TE activities in domesticated grapevines, likely a result of long-term vegetative propagation and intensive human care, which simultaneously promote TE proliferation and repress TE removal mechanisms such as recombination. The intron size expansion observed in domesticated grapevines provided an example of rapid plant genome evolution in response to artificial selection and propagation, and may shed light on the important genomic changes during domestication. In addition, the transcriptome approach used to gather intron size data significantly improved annotations of the V. vinifera genome.