Orthologous comparison in a gene-rich region among grasses reveals stability in the sugarcane polyploid genome

Modern sugarcane (Saccharum spp.) is an important grass that contributes 60% of the raw sugar produced worldwide and has a high biofuel production potential. It was created about a century ago through hybridization of two highly polyploid species, namely S. officinarum and S. spontaneum. We investigated genome dynamics in this highly polyploid context by analyzing two homoeologous sequences (97 and 126 kb) in a region that has already been studied in several cereals. Our findings indicate that the two Saccharum species diverged 1.5-2 million years ago from one another and 8-9 million years ago from sorghum. The two sugarcane homoeologous haplotypes show perfect colinearity as well as high gene structure conservation. Apart from the insertion of a few retrotransposable elements, high homology was also observed for the non-transcribed regions. Relative to sorghum, the sugarcane sequences displayed colinearity, with the exception of two genes present only in sorghum, and striking homology in most non-coding parts of the genome. The gene distribution highlighted high synteny and colinearity with rice, and partial colinearity with each homoeologous maize region, which became perfect when the sequences were combined. The haplotypes observed in sugarcane may thus closely represent the ancestral Andropogoneae haplotype. This analysis of sugarcane haplotype organization at the sequence level suggests that the high ploidy in sugarcane did not induce generalized reshaping of its genome, thus challenging the idea that polyploidy quickly induces generalized rearrangement of genomes. These results also confirm the view that sorghum is the model of choice for sugarcane.     

Three founding ancestral genomes involved in the origin of sugarcane

Background and AimsModern sugarcane cultivars (Saccharum spp.) are high polyploids, aneuploids (2n = ~12x = ~120) derived from interspecific hybridizations between the domesticated sweet species Saccharum officinarum and the wild species S. spontaneum.MethodsTo analyse the architecture and origin of such a complex genome, we analysed the sequences of all 12 hom(oe)ologous haplotypes (BAC clones) from two distinct genomic regions of a typical modern cultivar, as well as the corresponding sequence in Miscanthus sinense and Sorghum bicolor, and monitored their distribution among representatives of the Saccharum genus.Key ResultsThe diversity observed among haplotypes suggested the existence of three founding genomes (A, B, C) in modern cultivars, which diverged between 0.8 and 1.3 Mya. Two genomes (A, B) were contributed by S. officinarum; these were also found in its wild presumed ancestor S. robustum, and one genome (C) was contributed by S. spontaneum. These results suggest that S. officinarum and S. robustum are derived from interspecific hybridization between two unknown ancestors (A and B genomes). The A genome contributed most haplotypes (nine or ten) while the B and C genomes contributed one or two haplotypes in the regions analysed of this typical modern cultivar. Interspecific hybridizations likely involved accessions or gametes with distinct ploidy levels and/or were followed by a series of backcrosses with the A genome. The three founding genomes were found in all S. barberi, S. sinense and modern cultivars analysed. None of the analysed accessions contained only the A genome or the B genome, suggesting that representatives of these founding genomes remain to be discovered.ConclusionsThis evolutionary model, which combines interspecificity and high polyploidy, can explain the variable chromosome pairing affinity observed in Saccharum. It represents a major revision of the understanding of Saccharum diversity.     

Genomic organization of the S core region and the S flanking regions of a class-II S haplotype in Brassica rapa

The nucleotide sequence of an 86.4-kb region that includes the SP11, SRK, and SLG genes of Brassica rapa S-60 (a class-II S haplotype) was determined. In the sequenced region, 13 putative genes were found besides SP11-60, SRK-60, and SLG-60. Five of these sequences were isolated as cDNAs, five were homologues of known genes, cDNAs, or ORFs, and three are hypothetical ORFs. Based on their nucleotide sequences, however, some of them are thought to be non-functional. Two regions of colinearity between the class-II S-60 and Brassica class-I S haplotypes were identified, i.e., S flanking region 1 which shows partial colinearity of non-genic sequences and S flanking region 2 which shows a high level of colinearity. The observed colinearity made it possible to compare the order of SP-11, SRK, and SLG genes in the S locus between the five sequenced S haplotypes. It emerged that the order of SRK and SLG in class-II S-60 is the reverse of that in the four class-I S haplotypes reported so far, and the order of SP11, SRK and SLG is the opposite of that in the class-I haplotype S-910. The possible gene designated as SAN1 (S locus Anther-expressed Non-coding RNA like-1), which is located in the region between SP11-60 and SRK-60, has features reminiscent of genes for non-coding RNAs (ncRNAs), but no homologous sequences were found in the databases. This sequence is transcribed in anthers but not in stigmas or leaves. These features of the genomic structure of S-60 are discussed with special reference to the characteristics of class-II S haplotypes.     

A leucine tRNA gene adjacent to the QA gene cluster of Neurospora crassa.

A single tRNALeu gene has been localized and sequenced from Neurospora crassa. It is located only 375 bp from the qa gene cluster and it is the only tRNA or 5S rRNA gene within this cloned 37 kb region. The gene encodes a tRNALeu with the anti-codon AAG, and unlike the other nuclear eukaryotic tRNALeu (AAG) gene sequenced (from C. elegans), contains an intervening sequence of 27 bp. The Neurospora tRNALeu (AAG) is 84% and 73% homologous respectively to the C. elegans and bovine tRNALeu (AAG), and is 84% homologous to a Drosophila tRNALeu (CAA). However, it is only 65% homologous to a yeast tRNALeu (CAA) and there is little conservation of intervening sequences or V-loop regions. The gene hybridizes to at least 16 other DNA fragments in the Neurospora genome. Its expression does not seem to be linked to that of the qa genes.     

编码腈水解酶的大豆疫霉基因PsNIA的克隆与转录分析

包括大豆在内的许多植物都可以产生氰化物,对侵染的病原菌产生毒害作用而阻碍其进一步扩展。采用抑制性差减杂交(suppression subtractive hybridization,SSH)的方法,筛选到一个在大豆疫霉侵染早期上调表达的、编码腈水解酶的cDNA片段;克隆了该基因的全长序列,命名为PsNIA。Southern杂交结果显示,PsNIA在大豆疫霉基因组中只有1个拷贝。系统发育分析表明,PsNIA与绿脓杆菌Pseudomonas aeruginosa的腈水解酶的序列同源性最高,且该基因编码的氨基酸序列具有腈水解酶的保守结构域。RT-PCR分析表明,该基因在大豆疫霉侵染大豆12h时可以检测到转录。     

Antisense-Mediated Depletion of Tomato Chloroplast Omega-3 Fatty Acid Desaturase Enhances Thermal Tolerance

A chloroplast-localized tomato （Lycopersicon esculentum Mill.） ω-3 fatty acid desaturase gene （LeFADT） was isolated and characterized with regard to its sequence, response to various temperatures, and function in antisense transgenic tomato plants. The deduced amino acid sequence had four histidine-rich regions, of which three regions were highly conserved throughout the whole ω-3 fatty acid desaturasegene family. Southern blotting analysis showed that LeFAD7was encoded by a single copy gene and had two homologous genes in the tomato genome. Northern blot showed that LeFAD7 was expressed in all organs and was especially abundant in leaf tissue. Meanwhile, expression of LeFAD7 was induced by chilling stress （4 ℃）, but was inhibited by high temperature （45 ℃）, in leaves. Transgenic tomato plants were produced by integration of the antisense LeFAD7DNA under the control of a CaMV35S promoter into the genome. Antisense transgenic plants with lower 18 ： 3 content could maintain a higher maximal photochemical efficiency （Fv/Fm） and O2 evolution rate than wild-type plants. These results suggested that silence of the LeFAD7 gene alleviated high-temperature stress. There was also a correlation between the low content of 18 ： 3 resulting from silence of the LeFAD7 gene and tolerance to high-temperature stress.     

Detailed alignment of saccharum and sorghum chromosomes: comparative organization of closely related diploid and polyploid genomes.

The complex polyploid genomes of three Saccharum species have been aligned with the compact diploid genome of Sorghum (2n = 2x = 20). A set of 428 DNA probes from different Poaceae (grasses) detected 2460 loci in F1 progeny of the crosses Saccharum officinarum Green German x S. spontaneum IND 81-146, and S. spontaneum PIN 84-1 x S. officinarum Muntok Java. Thirty-one DNA probes detected 226 loci in S. officinarum LA Purple x S. robustum Molokai 5829. Genetic maps of the six Saccharum genotypes, including up to 72 linkage groups, were assembled into "homologous groups" based on parallel arrangements of duplicated loci. About 84% of the loci mapped by 242 common probes were homologous between Saccharum and Sorghum. Only one interchromosomal and two intrachromosomal rearrangements differentiated both S. officinarum and S. spontaneum from Sorghum, but 11 additional cases of chromosome structural polymorphism were found within Saccharum. Diploidization was advanced in S. robustum, incipient in S. officinarum, and absent in S. spontaneum, consistent with biogeographic data suggesting that S. robustum is the ancestor of S. officinarum, but raising new questions about the antiquity of S. spontaneum. The densely mapped Sorghum genome will be a valuable tool in ongoing molecular analysis of the complex Saccharum genome.     

Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.     

The Sugarcane Genome Challenge: Strategies for Sequencing a Highly Complex Genome

Sugarcane cultivars derive from interspecific hybrids obtained by crossing Saccharum officinarum and Saccharum spontaneum and provide feedstock used worldwide for sugar and biofuel production. The importance of sugarcane as a bioenergy feedstock has increased interest in the generation of new cultivars optimised for energy production. Cultivar improvement has relied largely on traditional breeding methods, which may be limited by the complexity of inheritance in interspecific polyploid hybrids, and the time-consuming process of selection of plants with desired agronomic traits. In this sense, molecular genetics can assist in the process of developing improved cultivars by generating molecular markers that can be used in the breeding process or by introducing new genes into the sugarcane genome. For meeting each of these, and additional goals, biotechnologists would benefit from a reference genome sequence of a sugarcane cultivar. The sugarcane genome poses challenges that have not been addressed in any prior sequencing project, due to its highly polyploid and aneuploid genome structure with a complete set of homeologous genes predicted to range from 10 to 12 copies (alleles) and to include representatives from each of two different species. Although sugarcane’s monoploid genome is about 1 Gb, its highly polymorphic nature represents another significant challenge for obtaining a genuine assembled monoploid genome. With a rich resource of expressed-sequence tag (EST) data in the public domain, the present article describes tools and strategies that may aid in the generation of a reference genome sequence.     

Enrichment of genomic DNA for polymorphism detection in a non-model highly polyploid crop plant

Large polyploid genomes of non-model species remain challenging targets for DNA polymorphism discovery despite the increasing throughput and continued reductions in cost of sequencing with new technologies. For these species especially, there remains a requirement to enrich genomic DNA to discover polymorphisms in regions of interest because of large genome size and to provide the sequence depth to enable estimation of copy number. Various methods of enriching DNA have been utilised, but some recent methods enable the efficient sampling of large regions (e.g. the exome). We have utilised one of these methods, solution-based hybridization (Agilent SureSelect), to capture regions of the genome of two sugarcane genotypes (one Saccharum officinarum and one Saccharum hybrid) based mainly on gene sequences from the close relative Sorghum bicolor. The capture probes span approximately 5.8?megabases (Mb). The enrichment over whole-genome shotgun sequencing was 10-11-fold for the two genotypes tested. This level of enrichment has important consequences for detecting single nucleotide polymorphisms (SNPs) from a single lane of Illumina (Genome Analyzer) sequence reads. The detection of polymorphisms was enabled by the depth of sequence at or near probe sites and enabled the detection of 270?000-280?000 SNPs within each genotype from a single lane of sequence using stringent detection parameters. The SNPs were present in 13?000-16?000 targeted genes, which would enable mapping of a large number of these chosen genes. SNP validation from 454 sequencing and between-genotype confirmations gave an 87%-91% validation rate.     

The identification and characterisation of alleles of sucrose phosphate synthase gene family III in sugarcane

Little is known about the extent of allelic diversity of genes in the complex polyploid, sugarcane. Using sucrose phosphate synthase (SPS) Gene (SPS) Family III as an example, we have amplified and sequenced a 400 nt region from this gene from two sugarcane lines that are parents of a mapping population. Ten single nucleotide polymorphisms (SNPs) were identified within the 400 nt region of which seven were present in both lines. In the elite commercial cultivar Q165^A, 10 sequence haplotypes were identified, with four haplotypes recovered at 9% or greater frequency. Based on SNP presence, two clusters of haplotypes were observed. In IJ76-514, a Saccharum officinarum accession, 8 haplotypes were identified with 4 haplotypes recovered at 13% or greater frequency. Again, two clusters of haplotypes were observed. The results suggest that there may be two SPS Gene Family III genes per genome in sugarcane, each with different numbers of different alleles. This suggestion is supported by sequencing results in an elite parental sorghum line, 403463-2-1, in which 4 haplotypes, corresponding to two broad types, were also identified. Primers were designed to the sugarcane SNPs and screened over bulked DNA from high and low Sucrose-containing progeny from a cross between Q165^A and IJ76-514. The SNP frequency did not vary in the two bulked DNA samples, suggesting that these SNPs from this SPS gene family are not associated with variation in sucrose content. Using an ecotilling approach, two of the SPS Gene Family III haplotypes were mapped to two different linkage groups in homology group 1 in Q165^A. Both haplotypes mapped near QTLs for increased sucrose content but were not themselves associated with any sugar-related trait.     

Recurrent deletions of puroindoline genes at the grain hardness locus in four independent lineages of polyploid wheat

Polyploidy is known to induce numerous genetic and epigenetic changes but little is known about their physiological bases. In wheat, grain texture is mainly determined by the Hardness (Ha) locus consisting of genes Puroindoline a (Pina) and b (Pinb). These genes are conserved in diploid progenitors but were deleted from the A and B genomes of tetraploid Triticum turgidum (AB). We now report the recurrent deletions of Pina-Pinb in other lineages of polyploid wheat. We analyzed the Ha haplotype structure in 90 diploid and 300 polyploid accessions of Triticum and Aegilops spp. Pin genes were conserved in all diploid species and deletion haplotypes were detected in all polyploid Triticum and most of the polyploid Aegilops spp. Two Pina-Pinb deletion haplotypes were found in hexaploid wheat (Triticum aestivum; ABD). Pina and Pinb were eliminated from the G genome, but maintained in the A genome of tetraploid Triticum timopheevii (AG). Subsequently, Pina and Pinb were deleted from the A genome but retained in the A(m) genome of hexaploid Triticum zhukovskyi (A(m)AG). Comparison of deletion breakpoints demonstrated that the Pina-Pinb deletion occurred independently and recurrently in the four polyploid wheat species. The implications of Pina-Pinb deletions for polyploid-driven evolution of gene and genome and its possible physiological significance are discussed.     

假单胞菌M18rsmA-突变株的构建及其对Plt和PCA合成的区别性调控作用

假单胞菌(Pseudomonas sp．)M18是促进植物生长的根际细菌,能产生吩嗪-1-羧酸(PCA)和藤黄绿菌素(Plt)两种不同的抗生素抑制植物病原菌,保护植物免受病害。运用PCR方法,从M18基因组中,扩增出rsmA基因部分片段,并以该片段为探针,从M18的基因组柯斯文库中筛出阳性克隆,切取带有rsmA基因及两侧序列的1．5kb片段,中间插入编码Km‘的DNA片段,获得rsmA^-体外突变体。运用同源重组剔除技术,构建了M18菌株的rsmA突变株M18R^-。突变株M18R^-生物合成Plt的能力比野生型M18提高4倍,但是,PCA产量仅为野生型的20％。研究结果表明,全局性调控基因rsmA可能通过不同的机制区别性地影响Plt和PCA的生物合成。     

Colinearity and gene density in grass genomes

Grasses are the single most important plant family in agriculture. In the past years, comparative genetic mapping has revealed conserved gene order (colinearity) among many grass species. Recently, the first studies at gene level have demonstrated that microcolinearity of genes is less conserved: small scale rearrangements and deletions complicate the microcolinearity between closely related species, such as sorghum and maize, but also between rice and other crop plants. In spite of these problems, rice remains the model plant for grasses as there is limited useful colinearity between Arabidopsis and grasses. However, studies in rice have to be complemented by more intensive genetic work on grass species with large genomes (maize, Triticeae). Gene-rich chromosomal regions in species with large genomes, such as wheat, have a high gene density and are ideal targets for partial genome sequencing.     

Identification and expression analysis of microRNAs and targets in the biofuel crop sugarcane

Background  

MicroRNAs (miRNAs) are small regulatory RNAs, some of which are conserved in diverse plant genomes. Therefore, computational identification and further experimental validation of miRNAs from non-model organisms is both feasible and instrumental for addressing miRNA-based gene regulation and evolution. Sugarcane (Saccharum spp.) is an important biofuel crop with publicly available expressed sequence tag and genomic survey sequence databases, but little is known about miRNAs and their targets in this highly polyploid species.     

Single-Strand Conformational Polymorphism of EST-SSRs: A Potential Tool for Diversity Analysis and Varietal Identification in Sugarcane

Sugarcane (Saccharum spp.) has a large complex polyploid genome. Assay of molecular variation in the expressed component of its genome has relevance to the analysis of genetic diversity, variety identification and introgression of agronomically useful genes present in different members of the Saccharum complex. The present study was designed to evaluate single-strand conformational polymorphism (SSCP) as a potential tool to detect genetic variation in the expressed sequence tag (EST) derived microsatellites. Twenty primer pairs obtained from EST libraries and one designed from soluble acid invertase gene sequence were used to characterise 21 clones belonging to four different Saccharum species and 22 sugarcane varieties/genotypes. All the markers, including the two, which were reported monomorphic even at the interspecific level in an automated fragment analysis system in a previous study, could be successfully converted into polymorphic ones using SSCP analysis. A broad range of variation could be revealed by this technique. The Saccharum spp. clones could be grouped into distinct clusters, confirming the species relationships postulated earlier using morphological, biochemical and molecular methods. The polymorphic markers could also differentiate all the 22 sugarcane varieties from each other. This is a first report that demonstrates the usefulness of SSCP technique, in obtaining polymorphic microsatellite markers developed from EST sequences for various genetic and breeding applications, in this polyploid species.     

Seventy million years of concerted evolution of a homoeologous chromosome pair, in parallel, in major Poaceae lineages

Whole genome duplication ~70 million years ago provided raw material for Poaceae (grass) diversification. Comparison of rice (Oryza sativa), sorghum (Sorghum bicolor), maize (Zea mays), and Brachypodium distachyon genomes revealed that one paleo-duplicated chromosome pair has experienced very different evolution than all the others. For tens of millions of years, the two chromosomes have experienced illegitimate recombination that has been temporally restricted in a stepwise manner, producing structural stratification in the chromosomes. These strata formed independently in different grass lineages, with their similarities (low sequence divergence between paleo-duplicated genes) preserved in parallel for millions of years since the divergence of these lineages. The pericentromeric region of this homeologous chromosome pair accounts for two-thirds of the gene content differences between the modern chromosomes. Both intriguing and perplexing is a distal chromosomal region with the greatest DNA similarity between surviving duplicated genes but also with the highest concentration of lineage-specific gene pairs found anywhere in these genomes and with a significantly elevated gene evolutionary rate. Intragenomic similarity near this chromosomal terminus may be important in hom(e)ologous chromosome pairing. Chromosome structural stratification, together with enrichment of autoimmune response-related (nucleotide binding site-leucine-rich repeat) genes and accelerated DNA rearrangement and gene loss, confer a striking resemblance of this grass chromosome pair to the sex chromosomes of other taxa.     

Contributions of domesticated plant studies to our understanding of plant evolution

BACKGROUND: Plant evolutionary theory has been greatly enriched by studies on crop species. Over the last century, important information has been generated on many aspects of population biology, speciation and polyploid genetics. SCOPE: Searches for quantitative trait loci (QTL) in crop species have uncovered numerous blocks of genes that have dramatic effects on adaptation, particularly during the domestication process. Many of these QTL have epistatic and pleiotropic effects making rapid evolutionary change possible. Most of the pioneering work on the molecular basis of self-incompatibility has been conducted on crop species, along with the sequencing of the phytopathogenic resistance genes (R genes) responsible for the 'gene-to-gene' relations of coevolution observed in host-pathogen relationships. Some of the better examples of co-adaptation and early acting inbreeding depression have also been elucidated in crops. Crop-wild progenitor interactions have provided rich opportunities to study the evolution of novel adaptations subsequent to hybridization. Most crop/wild F1 hybrids have reduced fitness, but in some instances the crop relatives have acquired genes that make them more efficient weeds through crop mimicry. Studies on autopolyploid alfalfa and potato have uncovered the means by which polyploid gametes are formed and have led to hypotheses about how multiallelic interactions are associated with fitness and self-fertility. Research on the cole crops and wheat has discovered that newly formed polyploids can undergo dramatic genome rearrangements that could lead to rapid evolutionary change. CONCLUSIONS: Many more important evolutionary discoveries are on the horizon, now that the whole genome sequence is available of the two major subspecies of rice Oryza sativa ssp. japonica and O. sativa ssp. indica. The rice sequence data can be used to study the origin of genes and gene families, track rates of sequence divergence over time, and provide hints about how genes evolve and generate products with novel biological properties. The rice sequence data has already been mined to show that transposable elements often carry fragments of cellular genes. This type of genome shuffling could play a role in creating novel, reorganized genes with new adaptive properties.     

Single-copy genes define a conserved order between rice and wheat for understanding differences caused by duplication, deletion, and transposition of genes

The high-quality rice genome sequence is serving as a reference for comparative genome analysis in crop plants, especially cereals. However, early comparisons with bread wheat showed complex patterns of conserved synteny (gene content) and colinearity (gene order). Here, we show the presence of ancient duplicated segments in the progenitor of wheat, which were first identified in the rice genome. We also show that single-copy (SC) rice genes, those representing unique matches with wheat expressed sequence tag (EST) unigene contigs in the whole rice genome, show more than twice the proportion of genes mapping to syntenic wheat chromosome as compared to the multicopy (MC) or duplicated rice genes. While 58.7% of the 1,244 mapped SC rice genes were located in single syntenic wheat chromosome groups, the remaining 41.3% were distributed randomly to the other six non-syntenic wheat groups. This could only be explained by a background dispersal of genes in the genome through transposition or other unknown mechanism. The breakdown of rice–wheat synteny due to such transpositions was much greater near the wheat centromeres. Furthermore, the SC rice genes revealed a conserved primordial gene order that gives clues to the origin of rice and wheat chromosomes from a common ancestor through polyploidy, aneuploidy, centromeric fusions, and translocations. Apart from the bin-mapped wheat EST contigs, we also compared 56,298 predicted rice genes with 39,813 wheat EST contigs assembled from 409,765 EST sequences and identified 7,241 SC rice gene homologs of wheat. Based on the conserved colinearity of 1,063 mapped SC rice genes across the bins of individual wheat chromosomes, we predicted the wheat bin location of 6,178 unmapped SC rice gene homologs and validated the location of 213 of these in the telomeric bins of 21 wheat chromosomes with 35.4% initial success. This opens up the possibility of directed mapping of a large number of conserved SC rice gene homologs in wheat. Overall, only 46.4% of these SC genes code for proteins with known functional domains; the remaining 53.6% have unknown function, and hence, represent an important, but yet, under explored category of genes. Electronic supplementary material Supplementary material is available for this article at and is accessible for authorized users.