Transcriptome analysis of the 2,4-dichlorophenoxyacetic acid (2,4-D)-tolerant cotton chromosome substitution line CS-B15sh and its susceptible parental lines G. hirsutum L. cv. Texas Marker-1 and G. barbadense L. cv. Pima 379

Loida M Perez; Ramil Mauleon; Mark A Arick 2nd; Zenaida V Magbanua; Daniel G Peterson; Jeffrey F D Dean; Te Ming Tseng

doi:10.3389/fpls.2022.910369

Transcriptome analysis of the 2,4-dichlorophenoxyacetic acid (2,4-D)-tolerant cotton chromosome substitution line CS-B15sh and its susceptible parental lines G. hirsutum L. cv. Texas Marker-1 and G. barbadense L. cv. Pima 379

Front Plant Sci. 2022 Aug 22:13:910369. doi: 10.3389/fpls.2022.910369. eCollection 2022.

Authors

Loida M Perez¹, Ramil Mauleon², Mark A Arick 2nd³, Zenaida V Magbanua³, Daniel G Peterson³, Jeffrey F D Dean¹, Te Ming Tseng⁴

Affiliations

¹ Department of Biochemistry, Molecular Biology, Entomology & Plant Pathology, Mississippi State University, Starkville, MS, United States.
² Faculty of Science and Engineering, Southern Cross University, East Lismore, NSW, Australia.
³ Institute for Genomics, Biocomputing and Biotechnology, Mississippi State University, Starkville, MS, United States.
⁴ Department of Plant and Soil Sciences, Mississippi State University, Starkville, MS, United States.

Abstract

The cotton chromosome substitution line, CS-B15sh, exhibits 41% lower injury from 2,4-D when applied at the field recommended rate of 1.12 kg ae ha^-1 (1×) than does Texas Marker-1 (TM-1). CS-B15sh was developed in the genetic background of Gossypium hirsutum L. cv TM-1 and has chromosome introgression on the short arm of chromosome 15 from Gossypium barbadense L. cv. Pima 379. In a previous experiment, we observed reduced translocation of [¹⁴C]2,4-D outside the treated leaf tissue in CS-B15sh, which contrasted with an increased translocation of the herbicide in the tissues above and below the treated leaf in TM-1. Our results indicate a potential 2,4-D tolerance mechanism in CS-B15sh involving altered movement of 2,4-D. Here, we used RNA sequencing (RNA-seq) to determine the differential expression of genes between 2,4-D-challenged and control plants of the tolerant (CS-B15sh) and susceptible lines (TM-1 and Pima 379). Several components of the 2,4-D/auxin-response pathway-including ubiquitin E3 ligase, PB1|AUX/IAA, ARF transcription factors, and F-box proteins of the SCF^TIR1/AFB complex-were upregulated with at least threefold higher expression in TM-1 compared with CS-B15sh, while both Pima 379 and TM-1 showed the same fold change expression for PB1|AUX/IAA mRNA. Some genes associated with herbicide metabolism, including flavin monooxygenase (Gohir.A01G174100) and FAD-linked oxidase (Gohir.D06G002600), exhibited at least a twofold increase in CS-B15sh than in TM-1 (the gene was not expressed in Pima 379), suggesting a potential relationship between the gene's expression and 2,4-D tolerance. It is interesting to note that glutathione S-transferase was differentially expressed in both CS-B15sh and Pima 379 but not in TM-1, while cytochrome P450 and other genes involved in the oxidation-reduction process were significantly expressed only in CS-B15sh in response to 2,4-D. Gene set enrichment analysis on the union DEGs of the three cotton genotypes revealed the depletion of transcripts involved in photosynthesis and enrichment of transcripts involved in ABA response and signaling.

Keywords: Illumina sequencing; RNA sequencing; abiotic stress tolerance; auxin response and signaling; herbicide metabolism; herbicide resistance mechanism; plant physiology; transcriptome analysis.