A non-adaptationist perspective on evolution of genomic complexity or the continued dethroning of man

A new, non-adaptationist theory of evolution of genomic complexity was recently proposed by Lynch and Conery. This concept holds that increase in complexity seen in eukaryotic genomes is a 'syndrome' caused by increase in genome entropy, which is inevitably triggered by reduction of population size. Here, I discuss the definitions of genomic entropy and complexity and the evidence supporting the entropic theory of genome complexity evolution, including new observations on concordant gain and loss of genes and introns in eukaryotic genomes. I further consider the far-reaching biological and philosophical implications of this theory.     

Rearrangement rate following the whole-genome duplication in teleosts

It is now clear that a whole-genome duplication (WGD) occurred at the base of the teleost fish lineage. Like the other anciently polyploid genomes investigated so far, teleost genomes now behave like diploids with chromosomes forming pairs at meiosis. The diploidization process is currently poorly understood. It is associated with many gene deletions, such that one of the duplicates is lost at most loci and has also been proposed to coincide with an increase in genomic instability. Here we ask whether WGD is a determinant of the genomic rearrangement rate in teleosts. We study variability of the rates of rearrangement along a vertebrate phylogenetic tree, composed of 3 tetrapods (human, chicken, and mouse) and 3 teleost fishes (zebrafish, Tetraodon, and Takifugu), whose complete genome sequences are available. We devise a simple parsimony method for counting rearrangements, which takes into account various methodological complications caused by the WGD and the subsequent gene losses. We show that there does appear to be an increase in rearrangement rate after WGD, but that there is also a great deal of additional variability in rearrangement rates across species.     

Comparison of genomic DNA sequences: solved and unsolved problems

MOTIVATION: The DNA sequences of entire genomes are being determined at a rapid rate. Whereas initial genome sequencing efforts were for organisms chosen to be widely spaced in the tree of life, there is a growing emphasis on projects to sequence a species that is sufficiently similar to an already-sequenced species to allow direct comparison of those two DNA sequences. This and other changes in genome sequencing strategies have created a strong need for new methods to compare genomic sequences. RESULTS: We sketch the current state of software for comparing genomic DNA sequences and outline research directions that we believe are likely to result in important advances in practice.     

Evolutionary force of AT‐rich repeats to trap genomic and episomal DNAs into the rice genome: lessons from endogenous pararetrovirus

In plant genomes, the incorporation of DNA segments is not a common method of artificial gene transfer. Nevertheless, various segments of pararetroviruses have been found in plant genomes in recent decades. The rice genome contains a number of segments of endogenous rice tungro bacilliform virus‐like sequences (ERTBVs), many of which are present between AT dinucleotide repeats (ATrs). Comparison of genomic sequences between two closely related rice subspecies, japonica and indica, allowed us to verify the preferential insertion of ERTBVs into ATrs. In addition to ERTBVs, the comparative analyses showed that ATrs occasionally incorporate repeat sequences including transposable elements, and a wide range of other sequences. Besides the known genomic sequences, the insertion sequences also represented DNAs of unclear origins together with ERTBVs, suggesting that ATrs have integrated episomal DNAs that would have been suspended in the nucleus. Such insertion DNAs might be trapped by ATrs in the genome in a host‐dependent manner. Conversely, other simple mono‐ and dinucleotide sequence repeats (SSR) were less frequently involved in insertion events relative to ATrs. Therefore, ATrs could be regarded as hot spots of double‐strand breaks that induce non‐homologous end joining. The insertions within ATrs occasionally generated new gene‐related sequences or involved structural modifications of existing genes. Likewise, in a comparison between Arabidopsis thaliana and Arabidopsis lyrata, the insertions preferred ATrs to other SSRs. Therefore ATrs in plant genomes could be considered as genomic dumping sites that have trapped various DNA molecules and may have exerted a powerful evolutionary force.     

Length and GC Content Variability of Introns among Teleostean Genomes in the Light of the Metabolic Rate Hypothesis

A comparative analysis of five teleostean genomes, namely zebrafish, medaka, three-spine stickleback, fugu and pufferfish was performed with the aim to highlight the nature of the forces driving both length and base composition of introns (i.e., bpi and GCi). An inter-genome approach using orthologous intronic sequences was carried out, analyzing independently both variables in pairwise comparisons. An average length shortening of introns was observed at increasing average GCi values. The result was not affected by masking transposable and repetitive elements harbored in the intronic sequences. The routine metabolic rate (mass specific temperature-corrected using the Boltzmann''s factor) was measured for each species. A significant correlation held between average differences of metabolic rate, length and GC content, while environmental temperature of fish habitat was not correlated with bpi and GCi. Analyzing the concomitant effect of both variables, i.e., bpi and GCi, at increasing genomic GC content, a decrease of bpi and an increase of GCi was observed for the significant majority of the intronic sequences (from ∼40% to ∼90%, in each pairwise comparison). The opposite event, concomitant increase of bpi and decrease of GCi, was counter selected (from <1% to ∼10%, in each pairwise comparison). The results further support the hypothesis that the metabolic rate plays a key role in shaping genome architecture and evolution of vertebrate genomes.     

MRD: a microsatellite repeats database for prokaryotic and eukaryotic genomes

MRD is a database system to access the microsatellite repeats information of genomes such as archea, eubacteria, and other eukaryotic genomes whose sequence information is available in public domains. MRD stores information about simple tandemly repeated k-mer sequences where k= 1 to 6, i.e. monomer to hexamer. The web interface allows the users to search for the repeat of their interest and to know about the association of the repeat with genes and genomic regions in the specific organism. The data contains the abundance and distribution of microsatellites in the coding and non-coding regions of the genome. The exact location of repeats with respect to genomic regions of interest (such as UTR, exon, intron or intergenic regions) whichever is applicable to organism is highlighted. MRD is available on the World Wide Web at and/or . The database is designed as an open-ended system to accommodate the microsatellite repeats information of other genomes whose complete sequences will be available in future through public domain.     

Amino Acid Usage Is Asymmetrically Biased in AT- and GC-Rich Microbial Genomes

Introduction

Genomic base composition ranges from less than 25% AT to more than 85% AT in prokaryotes. Since only a small fraction of prokaryotic genomes is not protein coding even a minor change in genomic base composition will induce profound protein changes. We examined how amino acid and codon frequencies were distributed in over 2000 microbial genomes and how these distributions were affected by base compositional changes. In addition, we wanted to know how genome-wide amino acid usage was biased in the different genomes and how changes to base composition and mutations affected this bias. To carry this out, we used a Generalized Additive Mixed-effects Model (GAMM) to explore non-linear associations and strong data dependences in closely related microbes; principal component analysis (PCA) was used to examine genomic amino acid- and codon frequencies, while the concept of relative entropy was used to analyze genomic mutation rates.

Results

We found that genomic amino acid frequencies carried a stronger phylogenetic signal than codon frequencies, but that this signal was weak compared to that of genomic %AT. Further, in contrast to codon usage bias (CUB), amino acid usage bias (AAUB) was differently distributed in AT- and GC-rich genomes in the sense that AT-rich genomes did not prefer specific amino acids over others to the same extent as GC-rich genomes. AAUB was also associated with relative entropy; genomes with low AAUB contained more random mutations as a consequence of relaxed purifying selection than genomes with higher AAUB.

Conclusion

Genomic base composition has a substantial effect on both amino acid- and codon frequencies in bacterial genomes. While phylogeny influenced amino acid usage more in GC-rich genomes, AT-content was driving amino acid usage in AT-rich genomes. We found the GAMM model to be an excellent tool to analyze the genomic data used in this study.     

A novel common triplet profile for GC-rich prokaryotic genomes

It has been reported that there is a majority triplet profile among genomes, which was considered as a reflection of general mechanisms of genome evolution (Albrecht-Buehler, 2007). However, there are actually, according to our further analysis and at least among prokaryotic genomes, two common triplet profiles: one is from low-GC content genomes; the other is from high-GC content genomes. Both common profiles would be direct reflections of GC content variations and strand symmetry of genomic sequences.     

A kinetic model for subtractive hybridization.

Nucleic acid sequences that differ in abundance between two populations (target sequences) can be cloned by multiple rounds of subtractive hybridization and amplification by PCR. These sequences can be cDNAs representing up-regulated mRNAs, or genomic DNAs from deletion mutants. We have derived an equation that describes the recovery of such sequences, and have used this to simulate the outcome of up to 10 rounds of subtractive hybridization and PCR amplification. When the model was tested by comparing its predictions with the published results from genomic and cDNA subtractions, the predictions of the model were generally in good agreement with the published data. We have modelled the outcomes of genomic subtractions, for a variety of genomes, and have used it to compare various strategies for enriching targets. The model predicts that for genomes of less than 5 x 10(8) bp, deletions of as small as 1 kbp should represent > 99% of the DNA after three to six rounds of hybridization (depending on the enrichment procedure). As genomes increase in size, the kinetics of hybridization become an important limiting factor. However, even for genomes as large as 3 x 10(9) bp, it should be possible to isolate deletions of 5 kbp using the appropriate conditions. These simulations suggest that such methods offer a realistic alternative to chromosome walking for identifying genomic deletions for which there are known phenotypes, thereby considerably reducing time and effort. For cDNA subtractive hybridization, the model predicts that after six rounds of hybridization, sequences that do not differ in abundance between the tester and driver populations (the background) will represent < 1% of the subtracted population, and even quite modestly upregulated cDNAs should be successfully enriched. Where several up-regulated cDNAs are present, the predicted final representation is dependent on both the initial abundance and the degree of up-regulation.     

Characterization of species-specific repeats in 613 prokaryotic species

Prokaryotes are in general believed to possess small, compactly organized genomes, with repetitive sequences forming only a small part of them. Nonetheless, many prokaryotic genomes in fact contain species-specific repeats (>85 bp long genomic sequences with less than 60% identity to other species) as we have previously demonstrated. However, it is not known at present how frequent such species-specific repeats are and what their functional roles in bacterial genomes may be. Therefore, we have conducted a comprehensive survey of prokaryotic species-specific repeats and characterized them to examine as to whether there are functional classes among different repeats or not and how they are mutually related to each other. Of the 613 distinct prokaryotic species analyzed, 97% were found to contain at least one species-specific repeats. It seems interesting to note that the species-specific repeats thus identified appear to be functionally variable in different genomes: in some genomes, they are mostly associated with duplicated protein-coding genes, whereas in some other genomes with rRNA and tRNA genes. Contrary to what may be expected, only one-fourth of the species-specific repeats were found to be associated with mobile genetic elements.     

Genetic Relationship and Genomic in situ Hybridization Analysis of the Three Genomes in Triticum aestivum

Common wheat ( Triticum　aestivum L.) is an allohexaploid, consisting of three different genomes (Au, B and D ) which are genetically closely related. Genomic DNA of the three possible genome donors, T. urartu Thum., Aegilops speltoides Tausch and Ae. tauschii Coss.,were employed as probes to hybridize with the diploid genomic DNA digested by Eco RⅠand Hin dⅢ respectively. Both the hybridization strength and band patterns among the genomes would be good indicators of genome relationships. Combining distr ibution data of some repetitive DNA sequences cloned from T. urartu in the three genomes, the authors draw a conclusion that Au and D are more closely related to each other than either one to the B genome. Genomic in situ hybridization (GISH) of T. aestivum cv. Chinese Spring with genomic DNA probes of the three diploid progenitors respectively indicated that the three genomes could be discriminated clearly via GISH. The signals on the chromosomes of Au and D genomes were even. However, when Ae. speltoides DNA was used as probe, there were very strong cross hybridization and the signals condensed on some areas of the metaphasic chromosomes. In the interphase nucleus, the chromatin of B genome dispersed on the same region and the signals on the homologous chromosomes distributed symmetrically. Rich repetitive DNA sequences in B genome, especially the tandem repetitives, perhaps take an important role for the formation of the special hybridization pattern. The main difference between B and the other two genomes probably is in the repetitive DNA sequences.     

A thermal denaturation study of genomic DNAs from North American minnows (Cyprinidae: Teleostei)

Base compositions and differential melting rate profiles of genomic DNAs from twenty species of North American cyprinid fishes were generated via thermal denaturation. Base pair composition expressed as % GC values ranged among the twenty species from 36.1–41.3%. This range is considerably broader than that observed at comparable taxonomic levels in other vertebrate groups. Both the range and average difference in base pair composition between species in the diverse and rapidly evolving genus Notropis were considerably greater than those between species in other North American cyprinid genera. This may indicate that genomic changes at the level of base pair composition are frequent and possibly important events in cyprinid evolution. Compositional heterogeneity and asymmetry values among the twenty species were uniform and low, respectively, suggesting that most of the species lacked DNA components in their genomes which differed substantially from their main-band DNAs in base pair composition. The melting rate profiles revealed a prominent and distinct heavy or GC-rich DNA component in the genomes of three species belonging to the subgenus Cyprinella of Notropis. These and other data suggest that the heavy melting component may reflect a large, comparatively GC-rich family of highly repeated or satellite DNA sequences common to all three genomes.     

Prediction of recognition sites for genomic replication of classical swine fever virus with information analysis

In order to explore the mechanism for the genomic replication of classical swine fever virus (CSFV), so as to make a basis for investigating its pathogenicity, an introduction of the information theory is presented in connection with the statistical mechanics, whence small-sample statistics appears naturally as a consequence of the Bayesian approach. Furthermore, a selection rule for identifying the pattern of a recognition site for an RNA-binding protein is proposed by means of the maximum entropy principle. Based on those, the information contents of 3'-untranslated regions (3'UTRs) of genomes of 20 CSFV strains and 5'-untranslated regions (5'UTRs) of genomes of 58 CSFV strains are analyzed with a computational algorithm in a reduction mode, and the 3'UTR sites of 20 strains and 5'UTR sites of 58 strains containing important motifs are extracted from the unaligned RNA sequences of unequal lengths. These sites, which have the patterns of sequence and structure similar to the putative cis elements related to the regulation of genomic replication, would be identified as the potential recognition sites in 3'UTRs and 5'UTRs for CSFV replicase responsible for classical swine fever virus genomic replication, and to some extent, this identification is supported by experimental evidence. Finally, information analysis allows a presumption to be made about the CSFV RNA replication initiation mechanism.     

Evolution of Germline-Limited Sequences in Two Populations of the Ciliate Chilodonella uncinata

Ciliates are microbial eukaryotes that separate their nuclear functions into a germline micronucleus and a somatic macronucleus. During development of the macronucleus the genome undergoes a series of reorganization events that includes the precise excision of intervening DNA. Here, we determine the architecture of four loci in the micronuclear and macronuclear genomes of the ciliate Chilodonella uncinata and compare the levels of variation in micronuclear-limited sequences to macronuclear destined sequences at two of these loci. We find that within a population, germline-limited sequences are evolving at the same rate as other putatively neutral sites, but between populations germline-limited sequences are accumulating mutations at a much faster rate than other sites. We also find evidence of macronuclear recombination and incomplete elimination of intervening DNA, which result in increased diversity in the macronuclear genome. Our results support the assertion that the unusual genomic features of ciliates can result in rapid and unpredicted patterns of diversification.     

真核生物DNA非编码区的组分分析

在全基因组水平上,用直方图、混沌表示灰度图、距离差异度和信息熵差异度四种方法,研究了拟南芥、线虫、果蝇的DNA内含子、基因间隔区DNA、外显子三种区域的核苷酸短序列组分及组分复杂度.结果表明：a.不同基因组之间,不管基因数目多少,用4种方法得到的外显子部分其组分复杂度都比较接近,而非编码区部分的组分复杂度却很大.这一点定量地说明了物种之间的复杂程度,主要不体现在编码区部分,而体现在非编码区部分.b.同一基因组中,内含子的核苷酸短序列组分复杂度都是相似的,外显子和intergenic DNA部分的组分复杂度也是相似的.c.内含子和intergenic DNA在转录、剪切、二级结构等方面有很大的不同,但它们在核苷酸短序列组分上的差异却很小,说明内含子和intergenic DNA在转录、剪切、二级结构上的不同并不通过核苷酸短序列组分来进行限制.     

Microbial genome evolution: sources of variability

Comparative genome analyses of close relatives have yielded exciting insight into the sources of microbial genome variability with respect to gene content, gene order and evolution of genes with unknown functions. The genomes of free-living bacteria often carry phages and repetitive sequences that mediate genomic rearrangements in contrast to the small genomes of obligate host-associated bacteria. This suggests that genomic stability correlates with the genomic content of repeated sequences and movable genetic elements, and thereby with bacterial lifestyle. Genes with unknown functions present in a single species tend to be shorter than conserved, functional genes, indicating that the fraction of unique genes in microbial genomes has been overestimated.     

Four basic symmetry types in the universal 7-cluster structure of microbial genomic sequences

Coding information is the main source of heterogeneity (non-randomness) in the sequences of microbial genomes. The heterogeneity corresponds to a cluster structure in triplet distributions of relatively short genomic fragments (200-400 bp). We found a universal 7-cluster structure in microbial genomic sequences and explained its properties. We show that codon usage of bacterial genomes is a multi-linear function of their genomic G+C-content with high accuracy. Based on the analysis of 143 completely sequenced bacterial genomes available in Genbank in August 2004, we show that there are four "pure" types of the 7-cluster structure observed. All 143 cluster animated 3D-scatters are collected in a database which is made available on our web-site (http://www.ihes.fr/~zinovyev/7clusters). The findings can be readily introduced into software for gene prediction, sequence alignment or microbial genomes classification.     

Isolation of A/D and C genome specific dispersed and clustered repetitive DNA sequences from Avena sativa.

DNA gel-blot and in situ hybridization with genome-specific repeated sequences have proven to be valuable tools in analyzing genome structure and relationships in species with complex allopolyploid genomes such as hexaploid oat (Avena sativa L., 2n = 6x = 42; AACCDD genome). In this report, we describe a systematic approach for isolating genome-, chromosome-, and region-specific repeated and low-copy DNA sequences from oat that can presumably be applied to any complex genome species. Genome-specific DNA sequences were first identified in a random set of A. sativa genomic DNA cosmid clones by gel-blot hybridization using labeled genomic DNA from different Avena species. Because no repetitive sequences were identified that could distinguish between the A and D gneomes, sequences specific to these two genomes are refereed to as A/D genome specific. A/D or C genome specific DNA subfragments were used as screening probes to identify additional genome-specific cosmid clones in the A. sativa genomic library. We identified clustered and dispersed repetitive DNA elements for the A/D and C genomes that could be used as cytogenetic markers for discrimination of the various oat chromosomes. Some analyzed cosmids appeared to be composed entirely of genome-specific elements, whereas others represented regions with genome- and non-specific repeated sequences with interspersed low-copy DNA sequences. Thus, genome-specific hybridization analysis of restriction digests of random and selected A. sativa cosmids also provides insight into the sequence organization of the oat genome.