Modeling transposable element dynamics with fragmentation equations

In the first case, we assume that full-length TEs replicate at a constant rate regardless of the number of full-length TEs present in the genome. While this assumption simplifies the underlying biological processes, it allows us to derive an explicit analytical form of the time-varying TE density, as well as the steady-state behavior, under both discrete and continuous formulations. Next, we take into account the potential deleterious effects of TEs by modeling TE replication with a logistic growth equation. Under this assumption, the number of actively replicating TEs in a genome is limited by a carrying capacity. While we are still able to derive to derive analytical forms for the time-varying TE density, for both discrete and continuous length formulations, these solutions are implicit. For all four proposed models, we prove existence and uniqueness of these solutions describing TE length distributions. We compare both models and note that the logistic and exponential models initially agree. Since most TEs have not reached carrying capacity, we use the explicit exponential solution to quantify rates of replication to mutations. We apply our model to present day annotated collections of TEs from the genomes of species of fruit-flies, birds, and primates to uncover quantitative relationships of TE dynamics. With the increasing availability of complete genomes, our method is likely to uncover relationships of biological drivers of genomic variation in many species.