Recombination: an underappreciated factor in the evolution of plant genomes

Chromosomal phylogeny and karyotype evolution in x=7 crucifer species (Brassicaceae)

Karyotype evolution in species with identical chromosome number but belonging to distinct phylogenetic clades is a long-standing question of plant biology, intractable by conventional cytogenetic techniques. Here, we apply comparative chromosome painting (CCP) to reconstruct karyotype evolution in eight species with x=7 (2n=14, 28) chromosomes from six Brassicaceae tribes. CCP data allowed us to reconstruct an ancestral Proto-Calepineae Karyotype (PCK; n=7) shared by all x=7 species analyzed. The PCK has been preserved in the tribes Calepineae, Conringieae, and Noccaeeae, whereas karyotypes of Eutremeae, Isatideae, and Sisymbrieae are characterized by an additional translocation. The inferred chromosomal phylogeny provided compelling evidence for a monophyletic origin of the x=7 tribes. Moreover, chromosomal data along with previously published gene phylogenies strongly suggest the PCK to represent an ancestral karyotype of the tribe Brassiceae prior to its tribe-specific whole-genome triplication. As the PCK shares five chromosomes and conserved associations of genomic blocks with the putative Ancestral Crucifer Karyotype (n=8) of crucifer Lineage I, we propose that both karyotypes descended from a common ancestor. A tentative origin of the PCK via chromosome number reduction from n=8 to n=7 is outlined. Comparative chromosome maps of two important model species, Noccaea caerulescens and Thellungiella halophila, and complete karyotypes of two purported autotetraploid Calepineae species (2n=4x=28) were reconstructed by CCP.

A single Banana streak virus integration event in the banana genome as the origin of infectious endogenous pararetrovirus

Sequencing of plant nuclear genomes reveals the widespread presence of integrated viral sequences known as endogenous pararetroviruses (EPRVs). Banana is one of the three plant species known to harbor infectious EPRVs. Musa balbisiana carries integrated copies of Banana streak virus (BSV), which are infectious by releasing virions in interspecific hybrids. Here, we analyze the organization of the EPRV of BSV Goldfinger (BSGfV) present in the wild diploid M. balbisiana cv. Pisang Klutuk Wulung (PKW) revealed by the study of Musa bacterial artificial chromosome resources and interspecific genetic cross. cv. PKW contains two similar EPRVs of BSGfV. Genotyping of these integrants and studies of their segregation pattern show an allelic insertion. Despite the fact that integrated BSGfV has undergone extensive rearrangement, both EPRVs contain the full-length viral genome. The high degree of sequence conservation between the integrated and episomal form of the virus indicates a recent integration event; however, only one allele is infectious. Analysis of BSGfV EPRV segregation among an F1 population from an interspecific genetic cross revealed that these EPRV sequences correspond to two alleles originating from a single integration event. We describe here for the first time the full genomic and genetic organization of the two EPRVs of BSGfV present in cv. PKW in response to the challenge facing both scientists and breeders to identify and generate genetic resources free from BSV. We discuss the consequences of this unique host-pathogen interaction in terms of genetic and genomic plant defenses versus strategies of infectious BSGfV EPRVs.

Sequence level analysis of recently duplicated regions in soybean [Glycine max (L.) Merr.] genome

A single recessive gene, rxp, on linkage group (LG) D2 controls bacterial leaf-pustule resistance in soybean. We identified two homoeologous contigs (GmA and GmA') composed of five bacterial artificial chromosomes (BACs) during the selection of BAC clones around Rxp region. With the recombinant inbred line population from the cross of Pureunkong and Jinpumkong 2, single-nucleotide polymorphism and simple sequence repeat marker genotyping were able to locate GmA' on LG A1. On the basis of information in the Soybean Breeders Toolbox and our results, parts of LG A1 and LG D2 share duplicated regions. Alignment and annotation revealed that many homoeologous regions contained kinases and proteins related to signal transduction pathway. Interestingly, inserted sequences from GmA and GmA' had homology with transposase and integrase. Estimation of evolutionary events revealed that speciation of soybean from Medicago and the recent divergence of two soybean homoeologous regions occurred at 60 and 12 million years ago, respectively. Distribution of synonymous substitution patterns, K(s), yielded a first secondary peak (mode K(s) = 0.10-0.15) followed by two smaller bulges were displayed between soybean homologous regions. Thus, diploidized paleopolyploidy of soybean genome was again supported by our study.

Characterization of chromosome ends on the basis of the structure of TrsA subtelomeric repeats in rice (Oryza sativa L.)

Subtelomeres contain species-specific repetitive sequences. We characterized rice chromosome ends on the basis of the structure of TrsA, a subtelomeric repetitive sequence of rice. Among the 24 chromosome arms, TrsA was arrayed in tandem on the ends of five: 5L, 6S, 8L, 9L, and 12L. TrsA sequences were arranged in discrete clusters of 3-106 copies in a chromosome-specific manner, instead of being distributed uniformly throughout the subtelomeric regions. The clusters were located at the distal-most end of the sequenced region in 5L, 6S, 8L, and 9L, but in 12L expressed genes were present distal to the clusters. Thus, rice subtelomeres are composed of discrete clusters of a TrsA-rich region and a gene-rich region with high transcriptional activity. Intra-chromosomal duplications have resulted in a striking degree of variation in the number and distribution of TrsAs, suggesting that the areas near the ends of the chromosomes are dynamic and variable.

No effect of recombination on the efficacy of natural selection in primates

Population genetic theory suggests that natural selection should be less effective in regions of low recombination, potentially leading to differences in rates of adaptation among recombination environments. To date, this prediction has mainly been tested in Drosophila, with somewhat conflicting results. We investigated the association between human recombination rates and adaptation in primates, by considering rates of protein evolution (measured by d(N)/d(S)) between human, chimpanzee, and rhesus macaque. We found no correlation between either broad- or fine-scale rates of recombination and rates of protein evolution, once GC content is taken into account. Moreover, genes in regions of very low recombination, which are expected to show the most pronounced reduction in the efficacy of selection, do not evolve at a different rate than other genes. Thus, there is no evidence for differences in the efficacy of selection across recombinational environments. An interesting implication is that indirect selection for recombination modifiers has probably been a weak force in primate evolution.

Evolutionary origin of the segmental duplication encompassing the wheat GLU-B1 locus encoding the overexpressed Bx7 (Bx7OE) high molecular weight glutenin subunit

Sequencing of a BAC clone encompassing the Glu-B1 locus in Glenlea, revealed a 10.3 Kb segmental duplication including the Bx7 gene and flanking an LTR retroelement. To better understand the evolution of this locus, two collections of wheat were surveyed. The first consisted of 96 diploid and tetraploid species accessions while the second consisted of 316 Triticum aestivum cultivars and landraces from 41 countries. The genotypes were first characterized by SDS-PAGE and a total of 40 of the 316 T. aestivum accessions were found to display the overexpressed Bx7 phenotype (Bx7OE). Three lines from the 96 diploid/tetraploid collection also displayed the stronger intensity staining characteristic of the Bx7(OE) subunit. The relative amounts of the Bx7 subunit to total HMW-GS were quantified by RP-HPLC for all Bx7OE accessions and a number of checks. The entire collection was assessed for the presence of four DNA markers namely an 18 bp indel of the coding region of Bx7 variant alleles, a 43 bp indel of the 5'-region and the left and right junctions of the LTR retrotransposon borders and the duplicated segment. All 43 accessions found to have the Bx7OE subunit by SDS-PAGE and RP-HPLC produced the four diagnostic PCR amplicons. None of the lines without the Bx7OE had the LTR retroelement/duplication genomic structure. However, the 18 and 43 bp indel were found in accessions other than Bx7OE. These results indicate that the overexpression of the Bx7 HMW-GS is likely the result of a single event, i.e., a gene duplication at the Glu-B1 locus mediated by the insertion of a retroelement. Also, the 18 and 43 bp indels pre-date the duplication event. Allelic variants Bx7*, Bx7 with and without 43 bp insert and Bx7OE were found in both tetraploid and hexaploid collections and shared the same genomic organization. Though the possibility of introgression from T. aestivum to T. turgidum cannot be ruled out, the three structural genomic changes of the B-genome taken together support the hypothesis of multiple polyploidization events involving different tetraploid progenitors.

MareyMap: an R-based tool with graphical interface for estimating recombination rates

Comparing genetic and physical maps (the so-called Marey map approach) is still the most widely used approach to estimate genome-wide recombination rates. Remarkably, there is no available bioinformatics tool specifically devoted to Marey map approach. Here, we developed such a tool called MareyMap based on GNU R and Tcl/Tk. MareyMap offers a user-friendly graphical interface and includes useful features, such as data cleaning process, sophisticated interpolation methods to estimate local rates, possibility of complex queries, various range of import-export files. Moreover, MareyMap comes with ready-to-use maps for human, Drosophila, Caenorhabditis elegans and Arabidopsis. MareyMap has been made so that it can be easily upgraded with new data and interpolation methods.

AVAILABILITY

http://pbil.univ-lyon1.fr/software/mareymap/.

Sex-specific crossover distributions and variations in interference level along Arabidopsis thaliana chromosome 4

In many species, sex-related differences in crossover (CO) rates have been described at chromosomal and regional levels. In this study, we determined the CO distribution along the entire Arabidopsis thaliana Chromosome 4 (18 Mb) in male and female meiosis, using high density genetic maps built on large backcross populations (44 markers, >1,300 plants). We observed dramatic differences between male and female map lengths that were calculated as 88 cM and 52 cM, respectively. This difference is remarkably parallel to that between the total synaptonemal complex lengths measured in male and female meiocytes by immunolabeling of ZYP1 (a component of the synaptonemal complex). Moreover, CO landscapes were clearly different: in particular, at both ends of the map, male CO rates were higher (up to 4-fold the mean value), whereas female CO rates were equal or even below the chromosomal average. This unique material gave us the opportunity to perform a detailed analysis of CO interference on Chromosome 4 in male and female meiosis. The number of COs per chromosome and the distances between them clearly departs from randomness. Strikingly, the interference level (measured by coincidence) varied significantly along the chromosome in male meiosis and was correlated to the physical distance between COs. The significance of this finding on the relevance of current CO interference models is discussed.

Meiotic crossovers between homologous chromosomes ensure their proper segregation to generate ultimately gametes. They also create new allelic combinations which contribute to the diversity of traits among individuals. In all eukaryotes, the number and the localization of crossovers along chromosomes are not random. In addition, crossovers are not independent of each other: the occurrence of a crossover lowers the probability that another crossover arises in its vicinity. The mechanism of this phenomenon, called “crossover interference,” is one of the most challenging puzzles that geneticists have been faced with in the last century. In this paper, we precisely described the distribution of crossovers along Chromosome 4 of the model plant species Arabidopsis thaliana, separately in male and female meiosis. Interestingly, we observed that crossovers are 1.7 more numerous in male than in female meiosis, and this increase is especially marked at the ends of the chromosome. Moreover, our results provide the first evidence that the level of interference along a chromosome is not a constant and is correlated with the physical distance between crossovers. These results shed new light on the determinism of crossover localization and could have important outcomes on the relevance of current models of crossover interference.

Plant domestication, a unique opportunity to identify the genetic basis of adaptation

Despite the fundamental role of plant domestication in human history and the critical importance of a relatively small number of crop plants to modern societies, we still know little about adaptation under domestication. Here we focus on efforts to identify the genes responsible for adaptation to domestication. We start from a historical perspective, arguing that Darwin's conceptualization of domestication and unconscious selection provides valuable insight into the evolutionary history of crops and also provides a framework to evaluate modern methods used to decipher the genetic mechanisms underlying phenotypic change. We then review these methods, framing the discussion in terms of the phenotype-genotype hierarchy. Top-down approaches, such as quantitative trait locus and linkage disequilibrium mapping, start with a phenotype of interest and use genetic analysis to identify candidate genes. Bottom-up approaches, alternatively, use population genetic analyses to identify potentially adaptive genes and then rely on standard bioinformatics and reverse genetic tools to connect selected genes to a phenotype. We discuss the successes, advantages, and challenges of each, but we conclude that bottom-up approaches to understanding domestication as an adaptive process hold greater promise both for the study of adaptation and as a means to identify genes that contribute to agronomically important traits.

Meiosis-driven genome variation in plants

Meiosis includes two successive divisions of the nucleus with one round of DNA replication and leads to the formation of gametes with half of the chromosomes of the mother cell during sexual reproduction. It provides a cytological basis for gametogenesis and nheritance in eukaryotes. Meiotic cell division is a complex and dynamic process that involves a number of molecular and cellular events, such as DNA and chromosome replication, chromosome pairing, synapsis and recombination, chromosome segregation, and cytokinesis. Meiosis maintains genome stability and integrity over sexual life cycles. On the other hand, meiosis generates genome variations in several ways. Variant meiotic recombination resulting from specific genome structures induces deletions, duplications, and other rearrangements within the genic and non-genic genomic regions and has been considered a major driving force for gene and genome evolution in nature. Meiotic abnormalities in chromosome segregation lead to chromosomally imbalanced gametes and aneuploidy. Meiotic restitution due to failure of the first or second meiotic division gives rise to unreduced gametes, which triggers polyploidization and genome expansion. This paper reviews research regarding meiosis-driven genome variation, including deletion and duplication of genomic regions, aneuploidy, and polyploidization, and discusses the effect of related meiotic events on genome variation and evolution in plants. Knowledge of various meiosis-driven genome variations provides insight into genome evolution and genetic variability in plants and facilitates plant genome research.

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.

An evolutionary view of human recombination

Recombination has essential functions in mammalian meiosis, which impose several constraints on the recombination process. However, recent studies have shown that, in spite of these roles, recombination rates vary tremendously among humans, and show marked differences between humans and closely related species. These findings provide important insights into the determinants of recombination rates and raise new questions about the selective pressures that affect recombination over different genomic scales, with implications for human genetics and evolutionary biology.