- 1Department of Medicine, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
- 2Department of Biochemistry, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
- 3Department of Morphology and Genetics, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
- 4Department of Psychobiology, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
Klinefelter syndrome (KS) displays a broad dysmorphological, endocrinological, and neuropsychological clinical spectrum. We hypothesized that the neurocognitive dysfunction present in KS relies on an imbalance in X-chromosome gene expression. Thus, the X-chromosome inactivation (XCI) pattern and neurocognitive X-linked gene expression were tested and correlated with intelligence quotient (IQ) scores. We evaluated 11 KS patients by (a) IQ assessment, (b) analyzing the XCI patterns using both HUMARA and ZDHHC15 gene assays, and (c) blood RT-qPCR to investigate seven X-linked genes related to neurocognitive development (GTPBP6, EIF2S3, ITM2A, HUWE1, KDM5C, GDI1, and VAMP7) and XIST in comparison with 14 (male and female) controls. Considering IQ 80 as the standard minimum reference, we verified that the variability in IQ scores in KS patients seemed to be associated with the XCI pattern. Seven individuals in the KS group presented a random X-inactivation (RXI) and lower average IQ than the four individuals who presented a skewed X-inactivation (SXI) pattern. The evaluation of gene expression showed higher GTPBP6 expression in KS patients with RXI than in controls (p = 0.0059). Interestingly, the expression of GTPBP6 in KS patients with SXI did not differ from that observed in controls. Therefore, our data suggest for the first time that GTPBP6 expression is negatively associated with full-scale IQ under the regulation of the type of XCI pattern. The SXI pattern may regulate GTPBP6 expression, thereby dampening the impairment in cognitive performance and playing a role in intelligence variability in individuals with KS, which warrants further mechanistic investigations.
Introduction
Klinefelter Syndrome (KS) is the most frequent chromosomal disorder in males, with a prevalence of one in 650 newborn males regarding the most common KS karyotype, 47,XXY (Samango-Sprouse et al., 2015, 2019). In addition, this syndrome is associated with endocrinological dysfunctions, in which the individuals present hypogonadism, small and firm testes, gynecomastia, infertility, and some with cognitive impairment (Haltrich, 2019). Although a broad spectrum of phenotypes has been described, KS diagnosis continues to be a challenge, with many patients being misdiagnosed or remain undiagnosed. Currently, karyotype analysis is the most common genetic test for KS diagnosis. Delays in KS diagnosis impair the treatment of patients and the associated comorbidities.
The phenotypic characteristics of patients with KS vary according to the age at diagnosis and different chromosome constitutions, mainly due to mosaicisms. Before puberty, only minor changes are noted, such as slightly reduced testicular volume or elongated lower limbs (Lanfranco et al., 2004; Winge et al., 2020). Many KS individuals present signs of developmental disorders such as dyslexia (Samango-Sprouse, 2001) and attention deficit hyperactivity disorder (Lee et al., 2011). There are also neuropsychiatric conditions, including depression, bipolar disorder, and schizophrenia (Vawter et al., 2007; Giagulli et al., 2019). Behavioral alterations or cognitive dysfunctions are usually detected only at school-age, contributing to the delayed diagnosis of KS.
The Wechsler adult intelligence scale (WAIS-III), used to evaluate overall intelligence, provides a detailed neurocognitive evaluation also for KS (Kaufman et al., 2015). This instrument evaluates the full-intelligence quotient (IQ), as well as verbal (verbal IQ—VIQ), and non-verbal (performance IQ—PIQ) quotients. Individuals with the 47,XXY karyotype usually show average or below-average intelligence performance, with IQs ranging from 80 to 109 (Rovet et al., 1995; Itti et al., 2006; Skakkebaek et al., 2014). The discrepancy between VIQ and PIQ skills, which are considered psychometrically significant variations between standardized scores according to an age range, has not been reported in the literature.
Due to the KS cytogenetic constitution, bearing two X chromosomes, KS patients undergo the X-chromosome inactivation (XCI) process (Raznahan et al., 2018). X-inactivation appears to occur as a mechanism of dosage compensation of the expression of X-linked genes when the cell has more than one X chromosome (Lyon, 1960). This molecular mechanism might play a role in KS phenotype variability (Belling et al., 2017; Raznahan et al., 2018). Moreover, most of the genes on the Y chromosome are homologous to genes on the X chromosome (Mangs and Morris, 2007). These genes are part of the pseudoautosomal regions (PAR1- short arm and PAR2—long arm) and are not subject to X inactivation. However, as widely known, aneuploid cells present unusual gene expression due to the presence of the extra chromosomes, affecting the nuclear architecture (Raznahan et al., 2018; Braun et al., 2019). In the same way, sex chromosome aneuploidies also affect the expression of X-linked genes that escape from XCI (Zhang et al., 2020). For these reasons, the expression pattern of seven genes associated with cognition (GTPBP6, EIF2S3, ITM2A, HUWE1, KDM5C, GDI1, and VAMP7)(Vawter et al., 2007; Gecz et al., 2009; Tuttelmann and Gromoll, 2010; Lagha et al., 2013; Piton et al., 2013; Belling et al., 2017; Gravholt et al., 2018) and XIST was chosen as a strategic analysis to understand, at least in part, the phenotypic variability in KS (Skakkebaek et al., 2018).
We aimed to assess whether the intelligence performance variability of KS patients is associated with the X inactivation pattern and expression profile of X-linked genes. Our study regarding XCI and IQ analyzed individuals presenting karyotype 47,XXY and considered them in statistical comparisons with male 46,XY and female 46,XX controls. We also describe the findings observed in individuals presenting the 49,XXXXY, and 48,XXYY karyotypes.
Subjects and Methods
Cytogenetic Analysis
We enrolled 27 subjects, 11 KS patients (mean age 31.7 years; SEM 4.3) and two polysomic individuals (49,XXXXY, P12; 48,XXXY, P13), from the endocrinology outpatient clinic of University’s Hospital, as indicated in Table 1, as well as seven female (mean age 37.3 years; SEM 6.0) and seven male (mean age 35.3 years; SEM 4.3) healthy controls. The 550-resolution G-banding chromosomal analysis at a minimum of 20 lymphocyte metaphases confirmed karyotypes. In addition, we performed assays to study the XCI pattern in the patients and the mothers from three patients (P3, P5, and P12) who agreed to participate in the study. All patients signed informed consents after a fully explanating of the purpose and nature of all procedures used, which the University Ethics Committee had approved (CEP/Unifesp: 1.774.406).
Table 1. Sociodemographic, clinical, and neuropsychological characterization and X-chromosome inactivation features observed in Klinefelter syndrome.
Intelligence Quotient Assessment
The intelligence quotient (IQ) of the patients was assessed using the Brazilian version of the Wechsler Adult Intelligence Scale—WAIS-III (Bowden et al., 2006). A trained neuropsychologist undertook the process. WAIS-III assesses the overall intelligence quotient (full-scale IQ), as well as intelligence associated with verbal (verbal IQ) and non-verbal (performance IQ) domains. Verbal IQ (VIQ) is a measure of verbal reasoning and acquired knowledge, including social rules. Performance IQ (PIQ), on the other hand, is established by tasks that measure non-verbal reasoning, spatial processing, and visual-motor integration (Lange, 2011). WAIS-III is composed of 13 subtests that measure specific cognitive abilities linked to four factor-based indices: verbal comprehension (vocabulary, similarities, and information), perceptual organization (picture completion, block design, and matrix reasoning), working memory (arithmetic, digit span, and letter-number sequencing), and processing speed (digit symbol, and symbol search). The main IQ and Wechsler scale indexes are expressed as standard scores (100 ± 15) and the subtests as scaled scores (10 ± 3). The IQ ranges and classifications are: ≥ 130: very superior; 120–129: superior; 110–119: high average; 90–109: average; 80–89: low average; 70–79: borderline; ≤ 69: extremely low.
RNA Isolation and Quantitative RT-PCR
The expression level of seven X-chromosome genes EIF2S3, GDI1, ITM2A, GTPBP6, HUWE1, KDM5C, VAMP7 (Supplementary Table 1), and XIST was quantified by quantitative RT-PCR in the 11 KS patients, 2 polysomic patients and 14 (male and female) controls. Total RNA was isolated from peripheral blood using Trizol® reagent, according to the manufacturer’s instructions (Invitrogen Corp., Carlsbad, CA). One microgram of total RNA from each sample was reverse transcribed into cDNA using oligo (dT) primers. cDNA was then diluted fivefold, and 1.0 μL aliquot of cDNA was used in a 12 μL PCR reactions containing SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA) and 10 μmol/L of each primer for target gene or reference gene (RPS8). PCR primers are listed in Supplementary Table 2.
RT-qPCR reactions were performed in triplicate, and the threshold cycle (Ct) was averaged (SD ≤ 1). The fold changes were calculated according to the comparative ΔΔCt method using RPS8 as reference gene, as described by Carvalheira et al. (2015). The male individuals were used as a reference group. Alternatively, for the correlation test, we represented gene expression as 1/ΔCt with no normalization to a control group. Higher 1/ΔCq values indicate higher gene expression.
Human Androgen Receptor Assay (HUMARA)
XCI assays were performed in 16 individuals, comprising 11 karyotyped 47,XXY KS patients (P1-11), two patients with karyotypes 49,XXXXY (P12) and 48,XXYY (P13), and three mothers of from the patients P3, P5 and P12. The HUMARA assay was performed according to the manufacturer’s protocol as previously described (Sisdelli et al., 2016). Two primer pairs were designed to the first exon of the androgen receptor (AR) gene, located at Xq11.2. However, when the results were non-informative for the polymorphic region of the AR gene, another primer pair was designed to amplify the polymorphic region of the zinc finger DHHC-type containing 15 (ZDHHC15) gene, located at Xq13.3, as previously described (Bertelsen et al., 2011). Primers are listed in Supplementary Table 2. All reactions were performed in triplicate, except for patient 12, which were carried out in quintuplicate.
The pattern of XCI was determined by the equation proposed by Bolduc et al. (2008), which is used to obtain the percentage of a digested allele with the greatest number of CAGn repeats. The undigested samples were used as both the control of the digested reactions and to calculate the number of CAGn repeats. The XCI was considered random when the ratio presented values between 50 and 65%, indeterminate between 66 and 73%, skewed when between 74 and 90%, and extremely skewed between 90 and 100% (Bolduc et al., 2008).
Statistical Analysis
Intellectual performance measures were analyzed descriptively. Intelligence quotient scores and relative gene expression levels meeting parametric criteria were expressed as means, and the two-tailed Student’s t-test was applied for two group comparisons, while one-way ANOVA followed by Bonferroni test was performed for multiple group comparisons. Non-parametric data were expressed as median, and Kruskal-Wallis followed by Dunn test was used for multiple group comparisons. Pearson’s correlation coefficient (r) was calculated to evaluate the correlation between two numeric variables. GraphPad© PRISM (Version 8.0; San Diego, CA) was used for all statistical analyses. A P-value of < 0.05 was considered statistically significant. Polysomic patients were not considered in the statistical comparison of gene expression between KS (47,XXY) patients and controls due to their phenotypic and molecular variability (Braun et al., 2019).
Results
Klinefelter Patients Present a Broad Intelligence Quotient Variability
Applying the Wechsler Adult Intelligence Scale with age-related assessment criteria, we verified that 11 KS patients with karyotype 47,XXY displayed total IQ values varying from 66 to 109, which correspond to extremely low to average levels a broad range of variation of the intelligence performance in our cohort. Sociodemographic, clinical and neuropsychological characterization, including the schooling, clinical comorbidity, VIQ (Vocabulary, Similarities, Arithmetic, Digit Span, Information and Comprehension), and PIQ (Figure Completion, Coding, Block Design, Matrix Reasoning, and Picture Arrangement), are described in Table 1. There was a difference of at least nine points between the PIQ and VIQ scores in five (P3, P4, P8, P10, and P11) individuals, with a better performance in the former. Regarding polysomic participants, the patient with the variant karyotype 49,XXXXY (P12) presented a severe intellectual disability, score < 1st percentile in the Raven’s Progressive Matrices Scale, corresponding to an IQ ≤ 55, while the patient with karyotype 48,XXYY (P13) showed a low average level IQ.
The Skewed X Inactivation Pattern in the Blood Is Associated With a Better Intelligence Quotient Score
To assess whether the XCI pattern is associated with the intelligence phenotype of KS patients, we used HUMARA and ZDHHC15 assays. Among the KS patients, seven (64%) demonstrated random X inactivation (RXI), while four (36%) presented skewed X inactivation (SXI) (Table 1). The mean ages at KS diagnosis and XCI pattern analysis were, respectively, 16.4 years (SEM 1.9) and 24.4 years (SEM 2.8) for patients with RXI and 38.8 years (SEM 7.9) and 44.5 (SEM 7.5) for patients with SXI pattern. In addition, the two polysomic patients with 49,XXXXY and 48,XXYY karyotypes showed RXI (P12) and SXI (P13), respectively (Table 1).
The mean IQ of the seven patients with RXI was at borderline levels (Mean = 76.86), while for the four patients with SXI, it was average (Mean = 93.75), demonstrating that in this cohort, SXI is present in KS patients with higher IQ values (Figure 1, p = 0.0451). Four patients with borderline IQ levels (P2, P3, P4, and P5), ranging from 70 and 77, presented RXI (Table 1). One patient (P1) with intellectual deficiency (IQ 66) presented RXI. On the other hand, four patients with low to average IQ levels (P8, P9, P10, and P11), ranging from 80 to 109, presented an SXI pattern. The polysomic patients also followed the association between XCI and IQ. P12 (49,XXXXY) presented an IQ ≤ 55 and RXI, while P13 (48,XXYY) presented an IQ of 80 and SXI (Table 1).
Figure 1. Intelligence quotient scores in Klinefelter syndrome according to the X-chromosome inactivation pattern. RXI patients (n = 7) show a lower IQ mean than SXI patients (n = 4). The mean IQ from patients with RXI was at borderline levels (Mean = 76.86), while in the SXI patients, it was average (Mean = 93.75). Data are presented as mean ± SEM. Statistical analysis used the t-test (*p = 0.0451). IQ, intelligence quotient. RXI, random X-chromosome inactivation; SXI, skewed X-chromosome inactivation.
Based on this broad IQ variability, we further tested whether seven X-chromosome genes, previously associated with neurocognitive function, would be altered in expression in KS patients and whether the expression of these genes was correlated with XCI patterns.
The X-Inactivation Pattern Interferes in GTPBP6 Expression, a Gene Related to Intellectual Performance
In respect of the general expression analysis of the eight selected X-linked genes, GTPBP6, EIF2S3, ITM2A, HUWE1, KDM5C, XIST, GDI1, and VAMP7, only XIST (Xq13.2) as expected, and the pseudoautosomal gene GTPBP6 (Xp22.33/Yp11.32) showed differential expression (Figure 2). XIST expression was higher in KS patients with RXI than in male controls (p = 0.0090, Figure 2A). KS patients with the RXI pattern presented higher GTPBP6 expression when compared with both male and female controls (p = 0.0059). Interestingly, the expression of GTPBP6 in KS patients with SXI did not differ from that observed for male and female controls (Figure 2B). Also, the four KS patients with the SXI pattern presented general expression level profiles similar to the male controls (Figure 2). The gene expression levels for the polysomic patients are presented in Supplementary Table 3 for comparison.
Figure 2. Neurocognitive-related X-chromosome linked gene expression analysis tested in Klinefelter Syndrome patients in comparison to healthy controls. Eight X-linked genes, (A) XIST, (B) GTPBP6, (C) EIF2S3, (D) KDM5C, (E) VAMP7, (F) GDI1, (G) ITM2A, and (H) HUWE1, were evaluated. Genes that escape from X-chromosome inactivation and are biallelically expressed are indicated in the gray area. Relative expression (fold change) is represented in relation to males (XY). Note that in (A), the XIST expression was higher in KS patients with RXI in comparison to male controls (p = 0.0090), while in (B), the GTPBP6 expression was higher in KS patients with RXI than in both female and male controls (p = 0.0059). In (A,C,D,E,F), data are presented as the median and the statistical test used was Kruskal-Wallis followed by the Dunn test. In (B,G,H), data are presented as the mean, and statistical analysis used was one-way ANOVA followed by the Bonferroni test. *p < 0.05; **p < 0.01. The number of individuals considered per group is indicated in Supplementary Table 3.
We further investigated whether GTPBP6 expression would correlate to intelligence assessment. Although Pearson’s correlation test suggested a negative correlation trend between GTPBP6 expression and full-scale IQ (r = –0.5045), it was not statistically significant (p = 0.0787) among all KS and polysomic patients (Figure 3A). However, when considering IQ 80 as the standard minimum reference, five patients (45%) with IQ < 80 presented the RXI pattern and showed higher GTPBP6 expression (Mean = 2.946) when compared with both female and male controls (p = 0.0012) (Figure 3B). On the other hand, all the SXI patients presented IQ ≥ 80, with four classified as SXI, and showed GTPBP6 expression levels that did not differ from controls (Figure 3B).
Figure 3. Intelligence quotient variability in association with the X-chromosome inactivation pattern and GTPBP6 gene expression evidenced in Klinefelter syndrome patients. (A) Pearson’s correlation between GTPBP6 expression and full-scale intelligence quotient in Klinefelter patients. GTPBP6 expression was calculated as 1/ΔCq using RPS8 as the reference gene. Higher 1/ΔCq values indicate greater GTPBP6 expression. Linear regression was calculated (slope = –0.01385). Klinefelter and polysomic patients with RXI and SXI are represented with white and black circles, respectively. Correlation r = –0.5045, p = 0.0787. Polysomic subjects are indicated. (B) Peripheral blood GTPBP6 gene expression analysis in KS patients with IQ below the standard average of 80 (n = 5) in comparison to both male (XY) and female (XX) controls. All KS patients with IQ < 80 presented a random X-inactivation pattern (RXI). Relative expression (fold change) is represented in relation to males. Note that GTPBP6 expression was higher in KS patients with IQ < 80 in comparison to both male and female controls (p = 0.0012). Data are presented as mean ± SEM and statistical analysis used One-way ANOVA followed by Bonferroni test (**p < 0.01). (C) Hypothetical model for understanding of the correlation between intelligence quotient, X-inactivation pattern and altered GTPBP6 expression in Klinefelter patients. This model suggests that the intelligence quotient in KS individuals may reflect the X-inactivation pattern. Also, when presenting RXI, the GTPBP6 expression may be elevated and the intelligence quotient is low. On the other hand, when presenting the SXI pattern, the GTPBP6 expression is low and the intelligence quotient is average. Our data indicate that the SXI pattern may downregulate GTPBP6 expression, leading to a better cognitive performance in KS patients. Color schemes for the XCI patterns show 50% of the cells with the same Xi representing RXI, and 75% of the cells with the same Xi representing SXI, following the classification in Bolduc et al. (2008). IQ, intelligence quotient; RXI, random X-inactivation; SXI, skewed X-inactivation.
Discussion
Although there has been significant progress in understanding genetics concerning neurodevelopmental disorders over the last two decades, little is known about the etiopathogenesis of cognitive impairment in 47,XXY patients. Typically, each second X-chromosome is silenced in half of the cells randomly, but this pattern may be skewed when preferential inactivation of an X-chromosome occurs. Therefore, we hypothesized that the X-inactivation pattern might interfere with cognitive phenotype development in KS. To address this topic, we assessed the XCI patterns as well as the expression of seven X-linked genes associated with neurocognitive function and their correlation with intelligence performance in KS patients.
Our data endorsed the variability in intelligence performance in the typical 47,XXY KS cohort, as reported previously (Samango-Sprouse, 2001; Tartaglia et al., 2010; Brandenburg-Goddard et al., 2014; Skakkebaek et al., 2014; Samango-Sprouse et al., 2018). We also verified that VIQ was significantly lower than PIQ in KS patients with average IQ, which has been known that lower VIQ is associated with language impairments, commonly observed in KS (Gravholt et al., 2018). Higher PIQ scores indicate that non-verbal skills (e.g., analyzing and synthesizing visual stimuli) are well developed. In other words, it seems that verbal tasks (word definition, verbal comprehension, concept formation, and the ability to define concepts orally) may be incredibly challenging for some KS patients. It is important to note that intellectual disability (IQ < 70) is rare in KS patients, and that verbal intellectual skills (VIQ) are usually more affected than non-verbal (PIQ) skills (Marco and Skuse, 2006; Skakkebaek et al., 2014; Melogno et al., 2018). A similar intellectual profile (VIQ < PIQ) is also observed in autism spectrum disorder (Ankenman et al., 2014).
Intellectual disability was detected in one and a borderline level IQ in four cases. Although intellectual variability in KS has already been described (Samango-Sprouse, 2001; Tartaglia et al., 2010; Skakkebaek et al., 2014; Brandenburg-Goddard et al., 2014; Samango-Sprouse et al., 2018), reports of disability are not usual. A high neurocognitive phenotypic variability implies that KS condition is frequently underdiagnosed. Consequently, hypogonadic individuals might not receive early hormonal replacement treatment, which may affect cognitive performance (Samango-Sprouse et al., 2013, 2015). It is also possible that the wide variation of IQ identified in our sample may be because the patients referred to our service, which is a tertiary reference center, usually present more complex clinical manifestations.
In KS, comorbidities such as an increased risk for obesity, diabetes, and seizures have been reported (Belling et al., 2017). A predisposition to metabolic diseases has been associated with a possible biological function of the methylation changes found in KS, as demonstrated by gene expression enrichment analysis (Skakkebaek et al., 2018). Diabetes, epilepsy, schizophrenia, or obesity was present in almost half (46%) of our patients. Some of these clinical conditions can negatively impact cognitive development and may have contributed to the intellectual variability (Al-Shazely and Al-Khaligy, 2014). A high rate of epilepsy has been reported in patients with other sex chromosome aneuploidies (Tatum et al., 1998; Verri et al., 2017). Noteworthy, P1 (47,XXY) and P12 (49,XXXXY) presented complex phenotypic variation besides epilepsy, including autoimmune disease (Churg-Strauss syndrome) and refractory controlling seizures (P1).
Many studies have attempted to establish criteria for early diagnosis of KS (Belling et al., 2017), with chromosome karyotype analysis being the best method of diagnostic confirmation so far. The phenotypic variability seems to result from both genetic and environmental factors, although the relative contribution of these events remains unclear. Regarding genetic influences, variability appears to be associated with a gene dosage effect due to the additional X chromosome. It can be the result of X-inactivation (Iitsuka et al., 2001), uniparental disomy of the sex chromosomes (Bruining et al., 2010), and a positive family history of learning disabilities (Samango-Sprouse et al., 2013, 2014).
As is widely known, XCI occurs in cells that contain more than one X chromosome, which increases the chance of survival of male aneuploid KS individuals compared to those with other forms of aneuploidy. XCI occurs due to the high expression of the non-coding RNA XIST, upregulated in 46,XX subjects (Vawter et al., 2007; Minks et al., 2008; Belling et al., 2017). This genetic phenomenon interferes with the X-linked genes’ expression and results in symptoms that may vary in individuals according to the XCI pattern. It has been suggested that SXI in KS patients occurs as a common genetic event (Iitsuka et al., 2001). However, it has also been shown that 10–50% of KS patients present SXI (Amos-Landgraf et al., 2006; Bojesen et al., 2011; Mehta et al., 2012; Kinjo et al., 2020). Our data show similar results to these studies. Seven patients (64%) presented the RXI pattern, and four (36%) the SXI pattern. It is known that 5–20% of seemingly healthy 46,XX normal subjects also display SXI (Belmont, 1996). It has been reported that females who carried out X-linked recessive lethal allele, or an allele that result in a proliferative advantage for cells, show skewing of X-inactivation (Naumova et al., 1996). These observations suggest that the XCI pattern in somatic cells may be influenced by the X-chromosome genetic background, promoting a cellular selection favoring the best genetic balance (Gartler and Sparkes, 1963). Another molecular mechanism involved in the XCI pattern is aging (Sharp et al., 2000; Amos-Landgraf et al., 2006). The prevalence of SXI is higher in older than in younger women (Sandovici et al., 2004). Patients with RXI pattern had a mean age at KS diagnosis and XCI pattern analysis lower than with SXI. These data may reflect that the RXI pattern could enhance the severity of the symptoms, while in KS patients with SXI, the clinical features could be mild, including those related to neurocognitive performance.
Another genetic feature related to the extra sex chromosome observed in the KS is the trisomy of pseudoautosomal genes (Zhang et al., 2013; Carrel and Brown, 2017; Balaton et al., 2018). These events would explain clinical symptoms and the susceptibility to comorbidities. The excess of gene expression from one or more X-linked genes that escape inactivation could be responsible for the cognitive phenotype in psychotic disorders diagnosed in KS (DeLisi et al., 2005; Wallentin et al., 2016), the reason by which we analyzed the expression of seven X-linked genes involved in neurodevelopment collected from peripheral blood cells. Indeed, a robust RNA-seq-based study has demonstrated that an XCI pattern from the blood is predictive of X-inactivation status in other tissues, such as the brain (Tukiainen et al., 2017), although these tissues have a different embryonic origin.
We demonstrated that GTP Binding Protein 6 (GTPBP6) gene presented differential expression among the other neurocognitive X-linked genes in KS studied group. GTPBP6 is expressed in the fetal and adult brain and has been negatively correlated with language functioning (Gianfrancesco et al., 1998; Vawter et al., 2007). Noteworthy, large-scale human exome sequencing data have found mutations in other GTPases (Rho family), causing X-linked intellectual disability (Piton et al., 2013). The overexpression of this gene was reported in lymphoblastoid cell lines derived from 11 males with KS, among 14 differentially expressed X-chromosome genes in a KS group (Vawter et al., 2007). Our results complement these data, demonstrating that higher GTPBP6 expression is observed in KS patients with RXI but not with SXI. This observation suggests that epigenetic events under the X-chromosome inactivation process may somehow interfere with the regulation of the genes involved in cognition in KS.
Moreover, our data suggested that GTPBP6 expression was inversely correlated with IQ scores in KS patients. On the other hand, KS patients with GTPBP6 expression similar to controls presented SXI patterns and higher IQ scores than those with RXI. This finding indicates that the SXI pattern somehow stabilizes GTPBP6 expression at levels similar to controls, which may function as a protective factor for cognitive development. Furthermore, it is known that in individuals with sex-chromosome aneuploidies even genes that escape X-inactivation are under the influence of the epigenetic events participating in this molecular mechanism, such as neighboring DNA methylation and changes in topologically associating domain (TAD) boundaries (Raznahan et al., 2018; Braun et al., 2019). Thus, these data corroborate the idea of GTPBP6 as a non-inactivated gene reflecting the XCI pattern. GTPBP6 is among the most sensitive genes to sex chromosome dosage (Raznahan et al., 2018), which may be reflected in the range of cellular and molecular strategies to regulate its expression in KS. These underlying molecular mechanisms in humans are still poorly understood and may involve complex signaling pathways that are yet to be investigated.
It is important to note that SXI may reflect selection against uniparental disomy and gene dosage, favoring cell survival and viability. Studies using murine models demonstrated a global increase in DNA methylation on the X chromosome from the paternal origin compared to the maternal one, implicating higher expression of X genes in XY than XX cells (Golden et al., 2019). In KS patients, neurocognitive functions were reported to be influenced by the parental origin of the extra X chromosome (Navarro-Cobos et al., 2020). Considering the data presented here, it is likely that the SXI pattern may be a molecular mechanism that guarantees cellular homeostasis. Thus, the effects of sex chromosome trisomy can be dampened, and, consequently, the genes present in the pseudoautosomal region, such as GTPBP6, would have their expression regulated by preferential inactivation of the X-chromosome.
Therefore, dysregulation of X-linked genes seems to cause an epigenetic imbalance that explains most of the changes observed in sex chromosome aneuploidies, particularly in KS. Our findings are summarized in Figure 3C, in which we propose a hypothetical model for understanding the correlation between X-inactivation pattern and altered GTPBP6 expression, being inversely related to intelligence performance. This model implicates the SXI pattern as a protective molecular strategy mitigating the deleterious effects of unbalanced gene expression related to the extra X chromosome. These data shed light on the contribution of the XCI pattern and GTPBP6 expression to the neurocognitive phenotype observed in KS patients and may help to elucidate the molecular mechanisms associated with neurodevelopment dysfunction in X-linked aneuploid disorders.
Data Availability Statement
The datasets analyzed for this study will be found in the Unifesp Repository (http://repositorio.unifesp.br/).
Ethics Statement
The studies involving human participants were reviewed and approved by the Universidade Federal de São Paulo (Unifesp) Ethics Committee N° 1.774.406. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
LS assisted patients and performed all neuropsychological tests. LF, AV, JS, and IK performed molecular and bioinformatics analysis. MM and GC reviewed cytogenetics and gene expression data. LS, LF, CM, GC, and MD designed the study, performed, and reviewed statistics. GC and MD equally conceived and coordinated the study, edited, and reviewed the manuscript. All authors discussed the data and contributed to writing the final manuscript.
Funding
This work was supported by grant 312746/2017-2 from the Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq and by grants 17/13513-7 and 18/03511-0 from the São Paulo Research Foundation—FAPESP. JS was granted FAPESP scholarships (2017/07053-3 and 2018/22763-0), and funding 2021/00684-3 to MD.
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank reported patient and their mothers for their cooperation in this study, Patricia Monteagudo, Andrea Mathez, Joao Ferro, and Julia Oliveira for the clinical follow-up of patients; Janete Cerutti and Rui Maciel for lab support. Conselho Nacional de Desenvolvimento Científico e Tecnológico—CNPq and São Paulo Research Foundation—FAPESP for financial sponsorship.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.724625/full#supplementary-material
References
Al-Shazely, S. M., and Al-Khaligy, H. (2014). Intelligence quotient in children with epilepsy. Curr. Res. Neurosci. 4, 10–17. doi: 10.3923/crn.2014.10.17
Amos-Landgraf, J. M., Cottle, A., Plenge, R. M., Friez, M., Schwartz, C. E., Longshore, J., et al. (2006). X chromosome-inactivation patterns of 1,005 phenotypically unaffected females. Am. J. Hum. Genet. 79, 493–499. doi: 10.1086/507565
Ankenman, K., Elgin, J., Sullivan, K., Vincent, L., and Bernier, R. (2014). Nonverbal and verbal cognitive discrepancy profiles in autism spectrum disorders: influence of age and gender. Am. J. Intellect. Dev. Disabil. 119, 84–99. doi: 10.1352/1944-7558-119.1.84
Balaton, B. P., Dixon-McDougall, T., Peeters, S. B., and Brown, C. J. (2018). The eXceptional nature of the X chromosome. Hum. Mol. Genet. 27, R242–R249.
Belling, K., Russo, F., Jensen, A. B., Dalgaard, M. D., Westergaard, D., Rajpert-De Meyts, E., et al. (2017). Klinefelter syndrome comorbidities linked to increased X chromosome gene dosage and altered protein interactome activity. Hum. Mol. Genet. 26, 1219–1229. doi: 10.1093/hmg/ddx014
Belmont, J. W. (1996). Genetic control of X inactivation and processes leading to X-inactivation skewing. Am. J. Hum. Genet. 58, 1101–1108.
Bertelsen, B., Tumer, Z., and Ravn, K. (2011). Three new loci for determining X chromosome inactivation patterns. J. Mol. Diagn. 13, 537–540. doi: 10.1016/j.jmoldx.2011.05.003
Bojesen, A., Hertz, J. M., and Gravholt, C. H. (2011). Genotype and phenotype in Klinefelter syndrome–impact of androgen receptor polymorphism and skewed X inactivation. Int. J. Androl. 34, e642–e648.
Bolduc, V., Chagnon, P., Provost, S., Dube, M. P., Belisle, C., Gingras, M., et al. (2008). No evidence that skewing of X chromosome inactivation patterns is transmitted to offspring in humans. J. Clin. Invest. 118, 333–341. doi: 10.1172/JCI33166
Bowden, S. C., Weiss, L. G., Holdnack, J. A., and Lloyd, D. (2006). Age-related invariance of abilities measured with the wechsler adult intelligence scale–III. Psychol. Assess. 18, 334–339. doi: 10.1037/1040-3590.18.3.334
Brandenburg-Goddard, M. N., van Rijn, S., Rombouts, S. A., Veer, I. M., and Swaab, H. (2014). A comparison of neural correlates underlying social cognition in Klinefelter syndrome and autism. Soc. Cogn. Affect. Neurosci. 9, 1926–1933. doi: 10.1093/scan/nst190
Braun, R., Ronquist, S., Wangsa, D., Chen, H., Anthuber, L., Gemoll, T., et al. (2019). Single chromosome aneuploidy induces genome-wide perturbation of nuclear organization and gene expression. Neoplasia 21, 401–412. doi: 10.1016/j.neo.2019.02.003
Bruining, H., van Rijn, S., Swaab, H., Giltay, J., Kates, W., Kas, M. J., et al. (2010). The parent-of-origin of the extra X chromosome may differentially affect psychopathology in Klinefelter syndrome. Biol. Psychiatry 68, 1156–1162. doi: 10.1016/j.biopsych.2010.08.034
Carrel, L., and Brown, C. J. (2017). When the Lyon(ized chromosome) roars: ongoing expression from an inactive X chromosome. Philos. Trans. R. Soc. Lond. B Biol. Sci. 372:20160355. doi: 10.1098/rstb.2016.0355
Carvalheira, G., Nozima, B. H., and Cerutti, J. M. (2015). MicroRNA-106b-mediated down-regulation of C1orf24 expression induces apoptosis and suppresses invasion of thyroid cancer. Oncotarget 6, 28357–28370. doi: 10.18632/oncotarget.4947
DeLisi, L. E., Maurizio, A. M., Svetina, C., Ardekani, B., Szulc, K., Nierenberg, J., et al. (2005). Klinefelter’s syndrome (XXY) as a genetic model for psychotic disorders. Am. J. Med. Genet. B: Neuropsychiatr. Genet. 135, 15–23. doi: 10.1002/ajmg.b.30163
Gartler, S. M., and Sparkes, R. S. (1963). The lyon-beutler hypothesis and isochromosome x patients with turner syndrome. Lancet 2:411. doi: 10.1016/s0140-6736(63)93096-4
Gecz, J., Shoubridge, C., and Corbett, M. (2009). The genetic landscape of intellectual disability arising from chromosome X. Trends Genet. 25, 308–316. doi: 10.1016/j.tig.2009.05.002
Giagulli, V. A., Campone, B., Castellana, M., Salzano, C., Fisher, A. D., de Angelis, C., et al. (2019). Neuropsychiatric aspects in men with Klinefelter syndrome. Endocr. Metab. Immune. Disord. Drug Targets 19, 109–115. doi: 10.2174/1871530318666180703160250
Gianfrancesco, F., Esposito, T., Montanini, L., Ciccodicola, A., Mumm, S., Mazzarella, R., et al. (1998). A novel pseudoautosomal gene encoding a putative GTP-binding protein resides in the vicinity of the Xp/Yp telomere. Hum. Mol. Genet. 7, 407–414. doi: 10.1093/hmg/7.3.407
Golden, L. C., Itoh, Y., Itoh, N., Iyengar, S., Coit, P., Salama, Y., et al. (2019). Parent-of-origin differences in DNA methylation of X chromosome genes in T lymphocytes. Proc. Natl. Acad. Sci. U.S.A. 116, 26779–26787. doi: 10.1073/pnas.1910072116
Gravholt, C. H., Chang, S., Wallentin, M., Fedder, J., Moore, P., and Skakkebaek, A. (2018). Klinefelter syndrome: integrating genetics. Neuropsychology, an endocrinology. Endocr. Rev. 39, 389–423. doi: 10.1210/er.2017-00212
Haltrich, I. (2019). Chromosomal aberrations with endocrine relevance (turner syndrome). Exp. Suppl. 111, 443–473. doi: 10.1007/978-3-030-25905-1_20
Iitsuka, Y., Bock, A., Nguyen, D. D., Samango-Sprouse, C. A., Simpson, J. L., and Bischoff, F. Z. (2001). Evidence of skewed X-chromosome inactivation in 47,XXY and 48,XXYY Klinefelter patients. Am. J. Med. Genet. 98, 25–31. doi: 10.1002/1096-8628(20010101)98:1<25::aid-ajmg1015>3.0.co;2-x
Itti, E., Gaw Gonzalo, I. T., Pawlikowska-Haddal, A., Boone, K. B., Mlikotic, A., Itti, L., et al. (2006). The structural brain correlates of cognitive deficits in adults with Klinefelter’s syndrome. J. Clin. Endocrinol. Metab. 91, 1423–1427. doi: 10.1210/jc.2005-1596
Kaufman, A. S., Raiford, S. E., and Coalson, D. L. (2015). Intelligent Testing with the WISC-V. New York, NY: Wiley, 832.
Kinjo, K., Yoshida, T., Kobori, Y., Okada, H., Suzuki, E., Ogata, T., et al. (2020). Random X chromosome inactivation in patients with Klinefelter syndrome. Mol. Cell. Pediatr. 7:1. doi: 10.1186/s40348-020-0093-x
Lagha, M., Mayeuf-Louchart, A., Chang, T., Montarras, D., Rocancourt, D., Zalc, A., et al. (2013). Itm2a is a Pax3 target gene, expressed at sites of skeletal muscle formation in vivo. PLoS One 8:e63143. doi: 10.1371/journal.pone.0063143
Lanfranco, F., Kamischke, A., Zitzmann, M., and Nieschlag, E. (2004). Klinefelter’s syndrome. Lancet 364, 273–283. doi: 10.1016/S0140-6736(04)16678-6
Lange, R. T. (2011). “Performance IQ,” in Encyclopedia of Clinical Neuropsychology, eds J. S. Kreutzer, J. De Luca, and B. Caplan (New York, NY: Springer), 1907–1908. doi: 10.1007/978-0-387-79948-3_1066
Lee, N. R., Wallace, G. L., Clasen, L. S., Lenroot, R. K., Blumenthal, J. D., White, S. L., et al. (2011). Executive function in young males with Klinefelter (XXY) syndrome with and without comorbid attention-deficit/hyperactivity disorder. J. Int. Neuropsychol. Soc. 17, 522–530. doi: 10.1017/s1355617711000312
Lyon, M. F. (1960). A further mutation of the mottled type in the house mouse. J. Hered. 51, 116–121. doi: 10.1093/oxfordjournals.jhered.a106965
Mangs, H. A., and Morris, B. J. (2007). The human pseudoautosomal region (par): origin, function and future. Curr. Genomics 8, 129–136. doi: 10.2174/138920207780368141
Marco, E. J., and Skuse, D. H. (2006). Autism-lessons from the X chromosome. Soc. Cogn. Affect. Neurosci. 1, 183–193. doi: 10.1093/scan/nsl028
Mehta, A., Malek-Jones, M., Bolyakov, A., Mielnik, A., Schlegel, P. N., and Paduch, D. A. (2012). Methylation-Specific PCR allows for fast diagnosis of X chromosome disomy and reveals skewed inactivation of the X chromosome in men with Klinefelter syndrome. J. Androl. 33, 955–962. doi: 10.2164/jandrol.111.016030
Melogno, S., Pinto, M. A., Orsolini, M., and Tarani, L. (2018). Beyond the literal meaning of words in children with Klinefelter syndrome: two case studies. Brain Sci. 8, 171. doi: 10.3390/brainsci8090171
Minks, J., Robinson, W. P., and Brown, C. J. (2008). A skewed view of X chromosome inactivation. J. Clin. Invest. 118, 20–23. doi: 10.1172/JCI34470
Naumova, A. K., Plenge, R. M., Bird, L. M., Leppert, M., Morgan, K., Willard, H. F., et al. (1996). Heritability of X chromosome–inactivation phenotype in a large family. Am. J. Hum. Genet. 58, 1111–1119.
Navarro-Cobos, M. J., Balaton, B. P., and Brown, C. J. (2020). Genes that escape from X-chromosome inactivation: potential contributors to Klinefelter syndrome. Am. J. Med. Genet. C Semin. Med. Genet. 184, 226–238. doi: 10.1002/ajmg.c.31800
Piton, A., Redin, C., and Mandel, J. L. (2013). XLID-causing mutations and associated genes challenged in light of data from large-scale human exome sequencing. Am. J. Hum. Genet. 93, 368–383. doi: 10.1016/j.ajhg.2013.06.013
Raznahan, A., Parikshak, N. N., Chandran, V., Blumenthal, J. D., Clasen, L. S., Alexander-Bloch, A. F., et al. (2018). Sex-chromosome dosage effects on gene expression in humans. Proc. Natl. Acad. Sci. U.S.A. 115, 7398–7403. doi: 10.1073/pnas.1802889115
Rovet, J., Netley, C., Bailey, J., Keenan, M., and Stewart, D. (1995). Intelligence and achievement in children with extra X aneuploidy: a longitudinal perspective. Am. J. Med. Genet. 60, 356–363. doi: 10.1002/ajmg.1320600503
Samango-Sprouse, C. (2001). Mental development in polysomy X Klinefelter syndrome (47,XXY; 48,XXXY): effects of incomplete X inactivation. Semin. Reprod. Med. 19, 193–202.
Samango-Sprouse, C., Stapleton, E., Chea, S., Lawson, P., Sadeghin, T., Cappello, C., et al. (2018). International investigation of neurocognitive and behavioral phenotype in 47,XXY (Klinefelter syndrome): predicting individual differences. Am. J. Med. Genet. A 176, 877–885. doi: 10.1002/ajmg.a.38621
Samango-Sprouse, C., Stapleton, E. J., Lawson, P., Mitchell, F., Sadeghin, T., Powell, S., et al. (2015). Positive effects of early androgen therapy on the behavioral phenotype of boys with 47,XXY. Am. J. Med. Genet C Semin. Med. Genet. 169, 150–157. doi: 10.1002/ajmg.c.31437
Samango-Sprouse, C. A., Counts, D. R., Tran, S. L., Lasutschinkow, P. C., Porter, G. F., and Gropman, A. L. (2019). Update on the clinical perspectives and care of the child with 47,xxy (Klinefelter syndrome). Appl. Clin. Genet. 12, 191–202. doi: 10.2147/tacg.s180450
Samango-Sprouse, C. A., Sadeghin, T., Mitchell, F. L., Dixon, T., Stapleton, E., Kingery, M., et al. (2013). Positive effects of short course androgen therapy on the neurodevelopmental outcome in boys with 47,XXY syndrome at 36 and 72 months of age. Am. J. Med. Genet. A 161A, 501–508. doi: 10.1002/ajmg.a.35769
Samango-Sprouse, C. A., Stapleton, E. J., Mitchell, F. L., Sadeghin, T., Donahue, T. P., and Gropman, A. L. (2014). Expanding the phenotypic profile of boys with 47, XXY: the impact of familial learning disabilities. Am. J. Med. Genet. A 164A, 1464–1469. doi: 10.1002/ajmg.a.36483
Sandovici, I., Naumova, A. K., Leppert, M., Linares, Y., and Sapienza, C. (2004). A longitudinal study of X-inactivation ratio in human females. Hum. Genet. 115, 387–392.
Sharp, A., Robinson, D., and Jacobs, P. (2000). Age-and tissue-specific variation of X chromosome inactivation ratios in normal women. Hum. Genet. 107, 343–349. doi: 10.1007/s004390000382
Sisdelli, L., Vidi, A. C., Moysés-Oliveira, M., Di Battista, A., Bortolai, A., Moretti-, Ferreira, D., et al. (2016). Incorporation of 5-ethynyl-2’-deoxyuridine (EdU) as a novel strategy for identification of the skewed X inactivation pattern in balanced and unbalanced X-rearrangements. Hum. Genet. 135, 185–192. doi: 10.1007/s00439-015-1622-x
Skakkebaek, A., Gravholt, C. H., Rasmussen, P. M., Bojesen, A., Jensen, J. S., Fedder, J., et al. (2014). Neuroanatomical correlates of Klinefelter syndrome studied in relation to the neuropsychological profile. Neuroimage Clin. 4, 1–9. doi: 10.1016/j.nicl.2013.10.013
Skakkebaek, A., Nielsen, M. M., Trolle, C., Vang, S., Hornshoj, H., Hedegaard, J., et al. (2018). DNA hypermethylation and differential gene expression associated with Klinefelter syndrome. Sci. Rep. 8:13740. doi: 10.1038/s41598-018-31780-0.
Tartaglia, N., Cordeiro, L., Howell, S., Wilson, R., and Janusz, J. (2010). The spectrum of the behavioral phenotype in boys and adolescents 47,XXY (Klinefelter syndrome). Pediatr. Endocrinol. Rev. 8, 151–159.
Tatum, W. O., Passaro, E. A., Elia, M., Guerrini, R., Gieron, M., and Genton, P. (1998). Seizures in Klinefelter’s syndrome. Pediatr. Neurol. 19, 275–278. doi: 10.1016/s0887-8994(98)00055-1
Tukiainen, T., Villani, A.-C., Yen, A., Rivas, M. A., Marshall, J. L., Satija, R., et al. (2017). Landscape of X chromosome inactivation across human tissues. Nature 550, 244–248. doi: 10.1038/nature24265
Tuttelmann, F., and Gromoll, J. (2010). Novel genetic aspects of Klinefelter’s syndrome. Mol Hum. Reprod. 16, 386–395. doi: 10.1093/molehr/gaq019.
Vawter, M. P., Harvey, P. D., and DeLisi, L. E. (2007). Dysregulation of X-linked gene expression in Klinefelter’s syndrome and association with verbal cognition. Am. J. Med. Genet. B Neuropsychiatr. Genet. 144B, 728–734. doi: 10.1002/ajmg.b.30454
Verri, A., Carmen, D., Anna, C., Federica, C., Anna, M., and Chiara, C. (2017). Variability in cognitive behavioral phenotypes in Klinefelter syndrome (KS) and other sex chromosomal aneuploidies (SCAs). Andrology 6, 1–9. doi: 10.4172/2167-0250.1000175
Wallentin, M., Skakkebaek, A., Bojesen, A., Fedder, J., Laurberg, P., Ostergaard, J. R., et al. (2016). Klinefelter syndrome has increased brain responses to auditory stimuli and motor output, but not to visual stimuli or Stroop adaptation. Neuroimage Clin. 11, 239–251. doi: 10.1016/j.nicl.2016.02.002
Winge, S. B., Soraggi, S., Schierup, M. H., Rajpert-De Meyts, E., and Almstrup, K. (2020). Integration and reanalysis of transcriptomics and methylomics data derived from blood and testis tissue of men with 47,XXY Klinefelter syndrome indicates the primary involvement of sertoli cells in the testicular pathogenesis. Am. J. Med. Genet. C Semin. Med. Genet. 184, 239–255. doi: 10.1002/ajmg.c.31793
Zhang, X., Hong, D., Ma, S., Ward, T., Ho, M., Pattni, R., et al. (2020). Integrated functional genomic analyses of Klinefelter and Turner syndromes reveal global network effects of altered X chromosome dosage. Proc. Natl. Acad. Sci. U.S.A. 117, 4864–4873. doi: 10.1073/pnas.1910003117
Zhang, Y., Castillo-Morales, A., Jiang, M., Zhu, Y., Hu, L., Urrutia, A. O., et al. (2013). Genes that escape X-inactivation in humans have high intraspecific variability in expression, are associated with mental impairment but are not slow evolving. Mol. Biol. Evol. 30, 2588–2601. doi: 10.1093/molbev/mst148
Keywords: X-chromosome inactivation, Klinefelter syndrome, X-linked genes, intelligence quotient, GTPBP6, gene expression
Citation: Simonetti L, Ferreira LGA, Vidi AC, Souza JSd, Kunii IS, Melaragno MI, Mello CBd, Carvalheira G and Dias da Silva MR (2021) Intelligence Quotient Variability in Klinefelter Syndrome Is Associated With GTPBP6 Expression Under Regulation of X-Chromosome Inactivation Pattern. Front. Genet. 12:724625. doi: 10.3389/fgene.2021.724625
Received: 15 June 2021; Accepted: 24 August 2021;
Published: 20 September 2021.
Edited by:
Enrique Medina-Acosta, State University of the North Fluminense Darcy Ribeiro, BrazilReviewed by:
Srimonta Gayen, Indian Institute of Science (IISc), IndiaHala El-Bassyouni, National Research Centre, Egypt
Cíntia Barros Santos-Rebouças, Rio de Janeiro State University, Brazil
Copyright © 2021 Simonetti, Ferreira, Vidi, Souza, Kunii, Melaragno, Mello, Carvalheira and Dias da Silva. 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: Magnus R. Dias da Silva, bXJkc2lsdmFAdW5pZmVzcC5icg==
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share senior authorship