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

The reference genome of Macropodus opercularis (the paradise fish)

Amongst fishes, zebrafish (Danio rerio) has gained popularity as a model system over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species has a highly complex behavioral repertoire and has been the subject of many ethological investigations but lacks genomic resources. Here we report the reference genome assembly of M. opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of 483,077,705 base pairs (~483 Mb) on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ~90% of them to orthogroups.

The reference genome of the paradise fish (Macropodus opercularis)

Over the decades, a small number of model species, each representative of a larger taxa, have dominated the field of biological research. Amongst fishes, zebrafish (Danio rerio) has gained popularity over most other species and while their value as a model is well documented, their usefulness is limited in certain fields of research such as behavior. By embracing other, less conventional experimental organisms, opportunities arise to gain broader insights into evolution and development, as well as studying behavioral aspects not available in current popular model systems. The anabantoid paradise fish (Macropodus opercularis), an "air-breather" species from Southeast Asia, has a highly complex behavioral repertoire and has been the subject of many ethological investigations, but lacks genomic resources. Here we report the reference genome assembly of Macropodus opercularis using long-read sequences at 150-fold coverage. The final assembly consisted of ≈483 Mb on 152 contigs. Within the assembled genome we identified and annotated 20,157 protein coding genes and assigned ≈90% of them to orthogroups. Completeness analysis showed that 98.5% of the Actinopterygii core gene set (ODB10) was present as a complete ortholog in our reference genome with a further 1.2 % being present in a fragmented form. Additionally, we cloned multiple genes important during early development and using newly developed in situ hybridization protocols, we showed that they have conserved expression patterns.

Advances in Vertebrate (Cyto)Genomics Shed New Light on Fish Compositional Genome Evolution

Cytogenetic and compositional studies considered fish genomes rather poor in guanine-cytosine content (GC%) because of a putative "sharp increase in genic GC% during the evolution of higher vertebrates". However, the available genomic data have not been exploited to confirm this viewpoint. In contrast, further misunderstandings in GC%, mostly of fish genomes, originated from a misapprehension of the current flood of data. Utilizing public databases, we calculated the GC% in animal genomes of three different, technically well-established fractions: DNA (entire genome), cDNA (complementary DNA), and cds (exons). Our results across chordates help set borders of GC% values that are still incorrect in literature and show: (i) fish in their immense diversity possess comparably GC-rich (or even GC-richer) genomes as higher vertebrates, and fish exons are GC-enriched among vertebrates; (ii) animal genomes generally show a GC-enrichment from the DNA, over cDNA, to the cds level (i.e., not only the higher vertebrates); (iii) fish and invertebrates show a broad(er) inter-quartile range in GC%, while avian and mammalian genomes are more constrained in their GC%. These results indicate no sharp increase in the GC% of genes during the transition to higher vertebrates, as stated and numerously repeated before. We present our results in 2D and 3D space to explore the compositional genome landscape and prepared an online platform to explore the AT/GC compositional genome evolution.

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

Intronic regions of the human coagulation factor VIII gene harboring transcription factor binding sites with a strong bias towards the short-interspersed elements

Increasing data show that intronic derived regulatory elements, such as transcription factor binding sites (TFBs), play key roles in gene regulation, and malfunction. Accordingly, characterizing the sequence context of the intronic regions of the human coagulation factor VIII (hFVIII) gene can be important. In this study, the intronic regions of the hFVIII gene were scrutinized based on in-silico methods. The results disclosed that these regions harbor a rich array of functional elements such as repetitive elements (REs), splicing sites, and transcription factor binding sites (TFBs). Among these elements, TFBs and REs showed a significant distribution and correlation to each other. This survey indicated that 31% of TFBs are localized in the intronic regions of the gene. Moreover, TFBs indicate a strong bias in the regions far from splice sites of introns with mapping to different REs. Accordingly, TFBs showed highly bias toward Short Interspersed Elements (SINEs), which in turn they covering about 12% of the total of REs. However, the distribution pattern of TFBs-REs showed different bias in the intronic regions, spatially into the Introns 13 and 25. The rich array of SINE-TFBs and CR1-TFBs were situated within 5′UTR of the gene that may be an important driving force for regulatory innovation of the hFVIII gene. Taken together, these data may lead to revealing intronic regions with the capacity to renewing gene regulatory networks of the hFVIII gene. On the other hand, these correlations might provide the novel idea for a new hypothesis of molecular evolution of the FVIII gene, and treatment of Hemophilia A which should be considered in future studies.

Hybrid genome assembly and annotation of Danionella translucida

Studying neuronal circuits at cellular resolution is very challenging in vertebrates due to the size and optical turbidity of their brains. Danionella translucida, a close relative of zebrafish, was recently introduced as a model organism for investigating neural network interactions in adult individuals. Danionella remains transparent throughout its life, has the smallest known vertebrate brain and possesses a rich repertoire of complex behaviours. Here we sequenced, assembled and annotated the Danionella translucida genome employing a hybrid Illumina/Nanopore read library as well as RNA-seq of embryonic, larval and adult mRNA. We achieved high assembly continuity using low-coverage long-read data and annotated a large fraction of the transcriptome. This dataset will pave the way for molecular research and targeted genetic manipulation of this novel model organism.

The Most Developmentally Truncated Fishes Show Extensive Hox Gene Loss and Miniaturized Genomes

The world’s smallest fishes belong to the genus Paedocypris. These miniature fishes are endemic to an extreme habitat: the peat swamp forests in Southeast Asia, characterized by highly acidic blackwater. This threatened habitat is home to a large array of fishes, including a number of miniaturized but also developmentally truncated species. Especially the genus Paedocypris is characterized by profound, organism-wide developmental truncation, resulting in sexually mature individuals of <8 mm in length with a larval phenotype. Here, we report on evolutionary simplification in the genomes of two species of the dwarf minnow genus Paedocypris using whole-genome sequencing. The two species feature unprecedented Hox gene loss and genome reduction in association with their massive developmental truncation. We also show how other genes involved in the development of musculature, nervous system, and skeleton have been lost in Paedocypris, mirroring its highly progenetic phenotype. Further, our analyses suggest two mechanisms responsible for the genome streamlining in Paedocypris in relation to other Cypriniformes: severe intron shortening and reduced repeat content. As the first report on the genomic sequence of a vertebrate species with organism-wide developmental truncation, the results of our work enhance our understanding of genome evolution and how genotypes are translated to phenotypes. In addition, as a naturally simplified system closely related to zebrafish, Paedocypris provides novel insights into vertebrate development.

Evolutionary divergence of 3' UTRs in cichlid fishes

Background

Post-transcriptional regulation is crucial for the control of eukaryotic gene expression and might contribute to adaptive divergence. The three prime untranslated regions (3’ UTRs), that are located downstream of protein-coding sequences, play important roles in post-transcriptional regulation. These regions contain functional elements that influence the fate of mRNAs and could be exceptionally important in groups such as rapidly evolving cichlid fishes.

Results

To examine cichlid 3’ UTR evolution, we 1) identified gene features in nine teleost genomes and 2) performed comparative analyses to assess evolutionary variation in length, functional motifs, and evolutionary rates of 3’ UTRs. In all nine teleost genomes, we found a smaller proportion of repetitive elements in 3’ UTRs than in the whole genome. We found that the 3’ UTRs in cichlids tend to be longer than those in non-cichlids, and this was associated, on average, with one more miRNA target per gene in cichlids. Moreover, we provided evidence that 3’ UTRs on average have evolved faster in cichlids than in non-cichlids. Finally, analyses of gene function suggested that both the top 5% longest and 5% most rapidly evolving 3’ UTRs in cichlids tended to be involved in ribosome-associated pathways and translation.

Conclusions

Our results reveal novel patterns of evolution in the 3’ UTRs of teleosts in general and cichlids in particular. The data suggest that 3’ UTRs might serve as important meta-regulators, regulators of other mechanisms governing post-transcriptional regulation, especially in groups like cichlids that have undergone extremely fast rates of phenotypic diversification and speciation.

Electronic supplementary material

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

Population genetic diversity in zebrafish lines

Toxicological and pharmacological researchers have seized upon the many benefits of zebrafish, including the short generation time, well-characterized development, and early maturation as clear embryos. A major difference from many model organisms is that standard husbandry practices in zebrafish are designed to maintain population diversity. While this diversity is attractive for translational applications in human and ecological health, it raises critical questions on how interindividual genetic variation might contribute to chemical exposure or disease susceptibility differences. Findings from pooled samples of zebrafish support this supposition of diversity yet cannot directly measure allele frequencies for reference versus alternate alleles. Using the Tanguay lab Tropical 5D zebrafish line (T5D), we performed whole genome sequencing on a large group (n = 276) of individual zebrafish embryos. Paired-end reads were collected on an Illumina 3000HT, then aligned to the most recent zebrafish reference genome (GRCz10). These data were used to compare observed population genetic variation across species (humans, mice, zebrafish), then across lines within zebrafish. We found more single nucleotide polymorphisms (SNPs) in T5D than have been reported in SNP databases for any of the WIK, TU, TL, or AB lines. We theorize that some subset of the novel SNPs may be shared with other zebrafish lines but have not been identified in other studies due to the limitations of capturing population diversity in pooled sequencing strategies. We establish T5D as a model that is representative of diversity levels within laboratory zebrafish lines and demonstrate that experimental design and analysis can exert major effects when characterizing genetic diversity in heterogeneous populations.

COGNATE: comparative gene annotation characterizer

BACKGROUND

The comparison of gene and genome structures across species has the potential to reveal major trends of genome evolution. However, such a comparative approach is currently hampered by a lack of standardization (e.g., Elliott TA, Gregory TR, Philos Trans Royal Soc B: Biol Sci 370:20140331, 2015). For example, testing the hypothesis that the total amount of coding sequences is a reliable measure of potential proteome diversity (Wang M, Kurland CG, Caetano-Anollés G, PNAS 108:11954, 2011) requires the application of standardized definitions of coding sequence and genes to create both comparable and comprehensive data sets and corresponding summary statistics. However, such standard definitions either do not exist or are not consistently applied. These circumstances call for a standard at the descriptive level using a minimum of parameters as well as an undeviating use of standardized terms, and for software that infers the required data under these strict definitions. The acquisition of a comprehensive, descriptive, and standardized set of parameters and summary statistics for genome publications and further analyses can thus greatly benefit from the availability of an easy to use standard tool.

RESULTS

We developed a new open-source command-line tool, COGNATE (Comparative Gene Annotation Characterizer), which uses a given genome assembly and its annotation of protein-coding genes for a detailed description of the respective gene and genome structure parameters. Additionally, we revised the standard definitions of gene and genome structures and provide the definitions used by COGNATE as a working draft suggestion for further reference. Complete parameter lists and summary statistics are inferred using this set of definitions to allow down-stream analyses and to provide an overview of the genome and gene repertoire characteristics. COGNATE is written in Perl and freely available at the ZFMK homepage ( https://www.zfmk.de/en/COGNATE ) and on github ( https://github.com/ZFMK/COGNATE ).

CONCLUSION

The tool COGNATE allows comparing genome assemblies and structural elements on multiples levels (e.g., scaffold or contig sequence, gene). It clearly enhances comparability between analyses. Thus, COGNATE can provide the important standardization of both genome and gene structure parameter disclosure as well as data acquisition for future comparative analyses. With the establishment of comprehensive descriptive standards and the extensive availability of genomes, an encompassing database will become possible.

Transcriptome analysis of pancreatic cells across distant species highlights novel important regulator genes

BACKGROUND

Defining the transcriptome and the genetic pathways of pancreatic cells is of great interest for elucidating the molecular attributes of pancreas disorders such as diabetes and cancer. As the function of the different pancreatic cell types has been maintained during vertebrate evolution, the comparison of their transcriptomes across distant vertebrate species is a means to pinpoint genes under strong evolutionary constraints due to their crucial function, which have therefore preserved their selective expression in these pancreatic cell types.

RESULTS

In this study, RNA-sequencing was performed on pancreatic alpha, beta, and delta endocrine cells as well as the acinar and ductal exocrine cells isolated from adult zebrafish transgenic lines. Comparison of these transcriptomes identified many novel markers, including transcription factors and signaling pathway components, specific for each cell type. By performing interspecies comparisons, we identified hundreds of genes with conserved enriched expression in endocrine and exocrine cells among human, mouse, and zebrafish. This list includes many genes known as crucial for pancreatic cell formation or function, but also pinpoints many factors whose pancreatic function is still unknown. A large set of endocrine-enriched genes can already be detected at early developmental stages as revealed by the transcriptomic profiling of embryonic endocrine cells, indicating a potential role in cell differentiation. The actual involvement of conserved endocrine genes in pancreatic cell differentiation was demonstrated in zebrafish for myt1b, whose invalidation leads to a reduction of alpha cells, and for cdx4, selectively expressed in endocrine delta cells and crucial for their specification. Intriguingly, comparison of the endocrine alpha and beta cell subtypes from human, mouse, and zebrafish reveals a much lower conservation of the transcriptomic signatures for these two endocrine cell subtypes compared to the signatures of pan-endocrine and exocrine cells. These data suggest that the identity of the alpha and beta cells relies on a few key factors, corroborating numerous examples of inter-conversion between these two endocrine cell subtypes.

CONCLUSION

This study highlights both evolutionary conserved and species-specific features that will help to unveil universal and fundamental regulatory pathways as well as pathways specific to human and laboratory animal models such as mouse and zebrafish.

The genetic map of goldfish (Carassius auratus) provided insights to the divergent genome evolutions in the Cyprinidae family

A high-density linkage map of goldfish (Carassius auratus) was constructed using RNA-sequencing. This map consists of 50 linkage groups with 8,521 SNP markers and an average resolution of 0.62 cM. Approximately 84% of markers are in protein-coding genes orthologous to zebrafish proteins. We performed comparative genome analysis between zebrafish and medaka, common carp, grass carp, and goldfish to study the genome evolution events in the Cyprinidae family. The comparison revealed large synteny blocks among Cyprinidae fish and we hypothesized that the Cyprinidae ancestor undergone many inter-chromosome rearrangements after speciation from teleost ancestor. The study also showed that goldfish genome had one more round of whole genome duplication (WGD) than zebrafish. Our results illustrated that most goldfish markers were orthologous to genes in common carp, which had four rounds of WGD. Growth-related regions and genes were identified by QTL analysis and association study. Function annotations of the associated genes suggested that they might regulate development and growth in goldfish. This first genetic map enables us to study the goldfish genome evolution and provides an important resource for selective breeding of goldfish.

Direct Visualization of DNA Replication Dynamics in Zebrafish Cells

Spatiotemporal regulation of DNA replication in the S-phase nucleus has been extensively studied in mammalian cells because it is tightly coupled with the regulation of other nuclear processes such as transcription. However, little is known about the replication dynamics in nonmammalian cells. Here, we analyzed the DNA replication processes of zebrafish (Danio rerio) cells through the direct visualization of replicating DNA in the nucleus and on DNA fiber molecules isolated from the nucleus. We found that zebrafish chromosomal DNA at the nuclear interior was replicated first, followed by replication of DNA at the nuclear periphery, which is reminiscent of the spatiotemporal regulation of mammalian DNA replication. However, the relative duration of interior DNA replication in zebrafish cells was longer compared to mammalian cells, possibly reflecting zebrafish-specific genomic organization. The rate of replication fork progression and ori-to-ori distance measured by the DNA combing technique were ∼ 1.4 kb/min and 100 kb, respectively, which are comparable to those in mammalian cells. To our knowledge, this is a first report that measures replication dynamics in zebrafish cells.

Rapid genomic DNA changes in allotetraploid fish hybrids

Rapid genomic change has been demonstrated in several allopolyploid plant systems; however, few studies focused on animals. We addressed this issue using an allotetraploid lineage (4nAT) of freshwater fish originally derived from the interspecific hybridization of red crucian carp (Carassius auratus red var., ♀, 2n=100) × common carp (Cyprinus carpio L., ♂, 2n=100). We constructed a bacterial artificial chromosome (BAC) library from allotetraploid hybrids in the 20th generation (F20) and sequenced 14 BAC clones representing a total of 592.126 kb, identified 11 functional genes and estimated the guanine-cytosine content (37.10%) and the proportion of repetitive elements (17.46%). The analysis of intron evolution using nine orthologous genes across a number of selected fish species detected a gain of 39 introns and a loss of 30 introns in the 4nAT lineage. A comparative study based on seven functional genes among 4nAT, diploid F1 hybrids (2nF1) (first generation of hybrids) and their original parents revealed that both hybrid types (2nF1 and 4nAT) not only inherited genomic DNA from their parents, but also demonstrated rapid genomic DNA changes (homoeologous recombination, parental DNA fragments loss and formation of novel genes). However, 4nAT presented more genomic variations compared with their parents than 2nF1. Interestingly, novel gene fragments were found for the iqca1 gene in both hybrid types. This study provided a preliminary genomic characterization of allotetraploid F20 hybrids and revealed evolutionary and functional genomic significance of allopolyploid animals.

Increased functional protein expression using nucleotide sequence features enriched in highly expressed genes in zebrafish

Many genetic manipulations are limited by difficulty in obtaining adequate levels of protein expression. Bioinformatic and experimental studies have identified nucleotide sequence features that may increase expression, however it is difficult to assess the relative influence of these features. Zebrafish embryos are rapidly injected with calibrated doses of mRNA, enabling the effects of multiple sequence changes to be compared in vivo. Using RNAseq and microarray data, we identified a set of genes that are highly expressed in zebrafish embryos and systematically analyzed for enrichment of sequence features correlated with levels of protein expression. We then tested enriched features by embryo microinjection and functional tests of multiple protein reporters. Codon selection, releasing factor recognition sequence and specific introns and 3′ untranslated regions each increased protein expression between 1.5- and 3-fold. These results suggested principles for increasing protein yield in zebrafish through biomolecular engineering. We implemented these principles for rational gene design in software for codon selection (CodonZ) and plasmid vectors incorporating the most active non-coding elements. Rational gene design thus significantly boosts expression in zebrafish, and a similar approach will likely elevate expression in other animal models.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Molecular cloning and functional characterization of zebrafish Slc4a3/Ae3 anion exchanger

The zebrafish genome encodes two slc4a1 genes, one expressed in erythroid tissues and the other in the HR (H(+)-ATPase-rich) type of embryonic skin ionocytes, and two slc4a2 genes, one in proximal pronephric duct and the other in several extrarenal tissues of the embryo. We now report cDNA cloning and functional characterization of zebrafish slc4a3/ae3 gene products. The single ae3 gene on chromosome 9 generates at least two low-abundance ae3 transcripts differing only in their 5'-untranslated regions and encoding a single definitive Ae3 polypeptide of 1170 amino acids. The 7 kb upstream of the apparent initiator Met in ae3 exon 3 comprises multiple diverse, mobile repeat elements which disrupt and appear to truncate the Ae3 N-terminal amino acid sequence that would otherwise align with brain Ae3 of other species. Embryonic ae3 mRNA expression was detected by whole mount in situ hybridization only in fin buds at 24-72 hpf, but was detectable by RT-PCR across a range of embryonic and adult tissues. Epitope-tagged Ae3 polypeptide was expressed at or near the surface of Xenopus oocytes, and mediated low rates of DIDS-sensitive (36)Cl(-)/Cl(-) exchange in influx and efflux assays. As previously reported for Ae2 polypeptides, (36)Cl(-) transport by Ae3 was inhibited by both extracellular and intracellular acidic pH, and stimulated by alkaline pH. However, zebrafish Ae3 differed from Ae2 polypeptides in its insensitivity to NH4Cl and to hypertonicity. We conclude that multiple repeat elements have disrupted the 5'-end of the zebrafish ae3 gene, associated with N-terminal truncation of the protein and reduced anion transport activity.

Identification, expression profiling of a grass carp TLR8 and its inhibition leading to the resistance to reovirus in CIK cells

TLR8 (toll-like receptor 8), a homolog of TLR3, TLR7 and TLR9 as prototypical intracellular members of TLR family, is generally associated with sensing single stranded RNA and plays a pivotal role in antiviral immune response. In this study, a TLR8 gene from grass carp Ctenopharyngodon idella (designated as CiTLR8) was obtained and characterized. The full-length cDNA of CiTLR8 was of 3766 bp. The open reading frame was of 3072 bp and encoded a polypeptide of 1023 amino acids, including seventeen LRR (leucine-rich repeat) motifs, one transmembrane domain and one TIR (toll/interleukin-1 receptor) domain. A single intron with the size of 839 bp was found on the neck of start codon (ATG). CiTLR8 mRNA was ubiquitously expressed in the 15 tested tissues and the expression level in gas bladder, spleen, brain, hindgut and trunk kidney tissues was high. Besides, the CiTLR8 expression in spleen and head kidney was significantly up-regulated and reached peak at 24 h post-injection of grass carp reovirus (GCRV). CiTLR8 transcription reached peak at 8 h and then declined below the normal level post-GCRV infection in the C. idella kidney (CIK) cell line; and it was rapidly and significantly down-regulated by the stimulation of the synthetic double-stranded RNA polyriboinosinic-polyribocytidylic acid sodium salt (poly I:C) in CIK cells in a dose and time-dependent manner. The inhibitor expression vectors were constructed and transfected into CIK cell line to obtain stably expressing shRNA targeting TLR8. In CiTLR8-knockdown cells, CiTLR7 transcript weakly increased, CiIFN-I mRNA was significantly down-regulated, and the expression of CiMyD88, CiIRF7 and CiMx1 scarcely changed. Post poly I:C stimulation, CiTLR8, CiTLR7 and CiMyD88 transcripts significantly increased, CiIRF7 was down-regulated after an initial phase of increase, and CiIFN-I and CiMx1 transcripts were up-regulated. After GCRV infection, the transcripts of CiTLR8, CiTLR7, CiMyD88 and CiIRF7 were up-regulated, but CiIFN-I and CiMx1 transcripts were inhibited. Nevertheless, cells transfected with pshTLR8 vectors had strong resistance against GCRV injection, suggesting CiTLR8 might play a negative role in antiviral immune response. These results collectively suggested that CiTLR8 was a novel member of TLR gene family, engaging in antiviral innate immune defense in C. idella, which laid a foundation for the further mechanism research of TLR8 in fishes.

Function-valued traits in evolution

Many biological characteristics of evolutionary interest are not scalar variables but continuous functions. Given a dataset of function-valued traits generated by evolution, we develop a practical, statistical approach to infer ancestral function-valued traits, and estimate the generative evolutionary process. We do this by combining dimension reduction and phylogenetic Gaussian process regression, a non-parametric procedure that explicitly accounts for known phylogenetic relationships. We test the performance of methods on simulated, function-valued data generated from a stochastic evolutionary model. The methods are applied assuming that only the phylogeny, and the function-valued traits of taxa at its tips are known. Our method is robust and applicable to a wide range of function-valued data, and also offers a phylogenetically aware method for estimating the autocorrelation of function-valued traits.

In with the old, in with the new: the promiscuity of the duplication process engenders diverse pathways for novel gene creation

The gene duplication process has exhibited far greater promiscuity in the creation of paralogs with novel exon-intron structures than anticipated even by Ohno. In this paper I explore the history of the field, from the neo-Darwinian synthesis through Ohno's formulation of the canonical model for the evolution of gene duplicates and culminating in the present genomic era. I delineate the major tenets of Ohno's model and discuss its failure to encapsulate the full complexity of the duplication process as revealed in the era of genomics. I discuss the diverse classes of paralogs originating from both DNA- and RNA-mediated duplication events and their evolutionary potential for assuming radically altered functions, as well as the degree to which they can function unconstrained from the pressure of gene conversion. Lastly, I explore theoretical population-genetic considerations of how the effective population size (N(e)) of a species may influence the probability of emergence of genes with radically altered functions.

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.