Comparative genome analysis reveals extensive conservation of genome organisation for Arabidopsis thaliana and Capsella rubella

Origin and diversity of Capsella bursa-pastoris from the genomic point of view

BACKGROUND

Capsella bursa-pastoris, a cosmopolitan weed of hybrid origin, is an emerging model object for the study of early consequences of polyploidy, being a fast growing annual and a close relative of Arabidopsis thaliana. The development of this model is hampered by the absence of a reference genome sequence.

RESULTS

We present here a subgenome-resolved chromosome-scale assembly and a genetic map of the genome of Capsella bursa-pastoris. It shows that the subgenomes are mostly colinear, with no massive deletions, insertions, or rearrangements in any of them. A subgenome-aware annotation reveals the lack of genome dominance-both subgenomes carry similar number of genes. While most chromosomes can be unambiguously recognized as derived from either paternal or maternal parent, we also found homeologous exchange between two chromosomes. It led to an emergence of two hybrid chromosomes; this event is shared between distant populations of C. bursa-pastoris. The whole-genome analysis of 119 samples belonging to C. bursa-pastoris and its parental species C. grandiflora/rubella and C. orientalis reveals introgression from C. orientalis but not from C. grandiflora/rubella.

CONCLUSIONS

C. bursa-pastoris does not show genome dominance. In the earliest stages of evolution of this species, a homeologous exchange occurred; its presence in all present-day populations of C. bursa-pastoris indicates on a single origin of this species. The evidence coming from whole-genome analysis challenges the current view that C. grandiflora/rubella was a direct progenitor of C. bursa-pastoris; we hypothesize that it was an extinct (or undiscovered) species sister to C. grandiflora/rubella.

Perspective Research Progress in Cold Responses of Capsella bursa-pastoris

Plants respond to cold stress by modulating biochemical pathways and array of molecular events. Plant morphology is also affected by the onset of cold conditions culminating at repression in growth as well as yield reduction. As a preventive measure, cascades of complex signal transduction pathways are employed that permit plants to endure freezing or chilling periods. The signaling pathways and related events are regulated by the plant hormonal activity. Recent investigations have provided a prospective understanding about plant response to cold stress by means of developmental pathways e.g., moderate growth involved in cold tolerance. Cold acclimation assays and bioinformatics analyses have revealed the role of potential transcription factors and expression of genes like CBF, COR in response to low temperature stress. Capsella bursa-pastoris is a considerable model plant system for evolutionary and developmental studies. On different occasions it has been proved that C. bursa-pastoris is more capable of tolerating cold than A. thaliana. But, the mechanism for enhanced low or freezing temperature tolerance is still not clear and demands intensive research. Additionally, identification and validation of cold responsive genes in this candidate plant species is imperative for plant stress physiology and molecular breeding studies to improve cold tolerance in crops. We have analyzed the role of different genes and hormones in regulating plant cold resistance with special reference to C. bursa-pastoris. Review of collected data displays potential ability of Capsella as model plant for improvement in cold stress regulation. Information is summarized on cold stress signaling by hormonal control which highlights the substantial achievements and designate gaps that still happen in our understanding.

Genome-Wide Analysis of the Glutathione S-Transferase Gene Family in Capsella rubella: Identification, Expression, and Biochemical Functions

Extensive subfunctionalization might explain why so many genes have been maintained after gene duplication, which provides the engine for gene family expansion. However, it is still a particular challenge to trace the evolutionary dynamics and features of functional divergences in a supergene family over the course of evolution. In this study, we identified 49 Glutathione S-transferase (GST) genes from the Capsella rubella, a close relative of Arabidopsis thaliana and a member of the mustard family. Capsella GSTs can be categorized into eight classes, with tau and phi GSTs being the most numerous. The expansion of the two classes mainly occurs through tandem gene duplication, which results in tandem-arrayed gene clusters on chromosomes. By integrating phylogenetic analysis, expression patterns, and biochemical functions of Capsella and Arabidopsis GSTs, functional divergence, both in gene expression and enzymatic properties, were clearly observed in paralogous gene pairs in Capsella (even the most recent duplicates), and orthologous GSTs in Arabidopsis/Capsella. This study provides functional evidence for the expansion and organization of a large gene family in closely related species.

Microsynteny and phylogenetic analysis of tandemly organised miRNA families across five members of Brassicaceae reveals complex retention and loss history

Plant genomes are characterized by the presence of large miRNA gene families which are few in number. The expansion of miRNA families is thought to be driven by gene and genome duplication. Some members of these miRNA gene families are tandemly arranged and their analysis is of interest because such organisation may indicate origin through tandem duplication and also to investigate whether some such tandem clusters have similar expression patterns, and whether these are regulated through a common set of cis-regulatory elements (eg. promoters and enhancers). As a first step, we undertake a comprehensive study using micro-synteny analyses of tandemly organised miRNA families across the Brassicaceae spanning an evolutionary time scale of ca. 45 million years, among Arabidopsis, Capsella, Brassica and Thellungiella species, to address the following questions: Are most miRNA gene families present as tandem clusters? To what extent are these tandem patterns retained? To what extent can family sizes be ascribed to genome duplication? Our analysis of thirteen tandemly organised miRNA families revealed that synteny is largely conserved among Arabidopsis thaliana, A. lyrata and Capsella rubella, which form a clade spanning approximately between 6.2-9.8 my (Acarkan et al., 2000) [1]. On the other hand, comparison of sequences from these species with Brassica rapa, B. oleracea and Thellungiella halophila, which form a separate clade spanning 31 my (Franzke et al., 2011)[2] reveals many differences. The latter clade reveals several paralogous duplications that probably resulted from whole genome duplication, as well as disrupted synteny. Phylogenetic analyses of precursor sequences generally support the history inferred from synteny analysis. Synteny and phylogenetic analysis of six members of the tandemly organised miR169 family suggest that the Brassicaceae ancestral state consisted of a "dimer as a unit" which may have undergone direct local duplication to retain the transcriptional orientation followed by lineage specific changes. MiR169, to the best of our knowledge, is one of the largest tandemly organised miRNA gene family across plant kingdom and further analysis should reveal the generality of this pattern of evolution. The conserved organisation of miR395A-B-C and miR395 D-E-F as two clusters on same chromosome/scaffold across A. thaliana, B. rapa and salsuginea demonstrates retention of the large chromosomal segment across the two lineages. MiRNA family miR845 was detected only in Arabidopsis species and Thellungiella indicating a complex loss and retention history. MiR447A-B family was only found in A. thaliana indicating that it is a species-specific gene family of recent origin.

Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

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

Self-incompatibility in Papaver: advances in integrating the signalling network

Self-fertilization, which results in reduced fitness of offspring, is a common problem in hermaphrodite angiosperms. To prevent this, many plants utilize SI (self-incompatibility), which is determined by the multi-allelic S-locus, that allows discrimination between self (incompatible) and non-self (compatible) pollen by the pistil. In poppy (Papaver rhoeas), the pistil S-determinant (PrsS) is a small secreted protein which interacts with the pollen S-determinant PrpS, a ~20 kDa novel transmembrane protein. Interaction of matching pollen and pistil S-determinants results in self-recognition, initiating a Ca2+-dependent signalling network in incompatible pollen. This triggers several downstream events, including alterations to the cytoskeleton, phosphorylation of sPPases (soluble inorganic pyrophosphatases) and an MAPK (mitogen-activated protein kinase), increases in ROS (reactive oxygen species) and nitric oxide (NO), and activation of several caspase-like activities. This results in the inhibition of pollen tube growth, prevention of self-fertilization and ultimately PCD (programmed cell death) in incompatible pollen. The present review focuses on our current understanding of the integration of these signals with their targets in the SI/PCD network. We also discuss our recent functional expression of PrpS in Arabidopsis thaliana pollen.

Conservation and functional element discovery in 20 angiosperm plant genomes

Here, we describe the construction of a phylogenetically deep, whole-genome alignment of 20 flowering plants, along with an analysis of plant genome conservation. Each included angiosperm genome was aligned to a reference genome, Arabidopsis thaliana, using the LASTZ/MULTIZ paradigm and tools from the University of California-Santa Cruz Genome Browser source code. In addition to the multiple alignment, we created a local genome browser displaying multiple tracks of newly generated genome annotation, as well as annotation sourced from published data of other research groups. An investigation into A. thaliana gene features present in the aligned A. lyrata genome revealed better conservation of start codons, stop codons, and splice sites within our alignments (51% of features from A. thaliana conserved without interruption in A. lyrata) when compared with previous publicly available plant pairwise alignments (34% of features conserved). The detailed view of conservation across angiosperms revealed not only high coding-sequence conservation but also a large set of previously uncharacterized intergenic conservation. From this, we annotated the collection of conserved features, revealing dozens of putative noncoding RNAs, including some with recorded small RNA expression. Comparing conservation between kingdoms revealed a faster decay of vertebrate genome features when compared with angiosperm genomes. Finally, conserved sequences were searched for folding RNA features, including but not limited to noncoding RNA (ncRNA) genes. Among these, we highlight a double hairpin in the 5'-untranslated region (5'-UTR) of the PRIN2 gene and a putative ncRNA with homology targeting the LAF3 protein.

Comparative mapping, genomic structure, and expression analysis of eight pseudo-response regulator genes in Brassica rapa

Circadian clocks regulate plant growth and development in response to environmental factors. In this function, clocks influence the adaptation of species to changes in location or climate. Circadian-clock genes have been subject of intense study in models such as Arabidopsis thaliana but the results may not necessarily reflect clock functions in species with polyploid genomes, such as Brassica species, that include multiple copies of clock-related genes. The triplicate genome of Brassica rapa retains high sequence-level co-linearity with Arabidopsis genomes. In B. rapa we had previously identified five orthologs of the five known Arabidopsis pseudo-response regulator (PRR) genes that are key regulators of the circadian clock in this species. Three of these B. rapa genes, BrPRR1, BrPPR5, and BrPPR7, are present in two copies each in the B. rapa genome, for a total of eight B. rapa PRR (BrPRR) orthologs. We have now determined sequences and expression characteristics of the eight BrPRR genes and mapped their positions in the B. rapa genome. Although both members of each paralogous pair exhibited the same expression pattern, some variation in their gene structures was apparent. The BrPRR genes are tightly linked to several flowering genes. The knowledge about genome location, copy number variation and structural diversity of these B. rapa clock genes will improve our understanding of clock-related functions in this important crop. This will facilitate the development of Brassica crops for optimal growth in new environments and under changing conditions.

'Missing link' species Capsella orientalis and Capsella thracica elucidate evolution of model plant genus Capsella (Brassicaceae)

To elucidate the evolutionary history of the genus Capsella, we included the hitherto poorly known species C. orientalis and C. thracica into our studies together with C. grandiflora, C. rubella and C. bursa-pastoris. We sequenced the ITS and four loci of noncoding cpDNA regions (trnL - F, rps16, trnH -psbA and trnQ -rps16). Sequence data were evaluated with parsimony and Bayesian analyses. Divergence time estimates were carried out with the software package BEAST. We also performed isozyme, cytological, morphological and biogeographic studies. Capsella orientalis (self-compatible, SC; 2n = 16) forms a clade (eastern lineage) with C. bursa-pastoris (SC; 2n = 32), which is a sister clade (western lineage) to C. grandiflora (self-incompatible, SI; 2n = 16) and C. rubella (SC; 2n = 16). Capsella bursa-pastoris is an autopolyploid species of multiple origin, whereas the Bulgarian endemic C. thracica (SC; 2n = 32) is allopolyploid and emerged from interspecific hybridization between C. bursa-pastoris and C. grandiflora. The common ancestor of the two lineages was diploid and SI, and its distribution ranged from eastern Europe to central Asia, predominantly confined to steppe-like habitats. Biogeographic dynamics during the Pleistocene caused geographic and genetic subdivisions within the common ancestor giving rise to the two extant lineages.

The Papaver self-incompatibility pollen S-determinant, PrpS, functions in Arabidopsis thaliana

Many angiosperms use specific interactions between pollen and pistil proteins as "self" recognition and/or rejection mechanisms to prevent self-fertilization. Self-incompatibility (SI) is encoded by a multiallelic S locus, comprising pollen and pistil S-determinants. In Papaver rhoeas, cognate pistil and pollen S-determinants, PrpS, a pollen-expressed transmembrane protein, and PrsS, a pistil-expressed secreted protein, interact to trigger a Ca(2+)-dependent signaling network, resulting in inhibition of pollen tube growth, cytoskeletal alterations, and programmed cell death (PCD) in incompatible pollen. We introduced the PrpS gene into Arabidopsis thaliana, a self-compatible model plant. Exposing transgenic A. thaliana pollen to recombinant Papaver PrsS protein triggered remarkably similar responses to those observed in incompatible Papaver pollen: S-specific inhibition and hallmark features of Papaver SI. Our findings demonstrate that Papaver PrpS is functional in a species with no SI system that diverged ~140 million years ago. This suggests that the Papaver SI system uses cellular targets that are, perhaps, common to all eudicots and that endogenous signaling components can be recruited to elicit a response that most likely never operated in this species. This will be of interest to biologists interested in the evolution of signaling networks in higher plants.

Comparative analysis of rosaceous genomes and the reconstruction of a putative ancestral genome for the family

BACKGROUND

Comparative genome mapping studies in Rosaceae have been conducted until now by aligning genetic maps within the same genus, or closely related genera and using a limited number of common markers. The growing body of genomics resources and sequence data for both Prunus and Fragaria permits detailed comparisons between these genera and the recently released Malus × domestica genome sequence.

RESULTS

We generated a comparative analysis using 806 molecular markers that are anchored genetically to the Prunus and/or Fragaria reference maps, and physically to the Malus genome sequence. Markers in common for Malus and Prunus, and Malus and Fragaria, respectively were 784 and 148. The correspondence between marker positions was high and conserved syntenic blocks were identified among the three genera in the Rosaceae. We reconstructed a proposed ancestral genome for the Rosaceae.

CONCLUSIONS

A genome containing nine chromosomes is the most likely candidate for the ancestral Rosaceae progenitor. The number of chromosomal translocations observed between the three genera investigated was low. However, the number of inversions identified among Malus and Prunus was much higher than any reported genome comparisons in plants, suggesting that small inversions have played an important role in the evolution of these two genera or of the Rosaceae.

Recent speciation of Capsella rubella from Capsella grandiflora, associated with loss of self-incompatibility and an extreme bottleneck

Flowering plants often prevent selfing through mechanisms of self-incompatibility (S.I.). The loss of S.I. has occurred many times independently, because it provides short-term advantages in situations where pollinators or mates are rare. The genus Capsella, which is closely related to Arabidopsis, contains a pair of closely related diploid species, the self-incompatible Capsella grandiflora and the self-compatible Capsella rubella. To elucidate the transition to selfing and its relationship to speciation of C. rubella, we have made use of comparative sequence information. Our analyses indicate that C. rubella separated from C. grandiflora recently ( approximately 30,000-50,000 years ago) and that breakdown of S.I. occurred at approximately the same time. Contrasting the nucleotide diversity patterns of the 2 species, we found that C. rubella has only 1 or 2 alleles at most loci, suggesting that it originated through an extreme population bottleneck. Our data are consistent with diploid speciation by a single, selfing individual, most likely living in Greece. The new species subsequently colonized the Mediterranean by Northern and Southern routes, at a time that also saw the spread of agriculture. The presence of phenotypic diversity within modern C. rubella suggests that this species will be an interesting model to understand divergence and adaptation, starting from very limited standing genetic variation.

Role of model plant species

The use of model or reference species has played a major role in furthering detailed understanding of mechanisms and processes in the plant kingdom over the past 25 years. Species which have been adopted as models for dicotyledons and monocotyledons include arabidopsis and rice and more recently brachy-podium,Such models are diploids, have few and small chromosomes, well developed genetics, rapid life cycles, are easily transformed and have extensive sets of technical resources and databases curated by international resource centres. The study of crop genomics today is deeply rooted in earlier studies on model species. Genomes of model species share reasonable genetic synteny with key crop plants which facilitates the discovery of genes and association of genes with phenotypes. While some mechanisms and processes are conserved across the plant kingdom and so can be revealed by studies on any model species,others have diverged during evolution and so are revealed by studying only a closely related model species.Examples of processes that are conserved across the plant kingdom and others that have diverged and therefore need to be understood by studying a more closely related model species are described.

Contrasting patterns of transposable-element insertion polymorphism and nucleotide diversity in autotetraploid and allotetraploid Arabidopsis species

It has been hypothesized that polyploidy permits the proliferation of transposable elements, due to both the masking of deleterious recessive mutations and the breakdown of host silencing mechanisms. We investigated the patterns of insertion polymorphism of an Ac-like transposable element and nucleotide diversity at 18 gene fragments in the allotetraploid Arabidopsis suecica and the autotetraploid A. arenosa. All identified insertions were fixed in A. suecica, and many were clearly inherited from the parental species A. thaliana or A. arenosa. These results are inconsistent with a rapid increase in transposition associated with hybrid breakdown but support the evidence from nucleotide polymorphism patterns of a recent single origin of this species leading to genomewide fixations of transposable elements. In contrast, most insertions were segregating at very low frequencies in A. arenosa samples, showing a significant departure from neutrality in favor of purifying selection, even when we account for population subdivision inferred from sequence variation. Patterns of nucleotide variation at reference genes are consistent with the TE results, showing evidence for higher effective population sizes in A. arenosa than in related diploid taxa but a near complete population bottleneck associated with the origins of A. suecica.

Germline transformation of Shepherd's purse (Capsella bursa-pastoris) by the 'floral dip' method as a tool for evolutionary and developmental biology

Capsella bursa-pastoris is an attractive model system for evolutionary and developmental biology. To facilitate future studies on gene function, the 'floral dip' method was adapted to achieve germline transformation of C. bursa-pastoris. The GFP and BASTA-resistance (BAR (r)) genes were used as markers for screening or selecting, respectively, putative transgenic C. bursa-pastoris plants and the beta-glucuronidase (GUS) gene as well as the GFP gene for monitoring transgene expression level. We tested two Agrobacterium strains, LBA4404 and GV3101, for their ability to transform C. bursa-pastoris. In contrast to Arabidopsis thaliana, for which both strains were able to transform different ecotypes, only GV3101 gave satisfactory transformation rates with C. bursa-pastoris. Furthermore, we evaluated the effects of different concentrations of sucrose and the surfactant Silwet L-77 on the efficiency to generate transgenic C. bursa-pastoris plants and identified an efficient medium containing 10% (w/v) sucrose and 0.02-0.05% (v/v) Silwet L-77. Using Southern hybridisation, we confirmed the integration of the marker gene in the plant genome and the stable heredity of the introduced genes in the next generation.

Transcriptional regulation of NADP-dependent malate dehydrogenase: comparative genetics and identification of DNA-binding proteins

The transcriptional regulation of NADP-malate dehydrogenase (NADP-MDH) was analyzed in Arabidopsis ecotypes and other Brassicaceae. The amount of transcript increased twofold after transfer into low temperature (12 degrees C) or high light (750 microE) in all species. Analysis of the genomic DNA reveals that the NADP-MDH gene (At5g58330 in A. thaliana) in Brassicaceae is located between two other genes (At5g58320 and At5g58340 in Arabidopsis), both encoded on the opposite DNA strand. No promoter elements were identified in 5' direction of the NADP-MDH gene, and the expression of NADP-MDH was not affected in knock-out plants carrying a DNA insert in the 5' region. A yeast-one hybrid approach yielded only three DNA-binding proteins for the 500-bp fragment located upstream of the ATG sequence, but 34 proteins for its coding region. However, in Chlamydomonas and in some Poaceae, which do not possess any genes within the 1200 bp upstream region, typical promoter elements were identified. Alignments of genomic DNA reveal that, in contrast to Poaceae, the introns are highly conserved within Brassicaceae. We conclude that in Brassicaceae the majority of regulatory elements are located within the coding region. The NADP-MDH gene of both families evolved from a common precursor, similar to the gene in Chlamydomonas. Changes in the selection pressure allowed the insertion of At5g58340 into the promoter region of a common ancestor. When the demand for transcriptional regulation increased, At5g58340 disappeared in Poaceae, and a promoter developed in the 5' region. In contrast, Brassicaceae maintained At5g58340 and shifted all regulatory elements into the coding region of NADP-MDH.

Production and genetic analysis of partial hybrids in intertribal crosses between Brassica species (B. rapa, B. napus) and Capsella bursa-pastoris

Capsella bursa-pastoris (L.) Medic (2n = 4x = 32) is a natural double-low (erucic acid < 1%, glucosinolates < 30 micromol/g) germplasm and shows high degree of resistance to Sclerotinia sclerotiorum. Hybridizations were carried out between two Brassica species viz. B. rapa (2n = 20) and B. napus (2n = 38) as female and C. bursa-pastoris as male parent to introduce these desirable traits into cultivated Brassica species. Majority of F(1) plants resembled female parents in morphology and only a few expressed some characters of male parent, including the white petals. Based on cytological observation of somatic cells, the F(1) plants were classified into five types: two types from the cross with B. rapa, type I had 2n = 27-29; type II had 2n = 20; three types from the crosses with B. napus, type III was haploids with 2n = 19; type IV had 2n = 29; type V had 2n = 38. One to two chromosomes of C. bursa-pastoris were detected in pollen mother cells (PMCs) of type I plant by genomic in situ hybridization (GISH), together with chromosomal segments in ovary cells and PMCs of some F1 plants. Amplified fragment length polymorphism (AFLP) bands specific for the male parent, novel for two parents and absent bands in Brassica parents were generated in majority of F1 plants, even in Brassica-types and haploids, indicating the introgressions at various levels from C. bursa-pastoris and genomic alterations following hybridization. Some Brassica-type progeny plants had reduced contents of erucic acid and glucosinolates associated with improved resistance to S. sclerotiorum. The cytological and molecular mechanisms behind these results are discussed.

Differential expression of genes important for adaptation in Capsella bursa-pastoris (Brassicaceae)

Understanding the genetic basis of natural variation is of primary interest for evolutionary studies of adaptation. In Capsella bursa-pastoris, a close relative of Arabidopsis (Arabidopsis thaliana), variation in flowering time is correlated with latitude, suggestive of an adaptation to photoperiod. To identify pathways regulating natural flowering time variation in C. bursa-pastoris, we have studied gene expression differences between two pairs of early- and late-flowering C. bursa-pastoris accessions and compared their response to vernalization. Using Arabidopsis microarrays, we found a large number of significant differences in gene expression between flowering ecotypes. The key flowering time gene FLOWERING LOCUS C (FLC) was not differentially expressed prior to vernalization. This result is in contrast to those in Arabidopsis, where most natural flowering time variation acts through FLC. However, the gibberellin and photoperiodic flowering pathways were significantly enriched for gene expression differences between early- and late-flowering C. bursa-pastoris. Gibberellin biosynthesis genes were down-regulated in late-flowering accessions, whereas circadian core genes in the photoperiodic pathway were differentially expressed between early- and late-flowering accessions. Detailed time-series experiments clearly demonstrated that the diurnal rhythm of CIRCADIAN CLOCK-ASSOCIATED1 (CCA1) and TIMING OF CAB EXPRESSION1 (TOC1) expression differed between flowering ecotypes, both under constant light and long-day conditions. Differential expression of flowering time genes was biologically validated in an independent pair of flowering ecotypes, suggesting a shared genetic basis or parallel evolution of similar regulatory differences. We conclude that genes involved in regulation of the circadian clock, such as CCA1 and TOC1, are strong candidates for the evolution of adaptive flowering time variation in C. bursa-pastoris.

Epigenetic mechanisms for breakdown of self-incompatibility in interspecific hybrids

As a major agent of rapid speciation, interspecific hybridization has played an important role in plant evolution. When hybridization involves species that exhibit self-incompatibility (SI), this prezygotic barrier to self-fertilization must be overcome or lost to allow selfing. How SI, a normally dominant trait, is lost in nascent hybrids is not known, however. Here we demonstrate that hybrid self-fertility can result from epigenetic changes in expression of the S-locus genes that determine specificity in the SI response. We analyzed loss of SI in synthetic hybrids produced by crossing self-fertile and self-incompatible species in each of two crucifer genera. We show that SI is lost in the stigmas of A. thaliana-lyrata hybrids and their neo-allotetraploid derivatives and in the pollen of C. rubella-grandiflora hybrids and their homoploid progenies. Aberrant processing of S-locus receptor kinase gene transcripts as detected in Arabidopsis hybrids and suppression of the S-locus cysteine-rich protein gene as observed in Capsella hybrids are two reversible mechanisms by which SI might break down upon interspecific hybridization to generate self-fertile hybrids in nature.

Evolution of the self-incompatibility system in the Brassicaceae: identification of S-locus receptor kinase (SRK) in self-incompatible Capsella grandiflora

Self-incompatibility (SI) has been well studied in the genera Brassica and Arabidopsis, which have become models for investigation into the SI system. To understand the evolution of the SI system in the Brassicaceae, comparative analyses of the S-locus in genera other than Brassica and Arabidopsis are necessary. We report the identification of six putative S-locus receptor kinase genes (SRK) in natural populations of Capsella grandiflora, an SI species from a genus which is closely related to Arabidopsis. These S-alleles display striking similarities to the Arabidopsis lyrata SRK alleles in sequence and structure. Our phylogenetic analysis supports the scenario of differing SI evolution along the two lineages (The Brassica lineage and Arabidopsis/Capsella lineage). Our results also argue that the ancestral S-locus lacked the SLG gene (S-locus glycoprotein) and that the diversification of S-alleles predates the separation of Arabidopsis and Capsella.

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.

Dynamic evolution at pericentromeres

Pericentromeres are exceptional genomic regions: in animals they contain extensive segmental duplications implicated in gene creation, and in plants they sustain rearrangements and insertions uncommon in euchromatin. To examine the mechanisms and patterns of plant pericentromere evolution, we compared pericentromere sequence from four Brassicaceae species separated by <15 million years (Myr). This flowering plant family is ideal for studying relationships between genome reorganization and pericentromere evolution-its members have undergone recent polyploidization and hybridization, with close relatives changing in genome size and chromosome number. Through sequence and hybridization analyses, we examined regions from Arabidopsis arenosa, Capsella rubella, and Olimarabidopsis pumila that are homologous to Arabidopsis thaliana pericentromeres (peri-CENs) III and V, and used FISH to demonstrate they have been maintained near centromere satellite arrays in each species. Sequence analysis revealed a set of highly conserved genes, yet we discovered substantial differences in intergenic length and species-specific changes in sequence content and gene density. We discovered that A. thaliana has undergone recent, significant expansions within its pericentromeres, in some cases measuring hundreds of kilobases; these findings are in marked contrast to euchromatic segments in these species that exhibit only minor length changes. While plant pericentromeres do contain some duplications, we did not find evidence of extensive segmental duplications, as has been documented in primates. Our data support a model in which plant pericentromeres may experience selective pressures distinct from euchromatin, tolerating rapid, dynamic changes in structure and sequence content, including large insertions of mobile elements, 5S rDNA arrays and pseudogenes.

Construction of a genetic linkage map of Thlaspi caerulescens and quantitative trait loci analysis of zinc accumulation

Zinc (Zn) hyperaccumulation seems to be a constitutive species-level trait in Thlaspi caerulescens. When compared under conditions of equal Zn availability, considerable variation in the degree of hyperaccumulation is observed among accessions originating from different soil types. This variation offers an excellent opportunity for further dissection of the genetics of this trait. A T. caerulescens intraspecific cross was made between a plant from a nonmetallicolous accession [Lellingen (LE)], characterized by relatively high Zn accumulation, and a plant from a calamine accession [La Calamine (LC)], characterized by relatively low Zn accumulation. Zinc accumulation in roots and shoots segregated in the F3 population. This population was used to construct an LE/LC amplified fragment length polymorphism (AFLP)-based genetic linkage map and to map quantitative trait loci (QTL) for Zn accumulation. Two QTL were identified for root Zn accumulation, with the trait-enhancing alleles being derived from each of the parents, and explaining 21.7 and 16.6% of the phenotypic variation observed in the mapping population. Future development of more markers, based on Arabidopsis orthologous genes localized in the QTL regions, will allow fine-mapping and map-based cloning of the genes underlying the QTL.

Comparative genome analyses of Arabidopsis spp.: inferring chromosomal rearrangement events in the evolutionary history of A. thaliana

Comparative genome analysis is a powerful tool that can facilitate the reconstruction of the evolutionary history of the genomes of modern-day species. The model plant Arabidopsis thaliana with its n = 5 genome is thought to be derived from an ancestral n = 8 genome. Pairwise comparative genome analyses of A. thaliana with polyploid and diploid Brassicaceae species have suggested that rapid genome evolution, manifested by chromosomal rearrangements and duplications, characterizes the polyploid, but not the diploid, lineages of this family. In this study, we constructed a low-density genetic linkage map of Arabidopsis lyrata ssp. lyrata (A. l. lyrata; n = 8, diploid), the closest known relative of A. thaliana (MRCA approximately 5 Mya), using A. thaliana-specific markers that resolve into the expected eight linkage groups. We then performed comparative Bayesian analyses using raw mapping data from this study and from a Capsella study to infer the number and nature of rearrangements that distinguish the n = 8 genomes of A. l. lyrata and Capsella from the n = 5 genome of A. thaliana. We conclude that there is strong statistical support in favor of the parsimony scenarios of 10 major chromosomal rearrangements separating these n = 8 genomes from A. thaliana. These chromosomal rearrangement events contribute to a rate of chromosomal evolution higher than previously reported in this lineage. We infer that at least seven of these events, common to both sets of data, are responsible for the change in karyotype and underlie genome reduction in A. thaliana.

Comparing the linkage maps of the close relatives Arabidopsis lyrata and A. thaliana

We have constructed a genetic map of Arabidopsis lyrata, a self-incompatible relative of the plant model species A. thaliana. A. lyrata is a diploid (n = 8) species that diverged from A. thaliana (n = 5) approximately 5 MYA. Mapping was conducted in a full-sib progeny of two unrelated F(1) hybrids between two European populations of A. lyrata ssp. petraea. We used the least-squares method of the Joinmap program for map construction. The gross chromosomal differences between the two species were most parsimoniously explained with three fusions, two reciprocal translocations, and one inversion. The total map length was 515 cM, and the distances were 12% larger than those between corresponding markers in the linkage map of A. thaliana. The 72 markers, consisting of microsatellites and gene-based markers, were spaced on average every 8 cM. Transmission ratio distortion was extensive, and most distortions were specific to each reciprocal cross, suggesting cytoplasmic interactions. We estimate locations and most probable genotype frequencies of transmission ratio distorting loci (TRDL) with a Bayesian method and discuss the possible reasons for the observed distortions.

Comparative mapping between Medicago sativa and Pisum sativum

Comparative genome analysis has been performed between alfalfa (Medicago sativa) and pea (Pisum sativum), species which represent two closely related tribes of the subfamily Papilionoideae with different basic chromosome numbers. The positions of genes on the most recent linkage map of diploid alfalfa were compared to those of homologous loci on the combined genetic map of pea to analyze the degree of co-linearity between their linkage groups. In addition to using unique genes, analysis of the map positions of multicopy (homologous) genes identified syntenic homologs (characterized by similar positions on the maps) and pinpointed the positions of non-syntenic homologs. The comparison revealed extensive conservation of gene order between alfalfa and pea. However, genetic rearrangements (due to breakage and reunion) were localized which can account for the difference in chromosome number (8 for alfalfa and 7 for pea). Based on these genetic events and our increasing knowledge of the genomic structure of pea, it was concluded that the difference in genome size between the two species (the pea genome is 5- to 10-fold larger than that of alfalfa) is not a consequence of genome duplication in pea. The high degree of synteny observed between pea and Medicago loci makes further map-based cloning of pea genes based on the genome resources now available for M. truncatula a promising strategy.

The Iccare web server: an attempt to merge sequence and mapping information for plant and animal species

The Iccare web server, http://genopole.toulouse.inra.fr/bioinfo/Iccare, provides a simple yet efficient tool for crude EST (expressed sequence tag) annotation specifically dedicated to comparative mapping approaches. Iccare uses all the EST and mRNA sequences from public databases for an organism of interest (query species) and compares them to all the transcripts of one reference organism (Homo sapiens or Arabidopsis thaliana). The results are displayed according to the location of the genes on the chromosomes of the reference organism. Gene structure information and sequence similarities are combined in a graphical representation in order to pinpoint the nature of the transcript query sequence. The user can subsequently design primers or probes for the purpose of physical or genetic mapping. In addition to the query organisms already available in Iccare, users can perform a tailor-made search with their own sequences against the animal or plant reference organism genes.

Recent progress in chromosome painting of Arabidopsis and related species

This paper reports on state-of-the-art achievements of chromosome painting in Arabidopsis thaliana (2n = 10). Arabidopsis chromosomes 1, 2 and 4 were painted using chromosome-specific BAC contigs. We consider technical aspects of the painting approach and document major applications, such as the tracing of Arabidopsis chromosomes as interphase chromosome territories and during mitotic and meiotic cell cycles as well as comparative chromosome painting in related species. This is the first report of successful interspecific chromosome painting in plants. The evolutionary history of chromosomes homeologous to Arabidopsis chromosome 4 was reconstructed by hybridization of chromosome-4-specific painting probes to karyotypes of Brassicaceae species with x = 8 chromosomes. Future perspectives of chromosome painting in A. thaliana and its wild relatives are outlined.

Duplicate marker loci can result in incorrect locus orders on linkage maps

Genetic linkage maps, constructed from multi-locus recombination data, are the basis for many applications of molecular markers. For the successful employment of a linkage map, it is essential that the linear order of loci on a chromosome is correct. The objectives of this theoretical study were to (1) investigate the occurrence of incorrect locus orders caused by duplicate marker loci, (2) develop a statistical test for the detection of duplicate markers, and (3) discuss the implications for practical applications of linkage maps. We derived conditions, under which incorrect locus orders do or do not occur with duplicate marker loci for the general case of n markers on a chromosome in a BC(1) mapping population. We further illustrated these conditions numerically for the special case of four markers. On the basis of the extent of segregation distortion, an exact test for the presence of duplicate marker loci was suggested and its power was investigated numerically. Incorrect locus orders caused by duplicate marker loci can (1) negatively affect the assignment of target genes to chromosome regions in a map-based cloning experiment, (2) hinder indirect selection for a favorable allele at a quantitative trait locus, and (3) decrease the efficiency of reducing the length of the chromosome segment attached to a target gene in marker-assisted backcrossing.

Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution

To study the evolutionary effects of polyploidy on plant gene functions, we analyzed functional genomics data for a large number of duplicated gene pairs formed by ancient polyploidy events in Arabidopsis thaliana. Genes retained in duplicate are not distributed evenly among Gene Ontology or Munich Information Center for Protein Sequences functional categories, which indicates a nonrandom process of gene loss. Genes involved in signal transduction and transcription have been preferentially retained, and those involved in DNA repair have been preferentially lost. Although the two members of each gene pair must originally have had identical transcription profiles, less than half of the pairs formed by the most recent polyploidy event still retain significantly correlated profiles. We identified several cases where groups of duplicated gene pairs have diverged in concert, forming two parallel networks, each containing one member of each gene pair. In these cases, the expression of each gene is strongly correlated with the other nonhomologous genes in its network but poorly correlated with its paralog in the other network. We also find that the rate of protein sequence evolution has been significantly asymmetric in >20% of duplicate pairs. Together, these results suggest that functional diversification of the surviving duplicated genes is a major feature of the long-term evolution of polyploids.

The Arabidopsis genome sequence as a tool for genome analysis in Brassicaceae. A comparison of the Arabidopsis and Capsella rubella genomes

The annotated Arabidopsis genome sequence was exploited as a tool for carrying out comparative analyses of the Arabidopsis and Capsella rubella genomes. Comparison of a set of random, short C. rubella sequences with the corresponding sequences in Arabidopsis revealed that aligned protein-coding exon sequences differ from aligned intron or intergenic sequences in respect to the degree of sequence identity and the frequency of small insertions/deletions. Molecular-mapped markers and expressed sequence tags derived from Arabidopsis were used for genetic mapping in a population derived from an interspecific cross between Capsella grandiflora and C. rubella. The resulting eight Capsella linkage groups were compared to the sequence maps of the five Arabidopsis chromosomes. Fourteen colinear segments spanning approximately 85% of the Arabidopsis chromosome sequence maps and 92% of the Capsella genetic linkage map were detected. Several fusions and fissions of chromosomal segments as well as large inversions account for the observed arrangement of the 14 colinear blocks in the analyzed genomes. In addition, evidence for small-scale deviations from genome colinearity was found. Colinearity between the Arabidopsis and Capsella genomes is more pronounced than has been previously reported for comparisons between Arabidopsis and different Brassica species.

Comparison of a Brassica oleracea genetic map with the genome of Arabidopsis thaliana

Brassica oleracea is closely related to the model plant, Arabidopsis thaliana. Despite this relationship, it has been difficult to both identify the most closely related segments between the genomes and determine the degree of genome replication within B. oleracea relative to A. thaliana. These difficulties have arisen in part because both species have replicated genomes, and the criteria used to identify orthologous regions between the genomes are often ambiguous. In this report, we compare the positions of sequenced Brassica loci with a known position on a B. oleracea genetic map to the positions of their putative orthologs within the A. thaliana genome. We use explicit criteria to distinguish orthologous from paralogous loci. In addition, we develop a conservative algorithm to identify collinear loci between the genomes and a permutation test to evaluate the significance of these regions. The algorithm identified 34 significant A. thaliana regions that are collinear with >28% of the B. oleracea genetic map. These regions have a mean of 3.3 markers spanning 2.1 Mbp of the A. thaliana genome and 2.5 cM of the B. oleracea genetic map. Our findings are consistent with the hypothesis that the B. oleracea genome has been highly rearranged since divergence from A. thaliana, likely as a result of polyploidization.

Syntenic relationships between Medicago truncatula and Arabidopsis reveal extensive divergence of genome organization

Arabidopsis and Medicago truncatula represent sister clades within the dicot subclass Rosidae. We used genetic map-based and bacterial artificial chromosome sequence-based approaches to estimate the level of synteny between the genomes of these model plant species. Mapping of 82 tentative orthologous gene pairs reveals a lack of extended macrosynteny between the two genomes, although marker collinearity is frequently observed over small genetic intervals. Divergence estimates based on non-synonymous nucleotide substitutions suggest that a majority of the genes under analysis have experienced duplication in Arabidopsis subsequent to divergence of the two genomes, potentially confounding synteny analysis. Moreover, in cases of localized synteny, genetically linked loci in M. truncatula often share multiple points of synteny with Arabidopsis; this latter observation is consistent with the large number of segmental duplications that compose the Arabidopsis genome. More detailed analysis, based on complete sequencing and annotation of three M. truncatula bacterial artificial chromosome contigs suggests that the two genomes are related by networks of microsynteny that are often highly degenerate. In some cases, the erosion of microsynteny could be ascribed to the selective gene loss from duplicated loci, whereas in other cases, it is due to the absence of close homologs of M. truncatula genes in Arabidopsis.

Sequence analysis of a 282-kilobase region surrounding the citrus Tristeza virus resistance gene (Ctv) locus in Poncirus trifoliata L. Raf

Citrus tristeza virus (CTV) is the major virus pathogen causing significant economic damage to citrus worldwide, and a single dominant gene, Ctv, provides broad spectrum resistance to CTV in Poncirus trifoliata L. Raf. Ctv was physically mapped to a 282-kb region using a P. trifoliata bacterial artificial chromosome library. This region was completely sequenced to about 8x coverage using a shotgun sequencing strategy and primer walking for gap closure. Sequence analysis predicts 22 putative genes, two mutator-like transposons and eight retrotransposons. This sequence analysis also revealed some interesting features of this region of the P. trifoliata genome: a disease resistance gene cluster with seven members and eight retrotransposons clustered in a 125-kb gene-poor region. Comparative sequence analysis suggests that six genes in the Ctv region have significant sequence similarity with their orthologs in bacterial artificial chromosome clones F7H2 and F21T11 from Arabidopsis chromosome I. However, the analysis of gene colinearity between P. trifoliata and Arabidopsis indicates that Arabidopsis genome sequence information may be of limited use for positional gene cloning in P. trifoliata and citrus. Analysis of candidate genes for Ctv is also discussed.

A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome

The Arabidopsis genome contains numerous large duplicated chromosomal segments, but the different approaches used in previous analyses led to different interpretations regarding the number and timing of ancestral large-scale duplication events. Here, using more appropriate methodology and a more recent version of the genome sequence annotation, we investigate the scale and timing of segmental duplications in Arabidopsis. We used protein sequence similarity searches to detect duplicated blocks in the genome, used the level of synonymous substitution between duplicated genes to estimate the relative ages of the blocks containing them, and analyzed the degree of overlap between adjacent duplicated blocks. We conclude that the Arabidopsis lineage underwent at least two distinct episodes of duplication. One was a polyploidy that occurred much more recently than estimated previously, before the Arabidopsis/Brassica rapa split and probably during the early emergence of the crucifer family (24-40 Mya). An older set of duplicated blocks was formed after the monocot/dicot divergence, and the relatively low level of overlap among these blocks indicates that at least some of them are remnants of a larger duplication such as a polyploidy or aneuploidy.

Rates and patterns of molecular evolution in inbred and outbred Arabidopsis

The evolution of self-fertilization is associated with a large reduction in the effective rate of recombination and a corresponding decline in effective population size. If many spontaneous mutations are slightly deleterious, this shift in the breeding system is expected to lead to a reduced efficacy of natural selection and genome-wide changes in the rates of molecular evolution. Here, we investigate the effects of the breeding system on molecular evolution in the highly self-fertilizing plant Arabidopsis thaliana by comparing its coding and noncoding genomic regions with those of its close outcrossing relative, the self-incompatible A. lyrata. More distantly related species in the Brassicaceae are used as outgroups to polarize the substitutions along each lineage. In contrast to expectations, no significant difference in the rates of protein evolution is observed between selfing and outcrossing Arabidopsis species. Similarly, no consistent overall difference in codon bias is observed between the species, although for low-biased genes A. lyrata shows significantly higher major codon usage. There is also evidence of intron size evolution in A. thaliana, which has consistently smaller introns than its outcrossing congener, potentially reflecting directional selection on intron size. The results are discussed in the context of heterogeneity in selection coefficients across loci and the effects of life history and population structure on rates of molecular evolution. Using estimates of substitution rates in coding regions and approximate estimates of divergence and generation times, the genomic deleterious mutation rate (U) for amino acid substitutions in Arabidopsis is estimated to be approximately 0.2-0.6 per generation.

Applied plant genomics: the secret is integration

Although concerted efforts to understand selected botanical models have been made, the resulting basic knowledge varies in its applicability to other diverse species including the major crops. Recent advances in high-throughput genomics are offering new avenues through which to exploit model systems for the study of botanical diversity, providing prospects for crop improvement. In particular, whole-genome sequencing has provided opportunities for the broader application of reverse genetics, expression profiling, and molecular mapping in diverse species.

Comparative genomics between rice and Arabidopsis shows scant collinearity in gene order

We have investigated possible collinearity between the genomes of rice and Arabidopsis by comparing 126 annotated and mapped rice BAC sequences (approximately 20 Mb of sequence) with the annotated and complete Arabidopsis genome (approximately 115 Mb). Although we were able to identify several regions in which gene order is preserved, they are relatively small, and are interrupted by noncollinear genes. Computer simulation showed that these microscale collinearities are above the expectation for a random process. On the other hand, the order of exons within homologous genes (<2.5 kb) was preserved, as expected.

Contrasting genome organisation: two regions of the Brassica oleracea genome compared with collinear regions of the Arabidopsis thaliana genome

Brassica crop species are of worldwide importance and are closely related to the model plant Arabidopsis thaliana for which the complete genome sequence has recently been established. We investigated collinearity of marker order by comparing two contrasting regions of the Brassica oleracea genome with homologous regions of A. thaliana. Although there is widespread replication of marker loci in both A. thaliana and B. oleracea, we found that a combination of genetic markers mapped in B. oleracea, including RFLPs, CAPS, and SSRs allowed comparison and interpretation of medium-scale chromosomal organisation and rearrangements. The interpretation of data was facilitated by hybridising probes onto the whole A. thaliana genome, as represented by BAC contigs. Twenty marker loci were sampled from the whole length of the shortest B. oleracea linkage group, O6, and 21 from a 30.4-cM section of the longest linkage group, O3. There is evidence of locus duplication on linkage group O6. Locus order is well conserved between a putative duplicated region of 10.5 cM and a discrete region comprising 25 cM of A. thaliana chromosome I. This was supported by evidence from seven paralogous loci, three of which were duplicated in a 30.6-cM region of linkage group O6. The pattern of locus order for the remainder of linkage group O6 and the sampled section of linkage group O3 was more complex when compared with the A. thaliana genome. Although there was some conservation of locus order between markers on linkage group O3 and approximately 9 cM of A. thaliana chromosome I, this was superimposed upon a complex pattern of additional loci that were replicated in both A. thaliana and B. oleracea. The results are discussed in the context of the ability to use collinear information to assist map-based cloning.Key words: comparative mapping, BAC, physical contig, MADS box.

Multiple losses and transfers to the nucleus of two mitochondrial succinate dehydrogenase genes during angiosperm evolution

Unlike in animals, the functional transfer of mitochondrial genes to the nucleus is an ongoing process in plants. All but one of the previously reported transfers in angiosperms involve ribosomal protein genes. Here we report frequent transfer of two respiratory genes, sdh3 and sdh4 (encoding subunits 3 and 4 of succinate dehydrogenase), and we also show that these genes are present and expressed in the mitochondria of diverse angiosperms. Southern hybridization surveys reveal that sdh3 and sdh4 have been lost from the mitochondrion about 40 and 19 times, respectively, among the 280 angiosperm genera examined. Transferred, functional copies of sdh3 and sdh4 were characterized from the nucleus in four and three angiosperm families, respectively. The mitochondrial targeting presequences of two sdh3 genes are derived from preexisting genes for anciently transferred mitochondrial proteins. On the basis of the unique presequences of the nuclear genes and the recent mitochondrial gene losses, we infer that each of the seven nuclear sdh3 and sdh4 genes was derived from a separate transfer to the nucleus. These results strengthen the hypothesis that angiosperms are experiencing a recent evolutionary surge of mitochondrial gene transfer to the nucleus and reveal that this surge includes certain respiratory genes in addition to ribosomal protein genes.

Genome sequences and evolutionary biology, a two-way interaction

Complete genome sequences are accumulating rapidly, culminating with the announcement of the human genome sequence in February 2001. In addition to cataloguing the diversity of genes and other sequences, genome sequences will provide the first detailed and complete data on gene families and genome organization, including data on evolutionary changes. Reciprocally, evolutionary biology will make important contributions to the efforts to understand functions of genes and other sequences in genomes. Large-scale, detailed and unbiased comparisons between species will illuminate the evolution of genes and genomes, and population genetics methods will enable detection of functionally important genes or sequences, including sequences that have been involved in adaptive changes.

Duplicate and diverge: the evolution of plant genome microstructure

The use of positional approaches for the isolation of genes from most crop species is difficult due to the large size of their genomes. If the order of genes in segments of the genomes is similar in different plants, it might be feasible to use smaller genomes as templates upon which to base strategies for the positional cloning of genes from other species. Comparative genetic mapping, using markers such as restriction-fragment length polymorphisms, has revealed extensive conservation of long-range genome organization (macrostructure) between related species. But is the organization of the tens or hundreds of genes between the genetic markers also conserved? Recent results suggest that the fine-scale structure (microstructure) of plant genomes is more dynamic than previously assumed from investigations of the macrostructure.

Arabidopsis genome sequence as a tool for functional genomics in tomato

Tomato is a well-established model organism for studying many biological processes including resistance and susceptibility to pathogens and the development and ripening of fleshy fruits. The availability of the complete Arabidopsis genome sequence will expedite map-based cloning in tomato on the basis of chromosomal synteny between the two species, and will facilitate the functional analysis of tomato genes.

Analysis of the genome sequence of the flowering plant Arabidopsis thaliana

The flowering plant Arabidopsis thaliana is an important model system for identifying genes and determining their functions. Here we report the analysis of the genomic sequence of Arabidopsis. The sequenced regions cover 115.4 megabases of the 125-megabase genome and extend into centromeric regions. The evolution of Arabidopsis involved a whole-genome duplication, followed by subsequent gene loss and extensive local gene duplications, giving rise to a dynamic genome enriched by lateral gene transfer from a cyanobacterial-like ancestor of the plastid. The genome contains 25,498 genes encoding proteins from 11,000 families, similar to the functional diversity of Drosophila and Caenorhabditis elegans--the other sequenced multicellular eukaryotes. Arabidopsis has many families of new proteins but also lacks several common protein families, indicating that the sets of common proteins have undergone differential expansion and contraction in the three multicellular eukaryotes. This is the first complete genome sequence of a plant and provides the foundations for more comprehensive comparison of conserved processes in all eukaryotes, identifying a wide range of plant-specific gene functions and establishing rapid systematic ways to identify genes for crop improvement.