Euchromatin and pericentromeric heterochromatin: comparative composition in the tomato genome

Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity

Soil salinization is a serious problem for cultivation of rice, as among cereals rice is the most salt sensitive crop, and more than 40% of the total agricultural land amounting to approximately 80 million ha the world over is salt affected. Salinity affects a plant in a varieties of ways, including ion toxicity, osmotic stress and oxidative damage. Since miRNAs occupy the top place in biochemical events determining a trait, understanding their role in salt tolerance is highly desirable, which may allow introduction of the trait in the rice cultivars of choice through biotechnological interventions. High throughput sequencing of sRNAs in the root and shoot tissues of the seedlings of the control and NaCl treated Pokkali, a salt-tolerant rice variety, identified 75 conserved miRNAs and mapped 200 sRNAs to the rice genome as novel miRNAs. Expression of nine novel miRNAs and two conserved miRNAs were confirmed by Northern blotting. Several of both conserved and novel miRNAs that expressed differentially in root and/or shoot tissues targeted transcription factors like AP2/EREBP domain protein, ARF, NAC, MYB, NF-YA, HD-Zip III, TCP and SBP reported to be involved in salt tolerance or in abiotic stress tolerance in general. Most of the novel miRNAs expressed in the salt tolerant wild rice Oryza coarctata, suggesting conservation of miRNAs in taxonomically related species. One of the novel miRNAs, osa-miR12477, also targeted L-ascorbate oxidase (LAO), indicating build-up of oxidative stress in the plant upon salt treatment, which was confirmed by DAB staining. Thus, salt tolerance might involve miRNA-mediated regulation of 1) cellular abundance of the hormone signaling components like EREBP and ARF, 2) synthesis of abiotic stress related transcription factors, and 3) antioxidative component like LAO for mitigation of oxidative damage. The study clearly indicated importance of osa-miR12477 regulated expression of LAO in salt tolerance in the plant.

LTR-TEs abundance, timing and mobility in Solanum commersonii and S. tuberosum genomes following cold-stress conditions

Copia/Ale is the youngest lineage in both Solanum tuberosum and S. commersonii. Within it, we identified nightshade, a new LTR element active in the cultivated potato. From an evolutionary perspective, long-terminal repeat retrotransposons (LTR-RT) activity during stress may be viewed as a mean by which organisms can keep up rates of genetic adaptation to changing conditions. Potato is one of the most important crop consumed worldwide, but studies on LTR-RT characterization are still lacking. Here, we assessed the abundance, insertion time and activity of LTR-RTs in both cultivated Solanum tuberosum and its cold-tolerant wild relative S. commersonii genomes. Gypsy elements were more abundant than Copia ones, suggesting that the former was somehow more successful in colonizing potato genomes. However, Copia elements, and in particular, the Ale lineage, are younger than Gypsy ones, since their insertion time was in average ~ 2 Mya. Due to the ability of LTR-RTs to be circularized by the host DNA repair mechanisms, we identified via mobilome-seq a Copia/Ale element (called nightshade, informal name used for potato family) active in S. tuberosum genome. Our analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in potato.

Environmental and epigenetic regulation of Rider retrotransposons in tomato

Transposable elements in crop plants are the powerful drivers of phenotypic variation that has been selected during domestication and breeding programs. In tomato, transpositions of the LTR (long terminal repeat) retrotransposon family Rider have contributed to various phenotypes of agronomical interest, such as fruit shape and colour. However, the mechanisms regulating Rider activity are largely unknown. We have developed a bioinformatics pipeline for the functional annotation of retrotransposons containing LTRs and defined all full-length Rider elements in the tomato genome. Subsequently, we showed that accumulation of Rider transcripts and transposition intermediates in the form of extrachromosomal DNA is triggered by drought stress and relies on abscisic acid signalling. We provide evidence that residual activity of Rider is controlled by epigenetic mechanisms involving siRNAs and the RNA-dependent DNA methylation pathway. Finally, we demonstrate the broad distribution of Rider-like elements in other plant species, including crops. Our work identifies Rider as an environment-responsive element and a potential source of genetic and epigenetic variation in plants.

Transposons are major constituents of plant genomes and represent a powerful source of internal genetic and epigenetic variation. For example, domestication of maize has been facilitated by a dramatic change in plant architecture, the consequence of a transposition event. Insertion of transposons near genes often confers quantitative phenotypic variation linked to changes in transcriptional patterns, as documented for blood oranges and grapes. In tomato, the most widely grown fruit crop and model for fleshy fruit biology, occurrences of several beneficial traits related to fruit shape and plant architecture are due to the activity of the transposon family Rider. While Rider represents a unique endogenous source of genetic and epigenetic variation, mechanisms regulating Rider activity remain unexplored. By achieving experimentally-controlled activation of the Rider family, we shed light on the regulation of these transposons by drought stress, signalling by phytohormones, as well as epigenetic pathways. Furthermore, we reveal the presence of Rider-like elements in other economically important crops such as rapeseed, beetroot and quinoa. This suggests that drought-inducible Rider activation could be further harnessed to generate genetic and epigenetic variation for crop breeding, and highlights the potential of transposon-directed mutagenesis for crop improvement.

Genome-wide nucleotide patterns and potential mechanisms of genome divergence following domestication in maize and soybean

BACKGROUND

Plant domestication provides a unique model to study genome evolution. Many studies have been conducted to examine genes, genetic diversity, genome structure, and epigenome changes associated with domestication. Interestingly, domesticated accessions have significantly higher [A] and [T] values across genome-wide polymorphic sites than accessions sampled from the corresponding progenitor species. However, the relative contributions of different genomic regions to this genome divergence pattern and underlying mechanisms have not been well characterized.

RESULTS

Here, we investigate the genome-wide base-composition patterns by analyzing millions of SNPs segregating among 100 accessions from a teosinte-maize comparison set and among 302 accessions from a wild-domesticated soybean comparison set. We show that non-genic part of the genome has a greater contribution than genic SNPs to the [AT]-increase observed between wild and domesticated accessions in maize and soybean. The separation between wild and domesticated accessions in [AT] values is significantly enlarged in non-genic and pericentromeric regions. Motif frequency and sequence context analyses show the motifs (PyCG) related to solar-UV signature are enriched in these regions, particularly when they are methylated. Additional analysis using population-private SNPs also implicates the role of these motifs in relatively recent mutations. With base-composition across polymorphic sites as a genome phenotype, genome scans identify a set of putative candidate genes involved in UV damage repair pathways.

CONCLUSIONS

The [AT]-increase is more pronounced in genomic regions that are non-genic, pericentromeric, transposable elements; methylated; and with low recombination. Our findings establish important links among UV radiation, mutation, DNA repair, methylation, and genome evolution.

Comparative analysis of repetitive sequences among species from the potato and the tomato clades

BACKGROUND AND AIMS

The genus Solanum includes important vegetable crops and their wild relatives. Introgression of their useful traits into elite cultivars requires effective recombination between hom(e)ologues, which is partially determined by genome sequence differentiation. In this study we compared the repetitive genome fractions of wild and cultivated species of the potato and tomato clades in a phylogenetic context.

METHODS

Genome skimming followed by a clustering approach was used as implemented in the RepeatExplorer pipeline. Repeat classes were annotated and the sequences of their main domains were compared.

KEY RESULTS

Repeat abundance and genome size were correlated and the larger genomes of species in the tomato clade were found to contain a higher proportion of unclassified elements. Families and lineages of repetitive elements were largely conserved between the clades, but their relative proportions differed. The most abundant repeats were Ty3/Gypsy elements. Striking differences in abundance were found in the highly dynamic Ty3/Gypsy Chromoviruses and Ty1/Copia Tork elements. Within the potato clade, early branching Solanum cardiophyllum showed a divergent repeat profile. There were also contrasts between cultivated and wild potatoes, mostly due to satellite amplification in the cultivated species. Interspersed repeat profiles were very similar among potatoes. The repeat profile of Solanum etuberosum was more similar to that of the potato clade.

CONCLUSIONS

The repeat profiles in Solanum seem to be very similar despite genome differentiation at the level of collinearity. Removal of transposable elements by unequal recombination may have been responsible for structural rearrangements across the tomato clade. Sequence variability in the tomato clade is congruent with clade-specific amplification of repeats after its divergence from S. etuberosum and potatoes. The low differentiation among potato and its wild relatives at the level of interspersed repeats may explain the difficulty in discriminating their genomes by genomic in situ hybridization techniques.

Retrotranspositional landscape of Asian rice revealed by 3000 genomes

The recent release of genomic sequences for 3000 rice varieties provides access to the genetic diversity at species level for this crop. We take advantage of this resource to unravel some features of the retrotranspositional landscape of rice. We develop software TRACKPOSON specifically for the detection of transposable elements insertion polymorphisms (TIPs) from large datasets. We apply this tool to 32 families of retrotransposons and identify more than 50,000 TIPs in the 3000 rice genomes. Most polymorphisms are found at very low frequency, suggesting that they may have occurred recently in agro. A genome-wide association study shows that these activations in rice may be triggered by external stimuli, rather than by the alteration of genetic factors involved in transposable element silencing pathways. Finally, the TIPs dataset is used to trace the origin of rice domestication. Our results suggest that rice originated from three distinct domestication events.

Low-coverage resequencing detects meiotic recombination pattern and features in tomato RILs

Traditional plant breeding relies on meiotic recombination for mixing of parental alleles to create novel allele combinations. Detailed analysis of recombination patterns in model organisms shows that recombination is tightly regulated within the genome, but frequencies vary extensively along chromosomes. Despite being a model organism for fruit developmental studies, high-resolution recombination patterns are lacking in tomato. In this study, we developed a novel methodology to use low-coverage resequencing to identify genome-wide recombination patterns and applied this methodology on 60 tomato Recombinant Inbred Lines (RILs). Our methodology identifies polymorphic markers from the low-coverage resequencing population data and utilizes the same data to locate the recombination breakpoints in individuals by using a variable sliding window. We identified 1,445 recombination sites comprising 112 recombination prone regions enriched for AT-rich DNA motifs. Furthermore, the recombination prone regions in tomato preferably occurred in gene promoters over intergenic regions, an observation consistent with Arabidopsis thaliana, Zea mays and Mimulus guttatus. Overall, our cost effective method and findings enhance the understanding of meiotic recombination in tomato and suggest evolutionarily conserved recombination associated genomic features.

A collection of enhancer trap insertional mutants for functional genomics in tomato

With the completion of genome sequencing projects, the next challenge is to close the gap between gene annotation and gene functional assignment. Genomic tools to identify gene functions are based on the analysis of phenotypic variations between a wild type and its mutant; hence, mutant collections are a valuable resource. In this sense, T-DNA collections allow for an easy and straightforward identification of the tagged gene, serving as the basis of both forward and reverse genetic strategies. This study reports on the phenotypic and molecular characterization of an enhancer trap T-DNA collection in tomato (Solanum lycopersicum L.), which has been produced by Agrobacterium-mediated transformation using a binary vector bearing a minimal promoter fused to the uidA reporter gene. Two genes have been isolated from different T-DNA mutants, one of these genes codes for a UTP-glucose-1-phosphate uridylyltransferase involved in programmed cell death and leaf development, which means a novel gene function reported in tomato. Together, our results support that enhancer trapping is a powerful tool to identify novel genes and regulatory elements in tomato and that this T-DNA mutant collection represents a highly valuable resource for functional analyses in this fleshy-fruited model species.

Reverse transcriptase sequences from mulberry LTR retrotransposons: characterization analysis

Copia and Gypsy play important roles in structural, functional and evolutionary dynamics of plant genomes. In this study, a total of 106 and 101, Copia and Gypsy reverse transcriptase (rt) were amplified respectively in the Morus notabilis genome using degenerate primers. All sequences exhibited high levels of heterogeneity, were rich in AT and possessed higher sequence divergence of Copia rt in comparison to Gypsy rt. Two reasons are likely to account for this phenomenon: a) these elements often experience deletions or fragmentation by illegitimate or unequal homologous recombination in the transposition process; b) strong purifying selective pressure drives the evolution of these elements through “selective silencing” with random mutation and eventual deletion from the host genome. Interestingly, mulberry rt clustered with other rt from distantly related taxa according to the phylogenetic analysis. This phenomenon did not result from horizontal transposable element transfer. Results obtained from fluorescence in situ hybridization revealed that most of the hybridization signals were preferentially concentrated in pericentromeric and distal regions of chromosomes, and these elements may play important roles in the regions in which they are found. Results of this study support the continued pursuit of further functional studies of Copia and Gypsy in the mulberry genome.

Development and Genetic Characterization of Advanced Backcross Materials and An Introgression Line Population of Solanum incanum in a S. melongena Background

Advanced backcrosses (ABs) and introgression lines (ILs) of eggplant (Solanum melongena) can speed up genetics and genomics studies and breeding in this crop. We have developed the first full set of ABs and ILs in eggplant using Solanum incanum, a wild eggplant that has a relatively high tolerance to drought, as a donor parent. The development of these ABs and IL eggplant populations had a low efficiency in the early stages, because of the lack of molecular markers and genomic tools. However, this dramatically improved after performing genotyping-by-sequencing in the first round of selfing, followed by high-resolution-melting single nucleotide polymorphism genotyping in subsequent selection steps. A set of 73 selected ABs covered 99% of the S. incanum genome, while 25 fixed immortal ILs, each carrying a single introgressed fragment in homozygosis, altogether spanned 61.7% of the S. incanum genome. The introgressed size fragment in the ILs contained between 0.1 and 10.9% of the S. incanum genome, with a mean value of 4.3%. Sixty-eight candidate genes involved in drought tolerance were identified in the set of ILs. This first set of ABs and ILs of eggplant will be extremely useful for the genetic dissection of traits of interest for eggplant, and represents an elite material for introduction into the breeding pipelines for developing new eggplant cultivars adapted to the challenges posed by the climate-change scenario.

Condensation patterns of prophase/prometaphase chromosome are correlated with H4K5 histone acetylation and genomic DNA contents in plants

Mitotic prophase chromosome condensation plays an essential role in nuclear division being therefore regulated by highly conserved mechanisms. However, degrees of chromatin condensation in prophase-prometaphase cells may vary along the chromosomes resulting in specific condensation patterns. We examined different condensation patterns (CPs) of prophase and prometaphase chromosomes and investigated their relationship with genome size and distribution of histone H4 acetylated at lysine 5 (H4K5ac) in 17 plant species. Our results showed that most species with small genomes (2C < 5 pg) (Arachis pusilla, Bixa orellana, Costus spiralis, Eleutherine bulbosa, Indigofera campestris, Phaseolus lunatus, P. vulgaris, Poncirus trifoliata, and Solanum lycopersicum) displayed prophase chromosomes with late condensing terminal regions that were highly enriched in H4K5ac, and early condensing regions with apparently non-acetylated proximal chromatin. The species with large genomes (Allium cepa, Callisia repens, Araucaria angustifolia and Nothoscordum pulchellum) displayed uniformly condensed and acetylated prophase/prometaphase chromosomes. Three species with small genomes (Eleocharis geniculata, Rhynchospora pubera, and R. tenuis) displayed CP and H4K5ac labeling patterns similar to species with large genomes, whereas a forth species (Emilia sonchifolia) exhibited a gradual chromosome labeling, being more acetylated in the terminal regions and less acetylated in the proximal ones. The nucleolus organizer chromatin was the only chromosomal region that in prometaphase or metaphase could be hyperacetylated, hypoacetylated or non-acetylated, depending on the species. Our data indicate that the CP of a plant chromosome complement is influenced but not exclusively determined by nuclear and chromosomal DNA contents, whereas the CP of individual chromosomes is clearly correlated with H4K5ac distribution.

Diversity, distribution and dynamics of full-length Copia and Gypsy LTR retroelements in Solanum lycopersicum

Transposable elements are the most abundant components of plant genomes and can dramatically induce genetic changes and impact genome evolution. In the recently sequenced genome of tomato (Solanum lycopersicum), the estimated fraction of elements corresponding to retrotransposons is nearly 62%. Given that tomato is one of the most important vegetable crop cultivated and consumed worldwide, understanding retrotransposon dynamics can provide insight into its evolution and domestication processes. In this study, we performed a genome-wide in silico search of full-length LTR retroelements in the tomato nuclear genome and annotated 736 full-length Gypsy and Copia retroelements. The dispersion level across the 12 chromosomes, the diversity and tissue-specific expression of those elements were estimated. Phylogenetic analysis based on the retrotranscriptase region revealed the presence of 12 major lineages of LTR retroelements in the tomato genome. We identified 97 families, of which 77 and 20 belong to the superfamilies Copia and Gypsy, respectively. Each retroelement family was characterized according to their element size, relative frequencies and insertion time. These analyses represent a valuable resource for comparative genomics within the Solanaceae, transposon-tagging and for the design of cultivar-specific molecular markers in tomato.

Two Loci Contribute Epistastically to Heterospecific Pollen Rejection, a Postmating Isolating Barrier Between Species

Recognition and rejection of heterospecific male gametes occurs in a broad range of taxa, although the complexity of mechanisms underlying these components of postmating cryptic female choice is poorly understood. In plants, the arena for postmating interactions is the female reproductive tract (pistil), within which heterospecific pollen tube growth can be arrested via active molecular recognition and rejection. Unilateral incompatibility (UI) is one such postmating barrier in which pollen arrest occurs in only one direction of an interspecific cross. We investigated the genetic basis of pistil-side UI between Solanum species, with the specific goal of understanding the role and magnitude of epistasis between UI QTL. Using heterospecific introgression lines (ILs) between Solanum pennellii and S. lycopersicum, we assessed the individual and pairwise effects of three chromosomal regions (ui1.1, ui3.1, and ui12.1) previously associated with interspecific UI among Solanum species. Specifically, we generated double introgression (‘pyramided’) genotypes that combined ui12.1 with each of ui1.1 and ui3.1, and assessed the strength of UI pollen rejection in the pyramided lines, compared to single introgression genotypes. We found that none of the three QTL individually showed UI rejection phenotypes, but lines combining ui3.1 and ui12.1 showed significant pistil-side pollen rejection. Furthermore, double ILs (DILs) that combined different chromosomal regions overlapping ui3.1 differed significantly in their rate of UI, consistent with at least two genetic factors on chromosome three contributing quantitatively to interspecific pollen rejection. Together, our data indicate that loci on both chromosomes 3 and 12 are jointly required for the expression of UI between S. pennellii and S. lycopersicum, suggesting that coordinated molecular interactions among a relatively few loci underlie the expression of this postmating prezygotic barrier. In addition, in conjunction with previous data, at least one of these loci appears to also contribute to conspecific self-incompatibility (SI), consistent with a partially shared genetic basis between inter- and intraspecific mechanisms of postmating prezygotic female choice.

Collinearity between potato (Solanum tuberosum L.) and wild relatives assessed by comparative cytogenetic mapping

A major bottleneck to introgressive hybridization is the lack of genome collinearity between the donor (alien) genome and the recipient crop genome. Structural differences between the homeologs may create unbalanced segregation of chromosomes or cause linkage drag. To assess large-scale collinearity between potato and two of its wild relatives (Solanum commersonii and Solanum chacoense), we used BAC-FISH mapping of sequences with known positions on the RH potato map. BAC probes could successfully be hybridized to the S. commersonii and S. chachoense pachytene chromosomes, confirming their correspondence with linkage groups in RH potato. Our study shows that the order of BAC signals is conserved. Distances between BAC signals were quantified and compared; some differences found suggest either small-scale rearrangements or reduction/amplification of repeats. We conclude that S. commersonii and S. chacoense are collinear with cultivated Solanum tuberosum on the whole chromosome scale, making these amenable species for efficient introgressive hybridization breeding.

Advancing Eucalyptus Genomics: Cytogenomics Reveals Conservation of Eucalyptus Genomes

The genus Eucalyptus encloses several species with high ecological and economic value, being the subgenus Symphyomyrtus one of the most important. Species such as E. grandis and E. globulus are well characterized at the molecular level but knowledge regarding genome and chromosome organization is very scarce. Here we characterized and compared the karyotypes of three economically important species, E. grandis, E. globulus, and E. calmadulensis, and three with ecological relevance, E. pulverulenta, E. cornuta, and E. occidentalis, through an integrative approach including genome size estimation, fluorochrome banding, rDNA FISH, and BAC landing comprising genes involved in lignin biosynthesis. All karyotypes show a high degree of conservation with pericentromeric 35S and 5S rDNA loci in the first and third pairs, respectively. GC-rich heterochromatin was restricted to the 35S rDNA locus while the AT-rich heterochromatin pattern was species-specific. The slight differences in karyotype formulas and distribution of AT-rich heterochromatin, along with genome sizes estimations, support the idea of Eucalyptus genome evolution by local expansions of heterochromatin clusters. The unusual co-localization of both rDNA with AT-rich heterochromatin was attributed mainly to the presence of silent transposable elements in those loci. The cinnamoyl CoA reductase gene (CCR1) previously assessed to linkage group 10 (LG10) was clearly localized distally at the long arm of chromosome 9 establishing an unexpected correlation between the cytogenetic chromosome 9 and the LG10. Our work is novel and contributes to the understanding of Eucalyptus genome organization which is essential to develop successful advanced breeding strategies for this genus.

Highly distinct chromosomal structures in cowpea (Vigna unguiculata), as revealed by molecular cytogenetic analysis

Cowpea (Vigna unguiculata (L.) Walp) is an important legume, particularly in developing countries. However, little is known about its genome or chromosome structure. We used molecular cytogenetics to characterize the structure of pachytene chromosomes to advance our knowledge of chromosome and genome organization of cowpea. Our data showed that cowpea has highly distinct chromosomal structures that are cytologically visible as brightly DAPI-stained heterochromatic regions. Analysis of the repetitive fraction of the cowpea genome present at centromeric and pericentromeric regions confirmed that two retrotransposons are major components of pericentromeric regions and that a 455-bp tandem repeat is found at seven out of 11 centromere pairs in cowpea. These repeats likely evolved after the divergence of cowpea from common bean and form chromosomal structure unique to cowpea. The integration of cowpea genetic and physical chromosome maps reveals potential regions of suppressed recombination due to condensed heterochromatin and a lack of pairing in a few chromosomal termini. This study provides fundamental knowledge on cowpea chromosome structure and molecular cytogenetics tools for further chromosome studies.

Multiple QTL for horticultural traits and quantitative resistance to Phytophthora infestans linked on Solanum habrochaites chromosome 11

Previously, a Phytophthora infestans resistance QTL from Solanum habrochaites chromosome 11 was introgressed into cultivated tomato (S. lycopersicum). Fine mapping of this resistance QTL using near-isogenic lines (NILs) revealed some co-located QTL with undesirable effects on plant size, canopy density, and fruit size traits. Subsequently, higher-resolution mapping with sub-NILs detected multiple P. infestans resistance QTL within this 9.4-cM region of chromosome 11. In our present study, these same sub-NILs were also evaluated for 17 horticultural traits, including yield, maturity, fruit size and shape, fruit quality, and plant architecture traits in replicated field experiments over 2 years. The horticultural trait QTL originally detected by fine mapping each fractionated into two or more QTL at higher resolution. A total of 34 QTL were detected across all traits, with 14% exhibiting significant QTL × environment interactions (QTL × E). QTL for many traits were co-located, suggesting either pleiotropic effects or tight linkage among genes controlling these traits. Recombination in the pericentromeric region of the introgression between markers TG147 and At4g10050 was suppressed to approximately 29.7 Mbp per cM, relative to the genomewide average of 750 kbp per cM. The genetic architecture of many of the horticultural and P. infestans resistance traits that mapped within this chromosome 11 S. habrochaites region is complex. Complicating factors included fractionation of QTL, pleiotropy or tight linkage of QTL for multiple traits, pericentromeric chromosomal location(s), and/or QTL × E. High-resolution mapping of QTL in this region would be needed to determine which specific target QTL could be useful in breeding cultivated tomato.

Young but not relatively old retrotransposons are preferentially located in gene-rich euchromatic regions in tomato (Solanum lycopersicum) plants

Long terminal repeat (LTR) retrotransposons are the major DNA components of flowering plants. They are generally enriched in pericentromeric heterochromatin regions of their host genomes, which could result from the preferential insertion of LTR retrotransposons and the low effectiveness of purifying selection in these regions. To estimate the relative importance of the actions of these two factors on their distribution pattern, the LTR retrotransposons in Solanum lycopersicum (tomato) plants were characterized at the genome level, and then the distribution of young elements was compared with that of relatively old elements. The current data show that old elements are mainly located in recombination-suppressed heterochromatin regions, and that young elements are preferentially located in the gene-rich euchromatic regions. Further analysis showed a negative correlation between the insertion time of LTR retrotransposons and the recombination rate. The data also showed there to be more solo LTRs in genic regions than in intergenic regions or in regions close to genes. These observations indicate that, unlike in many other plant genomes, the current LTR retrotransposons in tomatoes have a tendency to be preferentially located into euchromatic regions, probably caused by their severe suppression of activities in heterochromatic regions. These elements are apt to be maintained in heterochromatin regions, probably as a consequence of the pericentromeric effect in tomatoes. These results also indicate that local recombination rates and intensities of purifying selection in different genomic regions are largely responsible for structural variation and non-random distribution of LTR retrotransposons in tomato plants.

Variation block-based genomics method for crop plants

BACKGROUND

In contrast with wild species, cultivated crop genomes consist of reshuffled recombination blocks, which occurred by crossing and selection processes. Accordingly, recombination block-based genomics analysis can be an effective approach for the screening of target loci for agricultural traits.

RESULTS

We propose the variation block method, which is a three-step process for recombination block detection and comparison. The first step is to detect variations by comparing the short-read DNA sequences of the cultivar to the reference genome of the target crop. Next, sequence blocks with variation patterns are examined and defined. The boundaries between the variation-containing sequence blocks are regarded as recombination sites. All the assumed recombination sites in the cultivar set are used to split the genomes, and the resulting sequence regions are termed variation blocks. Finally, the genomes are compared using the variation blocks. The variation block method identified recurring recombination blocks accurately and successfully represented block-level diversities in the publicly available genomes of 31 soybean and 23 rice accessions. The practicality of this approach was demonstrated by the identification of a putative locus determining soybean hilum color.

CONCLUSIONS

We suggest that the variation block method is an efficient genomics method for the recombination block-level comparison of crop genomes. We expect that this method will facilitate the development of crop genomics by bringing genomics technologies to the field of crop breeding.

Fluorescence in situ hybridization and optical mapping to correct scaffold arrangement in the tomato genome

The order and orientation (arrangement) of all 91 sequenced scaffolds in the 12 pseudomolecules of the recently published tomato (Solanum lycopersicum, 2n = 2x = 24) genome sequence were positioned based on marker order in a high-density linkage map. Here, we report the arrangement of these scaffolds determined by two independent physical methods, bacterial artificial chromosome–fluorescence in situ hybridization (BAC-FISH) and optical mapping. By localizing BACs at the ends of scaffolds to spreads of tomato synaptonemal complexes (pachytene chromosomes), we showed that 45 scaffolds, representing one-third of the tomato genome, were arranged differently than predicted by the linkage map. These scaffolds occur mostly in pericentric heterochromatin where 77% of the tomato genome is located and where linkage mapping is less accurate due to reduced crossing over. Although useful for only part of the genome, optical mapping results were in complete agreement with scaffold arrangement by FISH but often disagreed with scaffold arrangement based on the linkage map. The scaffold arrangement based on FISH and optical mapping changes the positions of hundreds of markers in the linkage map, especially in heterochromatin. These results suggest that similar errors exist in pseudomolecules from other large genomes that have been assembled using only linkage maps to predict scaffold arrangement, and these errors can be corrected using FISH and/or optical mapping. Of note, BAC-FISH also permits estimates of the sizes of gaps between scaffolds, and unanchored BACs are often visualized by FISH in gaps between scaffolds and thus represent starting points for filling these gaps.

Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species

Hot pepper (Capsicum annuum), one of the oldest domesticated crops in the Americas, is the most widely grown spice crop in the world. We report whole-genome sequencing and assembly of the hot pepper (Mexican landrace of Capsicum annuum cv. CM334) at 186.6× coverage. We also report resequencing of two cultivated peppers and de novo sequencing of the wild species Capsicum chinense. The genome size of the hot pepper was approximately fourfold larger than that of its close relative tomato, and the genome showed an accumulation of Gypsy and Caulimoviridae family elements. Integrative genomic and transcriptomic analyses suggested that change in gene expression and neofunctionalization of capsaicin synthase have shaped capsaicinoid biosynthesis. We found differential molecular patterns of ripening regulators and ethylene synthesis in hot pepper and tomato. The reference genome will serve as a platform for improving the nutritional and medicinal values of Capsicum species.

Sequence-based SSR marker development and their application in defining the Introgressions of LA0716 (Solanum pennellii) in the background of cv. M82 (Solanum lycopersicum)

The introgression lines (ILs) from cv. M82 (Solanum lycopersicum) × LA0716 (S. pennellii) have been proven to be exceptionally useful for genetic analysis and gene cloning. The introgressions were originally defined by RFLP markers at their development. The objectives of this study are to develop polymorphic SSR markers, and to re-define the DNA introgression from LA0716 in the ILs. Tomato sequence data was scanned by software to generate SSR markers. In total, 829 SSRs, which could be robustly amplified by PCR, were developed. Among them, 658 SSRs were dinucleotide repeats, 162 were trinucleotide repeats, and nine were tetranucleotide repeats. The 829 SSRs together with 96 published RFLPs were integrated into the physical linkage map of S. lycopersicum. Introgressions of DNA fragments from LA0716 were re-defined among the 75 ILs using the newly developed SSRs. A specific introgression of DNA fragment from LA0716 was identified in 72 ILs as described previously by RFLP, whereas the specific DNA introgression described previously were not detected in the ILs LA4035, LA4059 and LA4091. The physical location of each investigated DNA introgression was finely determined by SSR mapping. Among the 72 ILs, eight ILs showed a shorter and three ILs (IL3-2, IL12-3 and IL12-3-1) revealed a longer DNA introgression than that framed by RFLPs. Furthermore, 54 previously undefined segments were found in 21 ILs, ranging from 1 to 11 DNA introgressions per IL. Generally, the newly developed SSRs provide additional markers for genetic studies of tomatoes, and the fine definition of DNA introgressions from LA0716 would facilitate the use of the ILs for genetic analysis and gene cloning.

Mapping pachytene chromosomes of coffee using a modified protocol for fluorescence in situ hybridization

Fluorescence in situ hybridization (FISH) is the most direct method for physically mapping DNA sequences on chromosomes. Fluorescence in situ hybridization mapping of meiotic chromosomes during the pachytene stage is an important tool in plant cytogenetics, because it provides high-resolution measurements of physical distances. Fluorescence in situ hybridization mapping of coffee pachytene chromosomes offers significant advantages compared with FISH mapping of somatic chromosomes, because pachytene chromosomes are 30 times longer and provide additional cytological markers. However, the application of this technique to pachytene chromosomes has been complicated by problems in making preparations of meiotic chromosomes and by difficulties in the application of standard FISH protocols. We have been able to overcome most of these obstacles in applying the FISH technique to the pachytene chromosomes of coffee plants. Digesting the external callose layer surrounding the pollen mother cells (PMCs) in conjunction with other procedures permitted suitable pachytene chromosomes to be obtained by increasing cell permeability, which allowed the probe sequences to enter the cells. For the first time, hybridization signals were registered on coffee pachytene chromosomes using the FISH technique with a repetitive sequence as a probe. We obtained slides on which 80 % of the PMCs had hybridization signals, resulting in FISH labelling with high efficiency. The procedure does not seem to be dependent on the genotype, because hybridization signals were detected in genetically different coffee plants. These findings enhance the possibilities for high-resolution physical mapping of coffee chromosomes.

Isolation and characterization of reverse transcriptase fragments of LTR retrotransposons from the genome of Chenopodium quinoa (Amaranthaceae)

High heterogeneity was observed among conserved domains of reverse transcriptase ( rt ) isolated from quinoa. Only one Ty1- copia rt was highly amplified. Reverse transcriptase sequences were located predominantly in pericentromeric region of quinoa chromosomes. The heterogeneity, genomic abundance, and chromosomal distribution of reverse transcriptase (rt)-coding fragments of Ty1-copia and Ty3-gypsy long terminal repeat retrotransposons were analyzed in the Chenopodium quinoa genome. Conserved domains of the rt gene were amplified and characterized using degenerate oligonucleotide primer pairs. Sequence analyses indicated that half of Ty1-copia rt (51 %) and 39 % of Ty3-gypsy rt fragments contained intact reading frames. High heterogeneity among rt sequences was observed for both Ty1-copia and Ty3-gypsy rt amplicons, with Ty1-copia more heterogeneous than Ty3-gypsy. Most of the isolated rt fragments were present in quinoa genome in low copy numbers, with only one highly amplified Ty1-copia rt sequence family. The gypsy-like RNase H fragments co-amplified with Ty1-copia-degenerate primers were shown to be highly amplified in the quinoa genome indicating either higher abundance of some gypsy families of which rt domains could not be amplified, or independent evolution of this gypsy-region in quinoa. Both Ty1-copia and Ty3-gypsy retrotransposons were preferentially located in pericentromeric heterochromatin of quinoa chromosomes. Phylogenetic analyses of newly amplified rt fragments together with well-characterized retrotransposon families from other organisms allowed identification of major lineages of retroelements in the genome of quinoa and provided preliminary insight into their evolutionary dynamics.

TARE1, a mutated Copia-like LTR retrotransposon followed by recent massive amplification in tomato

Long terminal repeat retrotransposons (LTR-RTs) are the major DNA components in flowering plants. Most LTR-RTs contain dinucleotides ‘TG’ and ‘CA’ at the ends of the two LTRs. Here we report the structure, evolution, and propensity of a tomato atypical retrotransposon element (TARE1) with both LTRs starting as ‘TA’. This family is also characterized by high copy numbers (354 copies), short LTR size (194 bp), extremely low ratio of solo LTRs to intact elements (0.05∶1), recent insertion (most within 0.75∼1.75 million years, Mys), and enrichment in pericentromeric region. The majority (83%) of the TARE1 elements are shared between S. lycopersicum and its wild relative S. pimpinellifolium, but none of them are found in potato. In the present study, we used shared LTR-RTs as molecular markers and estimated the divergence time between S. lycopersicum and S. pimpinellifolium to be <0.5 Mys. Phylogenetic analysis showed that the TARE1 elements, together with two closely related families, TARE2 and TGRE1, have formed a sub-lineage belonging to a Copia-like Ale lineage. Although TARE1 and TARE2 shared similar structural characteristics, the timing, scale, and activity of their amplification were found to be substantially different. We further propose a model wherein a single mutation from ‘G’ to ‘A’ in 3′ LTR followed by amplification is responsible for the origin of TARE1, thus providing evidence that the proliferation of a spontaneous mutation can be mediated by the amplification of LTR-RTs at the level of RNA.

First insights into the large genome of Epimedium sagittatum (Sieb. et Zucc) Maxim, a Chinese Ttaditional medicinal plant

Epimedium sagittatum (Sieb. et Zucc) Maxim is a member of the Berberidaceae family of basal eudicot plants, widely distributed and used as a traditional medicinal plant in China for therapeutic effects on many diseases with a long history. Recent data shows that E. sagittatum has a relatively large genome, with a haploid genome size of ~4496 Mbp, divided into a small number of only 12 diploid chromosomes (2n = 2x = 12). However, little is known about Epimedium genome structure and composition. Here we present the analysis of 691 kb of high-quality genomic sequence derived from 672 randomly selected plasmid clones of E. sagittatum genomic DNA, representing ~0.0154% of the genome. The sampled sequences comprised at least 78.41% repetitive DNA elements and 2.51% confirmed annotated gene sequences, with a total GC% content of 39%. Retrotransposons represented the major class of transposable element (TE) repeats identified (65.37% of all TE repeats), particularly LTR (Long Terminal Repeat) retrotransposons (52.27% of all TE repeats). Chromosome analysis and Fluorescence in situ Hybridization of Gypsy-Ty3 retrotransposons were performed to survey the E. sagittatum genome at the cytological level. Our data provide the first insights into the composition and structure of the E. sagittatum genome, and will facilitate the functional genomic analysis of this valuable medicinal plant.

Identification and QTL mapping of whitefly resistance components in Solanum galapagense

Solanum galapagense is closely related to the cultivated tomato and can show a very good resistance towards whitefly. A segregating population resulting from a cross between the cultivated tomato and a whitefly resistant S. galapagense was created and used for mapping whitefly resistance and related traits, which made it possible to study the genetic basis of the resistance. Quantitative trait loci (QTL) for adult survival co-localized with type IV trichome characteristics (presence, density, gland longevity and gland size). A major QTL (Wf-1) was found for adult survival and trichome characters on Chromosome 2. This QTL explained 54.1 % of the variation in adult survival and 81.5 % of the occurrence of type IV trichomes. A minor QTL (Wf-2) for adult survival and trichome characters was identified on Chromosome 9. The major QTL was confirmed in F3 populations. Comprehensive metabolomics, based on GCMS profiling, revealed that 16 metabolites segregating in the F2 mapping population were associated with Wf-1 and/or Wf-2. Analysis of the 10 most resistant and susceptible F2 genotypes by LCMS showed that several acyl sugars were present in significantly higher concentration in the whitefly resistant genotypes, suggesting a role for these components in the resistance as well. Our results show that whitefly resistance in S. galapagense seems to inherit relatively simple compared to whitefly resistance from other sources and this offers great prospects for resistance breeding as well as elucidating the underlying molecular mechanism(s) of the resistance.

Mapping of Micro-Tom BAC-End Sequences to the Reference Tomato Genome Reveals Possible Genome Rearrangements and Polymorphisms

A total of 93,682 BAC-end sequences (BESs) were generated from a dwarf model tomato, cv. Micro-Tom. After removing repetitive sequences, the BESs were similarity searched against the reference tomato genome of a standard cultivar, "Heinz 1706." By referring to the "Heinz 1706" physical map and by eliminating redundant or nonsignificant hits, 28,804 "unique pair ends" and 8,263 "unique ends" were selected to construct hypothetical BAC contigs. The total physical length of the BAC contigs was 495, 833, 423 bp, covering 65.3% of the entire genome. The average coverage of euchromatin and heterochromatin was 58.9% and 67.3%, respectively. From this analysis, two possible genome rearrangements were identified: one in chromosome 2 (inversion) and the other in chromosome 3 (inversion and translocation). Polymorphisms (SNPs and Indels) between the two cultivars were identified from the BLAST alignments. As a result, 171,792 polymorphisms were mapped on 12 chromosomes. Among these, 30,930 polymorphisms were found in euchromatin (1 per 3,565 bp) and 140,862 were found in heterochromatin (1 per 2,737 bp). The average polymorphism density in the genome was 1 polymorphism per 2,886 bp. To facilitate the use of these data in Micro-Tom research, the BAC contig and polymorphism information are available in the TOMATOMICS database.

Evolution of the large genome in Capsicum annuum occurred through accumulation of single-type long terminal repeat retrotransposons and their derivatives

Although plant genome sizes are extremely diverse, the mechanism underlying the expansion of huge genomes that did not experience whole-genome duplication has not been elucidated. The pepper, Capsicum annuum, is an excellent model for studies of genome expansion due to its large genome size (2700 Mb) and the absence of whole genome duplication. As most of the pepper genome structure has been identified as constitutive heterochromatin, we investigated the evolution of this region in detail. Our findings show that the constitutive heterochromatin in pepper was actively expanded 20.0-7.5 million years ago through a massive accumulation of single-type Ty3/Gypsy-like elements that belong to the Del subgroup. Interestingly, derivatives of the Del elements, such as non-autonomous long terminal repeat retrotransposons and long-unit tandem repeats, played important roles in the expansion of constitutive heterochromatic regions. This expansion occurred not only in the existing heterochromatic regions but also into the euchromatic regions. Furthermore, our results revealed a repeat of unit length 18-24 kb. This repeat was found not only in the pepper genome but also in the other solanaceous species, such as potato and tomato. These results represent a characteristic mechanism for large genome evolution in plants.

Euchromatic and heterochromatic compositional properties emerging from the analysis of Solanum lycopersicum BAC sequences

The consortium responsible for the sequencing of the tomato (Solanum lycopersicum) genome initially focused on the sequencing of the euchromatic regions using a BAC-by-BAC strategy. We analyzed the compositional features of the whole collection of BAC sequences publically available. This analysis highlights specific peculiarities of heterochromatic and euchromatic BACs, in particular: the whole BAC collection has i) a large variability in repeat and gene content, ii) a positive and significant correlation of LTR retrotransposons of the Gypsy class with the repeat content and iii) the preferential location of the SINEs (short interspersed nuclear elements) in BAC sequences showing a low repeat content. Our results point out a typical design of the tomato chromosomes and pave the way for further investigations on the relationship between DNA primary structure and chromatin organization in Solanaceae genomes.

Minimizing the unpredictability of transgene expression in plants: the role of genetic insulators

The genetic transformation of plants has become a necessary tool for fundamental plant biology research, as well as the generation of engineered plants exhibiting improved agronomic and industrial traits. However, this technology is significantly hindered by the fact that transgene expression is often highly variable amongst independent transgenic lines. Two of the major contributing factors to this type of inconsistency are inappropriate enhancer-promoter interactions and chromosomal position effects, which frequently result in mis-expression or silencing of the transgene, respectively. Since the precise, often tissue-specific, expression of the transgene(s) of interest is often a necessity for the successful generation of transgenic plants, these undesirable side effects have the potential to pose a major challenge for the genetic engineering of these organisms. In this review, we discuss strategies for improving foreign gene expression in plants via the inclusion of enhancer-blocking insulators, which function to impede enhancer-promoter communication, and barrier insulators, which block the spread of heterochromatin, in transgenic constructs. While a complete understanding of these elements remains elusive, recent studies regarding their use in genetically engineered plants indicate that they hold great promise for the improvement of transgene expression, and thus the future of plant biotechnology.

Comparative analysis of pepper and tomato reveals euchromatin expansion of pepper genome caused by differential accumulation of Ty3/Gypsy-like elements

Background

Among the Solanaceae plants, the pepper genome is three times larger than that of tomato. Although the gene repertoire and gene order of both species are well conserved, the cause of the genome-size difference is not known. To determine the causes for the expansion of pepper euchromatic regions, we compared the pepper genome to that of tomato.

Results

For sequence-level analysis, we generated 35.6 Mb of pepper genomic sequences from euchromatin enriched 1,245 pepper BAC clones. The comparative analysis of orthologous gene-rich regions between both species revealed insertion of transposons exclusively in the pepper sequences, maintaining the gene order and content. The most common type of the transposon found was the LTR retrotransposon. Phylogenetic comparison of the LTR retrotransposons revealed that two groups of Ty3/Gypsy-like elements (Tat and Athila) were overly accumulated in the pepper genome. The FISH analysis of the pepper Tat elements showed a random distribution in heterochromatic and euchromatic regions, whereas the tomato Tat elements showed heterochromatin-preferential accumulation.

Conclusions

Compared to tomato pepper euchromatin doubled its size by differential accumulation of a specific group of Ty3/Gypsy-like elements. Our results could provide an insight on the mechanism of genome evolution in the Solanaceae family.

SNP discovery and linkage map construction in cultivated tomato

Few intraspecific genetic linkage maps have been reported for cultivated tomato, mainly because genetic diversity within Solanum lycopersicum is much less than that between tomato species. Single nucleotide polymorphisms (SNPs), the most abundant source of genomic variation, are the most promising source of polymorphisms for the construction of linkage maps for closely related intraspecific lines. In this study, we developed SNP markers based on expressed sequence tags for the construction of intraspecific linkage maps in tomato. Out of the 5607 SNP positions detected through in silico analysis, 1536 were selected for high-throughput genotyping of two mapping populations derived from crosses between 'Micro-Tom' and either 'Ailsa Craig' or 'M82'. A total of 1137 markers, including 793 out of the 1338 successfully genotyped SNPs, along with 344 simple sequence repeat and intronic polymorphism markers, were mapped onto two linkage maps, which covered 1467.8 and 1422.7 cM, respectively. The SNP markers developed were then screened against cultivated tomato lines in order to estimate the transferability of these SNPs to other breeding materials. The molecular markers and linkage maps represent a milestone in the genomics and genetics, and are the first step toward molecular breeding of cultivated tomato. Information on the DNA markers, linkage maps, and SNP genotypes for these tomato lines is available at http://www.kazusa.or.jp/tomato/.

An interspecific linkage map of SSR and intronic polymorphism markers in tomato

Despite the collection and availability of abundant tomato genome sequences, PCR-based markers adapted to large scale analysis have not been developed in tomato species. Therefore, using public genome sequence data in tomato, we developed three types of DNA markers: expressed sequence tag (EST)-derived simple sequence repeat (SSR) markers (TES markers), genome-derived SSR markers (TGS markers) and EST-derived intronic polymorphism markers (TEI markers). A total of 2,047 TES, 3,510 TGS and 674 TEI markers were established and used in the polymorphic analysis of a cultivated tomato (Solanum lycopersicum) 'LA925' and its wild relative Solanum pennellii 'LA716', parents of the Tomato-EXPEN 2000 mapping population. The polymorphic ratios between parents revealed by the TES, TGS and TEI markers were 37.3, 22.6 and 80.0%, respectively. Those showing polymorphisms were used to genotype the Tomato-EXPEN 2000 mapping population, and a high-density genetic linkage map composed of 1,433 new and 683 existing marker loci was constructed on 12 chromosomes, covering 1,503.1 cM. In the present map, 48% of the mapped TGS loci were located within heterochromatic regions, while 18 and 21% of TES and TEI loci, respectively, were located in heterochromatin. The large number of SSR and SNP markers developed in this study provide easily handling genomic tools for molecular breeding in tomato. Information on the DNA markers developed in this study is available at http://www.kazusa.or.jp/tomato/.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

Review of the Application of Modern Cytogenetic Methods (FISH/GISH) to the Study of Reticulation (Polyploidy/Hybridisation)

The convergence of distinct lineages upon interspecific hybridisation, including when accompanied by increases in ploidy (allopolyploidy), is a driving force in the origin of many plant species. In plant breeding too, both interspecific hybridisation and allopolyploidy are important because they facilitate introgression of alien DNA into breeding lines enabling the introduction of novel characters. Here we review how fluorescence in situ hybridisation (FISH) and genomic in situ hybridisation (GISH) have been applied to: 1) studies of interspecific hybridisation and polyploidy in nature, 2) analyses of phylogenetic relationships between species, 3) genetic mapping and 4) analysis of plant breeding materials. We also review how FISH is poised to take advantage of nextgeneration sequencing (NGS) technologies, helping the rapid characterisation of the repetitive fractions of a genome in natural populations and agricultural plants.

ORTom: a multi-species approach based on conserved co-expression to identify putative functional relationships among genes in tomato

Co-expressed genes are often expected to be functionally related and many bioinformatics approaches based on co-expression have been developed to infer their biological role. However, such annotations may be unreliable, whereas the evolutionary conservation of gene co-expression among species may form a basis for more confident predictions. The huge amount of expression data (microarrays, SAGE, ESTs) has already allowed functional studies based on conserved co-expression in animals. Up to now, the implementation of analogous tools for plants has been strongly limited probably by the paucity and heterogeneity of data. Here we present ORTom, a tomato-centred EST data-mining approach based on conserved co-expression in the Solanaceae family. ORTom can be used to predict functional relationships among genes and to prioritize candidate genes for targeted studies. The method consists in ranking ESTs co-expressed with a gene of interest according to the level of expression pattern conservation in phylogenetically-related plants (potato, tobacco and pepper) to obtain lists of putative functionally-related genes. The lists are then analyzed for Gene Ontology keyword enrichment. The web server ORTom has been implemented to make the results publicly-available and searchable. Few biological examples on how the tool can be used are presented.

Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4-1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.

A lineage-specific centromere retrotransposon in Oryza brachyantha

Most eukaryotic centromeres contain large quantities of repetitive DNA, such as satellite repeats and retrotransposons. Unlike most transposons in plant genomes, the centromeric retrotransposon (CR) family is conserved over long evolutionary periods among a majority of the grass species. CR elements are highly concentrated in centromeres, and are likely to play a role in centromere function. In order to study centromere evolution in the Oryza (rice) genus, we sequenced the orthologous region to centromere 8 of Oryza sativa from a related species, Oryza brachyantha. We found that O. brachyantha does not have the canonical CRR (CR of rice) found in the centromeres of all other Oryza species. Instead, a new Ty3-gypsy (Metaviridae) retroelement (FRetro3) was found to colonize the centromeres of this species. This retroelement is found in high copy numbers in the O. brachyantha genome, but not in other Oryza genomes, and based on the dating of long terminal repeats (LTRs) of FRetro3 it was amplified in the genome in the last few million years. Interestingly, there is a high level of removal of FRetro3 based on solo-LTRs to full-length elements, and this rapid turnover may have played a role in the replacement of the canonical CRR with the new element by active deletion. Comparison with previously described ChIP cloning data revealed that FRetro3 is found in CENH3-associated chromatin sequences. Thus, within a single lineage of the Oryza genus, the canonical component of grass centromeres has been replaced with a new retrotransposon that has all the hallmarks of a centromeric retroelement.

Genome organization of the tomato sun locus and characterization of the unusual retrotransposon Rider

DNA sequences provide useful insights into genome structure and organization as well as evolution of species. We report on a detailed analysis of the locus surrounding the tomato (Solanum lycopersicum) fruit-shape gene SUN to determine the driving force and genome environment that foster the appearance of novel phenotypes. The gene density at the sun locus is similar to that described in other euchromatic portions of the tomato genome despite the relatively high number of transposable elements. Genes at the sun locus include protein-coding as well as RNA genes, are small in size, and belong to families that were duplicated at the locus an estimated 5-74 million years ago. In general, the DNA transposons at the sun locus are older than the RNA transposons, and their insertion pre-dates the speciation of S. lycopersicum and S. pimpinellifolium. Gene redundancy and large intergenic regions may explain the tolerance of the sun locus to frequent rearrangements and transpositions. The most recent transposition event at the sun locus involved Rider, a recently discovered high-copy retrotransposon. Rider probably arose early during the speciation of tomato. The element inserts into or near to genes and may still be active, which are unusual features for a high-copy element. Rider full-length and read-through transcripts past the typical transcription termination stop are detected, and the latter are required for mobilizing nearby sequences. Rider activity has resulted in an altered phenotype in three known cases, and may therefore have played an important role in tomato evolution and domestication.

Assignment of genetic linkage maps to diploid Solanum tuberosum pachytene chromosomes by BAC-FISH technology

A cytogenetic map has been developed for diploid potato (Solanum tuberosum), in which the arms of the 12 potato bivalents can be identified in pachytene complements using multicolor fluorescence in situ hybridization (FISH) with a set of 60 genetically anchored bacterial artificial chromosome (BAC) clones from the RHPOTKEY BAC library. This diagnostic set of selected BACs (five per chromosome) hybridizes to euchromatic regions and corresponds to well-defined loci in the ultradense genetic map, and with these probes a new detailed and reliable pachytene karyotype could be established. Chromosome size has been estimated both from microscopic length measurements and from 4′,6-diamidino-2-phenylindole fluorescence-based DNA content measurements. In both approaches, chromosome 1 is the largest (100–115 Mb) and chromosome 11 the smallest (49–53 Mb). Detailed measurements of mega-base-pair to micrometer ratios have been obtained for chromosome 5, with average values of 1.07 Mb/μm for euchromatin and 3.67 Mb/μm for heterochromatin. In addition, our FISH results helped to solve two discrepancies in the potato genetic map related to chromosomes 8 and 12. Finally, we discuss the significance of the potato cytogenetic map for extended FISH studies in potato and related Solanaceae, which will be especially beneficial for the potato genome-sequencing project.

Genome analysis and genetic enhancement of tomato

The Solanaceae is an important family of vegetable crops, ornamentals and medicinal plants. Tomato has served as a model member of this family largely because of its enriched cytogenetic, genetic, as well as physical, maps. Mapping has helped in cloning several genes of importance such as Pto, responsible for resistance against bacterial speck disease, Mi-1.2 for resistance against nematodes, and fw2.2 QTL for fruit weight. A high-throughput genome-sequencing program has been initiated by an international consortium of 10 countries. Since heterochromatin has been found to be concentrated near centromeres, the consortium is focusing on sequencing only the gene-rich euchromatic region. Genomes of the members of Solanaceae show a significant degree of synteny, suggesting that the tomato genome sequence would help in the cloning of genes for important traits from other Solanaceae members as well. ESTs from a large number of cDNA libraries have been sequenced, and microarray chips, in conjunction with wide array of ripening mutants, have contributed immensely to the understanding of the fruit-ripening phenomenon. Work on the analysis of the tomato proteome has also been initiated. Transgenic tomato plants with improved abiotic stress tolerance, disease resistance and insect resistance, have been developed. Attempts have also been made to develop tomato as a bioreactor for various pharmaceutical proteins. However, control of fruit quality and ripening remains an active and challenging area of research. Such efforts should pave the way to improve not only tomato, but also other solanaceous crops.

The genomic organization of Ty3/gypsy-like retrotransposons in Helianthus (Asteraceae) homoploid hybrid species

The origin of new diploid, or homoploid, hybrid species is associated with rapid genomic restructuring in the hybrid neospecies. This mode of speciation has been best characterized in wild sunflower species in the genus Helianthus, where three homoploid hybrid species (H. anomalus, H. deserticola, and H. paradoxus) have independently arisen via ancient hybridization events between the same two parental species (H. annuus and H. petiolaris). Most previous work examining genomic restructuring in these sunflower hybrid species has focused on chromosomal rearrangements. However, the origin of all three homoploid hybrid sunflower species also is associated with massive proliferation events of Ty3/gypsy-like retrotransposons in the hybrid species' genomes. We compared the genomic organization of these elements in the parent species and two of the homoploid hybrid species using fluorescence in situ hybridization (FISH). We found a significant expansion of Ty3/gypsy-like retrotransposons confined to the pericentromeric regions of two hybrid sunflower species, H. deserticola and H. paradoxus. In contrast, we detected no significant increase in the frequency or extent of dispersed retrotransposon populations in the hybrid species within the resolution limits of our assay. We discuss the potential role that transposable element proliferation and localization plays in the evolution of homoploid hybrid species.

A physical map of the papaya genome with integrated genetic map and genome sequence

BACKGROUND

Papaya is a major fruit crop in tropical and subtropical regions worldwide and has primitive sex chromosomes controlling sex determination in this trioecious species. The papaya genome was recently sequenced because of its agricultural importance, unique biological features, and successful application of transgenic papaya for resistance to papaya ringspot virus. As a part of the genome sequencing project, we constructed a BAC-based physical map using a high information-content fingerprinting approach to assist whole genome shotgun sequence assembly.

RESULTS

The physical map consists of 963 contigs, representing 9.4x genome equivalents, and was integrated with the genetic map and genome sequence using BAC end sequences and a sequence-tagged high-density genetic map. The estimated genome coverage of the physical map is about 95.8%, while 72.4% of the genome was aligned to the genetic map. A total of 1,181 high quality overgo (overlapping oligonucleotide) probes representing conserved sequences in Arabidopsis and genetically mapped loci in Brassica were anchored on the physical map, which provides a foundation for comparative genomics in the Brassicales. The integrated genetic and physical map aligned with the genome sequence revealed recombination hotspots as well as regions suppressed for recombination across the genome, particularly on the recently evolved sex chromosomes. Suppression of recombination spread to the adjacent region of the male specific region of the Y chromosome (MSY), and recombination rates were recovered gradually and then exceeded the genome average. Recombination hotspots were observed at about 10 Mb away on both sides of the MSY, showing 7-fold increase compared with the genome wide average, demonstrating the dynamics of recombination of the sex chromosomes.

CONCLUSION

A BAC-based physical map of papaya was constructed and integrated with the genetic map and genome sequence. The integrated map facilitated the draft genome assembly, and is a valuable resource for comparative genomics and map-based cloning of agronomically and economically important genes and for sex chromosome research.

Different chromatin fractions of tomato (Solanum lycopersicum L.) and related species

Conventional chromosome staining has suggested that more than 75% of the tomato chromosomes are constituted by heterochromatin. In order to determine whether more deeply stained proximal regions are classic heterochromatin, the distributions of C-bands and chromomycin A(3) (CMA) bands, and the prophase condensation patterns, were analysed in tomato. In this and most other species of the tomato clade, the 5S and 45S rDNA sites were also localised. In tomato, CMA banding was similar to C-banding. After conventional staining, all species displayed large condensed heteropycnotic regions that did not correspond to C-bands or CMA bands. Analyses of the CMA banded karyotypes revealed a low heterochromatin content. Around 12-17% of the chromatin of tomato was CMA(+) and 1/4 to 1/5 of this heterochromatin corresponded to 45S rDNA. In other species, the CMA(+) heterochromatin showed extensive variation (8-35%), but was never near the values found in the literature for tomato. These data suggest the existence of three principal fractions of chromatin in tomato and related species: the late condensed euchromatin corresponding to the terminal regions of the chromosomes, the precocious condensed euchromatin that occupies the major part of the chromosomes and the constitutive heterochromatin that represents those regions revealed by C-bands.

Solanum lycopersicum cv. Heinz 1706 chromosome 6: distribution and abundance of genes and retrotransposable elements

We studied the physical and genetic organization of chromosome 6 of tomato (Solanum lycopersicum) cv. Heinz 1706 by combining bacterial artificial chromosome (BAC) sequence analysis, high-information-content fingerprinting, genetic analysis, and BAC-fluorescent in situ hybridization (FISH) mapping data. The chromosome positions of 81 anchored seed and extension BACs corresponded in most cases with the linear marker order on the high-density EXPEN 2000 linkage map. We assembled 25 BAC contigs and eight singleton BACs spanning 2.0 Mb of the short-arm euchromatin, 1.8 Mb of the pericentromeric heterochromatin and 6.9 Mb of the long-arm euchromatin. Sequence data were combined with their corresponding genetic and pachytene chromosome positions into an integrated map that covers approximately a third of the chromosome 6 euchromatin and a small part of the pericentromeric heterochromatin. We then compared physical length (Mb), genetic (cM) and chromosome distances (microm) for determining gap sizes between contigs, revealing relative hot and cold spots of recombination. Through sequence annotation we identified several clusters of functionally related genes and an uneven distribution of both gene and repeat sequences between heterochromatin and euchromatin domains. Although a greater number of the non-transposon genes were located in the euchromatin, the highly repetitive (22.4%) pericentromeric heterochromatin displayed an unexpectedly high gene content of one gene per 36.7 kb. Surprisingly, the short-arm euchromatin was relatively rich in repeats as well, with a repeat content of 13.4%, yet the ratio of Ty3/Gypsy and Ty1/Copia retrotransposable elements across the chromosome clearly distinguished euchromatin (2:3) from heterochromatin (3:2).