Zhongli Hu

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

BackgroundAmorphophallus 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.ResultscDNAs 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.ConclusionA 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.

QTL Mapping of Combining Ability and Heterosis of Agronomic Traits in Rice Backcross Recombinant Inbred Lines and Hybrid Crosses

Background Combining ability effects are very effective genetic parameters in deciding the next phase of breeding programs. Although some breeding strategies on the basis of evaluating combining ability have been utilized extensively in hybrid breeding, little is known about the genetic basis of combining ability. Combining ability is a complex trait that is controlled by polygenes. With the advent and development of molecular markers, it is feasible to evaluate the genetic bases of combining ability and heterosis of elite rice hybrids through QTL analysis. Methodology/Principal Findings In the present study, we first developed a QTL-mapping method for dissecting combining ability and heterosis of agronomic traits. With three testcross populations and a BCRIL population in rice, biometric and QTL analyses were conducted for ten agronomic traits. The significance of general combining ability and special combining ability for most of the traits indicated the importance of both additive and non-additive effects on expression levels. A large number of additive effect QTLs associated with performance per se of BCRIL and general combining ability, and dominant effect QTLs associated with special combining ability and heterosis were identified for the ten traits. Conclusions/Significance The combining ability of agronomic traits could be analyzed by the QTL mapping method. The characteristics revealed by the QTLs for combining ability of agronomic traits were similar with those by multitudinous QTLs for agronomic traits with performance per se of BCRIL. Several QTLs (1–6 in this study) were identified for each trait for combining ability. It demonstrated that some of the QTLs were pleiotropic or linked tightly with each other. The identification of QTLs responsible for combining ability and heterosis in the present study provides valuable information for dissecting genetic basis of combining ability.

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.

Development and characterisation of EST-SSR markers by transcriptome sequencing in taro ( Colocasia esculenta (L.) Schoot)

Taro (Colocasia esculenta) is an important crop with a long history of cultivation. In this study 5278 SSRs were identified in taro transcriptome data. A total of 2858 primer pairs were designed for marker development. 100 primers were randomly selected and synthesized. Among them, 72 primer pairs were successfully amplified and 62 were polymorphic in taro accessions. The number of alleles ranged from 2 to 14 for each different polymorphic locus and the polymorphism information content valued ranged from 0.01 to 0.82. The phylogenetic tree was also constructed to analyse the genetic diversity in 68 taro accessions. The large number of taro SSR markers developed in the present study will be useful in the researches of genetic diversity, germplasm characterization and molecular breeding etc.

Genomewide mapping reveals a combination of different genetic effects causing the genetic basis of heterosis in two elite rice hybrids

North Carolina design III (NCIII) is one of the most powerful and widely used mating designs for understanding the genetic basis of heterosis. However, the quantitative trait mapping (QTL) conducted in previous studies with this design was mainly based on analysis of variance (ANOVA), composite interval or multiple interval mapping methods. These methodologies could not investigate all kinds of genetic effects, especially epistatic effects, simultaneously on the whole genome. In this study, with a statistical method for mapping epistatic QTL associated with heterosis using the recombinant inbred line (RIL)-based NCIII design, we conducted QTL mapping for nine agronomic traits of two elite hybrids to characterize the mode of gene action contributing to heterosis on a whole genomewide scale. In total, 23 main-effect QTL (M-QTL) and 23 digenic interactions in IJ (indica ×japonica) hybrids, 11 M-QTL and 82 digenic interactions in II (indica ×indica) hybrid QTLs were identified in the present study. The variation explained by individual M-QTL or interactions ranged from 2.3 to 11.0%. The number of digenic interactions and the total variation explained by interactions of each trait were larger than those of M-QTL. The augmented genetic effect ratio of most M-QTL and digenic interactions in (L 1–L 2) data of two backcross populations (L 1 and L 2) showed complete dominance or overdominance, and in (L 1+L2) data showed an additive effect. Our results indicated that the dominance, overdominance and epistatic effect were important in conditioning the genetic basis of heterosis of the two elite hybrids. The relative contributions of the genetic components varied with traits and the genetic basis of the two hybrids was different.

QTL mapping for combining ability in different population-based NCII designs: a simulation study

The NCII design (North Carolina mating design II) has been widely applied in studies of combining ability and heterosis. The objective of our research was to estimate how different base populations, sample sizes, testcross numbers and heritability influence QTL analyses of combining ability and heterosis. A series of Monte Carlo simulation experiments with QTL mapping were then conducted for the base population performance, testcross population phenotypic values and the general combining ability (GCA), specific combining ability (SCA) and Hmp (midparental heterosis) datasets. The results indicated that: (i) increasing the number of testers did not necessarily enhance the QTL detection power for GCA, but it was significantly related to the QTL effect. (ii) The QTLs identified in the base population may be different from those from GCA dataset. Similar phenomena can be seen from QTL detected in SCA and Hmp datasets. (iii) The QTL detection power for GCA ranked in the order of DH(RIL) based > F2 based > BC based NCII design, when the heritability was low. The recombinant inbred lines (RILs) (or DHs) allows more recombination and offers higher mapping resolution than other populations. Further, their testcross progeny can be repeatedly generated and phenotyped. Thus, RIL based (or DH based) NCII design was highly recommend for combining ability QTL analysis. Our results expect to facilitate selecting elite parental lines with high combining ability and for geneticists to research the genetic basis of combining ability.

Nuclear DNA content in Miscanthus sp. and the geographical variation pattern in Miscanthus lutarioriparius

The genome sizes of five Miscanthus species, including 79 accessions of M. lutarioriparius, 8 of M. floridulus, 6 of M. sacchariflorus, 7 of M. sinensis, and 4 of M. × giganteus were examined using flow cytometry. The overall average nuclear DNA content were 4.256 ± 0.6 pg/2C in M. lutarioriparius, 5.175 ± 0.3 pg/2C in M. floridulus, 3.956 ± 0.2 pg/2C in M. sacchariflorus, 5.272 ± 0.2 pg/2C in M. sinensis, and 6.932 ± 0.1 pg/2C in M. × giganteus. Interspecific variation was found at the diploid level, suggesting that DNA content might be a parameter that can be used to differentiate the species. Tetraploid populations were found in M. lutarioriparius, M. sacchariflorus, and M. sinensis, and their DNA content were 8.34 ± 1.2, 8.52, and 8.355 pg, respectively. The association between the DNA content of M. lutarioriparius, collected from representative ranges across the Yangtze River, and its geographic distribution was statistically analyzed. A consistent pattern of DNA content variation in 79 M. lutarioriparius accessions across its entire geographic range was found in this study. Along the Yangtze River, the DNA content of M. lutarioriparius tended to increase from the upstream to the downstream areas, and almost all tetraploids gathered in the upstream area extended to coastal regions.

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.

High embryogenic ability and regeneration from floral axis of Amorphophallus konjac (Araceae)

Abstract Amorphophallus konjac (Araceae) a perennial herb, it has high medicinal and industrial value. In this study, a simple and efficient system for direct somatic embryogenesis and plantlet regeneration of Amorphophallus konjac was developed. The floral axis was used as the experimental material. The primary callus, developed from the floral axis grown on Murashige and Skoog (MS) medium supplemented with different hormone combination at different concentrations. The highest rate of embryogenic callus formation was observed on the MS medium containing 9.04 µM 2, 4-dichlorophenoxyacetic acid (2, 4-D) and 5.37 µM naphthalene acetic acid (NAA). The maximum induction rate was 79.8%, and the embryogenic calli were able to subculture on a medium containing similar hormone combination for over 1 year. The calli were also placed on different media for regeneration and it produced complete plants with shoots and root systems simultaneously. The highest differentiation rate of the embryogenic calli grown on differentiation medium supplemented with 8.88 µM 6-benzylaminopurine (6-BA) and 5.37 µM NAA was 95.6%. Flow cytometry analysis showed no ploidy variation in all the regenerate plantlets.

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.