Identification of genome constitution of Oryza malampuzhaensis, O. minuta,and O. punctata by multicolor genomic in situ hybridization

Multicolor genomic in situ hybridization (McGISH) was applied to identify the genomic constitution of three tetraploid species (2n = 4x = 48) in the Oryza officinalis complex of the genus Oryza, i.e. Oryza malam-puzhaensis, Oryza minuta, and Oryza punctata. The genomic probes used were from three diploids, i.e. Oryza officinalis (CC), Oryza eichingeri (CC) and Oryza punctata (BB), respectively. The results indicated that all three tetraploids are allotetraploid with the genomic constitution of BBCC, and among them the genome constitution of O. malampuzhaensis was verified for the first time. Restoration of the independent taxonomic status of O. malampuzhaensis is suggested. One pair of satellite chromosomes belonging to the B genome was identified in O. malampuzhaensis, but no such satellite chromosomes were found in either O. minuta or the tetraploid O. punctata. The average chromosome length of the C genome was found to be slightly larger than that of the B-genome chromosomes of O. minuta, but not in the tetraploids O. punctata and O. malampuzhaensis. McGISH also revealed that the B genome of O. minuta and the B genome of diploid O. punctata showed clear differentiation from each other. Therefore, the suggestion was proposed that the B genome in diploid O. punctata was not the source of the B genome of O. minuta. The present results proved that multicolor GISH had high resolution in identifying the genomic constitution of polyploid Oryza species. Received: 14 February 2000 / Accepted: 13 November 2000     

Genetic diversity and evolutionary relationships ofOryza species with the B- and C-genomes as revealed by SSR markers

Genetic diversity and evolutionary relationships of 72 accessions representing six species with the B-, C-, and BC-genomes in the genusOryza were investigated by seven microsatellite markers. Of four diploid species,Oryza officinalis maintained the highest diversity (P=71.4%, He=0.565), followed by Oryza eichingeri (P=57.1%, He=0.376), Oryzapunctata (P=57.1%, He=0.272) and Oryza rhizomatis (P=42.9%, He=0.222). In comparison, a higher level of genetic diversity was revealed in the tetraploid (P=71.4%, He=0.461-0.637). UPGMA dendrograms based on genetic distance revealed an obvious genetic differentiation between Asian and African races ofO. eichingeri. Three BBCC species clustered with different accessions of the diploidO. punctata, suggestive of their multiple origins. The results inferred from the dendrogram suggested that diploid species,O. officinalis and AfricanO. eichingeri might be the C-genome donors for tetraploid species,Oryza minuta andO. punctata, respectively, while the C-genome ancestor ofOryza malampuzhaensis seemed to be eitherO. rhizomatis or the Sri LankanO. eichingeri species. The genetic relationship among the CC and BBCC species further indicated that the tetraploid species with the BC-genome have originated independently, at least three times in history. In addition, we have demonstrated successful cross-species amplification of seven rice SSR loci acrossOryza species with B-and C-genomes.     

Singular value decomposition-based regression identifies activation of endogenous signaling pathways in vivo

Background

The physical organization and chromosomal localization of genes within genomes is known to play an important role in their function. Most genes arise by duplication and move along the genome by random shuffling of DNA segments. Higher order structuring of the genome occurs in eukaryotes, where groups of physically linked genes are co-expressed. However, the contribution of gene duplication to gene order has not been analyzed in detail, as it is believed that co-expression due to recent duplicates would obscure other domains of co-expression.

Results

We have catalogued ordered duplicated genes in Drosophila melanogaster, and found that one in five of all genes is organized as tandem arrays. Furthermore, among arrays that have been spatially conserved over longer periods than would be expected on the basis of random shuffling, a disproportionate number contain genes encoding developmental regulators. Using in situ gene expression data for more than half of the Drosophila genome, we find that genes in these conserved clusters are co-expressed to a much higher extent than other duplicated genes.

Conclusions

These results reveal the existence of functional constraints in insects that retain copies of genes encoding developmental and regulatory proteins as neighbors, allowing their co-expression. This co-expression may be the result of shared cis-regulatory elements or a shared need for a specific chromatin structure. Our results highlight the association between genome architecture and the gene regulatory networks involved in the construction of the body plan.     

Phylogenetic analysis of Oryza species, based on simple sequence repeats and their flanking nucleotide sequences from the mitochondrial and chloroplast genomes

Simple sequence repeats (SSR) and their flanking regions in the mitochondrial and chloroplast genomes were sequenced in order to reveal DNA sequence variation. This information was used to gain new insights into phylogenetic relationships among species in the genus Oryza. Seven mitochondrial and five chloroplast SSR loci equal to or longer than ten mononucleotide repeats were chosen from known rice mitochondrial and chloroplast genome sequences. A total of 50 accessions of Oryza that represented six different diploid genomes and three different allopolyploid genomes of Oryza species were analyzed. Many base substitutions and deletions/insertions were identified in the SSR loci as well as their flanking regions. Of mononucleotide SSR, G (or C) repeats were more variable than A (or T) repeats. Results obtained by chloroplast and mitochondrial SSR analyses showed similar phylogenetic relationships among species, although chloroplast SSR were more informative because of their higher sequence diversity. The CC genome is suggested to be the maternal parent for the two BBCC genome species (O. punctata and O. minuta) and the CCDD species O. latifolia, based on the high level of sequence conservation between the diploid CC genome species and these allotetraploid species. This is the first report of phylogenetic analysis among plant species, based on mitochondrial and chloroplast SSR and their flanking sequences.     

A Tourist element in the 5′-flanking region of the catalase gene CatA reveals evolutionary relationships among Oryza species with various genome types

Tourist-OsaCatA, a transposable element, was found in the 5′-flanking region of the rice gene CatA. The characteristics of this element are similar to those of the other Tourist elements so far found in Oryza sativa. PCR and sequence analyses of 37 accessions of 18 species revealed that all the Oryza species examined, except for one accession, have either a full-length or a partial Tourist element at this locus. Unlike the Tourist elements previously reported, this Tourist element is found in all four Oryza species complexes in the Oryzeae tribe. All AA genome Oryza species, except O.?longistaminata, contain the full-length Tourist element. O. longistaminata and the species of the O. officinalis, O. meyeriana and O.?ridleyi complexes contain the partial element. A phylogenetic tree of Oryza species based on the nucleotide sequences of these Tourist elements was constructed. The O.?longistaminata accessions were placed near the neighboring cluster of the officinalis complex. We propose that the ancestor of O.?longistaminata and that of other species with the AA genome diverged, and the ancestor(s) of the O.?officinalis, O.?ridleyi and O.?meyeriana complexes then diverged from the ancestor of O.?longistaminata in the course of the evolution of the Oryza species. The Tourist elements associated with CatA and its orthologs thus provide useful tools for examining evolutionary relationships among Oryza species.     

Polymorphism and phylogenetic relationships among species in the genus Oryza as determined by analysis of nuclear RFLPs

Summary Ninety-three accessions representing 21 species from the genus Oryza were examined for restriction fragment length polymorphism. The majority (78%) of the accessions, for which five individuals were tested, were found to be monomorphic. Most of the polymorphic accessions segregated for only one or two probes and appeared to be mixed pure lines. For most of the Oryza species tested, the majority of the genetic variation (83%) was found between accessions from different species with only 17% between accessions within species. Tetraploid species were found to have, on average, nearly 50% more alleles (unique fragments) per individual than diploid species reflecting the allopolyploid nature of their genomes.Classification of Oryza species based on RFLPs matches remarkably well previous classifications based on morphology, hybridization and isozymes. In the current study, four species complexes could be identified corresponding to those proposed by Vaughan (1989): the O. ridleyi complex, the O. meyeriana complex, the O. officinalis complex and the O. sativa complex. Within the O. sativa complex, accessions of O. rufipogon from Asia (including O. nivara) and perennial forms of O. rufipogon from Australia clustered together with accessions of cultivated rice O. sativa. Surprisingly, indica and japonica (the two major subspecies of cultivated rice) showed closer affinity with different accessions of wild O. Rufipogon than to each other, supporting a hypothesis of independent domestication events for these two types of rice. Australian annual wild rice O. meridionalis (previously classified as O. rufipogon) was clearly distinct from all other O. rufipogon accessions supporting its recent reclassification as O. meridionalis (Ng et al. 1981). Using genetic relatedness as a criterion, it was possible to identify the closest living diploid relatives of the currently known tetraploid rice species. Results from these analyses suggest that BBCC tetraploids (O. malampuzhaensis, O. punctata and O. minuta) are either of independent origins or have experienced introgression from sympatric C-genome diploid rice species. CCDD tetraploid species from America (O. latifolia, O. alta and O. grandiglumis) may be of ancient origin since they show a closer affinity to each other than to any known diploid species. Their closest living diploid relatives belong to C genome (O. eichingeri) and E genome (O. Australiensis) species. Comparisons among African, Australian and Asian rice species suggest that Oryza species in Africa and Australia are of polyphyletic origin and probably migrated to these regions at different times in the past.Finally, on a practical note, the majority of probes used in this study detected polymorphism between cultivated rice and its wild relatives. Hence, RFLP markers and maps based on such markers are likely to be very useful in monitoring and aiding introgression of genes from wild rice into modern cultivars.     

Genome evolution in Oryza allopolyploids of various ages: Insights into the process of diploidization

The prevalence and recurrence of whole-genome duplication in plants and its major role in evolution have been well recognized. Despite great efforts, many aspects of genome evolution, particularly the temporal progression of genomic responses to allopolyploidy and the underlying mechanisms, remain poorly understood. The rice genus Oryza consists of both recently formed and older allopolyploid species, representing an attractive system for studying the genome evolution after allopolyploidy. In this study, through screening BAC libraries and sequencing and annotating the targeted BAC clones, we generated orthologous genomic sequences surrounding the DEP1 locus, a major grain yield QTL in cultivated rice, from four Oryza polyploids of various ages and their likely diploid genome donors or close relatives. Based on sequenced DEP1 region and published data from three other genomic regions, we investigated the temporal evolutionary dynamics of four polyploid genomes at both genetic and expression levels. In the recently formed BBCC polyploid, Oryza minuta, genome dominance was not observed and its short-term responses to allopolyploidy are mainly manifested as a high proportion of homoeologous gene pairs showing unequal expression. This could partly be explained by parental legacy, rewiring of divergent regulatory networks and epigenetic modulation. Moreover, we detected an ongoing diploidization process in this genus, and suggest that the expression divergence driven by changes of selective constraint probably plays a big role in the long-term diploidization. These findings add novel insights into our understanding of genome evolution after allopolyploidy, and could facilitate crop improvements through hybridization and polyploidization.     

Dynamic Evolution of Oryza Genomes Is Revealed by Comparative Genomic Analysis of a Genus-Wide Vertical Data Set

Oryza (23 species; 10 genome types) contains the world's most important food crop — rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past ~15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.     

Extensive genomic diversity and selective conservation of virulence-determinants in enterohemorrhagic Escherichia coli strains of O157 and non-O157 serotypes

Background

Enterohemorrhagic Escherichia coli (EHEC) O157 causes severe food-borne illness in humans. The chromosome of O157 consists of 4.1 Mb backbone sequences shared by benign E. coli K-12, and 1.4 Mb O157-specific sequences encoding many virulence determinants, such as Shiga toxin genes (stx genes) and the locus of enterocyte effacement (LEE). Non-O157 EHECs belonging to distinct clonal lineages from O157 also cause similar illness in humans. According to the 'parallel' evolution model, they have independently acquired the major virulence determinants, the stx genes and LEE. However, the genomic differences between O157 and non-O157 EHECs have not yet been systematically analyzed.

Results

Using microarray and whole genome PCR scanning analyses, we performed a whole genome comparison of 20 EHEC strains of O26, O111, and O103 serotypes with O157. In non-O157 EHEC strains, although genome sizes were similar with or rather larger than O157 and the backbone regions were well conserved, O157-specific regions were very poorly conserved. Around only 20% of the O157-specific genes were fully conserved in each non-O157 serotype. However, the non-O157 EHECs contained a significant number of virulence genes that are found on prophages and plasmids in O157, and also multiple prophages similar to, but significantly divergent from, those in O157.

Conclusion

Although O157 and non-O157 EHECs have independently acquired a huge amount of serotype- or strain-specific genes by lateral gene transfer, they share an unexpectedly large number of virulence genes. Independent infections of similar but distinct bacteriophages carrying these virulence determinants are deeply involved in the evolution of O157 and non-O157 EHECs.     

Two C or not two C: recurrent disruption of Zn-ribbons, gene duplication, lineage-specific gene loss, and horizontal gene transfer in evolution of bacterial ribosomal proteins

Background

Ribosomal proteins are encoded in all genomes of cellular life forms and are, generally, well conserved during evolution. In prokaryotes, the genes for most ribosomal proteins are clustered in several highly conserved operons, which ensures efficient co-regulation of their expression. Duplications of ribosomal-protein genes are infrequent, and given their coordinated expression and functioning, it is generally assumed that ribosomal-protein genes are unlikely to undergo horizontal transfer. However, with the accumulation of numerous complete genome sequences of prokaryotes, several paralogous pairs of ribosomal protein genes have been identified. Here we analyze all such cases and attempt to reconstruct the evolutionary history of these ribosomal proteins.

Results

Complete bacterial genomes were searched for duplications of ribosomal proteins. Ribosomal proteins L36, L33, L31, S14 are each duplicated in several bacterial genomes and ribosomal proteins L11, L28, L7/L12, S1, S15, S18 are so far duplicated in only one genome each. Sequence analysis of the four ribosomal proteins, for which paralogs were detected in several genomes, two of the ribosomal proteins duplicated in one genome (L28 and S18), and the ribosomal protein L32 showed that each of them comes in two distinct versions. One form contains a predicted metal-binding Zn-ribbon that consists of four conserved cysteines (in some cases replaced by histidines), whereas, in the second form, these metal-chelating residues are completely or partially replaced. Typically, genomes containing paralogous genes for these ribosomal proteins encode both versions, designated C+ and C-, respectively. Analysis of phylogenetic trees for these seven ribosomal proteins, combined with comparison of genomic contexts for the respective genes, indicates that in most, if not all cases, their evolution involved a duplication of the ancestral C+ form early in bacterial evolution, with subsequent alternative loss of the C+ and C- forms in different lineages. Additionally, evidence was obtained for a role of horizontal gene transfer in the evolution of these ribosomal proteins, with multiple cases of gene displacement 'in situ', that is, without a change of the gene order in the recipient genome.

Conclusions

A more complex picture of evolution of bacterial ribosomal proteins than previously suspected is emerging from these results, with major contributions of lineage-specific gene loss and horizontal gene transfer. The recurrent theme of emergence and disruption of Zn-ribbons in bacterial ribosomal proteins awaits a functional interpretation.     

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Duplicated genes frequently experience asymmetric rates of sequence evolution. Relaxed selective constraints and positive selection have both been invoked to explain the observation that one paralog within a gene-duplicate pair exhibits an accelerated rate of sequence evolution. In the majority of studies where asymmetric divergence has been established, there is no indication as to which gene copy, ancestral or derived, is evolving more rapidly. In this study we investigated the effect of local synteny (gene-neighborhood conservation) and codon usage on the sequence evolution of gene duplicates in the S. cerevisiae genome. We further distinguish the gene duplicates into those that originated from a whole-genome duplication (WGD) event (ohnologs) versus small-scale duplications (SSD) to determine if there exist any differences in their patterns of sequence evolution.

Results

For SSD pairs, the derived copy evolves faster than the ancestral copy. However, there is no relationship between rate asymmetry and synteny conservation (ancestral-like versus derived-like) in ohnologs. mRNA abundance and optimal codon usage as measured by the CAI is lower in the derived SSD copies relative to ancestral paralogs. Moreover, in the case of ohnologs, the faster-evolving copy has lower CAI and lowered expression.

Conclusions

Together, these results suggest that relaxation of selection for codon usage and gene expression contribute to rate asymmetry in the evolution of duplicated genes and that in SSD pairs, the relaxation of selection stems from the loss of ancestral regulatory information in the derived copy.     

Molecular diversity in the primary and secondary gene pools of genus Oryza

The objective of the present investigation was to assess the genetic relationships among the species of Oryza that belong to the primary gene pool (sativa complex) and the secondary gene pool (officinalis complex) using three marker systems such as RAPDs, ISSRs and SSRs. A total of 432 clear and reproducible bands were amplified from 18 RAPD primers; 113 bands were detected from 8 ISSR primers and 78 alleles were found to be amplified across the Oryza species from 13 SSR primer pairs. All the three dendrograms constructed, using UPGMA from the genetic similarity matrices based on the three marker data sets, were similar in their groupings. In all the three trees, two accessions of Oryza sativa formed an exclusive group indicating its genomic differentiation from its wild ancestors through the process of domestication. Distinctness between the wild species of the sativa and officinalis complexes was evident in all the trees derived from different markers. The groupings obtained among the species of the sativa complex were in perfect concordance with the species relationships established through classical crossability and cytogenetic analysis. This study has brought out some information on the species relationship between the diploid and tetraploid genomes of the officinalis complex possessing BB, CC and DD genomes. The higher level of similarity observed between the species possessing C and D genomes supports the view of many earlier authors that these two genomes might have originated from a single hybridization event. The results of this study also show that the diploid species possessing C genomes such as Oryza officinalis, Oryza rhizomatis and Oryza eichingeri are distinct from their allotetraploid counterparts possessing BBCC and CCDD genomes indicating a wider genomic differentiation in their evolutionary process.     

Genomic and gene regulatory signatures of cryptozoic adaptation: Loss of blue sensitive photoreceptors through expansion of long wavelength-opsin expression in the red flour beetle Tribolium castaneum

Background

Hox genes are expressed in specific domains along the anterior posterior body axis and define the regional identity. In most animals these genes are organized in a single cluster in the genome and the order of the genes in the cluster is correlated with the anterior to posterior expression of the genes in the embryo. The conserved order of the various Hox gene orthologs in the cluster among most bilaterians implies that such a Hox cluster was present in their last common ancestor. Vertebrates are the only metazoans so far that have been shown to contain duplicated Hox clusters, while all other bilaterians seem to possess only a single cluster.

Results

We here show that at least three Hox genes of the spider Cupiennius salei are present as two copies in this spider. In addition to the previously described duplicated Ultrabithorax gene, we here present sequence and expression data of a second Deformed gene, and of two Sex comb reduced genes. In addition, we describe the sequence and expression of the Cupiennius proboscipedia gene. The spider Cupiennius salei is the first chelicerate for which orthologs of all ten classes of arthropod Hox genes have been described. The posterior expression boundary of all anterior Hox genes is at the tagma border of the prosoma and opisthosoma, while the posterior boundary of the posterior Hox genes is at the posterior end of the embryo.

Conclusion

The presence of at least three duplicated Hox genes points to a major duplication event in the lineage to this spider, perhaps even of the complete Hox cluster as has taken place in the lineage to the vertebrates. The combined data of all Cupiennius Hox genes reveal the existence of two distinct posterior expression boundaries that correspond to morphological tagmata boundaries.     

A physical map for the Amborella trichopoda genome sheds light on the evolution of angiosperm genome structure

Background

Recent phylogenetic analyses have identified Amborella trichopoda, an understory tree species endemic to the forests of New Caledonia, as sister to a clade including all other known flowering plant species. The Amborella genome is a unique reference for understanding the evolution of angiosperm genomes because it can serve as an outgroup to root comparative analyses. A physical map, BAC end sequences and sample shotgun sequences provide a first view of the 870 Mbp Amborella genome.

Results

Analysis of Amborella BAC ends sequenced from each contig suggests that the density of long terminal repeat retrotransposons is negatively correlated with that of protein coding genes. Syntenic, presumably ancestral, gene blocks were identified in comparisons of the Amborella BAC contigs and the sequenced Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera and Oryza sativa genomes. Parsimony mapping of the loss of synteny corroborates previous analyses suggesting that the rate of structural change has been more rapid on lineages leading to Arabidopsis and Oryza compared with lineages leading to Populus and Vitis. The gamma paleohexiploidy event identified in the Arabidopsis, Populus and Vitis genomes is shown to have occurred after the divergence of all other known angiosperms from the lineage leading to Amborella.

Conclusions

When placed in the context of a physical map, BAC end sequences representing just 5.4% of the Amborella genome have facilitated reconstruction of gene blocks that existed in the last common ancestor of all flowering plants. The Amborella genome is an invaluable reference for inferences concerning the ancestral angiosperm and subsequent genome evolution.     

Two new genomes in the Oryza complex identified on the basis of molecular divergence analysis using total genomic DNA hybridization

The genus Oryza to which cultivated rice belongs has 24 species (2n?=?24 or 48), representing seven genomes (AA, BB, CC, EE, FF, BBCC and CCDD). The genomic constitution of five of these species is unknown. These five species have been grouped into two species complexes, the tetraploid ridleyi complex (O. ridleyi, O.?longiglumis) and the diploid meyeriana complex (O.?granulata, O. meyeriana, O. indandamanica). To evaluate the genomic structure of these species in terms of divergence at the molecular level vis-à-vis other known genomes of Oryza, we used the total genomic DNA hybridization approach. Total genomic DNA (after restriction digestion) of 79 accessions of 23 Oryza species, 6 related genera, 5 outgroup taxa (2 monocots, 3 dicots) and 6 F₁s and BC₁s derived from crosses of O.?sativa with wild species were hybridized individually with ³²P-labeled total genomic DNA from 12 Oryza species: O. ridleyi, O.?longiglumis, O. granulata, O.?meyeriana, O. brachyantha, O. punctata, O. officinalis, O. eichingeri, O. alta, O. latifolia, O. australiensis, and O.?sativa. The labeled genomic DNAs representing the ridleyi and meyeriana complexes cross-hybridized best to all the accessions of their respective species, less to those representing other genomes of Oryza and related genera, and least to outgroup taxa. In general, the hybridization differential measured in terms of signal intensities was >50-fold under conditions that permit detection of 70–75% homologous sequences, both in the presence and in the absence of O. sativa DNA as competitor. In contrast, when total DNAs representing other Oryza genomes were used as probes, species of the O.?ridleyi and O.?meyeriana complexes did not show any significant cross-hybridization (<5%). These results demonstrate that the genome(s) of both of these complexes are highly diverged and distinct from all other known genomes of Oryza. We, therefore, propose new genomic designations for these two species complexes: GG for the diploid O. meyeriana complex and HHJJ for the allotetraploid O. ridleyi complex. The results also suggest that the uniqueness of these genomes is not restricted to species-specific highly repetitive DNA sequences, but also applies to dispersed sequences present in single or low to moderate copy numbers. Furthermore these appear to share relatively more genome-specific repeat sequences between themselves than with other genomes of rice. The study also demonstrates the potential of total genomic DNA hybridization as a simple but powerful tool, complementary to existing approaches, for ascertaining the genomic makeup of an organism.