Effective identification of varieties by nucleotide polymorphisms and its application for essentially derived variety identification in rice.

BACKGROUND: Plant variety identification is the one most important of agricultural systems. Development of DNA marker profiles of released varieties to compare with candidate variety or future variety is required. However, strictly speaking, scientists did not use most existing variety identification techniques for "identification" but for "distinction of a limited number of cultivars," of which generalization ability always not be well estimated. Because many varieties have similar genetic backgrounds, even some essentially derived varieties (EDVs) are involved, which brings difficulties for identification and breeding progress. A fast, accurate variety identification method, which also has good performance on EDV determination, needs to be developed. RESULTS: In this study, with the strategy of "Divide and Conquer," a variety identification method Conditional Random Selection (CRS) method based on SNP of the whole genome of 3024 rice varieties was developed and be applied in essentially derived variety (EDV) identification of rice. CRS is a fast, efficient, and automated variety identification method. Meanwhile, in practical, with the optimal threshold of identity score searched in this study, the set of SNP (including 390 SNPs) showed optimal performance on EDV and non-EDV identification in two independent testing datasets. CONCLUSION: This approach first selected a minimal set of SNPs to discriminate non-EDVs in the 3000 Rice Genome Project, then united several simplified SNP sets to improve its generalization ability for EDV and non-EDV identification in testing datasets. The results suggested that the CRS method outperformed traditional feature selection methods. Furthermore, it provides a new way to screen out core SNP loci from the whole genome for DNA fingerprinting of crop varieties and be useful for crop breeding.

Identification and Allele Combination Analysis of Rice Grain Shape-Related Genes by Genome-Wide Association Study.

Grain shape is an important agronomic character of rice, which affects the appearance, processing, and the edible quality. Screening and identifying more new genes associated with grain shape is beneficial to further understanding the genetic basis of grain shape and provides more gene resources for genetic breeding. This study has a natural population containing 623 indica rice cultivars. Genome-wide association studies/GWAS of several traits related to grain shape (grain length/GL, grain width/GW, grain length to width ratio/GLWR, grain circumferences/GC, and grain size/grain area/GS) were conducted by combining phenotypic data from four environments and the second-generation resequencing data, which have identified 39 important Quantitative trait locus/QTLs. We analyzed the 39 QTLs using three methods: gene-based association analysis, haplotype analysis, and functional annotation and identified three cloned genes (GS3, GW5, OsDER1) and seven new candidate genes in the candidate interval. At the same time, to effectively utilize the genes in the grain shape-related gene bank, we have also analyzed the allelic combinations of the three cloned genes. Finally, the extreme allele combination corresponding to each trait was found through statistical analysis. This study's novel candidate genes and allele combinations will provide a valuable reference for future breeding work.

Molecular cloning and characterization of two manganese superoxide dismutases from Miscanthus × giganteus.

KEY MESSAGE: Six MnSOD genes were isolated from five Miscanthus species. MgMnSOD1 functions in mitochondria and MgMnSOD1 seems to be the main MnSOD gene involved in stress response of M. × giganteus. Miscanthus × giganteus is a promising biomass energy crop with advantages of vigorous growth, high yield, low fertilizer and pesticide inputs. However, poor overwinter ability limits its widespread cultivation. Moreover, narrow genetic base may increase the risk of susceptibility to diseases and pests. Manganese superoxide dismutase (MnSOD), an important antioxidant enzyme involved in stress tolerance is able to protect plant cells from accumulated reactive oxygen species by converting superoxide to peroxide and oxygen. In many plants, overexpression of MnSOD has shown the ability to enhance the resistance to various stresses. This article describes the studies performed in an attempt to elucidate the molecular and enzymatic properties of MnSODs in M. × giganteus. MnSOD genes from M. × giganteus (MgMnSOD1, MgMnSOD2), M. lutarioriparia (MlMnSOD), M. sacchariflora (MsaMnSOD), M. sinensis (MsiMnSOD), and M. floridulus (MfMnSOD) were cloned and sequenced. The sequence analysis and expression patterns of MgMnSOD1 and MgMnSOD2 suggest that they were orthologous genes which were inherited from the two parents, M. sacchariflora and M. sinensis, respectively. In addition, MgMnSOD1 is predicted to be the main MnSOD gene involved in stress response of M. × giganteus. The activity of purified recombinant MgMnSOD1 was 1854.79 ± 39.98 U mg(-1) (mean ± SD). Further enzymatic assays revealed that the protein exhibited an outstanding thermal stability. MgMnSOD1 is predicted to be targeted to mitochondria and involved in removing the superoxide radical generated by respiration. The presence and sequences of other SOD isozymes transcripts were also investigated in this study.

Cloning a glutathione peroxidase gene from Nelumbo nucifera and enhanced salt tolerance by overexpressing in rice.

A full-length cDNA clone encoding an 866 bp-length glutathione peroxidase protein (NnGPX) was isolated from lotus (Nelumbo nucifera L.). The deduced amino acid sequence of the NnGPX gene had significant homology with ATGPX6. A 3D structural model of the NnGPX was constructed by homology modeling. The cloned NnGPX gene was expressed in Escherichia coli, and a fusion protein of about 40 kDa was detected after isopropyl thiogalactoside induction. Under different concentrations of Na2SeO3 treatments, NnGPX was found to be an enzyme that does not contain selenium. Real-time PCR analysis showed that the NnGPX gene was expressed in all organs of lotus, and its high expression mainly occurred in organs with active metabolisms. NnGPX transcript increased remarkably in response to cold, heat, mechanical damage, and salt treatment. Subsequently, the NnGPX gene was introduced in Oryza sativa cv. Yuetai B. PCR results verified the integration of this gene into the genome of rice and reverse transcription-PCR verified that this gene had been expressed in transgenic rice. The transgenic plants were significantly more tolerant to salt stress compared with the wild-type.

Development of microsatellite markers by transcriptome sequencing in two species of Amorphophallus (Araceae).

BACKGROUND: Amorphophallus is a genus of perennial plants widely distributed in the tropics or subtropics of West Africa and South Asia. Its corms contain a high level of water-soluble glucomannan; therefore, it has long been used as a medicinal herb and food source. Genetic studies of Amorphophallus have been hindered by a lack of genetic markers. A large number of molecular markers are required for genetic diversity study and improving disease resistance in Amorphophallus. Here, we report large scale of transcriptome sequencing of two species: Amorphophallus konjac and Amorphophallus bulbifer using deep sequencing technology, and microsatellite (SSR) markers were identified based on these transcriptome sequences. RESULTS: cDNAs of A. konjac and A. bulbifer were sequenced using Illumina HiSeq™ 2000 sequencing technology. A total of 135,822 non-redundant unigenes were assembled from about 9.66 gigabases, and 19,596 SSRs were identified in 16,027 non-redundant unigenes. Di-nucleotide SSRs were the most abundant motif (61.6%), followed by tri- (30.3%), tetra- (5.6%), penta- (1.5%), and hexa-nucleotides (1%) repeats. The top di- and tri-nucleotide repeat motifs included AG/CT (45.2%) and AGG/CCT (7.1%), respectively. A total of 10,754 primer pairs were designed for marker development. Of these, 320 primers were synthesized and used for validation of amplification and assessment of polymorphisms in 25 individual plants. The total of 275 primer pairs yielded PCR amplification products, of which 205 were polymorphic. The number of alleles ranged from 2 to 14 and the polymorphism information content valued ranged from 0.10 to 0.90. Genetic diversity analysis was done using 177 highly polymorphic SSR markers. A phenogram based on Jaccard's similarity coefficients was constructed, which showed a distinct cluster of 25 Amorphophallus individuals. CONCLUSION: A total of 10,754 SSR markers have been identified in Amorphophallus using transcriptome sequencing. One hundred and seventy-seven polymorphic markers were successfully validated in 25 individuals. The large number of genetic markers developed in the present study should contribute greatly to research into genetic diversity and germplasm characterization in Amorphophallus.

Joint analysis of phenotype-effect-generation identifies loci associated with grain quality traits in rice hybrids.

Genetic improvement of grain quality is more challenging in hybrid rice than in inbred rice due to additional nonadditive effects such as dominance. Here, we describe a pipeline developed for joint analysis of phenotypes, effects, and generations (JPEG). As a demonstration, we analyze 12 grain quality traits of 113 inbred lines (male parents), five tester lines (female parents), and 565 (113×5) of their hybrids. We sequence the parents for single nucleotide polymorphisms calling and infer the genotypes of the hybrids. Genome-wide association studies with JPEG identify 128 loci associated with at least one of the 12 traits, including 44, 97, and 13 loci with additive effects, dominant effects, and both additive and dominant effects, respectively. These loci together explain more than 30% of the genetic variation in hybrid performance for each of the traits. The JEPG statistical pipeline can help to identify superior crosses for breeding rice hybrids with improved grain quality.

Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans).

Lotus (Nelumbo Adans.), a relict plant, is the testimony of long-term sustained ecological success, but the underlying genetic changes related to its survival strategy remains unclear. Here, we assembled the high-quality lotus genome, investigated genome variation of lotus mutation accumulation (MA) lines and reconstructed the demographic history of wild Asian lotus, respectively. We identified and validated 43 base substitutions fixed in MA lines, implying a spontaneous mutation rate of 1.4 × 10-9 base/generation in lotus shoot stem cells. The past history of lotus revealed that the ancestors of lotus in eastern and southern Asia could be traced back ~20 million years ago (Mya) and experienced twice significant bottlenecks and population splits. We further identified the selected genes among three lotus groups in different habitats, suggesting that 453 genes between tropical and temperate group and 410 genes between two subgroups from Northeastern China and the Yangtze River - Yellow River Basin might play important roles in natural selection in lotus's adaptation and resilience. Our findings not only improve an understanding of the lotus evolutionary history and the genetic basis of its survival advantages, but also provide valuable data for addressing various questions in evolution and protection for the relict plants.

Cloning and expression of one chloroplastic ascorbate peroxidase gene from Nelumbo nucifera.

A novel ascorbate peroxidase (APX) cDNA was obtained from Nelumbo nucifera (Elian). The phylogenetic analysis indicated that N. nucifera APX grouped together with chloroplastic APX of high plants. The recombinant protein expressed by PET-30a vector showed APX activity (0.04 mM ascorbate min(-1) mg(-1) protein). The APX mRNA was expressed in young leaves, roots, terminal buds, and leafstalks. Synergistic expression of N. nucifera APX and MnSOD mRNA was indicated in the short-term response to mechanical wounding.

Genome-wide association study and Mendelian randomization analysis provide insights for improving rice yield potential.

Rice yield per plant has a complex genetic architecture, which is mainly determined by its three component traits: the number of grains per panicle (GPP), kilo-grain weight (KGW), and tillers per plant (TP). Exploring ideotype breeding based on selection for genetically less complex component traits is an alternative route for further improving rice production. To understand the genetic basis of the relationship between rice yield and component traits, we investigated the four traits of two rice hybrid populations (575 + 1495 F1) in different environments and conducted meta-analyses of genome-wide association study (meta-GWAS). In total, 3589 significant loci for three components traits were detected, while only 3 loci for yield were detected. It indicated that rice yield is mainly controlled by minor-effect loci and hardly to be identified. Selecting quantitative trait locus/gene affected component traits to further enhance yield is recommended. Mendelian randomization design is adopted to investigate the genetic effects of loci on yield through component traits and estimate the genetic relationship between rice yield and its component traits by these loci. The loci for GPP or TP mainly had a positive genetic effect on yield, but the loci for KGW with different direction effects (positive effect or negative effect). Additionally, TP (Beta = 1.865) has a greater effect on yield than KGW (Beta = 1.016) and GPP (Beta = 0.086). Five significant loci for component traits that had an indirect effect on yield were identified. Pyramiding superior alleles of the five loci revealed improved yield. A combination of direct and indirect effects may better contribute to the yield potential of rice. Our findings provided a rationale for using component traits as indirect indices to enhanced rice yield, which will be helpful for further understanding the genetic basis of yield and provide valuable information for improving rice yield potential.

Genome-Wide Association Study Identified Novel Candidate Loci/Genes Affecting Lodging Resistance in Rice.

Lodging reduces rice yield, but increasing lodging resistance (LR) usually limits yield potential. Stem strength and leaf type are major traits related to LR and yield, respectively. Hence, understanding the genetic basis of stem strength and leaf type is of help to reduce lodging and increase yield in LR breeding. Here, we carried out an association analysis to identify quantitative trait locus (QTLs) affecting stem strength-related traits (internode length/IL, stem wall thickness/SWT, stem outer diameter/SOD, and stem inner diameter/SID) and leaf type-associated traits (Flag leaf length/FLL, Flag leaf angle/FLA, Flag leaf width/FLW, leaf-rolling/LFR and SPAD/Soil, and plant analyzer development) using a diverse panel of 550 accessions and evaluated over two years. Genome-wide association study (GWAS) using 4,076,837 high-quality single-nucleotide polymorphisms (SNPs) identified 89 QTLs for the nine traits. Next, through "gene-based association analysis, haplotype analysis, and functional annotation", the scope was narrowed down step by step. Finally, we identified 21 candidate genes in 9 important QTLs that included four reported genes (TUT1, OsCCC1, CFL1, and ACL-D), and seventeen novel candidate genes. Introgression of alleles, which are beneficial for both stem strength and leaf type, or pyramiding stem strength alleles and leaf type alleles, can be employed for LR breeding. All in all, the experimental data and the identified candidate genes in this study provide a useful reference for the genetic improvement of rice LR.

Development and identification of three functional markers associated with starch content in lotus (Nelumbo nucifera).

It have been significantly demonstrated that Hexokinase (HXK), Granule-bound starch synthase (GBSS) and ADP-glucose pyrophosphorylase (AGPase) are three critical enzymes in the starch biosynthetic pathway and are related to starch (amylose, amylopectin and total starch) content in lotus. It is important to develop functional markers in marker-assisted selection of lotus breeding. So far there have been few reports about lotus functional markers. In this study, based on insertion-deletions (INDELs) and single-nucleotide polymorphisms (SNPs), we developed three functional markers, FMHXK-E1, FMGBSS-I8 and FMAGPL-I1. FMHXK-E1 was developed based on polymorphisms of two haplotypes of NnHXK. 26 lotus cultivars that the 320-bp fragment presented in NnHXK had a lower content of amylose and a higher content of amylopectin. FMGBSS-I8 was developed based on polymorphisms of two haplotypes of NnGBSS. The group containing 32 lotus cultivars with the 210-bp fragment had less amylose content and more amylopectin content. FMAGPL-I1 was developed based on polymorphisms of two haplotypes of NnAGPL (ADP-glucose pyrophosphorylase large subunit gene). The group containing 40 lotus cultivars with the 362-bp fragment had less amylopectin, total starch content and more amylose content. According to the study, FMHXK-E1, FMGBSS-I8 and FMAGPL-I1 are closely related to lotus starch content. It could be provided research basis for molecular assisted selection of lotus starch content improve breeding efficiency.

Genetic, transcriptional, and regulatory landscape of monolignol biosynthesis pathway in Miscanthus × giganteus.

BACKGROUND: Miscanthus × giganteus is widely recognized as a promising lignocellulosic biomass crop due to its advantages of high biomass production, low environmental impacts, and the potential to be cultivated on marginal land. However, the high costs of bioethanol production still limit the current commercialization of lignocellulosic bioethanol. The lignin in the cell wall and its by-products released in the pretreatment step is the main component inhibiting the enzymatic reactions in the saccharification and fermentation processes. Hence, genetic modification of the genes involved in lignin biosynthesis could be a feasible strategy to overcome this barrier by manipulating the lignin content and composition of M. × giganteus. For this purpose, the essential knowledge of these genes and understanding the underlying regulatory mechanisms in M. × giganteus is required. RESULTS: In this study, MgPAL1, MgPAL5, Mg4CL1, Mg4CL3, MgHCT1, MgHCT2, MgC3'H1, MgCCoAOMT1, MgCCoAOMT3, MgCCR1, MgCCR2, MgF5H, MgCOMT, and MgCAD were identified as the major monolignol biosynthetic genes in M. × giganteus based on genetic and transcriptional evidence. Among them, 12 genes were cloned and sequenced. By combining transcription factor binding site prediction and expression correlation analysis, MYB46, MYB61, MYB63, WRKY24, WRKY35, WRKY12, ERF021, ERF058, and ERF017 were inferred to regulate the expression of these genes directly. On the basis of these results, an integrated model was summarized to depict the monolignol biosynthesis pathway and the underlying regulatory mechanism in M. × giganteus. CONCLUSIONS: This study provides a list of potential gene targets for genetic improvement of lignocellulosic biomass quality of M. × giganteus, and reveals the genetic, transcriptional, and regulatory landscape of the monolignol biosynthesis pathway in M. × giganteus.

Genetic Diversity and DNA Fingerprints of Three Important Aquatic Vegetables by EST-SSR Markers.

Twenty-two sacred lotus (Nelumbo nucifera), 46 taros (Colocasia esculenta) and 10 arrowheads (Sagittaria trifolia) were used as materials and combined with EST-SSR (expressed sequence tag-simple sequence repeats) primers developed by our laboratory. Core primers were screened from a large number of primers that were able to distinguish all materials with a high frequency of polymorphisms. Six pairs, twenty pairs and three pairs of core primers were screened from sacred lotus, taro, and arrowhead, respectively. The SSR fingerprints of these three important aquatic vegetables, producing 17-, 87- and 14-bit binary molecular identity cards, respectively, were separately determined by using the core primers. Since there were few core primers of sacred lotus and arrowhead, 3 and 9 primer pairs with higher polymorphic information content (PIC), respectively, were selected as candidate primers. These core and candidate primers were used to identify the purities of No.36 space lotus, Shandong 8502 taro and Wuhan arrowhead, which were 93.3% (84/90), 98.9% (89/90) and 100.0% (90/90), respectively. The fingerprints, displayed as binary molecular identification cards of three important aquatic vegetables, were obtained, and their purity was successfully determined with EST-SSR labeling technology. Phylogenetic trees were also constructed to analyze the genetic diversity of 22 sacred lotus, 46 taros and 10 arrowheads. This study classifies and identifies germplasm resources and is an important reference to test the authenticity and variety purity of other aquatic vegetables in the future.

Genomic selection of agronomic traits in hybrid rice using an NCII population.

BACKGROUND: Hybrid breeding is an effective tool to improve yield in rice, while parental selection remains the key and difficult issue. Genomic selection (GS) provides opportunities to predict the performance of hybrids before phenotypes are measured. However, the application of GS is influenced by several genetic and statistical factors. Here, we used a rice North Carolina II (NC II) population constructed by crossing 115 rice varieties with five male sterile lines as a model to evaluate effects of statistical methods, heritability, marker density and training population size on prediction for hybrid performance. RESULTS: From the comparison of six GS methods, we found that predictabilities for different methods are significantly different, with genomic best linear unbiased prediction (GBLUP) and least absolute shrinkage and selection operation (LASSO) being the best, support vector machine (SVM) and partial least square (PLS) being the worst. The marker density has lower influence on predicting rice hybrid performance compared with the size of training population. Additionally, we used the 575 (115 × 5) hybrid rice as a training population to predict eight agronomic traits of all hybrids derived from 120 (115 + 5) rice varieties each mating with 3023 rice accessions from the 3000 rice genomes project (3 K RGP). Of the 362,760 potential hybrids, selection of the top 100 predicted hybrids would lead to 35.5%, 23.25%, 30.21%, 42.87%, 61.80%, 75.83%, 19.24% and 36.12% increase in grain yield per plant, thousand-grain weight, panicle number per plant, plant height, secondary branch number, grain number per panicle, panicle length and primary branch number, respectively. CONCLUSIONS: This study evaluated the factors affecting predictabilities for hybrid prediction and demonstrated the implementation of GS to predict hybrid performance of rice. Our results suggest that GS could enable the rapid selection of superior hybrids, thus increasing the efficiency of rice hybrid breeding.

Integration of transcriptome and proteome analyses reveal molecular mechanisms for formation of replant disease in Nelumbo nucifera.

The normal growth of Nelumbo nucifera, a widely planted aquatic crop in Asia, was severely ruined by replant disease. The mechanism of replant disease was still unknown in aquatic crops. Complementary transcriptomic and proteomic analyses were performed by comparing seedings of first-year planting (FP) and consecutive planting (CP). 9810 differentially expressed genes (DEGs) were identified between FP and CP. Additionally, 975 differentially expressed proteins (DEPs) were obtained. The correlation of proteome and transcriptome illustrated phenylpropanoid biosynthesis, flavonoid biosynthesis, metabolic pathways, and MAPK signaling pathways were significantly activated. Peroxidase, determined as one of the key proteins in replant disease of N. nucifera, was phylogenetically analyzed. A new depiction of the molecular mechanism causing replant disease in N. nucifera was illustrated. A consecutive monoculture stimulated the generation of reactive oxygen species (ROS) and ethylene, altered the metabolic balance of lignin and flavonoid, and attenuated the activity of antioxidant enzymes through DNA methylation. Therefore, the accumulation of autotoxic allelochemicals and the deficiency of antioxidant enzymes unavoidably suppressed the normal growth and development of replanted N. nucifera.

Transcriptomics and proteomics reveal genetic and biological basis of superior biomass crop Miscanthus.

Miscanthus is a rhizomatous C4 grass which is considered as potential high-yielding energy crop with the low-nutrient requirements, high water-use efficiency, and capability of C mitigation. To better understand the genetic basis, an integrative analysis of the transcriptome and proteome was performed to identify important genes and pathways involved in Miscanthus leaves. At the transcript level, 64,663 transcripts in M. lutarioriparius, 97,043 in M. sacchariflorus, 97,043 in M. sinensis, 67,323 in M. floridulus and 70,021 in M. × giganteus were detected by an RNA sequencing approach. At the protein level, 1964 peptide-represented proteins were identified and 1933 proteins differed by 1.5-fold or more in their relative abundance, as indicated by iTRAQ (isobaric tags for relative and absolute quantitation) analysis. Phylogenies were constructed from the nearly taxa of Miscanthus. A large number of genes closely related to biomass production were found. And SSR markers and their corresponding primers were derived from Miscanthus transcripts and 90% of them were successfully detected by PCR amplification among Miacanthus species. These similarities and variations on the transcriptional and proteomic level between Miscanthus species will serve as a resource for research in Miscanthus and other lignocellulose crops.

Pedigree-based genome re-sequencing reveals genetic variation patterns of elite backbone varieties during modern rice improvement.

Rice breeding has achieved great productivity improvements by semi-dwarf varieties and hybrid vigour. Due to poor understanding of genetic basis of elite backbone varieties, the continuous increasing in rice yield still faces great challenges. Here, 52 elite rice varieties from three historical representative pedigrees were re-sequenced with 10.1× depth on average, and ~6.5 million single nucleotide polymorphisms (SNPs) were obtained. We identified thousands of low-diversity genomic regions and 0-diversity genes during breeding. Using pedigree information, we also traced SNP transmission patterns and observed breeding signatures in pedigree. These regions included the larger number of key well-known functional genes. Besides, 35 regions spanning 0.16% of the rice gnome had been differentially selected between conventional and restorer pedigrees. These genes identified here will be useful to the further pedigree breeding. Our study provides insights into the genetic basis of backbone varieties and will have immediate implications for performing genome-wide breeding by design.

Molecular Cloning and Expression Analysis of a Catalase Gene (NnCAT) from Nelumbo nucifera.

Rapid amplification cDNA end (RACE) assay was established to achieve the complete cDNA sequence of a catalase gene (NnCAT) from Nelumbo nucifera. The obtained full-length cDNA was 1666 bp in size and contained a 1476-bp open reading frame. The 3D structural model of NnCAT was constructed by homology modeling. The putative NnCAT possessed all the main characteristic amino acid residues and motifs of catalase (CAT) protein family, and the phylogenetic analysis revealed that NnCAT grouped together with high plants. Moreover, recombinant NnCAT showed the CAT activity (758 U/mg) at room temperature, holding high activity during temperature range of 20-50 °C, then the optimal pH of recombinant protein was assessed from pH 4 to pH 11. Additionally, real-time PCR assay demonstrated that NnCAT mRNA was expressed in various tissues of N. nucifera, with the highest expression in young leaf and lowest level in the root, and mRNA level of NnCAT was significantly augmented in response to short-time mechanical wounding. Different expression pattern of NnCAT gene suggested that NnCAT probably played a defensive role in the initial stages of oxidative stress, regulating the level of reactive oxygen species (ROS) by extracellular stimuli such as short-time mechanical wounding.

De novo transcriptome and small RNA analyses of two amorphophallus species.

Konjac is one of the most important glucomannan crops worldwide. The breeding and genomic researches are largely limited by the genetic basis of Amorphophallus. In this study, the transcriptomes of A. konjac and A. bulbifer were constructed using a high-throughput Illumina sequencing platform. All 108,651 unigenes with average lengths of 430 nt in A. konjac and 119,678 unigenes with average lengths of 439 nt were generated from 54,986,020 reads and 52,334,098 reads after filtering and assembly, respectively. A total of 54,453 transcripts in A. konjac and 55,525 in A. bulbifier were annotated by comparison with Nr, Swiss-Prot, KEGG, and COG databases after removing exogenous contaminated sequences. A total of 80,332 transcripts differentially expressed between A. konjac and A. bulbifer. The majority of the genes that are associated with konjac glucomannan biosynthetic pathway were identified. Besides, the small RNAs in A. konjac leaves were also obtained by deep sequencing technology. All of 5,499,903 sequences of small RNAs were obtained with the length range between 18 and 30 nt. The potential targets for the miRNAs were also predicted according to the konjac transcripts. Our study provides a systematic overview of the konjac glucomannan biosynthesis genes that are involved in konjac leaves and should facilitate further understanding of the crucial roles of carbohydrate synthesis and other important metabolism pathways in Amorphophallus.