Identification of ceRNA and candidate genes related to fertility conversion of TCMS line YS3038 in wheat

Previous studies have indicated that noncoding RNAs are important factors in gene functions. To explore the mechanism of male sterility of YS3038, the sterile genes were mapped, and based on previous work, the expression of long noncoding RNAs (lncRNAs), circular RNAs (circRNAs), and their target genes was studied. Weighted gene coexpression network analysis (WGCNA) and competitive endogenous RNA (ceRNA) analysis were further performed for differentially expressed noncoding RNAs and target genes. At last, the candidate genes were silenced by barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) to prove their function. The sterile genes were mapped on chromosomes 1B and 6B based on chip mix pool analysis, and one major effect QTL (27.3190% variation) was found based on SSR primers. The WGCNA analysis revealed that the dark turquoise and steel blue modules were highly correlated with anther development and fertility conversion, respectively. The ceRNA analysis showed that a total of 184 RNAs interacted with each other, including 115 mRNAs, 55 microRNAs (miRNAs), eight circRNAs, and six lncRNAs. Finally, the seed setting rate of the plant was significantly decreased after fatty acyl-CoA reductase 5 silencing. This study provides breeders with a new option for the development of thermosensitive cytoplasmic male-sterile (TCMS) wheat lines, which will favor the sustainable development of two-line hybrid wheat.