Enhanced bovine genome annotation through integration of transcriptomics and epi-transcriptomics datasets facilitates genomic biology.

BACKGROUND: The accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues. RESULTS: A total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network. CONCLUSIONS: These validated results show significant improvement over current bovine genome annotations.

Evaluation of indicators of acute emotional states in dogs.

A complete assessment of animal welfare requires not just an understanding of negative emotional states, such as fear and anxiety, but also of positive states, such as calmness and happiness. However, few studies have identified accurate and reliable indicators of positive emotional states in dogs. This study aimed to identify parameters that may serve as indicators of short-term emotional states in dogs. Using a cross-over design, 60 dogs living at a research facility were exposed to six different 10-min scenarios expected to elicit responses varying in emotional valence and arousal. A range of behavioural and physiological parameters were collected and their relationship to anticipated emotional valence and arousal was analysed using linear and logistic mixed models. Cortisol, adrenocorticotropic hormone, heart rate variability, panting, whining, and body shake all demonstrated significant differences based on arousal levels, but only within negative valence scenarios. Scores from a qualitative behavioural assessment (QBA) were associated with both emotional valence and arousal and were considered the best indicator of positive valence. Activity, ear temperature, and sitting were associated with positive high arousal, although this may have been influenced by differing levels of movement induced during these scenarios. Meanwhile, heart rate, secretory immunoglobulin A, standing and lying all showed similar changes associated with arousal for both positive and negative valence scenarios. This study provides a critical first step towards identifying evidence-based indicators of short-term emotional states in dogs, while highlighting considerations that should be made when employing these parameters, including the influence of coder bias, food provision, exercise, and external temperature. Overall, it is recommended future dog emotion and welfare research use a combination of parameters including indicators of both emotional valence and arousal.

Current challenges and future of agricultural genomes to phenomes in the USA.

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

AgAnimalGenomes: browsers for viewing and manually annotating farm animal genomes.

Current genome sequencing technologies have made it possible to generate highly contiguous genome assemblies for non-model animal species. Despite advances in genome assembly methods, there is still room for improvement in the delineation of specific gene features in the genomes. Here we present genome visualization and annotation tools to support seven livestock species (bovine, chicken, goat, horse, pig, sheep, and water buffalo), available in a new resource called AgAnimalGenomes. In addition to supporting the manual refinement of gene models, these browsers provide visualization tracks for hundreds of RNAseq experiments, as well as data generated by the Functional Annotation of Animal Genomes (FAANG) Consortium. For species with predicted gene sets from both Ensembl and RefSeq, the browsers provide special tracks showing the thousands of protein-coding genes that disagree across the two gene sources, serving as a valuable resource to alert researchers to gene model issues that may affect data interpretation. We describe the data and search methods available in the new genome browsers and how to use the provided tools to edit and create new gene models.

Invited review: Reliability computation from the animal model era to the single-step genomic model era.

The calculation of exact reliabilities involving the inversion of mixed model equations poses a heavy computational challenge when the system of equations is large. This has prompted the development of different approximation methods. We give an overview of the various methods and computational approaches in calculating reliability from the era before the animal model to the era of single-step genomic models. The different methods are discussed in terms of modeling, development, and applicability in large dairy cattle populations. The paper also describes the problems faced in reliability computation. Many details dispersed throughout the literature are presented in this paper. It is clear that a universal solution applicable to every model and input data may not be possible, but we point out several efficient and accurate algorithms developed recently for a variety of very large genomic evaluations.

Safely resuming neglected tropical disease control activities during COVID-19: Perspectives from Nigeria and Guinea.

Since its early spread in early 2020, the disease caused by the novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Coronavirus Disease 2019 (COVID-19) has caused mass disruptions to health services. These have included interruptions to programs that aimed to prevent, control, and eliminate neglected tropical diseases (NTDs). In March 2020, the World Health Organization (WHO) released interim guidelines recommending the temporary cessation of mass drug administration (MDA), community-based surveys, and case detection, while encouraging continuation of morbidity management and vector control where possible. Over the course of the following months, national programs and implementing partners contributed to COVID-19 response efforts, while also beginning to plan for resumption of NTD control activities. To understand the challenges, opportunities, and recommendations for maximizing continuity of disease control during public health emergencies, we sought perspectives from Nigeria and Guinea on the process of restarting NTD control efforts during the COVID-19 pandemic. Through semistructured interviews with individuals involved with NTD control at the local and national levels, we identified key themes and common perspectives between the 2 countries, as well as observations that were specific to each. Overall, interviewees stressed the challenges posed by COVID-19 interruptions, particularly with respect to delays to activities and related knock-on impacts, such as drug expiry and prolonged elimination timelines, as well as concerns related to funding. However, respondents in both countries also highlighted the benefits of a formal risk assessment approach, particularly in terms of encouraging information sharing and increasing coordination and advocacy. Recommendations included ensuring greater availability of historical data to allow better monitoring of how future emergencies affect NTD control progress; continuing to use risk assessment approaches in the future; and identifying mechanisms for sharing lessons learned and innovations between countries as a means of advancing postpandemic health systems and disease control capacity strengthening.

Characterizing Genetic Regulatory Elements in Ovine Tissues.

The Ovine Functional Annotation of Animal Genomes (FAANG) project, part of the broader livestock species FAANG initiative, aims to identify and characterize gene regulatory elements in domestic sheep. Regulatory element annotation is essential for identifying genetic variants that affect health and production traits in this important agricultural species, as greater than 90% of variants underlying genetic effects are estimated to lie outside of transcribed regions. Histone modifications that distinguish active or repressed chromatin states, CTCF binding, and DNA methylation were used to characterize regulatory elements in liver, spleen, and cerebellum tissues from four yearling sheep. Chromatin immunoprecipitation with sequencing (ChIP-seq) was performed for H3K4me3, H3K27ac, H3K4me1, H3K27me3, and CTCF. Nine chromatin states including active promoters, active enhancers, poised enhancers, repressed enhancers, and insulators were characterized in each tissue using ChromHMM. Whole-genome bisulfite sequencing (WGBS) was performed to determine the complement of whole-genome DNA methylation with the ChIP-seq data. Hypermethylated and hypomethylated regions were identified across tissues, and these locations were compared with chromatin states to better distinguish and validate regulatory elements in these tissues. Interestingly, chromatin states with the poised enhancer mark H3K4me1 in the spleen and cerebellum and CTCF in the liver displayed the greatest number of hypermethylated sites. Not surprisingly, active enhancers in the liver and spleen, and promoters in the cerebellum, displayed the greatest number of hypomethylated sites. Overall, chromatin states defined by histone marks and CTCF occupied approximately 22% of the genome in all three tissues. Furthermore, the liver and spleen displayed in common the greatest percent of active promoter (65%) and active enhancer (81%) states, and the liver and cerebellum displayed in common the greatest percent of poised enhancer (53%), repressed enhancer (68%), hypermethylated sites (75%), and hypomethylated sites (73%). In addition, both known and de novo CTCF-binding motifs were identified in all three tissues, with the highest number of unique motifs identified in the cerebellum. In summary, this study has identified the regulatory regions of genes in three tissues that play key roles in defining health and economically important traits and has set the precedent for the characterization of regulatory elements in ovine tissues using the Rambouillet reference genome.

Insecticide exposure affects intergenerational patterns of DNA methylation in the Colorado potato beetle, Leptinotarsa decemlineata.

Insecticide use is pervasive as a selective force in modern agroecosystems. Insect herbivores exposed to these insecticides have been able to rapidly evolve resistance to them, but how they are able to do so is poorly understood. One possible but largely unexplored explanation is that exposure to sublethal doses of insecticides may alter epigenetic patterns that are heritable. For instance, epigenetic mechanisms, such as DNA methylation that modifies gene expression without changing the underlying genetic code, may facilitate the emergence of resistant phenotypes in complex ways. We assessed the effects of sublethal insecticide exposure, with the neonicotinoid imidacloprid, on DNA methylation in the Colorado potato beetle, Leptinotarsa decemlineata, examining both global changes in DNA methylation and specific changes found within genes and transposable elements. We found that exposure to insecticide led to decreases in global DNA methylation for parent and F2 generations and that many of the sites of changes in methylation are found within genes associated with insecticide resistance, such as cytochrome P450s, or within transposable elements. Exposure to sublethal doses of insecticide caused heritable changes in DNA methylation in an agricultural insect herbivore. Therefore, epigenetics may play a role in insecticide resistance, highlighting a fundamental mechanism of evolution while informing how we might better coexist with insect species in agroecosystems.

Global Fund contributions to health security in ten countries, 2014-20: mapping synergies between vertical disease programmes and capacities for preventing, detecting, and responding to public health emergencies.

BACKGROUND: The Global Fund to Fight AIDS, Tuberculosis and Malaria is a robust vertical global health programme. The extent to which vertical programmes financially support health security has not been investigated. We, therefore, endeavoured to quantify the extent to which the budgets of this vertical programme support health security. We believe this is a crucial area of work as the global community works to combine resources for COVID-19 response and future pandemic preparedness. METHODS: We examined budgets for work in Kenya, Uganda, Vietnam, Democratic Republic of the Congo, Guatemala, Guinea, India, Indonesia, Nigeria, and Sierra Leone from January, 2014 to December, 2020. These ten countries were selected because of the robustness of investments and the availability of data. Using the International Health Regulations Joint External Evaluation (JEE) tool as a framework, we mapped budget line items to health security capacities. Two researchers independently reviewed each budget and mapped items to the JEE. Budgets were then jointly reviewed until a consensus was reached regarding if an item supported health security directly, indirectly, or not at all. The budgets for the study countries were inputted into a single Microsoft Excel spreadsheet and line items that mapped to JEE indicators were scaled up to their respective JEE capacity. Descriptive analyses were then done to determine the total amount of money budgeted for activities that support health security, how much was budgeted for each JEE capacity, and how much of the support was direct or indirect. FINDINGS: The research team reviewed 37 budgets. Budgets totalled US$6 927 284 966, and $2 562 063 054 (37·0%) of this mapped to JEE capacities. $1 330 942 712 (19·2%) mapped directly to JEE capacities and $1 231 120 342 (17·8%) mapped indirectly to JEE capacities. Laboratory systems, antimicrobial resistance, and the deployment of medical countermeasures and personnel received the most overall budgetary support; laboratory systems, antimicrobial resistance, and workforce development received the greatest amount of direct budgetary support. INTERPRETATION: Over one-third of the Global Fund's work also supports health security and the organisation has budgeted more than $2 500 000 000 for activities that support health security in ten countries since 2014. Although these funds were not budgeted specifically for health security purposes, recognising how vertical programmes can synergistically support other global health efforts has important implications for policy related to health systems strengthening. FUNDING: Resolve to Save Lives: An Initiative of Vital Strategies.

Whole Genome Methylation Analysis Reveals Role of DNA Methylation in Cow's Ileal and Ileal Lymph Node Responses to Mycobacterium avium subsp. paratuberculosis Infection.

Johne's Disease (JD), caused by Mycobacterium avium subsp paratuberculosis (MAP), is an incurable disease of ruminants and other animal species and is characterized by an imbalance of gut immunity. The role of MAP infection on the epigenetic modeling of gut immunity during the progression of JD is still unknown. This study investigated the DNA methylation patterns in ileal (IL) and ileal lymph node (ILLN) tissues from cows diagnosed with persistent subclinical MAP infection over a one to 4 years period. DNA samples from IL and ILLN tissues from cows negative (MAPneg) (n = 3) or positive for MAP infection (MAPinf) (n = 4) were subjected to whole genome bisulfite sequencing. A total of 11,263 and 62,459 differentially methylated cytosines (DMCs), and 1259 and 8086 differentially methylated regions (DMRs) (FDR<0.1) were found between MAPinf and MAPneg IL and ILLN tissues, respectively. The DMRs were found on 394 genes (denoted DMR genes) in the IL and on 1305 genes in the ILLN. DMR genes with hypermethylated promoters/5'UTR [3 (IL) and 88 (ILLN)] or hypomethylated promoters/5'UTR [10 (IL) and 25 (ILLN)] and having multiple functions including response to stimulus/immune response (BLK, BTC, CCL21, AVPR1A, CHRNG, GABRA4, TDGF1), cellular processes (H2AC20, TEX101, GLA, NCKAP5L, RBM27, SLC18A1, H2AC20BARHL2, NLGN3, SUV39H1, GABRA4, PPA1, UBE2D2) and metabolic processes (GSTO2, H2AC20, SUV39H1, PPA1, UBE2D2) are potential DNA methylation candidate genes of MAP infection. The ILLN DMR genes were enriched for more biological process (BP) gene ontology (GO) terms (n = 374), most of which were related to cellular processes (27.6%), biological regulation (16.6%), metabolic processes (15.4%) and response to stimulus/immune response (8.2%) compared to 75 BP GO terms (related to cellular processes, metabolic processes and transport, and system development) enriched for IL DMR genes. ILLN DMR genes were enriched for more pathways (n = 47) including 13 disease pathways compared with 36 enriched pathways, including 7 disease/immune pathways for IL DMR genes. In conclusion, the results show tissue specific responses to MAP infection with more epigenetic changes (DMCs and DMRs) in the ILLN than in the IL tissue, suggesting that the ILLN and immune processes were more responsive to regulation by methylation of DNA relative to IL tissue. Our data is the first to demonstrate a potential role for DNA methylation in the pathogenesis of MAP infection in dairy cattle.

De novo assembly of the cattle reference genome with single-molecule sequencing.

BACKGROUND: Major advances in selection progress for cattle have been made following the introduction of genomic tools over the past 10-12 years. These tools depend upon the Bos taurus reference genome (UMD3.1.1), which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. RESULTS: We present the new reference genome for cattle, ARS-UCD1.2, based on the same animal as the original to facilitate transfer and interpretation of results obtained from the earlier version, but applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly includes 2.7 Gb and is >250× more continuous than the original assembly, with contig N50 >25 Mb and L50 of 32. We also greatly expanded supporting RNA-based data for annotation that identifies 30,396 total genes (21,039 protein coding). The new reference assembly is accessible in annotated form for public use. CONCLUSIONS: We demonstrate that improved continuity of assembled sequence warrants the adoption of ARS-UCD1.2 as the new cattle reference genome and that increased assembly accuracy will benefit future research on this species.

A novel understanding of global DNA methylation in bobcat (Lynx rufus).

Epigenetic mechanisms may provide a novel prospective of bobcat (Lynx rufus) adaptation to habitat loss/fragmentation. Previous research has focused on bobcat behavior and genetics, but epigenetics has not been studied in bobcat. The aim of this study was to determine the quantity of global DNA methylation in the liver of 30 bobcats. DNA was extracted from liver samples obtained from the Vermont Fish and Wildlife Department. The percent of global DNA methylation was quantified and calculated using the MethylFlashTM Methylated DNA 5-mC Quantification Kit from Epigentek (Farmingdale, NY, USA). Age, sex, and carcass weight data were collected at sampling and analyzed with percent of global DNA methylation. Global DNA methylation was found to range from 0.46% to 2.76%. Age ranged from <1 to 12 years old and weight ranged from 3.18 to 13.61 kg. Further analysis of differential methylation may provide insight into novel means of bobcat conservation within different regions of Vermont. These results reinforce the need for genome-wide epigenetic studies in conservation biology.

Identification of bovine CpG SNPs as potential targets for epigenetic regulation via DNA methylation.

Methylation patterns established and maintained at CpG sites may be altered by single nucleotide polymorphisms (SNPs) within these sites and may affect the regulation of nearby genes. Our aims were to: 1) identify and generate a database of SNPs potentially subject to epigenetic control by DNA methylation via their involvement in creating, removing or displacing CpG sites (meSNPs), and; 2) investigate the association of these meSNPs with CpG islands (CGIs), and with methylation profiles of DNA extracted from tissues from cattle with divergent feed efficiencies detected using MIRA-Seq. Using the variant annotation for 56,969,697 SNPs identified in Run5 of the 1000 Bull Genomes Project and the UMD3.1.1 bovine reference genome sequence assembly, we identified and classified 12,836,763 meSNPs according to the nature of variation created at CpGs. The majority of the meSNPs were located in intergenic regions (68%) or introns (26.3%). We found an enrichment (p<0.01) of meSNPs located in CGIs relative to the genome as a whole, and also in differentially methylated sequences in tissues from animals divergent for feed efficiency. Seven meSNPs, located in differentially methylated regions, were fixed for methylation site creating (MSC) or destroying (MSD) alleles in the differentially methylated genomic sequences of animals differing in feed efficiency. These meSNPs may be mechanistically responsible for creating or deleting methylation targets responsible for the differential expression of genes underlying differences in feed efficiency. Our methyl SNP database (dbmeSNP) is useful for identifying potentially functional "epigenetic polymorphisms" underlying variation in bovine phenotypes.

Global DNA Methylation in the Limbic System of Cattle.

To elucidate the extent to which DNA methylation varies across multiple tissues in the brain and between animals, we have quantified global DNA methylation in tissues comprising the limbic system for six Red Angus x Simmental steers. Global DNA methylation was measured in nine regions of the bovine brain: amygdala, the bed nucleus of the stria terminalis, cingulate gyrus, dorsal raphe, hippocampus, hypothalamus, nucleus accumbens, periaqueductal gray and prefrontal cortex. DNA methylation varies among animals for each tissue type and varies among tissue types for each animal. The highest amounts of DNA methylation were found in the amygdala, cingulate gyrus and dorsal raphe, while the bed nucleus of the stria terminalis, nucleus accumbens and periaqueductal gray had the lowest amounts of DNA methylation. A heatmap sorted by k-means clustering was generated to graphically display percent DNA methylation in relation to tissue type and animal number. This is the first study to report measures of DNA methylation in the limbic system of the bovine brain and can be used to inform the cattle genomics community of expected variation in cattle brain methylation.

A Vision for Development and Utilization of High-Throughput Phenotyping and Big Data Analytics in Livestock.

Automated high-throughput phenotyping with sensors, imaging, and other on-farm technologies has resulted in a flood of data that are largely under-utilized. Drastic cost reductions in sequencing and other omics technology have also facilitated the ability for deep phenotyping of livestock at the molecular level. These advances have brought the animal sciences to a cross-roads in data science where increased training is needed to manage, record, and analyze data to generate knowledge and advances in Agriscience related disciplines. This paper describes the opportunities and challenges in using high-throughput phenotyping, "big data," analytics, and related technologies in the livestock industry based on discussions at the Livestock High-Throughput Phenotyping and Big Data Analytics meeting, held in November 2017 (see: https://www.animalgenome.org/bioinfo/community/workshops/2017/). Critical needs for investments in infrastructure for people (e.g., "big data" training), data (e.g., data transfer, management, and analytics), and technology (e.g., development of low cost sensors) were defined by this group. Though some subgroups of animal science have extensive experience in predictive modeling, cross-training in computer science, statistics, and related disciplines are needed to use big data for diverse applications in the field. Extensive opportunities exist for public and private entities to harness big data to develop valuable research knowledge and products to the benefit of society under the increased demands for food in a rapidly growing population.

An Evaluation Tool for National-Level Pandemic Influenza Planning.

The historical threat of pandemic influenza and the circulation of novel influenza viruses have led countries to strengthen their efforts in pandemic influenza preparedness planning. A cornerstone of these efforts is the creation of a comprehensive national plan that addresses all of the capacities required to prevent, detect, and respond to novel influenza outbreaks. In 2017 and 2018, the World Health Organization issued updated guidance for national pandemic planning efforts, based on lessons learned from the 2009 influenza A (H1N1) outbreak, the International Health Regulations (World Health Organization, 2005a), and other developments in health security. We have created a tool to assess national-level plans based on these updated guidelines. This tool will allow for countries to identify both strengths and weaknesses in their national plans, identify capacities and sectors that require improvement, and to help frame the updating or drafting of plans in line with the most updated guidance.

Meiotic Recombination Differences in Rams from Three Breeds of Sheep in the US.

Meiotic recombination is an important contributor to genetic variation and ensures proper chromosome segregation during gametogenesis. Previous studies suggest that at least 1 crossover (CO) per chromosome arm is important to avoid mis-segregation. While the total number of COs per spermatocyte is known to differ in mice, this is only beginning to be evaluated in sheep. This study used a cytogenetic approach to quantify and compare the number of COs per spermatocyte in rams from 3 breeds of sheep: Suffolk, Icelandic, and Targhee. In total, 2,758 spermatocytes and over 170,000 COs were examined. Suffolk rams exhibited the lowest mean number of COs (61.1 ± 0.15) compared to Icelandic (63.5 ± 0.27) and Targhee (65.9 ± 0.26) rams. Significant differences in the number of COs per spermatocyte were observed between Suffolk, Icelandic, and Targhee breeds as well as within each breed. Additionally, the number and location of COs were characterized for homologous chromosomes in a subset of spermatocytes for each ram. A positive correlation was identified between the number of COs and the length of the homologous chromosome pair. Suffolk and Icelandic rams exhibited up to 7 COs per chromosome, while Targhee rams exhibited up to 9. Further, distinct CO location preferences on homologous chromosome pairs with 1, 2, 3, and 4 COs were observed in all 3 breeds. These data in sheep will aid in elucidating the mechanism of mammalian meiotic recombination, an important contributor to genetic diversity.

Transgenerational effects of insecticides-implications for rapid pest evolution in agroecosystems.

Although pesticides are a major selective force in driving the evolution of insect pests, the evolutionary processes that give rise to insecticide resistance remain poorly understood. Insecticide resistance has been widely observed to increase with frequent and intense insecticide exposure, but can be lost following the relaxation of insecticide use. One possible but rarely explored explanation is that insecticide resistance may be associated with epigenetic modifications, which influence the patterning of gene expression without changing underlying DNA sequence. Epigenetic modifications such as DNA methylation, histone modifications, and small RNAs have been observed to be heritable in arthropods, but their role in the context of rapid evolution of insecticide resistance remain poorly understood. Here, we discuss evidence supporting how: firstly, insecticide-induced effects can be transgenerationally inherited; secondly, epigenetic modifications are heritable; and thirdly, epigenetic modifications are responsive to pesticide and xenobiotic stress. Therefore, pesticides may drive the evolution of resistance via epigenetic processes. Moreover, insect pests primed by pesticides may be more tolerant of other stress, further enhancing their success in adapting to agroecosystems. Resolving the role of epigenetic modifications in the rapid evolution of insect pests has the potential to lead to new approaches for integrated pest management as well as improve our understanding of how anthropogenic stress may drive the evolution of insect pests.