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

Chromosomal localization of rDNA genes and genomic organization of 5S rDNA in Oreochromis mossambicus, O. urolepis hornorum and their hybrid

In this study, classical and molecular cytogenetic analyses were performed in tilapia fishes, Oreochromis mossambicus (XX/XY sex determination system), O. urolepis hornorum (WZ/ZZ sex determination system) and their hybrid by crossing O. mossambicus female x O. u. hornorum male. An identical karyotype ((2n = 44, NF (total number of chromosomal arms) = 50) was obtained from three examined tilapia samples. Genomic organization analysis of 5S rDNA revealed two different types of 5S rDNA sequences, 5S type I and 5S type II. Moreover, fluorescence in situ hybridization (FISH) with 5S rDNA probes showed six positive fluorescence signals on six chromosomes of all the analysed metaphases from the three tilapia samples. Subsequently, 45S rDNA probes were also prepared, and six positive fluorescence signals were observed on three chromosome pairs in all analysed metaphases of the three tilapia samples. The correlation between 45 rDNA localization and nucleolar organizer regions (NORs) was confirmed by silver nitrate staining in tilapia fishes. Further, different chromosomal localizations of 5S rDNA and 45S rDNA were verified by two different colour FISH probes. Briefly, the current data provide an insights for hybridization projects and breeding improvement of tilapias.

Chromosomal evolution in the plant family Solanaceae

BACKGROUND

Over the past decades, extensive comparative mapping research has been performed in the plant family Solanaceae. The recent identification of a large set of single-copy conserved orthologous (COSII) markers has greatly accelerated comparative mapping studies among major solanaceous species including tomato, potato, eggplant, pepper and diploid Nicotiana species (as well as tetraploid tobacco). The large amount of comparative data now available for these species provides the opportunity to describe the overall patterns of chromosomal evolution in this important plant family. The results of this investigation are described herein.

RESULTS

We combined data from multiple COSII studies, and other comparative mapping studies performed in tomato, potato, eggplant, pepper and diploid Nicotiana species, to deduce the features and outcomes of chromosomal evolution in the Solanaceae over the past 30 million years. This includes estimating the rates and timing of chromosomal changes (inversions and translocations) as well as deducing the age of ancestral progenitor species and predicting their genome configurations.

CONCLUSIONS

The Solanaceae has experienced chromosomal changes at a modest rate compared with other families and the rates are likely conserved across different lineages of the family. Chromosomal inversions occur at a consistently higher rate than do translocations. Further, we find evidences for non-random positioning of the chromosomal rearrangement breakpoints. This finding is consistent with the similar finding in mammals, where hot spots for chromosomal breakages have apparently played a significant role in shaping genome evolution. Finally, by utilizing multiple genome comparisons we were able to reconstruct the most likely genome configuration for a number of now-extinct progenitor species that gave rise to the extant solanaceous species used in this research. The results from this study provide the first broad overview of chromosomal evolution in the family Solanaceae, and one of the most detailed thus far for any family of plants.

Chromosomal evolution in the plant family Solanaceae

Background

Over the past decades, extensive comparative mapping research has been performed in the plant family Solanaceae. The recent identification of a large set of single-copy conserved orthologous (COSII) markers has greatly accelerated comparative mapping studies among major solanaceous species including tomato, potato, eggplant, pepper and diploid Nicotiana species (as well as tetraploid tobacco). The large amount of comparative data now available for these species provides the opportunity to describe the overall patterns of chromosomal evolution in this important plant family. The results of this investigation are described herein.

Results

We combined data from multiple COSII studies, and other comparative mapping studies performed in tomato, potato, eggplant, pepper and diploid Nicotiana species, to deduce the features and outcomes of chromosomal evolution in the Solanaceae over the past 30 million years. This includes estimating the rates and timing of chromosomal changes (inversions and translocations) as well as deducing the age of ancestral progenitor species and predicting their genome configurations.

Conclusions

The Solanaceae has experienced chromosomal changes at a modest rate compared with other families and the rates are likely conserved across different lineages of the family. Chromosomal inversions occur at a consistently higher rate than do translocations. Further, we find evidences for non-random positioning of the chromosomal rearrangement breakpoints. This finding is consistent with the similar finding in mammals, where hot spots for chromosomal breakages have apparently played a significant role in shaping genome evolution. Finally, by utilizing multiple genome comparisons we were able to reconstruct the most likely genome configuration for a number of now-extinct progenitor species that gave rise to the extant solanaceous species used in this research. The results from this study provide the first broad overview of chromosomal evolution in the family Solanaceae, and one of the most detailed thus far for any family of plants.

Distinct, ecotype-specific genome and proteome signatures in the marine cyanobacteria Prochlorococcus

BACKGROUND

The marine cyanobacterium Prochlorococcus marinus, having multiple ecotypes of distinct genotypic/phenotypic traits and being the first documented example of genome shrinkage in free-living organisms, offers an ideal system for studying niche-driven molecular micro-diversity in closely related microbes. The present study, through an extensive comparative analysis of various genomic/proteomic features of 6 high light (HL) and 6 low light (LL) adapted strains, makes an attempt to identify molecular determinants associated with their vertical niche partitioning.

RESULTS

Pronounced strand-specific asymmetry in synonymous codon usage is observed exclusively in LL strains. Distinct dinucleotide abundance profiles are exhibited by 2 LL strains with larger genomes and G+C-content approximately 50% (group LLa), 4 LL strains having reduced genomes and G+C-content approximately 35-37% (group LLb), and 6 HL strains. Taking into account the emergence of LLa, LLb and HL strains (based on 16S rRNA phylogeny), a gradual increase in average aromaticity, pI values and beta- & coil-forming propensities and a decrease in mean hydrophobicity, instability indices and helix-forming propensities of core proteins are observed. Greater variations in orthologous gene repertoire are found between LLa and LLb strains, while higher number of positively selected genes exist between LL and HL strains.

CONCLUSION

Strains of different Prochlorococcus groups are characterized by distinct compositional, physicochemical and structural traits that are not mere remnants of a continuous genetic drift, but are potential outcomes of a grand scheme of niche-oriented stepwise diversification, that might have driven them chronologically towards greater stability/fidelity and invoked upon them a special ability to inhabit diverse oceanic environments.

Comparative analysis of the Brassica oleracea genetic map and the Arabidopsis thaliana genome

We further investigated genome macrosynteny for Brassica species and Arabidopsis thaliana. This work aimed at comparative map construction for B. oleracea and A. thaliana chromosomes based on 160 known A. thaliana probes: 147 expressed sequence tags (ESTs) and 13 full-length cDNA clones. Based on an in silico study of the A. thaliana genome, most of the selected ESTs (83%) represented unique or low-copy genes. We identified conserved segments by the visual inspection of comparative data with a priori assumptions, and established their significance with the LineUp algorithm. Evaluation of the number of B. oleracea gene copies per A. thaliana EST revealed a fixed upward trend. We established a segregation distortion pattern for all genetic loci, with particular consideration of the type of selection (gametic or zygotic), and discuss its possible impact on genetic map construction. Consistent with previous reports, we found evidence for numerous chromosome rearrangements and the genome fragment replication of B. oleracea that have taken place since the divergence of the two species. Also, we found that over 54% of the B. oleracea genome is covered by 24 segments conserved with the A. thaliana genome. The average conserved segment is composed of 5 loci covering 19.3 cM in the B. oleracea genetic map and 2.42 Mb in the A. thaliana physical map. We have also attempted to use a unified system of conserved blocks (previously described) to verify our results and perform a comprehensive comparison with other Brassica species.

Contrasting evolutionary dynamics between angiosperm and mammalian genomes

Continuing advances in genomics are revealing substantial differences between genomes of major eukaryotic lineages. Because most data (in terms of depth and phylogenetic breadth) are available for angiosperms and mammals, we explore differences between these groups and show that angiosperms have less highly compartmentalized and more diverse genomes than mammals. In considering the causes of these differences, four mechanisms are highlighted: polyploidy, recombination, retrotransposition and genome silencing, which have different modes and time scales of activity. Angiosperm genomes are evolutionarily more dynamic and labile, whereas mammalian genomes are more stable at both the sequence and chromosome level. We suggest that fundamentally different life strategies and development feedback on the genome exist, influencing dynamics and evolutionary trajectories at all levels from the gene to the genome.

Footprints of inversions at present and past pseudoautosomal boundaries in human sex chromosomes

The human sex chromosomes have stopped recombining gradually, which has left five evolutionary strata on the X chromosome. Y inversions are thought to have suppressed X–Y recombination but clear evidence is missing. Here, we looked for such evidence by focusing on a region—the X-added region (XAR)—that includes the pseudoautosomal region and the most recent strata 3 to 5. We estimated and analyzed the whole set of parsimonious scenarios of Y inversions given the gene order in XAR and its Y homolog. Comparing these to scenarios for simulated sequences suggests that the strata 4 and 5 were formed by Y inversions. By comparing the X and Y DNA sequences, we found clear evidence of two Y inversions associated with duplications that coincide with the boundaries of strata 4 and 5. Divergence between duplicates is in agreement with the timing of strata 4 and 5 formation. These duplicates show a complex pattern of gene conversion that resembles the pattern previously found for AMELXY, a stratum 3 locus. This suggests that this locus—despite AMELY being unbroken—was possibly involved in a Y inversion that formed stratum 3. However, no clear evidence supporting the formation of stratum 3 by a Y inversion was found, probably because this stratum is too old for such an inversion to be detectable. Our results strongly support the view that the most recent human strata have arisen by Y inversions and suggest that inversions have played a major role in the differentiation of our sex chromosomes.

Genome sequence comparison of Col and Ler lines reveals the dynamic nature of Arabidopsis chromosomes

Large differences in plant genome sizes are mainly due to numerous events of insertions or deletions (indels). The balance between these events determines the evolutionary direction of genome changes. To address the question of what phenomena trigger these alterations, we compared the genomic sequences of two Arabidopsis thaliana lines, Columbia (Col) and Landsberg erecta (Ler). Based on the resulting alignments large indels (>100 bp) within these two genomes were analysed. There are ∼8500 large indels accounting for the differences between the two genomes. The genetic basis of their origin was distinguished as three main categories: unequal recombination (Urec)-derived, illegitimate recombination (Illrec)-derived and transposable elements (TE)-derived. A detailed study of their distribution and size variation along chromosomes, together with a correlation analyses, allowed us to demonstrate the impact of particular recombination-based mechanisms on the plant genome evolution. The results show that unequal recombination is not efficient in the removal of TEs within the pericentromeric regions. Moreover, we discovered an unexpectedly high influence of large indels on gene evolution pointing out significant differences between the various gene families. For the first time, we present convincing evidence that somatic events do play an important role in plant genome evolution.

The allotetraploid Arabidopsis thaliana-Arabidopsis lyrata subsp. petraea as an alternative model system for the study of polyploidy in plants

Polyploidy is known to be common in plants and recent work has focused on the rapid changes in genome structure and expression that occur upon polyploidization. In Arabidopsis, much of this work has been done on a synthetic allotetraploid obtained by crossing a tetraploid Arabidopsis thaliana (2n = 4x = 20) with A. arenosa (2n = 4x = 32). To explore an alternative route to polyploidy in this model species, we have developed a synthetic allopolyploid by crossing two diploid species: A. thaliana (2n = 2x = 10) and Arabidopsis lyrata subsp. petraea (2n = 2x = 16). F(1) hybrids were easy to obtain and phenotypically more similar to A. lyrata. Spontaneous chromosome doubling events occurred in about 25% of the F(1)s, thus restoring fertility. The resulting allotetraploids (2n = 26) exhibited many genomic changes typically reported upon polyploidization. Nucleolar dominance was observed as only the A. lyrata rDNA loci were expressed in the F(1) and allotetraploids. Changes in the degree of methylation were observed at almost 25% of the loci examined by MSAP analysis. Finally, structural genomic alterations did occur as a large deletion covering a significant portion of the upper arm of chromosome II was detected but no evidence of increased mobility of transposons was obtained. Such allotetraploids derived from two parents with sequenced (or soon to be sequenced) genomes offer much promise in elucidating the various changes that occur in newly synthesized polyploids.

Polyploidy and angiosperm diversification

Polyploidy has long been recognized as a major force in angiosperm evolution. Recent genomic investigations not only indicate that polyploidy is ubiquitous among angiosperms, but also suggest several ancient genome-doubling events. These include ancient whole genome duplication (WGD) events in basal angiosperm lineages, as well as a proposed paleohexaploid event that may have occurred close to the eudicot divergence. However, there is currently no evidence for WGD in Amborella, the putative sister species to other extant angiosperms. The question is no longer "What proportion of angiosperms are polyploid?", but "How many episodes of polyploidy characterize any given lineage?" New algorithms provide promise that ancestral genomes can be reconstructed for deep divergences (e.g., it may be possible to reconstruct the ancestral eudicot or even the ancestral angiosperm genome). Comparisons of diversification rates suggest that genome doubling may have led to a dramatic increase in species richness in several angiosperm lineages, including Poaceae, Solanaceae, Fabaceae, and Brassicaceae. However, additional genomic studies are needed to pinpoint the exact phylogenetic placement of the ancient polyploidy events within these lineages and to determine when novel genes resulting from polyploidy have enabled adaptive radiations.

Multilocus patterns of nucleotide diversity, population structure and linkage disequilibrium in Boechera stricta, a wild relative of Arabidopsis

Information about polymorphism, population structure, and linkage disequilibrium (LD) is crucial for association studies of complex trait variation. However, most genomewide studies have focused on model systems, with very few analyses of undisturbed natural populations. Here, we sequenced 86 mapped nuclear loci for a sample of 46 genotypes of Boechera stricta and two individuals of B. holboellii, both wild relatives of Arabidopsis. Isolation by distance was significant across the species range of B. stricta, and three geographic groups were identified by structure analysis, principal coordinates analysis, and distance-based phylogeny analyses. The allele frequency spectrum indicated a genomewide deviation from an equilibrium neutral model, with silent nucleotide diversity averaging 0.004. LD decayed rapidly, declining to background levels in approximately 10 kb or less. For tightly linked SNPs separated by <1 kb, LD was dependent on the reference population. LD was lower in the specieswide sample than within populations, suggesting that low levels of LD found in inbreeding species such as B. stricta, Arabidopsis thaliana, and barley may result from broad geographic sampling that spans heterogeneous genetic groups. Finally, analyses also showed that inbreeding B. stricta and A. thaliana have approximately 45% higher recombination per kilobase than outcrossing A. lyrata.

Many or most genes in Arabidopsis transposed after the origin of the order Brassicales

Previous to this work, typical genes were thought to move from one position to another infrequently. On the contrary, we now estimate that between one-fourth and three-fourths of the genes in Arabidopsis transposed in the Brassicales. We used the CoGe comparative genomics system to perform and visualize multiple orthologous chromosomal alignments. Using this tool, we found large differences between different categories of genes. Ten of the gene families examined, including genes in most transcription factor families, exhibited a median frequency of 5% transposed genes. In contrast, other gene families were composed largely of transposed genes: NB-LRR disease-resistance genes, genes encoding MADS-box and B3 transcription factors, and genes encoding F-box proteins. A unique method involving transposition-rich regions of genome allowed us to obtain an indirect estimate of the positional stability of the average gene. The observed differences between gene families raise important questions concerning the causes and consequences of gene transposition.

Synteny conservation between two distantly-related Rosaceae genomes: Prunus (the stone fruits) and Fragaria (the strawberry)

BACKGROUND

The Rosaceae encompass a large number of economically-important diploid and polyploid fruit and ornamental species in many different genera. The basic chromosome numbers of these genera are x = 7, 8 and 9 and all have compact and relatively similar genome sizes. Comparative mapping between distantly-related genera has been performed to a limited extent in the Rosaceae including a comparison between Malus (subfamily Maloideae) and Prunus (subfamily Prunoideae); however no data has been published to date comparing Malus or Prunus to a member of the subfamily Rosoideae. In this paper we compare the genome of Fragaria, a member of the Rosoideae, to Prunus, a member of the Prunoideae.

RESULTS

The diploid genomes of Prunus (2n = 2x = 16) and Fragaria (2n = 2x = 14) were compared through the mapping of 71 anchor markers - 40 restriction fragment length polymorphisms (RFLPs), 29 indels or single nucleotide polymorphisms (SNPs) derived from expressed sequence tags (ESTs) and two simple-sequence repeats (SSRs) - on the reference maps of both genera. These markers provided good coverage of the Prunus (78%) and Fragaria (78%) genomes, with maximum gaps and average densities of 22 cM and 7.3 cM/marker in Prunus and 32 cM and 8.0 cM/marker in Fragaria.

CONCLUSION

Our results indicate a clear pattern of synteny, with most markers of each chromosome of one of these species mapping to one or two chromosomes of the other. A large number of rearrangements (36), most of which produced by inversions (27) and the rest (9) by translocations or fission/fusion events could also be inferred. We have provided the first framework for the comparison of the position of genes or DNA sequences of these two economically valuable and yet distantly-related genera of the Rosaceae.

Synteny conservation between two distantly-related Rosaceae genomes: Prunus (the stone fruits) and Fragaria (the strawberry)

BACKGROUND

The Rosaceae encompass a large number of economically-important diploid and polyploid fruit and ornamental species in many different genera. The basic chromosome numbers of these genera are x = 7, 8 and 9 and all have compact and relatively similar genome sizes. Comparative mapping between distantly-related genera has been performed to a limited extent in the Rosaceae including a comparison between Malus (subfamily Maloideae) and Prunus (subfamily Prunoideae); however no data has been published to date comparing Malus or Prunus to a member of the subfamily Rosoideae. In this paper we compare the genome of Fragaria, a member of the Rosoideae, to Prunus, a member of the Prunoideae.

RESULTS

The diploid genomes of Prunus (2n = 2x = 16) and Fragaria (2n = 2x = 14) were compared through the mapping of 71 anchor markers - 40 restriction fragment length polymorphisms (RFLPs), 29 indels or single nucleotide polymorphisms (SNPs) derived from expressed sequence tags (ESTs) and two simple-sequence repeats (SSRs) - on the reference maps of both genera. These markers provided good coverage of the Prunus (78%) and Fragaria (78%) genomes, with maximum gaps and average densities of 22 cM and 7.3 cM/marker in Prunus and 32 cM and 8.0 cM/marker in Fragaria.

CONCLUSION

Our results indicate a clear pattern of synteny, with most markers of each chromosome of one of these species mapping to one or two chromosomes of the other. A large number of rearrangements (36), most of which produced by inversions (27) and the rest (9) by translocations or fission/fusion events could also be inferred. We have provided the first framework for the comparison of the position of genes or DNA sequences of these two economically valuable and yet distantly-related genera of the Rosaceae.

High DNA sequence diversity in pericentromeric genes of the plant Arabidopsis lyrata

Differences in neutral diversity at different loci are predicted to arise due to differences in mutation rates and from the "hitchhiking" effects of natural selection. Consistent with hitchhiking models, Drosophila melanogaster chromosome regions with very low recombination have unusually low nucleotide diversity. We compared levels of diversity from five pericentromeric regions with regions of normal recombination in Arabidopsis lyrata, an outcrossing close relative of the highly selfing A. thaliana. In contrast with the accepted theoretical prediction, and the pattern in Drosophila, we found generally high diversity in pericentromeric genes, which is consistent with the observation in A. thaliana. Our data rule out balancing selection in the pericentromeric regions, suggesting that hitchhiking is more strongly reducing diversity in the chromosome arms than the pericentromere regions.

Genetic causes and consequences of the breakdown of self-incompatibility: case studies in the Brassicaceae

The genetic consequences of inbreeding is a subject that has received thorough theoretical attention and has been of interest to empirical biologists since the time of Darwin. Particularly for species with genetically controlled mechanisms to promote outcrossing (self-incompatibility or SI systems), it is expected that high levels of genetic load should accumulate through sheltering of deleterious recessive mutations. Nevertheless, transitions to selfing are common across angiosperms, which suggests that the potentially negative consequences of reduced heterozygosity and genetic diversity are balanced by other factors, such as reproductive assurance. This mini-review focuses on empirical research in the Brassicaceae to emphasize some of the genetic consequences of shifts to inbreeding in terms of mechanisms for loss of SI, changes in genetic diversity following loss of SI, and inbreeding depression in relation to outcrossing history. Despite the long history of theoretical attention, there are still some surprisingly large gaps in our understanding in each of these areas. Rather than providing a complete overview, examples are drawn predominantly from published and emerging data from Arabidopsis thaliana and its relatives to highlight recent progress and remaining questions. We are currently on the brink of major breakthroughs in understanding due both to advances in sequencing technology and a shift in focus from crop plants to natural populations, where critical factors such as population structure, phylogeography, demographic history, partial compatibility and individual variation can be taken into account when investigating the nature of the selective forces regulating mating system evolution.

Selection on amino acid substitutions in Arabidopsis

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.

Comparative mapping of Brassica juncea and Arabidopsis thaliana using Intron Polymorphism (IP) markers: homoeologous relationships, diversification and evolution of the A, B and C Brassica genomes

Background

Extensive mapping efforts are currently underway for the establishment of comparative genomics between the model plant, Arabidopsis thaliana and various Brassica species. Most of these studies have deployed RFLP markers, the use of which is a laborious and time-consuming process. We therefore tested the efficacy of PCR-based Intron Polymorphism (IP) markers to analyze genome-wide synteny between the oilseed crop, Brassica juncea (AABB genome) and A. thaliana and analyzed the arrangement of 24 (previously described) genomic block segments in the A, B and C Brassica genomes to study the evolutionary events contributing to karyotype variations in the three diploid Brassica genomes.

Results

IP markers were highly efficient and generated easily discernable polymorphisms on agarose gels. Comparative analysis of the segmental organization of the A and B genomes of B. juncea (present study) with the A and B genomes of B. napus and B. nigra respectively (described earlier), revealed a high degree of colinearity suggesting minimal macro-level changes after polyploidization. The ancestral block arrangements that remained unaltered during evolution and the karyotype rearrangements that originated in the Oleracea lineage after its divergence from Rapa lineage were identified. Genomic rearrangements leading to the gain or loss of one chromosome each between the A-B and A-C lineages were deciphered. Complete homoeology in terms of block organization was found between three linkage groups (LG) each for the A-B and A-C genomes. Based on the homoeology shared between the A, B and C genomes, a new nomenclature for the B genome LGs was assigned to establish uniformity in the international Brassica LG nomenclature code.

Conclusion

IP markers were highly effective in generating comparative relationships between Arabidopsis and various Brassica species. Comparative genomics between the three Brassica lineages established the major rearrangements, translocations and fusions pivotal to karyotype diversification between the A, B and C genomes of Brassica species. The inter-relationships established between the Brassica lineages vis-à-vis Arabidopsis would facilitate the identification and isolation of candidate genes contributing to traits of agronomic value in crop Brassicas and the development of unified tools for Brassica genomics.

Construction, alignment and analysis of twelve framework physical maps that represent the ten genome types of the genus Oryza

Bacterial artificial chromosome (BAC) fingerprint and end-sequenced physical maps representing the ten genome types of Oryza are presented

We describe the establishment and analysis of a genus-wide comparative framework composed of 12 bacterial artificial chromosome fingerprint and end-sequenced physical maps representing the 10 genome types of Oryza aligned to the O. sativa ssp. japonica reference genome sequence. Over 932 Mb of end sequence was analyzed for repeats, simple sequence repeats, miRNA and single nucleotide variations, providing the most extensive analysis of Oryza sequence to date.

Hitch-hiking to a locus under balancing selection: high sequence diversity and low population subdivision at the S-locus genomic region inArabidopsis halleri

Hitch-hiking to a site under balancing selection is expected to produce a local increase in nucleotide polymorphism and a decrease in population differentiation compared with the background genomic level, but empirical evidence supporting these predictions is scarce. We surveyed molecular diversity at four genes flanking the region controlling self-incompatibility (the S-locus) in samples from six populations of the herbaceous plantArabidopsis halleri, and compared their polymorphism with sequences from five control genes unlinked to the S-locus. As a preliminary verification, the S-locus flanking genes were shown to co-segregate withSRK, the gene involved in the self-incompatibility reaction at the pistil level. In agreement with theory, our results demonstrated a significant peak of nucleotide diversity around the S-locus as well as a significant decrease in population genetic structure in the S-locus region compared with both control genes and a set of seven unlinked microsatellite markers. This is consistent with the theoretical expectation that balancing selection is increasing the effective migration rate in subdivided populations. Although only four S-locus flanking genes were investigated, our results suggest that these two signatures of the hitch-hiking effect are localized in a very narrow genomic region.

Arabidopsis: a model genus for speciation

What genetic and epigenetic changes underlie adaptation and divergence? Arabidopsis thaliana and its relatives are increasingly being employed to address such central questions of evolutionary biology. For example, comparative, genomic and classical genetic approaches are revealing mechanisms underlying processes relevant to speciation, including mating system evolution, the effects of ploidy and other chromosomal differences, and the roles that specific genes might play in Dobzhansky-Muller type incompatibilities. The considerable body of knowledge and resources available for A. thaliana and improvements in tools and technology applied to its close relatives are opening doors for combining experimental and comparative analyses to elucidate fundamental mechanisms of evolution.

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.

Linkage maps for Arabidopsis lyrata subsp. lyrata and Arabidopsis lyrata subsp. petraea combining anonymous and Arabidopsis thaliana-derived markers

Arabidopsis lyrata, a close relative of the model plant Arabidopsis thaliana, is 1 of a few plant species for which the genome is to be entirely sequenced, which promises to yield important insights into genome evolution. Only 2 sparse linkage maps have been published, and these were based solely on markers derived from the A. thaliana genome. Because the genome of A. lyrata is practically twice as large as that of A. thaliana, the extent of map coverage of the A. lyrata genome remains uncertain. In this study, a 2-way pseudo-testcross strategy was used to construct genetic linkage maps of A. lyrata subsp. petraea and A. lyrata subsp. lyrata, using simple sequence repeat (SSR) and cleaved amplified polymorphic sequence (CAPS) markers from the A. thaliana genome, and anonymous amplified fragment length polymorphism (AFLP) markers that could potentially uncover regions unique to the A. lyrata genome. The SSR and CAPS markers largely confirmed the relationships between linkage groups in A. lyrata and A. thaliana. AFLP markers slightly increased the coverage of the A. lyrata maps, but mostly increased marker density on the linkage groups. We noted a much lower level of polymorphism and a greater segregation distortion in A. lyrata subsp. lyrata markers. The implications of these findings for the sequencing of the A. lyrata genome are discussed.

Comparative physical mapping between Oryza sativa (AA genome type) and O. punctata (BB genome type)

A comparative physical map of the AA genome (Oryza sativa) and the BB genome (O. punctata) was constructed by aligning a physical map of O. punctata, deduced from 63,942 BAC end sequences (BESs) and 34,224 fingerprints, onto the O. sativa genome sequence. The level of conservation of each chromosome between the two species was determined by calculating a ratio of BES alignments. The alignment result suggests more divergence of intergenic and repeat regions in comparison to gene-rich regions. Further, this characteristic enabled localization of heterochromatic and euchromatic regions for each chromosome of both species. The alignment identified 16 locations containing expansions, contractions, inversions, and transpositions. By aligning 40% of the punctata BES on the map, 87% of the punctata FPC map covered 98% of the O. sativa genome sequence. The genome size of O. punctata was estimated to be 8% larger than that of O. sativa with individual chromosome differences of 1.5-16.5%. The sum of expansions and contractions observed in regions >500 kb were similar, suggesting that most of the contractions/expansions contributing to the genome size difference between the two species are small, thus preserving the macro-collinearity between these species, which diverged approximately 2 million years ago.

Comparative genetic mapping in Boechera stricta, a close relative of Arabidopsis

The angiosperm family Brassicaceae contains both the research model Arabidopsis (Arabidopsis thaliana) and the agricultural genus Brassica. Comparative genomics in the Brassicaceae has largely focused on direct comparisons between Arabidopsis and the species of interest. However, the reduced genome size and chromosome number (n = 5) of Arabidopsis complicates comparisons. Arabidopsis shows extensive genome and chromosome reshuffling compared to its close relatives Arabidopsis lyrata and Capsella rubella, both with n = 8. To facilitate comparative genomics across the Brassicaceae we recently outlined a system of 24 conserved chromosomal blocks based on their positions in an ancestral karyotype of n = 8, rather than by their position in Arabidopsis. In this report we use this system as a tool to understand genome structure and evolution in Boechera stricta (n = 7). B. stricta is a diploid, sexual, and highly self-fertilizing species occurring in mostly montane regions of western North America. We have created an F(2) genetic map of B. stricta based on 192 individuals scored at 196 microsatellite and candidate gene loci. Single-nucleotide polymorphism genotyping of 94 of the loci was done simultaneously using an Illumina bead array. The total map length is 725.8 cM, with an average marker spacing of 3.9 cM. There are no gaps greater than 19.3 cM. The chromosomal reduction from n = 8 to n = 7 and other genomic changes in B. stricta likely involved a pericentric inversion, a chromosomal fusion, and two reciprocal translocations that are easily visualized using the genomic blocks. Our genetic map will facilitate the analysis of ecologically relevant quantitative variation in Boechera. Sequence data from this article can be found in the GenBank/EMBL data libraries under accession numbers DU 667459 to DU 708532.

The genetic basis of zinc tolerance in the metallophyte Arabidopsis halleri ssp. halleri (Brassicaceae): an analysis of quantitative trait loci

The species Arabidopsis halleri, an emerging model for the study of heavy metal tolerance and accumulation in plants, has evolved a high level of constitutive zinc tolerance. Mapping of quantitative trait loci (QTL) was used to investigate the genetic architecture of zinc tolerance in this species. A first-generation backcross progeny of A. halleri ssp. halleri from a highly contaminated industrial site and its nontolerant relative A. lyrata ssp. petraea was produced and used for QTL mapping of zinc tolerance. A genetic map covering most of the A. halleri genome was constructed using 85 markers. Among these markers, 65 were anchored in A. thaliana and revealed high synteny with other Arabidopsis genomes. Three QTL of comparable magnitude on three different linkage groups were identified. At all QTL positions zinc tolerance was enhanced by A. halleri alleles, indicating directional selection for higher zinc tolerance in this species. The two-LOD support intervals associated with these QTL cover 24, 4, and 13 cM. The importance of each of these three regions is emphasized by their colocalization with HMA4, MTP1-A, and MTP1-B, respectively, three genes well known to be involved in metal homeostasis and tolerance in plants.

Transcription-related mutations and GC content drive variation in nucleotide substitution rates across the genomes of Arabidopsis thaliana and Arabidopsis lyrata

Background

There has been remarkably little study of nucleotide substitution rate variation among plant nuclear genes, in part because orthology is difficult to establish. Orthology is even more problematic for intergenic regions of plant nuclear genomes, because plant genomes generally harbor a wealth of repetitive DNA. In theory orthologous intergenic data is valuable for studying rate variation because nucleotide substitutions in these regions should be under little selective constraint compared to coding regions. As a result, evolutionary rates in intergenic regions may more accurately reflect genomic features, like recombination and GC content, that contribute to nucleotide substitution.

Results

We generated a set of 66 intergenic sequences in Arabidopsis lyrata, a close relative of Arabidopsis thaliana. The intergenic regions included transposable element (TE) remnants and regions flanking the TEs. We verified orthology of these amplified regions both by comparison of existing A. lyrata – A. thaliana genetic maps and by using molecular features. We compared substitution rates among the 66 intergenic loci, which exhibit ~5-fold rate variation, and compared intergenic rates to a set of 64 orthologous coding sequences. Our chief observations were that the average rate of nucleotide substitution is slower in intergenic regions than in synonymous sites, that rate variation in both intergenic and coding regions correlate with GC content, that GC content alone is not sufficient to explain differences in rates between intergenic and coding regions, and that rates of evolution in intergenic regions correlate negatively with gene density.

Conclusion

Our observations indicated that mutation rates vary among genomics regions as a function of base composition, suggesting that previous observations of "selective constraint" on non-coding regions could more accurately be attributed to a GC effect instead of selection. The negative correlation between nucleotide substitution rate and gene density provides a potential neutral explanation for a previously documented correlation between gene density and polymorphism levels within A. thaliana. Finally, we discuss potential forces that could contribute to rapid synonymous rates, and provide evidence to suggest that transcription-related mutation contributes to rate differences between intergenic and synonymous sites.

Evolution and genetic differentiation among relatives of Arabidopsis thaliana

Arabidopsis thaliana is one of the most intensively studied plant species. More recently, information is accumulating about its closest relatives, the former genus Cardaminopsis. A. thaliana diverged from these relatives, actually treated within three major lineages (Arabidopsis lyrata, Arabidopsis halleri, and Arabidopsis arenosa), approximately 5 mya. Significant karyotype evolution in A. thaliana with base chromosome number reduction from x=8 to x=5 might indicate and favor effective genetic isolation from these other species, although hybrids are occurring naturally and have been also constituted under controlled conditions. We tested the evolutionary significance to separate the x=5 from the x=8 lineage using DNA sequence data from the plastome and the nuclear ribosomal DNA based on an extensive, representative worldwide sampling of nearly all taxonomic entities. We conclude that (i) A. thaliana is clearly separated phylogenetically from the x=8 lineage, (ii) five major lineages outside A. thaliana can be identified (A. lyrata, A. arenosa, A. halleri, Arabidopsis croatica, and Arabidopsis pedemontana) together with Arabidopsis cebennensis, and (iii) centers of genetic and morphological diversity are mostly in congruence and are located close to the Balkans in Austria and Slovakia outside glaciated and permafrost regions with few notable exceptions.

Dependence of paracentric inversion rate on tract length

Background

We develop a Bayesian method based on MCMC for estimating the relative rates of pericentric and paracentric inversions from marker data from two species. The method also allows estimation of the distribution of inversion tract lengths.

Results

We apply the method to data from Drosophila melanogaster and D. yakuba. We find that pericentric inversions occur at a much lower rate compared to paracentric inversions. The average paracentric inversion tract length is approx. 4.8 Mb with small inversions being more frequent than large inversions.

If the two breakpoints defining a paracentric inversion tract are uniformly and independently distributed over chromosome arms there will be more short tract-length inversions than long; we find an even greater preponderance of short tract lengths than this would predict. Thus there appears to be a correlation between the positions of breakpoints which favors shorter tract lengths.

Conclusion

The method developed in this paper provides the first statistical estimator for estimating the distribution of inversion tract lengths from marker data. Application of this method for a number of data sets may help elucidate the relationship between the length of an inversion and the chance that it will get accepted.

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.

A quantitative trait loci analysis of zinc hyperaccumulation in Arabidopsis halleri

The mechanisms of metal hyperaccumulation are still not understood, so we conducted a quantitative trait locus (QTL) analysis of zinc (Zn) hyperaccumulation in Arabidopsis halleri, in a cross between this and its sister species, A. petraea, in order to determine the number and approximate location of the genomic regions significantly contributing to this adaptation. An F2 cross between the two species was made, and the leaf Zn concentration of 92 individuals was measured at both low (10 microm) and high (100 microm) Zn concentrations. Twenty-five markers were established that were distributed on all of the eight chromosomes. Mapping of the markers established that they were essentially collinear with previous studies. QTLs exceeding a logarithm to the base 10 of the odds (LOD) value of 3 were found on chromosomes 4 (low Zn), 6 (high Zn) and 7 (both high and low Zn). Evidence for a QTL on chromosome 3 (low Zn) was also found. This analysis validates a previously used method of QTL analysis, based on microarray analysis of segregating families. Genes that have altered during the evolution of this character should also be QTL: this analysis calls into question a number of candidate genes from consideration as such primary genes because they do not appear to be associated with QTLs.

Chromosome arrangement and nuclear architecture but not centromeric sequences are conserved between Arabidopsis thaliana and Arabidopsis lyrata

In contrast to the situation described for mammals and Drosophila, chromosome territory (CT) arrangement and somatic homologous pairing in interphase nuclei of Arabidopsis thaliana (n = 5) are predominantly random except for a more frequent association of the chromosomes bearing a homologous nucleolus organizer region. To find out whether this chromosome arrangement is also characteristic for other species of the genus Arabidopsis, we investigated Arabidopsis lyrata ssp. lyrata (n = 8), one of the closest relatives of A. thaliana. First, we determined the size of each chromosome and chromosome arm, the sequence type of centromeric repeats and their distribution between individual centromeres and the position of the 5S/45S rDNA arrays in A. lyrata. Then we demonstrated that CT arrangement, homologous pairing and sister chromatid alignment of distinct euchromatic and/or heterochromatic regions within A. lyrata interphase nuclei are similar to that in A. thaliana nuclei. Thus, the arrangement of interphase chromosomes appears to be conserved between both taxa that diverged about 5 million years ago. Since the chromosomes of A. lyrata resemble those of the presumed ancestral karyotype, a similar arrangement of interphase chromosomes is also to be expected for other closely related diploid species of the Brassicaceae family.

Triploid origin of the gibel carp as revealed by 5S rDNA localization and chromosome painting

5S ribosomal DNA (rDNA) was isolated and sequenced from the gibel carp Carassius auratus gibelio with 162 chromosomes and crucian carp Carassius auratus with 100 chromosomes, and fluorescent probes for chromosome localization were prepared to ascertain the ploidy origin and evolutionary relationship between the two species. Using fluorescence in-situ hybridization (FISH), major 5S rDNA signals were localized to the short arms of three subtelocentric chromosomes in the gibel carp and to the short arms of two subtelocentrics in the crucian carp. In addition, some minor signals were detected on other chromosomes of both species. Simultaneously, six chromosomes were microdissected from the gibel carp metaphase spreads using glass needles, and the isolated chromosomes were amplified in vitro by degenerate oligonucleotide primed-polymerase chain reaction (DOP-PCR). Significantly, when the DOP-PCR-generated probes prepared from each single chromosome were hybridized, three same-sized chromosomes were painted in each gibel carp metaphase, whereas only two painted chromosomes were observed in each crucian carp metaphase spread. The data indicate that gibel carp is of triploid origin in comparison with diploid crucian carp.

The ABC's of comparative genomics in the Brassicaceae: building blocks of crucifer genomes

In this review we summarize recent advances in our understanding of phylogenetics, polyploidization and comparative genomics in the family Brassicaceae. These findings pave the way for a unified comparative genomic framework. We integrate several of these findings into a simple system of 24 conserved chromosomal blocks (labeled A-X). The naming, order, orientation and color-coding of these blocks are based on their positions in a proposed ancestral karyotype (n=8), rather than by their position in the reduced genome of Arabidopsis thaliana (n=5). We show how these crucifer building blocks can be rearranged to model the genome structures of A. thaliana, Arabidopsis lyrata, Capsella rubella and Brassica rapa. A framework for comparison between species is timely because several crucifer genome-sequencing projects are underway.

Comparative genome mapping among Picea glauca, P. mariana x P. rubens and P. abies, and correspondence with other Pinaceae

A composite linkage map was constructed from four individual maps for the conifer Picea glauca (Moench) Voss, from anonymous and gene-specific markfers (714 AFLPs, 38 SSRs, and 53 ESTPs). A total of 12 linkage groups were delineated with an average marker density of 2.7 cM. Macro-synteny and macro-colinearity comparisons with two other composite linkage maps developed for the species complex P. mariana (Mill.) B.S.P. x P. rubens Sarg., and for P. abies (L.) Karst. revealed an identical number of linkage groups and a remarkable conservation of the gene content and gene order of linkage groups over the million years since the split between these taxa. Identical gene order among taxa was observed for 10 of the 12 assembled composite linkage groups. The discovery of one breakdown in synteny between P. glauca and the other two taxa indicated the occurrence of an inter-chromosomal rearrangement involving an insertional translocation. Analysis of marker colinearity also revealed a putative segmental duplication. The combined information from these three Picea genomes validated and improved large-scale genome comparisons at the inter-generic level in the family Pinaceae by allowing for the identification of 11 homoeologous linkage groups between Picea and Pinus, and nine such groups between Picea and Pseudotsuga menziesii. Notably, the analysis of synteny among the three genera revealed a putative case of chromosomal fission and an inter-chromosomal rearrangement in the genome of P. menziesii. Both of these changes are inter-connected, indicating much instability in this part of the P. menziesii genome. Overall, the macro-structure of the Pinaceae genome was well conserved, which is notable given the Cretaceous origin of its main lineages.

Testing chromosomal phylogenies and inversion breakpoint reuse in Drosophila

A combination of cytogenetic and bioinformatic procedures was used to test the chromosomal phylogeny relating Drosophila buzzatii with D. repleta. Chromosomes X and 2, harboring most of the inversions fixed between these two species, were analyzed. First, chromosomal segments conserved during the divergence of the two species were identified by comparative in situ hybridization to the D. repleta chromosomes of 180 BAC clones from a BAC-based physical map of the D. buzzatii genome. These conserved segments were precisely delimited with the aid of clones containing inversion breakpoints. Then GRIMM software was used to estimate the minimum number of rearrangements necessary to transform one genome into the other and identify all possible rearrangement scenarios. Finally, the most plausible inversion trajectory was tested by hybridizing 12 breakpoint-bearing BAC clones to the chromosomes of seven other species in the repleta group. The results show that chromosomes X and 2 of D. buzzatii and D. repleta differ by 12 paracentric inversions. Nine of them are fixed in chromosome 2 and entail two breakpoint reuses. Our results also show that the cytological relationship between D. repleta and D. mercatorum is closer than that between D. repleta and D. peninsularis, and we propose that the phylogenetic relationships in this lineage of the repleta group be reconsidered. We also estimated the rate of rearrangement between D. repleta and D. buzzatii and conclude that rates within the genus Drosophila vary substantially between lineages, even within a single species group.

Duplication of centromeric histone H3 (HTR12) gene in Arabidopsis halleri and A. lyrata, plant species with multiple centromeric satellite sequences

Arabidopsis halleri and lyrata have three different major centromeric satellite sequences, a unique finding for a diploid Arabidopsis species. Since centromeric histones coevolve with centromeric satellites, these proteins would be predicted to show signs of selection when new centromere satellites have recently arisen. We isolated centromeric protein genes from A. halleri and lyrata and found that one of them, HTR12 (CENP-A), is duplicated, while CENP-C is not. Phylogenetic analysis indicates that the HTR12 duplication occurred after these species diverged from A. thaliana. Genetic mapping shows that HTR12 copy B has the same genomic location as the A. thaliana gene; the other copy (A, at the other end of the same chromosome) is probably the new copy. To test for selection since the duplication, we surveyed diversity at both HTR12 loci within A. lyrata. Overall, there is no strong evidence for an "evolutionary arms race" causing multiple replacement substitutions. The A. lyrata HTR12B sequences fall into three classes of haplotypes, apparently maintained for a long time, but they all encode the same amino acid sequence. In contrast, HTR12A has low diversity, but many variants are amino acid replacements, possibly due to independent selective sweeps within populations of the species.

Zinc-dependent global transcriptional control, transcriptional deregulation, and higher gene copy number for genes in metal homeostasis of the hyperaccumulator Arabidopsis halleri

The metal hyperaccumulator Arabidopsis halleri exhibits naturally selected zinc (Zn) and cadmium (Cd) hypertolerance and accumulates extraordinarily high Zn concentrations in its leaves. With these extreme physiological traits, A. halleri phylogenetically belongs to the sister clade of Arabidopsis thaliana. Using a combination of genome-wide cross species microarray analysis and real-time reverse transcription-PCR, a set of candidate genes is identified for Zn hyperaccumulation, Zn and Cd hypertolerance, and the adjustment of micronutrient homeostasis in A. halleri. Eighteen putative metal homeostasis genes are newly identified to be more highly expressed in A. halleri than in A. thaliana, and 11 previously identified candidate genes are confirmed. The encoded proteins include HMA4, known to contribute to root-shoot transport of Zn in A. thaliana. Expression of either AtHMA4 or AhHMA4 confers cellular Zn and Cd tolerance to yeast (Saccharomyces cerevisiae). Among further newly implicated proteins are IRT3 and ZIP10, which have been proposed to contribute to cytoplasmic Zn influx, and FRD3 required for iron partitioning in A. thaliana. In A. halleri, the presence of more than a single genomic copy is a hallmark of several highly expressed candidate genes with possible roles in metal hyperaccumulation and metal hypertolerance. Both A. halleri and A. thaliana exert tight regulatory control over Zn homeostasis at the transcript level. Zn hyperaccumulation in A. halleri involves enhanced partitioning of Zn from roots into shoots. The transcriptional regulation of marker genes suggests that in the steady state, A. halleri roots, but not the shoots, act as physiologically Zn deficient under conditions of moderate Zn supply.

Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.

Genome evolution in Arabidopsis/Brassica: conservation and divergence of ancient rearranged segments and their breakpoints

Since the tetraploidization of the Arabidopsis thaliana ancestor 30-35 million years ago (Mya), a wave of chromosomal rearrangements have modified its genome architecture. The dynamics of this process is unknown, as it has so far been impossible to date individual rearrangement events. In this paper, we present evidence demonstrating that the majority of rearrangements occurred before the Arabidopsis-Brassica split 20-24 Mya, and that the segmental architecture of the A. thaliana genome is predominantly conserved in Brassica. This finding is based on the conservation of four rearrangement breakpoints analysed by fluorescence in situ hybridization (FISH) and RFLP mapping of three A. thaliana chromosomal regions. For this purpose, 95 Arabidopsis bacterial artificial chromosomes (BACs) spanning a total of 8.25 Mb and 81 genetic loci for 36 marker genes were studied in the Brassica oleracea genome. All the regions under study were triplicated in the B. oleracea genome, confirming the hypothesis of Brassica ancestral genome triplication. However, whilst one of the breakpoints was conserved at one locus, it was not at the two others. Further comparison of their organization may indicate that the evolution of the hexaploid Brassica progenitor proceeded by several events, separated in time. Genetic mapping and reprobing with rDNA allowed assignment of the regions to particular Brassica chromosomes. Based on this study of regional organization and evolution, a new insight into polyploidization/diploidization cycles is proposed.

History, protohistory and prehistory of the Arabidopsis thaliana chromosome complement

We propose an evolutionary scenario that could have shaped the modern Arabidopsis thaliana genome, which began with the reduction in chromosome number from n=8 to n=5 in the past 4 million to 5 million years as a result of chromosome fusion. The scenario also includes three ancient polyploidizations: the most recent occurred in an early Brassicaceae with n=4 chromosomes 24 million to 40 million years ago. The two other polyploidizations occurred after the emergence of the Eudicots and the Angiosperms, respectively. Angiosperm evolution includes recurrent cycles of genome duplication and gene and chromosome reorganizations.

Balancing selection and low recombination affect diversity near the self-incompatibility loci of the plant Arabidopsis lyrata

The self-incompatibility (S-) locus region of plants in the Brassica family is a small genome region. In Arabidopsis lyrata, the S-genes, SRK and SCR, encode the functional female and pollen recognition proteins, which must be coadapted to maintain correct associations between the two component genes, and thus self-incompatibility (SI). Recombinants would be self-compatible and thus probably disadvantageous in self-incompatible species. Therefore, tight linkage between the two genes in incompatibility systems is predicted to evolve to avoid producing such recombinant haplotypes. The evolution of low recombination in S-locus regions has not been rigorously tested. To test whether these regions' per-nucleotide recombination rates differ from those elsewhere in the genome, and to investigate whether the A. lyrata S-loci have the predicted effect on diversity in their immediate genome region, we studied diversity in genes that are linked to the S-loci but are not involved in incompatibility and are not under balancing selection. Compared with other A. lyrata loci, genes linked to the S-loci have extraordinarily high polymorphism. Our estimated recombination in this region, from fitting a model of the effects of S-allele polymorphism on linked neutral sites, supports the hypothesis of locally suppressed recombination around the S-locus.

Construction of a watermelon BAC library and identification of SSRs anchored to melon or Arabidopsis genomes

A bacterial artificial chromosome (BAC) library was constructed for watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus) with an average insert-size of 106 kb, providing 21 haploid genome equivalents. The library was used to identify BAC clones that are anchored to probes evenly distributed on the genomes of melon or Arabidopsis. Twenty eight probes (representing 66% of the tested probes) from melon and 30 probes (65%) from Arabidopsis identified positive BAC clones. Two methods were implemented to identify SSRs from the positively hybridizing BAC clones. First, analysis of BAC end sequences revealed 37 SSRs. For the second method, pooled DNA of BACs identified by the melon probes was used to develop a shotgun library. The library was then screened with synthetic SSR oligonucleotides by hybridization. Sequence analysis of positively hybridizing shotgun clones revealed 142 different SSRs. Thirty eight SSRs were characterized using three watermelon cultivars, five plant introduction (PI) accessions of C. lanatus var lanatus and four PIs of C. lanatus var citroides. Of these, 36 (95%) were found to be polymorphic with up to six alleles per marker. Polymorphism information content values for polymorphic markers varied between 0.22 and 0.79 with an average of 0.53. The methods described herein will be valuable for the construction of a watermelon linkage map with SSRs evenly distributed on its genome that is anchored to the genomes of melon and Arabidopsis.

Centromere locations and associated chromosome rearrangements in Arabidopsis lyrata and A. thaliana

We analyzed linkage and chromosomal positions of genes in A. lyrata ssp. petraea that are located near the centromere (CEN) regions of A. thaliana, using at least two genes from the short and long arms of each chromosome. In our map, genes from all 10 A. thaliana chromosome arms are also tightly linked in A. lyrata. Genes from the regions on the two sides of CEN5 have distant map localizations in A. lyrata (genes on the A. thaliana short-arm genes are on linkage group AL6, and long-arm genes are on AL7), but genes from the other four A. thaliana centromere regions remain closely linked in A. lyrata. The observation of complete linkage between short- and long-arm centromere genes, but not between genes in other genome regions that are separated by similar physical distances, suggests that crossing-over frequencies near the A. lyrata ssp. petraea centromere regions are low, as in A. thaliana. Thus, the centromere positions appear to be conserved between A. thaliana and A. lyrata, even though three centromeres have been lost in A. thaliana, and the core satellite sequences in the two species are very different. We can now definitively identify the three centromeres that were eliminated in the fusions that formed the A. thaliana chromosomes. However, we cannot tell whether genes were lost along with these centromeres, because such genes are absent from the A. thaliana genome, which is the sole source of markers for our mapping.

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.

Stability of large segmental duplications in the yeast genome

The high level of gene redundancy that characterizes eukaryotic genomes results in part from segmental duplications. Spontaneous duplications of large chromosomal segments have been experimentally demonstrated in yeast. However, the dynamics of inheritance of such structures and their eventual fixation in populations remain largely unsolved. We analyzed the stability of a vast panel of large segmental duplications in Saccharomyces cerevisiae (from 41 kb for the smallest to 268 kb for the largest). We monitored the stability of three different types of interchromosomal duplications as well as that of three intrachromosomal direct tandem duplications. In the absence of any selective advantage associated with the presence of the duplication, we show that a duplicated segment internally translocated within a natural chromosome is stably inherited both mitotically and meiotically. By contrast, large duplications carried by a supernumerary chromosome are highly unstable. Duplications translocated into subtelomeric regions are lost at variable rates depending on the location of the insertion sites. Direct tandem duplications are lost by unequal crossing over, both mitotically and meiotically, at a frequency proportional to their sizes. These results show that most of the duplicated structures present an intrinsic level of instability. However, translocation within another chromosome significantly stabilizes a duplicated segment, increasing its chance to get fixed in a population even in the absence of any immediate selective advantage conferred by the duplicated genes.

Mechanisms of chromosome number reduction in Arabidopsis thaliana and related Brassicaceae species

Evolution of chromosome complements can be resolved by genome sequencing, comparative genetic mapping, and comparative chromosome painting. Previously, comparison of genetic maps and gene-based phylogenies suggested that the karyotypes of Arabidopsis thaliana (n = 5) and of related species with six or seven chromosome pairs were derived from an ancestral karyotype with eight chromosome pairs. To test this hypothesis, we applied multicolor chromosome painting using contiguous bacterial artificial chromosome pools of A. thaliana arranged according to the genetic maps of Arabidopsis lyrata and Capsella rubella (both n = 8) to A. thaliana, A. lyrata, Neslia paniculata, Turritis glabra, and Hornungia alpina. This approach allowed us to map the A. lyrata centromeres as a prerequisite to defining a putative ancestral karyotype (n = 8) and to elucidate the evolutionary mechanisms that shaped the karyotype of A. thaliana and its relatives. We conclude that chromosome "fusions" in A. thaliana resulted from (i) generation of acrocentric chromosomes by pericentric inversions, (ii) reciprocal translocation between two chromosomes (one or both acrocentric), and (iii) elimination of a minichromosome that arose in addition to the "fusion chromosome." Comparative chromosome painting applied to N. paniculata (n = 7), T. glabra (n = 6), and H. alpina (n = 6), for which genetic maps are not available, revealed chromosomal colinearity between all species tested and allowed us to reconstruct the evolution of their chromosomes from a putative ancestral karyotype (n = 8). Although involving different ancestral chromosomes, chromosome number reduction followed similar routes as found within the genus Arabidopsis.