Similar protocols

To access compelling stats and trends, optimize your time and resources and pinpoint new correlations, you will need to subscribe to our premium service.

Subscribe

Pipeline publication

[…] sailles.inra.fr/Seq-Repository/RNA-Seq). This dataset covers five different organs (root, leaf, stem, spike, grain) at three developmental stages, each in two replicates. We used the HC gene predictions, version 2.1 (ftp://ftpmips.helmholtz-muenchen.de/plants/wheat/IWGSC/genePrediction_v2.1/) as reference for mapping the RNA-Seq reads. The reference was prepared by extracting the genomic sequence for each of the predicted genes with up to 2kb upstream and downstream bases whenever available from the corresponding IWGSC-CSS contig. The coordinates in the corresponding GTF/GFF transcript annotations file were adjusted accordingly. RNA-Seq reads were quality-, adapter- and length-trimmed using Trimmomatic (), version 0.30, with a custom list of adapter sequences and the following settings: ‘ILLUMINACLIP:adapters.fa:1:6:6 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:6 MINLEN:60’. After indexing the reference using Bowtie2 () version 2.2.1, trimmed reads were aligned to the reference using TopHat () version 2.2.1, not allowing any mismatches or indels. Paired reads were required to map concordantly (--no-discordant setting) to the same (--no-mixed setting) reference sequence. BAM files for the biological replicates were merged, apart from a single tissue/stage sample (spike_Z39) where only one sample was available. Expression was quantified by Cufflinks () version 2.1.1 utilizing (through the --GTF option) the adjusted version of the reference transcript annotations provided with the HC gene predictions. The remaining settings were left at their defaults, except for: --max-multiread-fraction 1, --frag-len-std-dev 50 --max-intron-length 5 000. FPKMs (fragments per kilobase of exon per million fragments mapped) per gene (rather than per isoform) were extracted and aggregated in tabular form. Functional annotation was obtained by comparison of the HC gene predictions to MSU Rice Genome, Release 7, and Arabidopsis TAIR using BlastX (expect-value <10−10)., For each genotype a total of 21 daytime whole-plant TR measurements (seven replicated three times) were carried out in the 1.1 to 3.7 kPa range simultaneously for the entire population (). The corresponding lowest (TR1.1) and highest (TR3.7) TR values were significantly variable between genotypes, reflecting a highly transgressive segregation in this population (). The genetic variability of TR3.7 was large, ranging between 94.1 and 236.7mg H2O m−2 s−1, dependent on the genotype. This variability was confirmed for TR values measured at the lowest VPD (TR1.1), which ranged from 20.3 to 55.1mg H2O m−2 s−1 kPa−1. This variation had a strong genetic basis as ind […]

Pipeline specifications

Software tools Trimmomatic, Bowtie2, TopHat, Cufflinks, BLASTX
Organisms Triticum aestivum