Analysis of 90 Mb of the potato genome reveals conservation of gene structures and order with tomato but divergence in repetitive sequence composition

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Characterization of two group III potato CDPKs, StCDPK22 and StCDPK24, that contain three EF-Hand motifs in their CLDs

Four members of the potato (Solanum tuberosum L.) calcium-dependent protein kinase (CDPK) family StCDPK22/23/24 and StCDPK27, present three functional EF-hands motifs in their calmodulin-like domain (CLD). StCDPK22/23/24 are clustered in clade III-b1 with tomato and Arabidopsis CDPKs that lack the first EF-hand motif, while StCDPK27 is clustered in clade III-b3 with CDPKs that lack EF-hand 2. Members of each clade share similar intron-exon structures and acylation profiles. 3D model predictions suggested that StCDPK22 and StCDPK24 are active kinases that undergo a conformational switch in the presence of Ca²⁺ even when lacking one functional EF-hand motif; however, assays performed with recombinant proteins indicated that StCDPK24:6xHis was active in all the conditions tested, and its activity was enhanced in the presence of Ca²⁺, but StCDPK22:6xHis had scarce or null activity. Both kinases share with AtCPK8 the same autophosphorylation pattern in the autoinhibitory (AD) and C-terminal variable (CTV) domains, suggesting that it could be a characteristic of clade III-b1. RT-qPCR analysis revealed that StCDPK22 is mainly expressed in early stages of tuberization, but not limited to, while StCDPK24 expression is more ubiquitous. In silico analysis predicted several abiotic stress-responsive elements in its promoters. Accordingly, StCDPK24 expression peaked at 10 h in in vitro plants exposed to salt shock and then declined. Moreover, a significant increase was observed at 2 h in stems of salt-treated greenhouse plants, suggesting that this CDPK could participate in the early events of the signaling cascade triggered in response to salt.

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.

Development of Genomic Microsatellite Markers in Carthamus tinctorius L. (Safflower) Using Next Generation Sequencing and Assessment of Their Cross-Species Transferability and Utility for Diversity Analysis

BACKGROUND

Safflower (Carthamus tinctorius L.), an Asteraceae member, yields high quality edible oil rich in unsaturated fatty acids and is resilient to dry conditions. The crop holds tremendous potential for improvement through concerted molecular breeding programs due to the availability of significant genetic and phenotypic diversity. Genomic resources that could facilitate such breeding programs remain largely underdeveloped in the crop. The present study was initiated to develop a large set of novel microsatellite markers for safflower using next generation sequencing.

PRINCIPAL FINDINGS

Low throughput genome sequencing of safflower was performed using Illumina paired end technology providing ~3.5X coverage of the genome. Analysis of sequencing data allowed identification of 23,067 regions harboring perfect microsatellite loci. The safflower genome was found to be rich in dinucleotide repeats followed by tri-, tetra-, penta- and hexa-nucleotides. Primer pairs were designed for 5,716 novel microsatellite sequences with repeat length ≥ 20 bases and optimal flanking regions. A subset of 325 microsatellite loci was tested for amplification, of which 294 loci produced robust amplification. The validated primers were used for assessment of 23 safflower accessions belonging to diverse agro-climatic zones of the world leading to identification of 93 polymorphic primers (31.6%). The numbers of observed alleles at each locus ranged from two to four and mean polymorphism information content was found to be 0.3075. The polymorphic primers were tested for cross-species transferability on nine wild relatives of cultivated safflower. All primers except one showed amplification in at least two wild species while 25 primers amplified across all the nine species. The UPGMA dendrogram clustered C. tinctorius accessions and wild species separately into two major groups. The proposed progenitor species of safflower, C. oxyacantha and C. palaestinus were genetically closer to cultivated safflower and formed a distinct cluster. The cluster analysis also distinguished diploid and tetraploid wild species of safflower.

CONCLUSION

Next generation sequencing of safflower genome generated a large set of microsatellite markers. The novel markers developed in this study will add to the existing repertoire of markers and can be used for diversity analysis, synteny studies, construction of linkage maps and marker-assisted selection.

A Deluge of Complex Repeats: The Solanum Genome

Repetitive elements have lately emerged as key components of genome, performing varieties of roles. It has now become necessary to have an account of repeats for every genome to understand its dynamics and state. Recently, genomes of two major Solanaceae species, Solanum tuberosum and Solanum lycopersicum, were sequenced. These species are important crops having high commercial significance as well as value as model species. However, there is a reasonable gap in information about repetitive elements and their possible roles in genome regulation for these species. The present study was aimed at detailed identification and characterization of complex repetitive elements in these genomes, along with study of their possible functional associations as well as to assess possible transcriptionally active repetitive elements. In this study, it was found that ~50-60% of genomes of S. tuberosum and S. lycopersicum were composed of repetitive elements. It was also found that complex repetitive elements were associated with >95% of genes in both species. These two genomes are mostly composed of LTR retrotransposons. Two novel repeat families very similar to LTR/ERV1 and LINE/RTE-BovB have been reported for the first time. Active existence of complex repeats was estimated by measuring their transcriptional abundance using Next Generation Sequencing read data and Microarray platforms. A reasonable amount of regulatory components like transcription factor binding sites and miRNAs appear to be under the influence of these complex repetitive elements in these species, while several genes appeared to possess exonized repeats.

Homologues of potato chromosome 5 show variable collinearity in the euchromatin, but dramatic absence of sequence similarity in the pericentromeric heterochromatin

BACKGROUND

In flowering plants it has been shown that de novo genome assemblies of different species and genera show a significant drop in the proportion of alignable sequence. Within a plant species, however, it is assumed that different haplotypes of the same chromosome align well. In this paper we have compared three de novo assemblies of potato chromosome 5 and report on the sequence variation and the proportion of sequence that can be aligned.

RESULTS

For the diploid potato clone RH89-039-16 (RH) we produced two linkage phase controlled and haplotype-specific assemblies of chromosome 5 based on BAC-by-BAC sequencing, which were aligned to each other and compared to the 52 Mb chromosome 5 reference sequence of the doubled monoploid clone DM 1-3 516 R44 (DM). We identified 17.0 Mb of non-redundant sequence scaffolds derived from euchromatic regions of RH and 38.4 Mb from the pericentromeric heterochromatin. For 32.7 Mb of the RH sequences the correct position and order on chromosome 5 was determined, using genetic markers, fluorescence in situ hybridisation and alignment to the DM reference genome. This ordered fraction of the RH sequences is situated in the euchromatic arms and in the heterochromatin borders. In the euchromatic regions, the sequence collinearity between the three chromosomal homologs is good, but interruption of collinearity occurs at nine gene clusters. Towards and into the heterochromatin borders, absence of collinearity due to structural variation was more extensive and was caused by hemizygous and poorly aligning regions of up to 450 kb in length. In the most central heterochromatin, a total of 22.7 Mb sequence from both RH haplotypes remained unordered. These RH sequences have very few syntenic regions and represent a non-alignable region between the RH and DM heterochromatin haplotypes of chromosome 5.

CONCLUSIONS

Our results show that among homologous potato chromosomes large regions are present with dramatic loss of sequence collinearity. This stresses the need for more de novo reference assemblies in order to capture genome diversity in this crop. The discovery of three highly diverged pericentric heterochromatin haplotypes within one species is a novelty in plant genome analysis. The possible origin and cytogenetic implication of this heterochromatin haplotype diversity are discussed.

Transcriptome responses to Ralstonia solanacearum infection in the roots of the wild potato Solanum commersonii

BACKGROUND

Solanum commersonii is a wild potato species that exhibits high tolerance to both biotic and abiotic stresses and has been used as a source of genes for introgression into cultivated potato. Among the interesting features of S. commersonii is resistance to the bacterial wilt caused by Ralstonia solanacearum, one of the most devastating bacterial diseases of crops.

RESULTS

In this study, we used deep sequencing of S. commersonii RNA (RNA-seq) to analyze the below-ground plant transcriptional responses to R. solanacearum. While a majority of S. commersonii RNA-seq reads could be aligned to the Solanum tuberosum Group Phureja DM reference genome sequence, we identified 2,978 S. commersonii novel transcripts through assembly of unaligned S. commersonii RNA-seq reads. We also used RNA-seq to study gene expression in pathogen-challenged roots of S. commersonii accessions resistant (F118) and susceptible (F97) to the pathogen. Expression profiles obtained from read mapping to the S. tuberosum reference genome and the S. commersonii novel transcripts revealed a differential response to the pathogen in the two accessions, with 221 (F118) and 644 (F97) differentially expressed genes including S. commersonii novel transcripts in the resistant and susceptible genotypes. Interestingly, 22.6% of the F118 and 12.8% of the F97 differentially expressed genes had been previously identified as responsive to biotic stresses and half of those up-regulated in both accessions had been involved in plant pathogen responses. Finally, we compared two different methods to eliminate ribosomal RNA from the plant RNA samples in order to allow dual mapping of RNAseq reads to the host and pathogen genomes and provide insights on the advantages and limitations of each technique.

CONCLUSIONS

Our work catalogues the S. commersonii transcriptome and strengthens the notion that this species encodes specific genes that are differentially expressed to respond to bacterial wilt. In addition, a high proportion of S. commersonii-specific transcripts were altered by R. solanacearum only in F118 accession, while phythormone-related genes were highly induced in F97, suggesting a markedly different response to the pathogen in the two plant accessions studied.

Chromosomal organizations of major repeat families on potato (Solanum tuberosum) and further exploring in its sequenced genome

One of the most powerful technologies in unraveling the organization of a eukaryotic plant genome is high-resolution Fluorescent in situ hybridization of repeats and single copy DNA sequences on pachytene chromosomes. This technology allows the integration of physical mapping information with chromosomal positions, including centromeres, telomeres, nucleolar-organizing region, and euchromatin and heterochromatin. In this report, we established chromosomal positions of different repeat fractions of the potato genomic DNA (Cot100, Cot500 and Cot1000) on the chromosomes. We also analysed various repeat elements that are unique to potato including the moderately repetitive P5 and REP2 elements, where the REP2 is part of a larger Gypsy-type LTR retrotransposon and cover most chromosome regions, with some brighter fluorescing spots in the heterochromatin. The most abundant tandem repeat is the potato genomic repeat 1 that covers subtelomeric regions of most chromosome arms. Extensive multiple alignments of these repetitive sequences in the assembled RH89-039-16 potato BACs and the draft assembly of the DM1-3 516 R44 genome shed light on the conservation of these repeats within the potato genome. The consensus sequences thus obtained revealed the native complete transposable elements from which they were derived.

Genome-wide data (ChIP-seq) enabled identification of cell wall-related and aquaporin genes as targets of tomato ASR1, a drought stress-responsive transcription factor

BACKGROUND

Identifying the target genes of transcription factors is important for unraveling regulatory networks in all types of organisms. Our interest was precisely to uncover the spectrum of loci regulated by a widespread plant transcription factor involved in physiological adaptation to drought, a type of stress that plants have encountered since the colonization of land habitats 400 MYA. The regulator under study, named ASR1, is exclusive to the plant kingdom (albeit absent in Arabidopsis) and known to alleviate the stress caused by restricted water availability. As its target genes are still unknown despite the original cloning of Asr1 cDNA 20 years ago, we examined the tomato genome for specific loci interacting in vivo with this conspicuous protein.

RESULTS

We performed ChIP followed by high throughput DNA sequencing (ChIP-seq) on leaves from stressed tomato plants, using a high-quality anti-ASR1 antibody. In this way, we unraveled a novel repertoire of target genes, some of which are clearly involved in the response to drought stress. Many of the ASR1-enriched genomic loci we found encode enzymes involved in cell wall synthesis and remodeling as well as channels implicated in water and solute flux, such as aquaporins. In addition, we were able to determine a robust consensus ASR1-binding DNA motif.

CONCLUSIONS

The finding of cell wall synthesis and aquaporin genes as targets of ASR1 is consistent with their suggested role in the physiological adaptation of plants to water loss. The results gain insight into the environmental stress-sensing pathways leading to plant tolerance of drought.

Comparative genome analysis of Solanum lycopersicum and Solanum tuberosum

Solanum lycopersicum and Solanum tuberosum are agriculturally important crop species as they are rich sources of starch, protein, antioxidants, lycopene, beta-carotene, vitamin C, and fiber. The genomes of S. lycopersicum and S. tuberosum are currently available. However the linear strings of nucleotides that together comprise a genome sequence are of limited significance by themselves. Computational and bioinformatics approaches can be used to exploit the genomes for fundamental research for improving their varieties. The comparative genome analysis, Pfam analysis of predicted reviewed paralogous proteins was performed. It was found that S. lycopersicum proteins belong to more families, domains and clans in comparison with S. tuberosum. It was also found that mostly intergenic regions are conserved in two genomes followed by exons, intron and UTR. This can be exploited to predict regions between genomes that are similar to each other and to study the evolutionary relationship between two genomes, leading towards the development of disease resistance, stress tolerance and improved varieties of tomato.

Individual chromosome identification, chromosomal collinearity and genetic-physical integrated map in Gossypium darwinii and four D genome cotton species revealed by BAC-FISH

The study was conducted on the individual chromosome identification in Gossypium darwinii (A(d)D(d)), G. klotzschianum (D(3k)), G. davidsonii (D(3d)), G. armourianum (D(2-1)) and G. aridum (D(4)) using a multi-probe fluorescence of in situ hybridization (FISH) system. Comparative analysis on their genetic maps with that of physical maps was made as well. The FISH probes used contained two sets of bacterial artificial chromosome (BAC) clones [one is specific to 26 individual chromosomes from A and D subgenomes of G. hirsutum (A(h) and D(h)) while the other is a D genome centromere-specific BAC clone 150D24], 45S and 5S rDNA clones. The results showed that all A(d) chromosomes were marked with the 13 A(h) chromosome-specific BAC clones, whilst all D(d), D(3k), D(3d), D(2-1) and D(4) chromosomes and chromosomal arms were identified with the 13 D(h) chromosome-specific BAC clones and the D genome centromere-specific BAC. According to the homology within D subgenomes which are between A (D) genome and A (D) subgenome, the systematic nomenclature for individual chromosome in the five species was established. The chromosomes of A (D) subgenomes in G. darwinii were named as A(d)01-A(d)13 (D(d)01-D(d)13). The chromosomes in D(3k), D(3d), D(2-1) and D(4) were named as D(3k)01-D(3k)13, D(3d)01-D(3d)13, D(2-1)01-D(2-1)13 and D(4)01-D(4)13, respectively. Based on the successful identification for individual chromosomes, 45S and 5S rDNA were located to the special chromosomes and chromosomal arms for all five species. And there appeared chromosomal collinearity from the BAC clones among different species by comparing BACs positions, which suggested that the majority of chromosome segment homology may exist between D genomes and D subgenome. Moreover, as the genetic map and physical map were integrated, the orientations of genetic maps for D(d) and D genomes of diploid cotton were established. The orientations of some of chromosomes in genetic maps (D(d)03, D(d)04, D(d)06, D(d)09, D(d)10 and D(d)12) were found switched. The SSR marker in the middle of linkage group 04 was corrected at nearby the end of chromosome 04 by FISH. The study will be helpful to establish a theoretical basis using the wild gene bank to exploit more genes aiming for cotton breeding and will provide powerful evidences both for the evolution of Gossypium and assembling the sequences of the obtained and as well the on-going whole genome sequencing projects of cotton.

Integration of mate pair sequences to improve shotgun assemblies of flow-sorted chromosome arms of hexaploid wheat

BACKGROUND

The assembly of the bread wheat genome sequence is challenging due to allohexaploidy and extreme repeat content (>80%). Isolation of single chromosome arms by flow sorting can be used to overcome the polyploidy problem, but the repeat content cause extreme assembly fragmentation even at a single chromosome level. Long jump paired sequencing data (mate pairs) can help reduce assembly fragmentation by joining multiple contigs into single scaffolds. The aim of this work was to assess how mate pair data generated from multiple displacement amplified DNA of flow-sorted chromosomes affect assembly fragmentation of shotgun assemblies of the wheat chromosomes.

RESULTS

Three mate pair (MP) libraries (2 Kb, 3 Kb, and 5 Kb) were sequenced to a total coverage of 89x and 64x for the short and long arm of chromosome 7B, respectively. Scaffolding using SSPACE improved the 7B assembly contiguity and decreased gene space fragmentation, but the degree of improvement was greatly affected by scaffolding stringency applied. At the lowest stringency the assembly N50 increased by ~7 fold, while at the highest stringency N50 was only increased by ~1.5 fold. Furthermore, a strong positive correlation between estimated scaffold reliability and scaffold assembly stringency was observed. A 7BS scaffold assembly with reduced MP coverage proved that assembly contiguity was affected only to a small degree down to ~50% of the original coverage.

CONCLUSION

The effect of MP data integration into pair end shotgun assemblies of wheat chromosome was moderate; possibly due to poor contig assembly contiguity, the extreme repeat content of wheat, and the use of amplified chromosomal DNA for MP library construction.

Interstitial telomeric repeats are enriched in the centromeres of chromosomes in Solanum species

Interstitial telomeric repeats (ITRs) were reported in a number of animal and plant species. Most ITRs are organized as short tandem arrays and are likely evolutionary relics derived from chromosomal rearrangements and DNA repairs. However, megabase-sized ITR arrays were reported in Solanum species. Here, we report a fluorescence in situ hybridization (FISH) survey of ITRs in all representative diploid Solanum species, including potato, tomato, and eggplant. FISH revealed massive amplification of ITRs in the centromeric regions of chromosomes from the Solanum species containing the B and P genomes. A significant proportion of the ITR FISH signals was mapped within the primary constrictions of the pachytene chromosomes of Solanum pinnatisectum (B genome). In addition, some ITR sites overlapped with St49, a satellite repeat enriched in centromeric DNA sequences associated with CENH3 nucleosomes, in both A and B genome Solanum species. These results show that some ITR subfamilies have been amplified and invaded in the functional centromeres of chromosomes in Solanum species.

A hybrid BAC physical map of potato: a framework for sequencing a heterozygous genome

BACKGROUND

Potato is the world's third most important food crop, yet cultivar improvement and genomic research in general remain difficult because of the heterozygous and tetraploid nature of its genome. The development of physical map resources that can facilitate genomic analyses in potato has so far been very limited. Here we present the methods of construction and the general statistics of the first two genome-wide BAC physical maps of potato, which were made from the heterozygous diploid clone RH89-039-16 (RH).

RESULTS

First, a gel electrophoresis-based physical map was made by AFLP fingerprinting of 64478 BAC clones, which were aligned into 4150 contigs with an estimated total length of 1361 Mb. Screening of BAC pools, followed by the KeyMaps in silico anchoring procedure, identified 1725 AFLP markers in the physical map, and 1252 BAC contigs were anchored the ultradense potato genetic map. A second, sequence-tag-based physical map was constructed from 65919 whole genome profiling (WGP) BAC fingerprints and these were aligned into 3601 BAC contigs spanning 1396 Mb. The 39733 BAC clones that overlap between both physical maps provided anchors to 1127 contigs in the WGP physical map, and reduced the number of contigs to around 2800 in each map separately. Both physical maps were 1.64 times longer than the 850 Mb potato genome. Genome heterozygosity and incomplete merging of BAC contigs are two factors that can explain this map inflation. The contig information of both physical maps was united in a single table that describes hybrid potato physical map.

CONCLUSIONS

The AFLP physical map has already been used by the Potato Genome Sequencing Consortium for sequencing 10% of the heterozygous genome of clone RH on a BAC-by-BAC basis. By layering a new WGP physical map on top of the AFLP physical map, a genetically anchored genome-wide framework of 322434 sequence tags has been created. This reference framework can be used for anchoring and ordering of genomic sequences of clone RH (and other potato genotypes), and opens the possibility to finish sequencing of the RH genome in a more efficient way via high throughput next generation approaches.

Sequencing crop genomes: approaches and applications

Many challenges face plant scientists, in particular those working on crop production, such as a projected increase in population, decrease in water and arable land, changes in weather patterns and predictability. Advances in genome sequencing and resequencing can and should play a role in our response to meeting these challenges. However, several barriers prevent rapid and effective deployment of these tools to a wide variety of crops. Because of the complexity of crop genomes, de novo sequencing with next-generation sequencing technologies is a process fraught with difficulties that then create roadblocks to the utilization of these genome sequences for crop improvement. Collecting rapid and accurate phenotypes in crop plants is a hindrance to integrating genomics with crop improvement, and advances in informatics are needed to put these tools in the hands of the scientists on the ground.

Organization and evolution of subtelomeric satellite repeats in the potato genome

Subtelomeric domains immediately adjacent to telomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes. A common feature associated with subtelomeric regions in different eukaryotes is the presence of long arrays of tandemly repeated satellite sequences. However, studies on molecular organization and evolution of subtelomeric repeats are rare. We isolated two subtelomeric repeats, CL14 and CL34, from potato (Solanum tuberosum). The CL14 and CL34 repeats are organized as independent long arrays, up to 1-3 Mb, of 182 bp and 339 bp monomers, respectively. The CL14 and CL34 repeat arrays are directly connected with the telomeric repeats at some chromosomal ends. The CL14 repeat was detected at the subtelomeric regions among highly diverged Solanum species, including tomato (Solanum lycopersicum). In contrast, CL34 was only found in potato and its closely related species. Interestingly, the CL34 repeat array was always proximal to the telomeres when both CL14 and CL34 were found at the same chromosomal end. In addition, the CL34 repeat family showed more sequence variability among monomers compared with the CL14 repeat family. We conclude that the CL34 repeat family emerged recently from the subtelomeric regions of potato chromosomes and is rapidly evolving. These results provide further evidence that subtelomeric domains are among the most dynamic regions in eukaryotic genomes.

Spectral counting robust on high mass accuracy mass spectrometers

Mass spectrometry is central to shotgun proteomics, an application that seeks to quantify as much of the total protein complement of a biological sample as possible. The high mass accuracy, resolution, capacity and scan rate of modern mass spectrometers have greatly facilitated this endeavor. The sum of MS to MS/MS transitions in tandem mass spectrometry, the spectral count (SC), of a peptide has been shown to be a reliable estimate of its relative abundance. However, when using SCs, optimal MS configurations are crucial in order to maximize the number of low abundant proteins quantified while keeping the estimates for the highly abundant proteins within the linear dynamic range.In this study, LC/MS/MS analysis was performed using an LTQ-OrbiTrap on a sample containing many highly abundant proteins. Tuning the LTQ-OrbiTrap mass spectrometer to minimize redundant MS/MS acquisition and to maximize resolution of the proteome by accurately measured m/z ratios resulted in an appreciable increase in quantified low abundant proteins. An exclusion duration of 90 s and an exclusion width of 10 ppm were found best of those tested. The spectral count of individual proteins was found to be highly reproducible and protein abundance ratios were not affected by the different settings that were applied. We conclude that on a high mass accuracy instrument spectral counting is a robust measure of protein abundance even for samples containing many highly abundant proteins and that tuning dynamic exclusion parameters appreciably improves the number of proteins that can be reliably quantified.

Characterization of microsatellites and gene contents from genome shotgun sequences of mungbean (Vigna radiata (L.) Wilczek)

BACKGROUND

Mungbean is an important economical crop in Asia. However, genomic research has lagged behind other crop species due to the lack of polymorphic DNA markers found in this crop. The objective of this work is to develop and characterize microsatellite or simple sequence repeat (SSR) markers from genome shotgun sequencing of mungbean.

RESULT

We have generated and characterized a total of 470,024 genome shotgun sequences covering 100.5 Mb of the mungbean (Vigna radiata (L.) Wilczek) genome using 454 sequencing technology. We identified 1,493 SSR motifs that could be used as potential molecular markers. Among 192 tested primer pairs in 17 mungbean accessions, 60 loci revealed polymorphism with polymorphic information content (PIC) values ranging from 0.0555 to 0.6907 with an average of 0.2594. Majority of microsatellite markers were transferable in Vigna species, whereas transferability rates were only 22.90% and 24.43% in Phaseolus vulgaris and Glycine max, respectively. We also used 16 SSR loci to evaluate phylogenetic relationship of 35 genotypes of the Asian Vigna group. The genome survey sequences were further analyzed to search for gene content. The evidence suggested 1,542 gene fragments have been sequence tagged, that fell within intersected existing gene models and shared sequence homology with other proteins in the database. Furthermore, potential microRNAs that could regulate developmental stages and environmental responses were discovered from this dataset.

CONCLUSION

In this report, we provided evidence of generating remarkable levels of diverse microsatellite markers and gene content from high throughput genome shotgun sequencing of the mungbean genomic DNA. The markers could be used in germplasm analysis, accessing genetic diversity and linkage mapping of mungbean.

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.

Faithful inheritance of cytosine methylation patterns in repeated sequences of the allotetraploid tobacco correlates with the expression of DNA methyltransferase gene families from both parental genomes

The widespread occurrence of epigenetic alterations in allopolyploid species deserves scrutiny that DNA methylation systems may be perturbed by interspecies hybridization and polyploidization. Here we studied the genes involved in DNA methylation in Nicotiana tabacum (tobacco) allotetraploid containing S and T genomes inherited from Nicotiana sylvestris and Nicotiana tomentosiformis progenitors. To determine the inheritance of DNA methyltransferase genes and their expression patterns we examined three major DNA methyltransferase families (MET1, CMT3 and DRM) from tobacco and the progenitor species. Using Southern blot hybridization and PCR-based methods (genomic CAPS), we found that the parental loci of these gene families are retained in tobacco. Homoeologous expression was found in all tissues examined (leaf, root, flower) suggesting that DNA methyltransferase genes were probably not themselves targets of uniparental epigenetic silencing for over thousands of generations of allotetraploid evolution. The level of CG and CHG methylation of selected high-copy repeated sequences was similar and high in tobacco and its diploid progenitors. We speculate that natural selection might favor additive expression of parental DNA methyltransferase genes maintaining high levels of DNA methylation in tobacco, which has a repeat-rich heterochromatic genome.

Single nucleotide polymorphisms in the allene oxide synthase 2 gene are associated with field resistance to late blight in populations of tetraploid potato cultivars

The oomycete Phytophthora infestans causes late blight, the most relevant disease of potato (Solanum tuberosum) worldwide. Field resistance to late blight is a complex trait. When potatoes are cultivated under long day conditions in temperate climates, this resistance is correlated with late plant maturity, an undesirable characteristic. Identification of natural gene variation underlying late blight resistance not compromised by late maturity will facilitate the selection of resistant cultivars and give new insight in the mechanisms controlling quantitative pathogen resistance. We tested 24 candidate loci for association with field resistance to late blight and plant maturity in a population of 184 tetraploid potato individuals. The individuals were genotyped for 230 single nucleotide polymorphisms (SNPs) and 166 microsatellite alleles. For association analysis we used a mixed model, taking into account population structure, kinship, allele substitution and interaction effects of the marker alleles at a locus with four allele doses. Nine SNPs were associated with maturity corrected resistance (P < 0.001), which collectively explained 50% of the genetic variance of this trait. A major association was found at the StAOS2 locus encoding allene oxide synthase 2, a key enzyme in the biosynthesis of jasmonates, plant hormones that function in defense signaling. This finding supports StAOS2 as being one of the factors controlling natural variation of pathogen resistance.

Characterization of a deep-coverage carrot (Daucus carota L.) BAC library and initial analysis of BAC-end sequences

Carrot is the most economically important member of the Apiaceae family and a major source of provitamin A carotenoids in the human diet. However, carrot molecular resources are relatively underdeveloped, hampering a number of genetic studies. Here, we report on the synthesis and characterization of a bacterial artificial chromosome (BAC) library of carrot. The library is 17.3-fold redundant and consists of 92,160 clones with an average insert size of 121 kb. To provide an overview of the composition and organization of the carrot nuclear genome we generated and analyzed 2,696 BAC-end sequences (BES) from nearly 2,000 BACs, totaling 1.74 Mb of BES. This analysis revealed that 14% of the BES consists of known repetitive elements, with transposable elements representing more than 80% of this fraction. Eleven novel carrot repetitive elements were identified, covering 8.5% of the BES. Analysis of microsatellites showed a comparably low frequency for these elements in the carrot BES. Comparisons of the translated BES with protein databases indicated that approximately 10% of the carrot genome represents coding sequences. Moreover, among eight dicot species used for comparison purposes, carrot BES had highest homology to protein-coding sequences from tomato. This deep-coverage library will aid carrot breeding and genetics.

Chromatin structure and physical mapping of chromosome 6 of potato and comparative analyses with tomato

Potato (Solanum tuberosum) has the densest genetic linkage map and one of the earliest established cytogenetic maps among all plant species. However, there has been limited effort to integrate these maps. Here, we report fluorescence in situ hybridization (FISH) mapping of 30 genetic marker-anchored bacterial artificial chromosome (BAC) clones on the pachytene chromosome 6 of potato. The FISH mapping results allowed us to define the genetic positions of the centromere and the pericentromeric heterochromatin and to relate chromatin structure to the distribution of recombination along the chromosome. A drastic reduction of recombination was associated with the pericentromeric heterochromatin that accounts for approximately 28% of the physical length of the pachytene chromosome. The pachytene chromosomes 6 of potato and tomato (S. lycopersicum) share a similar morphology. However, distinct differences of heterochromatin distribution were observed between the two chromosomes. FISH mapping of several potato BACs on tomato pachytene chromosome 6 revealed an overall colinearity between the two chromosomes. A chromosome inversion was observed in the euchromatic region of the short arms. These results show that the potato and tomato genomes contain more chromosomal rearrangements than those reported previously on the basis of comparative genetic linkage mapping.