Cells in multicellular organisms are genetically heterogeneous owing to somatic mutations. The accumulation of somatic genetic variation in species undergoing asexual (or clonal) reproduction (termed modular species) may lead to phenotypic heterogeneity among modules. However, abundance and dynamics of somatic genetic variation under clonal growth, a widespread life history in nature, remain poorly understood. Here we show that branching events in a seagrass clone or genet leads to population bottlenecks at the cellular level and hence the evolution of genetically differentiated modules. Studying inter-module somatic genetic variation, we uncovered thousands of SNPs that segregated among modules. The strength of purifying selection on mosaic genetic variation was greater at the intra-module comparing with the inter-module level. Our study provides evidence for the operation of selection at multiple levels, of cell population and modules. Somatic genetic drift leads to the emergence of genetically unique modules; hence, modules in long-lived clonal species constitute an appropriate elementary level of selection and individuality.All multicellular species, from plants to humans, are genetic mosaics, owing to mitotic errors (somatic mutations) during growth and development 1,2 . In unitary species, the resulting intra-organismal genetic heterogeneity may lead to genomic conflict and is often associated with degenerative disease such as cancer 3 (but see ref 4 ). Somatic genetic variation may play a different, more positive role in species undergoing asexual (or clonal) reproduction, hereafter called modular species, featured by 65% of all plant families 5 and 35% of all animal phyla 6 . Modular species have a simple body plan, and often indeterminate growth, during which iterative units (modules) emerge by asexual proliferation through fission, budding or branching 7 . Modules originating from the same zygote collectively form the clone or genet 8 . Under modular organization, somatic genetic variation may segregate when new modules are and the base of the leaves was sequenced with an average coverage of 81x per module ( Supplementary Tables 2 & 3). Sequenced reads were mapped against the Z. marina reference genome 24 with an average mapping rate of 94%, noting that the reference genome was from a clonal meadow only 22 km from the samples used in the present study ( Supplementary Fig. 2). We focused first on fixed genetic variants that were present in all the cells within the module. A total of 38,831 fully heterozygous SNPs were detected after stringent filtering ( Supplementary Fig. 2) that had a different allele compared to the reference genome.The resulting multi-locus SNP genotypes were first used to assess whether or not all 24 sampled modules were members of the same genet, i.e. originated exclusively by vegetative propagation. Members of the same genet will harbor an excess of shared heterozygous loci because under mitosis, all SNPs of the ancestral zygote will remain heterozygous in offspring. The majori...