Abstract
The underlying molecular mechanisms responsible for the etiology of autism and its sex-biased prevalence remain largely elusive. Abnormally shortened telomeres have recently been associated with autism. We have previously shown that children with non-syndromic autism exhibit a sexually dimorphic pattern of relative telomere length (RTL). Only male children with autism have significantly shorter RTLs than the healthy controls and paired siblings. Autistic females have substantially longer RTLs than autistic males. Aberrantly high levels of oxidative stress plays a fundamental role in the pathophysiology of autism, and telomeres are thought to be susceptible to oxidative damage due to their high guanine-repeat content. Employing a quantitative PCR (qPCR)-based method, telomeric oxidized base lesions were measured using genomic DNA extracted from saliva samples, and levels of telomeric RNA transcripts know as TERRA were evaluated using reverse transcriptase qPCR technique. Our data show that the autistic children exhibit substantially higher levels of oxidative base lesions at their telomeres than the healthy controls and paired siblings. Intriguingly, despite having significantly longer RTLs, female children with autism have even higher levels of telomeric oxidized bases than their male counterparts. Furthermore, despite having significantly shorter RTLs, the male children with autism exhibit lower levels of TERRA expression from the short arms of chromosomes 17 and X/P compared to their individually-matched healthy controls. These findings open a fresh angle into autism. Abnormal TL and high levels of telomeric oxidized bases may serve as biomarkers for childhood autism.
Introduction
Childhood autism, a severe subgroup of autism spectrum disorders (ASD), is a lifelong pervasive neurodevelopmental disorder with an early onset characterized by deficits in reciprocal social interaction, impaired verbal and nonverbal communication, and restricted and repetitive behavior [1]. The number of boys diagnosed with autism is 4-5 times greater than that of girls [2]. Genetic factors play a major role in the etiology of autism, but the number of genomic variations characterized is overwhelmingly high and extremely heterogeneous with a high degree of variance in phenotypic presentations among different cases [3, 4]. Autism has a complex pathophysiology with a number of outstanding features and key roles for abnormal alterations in oxidative stress, immune and neuroimmune regulation and inflammation, metabolism, and mitochondrial function [5–7]. Despite heightened worldwide attention and highly significant advances in autism research in recent years, the cellular processes involved and the underlying molecular mechanisms responsible for the etiology and pathophysiology of autism remain largely elusive.
Two recent studies reported a strong association between abnormally shortened telomeres and autism implicating a potential role for telomere biology in the pathogenesis of autism [8, 9]. Telomeres are specialized and dynamic nucleoprotein structures at the terminal ends of the mammalian chromosomes composed of tandem repeats of double-stranded (ds) 5'-TTAGGG-3' extending for several thousand base pairs (typically 4-20 kb in humans) and 150-200 nucleotides of the so-called G-rich single-stranded (ss) overhang at the very end. The G-rich overhang is looped back in a large lariat structure with D- and T-loops. Telomeric sequences are associated with specific protein factors and enzyme complexes, which are also involved in maintaining telomere length (TL) and/or structure [10]. TL is an important factor in maintaining telomere structure and functions, which include protecting the chromosome ends as well as genome integrity [11]. Moreover, the specific factors associated with the telomeric sequences have a number of extra-telomeric functions [12]. Telomere shortening elevates dispersion of the telomere associated factors to chromosome internal loci [13], and critically short telomeres are a prelude to cellular senescence and apoptosis [14]. On the other hand, abnormally long telomeres are also indicative of unusual telomere structure and alternative telomere maintenance and function [10, 15].
Despite their classic characterization as constitutive heterochromatin and high levels of protection, telomeres are both susceptible to DNA damage and accessible to transcriptional machinery [13, 16, 17]. In addition to telomere erosion in replicating somatic cells, which lack the reverse transcriptase activity of telomerase, TL is dependent on a number of genetic plus internal and external environmental factors [18]. TL is impacted by oxidative stress induced damage and a variety of DNA-repair mechanisms under unique direction/influence of the protein factors and enzymes involved in telomere protection and maintenance [19, 20]. If not properly structured and/or protected, telomeres are particularly susceptible to oxidative damage because of their repetitive sequence and high consecutive guanine base content. In fact, 8-oxo-7,8-dihydrodeoxyguanosine (8-oxoG) is a major product of oxidative stress induced DNA damage [16, 21, 22]. Surprisingly, despite having a highly specialized heterochromatin structure and being devoid of classically-defined segments as genes, telomeres are transcribed into long non-coding RNAs referred to as TERRA [23]. TERRA molecules not only seem to be involved in telomere structure and TL regulation, but also capable of impacting the expression of sub-telomeric genes and a number of chromosome internal target genes, as well as a large network of interacting proteins [24]. TERRA expression levels depend on telomere length and structure [25, 26].
In line with earlier reports indicating key roles for oxidative stress and aberrant alterations in cellular redox status (i.e., anti-oxidant defense capacity) in autism [27–30], a comprehensive study from S. Jill James group provided convincing evidence on the elevated levels of both oxidized protein derivatives and genome-wide oxidized DNA adducts (including 8-oxo-dG) in children with autism but not their healthy control or paired sibling groups [5]. These observations were further confirmed by follow-up works from the same research lab by post-mortem analysis of the brain tissues from individual with autism and also a mouse model of autism [31, 32].
In our earlier report, we have presented evidence showing a sexually dimorphic pattern of average RTL status among children with non-syndromic autism. While autistic boys followed a homogenous pattern of abnormally shortened telomeres, there was no evidence of unusual telomere shortening in girls. In fact, in a sharp contrast to the autistic male children, some of the female children with autism had substantially longer average RTL than their matched healthy controls [33]. However, many aspects of telomere biology such as TL distribution among different chromosome ends, structure, transcriptional activity, and levels of telomeric oxidation remain unexplored. Here, by employing a quantitative PCR (qPCR)-based method previously described by Nathan O’Callaghan et al. [34], we examine the levels of telomeric oxidized base lesions among male and female autistic children. We also take a look at the levels of TERRA transcripts coming from certain chromosome ends that are of interest in autism. To our knowledge, this is the first report on the status of telomeric oxidation and TERRA transcription in ASD research.
Results
Subjects and Study Groups
Demographics and characteristics of the subjects including 24 children (14M:10F) with non-syndromic typical autism from independent simplex families, 10 apparently healthy paired siblings (8M:2F), and 24 sex-, age-, and location-matched healthy control children have been described in detail in our previous paper [33]. Clinically, all 14 male and 10 female autistic children were in “moderate to severe” category. Five out of 14 male children with autism and six out of 10 autistic female children were nonverbal. Eighteen (14M:4F) of the children with autism were drug naïve (i.e., first time diagnosis) at the time of sampling and four of the females had extended wash periods [33].
The levels of oxidized residues in telomeric regions of autistic children are significantly higher than the healthy controls
Comparison of the average numbers of telomeric FPG-sensitive sites (i.e., oxidized base lesions) between different groups, after adjustment for covariates, showed that children with autism had significantly higher levels of telomeric base oxidation than both the healthy controls (β = 0.495, p = 0.0002) and apparently healthy paired siblings (β = 0.410, p = 0.003) (Fig. 1A, left). In a gender stratified manner, both male and female autistic children showed substantially higher levels of telomeric oxidized base lesions compared to their sex-, age-, and location-matched healthy controls (β = 0.506, p = 0.008; and β = 0.598, p = 0.002, respectively) (Fig. 1A, middle and right). Both male and female groups of children with autism also had significantly higher concentrations of telomeric oxidation compared with their respective paired siblings (β = 0.375, p = 0.048; and β = 0.622, p = 0.009, respectively). However, the difference for the female autistic children may not be very reliable because of the small number of paired siblings. Interestingly, the female autistic children had considerably higher levels of telomeric oxidation compared to their male counterparts, although this difference was of marginal statistical significance (β = 0.424, p = 0.053). It is also interesting to note that taking the sibling group out of the analysis further signifies the differences in the levels of telomeric oxidized bases between the children with autism either in combination or as male or female groups compared with their respective healthy control groups (β = 0.583, p = 0.00001; β = 0.611, p = 0.001; and β = 0.764, p = 0.001, respectively).
Similar paterns of elevated levels of telomeric oxidized bases among autistic male and female children groups
Organizing the subjects as individually-matched groups, showed that with the exception of a very few cases, and in a sharp contrast to the sexually dimorphic patterns of average relative telomere lengths (RTLs) [33], both male and female autistic children exhibited a similar pattern of elevated telomeric oxidized bases compared with their individually-matched healthy controls and apparently healthy paired siblings (Fig. 1B and 1C). Also in contrast to the increased telomere attrition with age [33], the levels of telomeric oxidized base lesions remained rather constant in both autistic and control subjects with increasing age (Fig. 2). Similar to our observation on a lack of correlation between ADI-R total scores and RTL [33], there were no significant associations between the levels of telomeric FPG-sensitive sites and total ADI-R scores among the patients (Supplementary Material, Fig. S2) except for seemingly interesting contrasting trends for male vs. female nonverbal autistic children (Fig. S2, panel E).
Differential TERRA expression in male children with autism
We have previously reported a detailed evaluation of GAPDH and YWHAZ as appropriate combined reference genes for normalization of RT-qPCR data for saliva transcript level analysis in autistic male children [35]. We performed a similar systematic two-phased analysis among our autistic and contorl female samples. Although GAPDH and YWHAZ again were ranked as the best candidate reference genes among the panel of nine housekeeping genes examined, they did not pass the evaluation criteria. In fact, even a combination of 4 genes in the 2nd phase of analysis did not seem good enough for RT-qPCR data normalization (Supplementray Material, Fig. S3). Thus, we limited our investigation of telomeric RNA transcription activity to the male samples alone.
We initially attempted to measure total TERRA levels in our male children groups following a protocol originally reported by Holzmann and colleagues [36]. However, despite careful treatment of the RNA samples with DNase I (see Materials and Methods), this approach did not produce a reliable outcome due to low levels of TERRA expression and background signals coming from even very small amounts of residual genomic DNA in the no-reverse transcription (NRT) control reactions (data not shown). Surprisingly, Holzmann and colleagues had not used any NRT controls for analysis of total TERRA expression in their study. Therefore, we switched our focus to measure TERRA molecules transcribed from specific chromosome ends 17p and XpYp, which are of interest in autism [37].
As shown in Figure 3 (panels A and D), both 17p- and XpYp-TERRA expression levels showed noticeable down-regulation in male children with autism compared with the healthy controls (4.8- and 29.6- fold, respectively). The expression levels of 17p- and XpYp-TERRA in individually-matched male groups are shown in Figures 3B and 3E. It should be noted that individual TERRA expression levels do not follow the pattern for the average RTLs shown in the accompanying paper [33]. However, the male children with autism, showed a significant decrease in both 17p- and XpYp-TERRA expression levels when compared with their individually matched healthy control samples [t(7) = −3.59, p = 0.009; and t(4) = −3.97, p = 0.016, respectively].
Discussion
To our knowledge, this is the first evaluation of telomeric oxidation and transcription levels in children with autism and ASD in general. In our earlier research work, we showed that children with non-syndromic autism exhibited sexually dimorphic patterns of average RTL: while the male children with autism followed a homogeneous pattern of abnormally shortened telomeres, the autistic girls showed no sign of abnormal RTL shortening. Moreover, some female autistic cases had substantially longer RTLs than the healthy controls [33]. Here, we show that the levels of telomeric oxidized base lesions are substantially higher in autistic children than the paired siblings and healthy controls. Despite having longer telomeres, female autistic children show even higher levels of telomeric oxidation levels than their male counterparts. Also somewhat surprising, despite having shorter RTLs, autistic male children show a significant decrease in 17p- and XpYp-TERRA expression levels when compared with their individually-matched healthy controls. We will discuss some possible implications of these seemingly paradoxical observations below.
Because of their repetitive sequence and strings of GGG, telomeres are considered as “hot zones” for oxidative damage [16, 38, 39]. As outlined above, our data show substantially higher levels of telomeric oxidized base lesions in both male and female children with autism. Thus, based on the conventional view on the association of increase in oxidative stress with accelerated telomere shortening, one would expect that autistic children to exhibit abnormally shortened telomeres. While this prediction fits well with our observations for the autistic male children, it sharply contrasts with our data from the autistic females. Previously, we provided a detailed argument explaining how differences in sex biology and hormones could lay the foundations for sexual dimorphism in RTL between male and female children with autism [33]. How come then autistic girls with longer RTL have even higher levels of telomeric oxidized base lesions than their male counterparts? First, we would like to remind the reader that we do not know of the actual length distribution of different chromosome-specific telomeres. It is therefore possible that the longer average RTL in autistic girls to be the result of certain abnormally long chromosome ends not all telomeres. Second, we do not know the exact composition of the oxidized bases within the telomeres, which are composed of 25% G and 25% T for most of their length as dsDNA, followed by 50% G and 33.3% T within the small but crucial ssDNA overhang at the very end. It is true that the bacterial enzyme FPG has a preference for oxidized purines, in particular 8-oxoG and its derivatives, but it also acts on thymine glycol (Tg) and has apurinic/apyrimidinic lyase activity, introducing ssDNA breaks [40]. In other words, FPG-sensitive sites could easily contain oxidized bases other than 8-oxoG. Most notably, it has been shown that Tg residues, which are a common oxidative product of thymine [21, 41] have significant impacts on telomere length, stability, and maintenance [42–44]. Last but not least, we do not know the exact locations and origins of telomeric oxidized bases. Nor do we know the structures of the autistic telomeres and the type of protein factors interacting with the oxidized bases. The nature, distribution, and impact of different oxidized lesions at human telomeres remain largely elusive and poorly understood. Different combinations of oxidized bases and lesions at telomeres, interactions with the shelterin proteins and various DNA repair intermediates could result in contrasting telomere length and integrity outcomes [19, 20].
The most studied oxidized base 8-oxoG itself is prone to further oxidation and turning into a number of hydantoin derivatives [45]. It is also important to note, however, that even if we assume that 8-oxoG and its derivatives make up the majority of oxidized lesions at autistic telomeres, it is difficult to predict the outcome on telomere length and structure without knowing the composition and locations of these residues, their sources, the composition of other oxidized lesions, and interacting factors and chromatin context [46]. Telomeres are highly unique chromatin entities capable of forming alternate structures including T- and D-loops and G-quadruplexes (G4), interacting with specially coded histones and histone variants, sheltering proteins, and chromatin remodeling complexes. The impact and processing of the 8-oxoG and its derivatives and other oxidized lesions appear to be distinctly different at the telomeres depending on their origins, compositions, and contexts [19, 20]. Moreover, although 8-oxoG is conventionally regarded as a moderate mutagen and a biomarker of cellular oxidative stress, recent reports suggest that it functions as an “epigenetic-like” mark at chromosome internal sites, and its impact on chromatin structure and transcription is very much dependent on the genomic context within which it is situated and the interacting protein factors involved [46]. Of particular interest is the evidence showing 8-oxoG residues occurring at specific G-rich regulatory elements and their ability to increase transcription by promoter G4 formation mediated by the base excision repair machinery [46, 47]. Recent studies by Opresko, Myong, and colleagues have revealed some interesting insights into the structural impact and regulatory mechanisms of the two most common oxidized bases, 8-oxoG and Tg, within the telomeric ssDNA 3’-overhang on human telomerase activity and telomere elongation using tailored ssDNA synthesized oligonucleotides in human HEK-293 and HeLa cells [44, 48]. Interestingly, while insertion of the oxidized dGTP free nucleotides by telomerase inhibits its reverse transcriptase activity and leads to telomere shortening, the 8-oxoG residues that already exist in the 3’-overhang sequence disturb the G4 structures and simulate telomerase loading and activity leading to telomere elongation [48]. Also interestingly, in contrast to the in vitro data showing Tg acting as a strong block to DNA replication [49] and evidence suggesting that Tg residues cause telomere fragility in transgenic mice lacking the NTH1 glycosylase, the enzyme which removes Tg from telomeric sequences [43], Myong and colleagues have shown that Tg lesions in the G-rich telomeric ssDNA oligonucleotide actually stimulate telomerase binding and telomere elongation [44]. The status, compositions, and functions of oxidized lesions at the native mammalian telomeres under physiological conditions are simply not known.
A study by Jill James’s group has shown significantly higher levels of genomewide oxidative markers (including 8-oxo-dG) in the autism cerebellum and Brodmann area 22 [31], which is part of the superior temporal gyrus in the left hemisphere and plays a crucial role in speech processing [50]. We see no evidence of higher levels of telomeric oxidized bases among our nonverbal group of patients. Interestingly, Rose et al. also concluded that since they saw no correlation between the levels of oxidative biomarkers and age, this was a chronic condition in autism [31]. We, too, found no association between age and the level of telomeric oxidized bases among either the children with autism (Fig. 2B) or the healthy controls (Fig. 2A). Taken together, this brings up the possibility that the oxidative status of telomeric regions may have been set during early/critical period of neurodevelopment. Moreover, based on our data showing an inverse correlation between age and RTL [33], it can be concluded that the overall distribution of oxidized bases remains about even throughout most of the telomeric dsDNA.
Another seemingly paradoxical observation among our findings is that despite having significantly shorter average RTL in male children with autism (M.Pedram and colleagues, manuscript in preparation), they had lower levels of 17p- and XpYp-TERRA expression, which may look contrary to the expected inverse correlation between telomere length and TERRA expression [25]. However, in addition to not knowing the telomere length status at the short arms of chromosomes 17, X, and Y in our samples, we are also unaware of the structure and epigenetic status of these telomeres. It is possible that despite having shorter than normal average telomere length, the male children with autism have longer than usual telomeres at certain chromosome ends. Furthermore, having shorter TL does not necessarily translate into elevated levels of TERRA transcription. A very interesting and real-life example of this is demonstrated in Cockayne syndrome, a rare hereditary human disorder, marked by severe segmental premature aging and progressive neurological degeneration. Most Cockayne syndrome patients carry a deleterious mutation in the Cockayne syndrome group B (CSB) gene that encodes for a large nuclear protein implicated in base excision repair, transcription regulation, and chromatin remodeling [51]. It has been shown that CSB protein localizes to only a small subset of telomeres, and that CSB loss/deficit leads to telomere shortening as well as lower levels of TERRA transcription [52]. Another interesting case supporting non-random TL distribution and maintenance among individual chromosome ends has been demonstrated in patients with ICF (immunodeficiency, centromeric instability and facial anomalies) syndrome, type I [53]. It is also interesting to note that a recent study on mouse embryonic stem cells has revealed that in addition to binding to telomeres and stabilizing them in cis, TERRA molecules bind to a number of select genes and also interact with a sizable proteome including a number of very important telomeric proteins, epigenetic factors, and non-telomeric enzymes such as the RNA helicase ATRX and chromatin modifier CHD8 [24]. The gene encoding CHD8 is frequently mutated in ASD, and in fact, CHD8 appears to be a master regulator of ASD risk genes [4]. Interestingly, a recent study shows that a human CHD8 mutation leads to sexually dimorphic behavioral phenotypes (with male-preponderant deficits) and neuronal activity in transgenic mice [54].
Our observations on the sexually dimorphic patterns of telomere length in children with non-syndromic autism [33], as well as differential 17p- and XpYp-TERRA expression levels in autistic male children, could provide important insights for deciphering telomere dynamics and the underlying molecular mechanisms of childhood autism. So many interesting and important questions come to surface: What is the status of TL at specific chromosome ends in male vs. female children with autism? What is the composition of telomeric oxidized bases in male vs. female autistic children vs. that of healthy controls? What factors interact with these oxidized bases/telomeres? What are the actual impacts of 8-oxoG and Tg at native telomeres in autistic children vs. healthy controls? What is the status of DNA methylation in subtelomeric regions of male vs. female children with autism vs. that of healthy controls? What is the implication of differential chromosome-specific TERRA expressions in male children with autism? What is the status of specific TERRA expression levels in female children with autism? What would be the impact of differential TERRA expression in autistic male (and perhaps, also female) children vs. controls on telomere structure, subtelomeric gene expression, and various TERRA target genes and interacting proteins in trans (for example, the chromatin modifier CHD8)? Do TERRA molecules play a part in explaining sexual dimorphism and male-bias in autism? Answering these questions could provide important and highly valuable insights and clues for deciphering the underlying molecular mechanisms involved in the etiology of autism and also provide a better understanding of the role of telomere biology in the molecular mechanisms responsible for sex bias and shared pathophysiology in autism.
Materials and Methods
Subjects and Sampling Procedure
This work follows our previous study on children with non-syndromic autism including: 24 children (14 male and 10 female) with non-syndromic autism (61.50 ± 27.65 mo), 24 sex-, age-, and location-matched healthy control children (66.08 ± 28.08 mo) as well as 10 apparently healthy siblings (eight males and two females, 80.40 ± 55.44 mo) [33]. Briefly, children with autism were diagnosed by two qualified clinicians, a child neurologist and a specialized psychiatrist, and the diagnosis of autism was further confirmed by Autism Diagnostic Interview-Revised (ADI-R) [55]. The Majority (18 subjects) of the children with autism (including all of the males) were antipsychotic drug naïve at the time of diagnosis and saliva sampling and six of the autistic females with treatment history had extended wash out periods prior to sampling. The present study was approved by Zanjan University of Medical Sciences (ZUMS) Research Ethics Committee (ZUMS.REC.1392.97 and REC.1394.265), it was conducted in accordance with the Declaration of Helsinki, and all children had signed informed consent provided by their parents or caregivers.
Saliva sampling was done using the Oragene collection kits (DNA genoTeK Inc., Ottawa, Ontario, Canada) as described previously [35]. OG-250 saliva-collection kits were used for genomic DNA (gDNA) extraction, and RE-100 kits were used for RNA extraction. Collected saliva samples from different cities and provinces were kept at room temperature for no longer than a week and were subsequently transferred to the central laboratory at ZUMS for appropriate processing.
Telomeric 8-oxoG/Oxidized Base Measurement
The concentration of oxidized base lesions in telomeric regions was evaluated using a qPCR-based technique introduced by Nathan O’Callaghan and colleagues [34]. This method combines gDNA digestion with the formamidopyrimidine (FPG) enzyme [40], a bacterial DNA glycosylase that is involved in DNA repair through Base Excision Repair (BER) mechanism [56], and a variation of the qPCR telomere measurement technique originally devised by Richard Cawthon [57]. FPG preferably cuts the oxidized purines, in particular 8-oxoG [58]. The difference in PCR threshold cycles between DNA templates completely digested by FPG and undigested DNA is then calculated (i.e., ΔCq = Cqdigested - Cqundigested).
8-oxoG Standard Curve Preparation
Four 84-mer oligomers of 14 5'-TTAGGG-3' repeats (as forward oligomers) containing either zero, one, two, or four 8-oxoG bases each and also a reverse oligomer (an 84-mer sequence of 14 5'- CCCTAA-3' repeats) [34] (Supplementary Material, Table S1), were purchased from Jena Bioscience (Jena, Thuringia, Germany). To prepare duplex standard oligonucleotides, each of the forward oligomers were mixed with the reverse oligomer in an equal molar ratio separately, followed by incubation at 95°C for 10 min and then cooling at room temperature for 30 min. The enzymatic digestions were performed using FPG (New England BioLabs, Hitchin, UK) on each of the ds-oligonucleotides containing zero, one, two and four 8-oxoG bases (Tel 0, Tel 1, Tel 2, and Tel 4, respectively) in microtubes. 200 ng of each duplex standard oligomers were incubated with 8 units of FPG in 1X NEBuffer 1 and 1X BSA. In case of undigested oligomers as controls, FPG was replaced with ddH2O in the reaction tubes. The reactions were set up on ice, and all microtubes/reactions (digested and undigested) were then incubated at 37°C for 4.5 hr, followed by an enzyme inactivation step at 60°C for 10 min. To evaluate complete digestion, a small portion of each sample was run on a 2.5% agarose gel.
Following enzymatic digestion with FPG (New England BioLabs, Hitchin, UK), amplification of the oligomers was performed using ABI Step One Plus Real-Time PCR instrument (Applied Biosystems, Foster City, CA, USA). Primer sequences were as follows: teloF: 5′- CGGTTTGTTTGGGTTTGGGTTTGGGTTTGGGTTTGGGTT-3′ and teloR: 5′- GGCTTGCCTTACCCTTACCCTTACCCTTACCCTTACCCT-3′ [59]. qPCR reactions were set up in a final volume of 20 µL containing 40 pg of digested or undigested standard oligomers, 1X SYBR Green master mix, high ROX (Ampliqon, Odense, Denmark), and 100 nM of each primer. PCR amplifications were performed in triplicates with an initial enzyme activation step at 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec and 60°C for 30 sec. Values for ΔCq (Cq treated - Cq untreated) were then calculated and used for preparing the 8-oxoG standard curve (Supplementary Material, Fig. S1).
Measurement of Telomeric Oxidized Base Lesions in gDNA Samples
Genomic DNA (gDNA) was extracted from whole saliva samples, using the manufacturer’s protocol (OG-250; DNA genoTeK Inc., Ottawa, Ontario, Canada) with slight modification [33]. 240 ng of each gDNA sample was incubated with 12 units of FPG in 1X NEBuffer 1 and 1X BSA. To prepare undigested controls, the FPG enzyme replaced with an equivalent volume of ddH2O. Reactions were set up on ice and all microtubes/reactions (digested and undigested) were incubated at 37°C overnight followed by an enzyme inactivation step at 60°C for 10 min.
Amplification of telomeric regions in digested and undigested gDNA was performed as described in the previous section. PCR reactions were prepared in a final volume of 20 μL containing 20 ng of digested or undigested genomic DNA, 1X SYBR Green master mix, high ROX (Ampliqon, Odense, Denmark), 100 nM teloF, and 100 nM teloR primers. PCR amplifications were performed in triplicates with an initial enzyme activation step at 95°C for 10 min followed by 40 cycles of 95°C for 15 sec and 60°C for 30 sec. Values for ΔCq (Cq treated - Cq untreated) were then calculated. The number of oxidized base lesions (i.e., FPG-sensitive sites) per 1 kb of telomeric DNA were then determined using the 8-oxoG standard curve.
Telomeric RNA (TERRA) Primer Set Selection
In an attempt to determine total (i.e., not chromosome-specific) TERRA levels, we used the telomere primer sets originally introduced by Richard Cawthon [60]. Chromosome-specific TERRA primer sequences for the chromosome ends 17p and XpYp were selected from previous publications [23, 61]: 17p TERRA-F, 5′-CCACAACCCCACCAGAAAGA-3′; 17p TERRA-R, 5′- GCGCGTCCGGAGTTTG-3′; XpYp TERRA-F, 5′ - GAATCCACGGATTGCTTTGTGTACTT-3′; and XpYp TERRA-R, 5′ - CCTCAGCCTCTCAACCTGCTTGG-3′.
RNA Extraction and Analysis of TERRA Transcript Levels with RT-qPCR
Total RNA was extracted from whole saliva as described in detail previously [35]. Briefly, total RNA was reverse transcribed to Complementary DNA (cDNA) by using the PrimeScriptTM RT reagent Kit (TaKaRa Inc., Shiga, Japan) following the manufacturer’s instructions. To eliminate potential residual genomic DNA, each 500-ng RNA sample were treated with 0.5U RNase-free DNase I (1U/ µL) (Thermo Scientific, Waltham, MA, USA) prior to cDNA synthesis. TERRA transcripts originating from chromosome ends 17p and XpYp were amplified using an ABI Step One Plus Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). Each 20-µL PCR reaction for either 17p or XpYp TERRA amplification was composed of 25 ng and 50 ng of cDNA, respectively, 10 µL of Power SYBR Green master mix (Applied Biosystems, Foster City, CA, USA), and 150 nM of each forward and reverse primers. PCR amplifications were performed in triplicates with an initial enzyme activation step at 95°C for 10 min, followed by: 50 cycles of 95°C for 15 s, either 6°C (for XpYp TERRA) or 62°C (for 17p TERRA) for 30 s, and 72°C for 1 min.
Selection and validation of GAPDH and YWHAZ combination as the most stable and suitable reference genes for analysis of the salivary transcriptome in male children with non-syndromic autism has been described previously [35]. A similar analysis of the cDNAs derived from female saliva samples showed that although GAPDH and YWHAZ were still the best among a panel of nine candidate reference genes examined, they would not be sufficient for proper normalization of the RT-qPCR data obtained from female subjects (Supplementary Material, Fig. S3). Therefore, the RT-qPCR data were normalized by using the geometric mean of GAPDH and YWHAZ Cq values in male samples only. All reactions were run in triplicates. No-template (NTC) and no-reverse transcriptase (NRT) controls were included in each run.
Statistical Analysis
Distributions of telomeric oxidized base levels among different clinical groups were non-normal based on their respective histograms and/or Shapiro-Wilk tests. To achieve normal distribution, data points were first natural log transformed and then compared by using multiple linear regression. To compare the mean values of telomeric oxidized base levels among different clinical groups, groups were coded as dummy variable and incorporated into the regression analysis with adjustment for covariates including sample’s age, sex, maternal ages at birth, and birth weight. Spearman’s rho correlations were used to examine the association between telomeric oxidized base levels and ADI-R scores. To evaluate TERRA expression levels among different clinical groups, a One-sample t-test was used, and each affected subject was calibrated relative to his/her respective healthy neighboring control to minimize the effect of environment. To avoid skewness, relative expression level values were first transformed to log values to assure a normal distribution before performing the One-sample t-test.
All statistical tests were two-tailed, and the P-values below 0.05 (p < 0.05) were considered as statistically significant. Statistical analysis was performed using the SPSS software, ver. 22.0 (IBM Corp., Armonk, NY, U.S.A.). Scatter plots and bar graphs were made using GraphPad Prism, ver. 6.07 (GraphPad Software Inc. La Jolla, CA, U.S.A.).
Data Availability
The only additional relevant data referred to in the manuscript could be obtained from our earlier manuscript submitted to medRxiv shown below: Panahi, Y., Salasar Moghaddam, F., Babaei, K., Eftekhar, M., Shervin Badv, R., Eskandari, M.R., Vafaee-Shahi, M., Pezeshk, H. and Pedram, M. (2020) Sexual Dimorphism in Telomere Length in Childhood Autism. MedRxive, doi: 10.1101/2020.04.30.20074765.
Supplementary Material
A composite PDF document containing the following:
One Table (Table S1); Three Figures (Figs. S1-S3); and References (total of 3).
Funding
This work was supported by Zanjan University of Medical Sciences (ZUMS) grant numbers A-12-534-1, A-12-534-9, and it had the approval of the ZUMS Research Ethics Committee (ZUMS.REC.1392.97 and REC.1394.265).
* Please, note that ZUMS is an academic institution. The present work/research is conducted independent of the Iranian government agencies.
Conflict of interest
The authors declare no conflict of interest.
Acknowledgements
We would like to thank Reza Shervin Badv, M.D. (Dept. of Pediatric Neurology, School of Medicine, Tehran University of Medical Sciences), Mohammad Vafaee-Shahi, M.D. (Growth and Development Research Center, Institute of Endocrinology and Metabolism, Iran University of Medical Sciences) for the initial recruitment and clinical examination of some of the autistic children, and clinical follow-ups.
Footnotes
Additional Information: Statistical Expert: Hamid Pezeshk, D. Phil (author no. 5, see above/previous page)
Abbreviations
- 8-oxoG
- 8-oxo-dihydro-2'-deoxyguanosine
- ASD
- Autism spectrum disorders
- ADI-R
- Autism Diagnostic Interview-Revised
- BER
- Base Exision Repair
- Cq
- Quantification Cycle
- ds
- doublestranded
- DSM-IV-TR
- Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, Text Revision
- FPG
- Formamidopyridine DNA Glycosylase
- GAPDH
- Glyceraldehyde-3-phosphate dehydrogenase gene
- G4
- G-quadruplex
- ICD-10
- International Statistical Classification of Diseases and Related Health Problems version 10
- NRT
- No Reverse Transcriptase
- NTC
- No Template Control
- ROS
- Reactive Oxygen Species
- RT
- Reverse Transcriptase
- RTL
- Relative Telomere Length
- RT-qPCR
- RT Quantitative real-time PCR
- ss
- single-stranded
- TERRA
- Telomeric Repeat Containing RNA
- Tg
- Thymine glycol
- TL
- Telomere Length
- YWHAZ
- Tyrosine 3 monooxygenase activation, zeta plolypeptide gene