- Brief report
- Open access
- Published:
PeakMatcher facilitates updated Aedes aegypti embryonic cis-regulatory element map
Hereditas volume 158, Article number: 7 (2021)
Abstract
Background
The Aedes aegypti mosquito is a threat to human health across the globe. The A. aegypti genome was recently re-sequenced and re-assembled. Due to a combination of long-read PacBio and Hi-C sequencing, the AaegL5 assembly is chromosome complete and significantly improves the assembly in key areas such as the M/m sex-determining locus. Release of the updated genome assembly has precipitated the need to reprocess historical functional genomic data sets, including cis-regulatory element (CRE) maps that had previously been generated for A. aegypti.
Results
We re-processed and re-analyzed the A. aegypti whole embryo FAIRE seq data to create an updated embryonic CRE map for the AaegL5 genome. We validated that the new CRE map recapitulates key features of the original AaegL3 CRE map. Further, we built on the improved assembly in the M/m locus to analyze overlaps of open chromatin regions with genes. To support the validation, we created a new method (PeakMatcher) for matching peaks from the same experimental data set across genome assemblies.
Conclusion
Use of PeakMatcher software, which is available publicly under an open-source license, facilitated the release of an updated and validated CRE map, which is available through the NIH GEO. These findings demonstrate that PeakMatcher software will be a useful resource for validation and transferring of previous annotations to updated genome assemblies.
Introduction
The Aedes aegypti mosquito, vector of the viruses responsible for the yellow, dengue, chikungunya, and Zika fevers, is a significant threat to global human health. A. aegypti is widespread throughout Africa, the Americas, and Asia, threatening a large fraction of human populations [5].
Although originally sequenced in 2007, a chromosome-complete assembly of the A. aegypti genome was produced only recently. Matthews, et al. [18] combined long-read PacBio and chromosome conformation capture (Hi-C) sequencing technologies using a new de novo assembly method recently introduced by Dudchenko, et al. [7] to generate a chromosome-complete assembly (AaegL5) for A. aegypti. The updated assembly substantially reduced sequence (by ~ 100 mb) and gene duplication (1463 genes previously annotated as paralogs were collapsed). Accuracy and completeness of gene models were increased, resulting in 915 additional genes with > 80% coverage when compared to the respective orthologs of these genes in Drosophila melanogaster and the identification of additional chemosensory receptors. The new assembly also improved the fidelity of the M/m locus which will be important to unraveling the underlying mechanisms of sex determination in A. aegypti.
Cis-regulatory elements (CREs) and the tissue- and developmental-specific chromatin accessibility patterns associated with these regions control gene expression. Mapping CREs and monitoring differences in the state of the chromatin accessibility to these elements is critical to gaining a full and complete understanding of gene expression. For example, differences in chromatin accessibility are associated with Plasmodium falciparum infection [27] and immune response [23] in the malaria vector Anopheles gambiae, developmental stage in D. melanogaster [20], regulation of silk protein genes in Bombyx mori [33, 34].
In 2016, Behura et al. [3] performed FAIRE-Seq (Formaldehyde-Assisted Identification of Regulatory Elements combined with next-generation sequencing) [9] on whole embryos to map CREs in the A. aegypti genome. FAIRE-Seq identifies open chromatin regions enriched with CREs [9, 13, 20, 29]. Compared with other techniques, FAIRE-Seq requires relatively small amounts of raw genetic material, demonstrates low technical variability, and is associated with a relatively straightforward experimental protocol [19, 28, 30, 31]. A subset of FAIRE peaks was confirmed in vivo to demonstrate enhancer-like activity in D. melanogaster reporter assays [3, 21]. Subsequently, the FAIRE data were used to identify CREs associated with expression of olfactory receptor neurons in A. aegypti [21]. Despite the value of the FAIRE-seq data set and the AaegL5 genome assembly to the mosquito research community, the embryonic CRE map which had been generated could not be viewed without reprocessing the FAIRE-seq data set for the updated A. aegypti genome assembly.
We reprocessed the raw FAIRE sequencing data to create a de novo annotation of CREs in the AaegL5 assembly. As part of that effort, we developed and open-sourced a new tool, PeakMatcher, to match DNA enrichment assay peaks called from the same sequencing data across different genome assembly versions. We applied PeakMatcher to create a list of corresponding peaks from the AaegL3.4 and AaegL5 CRE maps. Using the peak mapping, we confirmed that 14 of 16 experimentally validated peaks in Behura, et al. [3] and [21] were reconstructed in the latest genome assembly (AaegL5, [18]). Our updated and validated CRE map is publicly available through the NIH GEO (GSE162150). We anticipate that PeakMatcher will be useful for validating DNA enrichment assay (e.g., ChIP-Seq, DNase-Seq, or STARR-Seq) annotations across two genome assemblies.
Results and discussion
FAIRE peaks properties consistent across Aedes aegypti genome versions
We called 128,307 FAIRE peaks for AaegL5 [18] and compared the peaks against the previously called peaks for AaegL3 [22] from Behura et al. [3]. Peaks were called using the same parameters (extent size of 550, no shifting, p-value < 0.01 with no FDR control, and the same estimated mappable genome size used in [3]) to reduce bias in the comparison. 23.1, 41.2, and 34.6% of the AaegL5 peaks were located on chromosomes 1, 2, and 3 respectively, while the remaining 1.1% peaks were located on non-chromosome contigs.
Summary statistics of FAIRE peaks were largely in agreement for both of the genome versions considered (see Table 1). The distributions of FAIRE peaks across AaegL5 and AaegL3 were consistent (see Fig. 1) even though AaegL5 is 106 Mbp smaller, in part because of reduced duplication [18]. We were able to call an additional 6712 FAIRE peaks for AaegL5 that covered an additional 4.8 Mbp of the genome. FAIRE peaks covered an additional 1.2 percentage points (pp) of the fraction of the total genome (see Fig. 1). Increased coverage was observed primarily in exons (5.6 pp) and TSSs (7.9 pp), while coverage of introns (0.9 pp) and intergenic regions (0.3 pp) did not change substantially (see Fig. 1). In turn, enrichments of exons and TSSs were higher in AaegL5 (> 2.5X and > 3.4, respectively; see Fig. 1). FAIRE enrichment of TSSs is consistent with observations in other species [9, 29].
PeakMatcher method
We created a new method and software package called PeakMatcher for matching peaks from the same experiment across different genome assemblies. PeakMatcher is targeted towards cases where the same DNA enrichment assay data is being compared across two assembly versions. BEDTools [25, 26] and similar tools focus on analysis of overlaps of coordinate intervals, but this is not appropriate when two different genome assembly versions have very different coordinate systems. PeakMatcher uses aligned sequencing reads that are common to peaks in both genomes to match the peaks. It is an alternative to approaches based on whole-genome alignments (e.g., using MUMmer [17] or liftOver [12]).
We initially attempted to match peaks across the two A. aegypti genome assemblies based on whole-genome alignments, but we were only able to successfully find matches for ~ 40% of the peaks (data not shown). As we did not have access to the original AaegL3 alignments, we re-generated the alignments and re-processed the peaks. We validated the re-processed AaegL3 peaks against the original peak list distributed by Behura et al. [3]. We then matched peaks across the two genome versions using our PeakMatcher method. PeakMatcher was able to match 73.7% of the 124,959 re-processed AaegL3 peaks to 78.9% of the 128,307 AaegL5 peaks. The distributions of matched and all AaegL5 peaks were similar; 22.7, 40.5, and 35.8% of the matched peaks were located on chromosomes 1, 2, and 3 respectively, while the remaining 1.1% peaks were located on non-chromosome contigs.
Methods for combining long read and/or chromosome capture sequencing with next-generation sequencing can generate vastly improved and chromosome-complete assemblies. Fueled by these recent leaps in capability, we anticipate that a number of genomes will be re-sequenced and re-assembled. Functional genomics data sets will need to be reprocessed for use with the reassembled genomes. We anticipate that PeakMatcher will be a useful resource for transfer and validation of previous annotations to updated assemblies.
Experimentally validated peaks matched across genome versions
Sixteen AaegL3 FAIRE peaks were previously demonstrated to act as enhancers in transgenic insect reporter assays [3, 21]; 14 of these 16 AaegL3 peaks were matched to peaks in AaegL5 by PeakMatcher (see Table 2). Adjacent genes were consistent for 12 of the 14 matched peaks; inconsistencies for at least two of the three remaining peaks were clearly attributable to differences in the assemblies considered here (data not shown).
FAIRE peak coverage of TSSs associates with gene expression
FAIRE peaks indicate nucleosome-depleted regions associated with regulatory activity. We compared the overlap of FAIRE peaks with 500 bp windows upstream (downstream for negative strand) of transcription start sites with genes indicated to be expressed by RNA-Seq [2]. TSSs for 6611 (65.5%) of the 10,089 genes with increased and 1484 (32.4%) of the 4587 genes with decreased differential expression were overlapped by FAIRE peaks. The differences in frequencies are statistically significant (p-value < 10− 100, χ2 test of independence). We concluded that the FAIRE peaks successfully identify active promoters, consistent with observations in other species [9].
M/m locus
For the first time, a high-quality assembly of the M/m locus is available through the AaegL5 genome of A. aegypti, in which a dominant M locus establishes the male sex (male genotype = M/m; female genotype = m/m [18];). We analyzed the FAIRE peaks in a ~ 1.5 Mb region (151.68–152.95 Mb) in the predicted M/m locus (see Table 3, [18]). Peaks were associated with the TSSs (within 2.5 kb upstream) of five long non-coding RNA (lncRNA) genes: AAEL020975, AAEL022711, AAEL024704, AAEL025015, and AAEL026346. FAIRE elements were also associated with two protein-coding genes, Nix and myo-sex. Nix, which encodes a male-determining factor that is necessary and sufficient to drive male-specific development in A. aegypti [11], overalapped with FAIRE peaks located at the TSS and within an intron of this gene. The myo-sex gene, which is required for A. aegypti male flight [1, 10], is associated with eight intronic FAIRE peaks. FAIRE elements associated with these M/m locus genes (Table 3) may function as CREs that regulate sex-specific gene expression in A. aegypti.
FAIRE-seq mapping of open chromatin associated with the M/m locus represents a first glimpse at regions that could potentially have different chromatin signatures in male and female mosquitoes. The FAIRE DNA sequenced by Behura et al. [3] was prepared from a combination of male and female embryos. It would be interesting to pursue separate FAIRE-seq analyses in male vs. female mosquitoes, particularly once complete sequence information is available for the entire M and m regions. This would permit a more detailed understanding of both genetic, as well as epigenetic, differences between the two sexes. It is also interesting to note the many FAIRE elements associated with lncRNA genes residing at the M/m locus (Table 3). Although the roles of these genes have not yet been described in A. aegypti, our detection of open chromatin within and flanking these genes suggests that lncRNAs could play critical roles in the regulation of sex-specific development and differentiation.
Conclusion
Next-generation sequencing (NGS) drove a massive increase in genome sequencing by lowering costs and increasing per-base calling quality. Unfortunately, many of the resulting assemblies were fragmented. NGS has recently been combined with long-read and chromosome conformation capture (e.g., Hi-C) sequencing to produce new, vastly improved genome assemblies [7]. Resequencing and reassembly of A. aegypti demonstrated the power of these techniques, and many genomes will likely be re-sequenced and re-assembled in the next few years, particularly through efforts such as the Earth BioGenome project [14]. The approaches developed here will be useful to other researchers facing similar needs to reprocess functional genomics data for these updated genome assemblies.
We updated whole-embryo A. aegypti FAIRE sequencing data for the new AaegL5 genome assembly and validated the reprocessed data to ensure consistency with the previous FAIRE peaks. The resulting FAIRE peak lists are available to the public through the NIH GEO (GSE162150), where they will continue to be available to the insect vector community. Additionally, we anticipate that PeakMatcher will be useful to other researchers who are validating re-annotated DNA enrichment assay (e.g., ChIP-Seq, DNase-Seq, or STARR-Seq) data on updated genome assemblies.
Methods
FAIRE peak processing
FAIRE-Sequencing of A. aegypti whole embryos was originally performed by Behura, et al. [3]. Raw FAIRE-Seq sequencing data were downloaded from the NIH SRA (SRR2530418, SRR25304189, and SRR25304120). The A. aegypti AaegL3.5 [22] and AaegL5.2 [18] genomes and gene sets retrieved from Vectorbase [8]. 382,611,452, 341,380,350, and 296,902,904 sequencing reads were cleaned with trimmomatic (ILLUMINACLIP:TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36) [4]; 340,532,066, 305,899,612, and 274,957,716 paired reads survived cleaning.
The paired reads were aligned to the two reference genomes using BWA backtrack [15], and filtered with SAMtools [16] to remove reads with mapping qualities < 10, unmapped reads, and secondary alignments of reads. 179,608,283, 152,618,302, and 148,862,656 (total of 52.2%) and 171,583,885, 145,715,076, and 142,034,403 (total of 49.9%) of the reads aligned to AaegL3 and AaegL5, respectively, survived filtering. 92.6 and 92.5% of the AaegL3- and AaegL5-aligned reads were properly paired. 57.0 and 56.9% of the AaegL3- and AaegL5-aligned reads were identified as duplicates by Picard tools [24]. (Note that we did not subsequently filter out duplicate reads since MACS already does so as part of its workflow.) Average read lengths were similar across the two sets of alignments; the average lengths were 97, 96, and 97 bp, respectively, for each biological replicate. Insert sizes of 134.6 ± 16.6, 141.2 ± 19.0, and 135.8 ± 16.9 (AaegL3) and 135.0 ± 17.3, 240.5 ± 873.2, and 136.2 ± 17.5 (AaegL5) bp were observed.
Peaks were called from pooled alignments of the three biological replicates using MACS2 [32] (−t rep1.bam rep2.bam rep3.bam --nomodel -p 0.01 --extsize 550 -g 1.24e9). No controls were used. We used the default single-ended mode (BAM) which uses the 5′ end tag, disabled the shifting model (−-nomodel), used a p-value cutoff of 0.01 (in place of false discovery rate control), a shift size of 550 bp (−-extsize 500 bp), and an estimated mappable genome size of 1240 Mb to match the parameters used in [3] and facilitate comparison.
Peak-genomic element overlap enrichment analysis
Gene and exon coordinates were extracted from GFF3 files. Intronic regions were calculated by subtracting exon intervals from gene intervals using BEDtools [25, 26] subtract (default parameters). TSSs were identified as 500 bp windows upstream (positive strand) or downstream (negative strand) of genes using custom Python scripts.
Peak and element intervals were sorted by chromosome and then starting position using “sort -k 1,1 -k2,2n.” Overlaps between FAIRE peaks and genomic elements were determined with “bedtools intersect -a elements.bed -b peaks.bed -f 0.1”. Enrichment in peak-element overlaps was calculated using the following equations:
fraction of genome = total element coverage (bp)/genome size (bp)
Matching peaks across genome versions
The PeakMatcher method identifies corresponding peaks called for two different genome assembly versions from the same DNA enrichment assay sequencing data. One genome (usually the older genome which has been previously validated) is defined as the source genome, and the second genome (usually the new genome to be validated) is defined to be the target genome.
PeakMatcher operates in two steps (see Fig. 2). First, peaks from the source genome are associated with aligned reads. PeakMatcher takes the peak lists (e.g., generated by MACS2) and a source genome BAM file as inputs. Each peak is represented as a tuple of four elements (chromosome, start position, end position, peak and name) and an interval tree of peaks is constructed for each molecule (e.g., chromosome, chromosome arm, or scaffold). Interval trees are a specialized type of search tree data structure that enable fast O (log N) overlap queries (e.g., return all intervals in the tree that overlap a given query interval) [6]. Every read in the BAM file is represented as a tuple of three elements (chromosome, start position, and end position) and used to query the tree for overlaps with the called peaks. A list of read-peak names pairs is output.
In the second step, peaks from the source genome are associated with corresponding peaks from the target genome. In the same fashion as the source genome, an interval tree is constructed from the target genome peak list, and the tree is queried to find peaks overlapping each read. A list of target peak-read pairs is generated, and the original list of source peak-read pairs is read in. The two lists are joined on the read names to generate a final list of “matched” peaks as peak-peak pairs.
AaegL3 FAIRE peaks were matched to AaegL5 FAIRE peaks using this method. Sixteen experimentally validated AaegL3 FAIRE peaks listed in Table 2 of Behura, et al. [3] and Table 1 of Mysore, et al. [21] were identified for further validation. The Vectorbase genome browser was used to manually inspect the regions around each peak pair to confirm the match by verifying that neighboring genes in AegL3 for each peak were present in AaegL5.
RNA-Seq analysis
RNA-Seq data were retrieved from Akbari, et al. [2]. Genes were partitioned into high and low differential expression groups using a FPKM threshold of 1. TSSs for each gene present were extracted as described above. Overlaps between TSSs and peaks were determined using the BEDtools intersect command: “bedtools intersect -a promoters.bed -b peaks.bed -wa -F 0.1 -u”. Statistical differences in promoter overlap counts were confirmed with a Chi-squared test.
M/m locus
We followed the procedure described in the genomic element overlap section above to analyze the M/m locus. We restricted our analysis to FAIRE peaks, introns, and TSSs overlapping the M/m locus region (1:151.68–152.95 Mb).
Availability of data and materials
The original FAIRE sequencing data from Behura, et al. [3] is available from the NIH SRA (SRR2530418, SRR25304189, and SRR25304120). We made our re-processed AaegL3 and AaegL5 FAIRE peak lists available through the NIH GEO (GSE162150). The PeakMatcher software is hosted under a separate GitHub repository (https://github.com/rnowling/peak-matcher). The software is licensed under the open-source Apache Software License v2.
Abbreviations
- AaegL3:
-
A. aegypti genome assembly version 3
- AaegL5:
-
A. aegypti genome assembly version 5
- bp:
-
Base pair
- CRE:
-
cis-regulatory element
- FAIRE-Seq:
-
Formaldehyde-Assisted Identification of Regulatory Elements combined with next-generation sequencing
- Hi-C:
-
High-throughput chromosome conformation sequencing
- lncRNA:
-
Long non-coding RNA
- pp:
-
Percentage point
- RNA-Seq:
-
RNA sequencing
- TSS:
-
Transcription start site
References
Aryan A, Anderson MAE, Biedler JK, Qia Y, Overcash JM, Naumenko AN, et al. Nix alone is sufficient to convert female Aedes aegypti into fertile males and myo-sex is needed for male flight. Proc Natl Acad Sci. 2020;117(30):17702–9.
Akbari OS, Antoshechkin I, Amrhein H, Williams B, Diloreto R, Sandler J, et al. The developmental Transcriptome of the mosquito Aedes aegypti, an invasive species and major Arbovirus vector. G3. 2013;3(9):1493–509.
Behura SK, Sarro J, Li P, Mysore K, Severson DW, Emrich SJ, et al. High-throughput Cis-regulatory element discovery in the vector mosquito Aedes aegypti. BMC Genomics. 2016;17(May):341.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Centers for Disease Control, htpp://cdc.gov. Accessed Nov 2020.
Cormen TH, Leiserson CE, Rivest RL, Stein C. Introduction to algorithms. 3rd ed. Cambridge: MIT Press; 2009.
Dudchenko O, Batra SS, Omer AD, Nyquist SK, Hoeger M, Durand NC, et al. De novo assembly of the Aedes aegypti genome using hi-C yields chromosome-length scaffolds. Science. 2017;356(6333):92–5.
Giraldo-Calderón GI, Emrich SJ, MacCallum RM, Maslen G, Dialynas E, Topalis P, et al. VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic Acids Res. 2015;43(Database issue):D707–13.
Giresi PG, Kim J, McDaniell RM, Iyer VR, Lieb JD. FAIRE (formaldehyde-assisted isolation of regulatory elements) isolates active regulatory elements from human chromatin. Genome Res. 2007;17(6):877–85.
Hall AB, Timoshevskiy VA, Sharakhova MV, Jiang X, Basu S, Anderson MAE, et al. Insights into the preservation of the Homomorphic sex-determining chromosome of Aedes aegypti from the discovery of a male-biased gene tightly linked to the M-locus. Genome Biol Evol. 2014;6:179–91.
Hall AB, Basu S, Jiang X, Qi Y, Timoshevskiy VA, Biedler JK, et al. SEX DETERMINATION. A male-determining factor in the mosquito Aedes aegypti. Science. 2015;348(6240):1268–70.
Kuhn RM, Haussler D, Kent WJ. The UCSC genome browser and associated tools. Brief Bioinform. 2013;14(2):144–61.
Lai Y-T, Deem KD, Borràs-Castells F, Sambrani N, Rudolf H, Suryamohan K, et al. Enhancer identification and activity evaluation in the red flour beetle, Tribolium castaneum. Development. 2018;145(7). https://doi.org/10.1242/dev.160663.
Lewin HA, Robinson GE, John Kress W, Baker WJ, Coddington J, Crandall KA, et al. Earth BioGenome project: sequencing life for the future of life. Proc Natl Acad Sci U S A. 2018;115(17):4325–33.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Marçais G, Delcher AL, Phillippy AM, Coston R, Salzberg SL, Zimin A. MUMmer4: a fast and versatile genome alignment system. PLoS Comput Biol. 2018;14(1):e1005944.
Matthews BJ, Dudchenko O, Kingan SB, Koren S, Antoshechkin I, Crawford JE, et al. Improved reference genome of Aedes aegypti informs Arbovirus vector control. Nature. 2018;563(7732):501–7.
McKay DJ. Using formaldehyde-assisted isolation of regulatory elements (FAIRE) to identify functional regulatory DNA in insect genomes. In: Brown SJ, Pfrender ME, editors. Insect genomics: methods and protocols. New York: Springer New York; 2019. p. 89–97.
McKay DJ, Lieb JD. A common set of DNA regulatory elements shapes Drosophila appendages. Dev Cell. 2013;27(3):306–18.
Mysore K, Li P, Duman-Scheel M. Identification of Aedes aegypti Cis-regulatory elements that promote gene expression in olfactory receptor neurons of distantly related dipteran insects. Parasit Vectors. 2018;11(1):406.
Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Zhijian Jake T, et al. Genome sequence of Aedes aegypti, a major Arbovirus vector. Science. 2007;316(5832):1718–23.
Pérez-Zamorano B, Rosas-Madrigal S, Lozano OAM, Méndez MC, Valverde-Garduño V. Identification of Cis-regulatory sequences reveals potential participation of Lola and Deaf1 transcription factors in Anopheles gambiae innate immune response. PLoS One. 2017;12(10):e0186435.
Broad Institute. 2018. “Picard Tools.” 2018. http://broadinstitute.github.io/picard/.
Quinlan, Aaron R. 2014. “BEDTools: The Swiss-Army Tool for Genome Feature Analysis.” Current Protocols in Bioinformatics / Editoral Board, Andreas D. Baxevanis... [et Al.] 47 (September): 11.12.1–34.
Quinlan AR, Hall IM. BEDTools: a flexible suite of Utilities for Comparing Genomic Features. Bioinformatics. 2010;26(6):841–2.
Ruiz JL, Yerbanga RS, Lefèvre T, Ouedraogo JB, Corces VG, Gómez-Díaz E. Chromatin changes in Anopheles gambiae induced by plasmodium falciparum infection. Epigenetics Chromatin. 2019;12(1):5.
Simon JM, Giresi PG, Davis IJ, Lieb JD. Using formaldehyde-assisted isolation of regulatory elements (FAIRE) to isolate active regulatory DNA. Nat Protoc. 2012;7(2):256–67.
Song L, Zhang Z, Grasfeder LL, Boyle AP, Giresi PG, Lee B-K, et al. Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity. Genome Res. 2011;21(10):1757–67.
Sun Y, Miao N, Sun T. Detect accessible chromatin using ATAC-sequencing, from principle to applications. Hereditas. 2019;156(August):29.
Tsompana M, Buck MJ. Chromatin accessibility: a window into the genome. Epigenetics Chromatin. 2014;7(1):33.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.
Zhang Q, Cheng T, Jin S, Guo Y, Wu Y, Liu D, et al. Genome-wide open chromatin regions and their effects on the regulation of silk protein genes in Bombyx mori. Sci Rep. 2017;7(1):12919.
Zhang Q, Cheng T, Sun Y, Wang Y, Feng T, Li X, et al. Synergism of open chromatin regions involved in regulating genes in Bombyx mori. Insect Biochem Mol Biol. 2019;110(July):10–8.
Acknowledgements
We would like to thank Joseph Sarro, Keshava Mysore, Ping Li, and David W. Severson for their contributions to the initial acquisition and analyses of the A. aegypti FAIRE-Seq data and Tao Liu for advice on the parameters for MACS2.
Funding
This material is based upon work supported by the National Science Foundation under Grant No. IIS-1947257 to RJN and NIH/NIAID Award R21 AI117145–01 to MDS.
Author information
Authors and Affiliations
Contributions
RJN contributed to the analysis and interpretation of the data, creation of software, writing and editing of the manuscript. MDS contributed to the concept and design of the study, acquisition, analysis and interpretation of the data, and writing and editing of the manuscript. SKB contributed to the acquisition, analysis, and interpretation of the data. SJE and MSH contributed to the analysis and interpretation of the data and editing of the manuscript. The authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
Molly Duman-Scheel was listed as an author on U.S. Patent Office Application No. 62/751,052, “Sex-linked RNAi Insecticide Materials and Methods.” This application did not impact study design or interpretation of the data. All other authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Nowling, R.J., Behura, S.K., Halfon, M.S. et al. PeakMatcher facilitates updated Aedes aegypti embryonic cis-regulatory element map. Hereditas 158, 7 (2021). https://doi.org/10.1186/s41065-021-00172-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s41065-021-00172-2