Discovering single nucleotide variants and indels from bulk and single-cell ATAC-seq

Genetic variants and de novo mutations in regulatory regions of the genome are typically discovered by whole-genome sequencing (WGS), however WGS is expensive and most WGS reads come from non-regulatory regions. The Assay for Transposase-Accessible Chromatin (ATAC-seq) generates reads from regulatory sequences and could potentially be used as a low-cost 'capture' method for regulatory variant discovery, but its use for this purpose has not been systematically evaluated. Here we apply seven variant callers to bulk and single-cell ATAC-seq data and evaluate their ability to identify single nucleotide variants (SNVs) and insertions/deletions (indels). In addition, we develop an ensemble classifier, VarCA, which combines features from individual variant callers to predict variants. The Genome Analysis Toolkit (GATK) is the best-performing individual caller with precision/recall on a bulk ATAC test dataset of 0.92/0.97 for SNVs and 0.87/0.82 for indels within ATAC-seq peak regions with at least 10 reads. On bulk ATAC-seq reads, VarCA achieves superior performance with precision/recall of 0.99/0.95 for SNVs and 0.93/0.80 for indels. On single-cell ATAC-seq reads, VarCA attains precision/recall of 0.98/0.94 for SNVs and 0.82/0.82 for indels. In summary, ATAC-seq reads can be used to accurately discover non-coding regulatory variants in the absence of whole-genome sequencing data and our ensemble method, VarCA, has the best overall performance.

Shoot and root biomass production in semi-arid shrublands exposed to long-term experimental N input

Anthropogenic nitrogen (N) deposition has affected the primary production of terrestrial ecosystems worldwide; however, ecosystem responses often vary over time because of transient responses, interactions between N, precipitation, and/or other nutrients, and changes in plant species composition. Here we report N-induced changes in above- and below-ground standing crop and production over an 11-year period for two semi-arid shrublands, chaparral and coastal sage scrub (CSS), of Southern California. Shrubs were exposed to 50 kgN ha^-1 in the fall of each year to simulate the accumulation of dry N deposition, and shoot and root biomass and leaf area index (LAI) were measured every 3 months to assess how biomass production responded to chronic, dry N inputs. N inputs significantly altered above- and below-ground standing crop, production, and LAI; however, N impacts varied over time. For chaparral, N inputs initially increased root production but suppressed shoot production; however, over time biomass partitioning reversed and plants exposed to N had significantly more shoot biomass. In CSS, N inputs caused aboveground production to increase only during wet years, and this interaction between added N and precipitation was due in part to a highly flexible growth response of CSS shrubs to increases in N and water availability and to a shift from slower-growing native shrubs to fast-growing introduced annuals. Together, these results indicate that long-term N inputs will lead to complex, spatially and temporally variable growth responses for these, and similar, Mediterranean-type shrublands.