Understanding how regulatory sequences control gene expression is fundamental to explain how phenotypes arise in health and disease. The functions of regulatory elements must ultimately be understood within their genomic environment and developmental- or tissue-specific contexts. Here, we used induced Cas9 expression and multiplexed guide RNAs in C. elegans to create hundreds of mutations in enhancers, promoters and 3′ UTRs of 16 genes in parallel. We then analyzed the resulting complex populations by either selecting for phenotypic traits or reporter expression, or by DNA- or RNA amplicon sequencing of bulk samples. We developed a software pipeline, crispr-DART, to analyze targeted sequencing and describe the characteristics of >12,000 dsDNA break-induced indel mutations. We also analyzed the
lin-41 3′ UTR and found that the two
let-7 miRNA binding sites can function independently and that one of the sites is more important for mRNA repression and phenotype. In summary, our approach enables highly parallelized functional analysis of regulatory sequences in vivo.