Monoamine oxidase B gene variants associated with attention deficit hyperactivity disorder in the Indo-Caucasoid population from West Bengal

Background Attention deficit hyperactivity disorder (ADHD) is characterized by symptoms of inattention, excessive motor activity and impulsivity detected mostly during childhood. These traits are known to be controlled by monoamine neurotransmitters, chiefly dopamine, serotonin and norepinephrine. Monoamine oxidase A (MAOA) and B (MAOB), two isoenzymes bound to the outer membrane of mitochondria, are involved in the degradation of monoamines and were explored for association with ADHD in different ethnic groups. In the present study, few exonic as well as intronic MAOB variants were analyzed in ADHD probands (N = 150) and ethnically matched controls (N = 150) recruited following the Diagnostic and Statistical Manual for Mental Disorders-4th edition (DSM-IV). Appropriate scales were used for measuring the behavioural attributes. Gene variants were analyzed by amplification of target sites followed by DNA sequencing and data obtained were analyzed by population based statistical methods. Results Out of 34 variants present in the analyzed sites, only seven functional variants, rs4824562, rs56220155, rs2283728, rs2283727, rs3027441, rs6324 and rs3027440, were found to be polymorphic. rs2283728 ‘C’ (P = 3.45e-006) and rs3027440 ‘T’ (P = 0.02) alleles showed higher frequencies in ADHD probands as compared to controls. rs56220155 ‘A’ (P = 0.04) allele and ‘GA’ (P = 0.04) genotype showed higher frequencies in the male and female ADHD probands respectively as compared to sex-matched controls. Analysis of pairwise linkage disequilibrium revealed striking differences between probands and controls. Haplotype analysis revealed significantly higher occurrence of different haplotypes in the ADHD probands while some haplotypes were detected in the controls only. Higher scores for conduct problems were found to be associated with rs56220155 ‘A’ (P = 0.05) allele in the male ADHD probands. Multifactor dimensionality reduction analysis showed independent as well as interactive effects of polymorphic variants which were more robust in the male probands. Conclusions Since all the polymorphic variants analyzed were functional, it may be inferred that MAOB gene variants are contributing to the etiology of ADHD in the Indo-Caucasoid population from eastern India which merits further in depth analysis. Electronic supplementary material The online version of this article (doi:10.1186/s12863-016-0401-6) contains supplementary material, which is available to authorized users.


Background
Attention deficit hyperactivity disorder (ADHD) is an etiologically complex behavioural disorder. Major symptoms include persistent age-inappropriate hyperactivity and impulsivity, sometimes in association with inattention [1], leading to impairments in academic performances as well as social life [2,3]. Worldwide ADHD is highly prevalent and boys are more prone to the disorder than girls [4,5]. In India also, ADHD is quite prevalent and diagnosed more frequently in boys than girls [6][7][8]. Co-morbidity with other psychiatric disorders is common and in such condition, impairment is more [9,10].
Being a multi-factorial genetic disorder with around 76 % heritability [11,12], genetics is believed to play significant role in the etiology of ADHD [13,14]. Neurotransmitters like dopamine, serotonin, and norepinephrine regulate all vital behavioural attributes and studies on candidate genes involved in the regulation of these neurotransmitters [15,16] revealed associations between altered dopaminergic transmission and behavioural as well as cognitive deficits in various populations [17,18].
Monoamine oxidase A (MAOA) and B (MAOB) are mitochondrial outer membrane bound isoenzymes, catalyzing oxidative deamination of neurotransmitters like dopamine, serotonin, norepinephrine, and other neuromodulators like benzylamine, phenylethylamine (PEA), tyramine and tryptamine in the brain as well as peripheral tissues [19,20]. The two isoenzymes differ in substrate specificity [20,21]; while MAOA preferentially oxidizes serotonin and norepinephrine, MAOB prefers benzylamine and PEA. MAOB activity in human increases with age [22] and is selectively inhibited by low concentration of deprenyl; but in high concentration, the selectivity is lost [23]. In the human brain, MAOB is the key enzyme degrading dopamine [24][25][26] and subcortical regions exhibit higher MAOB activity [22]. MAOB was hypothesized to control impulsivity, attention and vulnerability to ADHD by degrading dopamine [24,25], which is the major factor responsible for regulating behaviour and cognitive function [18]. Further, MAOB knockout mice showed high level of PEA in the brain as well as an increased reactivity to stress and other behavioural alterations [20]. A correlation between platelet MAOB activity and sensation seeking as well as impulsiveness have also been reported [27,28]. Platelet MAOB activity was used as a marker for psychic behaviour though it was not evident whether platelet MAOB activity was correlated with brain MAOB activity or not [29].
Several MAOA variants have showed association with ADHD in various populations, including the Indo-Caucasoid population [30,31]. MAOB gene variants have also been studied, though the numbers of variants investigated were few and the data obtained were inconsistent [24,28,[32][33][34][35][36][37]. Since genes encoding for MAOA and MAOB are located on the X-chromosome [24], we hypothesized that these genes may have a role in the sex bias of ADHD and our earlier study revealed a biased maternal transmission of MAOA variants to the male probands [30]. In this study, for the first time few MAOB variants were explored for association with ADHD and its associated phenotypic traits in the Indo-Caucasoid population.

Subject recruitment
Sample size of 150 was determined statistically [38] considering 8 % prevalence of ADHD in this population [6]. ADHD cases (126 males and 24 females) were recruited by child psychiatrist and clinical psychologist following the Diagnostic and Statistical Manual of Mental Disorders-4 th edition (DSM-IV) criteria [1]. 73.34 % of the recruited cases were of the combined subtype, while hyperactive-impulsive and inattentive subtypes were of 13.33 % each. Mean age of the ADHD cases was 7.69 ± 2.99 years (range 3 to 18 years). Psychological evaluation was done through -The revised Conners' Parent Rating Scale (CPRS-R) [39] and Wechsler Intelligence Scale for Children >5 years [40]/Developmental Screening Test [41] for children < 5 years for the inattentionhyperactivity level and intelligent quotient (IQ) status respectively. DSM-IV score for oppositional defiant disorder (ODD) and Parental Account of Children's Symptoms (PACS) score for conduct problems were also used for assessing the traits in ADHD probands. Patients with any other neuropsychiatric disorders, mental retardation (IQ ≤ 70) including Down syndrome and Fragile-X syndrome, pervasive developmental disorder were excluded from the study.
The control group comprised of 150 ethnically matched healthy individuals (125 males and 25 females) assessed by the same psychometric evaluation procedure. Mean age of the control individuals was 18.41 ± 8.78 years (range 3 to 28 years).

Genotyping
Peripheral blood collected from the study participants was used for genomic DNA preparation using the standard protocol [42]. The target regions (detailed in Additional file 1) were amplified via polymerase chain reaction using primers (provided in Additional file 2) designed in the lab using the Primer3 software [43]. Applied Biosystems 3130 Genetic Analyzer with 98.5 % base calling accuracy and Read Length of upto 950 bp was used for sequence analysis of the amplicons using Big Dye v 3.1 chemistry and Sequencing Analysis Software, v 5.2 (Additional file 2). Chromatograms were also analyzed manually and mis-spaced letters/double peaks were investigated carefully for genotyping. For identification of heterozygous SNPs, >25 % base calling was accepted. Function of polymorphic variants was analyzed in silico using the is-rSNP [44].

Statistical analyses of data
To test the Hardy-Weinberg equilibrium [46], genotypic counts of only female ADHD probands and ethnically matched female controls were used since the MAOB gene is X-linked [24] and it is still unclear whether the male hemizygotes and female homozygotes are equivalent or not [47]. Allelic and genotypic association analyses for individual polymorphism as well as haplotype analysis were carried out using the UNPHASED v 3.1.5 [48] and correction for multiple testing was done while running the UNPHASED at 1000-fold iteration. To examine genotypic association, only female cases and female controls were considered. To calculate the power of the tests showing significant association, Piface version 1.72 [49] was used. Online odds ratio calculator [50] was used to calculate the odds ratio (OR). Pairwise linkage disequilibrium (LD) between the variants was measured using the Haploview program version 4.2 [51].

Analysis of interaction between the variants
Interaction between the variants or epistasis was analyzed by the Multifactor Dimensionality Reduction (MDR) program [52] through a 4-step process using the case-control data set. In the first step, using filter methods, interesting polymorphisms were selected from the pool of possible candidates through entropy-based measures of information gain (IG) for each individual polymorphism/attribute (i.e. main effects) and each pairwise combination of attributes (i.e. two way interaction effects) [53]. In the second step, a new multilocus attribute, which capture interaction information is constructed using previously selected polymorphisms in conjunction with constructive induction algorithm. Thus multilocus genotypes were pooled into high-risk and low-risk groups, effectively reducing the dimensionality of the attributes from multiple dimensions to one dimension. In the third step, the new multilocus attribute constructed in the previous step was evaluated using a machine learning method (i.e. naive Bayes classifier, based on probability theory). In the final step, an interaction circle graph, using the entropy estimates from step 1, was depicted by the program. Interaction circle graph comprised of a node for each attribute (i.e. polymorphism) with pairwise connections between them. The percentage of entropy (i.e. information gain or IG) by each polymorphism was visualized on each node, while the IG for each pairwise combination of polymorphisms was visualized on each connection. Thus, the independent main effects of each polymorphism were quickly compared to the interaction effect between them. Positive entropy values indicated synergy, while negative entropy values indicated redundancy [53]. All these analyses were implemented in the open-source MDR software package version 2.0 beta 8.4.

Association of alleles with phenotypic traits
Based on the CPRS-R, 'T scores' , ranging between 38 and 90, were obtained for ADHD probands. DSM-IV scores (ranging between 0 and 36) for assessing ODD trait and PACS scores (ranging between 0 and 90) for assessing conduct problems were also obtained. Male probands were sub-grouped based on the presence/ absence of the derived allele for each variant and distribution pattern of behavioural scores in each of the two comparing groups was checked using the Kolmogorov-Smirnov normality test [54]. Equality of variances was also checked using two sample F-test [55]. Allelic association with behavioural scores was analyzed using the Student's t-test [56] only when the variables (i.e. behavioural scores) showed a normal distribution and variances were equal. In other conditions, nonparametric test such as Mann-Whitney test [56] was performed. As the number of female probands was limited, comparative analysis on behavioural scores and genotypes was not performed for this group. Chi-square (p-value)

Results
Out of 34 variants localized in the investigated regions, only seven, rs4824562, rs56220155, rs2283728, rs2283727, rs3027441, rs6324 and rs3027440, were found to be polymorphic in the studied population (Additional file 1). In silico analysis revealed that all these variants have potential regulatory function (Additional file 3).
Frequencies of rs4824562 'G' , rs3027441'C' , rs6324 'T' and rs3027440 'C' alleles in the IND population revealed significant differences as compared to several world populations, while the distribution pattern of rs4824562, rs6324 and rs3027440 matched with populations from the South Asia (Table 1). On the other hand, rs56220155 exhibited statistically significant difference even with populations from South Asia (Table 1, BEB and GIH, P   (Table 1). No information was available in the 1000 Genomes Project Phase 3 (32) database for rs2283728. Case-control comparative analysis revealed statistically significant higher occurrence of rs2283728 'C' (P = 1.21e-005; power = 99 %) and rs3027440 'T' (P = 0.04; power = 53 %) alleles in the probands ( Table 2, Additional file 4). Stratification based on gender revealed statistically significant higher occurrence of rs2283728 'C' (P = 3.45e-006; power = 99.63 %), rs3027440 'T' (P = 0.02; power = 66 %) and rs56220155 'A' (P = 0.04; power = 54 %) alleles in the male probands in comparison to the male controls ( Table 2, Additional file 4). All the variants followed the HWE in the female subjects (Additional file 5). Statistically significant higher occurrence of rs56220155 'GA' genotype (P = 0.04; OR = 3.92; 95 % confidence interval (CI) = 1.28-11.95; power = 62 %) was also observed in the female probands as compared to the female controls (Additional file 5). Rest of the investigated variants did not show any biased occurrence (Additional files 4 and 5).
Pairwise LD analyses showed strikingly varied patterns (Fig. 1, Additional file 6); LD of the ADHD group was significantly different from that of the control group (Fig. 1a, b). Subjects stratified on the basis of gender revealed significant difference between the male subjects and sex-matched controls (Fig. 1c, d). Striking differences in the pattern of the LDs were also noticed between female ADHD probands and female control group (Fig. 1e, f).
The C-T haplotype, formed between rs2283728 and rs3027440 having strong LD (D' 0.84, r 2 0.58), exhibited significant higher occurrence (P = 1.07e-008; OR = 3.6) in the probands (Table 3). Significant higher occurrence of the A-C-T haplotype of rs56220155-rs2283728-rs3027440 (P = 1.99e-009; OR = 4.52) was also noticed in the male probands (Table 3). Haplotype analysis for all the seven variants (rs4824562-rs56220155-rs2283728-rs2283727-rs3027441-rs6324-rs3027440) showed statistically significant higher occurrence of four haplotypes (A-A-C-C-T-C-T, A-G-C-C-T-C-T, A-G-T-A-C-T-C and G-A-C-C-T-C-T) in the probands (P ≤ 0.04; OR > 2) ( Table 3). Stratified analysis based on gender, also revealed statistically significant higher occurrence of three haplotypes (A-A-C-C-T-C-T, A-G-T-A-C-T-C and G-A-C-C-T-C-T) in the male probands (P ≤ 0.005; OR > 3) and nominally significant higher occurrence of the A-G-C-C-T-C-T haplotype in the female probands (P = 0.055; OR = 8.23) as compared to sex-matched controls ( Table 3).

Discussion
MAOB is a key enzyme in the human brain, modulating oxidation of dopamine [24][25][26] as well as benzylamine, PEA, tyramine, and tryptamine [19][20][21]. Previous genetic association studies on MAOB revealed inconsistent findings in different populations. To find out the role of MAOB variants in the etiology of eastern Indian ADHD probands, we used a four step approach. Initially, allelic and genotypic frequencies of MAOB variants were analyzed by population-based methods to identify risk variants in the ADHD probands. Then, LD between the studied variants was analyzed to understand whether these variants are working independently or in a pairwise clubbed manner. Next, we verified independent main effects and epistatic effects of MAOB variants using case-control data set. Finally, to identify the relevance of these gene variants in disease associated symptoms, we analyzed association between alleles and behavioural attributes.
Alleles/genotypes of the eastern Indian control population studied in the present investigation resembled the South Asian ancestral population for most of the variants. However, allelic frequencies of rs2283727 and rs3027440 differed from the Gujrati population while rs56220155 differed from both Bengali from Bangladesh and Gujrati Indian from Houston, principally due to an increase in the minor allele frequency. We have earlier reported vast difference in the frequency of alleles in the Indian population [57]. Whether this drift in allelic frequencies is conferring any specific advantage is a matter of conjecture at the moment and merits further investigation in a large number of samples belonging to each ethnic group.
Earlier investigators reported significant positive association of MAOB gene variants with ADHD in the Spanish probands [37], while in the Irish [33] and Czech [28] population no association was noticed. International Multicentre ADHD Gene project, with Caucasian subjects from 12 specialized centres in eight different countries, also failed to notice any association between ADHD and MAOB gene variants, including rs4824562, rs6324 and rs3027440 [34]. In the present investigation, out of 34 variants only 7 were identified to be polymorphic and in silico analysis revealed that all can potentially regulate MAOB transcription. Four variants, rs56220155, rs2283728, rs2283727 and rs3027441, were analyzed for the first time for association with ADHD. Variants like rs4824562 and rs3027440 were previously studied in the European Caucasoid probands [34], while rs6324 was studied in ADHD probands belonging to Han Chinese [24,35] as well as European Caucasoid [34] populations; in the Han Chinese population, significant positive [24] as well as negative associations [35] were reported for rs6324. This site failed to show any positive association with ADHD in the Indo-Caucasoid population.
rs2283728 and rs3027440 showed allelic as well as haplotypic associations with ADHD in the Indo-Caucasoid population. rs56220155 showed genotypic association in the female probands. All these three variants also showed association in the male ADHD probands. Haplotypes consisting of all the seven variants, including the above three, showed significant association with the disorder. It can be hypothesized from the present observation that these  b Population based analysis for female subjects. The graphical interaction model describe the percentage of entropy (i.e. information gain or IG) in case-control status that is explained by each factor (i.e. gene variant) or two-way interaction. Two-way interactions between factors are depicted by line accompanied by a percent of entropy explained by that interaction. Values inside large boxes on nodes indicate information gain (IG) of individual/independent main effect of each polymorphic variant, whereas values inside small boxes between nodes exemplify IG of pairwise combination/interactive effects of respective variants. Positive IG values, between the nodes, indicate the synergistic interactions; whereas negative IG values indicate the redundancy between the respective nodes/variants. Schematic coloration represents a continuum from synergy to redundancy. The red lines represent a high degree of synergy. The orange lines represent moderate synergy. The golden yellow lines accompanied by a positive percent of entropy represent minimal synergy, whereas the golden yellow lines accompanied by a negative percent of entropy represent minimal redundancy. Green line represents high redundancy three variants independently as well as in combinations may play an important role in ADHD. We have noticed a striking difference in the LD pattern of ADHD probands and controls; pairwise all variants were in strong LD in the ADHD probands as compared to the controls. In absence of ethnic differences, recruitment of related individuals, and consanguineous marriage, the observed difference in the LD pattern may suggest a lower rate of recombination between the studied variants in the probands which facilitates generation of risk haplotypes associated with the disease etiology. Epistasis analysis showed significant pairwise synergistic interactive effects of most of the variants in the male ADHD probands. In female probands the interactive effects were very less. However, the number of female probands was limited and further exploration on the matter is desired before reaching into any conclusion. Allele 'A' of rs56220155 was associated with high conduct problems, as measured by the PACS score, in the male ADHD probands. The 'A' allele of rs56220155 also showed statistically significant higher occurrence in the male ADHD probands in comparison to the male controls. It can be inferred from these observations that MAOB has a significant role in the etiology of ADHD.

Conclusions
Our investigation for the first time revealed association of rs56220155 and rs2283728 with ADHD. rs3027440, previously reported to have no association in the Euro-Caucasoid ADHD subjects [34], also revealed positive association in the Indo-Caucasoid population. rs2283727 and rs3027441, in strong LD with rs2283728 and rs6324 respectively, were investigated for the first time in ADHD probands and statistical analysis failed to show any association in the studied population. The observed difference in allelic/genotypic association in the present study could be attributed to difference in allelic frequencies since the IND population revealed a different allelic distribution pattern as compared to other ancestral ethnic groups from other parts of the world. Stratified analysis revealed gross difference in the LD pattern of male and female ADHD probands as compared to sexmatched controls possibly be due to absence of recombination between the sites in the probands, thus creating a block conferring risk of ADHD. We have also noticed higher frequencies of rs56220155 'A' , rs2283728 "C" and rs3027440 "T" alleles in the male probands. Male probands exhibiting conduct problems also showed higher frequency of rs56220155 'A'. Whether this is really a male specific effect is a matter of conjecture at the moment since the major limitation of the present study was the low number of female subjects investigated. The high odd's ratio (>2) observed for a few association analyses could also be attributed to the limitation in sample number. We may infer from the data obtained that further investigation on a large cohort of samples belonging to different ethnic groups is warranted to validate our observation in the Indo-Caucasoid population from the eastern India.