Modeling characterization of the vertical and temporal variability of environmental DNA in the mesopelagic ocean

Increasingly, researchers are using innovative methods to census marine life, including identification of environmental DNA (eDNA) left behind by organisms in the water column. However, little is understood about how eDNA is distributed in the ocean, given that organisms are mobile and that physical and biological processes can transport eDNA after release from a host. Particularly in the vast mesopelagic ocean where many species vertically migrate hundreds of meters diurnally, it is important to link the location at which eDNA was shed by a host organism to the location at which eDNA was collected in a water sample. Here, we present a one-dimensional mechanistic model to simulate the eDNA vertical distribution after its release and to compare the impact of key biological and physical parameters on the eDNA vertical and temporal distribution. The modeled vertical eDNA profiles allow us to quantify spatial and temporal variability in eDNA concentration and to identify the most important parameters to consider when interpreting eDNA signals. We find that the vertical displacement by advection, dispersion, and settling has limited influence on the eDNA distribution, and the depth at which eDNA is found is generally within tens of meters of the depth at which the eDNA was originally shed from the organism. Thus, using information about representative vertical migration patterns, eDNA concentration variability can be used to answer ecological questions about migrating organisms such as what depths species can be found in the daytime and nighttime and what percentage of individuals within a species diurnally migrate. These findings are critical both to advance the understanding of the vertical distribution of eDNA in the water column and to link eDNA detection to organism presence in the mesopelagic ocean as well as other aquatic environments.