Radish cytoplasmic male sterility (CMS) has been widely used for breeding in Raphanus and Brassica genera. However, the detailed regulation network of the male sterility remains to be determined. Our previous work has shown that the abnormalities in a CMS radish appeared shortly after the tetrad stage when microspores were malformed and the tapetal cells grew abnormally large. In this work, histological analysis shows that anthers are at the tetrad stage when the radish buds are about 1.5 mm in length. Furthermore, a high throughput RNA sequencing technology was employed to characterize the transcriptome of radish buds with length about 1.5 mm from two CMS lines possessing the CMS-inducing orf138 gene and corresponding near-isogenic maintainer lines. A total of 67,140 unigenes were functionally annotated. Functional terms for these genes are significantly enriched in 55 Gene Ontology (GO) groups and 323 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The transcriptome detected transcripts for 72 out of a total of 79 protein genes encoded in the chloroplast genome from radish. In contrast, the radish mitochondrial genome contains 34 protein genes, but only 16 protein transcripts were detected from the transcriptome. The transcriptome comparison between CMS and near-isogenic maintainer lines revealed 539 differentially expressed genes (DEGs), indicating that the false positive rate for comparative transcriptome profiling was clearly decreased using two groups of CMS/maintainer lines with different nuclear background. The level of 127 transcripts was increased and 412 transcripts were decreased in the CMS lines. No change in levels of transcripts except CMS-inducing orf138 was identified from the mitochondrial and chloroplast genomes. Some DEGs which would be associated with the CMS, encoding MYB and bHLH transcription factors, pentatricopeptide repeat (PPR) proteins, heat shock transcription factors (HSFs) and heat shock proteins (HSPs), are discussed. The transcriptome dataset and comparative analysis will provide an important resource for further understanding anther development, the CMS mechanism and to improve molecular breeding in radish.

Genome-wide association studies are a powerful approach for identifying genes related to complex traits in organisms, but are limited by the requirement for a reference genome sequence of the species under study. To circumvent this problem, we propose a transcriptome-referenced association study (TRAS) that utilizes a transcriptome generated by single-molecule long-read sequencing as a reference sequence to score population variation at both transcript sequence and expression levels. Candidate transcripts are identified when both scores are associated with a trait and their potential interactions are ascertained by expression quantitative trait loci analysis. Applying this method to characterize garlic clove shape traits in 102 landraces, we identified 22 candidate transcripts, most of which showed extensive interactions. Eight transcripts were long non-coding RNAs (lncRNAs), and the others were proteins involved mainly in carbohydrate metabolism, protein degradation, etc. TRAS, as an efficient tool for association study independent of a reference genome, extends the applicability of association studies to a broad range of species.

Radish (Raphanus sativus) is an important vegetable worldwide that exhibits different flesh and skin colors. The anthocyanins responsible for the red and purple coloring in radishes possess nutritional value and pharmaceutical potential. To explore the structural and regulatory networks related to anthocyanin biosynthesis and identify key genes, we performed comparative transcriptome analyses of the skin and flesh of six colored radish accessions. The transcript profiles showed that each accession had a species-specific transcript profile. For radish pigmentation accumulation, the expression levels of anthocyanin biosynthetic genes (RsTT4, RsC4H, RsTT7, RsCCOAMT, RsDFR, and RsLDOX) were significantly upregulated in the red- and purple-colored accessions, but were downregulated or absent in the white and black accessions. The correlation test, combined with metabolome (PCC > 0.95), revealed five structural genes (RsTT4, RsDFR, RsCCOAMT, RsF3H, and RsBG8L) and three transcription factors (RsTT8-1, RsTT8-2, and RsPAR1) to be significantly correlated with flavonoids in the skin of the taproot. Four structural genes (RsBG8L, RsDFR, RsCCOAMT, and RsLDOX) and nine transcription factors (RsTT8-1, RsTT8-2, RsMYB24L, RsbHLH57, RsPAR2L, RsbHLH113L, RsOGR3L, RsMYB24, and RsMYB34L) were found to be significantly correlated with metabolites in the flesh of the taproot. This study provides a foundation for future studies on the gene functions and genetic diversity of radish pigmentation and should aid in the cultivation of new valuable radish varieties.

The northern root-knot nematode, Meloidogyne hapla, is a biotrophic parasite that infects many crops and causes severe economic losses worldwide. Rapid and accurate detection of M. hapla is crucial for disease forecasting and control. We developed a recombinase polymerase amplification combined with a lateral flow dipstick (RPA-LFD) assay for rapid detection of M. hapla. The primers and probe were designed based on the effector gene 16D10 sequence and were highly specific to M. hapla. The RPA reaction was performed at a wide range of temperatures from 25 to 45°C within 5 to 25 min, and the amplicon was visualized directly on the LFD within 5 min. The detection limits of the RPA-LFD assay were 10⁻³ females and 10⁻² second-stage juveniles/0.5 g of soil, which was 10 times more sensitive than the conventional PCR assay. In addition, the RPA-LFD assay can detect M. hapla from infested plant roots and soil samples, and the entire detection process can be completed within 1.5 h. These results indicate that the RPA-LFD assay is a simple, rapid, specific, sensitive, and visual method that can be used for rapid detection of M. hapla in the field and in resource-limited conditions.

CCT genes play vital roles in flowering, plant growth, development, and response to abiotic stresses. Although they have been reported in many plants, the characterization and expression pattern of CCT genes is still limited in R. sativus. In this study, a total of 58 CCT genes were identified in R. sativus. Phylogenetic tree, gene structure, and conserved domains revealed that all CCT genes were classified into three groups: COL, CMF, and PRR. Genome-wide identification and evolutionary analysis showed that segmental duplication expanded the CCT gene families considerably, with the LF subgenome retaining more CCT genes. We observed strong purifying selection pressure for CCT genes. RsCCT genes showed tissue specificity, and some genes (such as RsCCT22, RsCCT36, RsCCT42 and RsCCT51) were highly expressed in flowers. Promoter cis-elements and RNA-seq data analysis showed that RsCCT genes could play roles in controlling flowering through the photoperiodic pathway and vernalization pathway. The expression profiles of RsCCT genes under Cd, Cr, Pb, and heat and salt stresses revealed that many RsCCT genes could respond to one or more abiotic stresses. Our findings could provide essential information for further studies on the function of RsCCT genes.