Genetic basis of sRNA quantitative variation analyzed using an experimental population derived from an elite rice hybrid
- Huazhong Agricultural University, China
Figures
Figure 1 with 5 supplements
The distribution of sRNAs and s-traits in different genomic regions across the IMF2 population.
(A) The distribution of sRNAs across the IMF2 population. (B) The distribution of sRNAs aligned to the genic region, 2-kb upstream, and 500-bp downstream of annotated genes, as well as intergenic …
Figure 1—figure supplement 1
The number of sRNA reads of different lengths in the raw sequencing data of the IMF2 population.
The number of sRNAs with specific length for all the IMF2s is shown as box plot.
Figure 1—figure supplement 2
Processing of raw sRNA reads to genome-mapped reads.
The number of sRNAs of specific length in the IMF2 population is shown as box plot.
Figure 1—figure supplement 3
Distribution of sRNAs in the upstream and downstream of genic regions.
The distribution of the distance between sRNAs and transcription start site (TSS) (2-kb upstream) (A) and transcription termination site (TTS) (500-bp downstream) (B).
Figure 1—figure supplement 4
The distribution of coefficients of variation (CV) for 165,797 s-traits.
https://doi.org/10.7554/eLife.03913.008
Figure 1—figure supplement 5
The distributions of the expression values across all the IMF2s for 9 randomly selected s-traits.
https://doi.org/10.7554/eLife.03913.009
Figure 2 with 5 supplements
Expression correlations between s-traits originating from the same genes.
(A) The distribution of expression correlations between s-traits originating from the same genes. (B) The distribution of expression correlations between s-traits originating from LOC_Os03g01360. (C)…
Figure 2—figure supplement 1
Expression correlations between s-traits of different sizes originating from the same genes.
The s-traits originated from same genes were grouped by their sizes. The expression correlations between s-traits from pairs of groups were calculated, and results with more than 2000 correlations …
Figure 2—figure supplement 2
Expression correlations between s-traits originating from different regions of the same mother genes.
(A) The distribution of expression correlations between s-traits originating from introns of the same mother genes. (B) The distribution of expression correlations between s-traits originating from …
Figure 2—figure supplement 3
The expression correlations between s-traits originating from LOC_Os03g01360.
The structure of LOC_Os03g01360 is shown by the bars in the left and bottom of the plot. Each point represents the expression correlation between a pair of s-traits.
Figure 2—figure supplement 4
Distribution of sRNA cluster sizes.
(A) Distribution of the sizes of all 80,362 sRNA clusters. (B) Distribution of the sizes of sRNA clusters not longer than 10 kb. (C) Distribution of the sizes of sRNA clusters not longer than 1 kb. …
Figure 2—figure supplement 5
Expression correlations between s-traits originating from the same sRNA clusters.
(A) The distribution of expression correlations between s-traits originating from the same sRNA clusters. (B) The distribution of expression correlations between s-traits originating from …
Figure 3 with 2 supplements
Correlations between the s-traits and the transcripts of their mother genes.
(A) Correlations between s-traits and e-traits derived from LOC_Os03g01360. (B) Correlations between s-traits and e-traits derived from LOC_Os07g01240. (C) Correlations between s-traits and e-traits …
Figure 3—figure supplement 1
Processing of mRNA sequencing data from raw sequencing reads to genome-mapped reads.
The read number of all IMF2s in each step is shown as box plot.
Figure 3—figure supplement 2
The distribution of mRNA sequencing reads in different genomic regions.
https://doi.org/10.7554/eLife.03913.018
Figure 4 with 8 supplements
sQTLs for the 66,649 s-traits.
The color key shows the LOD value. X-axis, the physical position of sQTLs along the genome of 12 chromosomes. Y-axis, the physical position of s-traits. QTLs with LOD value <5 are not included in …
Figure 4—figure supplement 1
The distribution of the sizes of 1568 bins of the genetic map.
https://doi.org/10.7554/eLife.03913.020
Figure 4—figure supplement 2
The genotypes of 98 IMF2s across the 1568 bins.
Each row represents an IMF2, while each column represents a bin. Different color indicates different genotype. The 12 chromosomes were separated by black vertical lines.
Figure 4—figure supplement 3
The LOD values and expression variations explained by sQTLs (R2) for local- and distant-sQTLs.
(A) The LOD values for local-sQTLs. (B) The LOD values for distant-sQTLs. (C) Expression variation explained by local-sQTL. (D) Expression variation explained by distant-sQTL.
Figure 4—figure supplement 4
scQTLs for the 20,049 sc-traits.
The color key shows the LOD value. X-axis shows the physical position of scQTLs along the genome. Y-axis shows the physical position of sc-traits. QTLs with LOD value <5 are not included in the …
Figure 4—figure supplement 5
The LOD values and expression variations explained by scQTLs (R2) for local- and distant-scQTLs.
(A) The LOD values for local-scQTLs. (B) The LOD values for distant-scQTLs. (C) Expression variation explained by local-scQTL. (D) Expression variation explained by distant-scQTL.
Figure 4—figure supplement 6
eQTLs for the 6123 e-traits.
The color key shows the LOD value. X-axis shows the physical position of eQTLs along the genome. Y-axis shows the physical position of e-traits. QTLs with LOD value <5 are not included in the …
Figure 4—figure supplement 7
The LOD values and expression variations explained by eQTLs (R2) for local- and distant-eQTLs.
(A) The LOD values for local-eQTLs. (B) The LOD values for distant-eQTLs. (C) Expression variation explained by local-eQTL. (D) Expression variation explained by distant-eQTL.
Figure 4—figure supplement 8
Examples of s-traits expressed in only one of the two parents (Zhenshan 97).
The expression value (reads per millions [RPM]) of each s-trait in the IMF2 population is shown as box plot. IMF2s with different genotypes as determined by the genotype of the bin harboring these …
Figure 5 with 2 supplements
The distribution of traits and QTLs in the 1568 bins.
(A) s-traits and sQTLs. (B) sc-traits and scQTLs. (C) e-traits and e-QTLs. The 1568 bins are arranged from left to right according to their genomic position. The width of bar represents the size of …
Figure 5—figure supplement 1
DNA sequence polymorphisms of sRNA biogenesis genes between Zhenshan 97 and Minghui 63.
(A) SNPs in sRNA biogenesis genes between Zhenshan 97 and Minghui 63. (B) Indels in sRNA biogenesis genes between Zhenshan 97 and Minghui 63.
Figure 5—figure supplement 2
Expression correlations between sRNA biogenesis genes in distant-sQTL hotspots and s-traits regulated by these distant-sQTLs.
(A) Expression correlations between s-traits regulated by Bin632–Bin637 and the e-traits in these bins. (B) Expression correlations between s-traits regulated by Bin632–Bin637 and the e-traits of OsR…
Figure 6
Genetic co-regulation between s-traits originating from the same mother genes.
(A) Expression correlations between pairs of s-traits originating from the same mother genes and both regulated by local-sQTLs. (B) Expression correlations between pairs of s-traits originating from …
Figure 7
Genetic co-regulation between s-traits and their mother genes.
(A) Expression correlations between s-traits and their mother genes that were both controlled by local-QTLs. (B) Expression correlations between s-traits regulated by distant-sQTLs and their mother …
Figure 8 with 1 supplement
The expression levels of sRNA, mRNA, and the DNA methylation levels of the hybrid and its parents in the region of LOC_Os03g01360 and LOC_Os07g01240.
(A) The structure of LOC_Os03g01360 is shown in the top panel, with black color representing exons and red color indicating introns and UTRs. The number of mRNA reads is shown in blue color. The …
Figure 8—figure supplement 1
The expression level (fragments per kb per million [FPKM]) of the three transcripts of LOC_Os03g01360 (A) and the two transcripts of LOC_Os07g01240 (B) in Zhenshan 97, Minghui 63, and the hybrid.
https://doi.org/10.7554/eLife.03913.034Tables
Table 1
Average number of reads obtained for the 104 libraries from flag leaves of the IMF2 population and the parental lines
| Quality filtering* | Size filtering† | NcRNA filtering‡ | Genome mapping§ | Total# | |
|---|---|---|---|---|---|
| Redundant reads | 22,864,067 | 22,318,928 | 19,968,566 | 9,816,933 | 1,020,961,021 |
| Distinct reads | 5,106,654 | 5,035,073 | 4,986,184 | 2,499,233 | 53,613,739 |
-
*
Reads after filtering out low-quality reads.
-
†
Reads of 18–26 nt in length.
-
‡
Reads after eliminating ones matching tRNAs, rRNAs, snRNAs, and snoRNAs.
-
§
Reads mapped to the SNP-replaced reference genomes of the parents with unique locations allowing no mismatch.
-
#
The total reads of 104 libraries used to identify s-traits.
