Comparative analysis of sucrose phosphate synthase (SPS) gene family between Saccharum officinarum and Saccharum spontaneum.

BACKGROUND: Sucrose phosphate synthase (SPS) genes play vital roles in sucrose production across various plant species. Modern sugarcane cultivar is derived from the hybridization between the high sugar content species Saccharum officinarum and the high stress tolerance species Saccharum spontaneum, generating one of the most complex genomes among all crops. The genomics of sugarcane SPS remains under-studied despite its profound impact on sugar yield. RESULTS: In the present study, 8 and 6 gene sequences for SPS were identified from the BAC libraries of S. officinarum and S. spontaneum, respectively. Phylogenetic analysis showed that SPSD was newly evolved in the lineage of Poaceae species with recently duplicated genes emerging from the SPSA clade. Molecular evolution analysis based on Ka/Ks ratios suggested that polyploidy reduced the selection pressure of SPS genes in Saccharum species. To explore the potential gene functions, the SPS expression patterns were analyzed based on RNA-seq and proteome dataset, and the sugar content was detected using metabolomics analysis. All the SPS members presented the trend of increasing expression in the sink-source transition along the developmental gradient of leaves, suggesting that the SPSs are involved in the photosynthesis in both Saccharum species as their function in dicots. Moreover, SPSs showed the higher expression in S. spontaneum and presented expressional preference between stem (SPSA) and leaf (SPSB) tissue, speculating they might be involved in the differentia of carbohydrate metabolism in these two Saccharum species, which required further verification from experiments. CONCLUSIONS: SPSA and SPSB genes presented relatively high expression and differential expression patterns between the two Saccharum species, indicating these two SPSs are important in the formation of regulatory networks and sucrose traits in the two Saccharum species. SPSB was suggested to be a major contributor to the sugar accumulation because it presented the highest expressional level and its expression positively correlated with sugar content. The recently duplicated SPSD2 presented divergent expression levels between the two Saccharum species and the relative protein content levels were highest in stem, supporting the neofunctionalization of the SPSD subfamily in Saccharum.

Comprehensively Characterizing the Cytological Features of Saccharum spontaneum by the Development of a Complete Set of Chromosome-Specific Oligo Probes.

Chromosome-specific identification is a powerful technique in the study of genome structure and evolution. However, there is no reliable cytogenetic marker to unambiguously identify each of the chromosomes in sugarcane (Saccharum spp., Poaceae), which has a complex genome with a high level of ploidy and heterozygosity. In this study, we developed a set of oligonucleotide (oligo)-based probes through bioinformatic design and massive synthetization. These probes produced a clear and bright single signal in each of the chromosomes and their eight homologous chromosomes in the ancient species Saccharum spontaneum (2n = 8x = 64). Thus, they can be used as reliable markers to robustly label each of the chromosomes in S. spontaneum. We then obtained the karyotype data and established a nomenclature based on chromosomal sizes for the eight chromosomes of the octoploid S. spontaneum. In addition, we also found that the 45S and 5S rDNAs demonstrated high copy number variations among different homologous chromosomes, indicating a rapid evolution of the highly repeated sequence after polyploidization. Our fluorescence in situ hybridization (FISH) assay also demonstrated that these probes could be used as cross-species markers between or within the genera of Sorghum and Saccharum. By comparing FISH analyses, we discovered that several chromosome rearrangement events occurred in S. spontaneum, which might have contributed to the basic chromosome number reduction from 10 in sorghum to 8 in sugarcane. Consistent identification of individual chromosomes makes molecular cytogenetic study possible in sugarcane and will facilitate fine chromosomal structure and karyotype evolution of the genus Saccharum.

Comparative Analysis of Homologous Sequences of Saccharum officinarum and Saccharum spontaneum Reveals Independent Polyploidization Events.

Sugarcane (Saccharum spp. hybrids) is an economically important crop widely grown in tropical and subtropical regions for sugar and ethanol production. However, the large genome size, high ploidy level, interspecific hybridization and aneuploidy make sugarcane one of the most complex genomes and have long hampered genome research in sugarcane. Modern sugarcane cultivars are derived from interspecific hybridization between S. officinarum and S. spontaneum with 80-90% of the genome from S. officinarum and 10-20% of the genome from S. spontaneum. We constructed bacterial artificial chromosome (BAC) libraries of S. officinarum variety LA Purple (2n = 8x = 80) and S. spontaneum haploid clone AP85-441 (2n = 4x = 32), and selected and sequenced 97 BAC clones from the two Saccharum BAC libraries. A total of 5,847,280 bp sequence from S. officinarum and 5,011,570 bp from S. spontaneum were assembled and 749 gene models were annotated in these BACs. A relatively higher gene density and lower repeat content were observed in S. spontaneum BACs than in S. officinarum BACs. Comparative analysis of syntenic regions revealed a high degree of collinearity in genic regions between Saccharum and Sorghum bicolor and between S. officinarum and S. spontaneum. In the syntenic regions, S. spontaneum showed expansion relative to S. officinarum, and both S. officinarum and S. spontaneum showed expansion relative to sorghum. Among the 75 full-length LTR retrotransposons identified in the Saccharum BACs, none of them are older than 2.6 mys and no full-length LTR elements are shared between S. officinarum and S. spontaneum. In addition, divergence time estimated using a LTR junction marker and a syntenic gene shared by 3 S. officinarum and 1 S. spontaneum BACs revealed that the S. spontaneum intergenic region was distant to those from the 3 homologous regions in S. officinarum. Our results suggested that S. officinarum and S. spontaneum experienced at least two rounds of independent polyploidization in each lineage after their divergence from a common ancestor.

Development and Applications of Chromosome-Specific Cytogenetic BAC-FISH Probes in S. spontaneum.

Saccharum spontaneum is a major Saccharum species that contributed to the origin of modern sugarcane cultivars, and due to a high degree of polyploidy is considered to be a plant species with one of the most complex genetics. Fluorescence in situ hybridization (FISH) is a powerful and widely used tool in genome studies. Here, we demonstrated that FISH based on bacterial artificial chromosome (BAC) clones can be used as a specific cytological marker to identify S. spontaneum individual chromosomes and study the relationship between S. spontaneum and other related species. We screened low-copy BACs as probes from the sequences of a high coverage of S. spontaneum BAC library based on BLAST search of the sorghum genome. In total, we isolated 49 positive BAC clones, and identified 27 BAC clones that can give specific signals on the S. spontaneum chromosomes. Of the 27 BAC probes, 18 were confirmed to be able to discriminate the eight basic chromosomes of S. spontaneum. Moreover, BAC-24, BAC-66, BAC-78, BAC-69, BAC-71, BAC-73, and BAC-77 probes were used to construct physical maps of chromosome 1 and chromosome 2 of S. spontaneum, which indicated synteny in Sb01 between S. spontaneum and sorghum. Furthermore, we found that BAC-14 and BAC-19 probes, corresponding to the sorghum chromosomes 2 and 8, respectively, localized to different arms of the same S. spontaneum chromosome, suggesting that there was an inter-chromosomal rearrangement event between S. spontaneum and sorghum. Our study provides the first set of chromosome-specific cytogenetic markers in Saccharum and is critical for future advances in cytogenetics and genome sequencing studies in Saccharum.

Papain-like cysteine proteases in Carica papaya: lineage-specific gene duplication and expansion.

BACKGROUND: Papain-like cysteine proteases (PLCPs), a large group of cysteine proteases structurally related to papain, play important roles in plant development, senescence, and defense responses. Papain, the first cysteine protease whose structure was determined by X-ray crystallography, plays a crucial role in protecting papaya from herbivorous insects. Except the four major PLCPs purified and characterized in papaya latex, the rest of the PLCPs in papaya genome are largely unknown. RESULTS: We identified 33 PLCP genes in papaya genome. Phylogenetic analysis clearly separated plant PLCP genes into nine subfamilies. PLCP genes are not equally distributed among the nine subfamilies and the number of PLCPs in each subfamily does not increase or decrease proportionally among the seven selected plant species. Papaya showed clear lineage-specific gene expansion in the subfamily III. Interestingly, all four major PLCPs purified from papaya latex, including papain, chymopapain, glycyl endopeptidase and caricain, were grouped into the lineage-specific expansion branch in the subfamily III. Mapping PLCP genes on chromosomes of five plant species revealed that lineage-specific expansions of PLCP genes were mostly derived from tandem duplications. We estimated divergence time of papaya PLCP genes of subfamily III. The major duplication events leading to lineage-specific expansion of papaya PLCP genes in subfamily III were estimated at 48 MYA, 34 MYA, and 16 MYA. The gene expression patterns of the papaya PLCP genes in different tissues were assessed by transcriptome sequencing and qRT-PCR. Most of the papaya PLCP genes of subfamily III expressed at high levels in leaf and green fruit tissues. CONCLUSIONS: Tandem duplications played the dominant role in affecting copy number of PLCPs in plants. Significant variations in size of the PLCP subfamilies among species may reflect genetic adaptation of plant species to different environments. The lineage-specific expansion of papaya PLCPs of subfamily III might have been promoted by the continuous reciprocal selective effects of herbivore attack and plant defense.

Evolution and expression of the fructokinase gene family in Saccharum.

BACKGROUND: Sugarcane is an important sugar crop contributing up to about 80% of the world sugar production. Efforts to characterize the genes involved in sugar metabolism at the molecular level are growing since increasing sugar content is a major goal in the breeding of new sugarcane varieties. Fructokinases (FRK) are the main fructose phosphorylating enzymes with high substrate specificity and affinity. RESULTS: In this study, by combining comparative genomics approaches with BAC resources, seven fructokinase genes were identified in S. spontaneum. Phylogenetic analysis based on representative monocotyledon and dicotyledon plant species suggested that the FRK gene family is ancient and its evolutionary history can be traced in duplicated descending order: SsFRK4, SsFRK6/SsFRK7,SsFRK5, SsFRK3 and SsFRK1/SsFRK2. Among the close orthologs, the number and position of exons in FRKs were conserved; in contrast, the size of introns varied among the paralogous FRKs in Saccharum. Genomic constraints were analyzed within the gene alleles and between S. spontaneum and Sorghum bicolor, and gene expression analysis was performed under drought stress and with exogenous applications of plant hormones. FRK1, which was under strong functional constraint selection, was conserved among the gene allelic haplotypes, and displayed dominant expression levels among the gene families in the control conditions, suggesting that FRK1 plays a major role in the phosphorylation of fructose. FRK3 and FRK5 were dramatically induced under drought stress, and FRK5 was also found to increase its expression levels in the mature stage of Saccharum. Similarly, FRK3 and FRK5 were induced in response to drought stress in Saccharum. FRK2 and FRK7 displayed lower expression levels than the other FRK family members; FRK2 was under strong genomic selection constraints whereas FRK7 was under neutral selection. FRK7 may have become functionally redundant in Saccharum through pseudogenization. FRK4 and FRK6 shared the most similar expression pattern: FRK4 was revealed to have higher expression levels in mature tissues than in premature tissues of Saccharum, and FRK6 presented a slight increase under drought stress. CONCLUSIONS: Our study presents a comprehensive genomic study of the entire FRK gene family in Saccharum, providing the foundations for approaches to characterize the molecular mechanism regulated by the SsFRK family in sugarcane.

Comparative structural analysis of Bru1 region homeologs in Saccharum spontaneum and S. officinarum.

BACKGROUND: Sugarcane is a major sugar and biofuel crop, but genomic research and molecular breeding have lagged behind other major crops due to the complexity of auto-allopolyploid genomes. Sugarcane cultivars are frequently aneuploid with chromosome number ranging from 100 to 130, consisting of 70-80 % S. officinarum, 10-20 % S. spontaneum, and 10 % recombinants between these two species. Analysis of a genomic region in the progenitor autoploid genomes of sugarcane hybrid cultivars will reveal the nature and divergence of homologous chromosomes. RESULTS: To investigate the origin and evolution of haplotypes in the Bru1 genomic regions in sugarcane cultivars, we identified two BAC clones from S. spontaneum and four from S. officinarum and compared to seven haplotype sequences from sugarcane hybrid R570. The results clarified the origin of seven homologous haplotypes in R570, four haplotypes originated from S. officinarum, two from S. spontaneum and one recombinant.. Retrotransposon insertions and sequences variations among the homologous haplotypes sequence divergence ranged from 18.2 % to 60.5 % with an average of 33.7 %. Gene content and gene structure were relatively well conserved among the homologous haplotypes. Exon splitting occurred in haplotypes of the hybrid genome but not in its progenitor genomes. Tajima's D analysis revealed that S. spontaneum hapotypes in the Bru1 genomic regions were under strong directional selection. Numerous inversions, deletions, insertions and translocations were found between haplotypes within each genome. CONCLUSIONS: This is the first comparison among haplotypes of a modern sugarcane hybrid and its two progenitors. Tajima's D results emphasized the crucial role of this fungal disease resistance gene for enhancing the fitness of this species and indicating that the brown rust resistance gene in R570 is from S. spontaneum. Species-specific InDel, sequences similarity and phylogenetic analysis of homologous genes can be used for identifying the origin of S. spontaneum and S. officinarum haplotype in Saccharum hybrids. Comparison of exon splitting among the homologous haplotypes suggested that the genome rearrangements in Saccharum hybrids after hybridization. The combined minimum difference at 19.5 % among homologous chromosomes in S. officinarum would be sufficient for proper genome assembly of this autopolyploid genome. Retrotransposon insertions and sequences variations among the homologous haplotypes sequence divergence may allow sequencing and assembling the autopolyploid Saccharum genomes and the auto-allopolyploid hybrid genomes using whole genome shotgun sequencing.

Structure, phylogeny, allelic haplotypes and expression of sucrose transporter gene families in Saccharum.

BACKGROUND: Sugarcane is an economically important crop contributing to about 80% of the world sugar production. Increasing efforts in molecular biological studies have been performed for improving the sugar yield and other relevant important agronomic traits. However, due to sugarcane's complicated genomes, it is still challenging to study the genetic basis of traits, such as sucrose accumulation. Sucrose transporters (SUTs) are critical for both phloem loading in source tissue and sucrose uptaking in sink tissue, and are considered to be the control points for regulating sucrose storage. However, no genomic study for sugarcane sucrose transporter (SsSUT) families has been reported up to date. RESULTS: By using comparative genomics and bacterial artificial chromosomes (BACs), six SUT genes were identified and characterized in S. spontaenum. Phylogenetic analyses revealed that the two pairs SsSUTs (SsSUT1/SsSUT3 and SsSUT5/SsSUT6) could be clustered together into two separate monocot specific SUT groups, while SsSUT2 and SsSUT4 were separated into the other two groups, with members from both dicot and monocot species. Gene structure comparison demonstrated that the number and position of exons/introns in SUTs were highly conserved among the close orthologs; in contrast, there were variations among the paralogous SUTs in Sacchuarm. Though with the high polyploidy level, gene allelic haplotype comparative analysis showed that the examined four SsSUT members exhibited conservations of gene structures and amino acid sequences among the allelic haplotypes accompanied by variations of intron sizes. Gene expression analyses were performed for tissues from seedlings under drought stress and mature plants of three Saccharum species (S.officinarnum, S.spotaneum and S.robustum). Both SUT1 and SUT4 expressed abundantly at different conditions. SUT2 had similar expression level in all of the examined tissues, but SUT3 was undetectable. Both of SUT5 and SUT6 had lower expression level than other gene member, and expressed stronger in source leaves and are likely to play roles in phloem loading. In the seeding plant leave under water stress, four genes SUT1, SUT2, SUT4 and SUT5 were detectable. In these detectable genes, SUT1 and SUT4 were down regulated, while, SUT2 and SUT5 were up regulated. CONCLUSIONS: In this study, we presented the first comprehensive genomic study for a whole gene family, the SUT family, in Saccharum. We speculated that there were six SUT members in the S. spotaneum genome. Out of the six members, SsSUTs, SsSUT5 and SsSUT6 were recent duplication genes accompanied by rapid evolution, while, SsSUT2 and SsSUT4 were the ancient members in the families. Despite the high polypoidy genome, functional redundancy may not exist among the SUTs allelic haplotypes supported by the evidence of strong purifying selection of the gene allele. SUT3 could be a low active member in the family because it is undetectable in our study, but it might not be a pseudogene because it harbored integrated gene structure. SUT1 and SUT4 were the main members for the sucrose transporter, while, these SUTs had sub-functional divergence in response to sucrose accumulation and plant development in Saccharum.

Cold Responsive Gene Expression Profiling of Sugarcane and Saccharum spontaneum with Functional Analysis of a Cold Inducible Saccharum Homolog of NOD26-Like Intrinsic Protein to Salt and Water Stress.

Transcriptome analysis of sugarcane hybrid CP72-1210 (cold susceptible) and Saccharum spontaneum TUS05-05 (cold tolerant) using Sugarcane Assembled Sequences (SAS) from SUCEST-FUN Database showed that a total of 35,340 and 34,698 SAS genes, respectively, were expressed before and after chilling stress. The analysis revealed that more than 600 genes are differentially expressed in each genotype after chilling stress. Blast2Go annotation revealed that the major difference in gene expression profiles between CP72-1210 and TUS05-05 after chilling stress are present in the genes related to the transmembrane transporter activity. To further investigate the relevance of transmembrane transporter activity against abiotic stress tolerance, a S. spontaneum homolog of a NOD26-like major intrinsic protein gene (SspNIP2) was selected for functional analysis, of which expression was induced after chilling stress in the cold tolerant TUS05-05. Quantitative real-time PCR showed that SspNIP2 expression was increased ~2.5 fold at 30 minutes after cold treatment and stayed induced throughout the 24 hours of cold treatment. The amino acid sequence analysis of the cloned SspNIP2 confirmed the presence of six transmembrane domains and two NPA (Asn-Pro-Ala) motifs, signature features of major intrinsic protein families. Amino acid analysis confirmed that four amino acids, comprising the ar/R (aromatic residue/arginine) region responsible for the substrate specificity among MIPs, are conserved among monocot silicon transporters and SspNIP2. Salinity stress test on SspNIP2 transgenic tobacco plants resulted in more vigorous transgenic lines than the non-transgenic tobacco plants, suggesting some degree of tolerance to salt stress conferred by SspNIP2. SspNIP2-transgenic plants, exposed to 2 weeks of water stress without irrigation, developed various degrees of water stress symptom. The water stress test confirmed that the SspNIP2 transgenic lines had lower evapotranspiration rates than non-transgenic lines, suggesting that SspNIP2 transgenic lines showed a slight tolerance to the early water stress compared to wild type plants.

Micro-collinearity and genome evolution in the vicinity of an ethylene receptor gene of cultivated diploid and allotetraploid coffee species (Coffea).

Arabica coffee (Coffea arabica L.) is a self-compatible perennial allotetraploid species (2n=4x=44), whereas Robusta coffee (C. canephora L.) is a self-incompatible perennial diploid species (2n=2x=22). C. arabica (C(a) C(a) E(a) E(a) ) is derived from a spontaneous hybridization between two closely related diploid coffee species, C. canephora (CC) and C. eugenioides (EE). To investigate the patterns and degree of DNA sequence divergence between the Arabica and Robusta coffee genomes, we identified orthologous bacterial artificial chromosomes (BACs) from C. arabica and C. canephora, and compared their sequences to trace their evolutionary history. Although a high level of sequence similarity was found between BACs from C. arabica and C. canephora, numerous chromosomal rearrangements were detected, including inversions, deletions and insertions. DNA sequence identity between C. arabica and C. canephora orthologous BACs ranged from 93.4% (between E(a) and C(a) ) to 94.6% (between C(a) and C). Analysis of eight orthologous gene pairs resulted in estimated ages of divergence between 0.046 and 0.665 million years, indicating a recent origin of the allotetraploid species C. arabica. Analysis of transposable elements revealed differential insertion events that contributed to the size increase in the C(a) sub-genome compared to its diploid relative. In particular, we showed that insertion of a Ty1-copia LTR retrotransposon occurred specifically in C. arabica, probably shortly after allopolyploid formation. The two sub-genomes of C. arabica, C(a) and E(a) , showed sufficient sequence differences, and a whole-genome shotgun approach could be suitable for sequencing the allotetraploid genome of C. arabica.

Isolating promoters of multigene family members from the polyploid sugarcane genome by PCR-based walking in BAC DNA.

The availability of a wider range of promoters for regulated expression in valuable transgenic crops would benefit functional genomics studies and current biotechnology programs aimed at improved productivity. Polymerase chain reaction (PCR)-based genome walking techniques are commonly used to isolate promoters or 5' flanking genomic regions adjacent to known cDNA sequences in genomes that are not yet completely sequenced. However, these techniques are problematic when applied directly to DNA isolated from crops with highly complex and large genomes. An adaptor ligation-mediated PCR-based BAC genome walking method is described here for the efficient isolation of promoters of multigene family members, such as the putative defense and fiber biosynthesis DIRIGENT genes that are abundant in the stem of sugarcane, a species with a highly polyploid genome. The advantage of this method is the efficient and specific amplification of the target promoter using BAC genomic DNA as template for the adaptor ligation-mediated PCR walking.

Microcollinearity between autopolyploid sugarcane and diploid sorghum genomes.

BACKGROUND: Sugarcane (Saccharum spp.) has become an increasingly important crop for its leading role in biofuel production. The high sugar content species S. officinarum is an octoploid without known diploid or tetraploid progenitors. Commercial sugarcane cultivars are hybrids between S. officinarum and wild species S. spontaneum with ploidy at approximately 12x. The complex autopolyploid sugarcane genome has not been characterized at the DNA sequence level. RESULTS: The microsynteny between sugarcane and sorghum was assessed by comparing 454 pyrosequences of 20 sugarcane bacterial artificial chromosomes (BACs) with sorghum sequences. These 20 BACs were selected by hybridization of 1961 single copy sorghum overgo probes to the sugarcane BAC library with one sugarcane BAC corresponding to each of the 20 sorghum chromosome arms. The genic regions of the sugarcane BACs shared an average of 95.2% sequence identity with sorghum, and the sorghum genome was used as a template to order sequence contigs covering 78.2% of the 20 BAC sequences. About 53.1% of the sugarcane BAC sequences are aligned with sorghum sequence. The unaligned regions contain non-coding and repetitive sequences. Within the aligned sequences, 209 genes were annotated in sugarcane and 202 in sorghum. Seventeen genes appeared to be sugarcane-specific and all validated by sugarcane ESTs, while 12 appeared sorghum-specific but only one validated by sorghum ESTs. Twelve of the 17 sugarcane-specific genes have no match in the non-redundant protein database in GenBank, perhaps encoding proteins for sugarcane-specific processes. The sorghum orthologous regions appeared to have expanded relative to sugarcane, mostly by the increase of retrotransposons. CONCLUSIONS: The sugarcane and sorghum genomes are mostly collinear in the genic regions, and the sorghum genome can be used as a template for assembling much of the genic DNA of the autopolyploid sugarcane genome. The comparable gene density between sugarcane BACs and corresponding sorghum sequences defied the notion that polyploidy species might have faster pace of gene loss due to the redundancy of multiple alleles at each locus.