A study on the genetic relationships of Avena taxa and the origins of hexaploid oat.

KEY MESSAGE: Using next-generation DNA sequencing, it was possible to clarify the genetic relationships of Avena species and deduce the likely pathway from which hexaploid oat was formed by sequential polyploidization events. A sequence-based diversity study was conducted on a representative sample of accessions from species in the genus Avena using genotyping-by-sequencing technology. The results show that all Avena taxa can be assigned to one of four major genetic clusters: Cluster 1 = all hexaploids including cultivated oat, Cluster 2 = AC genome tetraploids, Cluster 3 = C genome diploids, Cluster 4 = A genome diploid and tetraploids. No evidence was found for the existence of discrete B or D genomes. Through a series of experiments involving the creation of in silico polyploids, it was possible to deduce that hexaploid oat likely formed by the fusion of an ancestral diploid species from Cluster 3 (A. clauda, A. eriantha) with an ancestral diploid species from Cluster 4D (A. longiglumis, A. canariensis, A. wiestii) to create the ancestral tetraploid from Cluster 2 (A. magna, A. murphyi, A. insularis). Subsequently, that ancestral tetraploid fused again with another ancestral diploid from Cluster 4D to create hexaploid oat. Based on the geographic distribution of these species, it is hypothesized that both the tetraploidization and hexaploidization events may have occurred in the region of northwest Africa, followed by radiation of hexaploid oat to its current worldwide distribution. The results from this study shed light not only on the origins of this important grain crop, but also have implications for germplasm collection and utilization in oat breeding.

Selfing for the design of genomic selection experiments in biparental plant populations.

Self-fertilization (selfing) is commonly used for population development in plant breeding, and it is well established that selfing increases genetic variance between lines, thus increasing response to phenotypic selection. Furthermore, numerous studies have explored how selfing can be deployed to maximal benefit in the context of traditional plant breeding programs (Cornish in Heredity 65:201-211,1990a, Heredity 65:213-220,1990b; Liu et al. in Theor Appl Genet 109:370-376, 2004; Pooni and Jinks in Heredity 54:255-260, 1985). However, the impact of selfing on response to genomic selection has not been explored. In the current study we examined how selfing impacts the two key aspects of genomic selection-GEBV prediction (training) and selection response. We reach the following conclusions: (1) On average, selfing increases genomic selection gains by more than 70 %. (2) The gains in genomic selection response attributable to selfing hold over a wide range population sizes (100-500), heritabilities (0.2-0.8), and selection intensities (0.01-0.1). However, the benefits of selfing are dramatically reduced as the number of QTLs drops below 20. (3) The major cause of the improved response to genomic selection with selfing is through an increase in the occurrence of superior genotypes and not through improved GEBV predictions. While performance of the training population improves with selfing (especially with low heritability and small population sizes), the magnitude of these improvements is relatively small compared with improvements observed in the selection population. To illustrate the value of these insights, we propose a practical genomic selection scheme that substantially shortens the number of generations required to fully capture the benefits of selfing. Specifically, we provide simulation evidence that indicates the proposed scheme matches or exceeds the selection gains observed in advanced populations (i.e. F 8 and doubled haploid) across a broad range of heritability and QTL models. Without sacrificing selection gains, we also predict that fully inbred candidates for potential commercialization can be identified as early as the F 4 generation.

The impact of recombination on short-term selection gain in plant breeding experiments.

Recombination is a requirement for response to selection, but researchers still debate whether increasing recombination beyond normal levels will result in significant gains in short-term selection. We tested this hypothesis, in the context of plant breeding, through a series of simulation experiments comparing short-term selection response (≤20 cycles) between populations with normal levels of recombination and similar populations with unconstrained recombination (i.e., free recombination). We considered additive and epistatic models and examined a wide range of values for key design variables: selection cycles, QTL number, heritability, linkage phase, selection intensity and population size. With few exceptions, going from normal to unconstrained levels of recombination produced only modest gains in response to selection (≈11 % on average). We then asked how breeders might capture some of this theoretical gain by increasing recombination through either (1) extra rounds of mating or (2) selection of highly recombinant individuals via use of molecular markers/maps. All methods tested captured less than half of the potential gain, but our analysis indicates that the most effective method is to select for increased recombination and the trait simultaneously. This recommendation is based on evidence of a favorable interaction between trait selection and the impact of recombination on selection gains. Finally, we examined the relative contributions of the two components of meiotic recombination, chromosome assortment and crossing over, to short-term selection gain. Depending primarily on the presence of trait selection pressure, chromosome assortment alone accounted for 40-75 % of gain in response to short-term selection.

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.

COSII genetic maps of two diploid Nicotiana species provide a detailed picture of synteny with tomato and insights into chromosome evolution in tetraploid N. tabacum.

Using single-copy conserved ortholog set (COSII) and simple sequence repeat (SSR) markers, we have constructed two genetic maps for diploid Nicotiana species, N. tomentosiformis and N. acuminata, respectively. N. acuminata is phylogenetically closer to N. sylvestris than to N. tomentosiformis, the latter two of which are thought to contribute the S-genome and T-genome, respectively, to the allotetraploid tobacco (N. tabacum L., 2n = 48). A comparison of the two maps revealed a minimum of seven inversions and one translocation subsequent to the divergence of these two diploid species. Further, comparing the diploid maps with a dense tobacco map revealed that the tobacco genome experienced chromosomal rearrangements more frequently than its diploid relatives, supporting the notion of accelerated genome evolution in allotetraploids. Mapped COSII markers permitted the investigation of Nicotiana-tomato syntenic relationships. A minimum of 3 (and up to 10) inversions and 11 reciprocal translocations differentiate the tomato genome from that of the last common ancestor of N. tomentosiformis and N. acuminata. Nevertheless, the marker/gene order is well preserved in 25 conserved syntenic segments. Molecular dating based on COSII sequences suggested that tobacco was formed 1.0 MYA or later. In conclusion, these COSII and SSR markers link the cultivated tobacco map to those of wild diploid Nicotiana species and tomato, thus providing a platform for cross-reference of genetic and genomic information among them as well as other solanaceous species including potato, eggplant, pepper and the closely allied coffee (Rubiaceae). Therefore they will facilitate genetic research in the genus Nicotiana.

A COSII genetic map of the pepper genome provides a detailed picture of synteny with tomato and new insights into recent chromosome evolution in the genus Capsicum.

We report herein the development of a pepper genetic linkage map which comprises 299 orthologous markers between the pepper and tomato genomes (including 263 conserved ortholog set II or COSII markers). The expected position of additional 288 COSII markers was inferred in the pepper map via pepper-tomato synteny, bringing the total orthologous markers in the pepper genome to 587. While pepper maps have been previously reported, this is the first complete map in the sense that all markers could be placed in 12 linkage groups corresponding to the 12 chromosomes. The map presented herein is relevant to the genomes of cultivated C. annuum and wild C. annuum (as well as related Capsicum species) which differ by a reciprocal chromosome translocation. This map is also unique in that it is largely based on COSII markers, which permits the inference of a detailed syntenic relationship between the pepper and tomato genomes-shedding new light on chromosome evolution in the Solanaceae. Since divergence from their last common ancestor is approximately 20 million years ago, the two genomes have become differentiated by a minimum number of 19 inversions and 6 chromosome translocations, as well as numerous putative single gene transpositions. Nevertheless, the two genomes share 35 conserved syntenic segments (CSSs) within which gene/marker order is well preserved. The high resolution COSII synteny map described herein provides a platform for cross-reference of genetic and genomic information (including the tomato genome sequence) between pepper and tomato and therefore will facilitate both applied and basic research in pepper.

A detailed synteny map of the eggplant genome based on conserved ortholog set II (COSII) markers.

We report herein the mapping of 115 PCR-based orthologous markers, including 110 conserved ortholog set or COSII markers, on the reference RFLP map of eggplant. The result permitted inference of a detailed syntenic relationship between the eggplant and tomato genomes. Further, the position of additional 522 COSII markers was inferred in the eggplant map via eggplant-tomato synteny, bringing the total number of markers in the eggplant genome to 869. Since divergence from their last common ancestor approximately 12 million years ago, the eggplant and tomato genomes have become differentiated by a minimum number of 24 inversions and 5 chromosomal translocations, as well as a number of single gene transpositions possibly triggered by transposable elements. Nevertheless, the two genomes share 37 conserved syntenic segments (CSSs) within which gene/marker order is well preserved. The high-resolution COSII synteny map described herein provides a platform for cross-reference of genetic and genomic information (including the tomato genome sequence) between eggplant and tomato and therefore will facilitate both applied and basic research in eggplant.

EST, COSII, and arbitrary gene markers give similar estimates of nucleotide diversity in cultivated tomato (Solanum lycopersicum L.).

Because cultivated tomato (Solanum lycopersicum L.) is low in genetic diversity, public, verified single nucleotide polymorphism (SNP) markers within the species are in demand. To promote marker development we resequenced approximately 23 kb in a diverse set of 31 tomato lines including TA496. Three classes of markers were sampled: (1) 26 expressed-sequence tag (EST), all of which were predicted to be polymorphic based on TA496, (2) 14 conserved ortholog set II (COSII) or unigene, and (3) ten published sequences, composed of nine fruit quality genes and one anonymous RFLP marker. The latter two types contained mostly noncoding DNA. In total, 154 SNPs and 34 indels were observed. The distributions of nucleotide diversity estimates among marker types were not significantly different from each other. Ascertainment bias of SNPs was evaluated for the EST markers. Despite the fact that the EST markers were developed using SNP prediction within a sample consisting of only one TA496 allele and one additional allele, the majority of polymorphisms in the 26 EST markers were represented among the other 30 tomato lines. Fifteen EST markers with published SNPs were more closely examined for bias. Mean SNP diversity observations were not significantly different between the original discovery sample of two lines (53 SNPs) and the 31 line diversity panel (56 SNPs). Furthermore, TA496 shared its haplotype with at least one other line at 11 of the 15 markers. These data demonstrate that public EST databases and noncoding regions are a valuable source of unbiased SNP markers in tomato.

Natural variation in an ABC transporter gene associated with seed size evolution in tomato species.

Seed size is a key determinant of evolutionary fitness in plants and is a trait that often undergoes tremendous changes during crop domestication. Seed size is most often quantitatively inherited, and it has been shown that Sw4.1 is one of the most significant quantitative trait loci (QTLs) underlying the evolution of seed size in the genus Solanum-especially in species related to the cultivated tomato. Using a combination of genetic, developmental, molecular, and transgenic techniques, we have pinpointed the cause of the Sw4.1 QTL to a gene encoding an ABC transporter gene. This gene exerts its control on seed size, not through the maternal plant, but rather via gene expression in the developing zygote. Phenotypic effects of allelic variation at Sw4.1 are manifested early in seed development at stages corresponding to the rapid deposition of starch and lipids into the endospermic cells. Through synteny, we have identified the Arabidopsis Sw4.1 ortholog. Mutagenesis has revealed that this ortholog is associated with seed length variation and fatty acid deposition in seeds, raising the possibility that the ABC transporter may modulate seed size variation in other species. Transcription studies show that the ABC transporter gene is expressed not only in seeds, but also in other tissues (leaves and roots) and, thus, may perform functions in parts of the plants other than developing seeds. Cloning and characterization of the Sw4.1 QTL gives new insight into how plants change seed during evolution and may open future opportunities for modulating seed size in crop plants for human purposes.

Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae.

Comparative genomics is a powerful tool for gaining insight into genomic function and evolution. However, in plants, sequence data that would enable detailed comparisons of both coding and noncoding regions have been limited in availability. Here we report the generation and analysis of sequences for an unduplicated conserved syntenic segment (CSS) in the genomes of five members of the agriculturally important plant family Solanaceae. This CSS includes a 105-kb region of tomato chromosome 2 and orthologous regions of the potato, eggplant, pepper, and petunia genomes. With a total neutral divergence of 0.73-0.78 substitutions/site, these sequences are similar enough that most noncoding regions can be aligned, yet divergent enough to be informative about evolutionary dynamics and selective pressures. The CSS contains 17 distinct genes with generally conserved order and orientation, but with numerous small-scale differences between species. Our analysis indicates that the last common ancestor of these species lived approximately 27-36 million years ago, that more than one-third of short genomic segments (5-15 bp) are under selection, and that more than two-thirds of selected bases fall in noncoding regions. In addition, we identify genes under positive selection and analyze hundreds of conserved noncoding elements. This analysis provides a window into 30 million years of plant evolution in the absence of polyploidization.

Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication.

Plant domestication represents an accelerated form of evolution, resulting in exaggerated changes in the tissues and organs of greatest interest to humans (for example, seeds, roots and tubers). One of the most extreme cases has been the evolution of tomato fruit. Cultivated tomato plants produce fruit as much as 1,000 times larger than those of their wild progenitors. Quantitative trait mapping studies have shown that a relatively small number of genes were involved in this dramatic transition, and these genes control two processes: cell cycle and organ number determination. The key gene in the first process has been isolated and corresponds to fw2.2, a negative regulator of cell division. However, until now, nothing was known about the molecular basis of the second process. Here, we show that the second major step in the evolution of extreme fruit size was the result of a regulatory change of a YABBY-like transcription factor (fasciated) that controls carpel number during flower and/or fruit development.

The SGN comparative map viewer.

MOTIVATION: With the rapid accumulation of genetic data for a multitude of different species, the availability of intuitive comparative genomic tools becomes an important requirement for the research community. Here we describe a web-based comparative viewer for mapping data, including genetic, physical and cytological maps, that is part of the SGN website (http://sgn.cornell.edu/) but that can also be installed and adapted for other websites. In addition to viewing and comparing different maps stored in the SGN database, the viewer allows users to upload their own maps and compare them to other maps in the system. The viewer is implemented in object oriented Perl, with a simple extensible interface to write data adapters for other relational database schemas and flat file formats.

Changes in regulation of a transcription factor lead to autogamy in cultivated tomatoes.

We report the cloning of Style2.1, the major quantitative trait locus responsible for a key floral attribute (style length) associated with the evolution of self-pollination in cultivated tomatoes. The gene encodes a putative transcription factor that regulates cell elongation in developing styles. The transition from cross-pollination to self-pollination was accompanied, not by a change in the STYLE2.1 protein, but rather by a mutation in the Style2.1 promoter that results in a down-regulation of Style2.1 expression during flower development.

Resistance to Phytophthora infestans in Lycopersicon pennellii.

To determine if the desert tomato, Lycopersicon pennellii, possesses resistance to late blight, caused by Phytophthora infestans, two plant populations were analyzed. Resistance was identified through assessments of disease progress in an F2 mapping population (L. esculentum × L. pennellii) and in a series of introgression lines (L. pennellii into L. esculentum). Levels of resistance varied widely among individuals within each population. However, the response of individuals to different strains of P. infestans was consistent. In the mapping population, a quantitative trait locus (QTL) was detected near marker T1556 on chromosome 6. This QTL accounted for 25% of the phenotypic variance in the population. The occurrence of this QTL was confirmed from analysis of the introgression lines (ILs), where IL 6-2 (containing marker T1556) was the most resistant IL in 2002 and the second most resistant IL in 2001. The identification of an additional QTL for resistance to late blight in tomato will aid in the development of durable resistance to this devastating disease.

Combining bioinformatics and phylogenetics to identify large sets of single-copy orthologous genes (COSII) for comparative, evolutionary and systematic studies: a test case in the euasterid plant clade.

We report herein the application of a set of algorithms to identify a large number (2869) of single-copy orthologs (COSII), which are shared by most, if not all, euasterid plant species as well as the model species Arabidopsis. Alignments of the orthologous sequences across multiple species enabled the design of "universal PCR primers," which can be used to amplify the corresponding orthologs from a broad range of taxa, including those lacking any sequence databases. Functional annotation revealed that these conserved, single-copy orthologs encode a higher-than-expected frequency of proteins transported and utilized in organelles and a paucity of proteins associated with cell walls, protein kinases, transcription factors, and signal transduction. The enabling power of this new ortholog resource was demonstrated in phylogenetic studies, as well as in comparative mapping across the plant families tomato (family Solanaceae) and coffee (family Rubiaceae). The combined results of these studies provide compelling evidence that (1) the ancestral species that gave rise to the core euasterid families Solanaceae and Rubiaceae had a basic chromosome number of x=11 or 12.2) No whole-genome duplication event (i.e., polyploidization) occurred immediately prior to or after the radiation of either Solanaceae or Rubiaceae as has been recently suggested.

EST database for early flower development in California poppy (Eschscholzia californica Cham., Papaveraceae) tags over 6,000 genes from a basal eudicot.

The Floral Genome Project (FGP) selected California poppy (Eschscholzia californica Cham. ssp. Californica) to help identify new florally-expressed genes related to floral diversity in basal eudicots. A large, non-normalized cDNA library was constructed from premeiotic and meiotic floral buds and sequenced to generate a database of 9,079 high quality Expressed Sequence Tags (ESTs). These sequences clustered into 5,713 unigenes, including 1,414 contigs and 4,299 singletons. Homologs of genes regulating many aspects of flower development were identified, including those for organ identity and development, cell and tissue differentiation, cell cycle control, and secondary metabolism. Over 5% of the transcriptome consisted of homologs to known floral gene families. Most are the first representatives of their respective gene families in basal eudicots and their conservation suggests they are important for floral development and/or function. App. 10% of the transcripts encoded transcription factors and other regulatory genes, including nine genes from the seven major lineages of the important MADS-box family of developmental regulators. Homologs of alkaloid pathway genes were also recovered, providing opportunities to explore adaptive evolution in secondary products. Furthermore, comparison of the poppy ESTs with the Arabidopsis genome provided support for putative Arabidopsis genes that previously lacked annotation. Finally, over 1,800 unique sequences had no observable homology in the public databases. The California poppy EST database and library will help bridge our understanding of flower initiation and development among higher eudicot and monocot model plants and provide new opportunities for comparative analysis of gene families across angiosperm species.

FW2.2 and cell cycle control in developing tomato fruit: a possible example of gene co-option in the evolution of a novel organ.

fw2.2 is one of the few QTLs thus far isolated from plants and the first one known to control fruit size. While it has been established that FW2.2 is a regulator (either directly or indirectly) of cell division, FW2.2 does not share sequence homology to any protein of known function (Frary et al. Science 289:85-88, 2000; Cong et al. Proc Natl Acad Sci USA 99:13606-13611, 2002; Liu et al. Plant Physiol 132:292-299, 2003). Thus, the mechanism by which FW2.2 mediates cell division in developing fruit is currently unknown. In an effort to remedy this situation, a combination of yeast two-hybrid screens, in vitro binding assays and cell bombardment studies were performed. The results provide strong evidence that FW2.2 physically interacts at or near the plasma membrane with the regulatory (beta) subunit of a CKII kinase. CKII kinases are well-studied in both yeast and animals where they form part of cell cycle related signaling pathway. Thus while FW2.2 is a plant-specific protein and regulates cell division in a specialized plant organ (fruit), it appears to participate in a cell-cycle control signal transduction pathway that predates the divergence of single- and multi-cellular organisms. These results thus provide a glimpse into how ancient and conserved regulatory processes can be co-opted in the evolution of novel organs such as fruit.

CpG + CpNpG analysis of protein-coding sequences from tomato.

We develop codon-based models for simultaneously inferring the mutational effects of CpG and CpNpG methylation in coding regions. In a data set of 369 tomato genes, we show that there is very little effect of CpNpG methylation but a strong effect of CpG methylation affecting almost all genes. We further show that the CpNpG and CpG effects are largely uncorrelated. Our results suggest different roles of CpG and CpNpG methylation, with CpNpG methylation possibly playing a specialized role in defense against transposons and RNA viruses.

Coffee and tomato share common gene repertoires as revealed by deep sequencing of seed and cherry transcripts.

An EST database has been generated for coffee based on sequences from approximately 47,000 cDNA clones derived from five different stages/tissues, with a special focus on developing seeds. When computationally assembled, these sequences correspond to 13,175 unigenes, which were analyzed with respect to functional annotation, expression profile and evolution. Compared with Arabidopsis, the coffee unigenes encode a higher proportion of proteins related to protein modification/turnover and metabolism-an observation that may explain the high diversity of metabolites found in coffee and related species. Several gene families were found to be either expanded or unique to coffee when compared with Arabidopsis. A high proportion of these families encode proteins assigned to functions related to disease resistance. Such families may have expanded and evolved rapidly under the intense pathogen pressure experienced by a tropical, perennial species like coffee. Finally, the coffee gene repertoire was compared with that of Arabidopsis and Solanaceous species (e.g. tomato). Unlike Arabidopsis, tomato has a nearly perfect gene-for-gene match with coffee. These results are consistent with the facts that coffee and tomato have a similar genome size, chromosome karyotype (tomato, n=12; coffee n=11) and chromosome architecture. Moreover, both belong to the Asterid I clade of dicot plant families. Thus, the biology of coffee (family Rubiacaeae) and tomato (family Solanaceae) may be united into one common network of shared discoveries, resources and information.

Transcriptome and selected metabolite analyses reveal multiple points of ethylene control during tomato fruit development.

Transcriptome profiling via cDNA microarray analysis identified 869 genes that are differentially expressed in developing tomato (Solanum lycopersicum) pericarp. Parallel phenotypic and targeted metabolite comparisons were employed to inform the expression analysis. Transcript accumulation in tomato fruit was observed to be extensively coordinated and often completely dependent on ethylene. Mutation of an ethylene receptor (Never-ripe [Nr]), which reduces ethylene sensitivity and inhibits ripening, alters the expression of 37% of these 869 genes. Nr also influences fruit morphology, seed number, ascorbate accumulation, carotenoid biosynthesis, ethylene evolution, and the expression of many genes during fruit maturation, indicating that ethylene governs multiple aspects of development both prior to and during fruit ripening in tomato. Of the 869 genes identified, 628 share homology (E-value < or = 1 x 10(-10)) with known gene products or known protein domains. Of these 628 loci, 72 share homology with previously described signal transduction or transcription factors, suggesting complex regulatory control. These results demonstrate multiple points of ethylene regulatory control during tomato fruit development and provide new insights into the molecular basis of ethylene-mediated ripening.