Transposons are a major contributor to gene expression variability under selection in rice populations

Reviewer #1 (Public Review):

For many years it has been understood that transposable elements (TEs) are an important source of natural variation. This is because, in addition to simple knockouts of genes, TEs carry regulatory sequences that can, and sometimes do, affect the expression of genes near the TEs. However, because TEs can be difficult to map to reference genomes, they have generally not been used for trait mapping. Instead, single nucleotide polymorphisms are widely used because they are easy to detect when using short reads. However, improvements in sequencing technology, as well as an increased appreciation of the importance of TEs to both linked to favorable alleles and are more likely to be causing the changes that make those alleles beneficial in a given environment. Further, because TE activity can occur after bottlenecks, they can provide polymorphisms in the absence of variation in point mutations.

In this manuscript, the authors carefully examine insertion polymorphisms in rice and demonstrate linkage to differences in expression. To do this, they used expression quantitative trait locus (eQTL) GWAS using TIPs as genetic markers to examine variation in 208 rice accessions. Because they chose to focus on genes that were expressed in at least 10% of the accessions, presumably because more rare variants would end up lacking statistical power. This is an understandable decision, but it says that recent insertions, such as the MITE elements detailed in a previous paper, would not be included. Importantly, although TIPs associated with differentially expressed genes are far less common than SNPs' traditional eQTLs, there were a significant number of eQTLs that showed linkage to TIPs but not to QTL.

The authors then show that of the eQTLs associated with both TIPs and SNPs, TIPs are more tightly linked to the eQTL, and are more likely to be associated with a reduction in expression, with variation in the effects of various TEs families supporting that hypothesis. Here and throughout, however, the distance of the TEs could be an important variable. It is also worth noting the relative numbers in order to assess the claim in the title of the paper. The total number of eQTL-TIPs is ten-fold less than the number of eQTL-SNPs, and, of the eQTLs that have both, there are a significant number of eQTL-TIPs that are not more tightly linked to the expression differences than the eQTL.

The authors show that eQTL-TIPs are more likely to be in the promoter-proximal region, but this may be due to insertion bias, which is well documented in DNA-type elements. Here and throughout the authors are careful to state that the data is consistent with the hypothesis that TEs are the cause of the change, but do not claim that the data demonstrate that they are.

Throughout the rest of the manuscript, the authors systematically build the case for a causal role for TEs by showing, for instance, that eQTL-TIPs show much stronger evidence for selection, with increased expression being more likely to be selected than decreased expression. The authors provide examples of genes most likely to have been affected by TE insertions.

Overall, the authors build a convincing case for TEs being an important source of regulatory information. I don't have any issues with the analysis, but I am concerned about the sweeping claims made in the title. Once you get rid of eQTLs that could be altered by either SNPs or TIPs and include only those insertions that show strong evidence of selection, the number of genes is reduced to only 30. And even in those cases, the observed linkage is just that, not definitive evidence for the involvement of TEs. Although clearly beyond the scope of this analysis, transgenic constructs with the TEs present or removed, or even segregating families, would have been far more convincing.

The fact that many of the eQTL-TIPs were relatively old is interesting because it suggests that selection in domesticated rice was on pre-existing variation rather than new insertions. This may strengthen the argument because those older insertions are less likely to be purged due to negative effects on gene expression. Given that the sequence of these TEs is likely to have diverged from others in the same family, it would have been interesting to see if selection in favor of a regulatory function had caused these particular insertions to move away from more typical examples of the family.

Reviewer #2 (Public Review):

In this manuscript, Castanera et al. investigated how transposable elements (TEs) altered gene expression in rice and how these changes were selected during the domestication of rice. Using GWAS, the authors found many TE polymorphisms in the proximity of genes to be correlated to distinct gene expression patterns between O. sativa ssp. japonica and O. sativa ssp. indica and between two different growing conditions (wet and drought). Thereby, the authors found some evidence of positive selection on some TE polymorphisms that could have contributed to the evolution of the different rice subspecies. These findings are underlined by some examples, which illustrate how changes in the expression of some specific genes could have been advantageous under different conditions. In this work, the authors manage to show that TEs should not be ignored when investigating the domestication of rise as they could have played an important role in contributing to the genetic diversity that was selected. However, this study stops short of identifying causations as the used method, GWAS, can only identify promising correlations. Nevertheless, this study contributes interesting insights into the role TEs played during the evolution of rice and will be of interest to a broader audience interested in the role TEs played during the evolution of plants in general.