- 1Faculty of Nursing, University of Calgary, Calgary, AB, Canada
- 2Mathison Centre for Mental Health Research and Education, Hotchkiss Brain Institute, Cumming School of Medicine, University of Calgary, Calgary, AB, Canada
- 3Biological Sciences Department, University of Lethbridge, Lethbridge, AB, Canada
- 4Department of Neuroscience, University of Lethbridge, Lethbridge, AB, Canada
- 5Faculty of Veterinary Medicine & Hotchkiss Brain Institute, University of Calgary, Calgary, AB, Canada
- 6Department of Psychiatry, McGill University, Montreal, QC, Canada
- 7Departments of Psychiatry and Medical Genetics, University of Alberta, Edmonton, AB, Canada
- 8Department of Oncology, University of Alberta, Edmonton, AB, Canada
- 9Child Development Centre, University of Calgary, Calgary, AB, Canada
- 10Department of Physiology and Pharmacology, Cumming School of Medicine, University of Calgary, Calgary, AB, Canada
- 11Department of Obstetrics and Gynaecology, Cumming School of Medicine, University of Calgary, Calgary, AB, Canada
Exposure to a natural disaster in childhood can have serious, long-lasting consequences, impacting physical and mental health, development, and learning. Although many children experience negative effects after a disaster, the majority do not, and what differentiates these groups is not well understood. Some of the factors that influence disaster-related outcomes in the midst of adversity include parents’ mental health, the home environment, and socioeconomic status. Furthermore, genetics has also a role to play in how children respond to stressors. We had the opportunity to conduct a natural experiment of disaster recovery following the Alberta 2013 Flood. This paper presents the detailed protocol on prediction of resilience in Albertan families, and validation with cortisol data. In addition, data collection procedures, developing resiliency screening tools, candidate gene identification, genotyping, DNA methylation, and genomic analyses are described to achieve the research objectives. This study produced new knowledge by using pre- and post-disaster information on children’s health and development, including children’s genetics and responses to stress. This information has been identified as important to governments and other organizations invested in early child development. Our comprehensive research plan generates evidence that can be mobilized population-based approaches to improve child and family resiliency.
Introduction
Population-based cohort studies provide vital information to examine and identify etiological factors in the development of diseases important for public health (1). Cohort studies are not conducted specifically on a diseased population, and hence they provide a comprehensive understanding of environmental factors and their impact before the development of diseases. Undertaking prenatal and postnatal cohort studies is particularly important for investigating the origin of diseases believed to have their origins in the in utero environment (2–4). Negative life events can have an impact at a number of points in child development, including the impact of maternal stress in utero (5), and children experiencing stress due to natural disaster are more vulnerable to mental health conditions such as depression, anxiety, behavioral, and learning problems (6–8). These stressors may induce epigenetic changes, and also increase risk for later development of adverse mental health outcomes (9). Recently, a study has shown that high stress levels in mothers’ increase the likelihood of epigenetic changes in the child’s genome (10). The importance of cohort studies in linking early events to later life outcomes has been well established, and this may help in identifying potential risks at early stage, as well as in developing preventive measures for the well-being of individuals and communities (11, 12).
The Issue: Prenatal Stress and Early Adverse Childhood Experiences May Have Serious and Enduring Consequences
The long-term impact of adverse childhood experiences (ACEs) on physical and mental health into late adulthood is well-established (13). Children exposed to natural disasters have higher rates of depression (18–25%), anxiety (12–41%), post-traumatic stress disorder (PTSD) (5–57%) (14–16) behavioral disorders (29%), and developmental/learning problems (29%) (8) than those not exposed. The chronicity of symptoms in this group is particularly concerning, with one-third of exposed children continuing to report symptoms (8, 15) or meet Diagnostic and Statistical Manual of Mental Disorders (DSM) criteria for a mental illness 1 to 3 years post-disaster (17). Recent Canadian studies have shed critical light on the extent of the influence of prenatal exposure to natural stressors including asthma, obesity, and epigenetic changes from early childhood to 13 years (10, 18–25). The consequences of exposure of children to natural disasters are of great public health concern, exerting significant burden on health and education sectors.
Which Children Are At Risk?
The most studied child outcome, related to disaster exposure, is psychopathology. Few studies have explored risk factors associated with post-disaster outcomes, and they are limited by an inability to control for pre-disaster influences. However, disaster research consistently demonstrates the increased risk of post-disaster psychopathology in children of mothers with mental health problems (25–28). Indeed, studies have found that maternal mental health is a better predictor of post-disaster child mental health and behavior than disaster exposure (14, 29). These findings suggest that maternal mental health is a central factor in whether a child develops mental health problems post-disaster, and importantly, one of the potentially modifiable risk factor identified to date (29, 30). Results from Project Ice Storm indicate that prenatal stress from a natural disaster is linked to epigenetic effects in adolescents demonstrating an enduring impact on the child’s epigenome (10). Only a few studies have explored the role of low family resilience (Family Resilience Measure) (31) and family dysfunction (14, 32) as risk factors for children developing mental health and behavior problems post-disaster, but this evidence is equivocal. However, it is evident that family functioning worsens in disaster-affected families, and that these changes are clearly linked to poor parental mental health (14).
Genetics Play a Role in Child Resilience to Stress
Childhood trauma is a type of stressor. In general, stress has acute, delayed, and long-term effects on the body (33), and can be classified into “good stress,” “bearable stress,” and “toxic stress” (34). “Good stress” can be coped with by physiological mechanisms, encouraging healthy growth; “bearable stress” states may eventually be turned into homeostasis through successful interventions; whereas, “toxic stress,” which is characterized by the prolonged or frequent activation and dysregulation of stress response pathways, induces long-term changes and damage not only in the brain but also in the whole body (35, 36). For example, it was reported that psychological stress is associated with significantly higher oxidative stress and at least one decade of additional aging in women (37).
In response to stress, stimulation of the amygdala leads to activation of the hypothalamic–pituitary–adrenal (HPA) and sympathetic-adrenal medullary (SAM) axes (36, 38–40). This involves numerous mediators of stress, including: the neurotransmitters glutamate, noradrenaline, serotonin (41), and γ-aminobutyric acid (GABA); neuropeptides such as corticotropin-releasing factor (CRF) and the CRF family of peptides such as urocortins, orexin (42), neuropeptide Y (43), dynorphin (44), and oxytocin; and various stress-related hormones such as cortisol. These mediators bind to neuronal and non-neuronal receptors throughout the body with resultant downstream effects.
Dopamine has been associated with the stress response, and dopaminergic signaling plays important roles in fear extinction (especially in hippocampal area CA3 and the medial prefrontal cortex) (45) and is also relevant to reinforcement learning, motor control, motivation, and mood (46). Dopamine receptors belong to the G-protein coupled receptor (GPCR) family. Dopaminergic neurons can be stimulated or spontaneously activated in various patterns or modes of firing (such as bursts), with different modes being associated with different behaviors (47, 48). Indeed, alterations in the dopaminergic system have been linked to psychosis (49), anxiety disorder (50), depression (49), addiction (51), and autism (52).
The DRD2 gene on chromosome 11q23.2 encodes the dopamine D2 receptor (D2R), which is able to interact with other GPCRs (53). There is a single nucleotide polymorphism (SNP), rs1800497, a C > T substitution is located 10 kb downstream of DRD2 and leads to an amino acid substitution (p.Glu713Lys) in a gene (ANKK1) on the complementary DNA strand (54). The T/T (A1/A1) genotype has been associated with reduced D2R availability in the striatum (55), and with addictions (56), mood disorder (57), antisocial behavior (58), and gambling (59). Conversely, the C/C (A2/A2) genotype has been linked to greater vulnerability to depression after stressful life events (60).
The dopamine transporter, encoded by SLC6A3 (also known as DAT1) on chromosome 5p15.33, mediates the reuptake of dopamine from the synaptic cleft (61). Variants in DAT1 have been associated with attention deficit hyperactivity disorder (ADHD) (62), Tourette’s Syndrome (63), childhood Parkinsonism–dystonia (64), substance abuse (65), and eating disorders (66).
There is a variable number tandem repeat (VNTR) in the 3′ untranslated region of DAT1, varying from 3–11 repeats of a 40-basepair element, with the 9 and 10 repeats being the most common (67, 68). Ten-repeat carriers express higher levels of DAT1 than 9-repeat carriers (68). Under conditions of adversity such as prenatal smoking exposure (69) and psychosocial trauma (70), the DAT1 10-repeat has been associated with ADHD symptoms (62).
The catechol-O-methyltransferase enzyme encoded by the COMT gene on chromosome 22q11.2 (71). Low, intermediate, and high levels of COMT enzyme activity were first identified in red blood cells (72, 73). A SNP in COMT, rs4680 (G > A), results in a valine (Val) to methionine (Met) substitution in the enzyme (74). The Val variant is associated with 3–4 fold higher enzyme activity than the Met, and thus a lower dopamine level in the prefrontal cortex (74, 75). Colloquially, this has been termed the warrior (G: lower dopamine level, less vulnerable to stress) versus worrier (A: higher dopamine level, more vulnerable to stress) variant (76).
The serotonin transporter is encoded by SLC6A4 (also known as the 5-HTT gene), on chromosome 17q12. At 1 kb upstream is the 5-HTT linked polymorphic region (5-HTTLPR), consisting of a 20 to 23 base pair (bp) GC-rich region which is repeated 14 times in the short allele (S) and 16 times in the long allele (L). The short variant of the polymorphism is associated with reduced 5-HTT expression (77), with S/S genotype being associated with increased risk of depression and suicide attempts after stressful events (78). An A > G SNP (rs25531) in the 5-HTTLPR appears to be associated with altered expression levels, such that on the L allele (LG), the L then behaves like an S (79).
In utero exposure to natural adversities can lead to epigenome changes that can affect developmental and adaptation aspects in children. DNA methylation is a well-studied mechanism of epigenetic, which impacts utero stress (80). Furthermore, studies have shown that environmental factors can alter patterns of DNA methylation resulting into long term changes in gene function (81). Similarly, stress reactivity is typically measured by activation of the hypothalamic pituitary adrenal (HPA) system through salivary cortisol (82). Cortisol stress response can be triggered from the nature of stimulus and individual’s perception (83). Anger predicted rise in salivary cortisol and its levels correlated with anger intensity in children (84). Studies have found some similarities and differences in physiological stress responses in boys and girls (85). These differences are believed to result from differential biological mechanisms and physical and psychosocial expression of stress.
Above observations suggest that genetic analyses involving these genes, and identification of new genes associated with stress response, may help in improving our understanding of child and family resilience to natural disaster.
Resilience: Why Do Some Children Experience Serious Consequences Post-Disaster, and Others Do Not?
Resilience is the ability to adapt or “rebound from adversity when one’s ability to function has been impaired” (86). The most pressing unanswered question in child outcomes research has been why some children do and other do not experience adverse outcomes. The majority of child outcomes research examines risk without resiliency, thus identifying a broad group of children who may develop morbidities is essential. To unravel resilience and attached risks in a longitudinal cohort study, the child’s status before and after the adversity need to be identified in order to understand the child’s response to adversity. To date, contemporary worldwide disaster research exploring children’s outcomes has been unable to tackle this question due to unavailability of pre-disaster data.
Challenges of Service Provision Post-Disaster
Recent disasters show that child mental health services are unable to address the breadth of services required from prevention to specialized services post-disaster (16). In addition, the increased burden on the healthcare system is prolonged, lasting up to one-year post-disaster (87). While stepped care service provision models have been found to be effective post-disaster, they require identification of children “at risk,” the level of severity of symptoms, and decision pathways that create appropriate linkages between need and service. Mcdermott et al. (16) noted that in response to lack of service capacity, services require new processes to facilitate assessment of children and adolescents (16).
One in seven children requiring mental health services receives them post-disaster, and a recent chart review of children receiving mental healthcare services 3 years post-Katrina reported multiple missed opportunities for early identification and intervention (8). Furthermore, evidence following Hurricane Katrina demonstrated that mental health functioning deteriorated markedly at follow-up assessments for children who did not access services (88). Efficient use of post-disaster resources relies upon accurate identification of children with sub-optimal resilience and efficiently linking those children with effective and acceptable services appropriate to their needs (89).
There is a paucity of knowledge about the relationship between parental resilience and disaster-related outcomes in children (31, 32). Inherent in the identification of at-risk children and consistent with a population-based approach to resiliency is the need to identify the impact of families with low resilience post-adversity. Given that the cohort study has emerged as a popular study design for investigating the gene-environment interactions for public health (1), this design was used in our Prediction and Understanding of Resilience in Albertan Families: Longitudinal Study of Disaster Responses PURLS to describe the effects of natural disaster (Calgary flood 2013) stress on the All Our Families (AOF) in Calgary, Alberta, Canada. The aims of this study include: 1) address the gaps that exist in post-disaster outcomes due to unavailability of pre-disaster data 2) widely disseminate the child resilience screening and triage tools for future use and evaluation 3) validate the resilience screening tools using biological data and 4) promote and guide further research in early prevention/intervention aimed at improving child and family resilience and strengthening protective factors.
Methods
The All Our Families (AOF) Cohort Population
The cohort is representative of the pregnant and parenting population in a Canadian urban center (90). The AOF cohort provides longitudinal information on preconception and prenatal experiences, birth outcomes, maternal mental health during the perinatal and early childhood periods, and child health and development outcomes up to 5 years of age.
Pre-Disaster Data
Approximately 3,200 women are recruited from health care offices, communities, and through a city-wide, single provider laboratory service in Calgary, Alberta (Table 1). In brief, 78% of participants are Caucasian, 22%, non-Caucasian, mean age at delivery is 31 years, and 50% are pregnant with their first child. At the time of the flood (June 2013), ∼450 children in the AOF cohort are between the 12 and 24-month assessments, ∼950 are between their 24 and 36-month assessments, and ∼1,400 are between their 36-month and 5-year assessments. Women are recruited at <25 weeks’ gestation between August 2008 to December 2010. Women complete self-report questionnaire twice during pregnancy and at 4, 12, 24, and 36 months postpartum. The questionnaires include standardized tools and content specific items as described in Mcdonald et al. (90).
Post-Disaster Data
The PURLS resilience study is conducted on the participants that are recruited out of the AOF cohort. People who complete the 2013 flood experiences survey in the AOF cohort are invited to join PURLS study. The flood impact survey (Supplementary Materials), is administered 6 months post-flood. The flood impact questionnaire for children (Supplementary Tables 1, 2), was administered in 2015–2016. The average age of children is 3 years, and range varied form 1–4 years at the time of flood. Consent forms (Supplementary Table 3) are sent to AOF in an email through a web link, and the families who consent to participate in the survey are subsequently sent an electronic survey to determine child resilience.
Of the total 1,711 participants that complete the 2,013 flood survey and are contacted for participation, 23 are excluded due to ineligibility or disinterest in participation (Figure 1). Out of the remaining 1,688 participants, 469 (28%) people participate by completing all or some of the electronic resilience questionnaire. A total 285 participants also consent to complete the saliva collection part of the study, and 44 participate in a telephone interview about their family experiences of the flood and their child’s resilience (Figure 1).
Development of Resiliency Screening Tools
Resilience status is derived by constructing a variable by cross-classifying flood impact with children’s pre-flood (baseline) and post-flood functioning in order to identify baseline status and changes occurring following the flood. Pre- and post-flood functioning is assessed by widely used tools with demonstrated psychometric properties. Specifically, internalizing and externalizing behaviors are measured using the Brief Infant Toddler Social Emotional Assessment (BITSEA) (age 2), the National Longitudinal Survey of Children and Youth’s Behavior Scales (age 3), and the Behavioral Assessment for Children, Second Edition (BASC-2) (age 5), Devereux Student Strengths Assessment (DESSA) and the Devereux Early Childhood Assessment for Preschoolers Second Edition (DECA-P2). The Ages and Stages Questionnaire, 3rd Edition (ASQ-3) assesses children’s cognitive, emotional, social, and motor functioning at all ages. Using scoring algorithms and established clinical cut-offs for these measures, we identify patterns of pre-post functioning i.e. trajectory of functioning worsened (−1SD), stayed same or improved (+1SD), which is combined with flood impact level to identify children’s resilience status (Figure 2). Maternal reported child healthcare utilization is used to identify poor mental functioning (e.g. visits to psychologist).
Potential factors associated with child resilience status include child factors (child gender, gestational age, early developmental functioning, temperament), maternal factors (maternal adverse childhood experiences, maternal mental health, parenting competence, and stress), and family/community factors (e.g. SES, child care, social support). “Low” family resilience is defined as a decline in quality of partner relationship or parenting style (e.g., more angry, hostile) post-disaster compared to pre-disaster.
Descriptive statistics is used to describe the sample using frequencies, proportions, means (SD)/medians, IQR) as appropriate. Bivariate and multivariable analyses are conducted followed by Latent Transition Analysis (LTA) to describe how children’s resilience status changes over time. Classification and Regression Trees (CART) analysis is used to determine children’s resilience status according to different combinations of risk factors that were identified in previous analyses (e.g., LTA). These analyses are performed for each separate construct (e.g., behavior, health service utilization) and for composite low resilience, operationalized as classification as low resilient in any of the separate constructs. Separate analysis is performed for low and high risk exposure (impact) (91).
Using a regression coefficient-based scoring algorithm applied to the final predictive regression model for composite low resilience, a screening tool is developed for low child resilience post-disaster. Receiver Operating Characteristics (92) analysis is applied to determine the sensitivity, specificity, positive, and negative predictive values for different cut‐offs of the screening score, and established cut-offs for mild, moderate, and severe status.
Biological Sample Collection
All biological samples are obtained from November 2015 to November 2016. Completed kits are picked up from the home by a study employee and returned to the laboratory by mail. The collection include a total of 8 samples per child for 285 children: 3 samples per day for two days for cortisol analysis (n=6 samples), a DNA sample, and an RNA sample.
Each package contains a kit for DNA and RNA sample collection and sample collection instructions. Back in the laboratory, DNA/RNA samples are separated from saliva samples and immediately stored at −80°C until analyses. Preliminary analyses confirm that the collected DNA is viable basis for a range of genomic and epigenomic applications, ranging from single nucleotide polymorphism genotyping, to microarray analysis, to methylation analysis using a bisulphite protocol.
Using the second package, children are provided six whole saliva samples using age-appropriate Salimetrics collection tubes and child sorbettes (cotton swabs: salimetrics.com). Cortisol samples are collected at home just after awakening, 30 min later, and just before bedtime on two consecutive days. Instructions are to delay toothbrushing until after collection, and to securely hold one end of the swab and place the other end under the child’s tongue as per manufacturer’s instructions (Salimetrics LLC, Carlsbad CA).
Validation of Screening Tool Using Cortisol Data
The field of resilience suffers from limited translation of biological markers of adversity into clinical assessment or intervention. A vital component of this project is the use of such biological markers to validate the resilience screening tool. Distinguishing the biological signatures of children with low and high resilience is thus a complementary approach to validating the definition of the resilience phenotype identified through psychosocial and clinical data. The adrenal steroid, cortisol hormone plays a crucial role in the response to adverse life events, a biomarker of current stress (93–95). Based on Cao et al. (23), 285 children are an adequate sample size for the proposed biological analyses. Resiliency screening tools are used to identify resilient and non-resilient children. Biological samples collected from the AOF children are analyzed for cortisol biomarkers, which is tested for association with resiliency status to biologically validate the resiliency tool. Moreover, genetic and genomic analyses is conducted to identify genes and/or single nucleotide polymorphism (SNPs) in resilience versus non-resilient children.
Circadian Cortisol Pattern and Repeatability Across Days
Established methods in liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) (96) are used to quantify the salivary adrenal stress hormone cortisol. Cortisol is purchased from Steroids Inc (Newport, RI). Deuterium labeled internal standard: cortisol-d4, is obtained from CDN Isotopes Inc (Pointe-Claire, Quebec, Canada). HPLC grade methanol, Optima grade acetonitrile and Optima grade water are purchased from Fisher Scientific (Edmonton, AB, Canada). Ammonium fluoride (NH4F) is purchased from Sigma-Aldrich (Oakville, ON, Canada).
Sample Preparation Procedure
Saliva sample tubes are centrifuged for 20 min at 2,000 ×g, and then 1 ml of each sample is transferred into a 1.5 ml Eppendorf tube. Sample preparation is performed on a Hamilton MICROLAB STARLET workstation by diluting 75 µl of each sample or calibrator or quality control (QC) with 75 µl of protein precipitation (PPT) solution including the deuterated cortisol internal standard [9 mg/ml ZnSO4.7H2O Methanol/H2O solution (90/10 v/v)]. The mixture is vortexed for 30 s before being incubated for 20 min at 4°C. Samples are then centrifuged at 18,000 ×g for 15 min and 90 µl of supernatant is submitted for quantitation.
Lc-Esi/Mrm Analysis
An Agilent 1,200 binary liquid chromatography (LC) coupled to an AB SCIEX QTRAP® 5,500 mass spectrometer equipped with electrospray ionization (ESI) source is used. LC separation uses an Agilent ZORBAX Eclipse plus C18 column (100 × 2.1 mm, 1.8 µm particle size) at 40°C. Mobile phase A is H2O/acetonitrile (95/5, v/v, 2mM NH4F) and mobile phase B is 100% acetonitrile (2mM NH4F). The 12 min gradient is 15–70% B (0–6 min), 70–100% B (6–7 min), 100% B (7–8.5 min), 100–15% B (8.5–9 min) and is then held at 15% B for 3 min. Cortisol and cortisol-d4 are detected by positive LC-ESI/MRM, with qualifier and quantifier transitions and ionization conditions as listed in Supplementary Tables 4–8.
Quantitation Method
All steroids are quantified as the area ratio relative to the biologically identical internal standard, using a linear calibration curve with 1/× weighting (range 0.02 to 100 ng/ml). All calibration curves (11 runs) are expected to have R2 > 0.998. The lower limit of quantitation, defined as the lowest concentration that gives <20% CV and < ± 30% error will be 0.20 ng/ml. Quality control solutions at 5 and at 0.5 ng/ml are included in all runs. QC replicates are expected to have a CV 6–8%, with precision >90%.
Saliva Parameters
Any samples (expected to be <5% of approximately 1,800 samples) that fall below the LLOQ have their concentration imputed evenly within the range of 0.1 to 0.2 for subsequent analyses. Cortisol concentrations are used to calculate 12 parameters for each child: global minimum and maximum concentration across six samples, difference between awakening and 30 min after awakening on each day (plus mean, sd, and CV across the two days), and area under the curve for the day (plus mean, sd, and CV across the two days).
Linear regression model is employed for cortisol levels at awakening, morning, night, and awakening response (morning minus awakening) are used as an outcome along with resilience, child sex, flood impact as independent variables. This analysis show if resilience predicts cortisol levels in boys/girls, impacted/not impacted by flood.
Genetic, Genomic, and Transcriptomic Analyses for Candidate Genes
A. Genetic Analysis
Candidate Genes Associated With Stress
Candidate genetic markers are chosen by virtue of their involvement in the stress response by the hypothalamus–pituitary axis (HPA) and relevant mental health phenotypes (see Introduction section). Genotyping of the few candidate genes associated with stress or resilient is conducted as described below. Genotyping of DRD2 rs1800497 is conducted using a TaqMan® SNP Genotyping Assay on an Applied Biosystems ViiA™ 7 Real-Time PCR System (ThermoFisher Scientific, Canada, formerly Applied Biosystems by Life Technologies, Canada). Similarly, genotyping of DAT1 is conducted as previously described, with minor modifications (97), using primer sequences as follows: forward 5′-TGT GGT GTA GGG AAC GGC CTG AG, and reverse: 5′-CTT CCT GGT CAC GGC TCA AGG.
The COMT rs4680 SNP is genotyped using a TaqMan® SNP Genotyping Assay on an Applied Biosystems ViiA™ 7 Real-Time PCR System (ThermoFisher Scientific, Canada, formerly Applied Biosystems by Life Technologies, Canada) as previously described (98). Furthermore, 5-HTTLPR variant is genotyped as previously described, with minor modifications (99), with primer sequences of 5′-ATG CCA GCA CCT AAC CCC TAA TGT (forward) and 5′-GGA CCG CAA GGT GGG CGG GA (reverse).
B. Epigenetics
DNA methylation is measured using whole genome bisulfite sequencing, with validation of a subset of identified affected CG sites. DNA undergoes bisulfite conversion with the EpiMark Bisulfite Conversion Kit (NEB, #E3318S) followed by library generation using the NEBNext Ultra II DNA Library Kit and NEB methylated adaptor/primers as per manufacture instructions. Each library is sequenced using an Illumina NextSeq sequencer using 75 bp single-end reads. Analysis is conducted at the University of Lethbridge using open-source software packages that enables the quantitative high-resolution analysis of DNA methylation patterns from bisulfite sequencing data.
The analysis is done in two different version; relaxed and stringent. In case of the relaxed version, a CpG region is covered in a minimum 4 samples in each of the experimental groups. In case of the stringent version, the minimum number of samples requirement is increased to 10 in each of the experimental groups. The minimum coverage to estimate percent methylation for CpG position is set at 10 reads. The minimum length of each clustered CpG stretch is set at 1 CpG (i.e. the smallest possible segment of continuous methylation change is a single CpG).
In case of the relaxed analysis, over 20,000 genomic intervals in total could be examined. Intervals with adjusted p values <0.05 and methylation difference over 10% is considered significantly changed — differentially methylated regions (DMRs). The results are annotated with closest feature, distance to transcription start site (TSS), feature overlaps, gene symbols, entrez ids, and gene description and overlap with CpG islands and CpG shores. It appears that DMRs tend to be located around TSS, but only a minority of them overlap with promoter areas, and a very few overlaps with CpGs and CpG shores.
Furthermore, we determine if DNA methylation is globally altered in vulnerable versus resilient populations, and also identify loci that is specifically affected. Validation of differential methylation patterns is conducted at University of Alberta using bisulfite sequencing. Briefly, after conducting a bisulfite conversion step, PCR based amplification is performed. PCR product is cloned and sequenced, to reveal differences in methylation at each CG dinucleotide. For these validation studies, 10 regions are analyzed, which are chosen based on the extent of differential methylation and on the probability that the expression of a gene related to stress responses or immunity may be affected by the methylation. We further investigate the degree of correlation between the level of methylation in the child DNA at 5-year follow-up and the cortisol levels in the child.
C. Genomic Analysis
The Infinium Global Screening Array (comprising over 640,000 markers) is used with PsychArray add-on (50,000 markers specific for neuropsychiatric disorders) at Deltagenomics (Edmonton AB). Arrays are run in two batches. Validation of array data for two markers (DRD2 rs1800497 and COMT rs4680) is undertaken by comparison with the above TaqMan genotyping data.
Genome-wide association analysis (GWAS) is undertaken using resilience as the primary phenotype of interest, taking into account variables pertaining to the children (gender, gestational age, early development functioning, and temperament), mothers (adverse childhood experiences, maternal mental health), and families (family resilience, social economic status, child care, and social support).
D. Transcriptomic Analysis
Gene expression analysis is conducted using saliva samples and following next generation sequencing (NGS) workflow. RNA is extracted using the Oragene-RNA-RE-100 kit, with PrepIT,LTP (#PT-L2P-5), RE-100-L2N (#RE-L2N-5) according to manufacturer instructions (DNA Genotek, Ontario, Canada) and mRNA libraries are prepared using TruSeq Stranded mRNA Library Prep Kit (Cat. No. RS-122-2101) and following instructions included in the kit. RNA sequencing reactions are performed using flow cell allowing for high output 75 nt single-end reads on Illumina NextSeq500 following manufacturer’s instructions (Illumina, San Diego CA). Initial quality check is performed on the obtained raw sequencing data using FastQC software (v0.11.5) (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Adapter is trimmed and low quality reads are removed using Trim Galore with cutadapt (v1.8.dev8) (https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/). Subsequently, sequenced reads are mapped against GRCh37 (Ensemble build) reference genome using bowtie (v.2.1.0.0) genome assembly tool available in Tophat package (v.2.10) (100). Read count is performed with featureCounts (v.1.5.0-p1) (101), and statistical comparisons are conducted using DESeq2 (v1.18.1) in R computing environment (https://www.r-project.org/). Gene expression results that are obtained from NGS data is further validated with quantitative real time PCR (qRT-PCR) (102).
Results
The PURLS study was funded in 2015, and AOF recruitment and data collection was started in 2016. Physiologically, lack of resilience is most evident in changes to biomarkers related to the stress response, including immune function (cytokines; steroid hormones-cortisol) and metabolomics (trace elements and minerals). Cortisol fluctuations have been associated with exposure to disaster and trauma in children (10, 103). Ensuring the accuracy of the tool is a necessary step prior to implementing it in a disaster situation. A rigorous method of verifying the accuracy of the screening tool is to compare children identified with “low” and “high” resiliency through the screening tool against children identified through the cortisol biomarkers. Results of this study determine 1) the effectiveness of resiliency screening tool for early identification of children mental health, 2) capability of the triage decision algorithm that links the type of need and its severity identified through the resilience tool, 3) association of genetic modification to child resilience, and 4) prospects of employing results from this study in policy making for mental health.
Discussion
No other contemporary studies worldwide have been able to study child resilience related to disaster rigorously because of the lack of pre-disaster data. Based on existing gaps in the science of resiliency and its integration into practice and policy, the four key areas that are likely to have the greatest impact in terms of generating and translating evidence to practice and policy are: 1) Early identification of at-risk children. Currently, there is no method to identify children’s level of resilience post-disaster in order to triage children appropriately to effective services. For children to access needed services in a timely manner post-disaster, there is a need to develop a screening tool and decision algorithm that would accurately identify children’s risk level post-disaster and efficiently triage them into a stepped-care system. Part of this screen would necessarily entail identifying at-risk families that do not function or parent effectively post-disaster because of the adverse impact this is likely to have on child resiliency, parental capacity to support the child post-disaster, and child outcomes. This study provides needed clinical evidence to implement population-based approaches to improve child resiliency. 2) Understanding the experience of families impacted by disaster. Only one study to date has described parents’ views of and responses to their children’s needs post-disaster (104). There is a need to support the clinical component of early identification and triage of at-risk children with children’s and families’ experiences, responses and needs post-disaster. 3) At present, there are no Canadian guidelines or position statements related to early identification/intervention of resilience in children in order to support population-based prevention activities or post-disaster intervention. Consensus statements are critical for the organization of prevention and intervention strategies, particularly where evidence gaps exist.
Conclusion
To our knowledge, PURLS is the first study that provides information on pre- and post-disaster data on child resilience. This study contributes to our understanding of methods for identification of resilient and non-resilient children supported by biological information. Identification of genes involved in stress response, associating them to child resiliency and exploring genetic modifications under stress, may enhance our understanding of family and child development. Resiliency screening tool may identify children at risk and their type of need, and this information can be helpful in policy making (Figure 3).
Ethics Statement
This study was carried out in accordance with the recommendations of “Conjoint Health Research Ethics Board, University of Calgary” with signed consent form/consent implied by specific action (Electronic consent) for surveys and interviews from all subjects. The protocol was approved by the University of Calgary Conjoint Health Research Ethics Board (#REB16-0133).
Author Contributions
DK conceptualized the project. DK, IK, GM, KW-E, SK, KA, LP, SM, PA, ST, and DS designed the project and secured grant funding. MA drafted the manuscript and handled the authors/reviewer’s comments. MM, AW, and SJ provided text for the manuscript. All the authors read and approved the final manuscript.
Funding
PURLS research study was funded through Alberta Innovates, Grant #201400571. PA receives funding from the Alberta Innovates Translational Health Chair in Child and Adolescent Mental Health.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors acknowledge the contribution and support of “Prediction and Understanding of Resilience in Albertan Families: Longitudinal Study of Disaster Responses (PURLS)” participants and PURLS team members. The authors also acknowledge the contribution and support of All Our Families participants and team—All Our Families is funded through Alberta Innovates Interdisciplinary Team Grant #200700595, the Alberta Children’s Hospital Foundation, and the Max Bell Foundation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2019.00729/full#supplementary-material
References
1. Manolio TA, Bailey-Wilson JE, Collins FS. Genes, environment and the value of prospective cohort studies. Nat Rev Genet (2006) 7:812–20. doi: 10.1038/nrg1919
2. Barker DJ. Maternal nutrition, fetal nutrition, and disease in later life. Nutrition (1997) 13:807–13. doi: 10.1016/S0899-9007(97)00193-7
3. Gillman MW. Developmental origins of health and disease. N Engl J Med (2005) 353:1848–50. doi: 10.1056/NEJMe058187
4. Gluckman PD, Hanson MA, Cooper C, Thornburg KL. Effect of in utero and early-life conditions on adult health and disease. N Engl J Med (2008) 359:61–73. doi: 10.1056/NEJMra0708473
5. Kingston D, Tough S. Prenatal and postnatal maternal mental health and school-age child development: a systematic review. Matern Child Health J (2014) 18:1728–41. doi: 10.1007/s10995-013-1418-3
6. Neria Y, Nandi A, Galea S. Post-traumatic stress disorder following disasters: a systematic review. Psychol Med (2008) 38:467–80. doi: 10.1017/S0033291707001353
7. Kar N. Psychological impact of disasters on children: review of assessment and interventions. World J Pediatr (2009) 5:5–11. doi: 10.1007/s12519-009-0001-x
8. Olteanu A, Arnberger R, Grant R, Davis C, Abramson D, Asola J. Persistence of mental health needs among children affected by Hurricane Katrina in New Orleans. Prehosp Dis Med (2011) 26:3–6. doi: 10.1017/S1049023X10000099
9. Babenko O, Kovalchuk I, Metz GA. Stress-induced perinatal and transgenerational epigenetic programming of brain development and mental health. Neurosci Biobehav Rev (2015) 48:70–91. doi: 10.1016/j.neubiorev.2014.11.013
10. Cao-Lei L, Massart R, Suderman MJ, Machnes Z, Elgbeili G, Laplante DP, et al. DNA methylation signatures triggered by prenatal maternal stress exposure to a natural disaster: project Ice Storm. PLoS One (2014) 9:e107653. doi: 10.1371/journal.pone.0107653
11. Jaddoe VW, Witteman JC. Hypotheses on the fetal origins of adult diseases: contributions of epidemiological studies. Eur J Epidemiol (2006) 21:91–102. doi: 10.1007/s10654-005-5924-5
12. Golding J, Jones R, Brune MN, Pronczuk J. Why carry out a longitudinal birth survey? Paediatr Perinat Epidemiol (2009) 23 Suppl 1:1–14. doi: 10.1111/j.1365-3016.2008.01009.x
13. Logan-Greene P, Green S, Nurius PS, Longhi D. Distinct contributions of adverse childhood experiences and resilience resources: a cohort analysis of adult physical and mental health. Soc Work Health Care (2014) 53:776–97. doi: 10.1080/00981389.2014.944251
14. Mcfarlane AC. Posttraumatic phenomena in a longitudinal study of children following a natural disaster. J Am Acad Child Adolesc Psychiatry (1987) 26:764–9. doi: 10.1097/00004583-198709000-00025
15. Yelland C, Robinson P, Lock C, La Greca AM, Kokegei B, Ridgway V, et al. Bushfire impact on youth. J Trauma Stress (2010) 23:274–7. doi: 10.1002/jts.20521
16. Mcdermott BM, Cobham VE. A stepped-care model of post-disaster child and adolescent mental health service provision. Eur J Psychotraumatol (2014) 5. doi: 10.3402/ejpt.v5.24294
17. Ularntinon S, Piyasil V, Ketumarn P, Sitdhiraksa N, Pityaratstian N, Lerthattasilp T, et al. Assessment of psychopathological consequences in children at 3 years after tsunami disaster. J Med Assoc Thai (2008) 91 Suppl 3:S69–75.
18. King S, Laplante DP. The effects of prenatal maternal stress on children's cognitive development: Project Ice Storm. Stress Int J Biol Stress (2005) 8:35–45. doi: 10.1080/10253890500108391
19. Laplante DP, Brunet A, Schmitz N, Ciampi A, King S. Project ice storm: Prenatal maternal stress affects cognitive and linguistic functioning in 51/2-year-old children. J Am Acad Child and Adolesc Psychiatry (2008) 47:1063–72. doi: 10.1097/CHI.0b013e31817eec80
20. Dancause KN, Laplante DP, Fraser S, Brunet A, Ciampi A, Schmitz N, et al. Prenatal exposure to a natural disaster increases risk for obesity in 5(1/2)-year-old children. Pediatr Res (2012) 71:126–31. doi: 10.1038/pr.2011.18
21. King S, Dancause K, Turcotte-Tremblay AM, Veru F, Laplante DP. Using natural disasters to study the effects of prenatal maternal stress on child health and development. Birth Defects Res Part C-Embryo Today-Rev (2012) 96:273–88. doi: 10.1002/bdrc.21026
22. Dancause KN, Veru F, Andersen RE, Laplante DP, King S. Prenatal stress due to a natural disaster predicts insulin secretion in adolescence. Early Hum Dev (2013) 89:773–6. doi: 10.1016/j.earlhumdev.2013.06.006
23. Cao X, Laplante DP, Brunet A, Ciampi A, King S. Prenatal maternal stress affects motor function in 5(1/2)-year-old children: project ice storm. Dev Psychobiol (2014) 56:117–25. doi: 10.1002/dev.21085
24. Turcotte-Tremblay AM, Lim R, Laplante DP, Kobzik L, Brunet A, King S. Prenatal maternal stress predicts childhood asthma in girls: Project Ice Storm. Biomed Res Int (2014) doi: 10.1155/2014/201717
25. Walder DJ, Laplante DR, Sousa-Pires A, Veru F, Brunet A, King S. Prenatal maternal stress predicts autism traits in 61/2 year-old children: Project Ice Storm. Psychiatry Res (2014) 219:353–60. doi: 10.1016/j.psychres.2014.04.034
26. Swenson CC, Saylor CF, Powell MP, Stokes SJ, Foster KY, Belter RW. Impact of a natural disaster on preschool children: Adjustment 14 months after a hurricane. Am J Orthopsychiatry (1996) 66:122–30. doi: 10.1037/h0080162
27. Jones RT, Ribbe DP, Cunningham PB, Weddle JD, Langley AK. Psychological impact of fire disaster on children and their parents. Behav Modif (2002) 26:163–86. doi: 10.1177/0145445502026002003
28. Mirzamani M, Bolton D. PTSD symptoms of mothers following occurrence of a disaster affecting their children. Psychol Rep (2002) 90:431–8. doi: 10.2466/pr0.2002.90.2.431
29. Spell AW, Kelley ML, Wang J, Self-Brown S, Davidson KL, Pellegrin A, et al. The moderating effects of maternal psychopathology on children's adjustment post-Hurricane Katrina. J Clin Child Adolesc Psychol (2008) 37:553–63. doi: 10.1080/15374410802148210
30. Mcfarlane AC, Policansky SK, Irwin C. A longitudinal-study of the psychological morbidity in children due to a natural disaster. Psychol Med (1987) 17:727–38. doi: 10.1017/S0033291700025964
31. Mcdermott BM, Cobham VE, Berry H, Stallman HM. Vulnerability factors for disaster-induced child post-traumatic stress disorder: the case for low family resilience and previous mental illness. Aust N Z J Psychiatry (2010) 44:384–9. doi: 10.3109/00048670903489916
32. Mcdermott BM, Cobham VE. Family functioning in the aftermath of a natural disaster. BMC Psychiatry (2012) 12:55. doi: 10.1186/1471-244X-12-55
33. Joëls M, Baram TZ. The neuro-symphony of stress. Nat Rev Neurosci (2009) 10:459. doi: 10.1038/nrn2632
34. Shonkoff JP, Boyce WT, Mcewen BS. Neuroscience, molecular biology, and the childhood roots of health disparities: building a new framework for health promotion and disease prevention. JAMA (2009) 301:2252–9. doi: 10.1001/jama.2009.754
35. Franke HA. Toxic stress: effects, prevention and treatment. Children (2014) 1:390–402. doi: 10.3390/children1030390
36. Bucci M, Marques SS, Oh D, Harris NB. Toxic stress in children and adolescents. Adv Pediatr (2016) 63:403–28. doi: 10.1016/j.yapd.2016.04.002
37. Epel ES, Blackburn EH, Lin J, Dhabhar FS, Adler NE, Morrow JD, et al. Accelerated telomere shortening in response to life stress. Proc Natl Acad Sci USA (2004) 101:17312–5. doi: 10.1073/pnas.0407162101
38. Mcewen BS, Bowles NP, Gray JD, Hill MN, Hunter RG, Karatsoreos IN, et al. Mechanisms of stress in the brain. Nat Neurosci (2015) 18:1353. doi: 10.1038/nn.4086
39. Janak PH, Tye KM. From circuits to behaviour in the amygdala. Nature (2015) 517:284. doi: 10.1038/nature14188
40. Scheuer S, Wiggert N, Brückl TM, Awaloff Y, Uhr M, Lucae S, et al. Childhood abuse and depression in adulthood: The mediating role of allostatic load. Psychoneuroendocrino (2018) 94:134–42. doi: 10.1016/j.psyneuen.2018.04.020
41. Popoli M, Yan Z, Mcewen BS, Sanacora G. The stressed synapse: the impact of stress and glucocorticoids on glutamate transmission. Nat Rev Neurosci (2012) 13:22. doi: 10.1038/nrn3138
42. Berridge CW, España RA, Vittoz NM. Hypocretin/orexin in arousal and stress. Brain Res (2010) 1314:91–102. doi: 10.1016/j.brainres.2009.09.019
43. Witt SH, Buchmann AF, Blomeyer D, Nieratschker V, Treutlein J, Esser G, et al. An interaction between a neuropeptide Y gene polymorphism and early adversity modulates endocrine stress responses. Psychoneuroendocrino (2011) 36:1010–20. doi: 10.1016/j.psyneuen.2010.12.015
44. Bailey CR, Cordell E, Sobin SM, Neumeister A. Recent progress in understanding the pathophysiology of post-traumatic stress disorder. CNS Drugs (2013) 27:221–32. doi: 10.1007/s40263-013-0051-4
45. Wen J-L, Xue L, Wang R-H, Chen Z-X, Shi Y-W, Zhao H. Involvement of the dopaminergic system in the consolidation of fear conditioning in hippocampal CA3 subregion. Behav Brain Res (2015) 278:527–34. doi: 10.1016/j.bbr.2014.10.049
46. Watabe-Uchida M, Zhu L, Ogawa SK, Vamanrao A, Uchida N. Whole-brain mapping of direct inputs to midbrain dopamine neurons. Neuron (2012) 74:858–73. doi: 10.1016/j.neuron.2012.03.017
47. Kopin I. Catecholamine metabolism: basic aspects and clinical significance. Pharmacol Rev (1985) 37:333–64.
48. Overton P, Clark D. Burst firing in midbrain dopaminergic neurons. Brain Res Rev (1997) 25:312–34. doi: 10.1016/S0165-0173(97)00039-8
49. Grace AA. Dysregulation of the dopamine system in the pathophysiology of schizophrenia and depression. Nat Rev Neurosci (2016) 17:524. doi: 10.1038/nrn.2016.57
50. Meyer A, Hajcak G, Torpey-Newman DC, Kujawa A, Klein DN. Enhanced error-related brain activity in children predicts the onset of anxiety disorders between the ages of 6 and 9. J Abnorm Psychol (2015) 124:266. doi: 10.1037/abn0000044
51. Pascoli V, Terrier J, Hiver A, Lüscher C. Sufficiency of mesolimbic dopamine neuron stimulation for the progression to addiction. Neuron (2015) 88:1054–66. doi: 10.1016/j.neuron.2015.10.017
52. Damasio AR, Maurer RG. A neurological model for childhood autism. Arch Neurol (1978) 35:777–86. doi: 10.1001/archneur.1978.00500360001001
53. Rocheville M, Lange DC, Kumar U, Patel SC, Patel RC, Patel YC. Receptors for dopamine and somatostatin: formation of hetero-oligomers with enhanced functional activity. Science (2000) 288:154–7. doi: 10.1126/science.288.5463.154
54. Neville MJ, Johnstone EC, Walton RT. Identification and characterization of ANKK1: a novel kinase gene closely linked to DRD2 on chromosome band 11q23.1. Hum Mutat (2004) 23:540–5. doi: 10.1002/humu.20039
55. Pohjalainen T, Rinne JO, Nagren K, Lehikoinen P, Anttila K, Syvalahti EK, et al. The A1 allele of the human D2 dopamine receptor gene predicts low D2 receptor availability in healthy volunteers. Mol Psychiatry (1998) 3:256–60. doi: 10.1038/sj.mp.4000350
56. Blum K, Thanos PK, Gold MS. Dopamine and glucose, obesity, and reward deficiency syndrome. Front Psychol (2014) 5:919. doi: 10.3389/fpsyg.2014.00919
57. Zhang L, Hu L, Li X, Zhang J, Chen B. The DRD2 rs1800497 polymorphism increase the risk of mood disorder: evidence from an update meta-analysis. J Affect Disord (2014) 158:71–7. doi: 10.1016/j.jad.2014.01.015
58. Lu RB, Lee JF, Huang SY, Lee SY, Chang YH, Kuo PH, et al. Interaction between ALDH2* 1* 1 and DRD2/ANKK1 TaqI A1A1 genes may be associated with antisocial personality disorder not co-morbid with alcoholism. (2012) 17:865–74. doi: 10.1111/j.1369-1600.2010.00268.x
59. Gray JC, Mackillop J. Genetic basis of delay discounting in frequent gamblers: examination of a priori candidates and exploration of a panel of dopamine-related loci. Brain Behav (2014) 4:812–21. doi: 10.1002/brb3.284
60. Elovainio M, Jokela M, Kivimaki M, Pulkki-Raback L, Lehtimaki T, Airla N, et al. Genetic variants in the DRD2 gene moderate the relationship between stressful life events and depressive symptoms in adults: Cardiovascular risk in young Finns study. Psychosom Med (2007) 69:391–5. doi: 10.1097/psy.0b013e31806bf365
61. Vandenbergh DJ, Persico AM, Uhl GR. A human dopamine transporter cDNA predicts reduced glycosylation, displays a novel repetitive element and provides racially-dimorphic TaqI RFLPs. Mol Brain Res (1992) 15:161–6. doi: 10.1016/0169-328X(92)90165-8
62. Cornish KM, Manly T, Savage R, Swanson J, Morisano D, Butler N, et al.. Association of the dopamine transporter (DAT1) 10/10-repeat genotype with ADHD symptoms and response inhibition in a general population sample. Mol Psychiatry (2005) 10:686. doi: 10.1038/sj.mp.4001641
63. Yoon DY, Rippel CA, Kobets AJ, Morris CM, Lee JE, Williams PN, et al. Dopaminergic polymorphisms in Tourette syndrome: association with the DAT gene (SLC6A3). Am J Med Genet Part B: Neuropsychiatr Genet (2007) 144:605–10. doi: 10.1002/ajmg.b.30466
64. Kurian MA, Li Y, Zhen J, Meyer E, Hai N, Christen H-J, et al. Clinical and molecular characterisation of hereditary dopamine transporter deficiency syndrome: an observational cohort and experimental study. Lancet Neurol (2011) 10:54–62. doi: 10.1016/S1474-4422(10)70269-6
65. Sabol SZ, Nelson ML, Fisher C, Gunzerath L, Brody CL, Hu S, et al. A genetic association for cigarette smoking behavior. Health Psychol (1999) 18:7. doi: 10.1037/0278-6133.18.1.7
66. Shinohara M, Mizushima H, Hirano M, Shioe K, Nakazawa M, Hiejima Y, et al. Eating disorders with binge-eating behaviour are associated with the s allele of the 3′-UTR VNTR polymorphism of the dopamine transporter gene. J Psychiatry Neurosci (2004) 29:134.
67. Sano A, Kondoh K, Kakimoto Y, Kondo I. A 40-nucleotide repeat polymorphism in the human dopamine transporter gene. Hum Genet (1993) 91:405–6. doi: 10.1007/BF00217369
68. Rowe DC, Stever C, Gard JM, Cleveland HH, Sanders ML, Abramowitz A, et al. The relation of the dopamine transporter gene (DAT1) to symptoms of internalizing disorders in children. Behav Genet (1998) 28:215–25. doi: 10.1023/A:1021427314941
69. Becker K, El-Faddagh M, Schmidt MH, Esser G, Laucht M. Interaction of dopamine transporter genotype with prenatal smoke exposure on ADHD symptoms. J Pediatr (2008) 152:263–269. e261. doi: 10.1016/j.jpeds.2007.07.004
70. Laucht M, Skowronek MH, Becker K, Schmidt MH, Esser G, Schulze TG, et al. Interacting effects of the dopamine transporter gene and psychosocial adversity on attention-deficit/hyperactivity disorder symptoms among 15-year-olds from a high-risk community sample. Arch Gen Psychiatry (2007) 64:585–90. doi: 10.1001/archpsyc.64.5.585
71. Grossman MH, Emanuel BS, Budarf ML. Chromosomal mapping of the human catechol-O-methyltransferase gene to 22q11.1→q11.2. Genomics (1992) 12:822–5. doi: 10.1016/0888-7543(92)90316-K
72. Axelrod J, Cohn CK. Methyltransferase enzymes in red blood cells. J Pharmacol Exp Ther (1971) 176:650–4.
73. Weinshilboum RM, Raymond FA. Inheritance of low erythrocyte catechol-o-methyltransferase activity in man. Am J Hum Genet (1977) 29:125–35.
74. Lachman HM, Papolos DF, Saito T, Yu Y-M, Szumlanski CL, Weinshilboum RM. Human catechol-O-methyltransferase pharmacogenetics: description of a functional polymorphism and its potential application to neuropsychiatric disorders. Pharmacogenetics (1996) 6:243–50. doi: 10.1097/00008571-199606000-00007
75. Chen J, Lipska BK, Halim N, Ma QD, Matsumoto M, Melhem S, et al. Functional analysis of genetic variation in catechol-O-methyltransferase (COMT): effects on mRNA, protein, and enzyme activity in postmortem human brain. Am J Hum Genet (2004) 75:807–21. doi: 10.1086/425589
76. Stein DJ, Newman TK, Savitz J, Ramesar R. Warriors versus worriers: the role of COMT gene variants. CNS Spectr (2006) 11:745–8. doi: 10.1017/S1092852900014863
77. Lesch K-P, Bengel D, Heils A, Sabol SZ, Greenberg BD, Petri S, et al. Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science (1996) 274:1527–31. doi: 10.1126/science.274.5292.1527
78. Caspi A, Sugden K, Moffitt TE, Taylor A, Craig IW, Harrington H, et al. Influence of life stress on depression: moderation by a polymorphism in the 5-HTT gene. Science (2003) 301:386–9. doi: 10.1126/science.1083968
79. Kraft JB, Slager SL, Mcgrath PJ, Hamilton SP. Sequence analysis of the serotonin transporter and associations with antidepressant response. Biol Psychiatry (2005) 58:374–81. doi: 10.1016/j.biopsych.2005.04.048
80. Cao-Lei L, De Rooij SR, King S, Matthews SG, Metz G.a.S., Roseboom TJ, et al. Prenatal stress and epigenetics. Neurosci Biobehav Rev (2017) S0149-7634(16)30726-6. doi: 10.1016/j.neubiorev.2017.05.016
81. Turecki G, Meaney MJ. Effects of the social environment and stress on glucocorticoid receptor gene methylation: a systematic review. Biol Psychiatry (2016) 79:87–96. doi: 10.1016/j.biopsych.2014.11.022
82. Hellhammer DH, Wust S, Kudielka BM. Salivary cortisol as a biomarker in stress research. Psychoneuroendocrino (2009) 34:163–71. doi: 10.1016/j.psyneuen.2008.10.026
83. Kudielka BM, Hellhammer DH, Wust S. Why do we respond so differently? Reviewing determinants of human salivary cortisol responses to challenge. Psychoneuroendocrino (2009) 34:2–18. doi: 10.1016/j.psyneuen.2008.10.004
84. Martorell GA. Maternal and child adrenocortical responses as a function of attachment classification, temperament, and maternal caregiving attributions. Dissertation Abstracts International: Section B: The Sciences and Engineering (2002) 62(7–B):3401.
85. Susman EJ, Pajer K. Biology-Behavior Integration and Antisocial Behavior in Girls. In: Putallaz M, Bierman KL, editors. Duke series in child development and public policy. Aggression, antisocial behavior, and violence among girls: A developmental perspective. Guilford Publications (2004). p. 23–47.
86. Karatsoreos IN, Mcewen BS. Psychobiological allostasis: resistance, resilience and vulnerability. Trends Cogn Sci (2011) 15:576–84. doi: 10.1016/j.tics.2011.10.005
87. Rodriguez JJ, Kohn R. Use of mental health services among disaster survivors. Curr Opin Psychiatry (2008) 21:370–8. doi: 10.1097/YCO.0b013e328304d984
88. Gibbs L, Waters E, Bryant RA, Pattison P, Lusher D, Harms L, et al. Beyond bushfires: Community, resilience and recovery — a longitudinal mixed method study of the medium to long term impacts of bushfires on mental health and social connectedness. BMC Public Health (2013) 13:1036. doi: 10.1186/1471-2458-13-1036
89. Kulig J, Botey AP. Facing a wildfire: What did we learn about individual and community resilience? Nat Hazards (2016) 82:1919–29. doi: 10.1007/s11069-016-2277-1
90. Mcdonald SW, Lyon AW, Benzies KM, Mcneil DA, Lye SJ, Dolan SM, et al. The All Our Babies pregnancy cohort: design, methods, and participant characteristics. BMC Pregnancy Childbirth (2013) 13.(Suppl 1):S2. doi: 10.1186/1471-2393-13-S1-S2
91. Hopkins KD, Zubrick SR, Taylor CL. Resilience amongst Australian aboriginal youth: An ecological analysis of factors associated with psychosocial functioning in high and low family risk contexts. PloS One (2014) 9(7):e102820. doi: 10.1371/journal.pone.0102820
92. Yonelinas AP. Receiver-operating characteristics in recognition memory — evidence for a dual-process model. J Exp Psychol Learn Memory Cogn (1994) 20:1341–54. doi: 10.1037/0278-7393.20.6.1341
93. Harkness KL, Stewart JG, Wynne-Edwards KE. Cortisol reactivity to social stress in adolescents: Role of depression severity and child maltreatment. Psychoneuroendocrino (2011) 36:173–81. doi: 10.1016/j.psyneuen.2010.07.006
94. Harkness KL, Bagby RM, Stewart JG, Larocque CL, Mazurka R, Strauss JS, et al. Childhood emotional and sexual maltreatment moderate the relation of the serotonin transporter gene to stress generation. J Abnorm Psychol (2015) 124:275–87. doi: 10.1037/abn0000034
95. Mazurka R, Wynne-Edwards KE, Harkness KL. Stressful Life Events Prior to Depression Onset and the Cortisol Response to Stress in Youth with First Onset Versus Recurrent Depression. J Abnorm Child Psychol (2016) 44:1173–84. doi: 10.1007/s10802-015-0103-y
96. Koren L, Ng ES, Soma KK, Wynne-Edwards KE. Sample preparation and liquid chromatography-tandem mass spectrometry for multiple steroids in mammalian and avian circulation. PLoS One (2012) 7:e32496. doi: 10.1371/journal.pone.0032496
97. Vanness SH, Owens MJ, Kilts CD. The variable number of tandem repeats element in DAT1 regulates in vitro dopamine transporter density. BMC Genet. (2005) 6:55. doi: 10.1186/1471-2156-6-55
98. Lodhi RJ, Wang Y, Rossolatos D, Macintyre G, Bowker A, Crocker C, et al. Investigation of the COMT Val158Met variant association with age of onset of psychosis, adjusting for cannabis use. Brain Behav (2017) 7:e00850. doi: 10.1002/brb3.850
99. Nan J, Yuan H, Li K, Jin Y, Yu M. 5-HTT SS genotype is associated with the pro-nociceptive sensation by alcoholic sting. Cell Biochem Biophys (2014) 68:629–33. doi: 10.1007/s12013-013-9759-5
100. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc (2012) 7:562–78. doi: 10.1038/nprot.2012.016
101. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics (2014) 30:923–30. doi: 10.1201/b16589
102. Consortium SM-I. A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the sequencing quality control consortium. Nat Biotechnol (2014) 32:903–14. doi: 10.1038/nbt.2957
103. Mondelli V, Dazzan P, Hepgul N, Di Forti M, Aas M, D'Albenzio A. Abnormal cortisol levels during the day and cortisol awakening response in first-episode psychosis: the role of stress and of antipsychotic treatment. Schizophr Res (2010) 116(2–3):234–42. doi: 10.1016/j.schres.2009.08.013
Keywords: stress, resiliency, genetics, protocol, child development
Citation: Kingston D, Mughal MK, Arshad M, Kovalchuk I, Metz GAS, Wynne-Edwards K, King S, Jiang S, Postovit L, Wajid A, McDonald S, Slater DM, Tough SC, Aitchison K and Arnold P (2019) Prediction and Understanding of Resilience in Albertan Families: Longitudinal Study of Disaster Responses (PURLS) – Protocol. Front. Psychiatry 10:729. doi: 10.3389/fpsyt.2019.00729
Received: 07 May 2019; Accepted: 11 September 2019;
Published: 31 October 2019.
Edited by:
Konstantin Loganovsky, National Academy of Medical Science of Ukraine (NAMS Ukraine), UkraineCopyright © 2019 Kingston, Mughal, Arshad, Kovalchuk, Metz, Wynne-Edwards, King, Jiang, Postovit, Wajid, McDonald, Slater, Tough, Aitchison and Arnold. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dawn Kingston, ZGF3bi5raW5nc3RvbkB1Y2FsZ2FyeS5jYQ==