Cytogenetic and genomic characterization of a novel tall wheatgrass-derived Fhb7 allele integrated into wheat B genome.

KEY MESSAGE: We identified and integrated the novel FHB-resistant Fhb7The2 allele into wheat B genome and made it usable in both common and durum wheat breeding programs without yellow flour linkage drag. A novel tall wheatgrass-derived (Thinopyrum elongatum, genome EE) Fhb7 allele, designated Fhb7The2, was identified and integrated into the wheat B genome through a small 7B-7E translocation (7BS·7BL-7EL) involving the terminal regions of the long arms. Fhb7The2 conditions significant Type II resistance to Fusarium head blight (FHB) in wheat. Integration of Fhb7The2 into the wheat B genome makes this wild species-derived FHB resistance gene usable for breeding in both common and durum wheat. By contrast, other Fhb7 introgression lines involving wheat chromosome 7D can be utilized only in common wheat breeding programs, not in durum wheat. Additionally, we found that Fhb7The2 does not have the linkage drag of the yellow flour pigment gene that is tightly linked to the decaploid Th. ponticum-derived Fhb7 allele Fhb7Thp. This will further improve the utility of Fhb7The2 in wheat breeding. DNA sequence analysis identified 12 single nucleotide polymorphisms (SNPs) in Fhb7The2, Fhb7Thp, and another Th. elongatum-derived Fhb7 allele Fhb7The1, which led to seven amino acid conversions in Fhb7The2, Fhb7Thp, and Fhb7The1, respectively. However, no significant variation was observed in their predicted protein configuration as a glutathione transferase. Diagnostic DNA markers were developed specifically for Fhb7The2. The 7EL segment containing Fhb7The2 in the translocation chromosome 7BS·7BL-7EL exhibited a monogenic inheritance pattern in the wheat genetic background. This will enhance the efficacy of marker-assisted selection for Fhb7The2 introgression, pyramiding, and deployment in wheat germplasm and varieties.

Genomic compatibility and inheritance of hexaploid-derived Fusarium head blight resistance genes in durum wheat.

Hexaploid-derived resistance genes exhibit complex inheritance and expression patterns in tetraploid backgrounds. This study aimed to characterize the inheritance patterns and genomic compatibilities of hexaploid-derived Fusarium head blight (FHB) resistance genes in tetraploid durum wheat (Triticum durum Desf.). Evaluation of FHB resistance for F1 hybrids of hexaploid 'Sumai 3' crossed with tetraploid and hexaploid wheats indicated that Sumai 3-derived FHB resistance genes exhibit a dominant phenotypic effect seen only in hexaploid hybrids. Alternately, the hexaploid-derived FHB resistance genes from PI 277012 exhibited complete dominance in the crosses with both tetraploid and hexaploid wheat. FHB evaluation of the F1 hybrids of Sumai 3 and PI 277012 crossed with 'Langdon' (LDN)-'Chinese Spring' D-genome substitution lines suggested that chromosomes 2B, 3B, 4B, 5B, 6B, 3A, 4A, 6A, and 7A contain genes that suppress expression of the Sumai 3-derived FHB resistance, whereas chromosomes 4A, 6A, and 6B contain genes required for expression of PI 277012-derived FHB resistance. A wide range of segregation for FHB severity (10-90%) was observed in the F2 generation of Sumai 3 crossed with durum cultivars LDN and 'Divide', but the distribution of F3 families derived from the most resistant F2 segregants was skewed towards susceptibility. Similar segregation trends were observed in the crosses of PI 277012 with other durum wheats, whereby FHB resistance became slightly diluted over successive generations. These results suggest tetraploid durum wheat contains the unique alleles at multiple gene loci on different chromosomes that positively and/or negatively regulate the expression of hexaploid-derived FHB resistance genes, which complicate efforts to deploy these genes in durum breeding programs.

Dissection and physical mapping of wheat chromosome 7B by inducing meiotic recombination with its homoeologues in Aegilops speltoides and Thinopyrum elongatum.

KEY MESSAGE: We constructed a homoeologous recombination-based bin map of wheat chromosome 7B, providing a unique physical framework for further study of chromosome 7B and its homoeologues in wheat and its relatives. Homoeologous recombination leads to the dissection and diversification of the wheat genome. Advances in genome sequencing and genotyping have dramatically improved the efficacy and throughput of homoeologous recombination-based genome studies and alien introgression in wheat and its relatives. In this study, we aimed to physically dissect and map wheat chromosome 7B by inducing meiotic recombination of chromosome 7B with its homoeologues 7E in Thinopyrum elongatum and 7S in Aegilops speltoides. The special genotypes, which were double monosomic for chromosomes 7B' + 7E' or 7B' + 7S' and homozygous for the ph1b mutant, were produced to enhance 7B - 7E and 7B - 7S recombination. Chromosome-specific DNA markers were developed and used to pre-screen the large recombination populations for 7B - 7E and 7B - 7S recombinants. The DNA marker-mediated preselections were verified by fluorescent genomic in situ hybridization (GISH). In total, 29 7B - 7E and 61 7B - 7S recombinants and multiple chromosome aberrations were recovered and delineated by GISH and the wheat 90 K SNP assay. Integrated GISH and SNP analysis of the recombinants physically mapped the recombination breakpoints and partitioned wheat chromosome 7B into 44 bins with 523 SNPs assigned within. A composite bin map was constructed for chromosome 7B, showing the bin size and physical distribution of SNPs. This provides a unique physical framework for further study of chromosome 7B and its homoeologues. In addition, the 7B - 7E and 7B - 7S recombinants extend the genetic variability of wheat chromosome 7B and represent useful germplasm for wheat breeding. Thereby, this genomics-enabled chromosome engineering approach facilitates wheat genome study and enriches the gene pool of wheat improvement.

Partitioning and physical mapping of wheat chromosome 3B and its homoeologue 3E in Thinopyrum elongatum by inducing homoeologous recombination.

KEY MESSAGE: We performed homoeologous recombination-based partitioning and physical mapping of wheat chromosome 3B and Th. elongatum chromosome 3E, providing a unique physical framework of this homoeologous pair for genome studies. The wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and Thinopyrum elongatum (2n = 2x = 14, EE) genomes can be differentiated from each other by fluorescent genomic in situ hybridization (FGISH) as well as molecular markers. This has facilitated homoeologous recombination-based partitioning and engineering of their genomes for physical mapping and alien introgression. Here, we constructed a special wheat genotype, which was double monosomic for wheat chromosome 3B and Th. elongatum chromosome 3E and homozygous for the ph1b mutant, to induce 3B-3E homoeologous recombination. Totally, 81 3B-3E recombinants were recovered and detected in the primary, secondary, and tertiary homoeologous recombination cycles by FGISH. Comparing to the primary recombination, the secondary and tertiary recombination shifted toward the proximal regions due to the increase in homology between the pairing partners. The 3B-3E recombinants were genotyped by high-throughput wheat 90-K single nucleotide polymorphism (SNP) arrays and their recombination breakpoints physically mapped based on the FGISH patterns and SNP results. The 3B-3E recombination physically partitioned chromosome 3B into 38 bins, and 429 SNPs were assigned to the distinct bins. Integrative analysis of FGISH and SNP results led to the construction of a composite bin map for chromosome 3B. Additionally, we developed 22 SNP-derived semi-thermal asymmetric reverse PCR markers specific for chromosome 3E and constructed a comparative map of homoeologous chromosomes 3E, 3B, 3A, and 3D. In summary, this work provides a unique physical framework for further studies of the 3B-3E homoeologous pair and diversifies the wheat genome for wheat improvement.

Chromosome engineering-mediated introgression and molecular mapping of novel Aegilops speltoides-derived resistance genes for tan spot and Septoria nodorum blotch diseases in wheat.

KEY MESSAGE: We identified, mapped and introduced novel Aegilops speltoides-derived resistance genes for tan spot and SNB diseases into wheat, enhancing understanding and utilization of host resistance to both diseases in wheat. Tan spot and Septoria nodorum blotch (SNB) are two important fungal diseases of wheat. Resistance to these diseases is often observed as the lack of sensitivity to the necrotrophic effectors (NE) produced by the fungal pathogens and thus exhibits a recessive inheritance pattern. In this study, we identified novel genes for resistance to tan spot and SNB on Aegilops speltoides (2n = 2x = 14, genome SS) chromosome 2S. These genes confer dominant resistance in the wheat background, indicating a distinct NE-independent mechanism of resistance. Ae. speltoides chromosome 2S was engineered for resistance gene introgression and molecular mapping by inducing meiotic homoeologous recombination with wheat chromosome 2B. Twenty representative 2B-2S recombinants were evaluated for reaction to tan spot and SNB and were delineated by genomic in situ hybridization and high-throughput wheat 90 K SNP assay. The resistance genes physically mapped to the sub-telomeric region (~ 8 Mb) on the short arm of chromosome 2S and designated TsrAes1 for tan spot resistance and SnbAes1 for SNB resistance. In addition, we developed SNP-derived PCR markers closely linked to TsrAes1/SnbAes1 for marker-assisted selection in wheat breeding. TsrAes1 and SnbAes1 are the first set of NE-independent tan spot, and SNB resistance genes are identified from Ae. speltoides. The 2SS-2BS·2BL recombinants with minimal amounts of Ae. speltoides chromatin containing TsrAes1/SnbAes1 were produced for germplasm development, making the wild species-derived resistance genes usable in wheat breeding. This will strengthen and diversify resistance of wheat to tan spot and SNB and facilitate understanding of resistance to these two diseases.

Meiotic homoeologous recombination-based mapping of wheat chromosome 2B and its homoeologues in Aegilops speltoides and Thinopyrum elongatum.

KEY MESSAGE: We physically dissected and mapped wheat chromosome 2B and its homoeologues in Aegilops speltoides and Thinopyrum elongatum based on meiotic homoeologous recombination, providing a unique physical framework for genome studies. Common wheat has a large and complex genome with narrow genetic diversity and various degrees of recombination between the A, B, and D subgenomes. This has limited the homologous recombination-based genome studies in wheat. Here, we exploited meiotic homoeologous recombination for molecular mapping of wheat chromosome 2B and its homoeologue 2S from Aegilops speltoides and 2E from Thinopyrum elongatum. The 2B-2S and 2B-2E recombination was induced by the ph1b mutant, and recovered using molecular markers and fluorescent genomic in situ hybridization (FGISH). A total of 112 2B-2S and 87 2B-2E recombinants involving different chromosome regions were developed and physically delineated by FGISH. The 2B-2S and 2B-2E recombination hotspots mapped to the subterminal regions on both arms. Recombination hotspots with the highest recombination rates mapped to the short arms. Eighty-three 2B-2S and 67 2B-2E recombinants were genotyped using the wheat 90 K SNP arrays. Based on the genotyping results and FGISH patterns of the recombinants, chromosomes 2B, 2S, and 2E were partitioned into 93, 66, and 46 bins, respectively. In total, 1037 SNPs physically mapped onto distinct bins of these three homoeologous chromosomes. A homoeologous recombination-based bin map was constructed for chromosome 2B, providing a unique physical framework for genome studies in wheat and its relatives. Meiotic homoeologous recombination also facilitates gene introgression to diversify the wheat genome for germplasm development. Therefore, homoeologous recombination-based studies enhance understanding of the wheat genome and its homoeologous counterparts from wild grasses, and expand the genetic variability of the wheat genome.

Genome-wide characterization of differentially expressed genes provides insights into regulatory network of heat stress response in radish (Raphanus sativus L.).

Heat stress (HS) causes detrimental effects on plant morphology, physiology, and biochemistry that lead to drastic reduction in plant biomass production and economic yield worldwide. To date, little is known about HS-responsive genes involved in thermotolerance mechanism in radish. In this study, a total of 6600 differentially expressed genes (DEGs) from the control and Heat24 cDNA libraries of radish were isolated by high-throughput sequencing. With Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, some genes including MAPK, DREB, ERF, AP2, GST, Hsf, and Hsp were predominantly assigned in signal transductions, metabolic pathways, and biosynthesis and abiotic stress-responsive pathways. These pathways played significant roles in reducing stress-induced damages and enhancing heat tolerance in radish. Expression patterns of 24 candidate genes were validated by reverse-transcription quantitative PCR (RT-qPCR). Based mainly on the analysis of DEGs combining with the previous miRNAs analysis, the schematic model of HS-responsive regulatory network was proposed. To counter the effects of HS, a rapid response of the plasma membrane leads to the opening of specific calcium channels and cytoskeletal reorganization, after which HS-responsive genes are activated to repair damaged proteins and ultimately facilitate further enhancement of thermotolerance in radish. These results could provide fundamental insight into the regulatory network underlying heat tolerance in radish and facilitate further genetic manipulation of thermotolerance in root vegetable crops.

Differential proteomic analysis reveals sequential heat stress-responsive regulatory network in radish (Raphanus sativus L.) taproot.

MAIN CONCLUSION: Differential abundance protein species (DAPS) involved in reducing damage and enhancing thermotolerance in radish were firstly identified. Proteomic analysis and omics association analysis revealed a HS-responsive regulatory network in radish. Heat stress (HS) is a major destructive factor influencing radish production and supply in summer, for radish is a cool season vegetable crop being susceptible to high temperature. In this study, the proteome changes of radish taproots under 40 °C treatment at 0 h (Control), 12 h (Heat12) and 24 h (Heat24) were analyzed using iTRAQ (Isobaric Tag for Relative and Absolute Quantification) approach. In total, 2258 DAPS representing 1542 differentially accumulated uniprotein species which respond to HS were identified. A total of 604, 910 and 744 DAPS was detected in comparison of Control vs. Heat12, Control vs. Heat24, and Heat12 vs. Heat24, respectively. Gene ontology and pathway analysis showed that annexin, ubiquitin-conjugating enzyme, ATP synthase, heat shock protein (HSP) and other stress-related proteins were predominately enriched in signal transduction, stress and defense pathways, photosynthesis and energy metabolic pathways, working cooperatively to reduce stress-induced damage in radish. Based on iTRAQ combined with the transcriptomics analysis, a schematic model of a sequential HS-responsive regulatory network was proposed. The initial sensing of HS occurred at the plasma membrane, and then key components of stress signal transduction triggered heat-responsive genes in the plant protective metabolism to re-establish homeostasis and enhance thermotolerance. These results provide new insights into characteristics of HS-responsive DAPS and facilitate dissecting the molecular mechanisms underlying heat tolerance in radish and other root crops.

Molecular cytogenetic and genomic analyses reveal new insights into the origin of the wheat B genome.

KEY MESSAGE: This work pinpointed the goatgrass chromosomal segment in the wheat B genome using modern cytogenetic and genomic technologies, and provided novel insights into the origin of the wheat B genome. Wheat is a typical allopolyploid with three homoeologous subgenomes (A, B, and D). The donors of the subgenomes A and D had been identified, but not for the subgenome B. The goatgrass Aegilops speltoides (genome SS) has been controversially considered a possible candidate for the donor of the wheat B genome. However, the relationship of the Ae. speltoides S genome with the wheat B genome remains largely obscure. The present study assessed the homology of the B and S genomes using an integrative cytogenetic and genomic approach, and revealed the contribution of Ae. speltoides to the origin of the wheat B genome. We discovered noticeable homology between wheat chromosome 1B and Ae. speltoides chromosome 1S, but not between other chromosomes in the B and S genomes. An Ae. speltoides-originated segment spanning a genomic region of approximately 10.46 Mb was detected on the long arm of wheat chromosome 1B (1BL). The Ae. speltoides-originated segment on 1BL was found to co-evolve with the rest of the B genome. Evidently, Ae. speltoides had been involved in the origin of the wheat B genome, but should not be considered an exclusive donor of this genome. The wheat B genome might have a polyphyletic origin with multiple ancestors involved, including Ae. speltoides. These novel findings will facilitate genome studies in wheat and other polyploids.

Genome-wide characterization and expression profiling of NAC transcription factor genes under abiotic stresses in radish (Raphanus sativus L.).

NAC (NAM, no apical meristem; ATAF, Arabidopsis transcription activation factor and CUC, cup-shaped cotyledon) proteins are among the largest transcription factor (TF) families playing fundamental biological processes, including cell expansion and differentiation, and hormone signaling in response to biotic and abiotic stresses. In this study, 172 RsNACs comprising 17 membrane-bound members were identified from the whole radish genome. In total, 98 RsNAC genes were non-uniformly distributed across the nine radish chromosomes. In silico analysis revealed that expression patterns of several NAC genes were tissue-specific such as a preferential expression in roots and leaves. In addition, 21 representative NAC genes were selected to investigate their responses to heavy metals (HMs), salt, heat, drought and abscisic acid (ABA) stresses using real-time polymerase chain reaction (RT-qPCR). As a result, differential expressions among these genes were identified where RsNAC023 and RsNAC080 genes responded positively to all stresses except ABA, while RsNAC145 responded more actively to salt, heat and drought stresses compared with other genes. The results provides more valuable information and robust candidate genes for future functional analysis for improving abiotic stress tolerances in radish.

Coordinated Regulation of Anthocyanin Biosynthesis Genes Confers Varied Phenotypic and Spatial-Temporal Anthocyanin Accumulation in Radish (Raphanus sativus L.).

Anthocyanins are natural pigments that have important functions in plant growth and development. Radish taproots are rich in anthocyanins which confer different taproot colors and are potentially beneficial to human health. The crop differentially accumulates anthocyanin during various stages of growth, yet molecular mechanisms underlying this differential anthocyanin accumulation remains unknown. In the present study, transcriptome analysis was used to concisely identify putative genes involved in anthocyanin biosynthesis in radish. Spatial-temporal transcript expressions were then profiled in four color variant radish cultivars. From the total transcript sequences obtained through illumina sequencing, 102 assembled unigenes, and 20 candidate genes were identified to be involved in anthocyanin biosynthesis. Fifteen genomic sequences were isolated and sequenced from radish taproot. The length of these sequences was between 900 and 1,579 bp, and the unigene coverage to all of the corresponding cloned sequences was more than 93%. Gene structure analysis revealed that RsF3'H is intronless and anthocyanin biosynthesis genes (ABGs) bear asymmetrical exons, except RsSAM. Anthocyanin accumulation showed a gradual increase in the leaf of the red radish and the taproot of colored cultivars during development, with a rapid increase at 30 days after sowing (DAS), and the highest content at maturity. Spatial-temporal transcriptional analysis of 14 genes revealed detectable expressions of 12 ABGs in various tissues at different growth levels. The investigation of anthocyanin accumulation and gene expression in four color variant radish cultivars, at different stages of development, indicated that total anthocyanin correlated with transcript levels of ABGs, particularly RsUFGT, RsF3H, RsANS, RsCHS3 and RsF3'H1. Our results suggest that these candidate genes play key roles in phenotypic and spatial-temporal anthocyanin accumulation in radish through coordinated regulation and the major control point in anthocyanin biosynthesis in radish is RsUFGT. The present findings lend invaluable insights into anthocyanin biosynthesis and may facilitate genetic manipulation for enhanced anthocyanin content in radish.

Genome-wide characterization of the WRKY gene family in radish (Raphanus sativus L.) reveals its critical functions under different abiotic stresses.

KEY MESSAGE: The radish WRKY gene family was genome-widely identified and played critical roles in response to multiple abiotic stresses. The WRKY is among the largest transcription factors (TFs) associated with multiple biological activities for plant survival, including control response mechanisms against abiotic stresses such as heat, salinity, and heavy metals. Radish is an important root vegetable crop and therefore characterization and expression pattern investigation of WRKY transcription factors in radish is imperative. In the present study, 126 putative WRKY genes were retrieved from radish genome database. Protein sequence and annotation scrutiny confirmed that RsWRKY proteins possessed highly conserved domains and zinc finger motif. Based on phylogenetic analysis results, RsWRKYs candidate genes were divided into three groups (Group I, II and III) with the number 31, 74, and 20, respectively. Additionally, gene structure analysis revealed that intron-exon patterns of the WRKY genes are highly conserved in radish. Linkage map analysis indicated that RsWRKY genes were distributed with varying densities over nine linkage groups. Further, RT-qPCR analysis illustrated the significant variation of 36 RsWRKY genes under one or more abiotic stress treatments, implicating that they might be stress-responsive genes. In total, 126 WRKY TFs were identified from the R. sativus genome wherein, 35 of them showed abiotic stress-induced expression patterns. These results provide a genome-wide characterization of RsWRKY TFs and baseline for further functional dissection and molecular evolution investigation, specifically for improving abiotic stress resistances with an ultimate goal of increasing yield and quality of radish.

Identification of critical genes associated with lignin biosynthesis in radish (Raphanus sativus L.) by de novo transcriptome sequencing.

Radish is an important root vegetable crop with high nutritional, economic, and medicinal value. Lignin is an important secondary metabolite possessing a great effect on plant growth and product quality. To date, lignin biosynthesis-related genes have been identified in some important plant species. However, little information on characterization of critical genes involved in plant lignin biosynthesis is available in radish. In this study, a total of 71,148 transcripts sequences were obtained from radish root, of which 66 assembled unigenes and ten candidate genes were identified to be involved in lignin monolignol biosynthesis. Full-length cDNA sequences of seven randomly selected genes were isolated and sequenced from radish root, and the assembled unigenes covered more than 80% of their corresponding cDNA sequences. Moreover, the lignin content gradually accumulated in leaf during the developmental stages, and it increased from pre-cortex to cortex splitting stage, followed by a decrease at thickening stage and then increased at mature stage in root. RT-qPCR analysis revealed that all these genes except RsF5H exhibited relatively low expression level in root at thickening stage. The expression profiles of Rs4CL5, RsCCoAOMT1, and RsCOMT genes were consistent with the changes of root lignin content, implying that these candidate genes may play important roles in lignin formation in radish root. These findings would provide valuable information for identification of lignin biosynthesis-related genes and facilitate dissection of molecular mechanism underlying lignin biosynthesis in radish and other root vegetable crops.

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.).

BACKGROUND: Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish. RESULTS: In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions. CONCLUSIONS: These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.

Transcriptome Profiling of Taproot Reveals Complex Regulatory Networks during Taproot Thickening in Radish (Raphanus sativus L.).

Radish (Raphanus sativus L.) is one of the most important vegetable crops worldwide. Taproot thickening represents a critical developmental period that determines yield and quality in radish life cycle. To isolate differentially expressed genes (DGEs) involved in radish taproot thickening process and explore the molecular mechanism underlying taproot development, three cDNA libraries from radish taproot collected at pre-cortex splitting stage (L1), cortex splitting stage (L2), and expanding stage (L3) were constructed and sequenced by RNA-Seq technology. More than seven million clean reads were obtained from the three libraries, from which 4,717,617 (L1, 65.35%), 4,809,588 (L2, 68.24%) and 4,973,745 (L3, 69.45%) reads were matched to the radish reference genes, respectively. A total of 85,939 transcripts were generated from three libraries, from which 10,450, 12,325, and 7392 differentially expressed transcripts (DETs) were detected in L1 vs. L2, L1 vs. L3, and L2 vs. L3 comparisons, respectively. Gene Ontology and pathway analysis showed that many DEGs, including EXPA9, Cyclin, CaM, Syntaxin, MADS-box, SAUR, and CalS were involved in cell events, cell wall modification, regulation of plant hormone levels, signal transduction and metabolisms, which may relate to taproot thickening. Furthermore, the integrated analysis of mRNA-miRNA revealed that 43 miRNAs and 92 genes formed 114 miRNA-target mRNA pairs were co-expressed, and three miRNA-target regulatory networks of taproot were constructed from different libraries. Finally, the expression patterns of 16 selected genes were confirmed using RT-qPCR analysis. A hypothetical model of genetic regulatory network associated with taproot thickening in radish was put forward. The taproot formation of radish is mainly attributed to cell differentiation, division and expansion, which are regulated and promoted by certain specific signal transduction pathways and metabolism processes. These results could provide new insights into the complex molecular mechanism underlying taproot thickening and facilitate genetic improvement of taproot in radish.

Identification of microRNAs and Their Target Genes Explores miRNA-Mediated Regulatory Network of Cytoplasmic Male Sterility Occurrence during Anther Development in Radish (Raphanus sativus L.).

MicroRNAs (miRNAs) are a type of endogenous non-coding small RNAs that play critical roles in plant growth and developmental processes. Cytoplasmic male sterility (CMS) is typically a maternally inherited trait and widely used in plant heterosis utilization. However, the miRNA-mediated regulatory network of CMS occurrence during anther development remains largely unknown in radish. In this study, a comparative small RNAome sequencing was conducted in floral buds of CMS line 'WA' and its maintainer line 'WB' by high-throughput sequencing. A total of 162 known miRNAs belonging to 25 conserved and 24 non-conserved miRNA families were isolated and 27 potential novel miRNA families were identified for the first time in floral buds of radish. Of these miRNAs, 28 known and 14 potential novel miRNAs were differentially expressed during anther development. Several target genes for CMS occurrence-related miRNAs encode important transcription factors and functional proteins, which might be involved in multiple biological processes including auxin signaling pathways, signal transduction, miRNA target silencing, floral organ development, and organellar gene expression. Moreover, the expression patterns of several CMS occurrence-related miRNAs and their targets during three stages of anther development were validated by qRT-PCR. In addition, a potential miRNA-mediated regulatory network of CMS occurrence during anther development was firstly proposed in radish. These findings could contribute new insights into complex miRNA-mediated genetic regulatory network of CMS occurrence and advance our understanding of the roles of miRNAs during CMS occurrence and microspore formation in radish and other crops.

Comprehensive transcriptome-based characterization of differentially expressed genes involved in microsporogenesis of radish CMS line and its maintainer.

Microsporogenesis is an indispensable period for investigating microspore development and cytoplasmic male sterility (CMS) occurrence. Radish CMS line plays a critical role in elite F1 hybrid seed production and heterosis utilization. However, the molecular mechanisms of microspore development and CMS occurrence have not been thoroughly uncovered in radish. In this study, a comparative analysis of radish floral buds from a CMS line (NAU-WA) and its maintainer (NAU-WB) was conducted using next generation sequencing (NGS) technology. Digital gene expression (DGE) profiling revealed that 3504 genes were significantly differentially expressed between NAU-WA and NAU-WB library, among which 1910 were upregulated and 1594 were downregulated. Gene ontology (GO) analysis showed that these differentially expressed genes (DEGs) were mainly enriched in extracellular region, catalytic activity, and response to stimulus. KEGG enrichment analysis revealed that the DEGs were predominantly associated with flavonoid biosynthesis, glycolysis, and biosynthesis of secondary metabolites. Real-time quantitative PCR analysis showed that the expression profiles of 13 randomly selected DEGs were in high agreement with results from Illumina sequencing. Several candidate genes encoding ATP synthase, auxin response factor (ARF), transcription factors (TFs), chalcone synthase (CHS), and male sterility (MS) were responsible for microsporogenesis. Furthermore, a schematic diagram for functional interaction of DEGs from NAU-WA vs. NAU-WB library in radish plants was proposed. These results could provide new information on the dissection of the molecular mechanisms underlying microspore development and CMS occurrence in radish.

De novo Taproot Transcriptome Sequencing and Analysis of Major Genes Involved in Sucrose Metabolism in Radish (Raphanus sativus L.).

Radish (Raphanus sativus L.) is an important annual or biennial root vegetable crop. The fleshy taproot comprises the main edible portion of the plant with high nutrition and medical value. Molecular biology study of radish begun rather later, and lacks sufficient transcriptomic and genomic data in pubic databases for understanding of the molecular mechanism during the radish taproot formation. To develop a comprehensive overview of the 'NAU-YH' root transcriptome, a cDNA library, prepared from three equally mixed RNA of taproots at different developmental stages including pre-cortex splitting stage, cortex splitting stage, and expanding stage was sequenced using high-throughput Illumina RNA sequencing. From approximately 51 million clean reads, a total of 70,168 unigenes with a total length of 50.28 Mb, an average length of 717 bp and a N50 of 994 bp were obtained. In total, 63,991 (about 91.20% of the assembled unigenes) unigenes were successfully annotated in five public databases including NR, GO, COG, KEGG, and Nt. GO analysis revealed that the majority of these unigenes were predominately involved in basic physiological and metabolic processes, catalytic, binding, and cellular process. In addition, a total of 103 unigenes encoding eight enzymes involved in the sucrose metabolism related pathways were also identified by KEGG pathway analysis. Sucrose synthase (29 unigenes), invertase (17 unigenes), sucrose-phosphate synthase (16 unigenes), fructokinase (17 unigenes), and hexokinase (11 unigenes) ranked top five in these eight key enzymes. From which, two genes (RsSuSy1, RsSPS1) were validated by T-A cloning and sequenced, while the expression of six unigenes were profiled with RT-qPCR analysis. These results would be served as an important public reference platform to identify the related key genes during taproot thickening and facilitate the dissection of molecular mechanisms underlying taproot formation in radish.

Transcriptional identification and characterization of differentially expressed genes associated with embryogenesis in radish (Raphanus sativus L.).

Embryogenesis is an important component in the life cycle of most plant species. Due to the difficulty in embryo isolation, the global gene expression involved in plant embryogenesis, especially the early events following fertilization are largely unknown in radish. In this study, three cDNA libraries from ovules of radish before and after fertilization were sequenced using the Digital Gene Expression (DGE) tag profiling strategy. A total of 5,777 differentially expressed transcripts were detected based on pairwise comparison in the three libraries (0_DAP, 7_DAP and 15_DAP). Results from Gene Ontology (GO) and pathway enrichment analysis revealed that these differentially expressed genes (DEGs) were implicated in numerous life processes including embryo development and phytohormones biosynthesis. Notably, some genes encoding auxin response factor (ARF ), Leafy cotyledon1 (LEC1) and somatic embryogenesis receptor-like kinase (SERK ) known to be involved in radish embryogenesis were differentially expressed. The expression patterns of 30 genes including LEC1-2, AGL9, LRR, PKL and ARF8-1 were validated by qRT-PCR. Furthermore, the cooperation between miRNA and mRNA may play a pivotal role in the radish embryogenesis process. This is the first report on identification of DEGs profiles related to radish embryogenesis and seed development. These results could facilitate further dissection of the molecular mechanisms underlying embryogenesis and seed development in radish.

Transcriptome-based gene expression profiling identifies differentially expressed genes critical for salt stress response in radish (Raphanus sativus L.).

KEY MESSAGE: Transcriptome-based gene expression analysis identifies many critical salt-responsive genes in radish and facilitates further dissecting the molecular mechanism underlying salt stress response. Salt stress severely impacts plant growth and development. Radish, a moderately salt-sensitive vegetable crop, has been studied for decades towards the physiological and biochemical performances under salt stress. However, no systematic study on isolation and identification of genes involved in salt stress response has been performed in radish, and the molecular mechanism governing this process is still indistinct. Here, the RNA-Seq technique was applied to analyze the transcriptomic changes on radish roots treated with salt (200 mM NaCl) for 48 h in comparison with those cultured in normal condition. Totally 8709 differentially expressed genes (DEGs) including 3931 up- and 4778 down-regulated genes were identified. Functional annotation analysis indicated that many genes could be involved in several aspects of salt stress response including stress sensing and signal transduction, osmoregulation, ion homeostasis and ROS scavenging. The association analysis of salt-responsive genes and miRNAs exhibited that 36 miRNA-mRNA pairs had negative correlationship in expression trends. Reverse-transcription quantitative PCR (RT-qPCR) analysis revealed that the expression profiles of DEGs were in line with results from the RNA-Seq analysis. Furthermore, the putative model of DEGs and miRNA-mediated gene regulation was proposed to elucidate how radish sensed and responded to salt stress. This study represents the first comprehensive transcriptome-based gene expression profiling under salt stress in radish. The outcomes of this study could facilitate further dissecting the molecular mechanism underlying salt stress response and provide a valuable platform for further genetic improvement of salt tolerance in radish breeding programs.