A comparative phylogenetic approach for dating whole genome duplication events

Whole genome duplication of intra- and inter-chromosomes in the tomato genome

Whole genome duplication (WGD) events have been proven to occur in the evolutionary history of most angiosperms. Tomato is considered a model species of the Solanaceae family. In this study, we describe the details of the evolutionary process of the tomato genome by detecting collinearity blocks and dating the WGD events on the tree of life by combining two different methods: synonymous substitution rates (Ks) and phylogenetic trees. In total, 593 collinearity blocks were discovered out of 12 pseudo-chromosomes constructed. It was evident that chromosome 2 had experienced an intra-chromosomal duplication event. Major inter-chromosomal duplication occurred among all the pseudo-chromosome. We calculated the Ks value of these collinearity blocks. Two peaks of Ks distribution were found, corresponding to two WGD events occurring approximately 36-82 million years ago (MYA) and 148-205 MYA. Additionally, the results of phylogenetic trees suggested that the more recent WGD event may have occurred after the divergence of the rosid-asterid clade, but before the major diversification in Solanaceae. The older WGD event was shown to have occurred before the divergence of the rosid-asterid clade and after the divergence of rice-Arabidopsis (monocot-dicot).

Identification and characterization of full-length cDNAs in channel catfish (Ictalurus punctatus) and blue catfish (Ictalurus furcatus)

BACKGROUND

Genome annotation projects, gene functional studies, and phylogenetic analyses for a given organism all greatly benefit from access to a validated full-length cDNA resource. While increasingly common in model species, full-length cDNA resources in aquaculture species are scarce.

METHODOLOGY AND PRINCIPAL FINDINGS

Through in silico analysis of catfish (Ictalurus spp.) ESTs, a total of 10,037 channel catfish and 7,382 blue catfish cDNA clones were identified as potentially encoding full-length cDNAs. Of this set, a total of 1,169 channel catfish and 933 blue catfish full-length cDNA clones were selected for re-sequencing to provide additional coverage and ensure sequence accuracy. A total of 1,745 unique gene transcripts were identified from the full-length cDNA set, including 1,064 gene transcripts from channel catfish and 681 gene transcripts from blue catfish, with 416 transcripts shared between the two closely related species. Full-length sequence characteristics (ortholog conservation, UTR length, Kozak sequence, and conserved motifs) of the channel and blue catfish were examined in detail. Comparison of gene ontology composition between full-length cDNAs and all catfish ESTs revealed that the full-length cDNA set is representative of the gene diversity encoded in the catfish transcriptome.

CONCLUSIONS

This study describes the first catfish full-length cDNA set constructed from several cDNA libraries. The catfish full-length cDNA sequences, and data gleaned from sequence characteristics analysis, will be a valuable resource for ongoing catfish whole-genome sequencing and future gene-based studies of function and evolution in teleost fishes.

Locating large-scale gene duplication events through reconciled trees: implications for identifying ancient polyploidy events in plants

Recent analyses of plant genomic data have found extensive evidence of ancient whole genome duplication (or polyploidy) events, but there are many unresolved questions regarding the number and timing of such events in plant evolutionary history. We describe the first exact and efficient algorithm for the Episode Clustering problem, which, given a collection of rooted gene trees and a rooted species tree, seeks the minimum number of locations on the species tree of gene duplication events. Solving this problem allows one to place gene duplication events onto nodes of a given species tree and potentially detect large-scale gene duplication events. We examined the performance of an implementation of our algorithm using 85 plant gene trees that contain genes from a total of 136 plant taxa. We found evidence of large-scale gene duplication events in Populus, Gossypium, Poaceae, Asteraceae, Brassicaceae, Solanaceae, Fabaceae, and near the root of the eudicot clade that are consistent with previous genomic evidence. However, a lack of phylogenetic signal within the gene trees can produce erroneous evidence of large-scale duplication events, especially near the root of the species tree. Although the results of our algorithm should be interpreted cautiously, they provide hypotheses for precise locations of large-scale gene duplication events with data from relatively few gene trees and can complement other genomic approaches to provide a more comprehensive view of ancient large-scale gene duplication events.

Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition

Each mode of gene duplication (tandem, tetraploid, segmental, transpositional) retains genes in a biased manner. A reciprocal relationship exists between plant genes retained postpaleotetraploidy versus genes retained after an ancient tandem duplication. Among the models (C, neofunctionalization, balanced gene drive) and ideas that might explain this relationship, only balanced gene drive predicts reciprocity. The gene balance hypothesis explains that more "connected" genes--by protein-protein interactions in a heteromer, for example--are less likely to be retained as a tandem or transposed duplicate and are more likely to be retained postpaleotetraploidy; otherwise, selectively negative dosage effects are created. Biased duplicate retention is an instant and neutral by-product, a spandrel, of purifying selection. Balanced gene drive expanded plant gene families, including those encoding proteasomal proteins, protein kinases, motors, and transcription factors, with each paleotetraploidy, which could explain trends involving complexity. Balanced gene drive is a saltation mechanism in the mutationist tradition.

Gene structure dynamics and divergence of the polygalacturonase gene family of plants and fungus

Whole copies of the polygalacturonase (PG) genes from rice ( Oryza sativa subsp. japonica) and a filamentous fungus ( Aspergillus oryzae ) were isolated. The orthologs of the rice PGs were also retrieved from other plant species. The 106 plant PGs analyzed were divided into 5 clades, A, B, C, D, and E. The fungus PGs were classified into 3 clades, of which one formed a loose cluster with clade E of the plant PGs. Four domain motifs (I, II, III, IV) were identified in all PGs. Motifs II and III were split by introns such as G/DDC and CGPGHGIS/IGSLG, respectively. In plant PGs there were 446 introns in total and 3.98 introns per gene. Intron phase distribution was 65.5% for phase 0, 19.7% for phase 1, and 14.8% for phase 2 in plant PGs. In the PGs of A. oryzae there were 37 introns of phase 0 (59.5%), phase 1 (24.3%), and phase 2 (16.2%), with 2.47 introns per gene. The 5 clades of plant PGs were divided into 3 basic gene structure lineages. Intron positions and phases were conserved among the PGs in the first 2 lineages. The third lineage consisted of PGs of clade E, which also carried highly conserved introns at different positions from other PGs. Intron positions were not as highly conserved in fungus PGs as in plant PGs. The introns in the current PGs have been present since before the divergence of monocots from dicots. The results obtained show that differential losses of introns created gene diversity, which was followed by segmental and tandem duplication in plant PGs.

Comparison between a coffee single copy chromosomal region and Arabidopsis duplicated counterparts evidenced high level synteny between the coffee genome and the ancestral Arabidopsis genome

The Arabidopsis thaliana genome sequence provides a catalogue of reference genes that can be used for comparative analysis of other species thereby facilitating map-based cloning in economically important crops. We made use of a coffee bacterial artificial chromosome (BAC) contig linked to the S(H)3 leaf rust resistance gene to assess microsynteny between coffee (Coffea arabica L.) and Arabidopsis. Microsynteny was revealed and the matching counterparts to C. arabica contigs were seen to be scattered throughout four different syntenic segments of Arabidopsis on chromosomes (Ath) I, III, IV and V. Coffee BAC filter hybridizations were performed using coffee putative conserved orthologous sequences to Arabidopsis predicted genes located on the different Arabidopsis syntenic regions. The coffee BAC contig related to the S(H)3 region was successfully consolidated and later on validated by fingerprinting. Furthermore, the anchoring markers appeared in same order on the coffee BAC contigs and in all Arabidopsis segments with the exception of a single inversion on AtIII and AtIV Arabidopsis segments. However, the S(H)3 coffee region appears to be closer to the ancestral genome segment (before the divergence of Arabidopsis and coffee) than any of the duplicated counterparts in the present-day Arabidopsis genome. The genome duplication events at the origin of its Arabidopsis counterparts occurred most probably after the separation (i.e. 94 million years ago) of Euasterid (Coffee) and Eurosid (Arabidopsis).

Towards the ictalurid catfish transcriptome: generation and analysis of 31,215 catfish ESTs

Background

EST sequencing is one of the most efficient means for gene discovery and molecular marker development, and can be additionally utilized in both comparative genome analysis and evaluation of gene duplications. While much progress has been made in catfish genomics, large-scale EST resources have been lacking. The objectives of this project were to construct primary cDNA libraries, to conduct initial EST sequencing to generate catfish EST resources, and to obtain baseline information about highly expressed genes in various catfish organs to provide a guide for the production of normalized and subtracted cDNA libraries for large-scale transcriptome analysis in catfish.

Results

A total of 17 cDNA libraries were constructed including 12 from channel catfish (Ictalurus punctatus) and 5 from blue catfish (I. furcatus). A total of 31,215 ESTs, with average length of 778 bp, were generated including 20,451 from the channel catfish and 10,764 from blue catfish. Cluster analysis indicated that 73% of channel catfish and 67% of blue catfish ESTs were unique within the project. Over 53% and 50% of the channel catfish and blue catfish ESTs, respectively, had significant similarities to known genes. All ESTs have been deposited in GenBank. Evaluation of the catfish EST resources demonstrated their potential for molecular marker development, comparative genome analysis, and evaluation of ancient and recent gene duplications. Subtraction of abundantly expressed genes in a variety of catfish tissues, identified here, will allow the production of low-redundancy libraries for in-depth sequencing.

Conclusion

The sequencing of 31,215 ESTs from channel catfish and blue catfish has significantly increased the EST resources in catfish. The EST resources should provide the potential for microarray development, polymorphic marker identification, mapping, and comparative genome analysis.

Gene-balanced duplications, like tetraploidy, provide predictable drive to increase morphological complexity

Controversy surrounds the apparent rising maximums of morphological complexity during eukaryotic evolution, with organisms increasing the number and nestedness of developmental areas as evidenced by morphological elaborations reflecting area boundaries. No "predictable drive" to increase this sort of complexity has been reported. Recent genetic data and theory in the general area of gene dosage effects has engendered a robust "gene balance hypothesis," with a theoretical base that makes specific predictions as to gene content changes following different types of gene duplication. Genomic data from both chordate and angiosperm genomes fit these predictions: Each type of duplication provides a one-way injection of a biased set of genes into the gene pool. Tetraploidies and balanced segments inject bias for those genes whose products are the subunits of the most complex biological machines or cascades, like transcription factors (TFs) and proteasome core proteins. Most duplicate genes are removed after tetraploidy. Genic balance is maintained by not removing those genes that are dose-sensitive, which tends to leave duplicate "functional modules" as the indirect products (spandrels) of purifying selection. Functional modules are the likely precursors of coadapted gene complexes, a unit of natural selection. The result is a predictable drive mechanism where "drive" is used rigorously, as in "meiotic drive." Rising morphological gain is expected given a supply of duplicate functional modules. All flowering plants have survived at least three large-scale duplications/diploidizations over the last 300 million years (Myr). An equivalent period of tetraploidy and body plan evolution may have ended for animals 500 million years ago (Mya). We argue that "balanced gene drive" is a sufficient explanation for the trend that the maximums of morphological complexity have gone up, and not down, in both plant and animal eukaryotic lineages.

Following tetraploidy in an Arabidopsis ancestor, genes were removed preferentially from one homeolog leaving clusters enriched in dose-sensitive genes

Approximately 90% of Arabidopsis' unique gene content is found in syntenic blocks that were formed during the most recent whole-genome duplication. Within these blocks, 28.6% of the genes have a retained pair; the remaining genes have been lost from one of the homeologs. We create a minimized genome by condensing local duplications to one gene, removing transposons, and including only genes within blocks defined by retained pairs. We use a moving average of retained and non-retained genes to find clusters of retention and then identify the types of genes that appear in clusters at frequencies above expectations. Significant clusters of retention exist for almost all chromosomal segments. Detailed alignments show that, for 85% of the genome, one homeolog was preferentially (1.6x) targeted for fractionation. This homeolog fractionation bias suggests an epigenetic mechanism. We find that islands of retention contain "connected genes," those genes predicted-by the gene balance hypothesis-to be resistant to removal because the products they encode interact with other products in a dose-sensitive manner, creating a web of dependency. Gene families that are overrepresented in clusters include those encoding components of the proteasome/protein modification complexes, signal transduction machinery, ribosomes, and transcription factor complexes. Gene pair fractionation following polyploidy or segmental duplication leaves a genome enriched for "connected" genes. These clusters of duplicate genes may help explain the evolutionary origin of coregulated chromosomal regions and new functional modules.

Independent ancient polyploidy events in the sister families Brassicaceae and Cleomaceae

Recent studies have elucidated the ancient polyploid history of the Arabidopsis thaliana (Brassicaceae) genome. The studies concur that there was at least one polyploidy event occurring some 14.5 to 86 million years ago (Mya), possibly near the divergence of the Brassicaceae from its sister family, Cleomaceae. Using a comparative genomics approach, we asked whether this polyploidy event was unique to members of the Brassicaceae, shared with the Cleomaceae, or an independent polyploidy event in each lineage. We isolated and sequenced three genomic regions from diploid Cleome spinosa (Cleomaceae) that are each homoeologous to a duplicated region shared between At3 and At5, centered on the paralogs of SEPALLATA (SEP) and CONSTANS (CO). Phylogenetic reconstructions and analysis of synonymous substitution rates support the hypothesis that a genomic triplication in Cleome occurred independently of and more recently than the duplication event in the Brassicaceae. There is a strong correlation in the copy number (single versus duplicate) of individual genes, suggesting functionally consistent influences operating on gene copy number in these two independently evolving lineages. However, the amount of gene loss in Cleome is greater than in Arabidopsis. The genome of C. spinosa is only 1.9 times the size of A. thaliana, enabling comparative genome analysis of separate but related polyploidy events.

Genome histories clarify evolution of the expansin superfamily: new insights from the poplar genome and pine ESTs

Expansins comprise a superfamily of plant cell wall-loosening proteins that has been divided into four distinct families, EXPA, EXPB, EXLA and EXLB. In a recent analysis of Arabidopsis thaliana and Oryza sativa expansins, we proposed a further subdivision of the families into 17 clades, representing independent lineages in the last common ancestor of monocots and eudicots. This division was based on both traditional sequence-based phylogenetic trees and on position-based trees, in which genomic locations and dated segmental duplications were used to reconstruct gene phylogeny. In this article we review recent work concerning the patterns of expansin evolution in angiosperms and include additional insights gained from the genome of a second eudicot species, Populus trichocarpa, which includes at least 36 expansin genes. All of the previously proposed monocot-eudicot orthologous groups, but no additional ones, are represented in this species. The results also confirm that all of these clades are truly independent lineages. Furthermore, we have used position-based phylogeny to clarify the history of clades EXPA-II and EXPA-IV. Most of the growth of the expansin superfamily in the poplar lineage is likely due to a recent polyploidy event. Finally, some monocot-eudicot clades are shown to have diverged before the separation of the angiosperm and gymnosperm lineages.

Taking the first steps towards a standard for reporting on phylogenies: Minimum Information About a Phylogenetic Analysis (MIAPA)

In the eight years since phylogenomics was introduced as the intersection of genomics and phylogenetics, the field has provided fundamental insights into gene function, genome history and organismal relationships. The utility of phylogenomics is growing with the increase in the number and diversity of taxa for which whole genome and large transcriptome sequence sets are being generated. We assert that the synergy between genomic and phylogenetic perspectives in comparative biology would be enhanced by the development and refinement of minimal reporting standards for phylogenetic analyses. Encouraged by the development of the Minimum Information About a Microarray Experiment (MIAME) standard, we propose a similar roadmap for the development of a Minimal Information About a Phylogenetic Analysis (MIAPA) standard. Key in the successful development and implementation of such a standard will be broad participation by developers of phylogenetic analysis software, phylogenetic database developers, practitioners of phylogenomics, and journal editors.

Use of genomic history to improve phylogeny and understanding of births and deaths in a gene family

Polyploidy events have played an important role in the evolution of angiosperm genomes. Here, we demonstrate how genomic histories can increase phylogenetic resolution in a gene family, specifically the expansin superfamily of cell wall proteins. There are 36 expansins in Arabidopsis and 58 in rice. Traditional sequence-based phylogenetic trees yield poor resolution below the family level. To improve upon these analyses, we searched for gene colinearity (microsynteny) between Arabidopsis and rice genomic segments containing expansin genes. Multiple rounds of genome duplication and extensive gene loss have obscured synteny. However, by simultaneously aligning groups of up to 10 potentially orthologous segments from the two species, we traced the history of 49 out of 63 expansin-containing segments back to the ancestor of monocots and eudicots. Our results indicate that this ancestor had 15-17 expansin genes, each ancestral to an extant clade. Some clades have strikingly different growth patterns in the rice and Arabidopsis lineages, with more than half of all rice expansins arising from two ancestral genes. Segmental duplications, most of them part of polyploidy events, account for 12 out of 21 new expansin genes in Arabidopsis and 16 out of 44 in rice. Tandem duplications explain most of the rest. We were also able to estimate a minimum of 28 gene deaths in the Arabidopsis lineage and nine in rice. This analysis greatly clarifies expansin evolution since the last common ancestor of monocots and eudicots and the method should be broadly applicable to many other gene families.

Structure, evolution, and expression of the two invertase gene families of rice

Invertases catalyze the irreversible hydrolysis of sucrose to glucose and fructose. Plants contain two unrelated families of these enzymes: acid forms that derive from periplasmic invertases of eubacteria and are found in cell wall and vacuole, and neutral/alkaline forms evolved from the cytosolic invertases of cyanobacteria. Genomes of rice (Oryza sativa) and thale cress (Arabidopsis thaliana) contain multiple genes encoding these two families. Here for rice we identify the member genes of a cell-wall group (designated OsCIN1-9), a vacuolar group (OsVIN1-2), and two ancient neutral/alkaline groups: alpha (OsNIN1-4) and beta (OsNIN5-8). In Arabidopsis these groups contain six, two, four and five members, respectively. It is believed that the vacuolar group evolved from the cell-wall group. We provide evidence that the N-terminal signal peptide that directs cell-wall invertases co-translationally into the endoplasmic reticulum for secretion was replaced in the vacuolar group by a sequence similar to the complex N-terminal motif that targets alkaline phosphatase post-translationally to the vacuolar membrane of yeast. Since the last common ancestor of Arabidopsis and rice, the two invertase families evolved equally rapidly via gene duplication and gene loss, but the acid invertase family underwent approximately 10 events of intron loss compared with a single event of intron gain in the neutral/alkaline invertase family. Transcripts were detected for all rice invertase genes except OsCIN9. The acid invertase genes showed greater spatial and temporal diversity of expression than the neutral/alkaline genes.

Two ancient rounds of polyploidy in rice genome

An ancient genome duplication (PPP1) that predates divergence of the cereals has recently been recognized. We report here another potentially older large-scale duplication (PPP2) event that predates monocot-dicot divergence in the genome of rice (Oryza sativa L.), as inferred from the age distribution of pairs of duplicate genes based on recent genome data for rice. Our results suggest that paleopolyploidy was widespread and played an important role in the evolution of rice.

Computational approaches to unveiling ancient genome duplications

Recent analyses of complete genome sequences have revealed that many genomes have been duplicated in their evolutionary past. Such events have been associated with important biological transitions, major leaps in evolution and adaptive radiations of species. Here, we consider recently developed computational methods to detect such ancient large-scale gene duplication events. Several new approaches have been used to show that large-scale gene duplications are more common than previously thought.

Comparative genome analysis of monocots and dicots, toward characterization of angiosperm diversity

The importance of angiosperms to sustaining humanity by providing a wide range of 'ecosystem services' warrants increased exploration of their genomic diversity. The nearly completed sequences for two species representing the major angiosperm subclasses, specifically the dicot Arabidopsis thaliana and the monocot Oryza sativa, provide a foundation for comparative analysis across the angiosperms. The angiosperms also exemplify some challenges to be faced as genomics makes new inroads into describing biotic diversity, in particular polyploidy (genome-wide chromatin duplication), and much larger genome sizes than have been studied to date.

Ancient polyploidization predating divergence of the cereals, and its consequences for comparative genomics

Integration of structural genomic data from a largely assembled rice genome sequence, with phylogenetic analysis of sequence samples for many other taxa, suggests that a polyploidization event occurred approximately 70 million years ago, before the divergence of the major cereals from one another but after the divergence of the Poales from the Liliales and Zingiberales. Ancient polyploidization and subsequent "diploidization" (loss) of many duplicated gene copies has thus shaped the genomes of all Poaceae cereal, forage, and biomass crops. The Poaceae appear to have evolved as separate lineages for approximately 50 million years, or two-thirds of the time since the duplication event. Chromosomes that are predicted to be homoeologs resulting from this ancient duplication event account for a disproportionate share of incongruent loci found by comparison of the rice sequence to a detailed sorghum sequence-tagged site-based genetic map. Differential gene loss during diploidization may have contributed many of these incongruities. Such predicted homoeologs also account for a disproportionate share of duplicated sorghum loci, further supporting the hypothesis that the polyploidization event was common to sorghum and rice. Comparative gene orders along paleo-homoeologous chromosomal segments provide a means to make phylogenetic inferences about chromosome structural rearrangements that differentiate among the grasses. Superimposition of the timing of major duplication events on taxonomic relationships leads to improved understanding of comparative gene orders, enhancing the value of data from botanical models for crop improvement and for further exploration of genomic biodiversity. Additional ancient duplication events probably remain to be discovered in other angiosperm lineages.