For RNA-seq library preparation, total RNA was treated with
Ambion Turbo DNase (Thermo Fisher Scientific, AM1907) and quantified using the
Qubit RNA BR Assay Kit (Invitrogen, Q10210). Four biological replicates (
n = 4) were used for 0-min samples of Col-0,
camta3-1,
fer-4, and
piezo and 22-min samples of Col-0 and
piezo. Three biological replicates (
n = 3) were used for 22-min samples of
camta3-1 and
fer-4. Then, 500 ng of RNA was used for library preparation using Illumina
TruSeq mRNA (poly-A selection) and
TruSeq RNA UD Indexes for up to 96 samples (Illumina, 20022371). Samples were pooled and sequenced on a half Illumina
NovaSeq6000 S4 lane, 2 × 150 bp reads, incl Xp kit. Data were processed using demultiplexing and quality controlled with FastQC. Alignment of reads was performed against the TAIR11 annotation using STAR (
67 (
link)). On average, 32 million reads per sample were generated. Counts were assigned to genes using featureCounts (
68 (
link)), and analysis of DEGs was performed with DeSEQ2 (
69 (
link)). Transcripts were considered differentially expressed if padj < 0.05 and fold change ≥ 1.5 or fold change ≤ −1.5. Raw RNA-seq data files are available on ArrayExpress (accession number E-MTAB- 10920). GO enrichment analysis was performed using TF2Network and GOrilla (
70 (
link)).
Darwish E., Ghosh R., Ontiveros-Cisneros A., Tran H.C., Petersson M., De Milde L., Broda M., Goossens A., Van Moerkercke A., Khan K, & Van Aken O. (2022). Touch signaling and thigmomorphogenesis are regulated by complementary CAMTA3- and JA-dependent pathways. Science Advances, 8(20), eabm2091.