- Biomedical and Pharmaceutical Research Unit, Department of Biomedical Sciences, College of Health Sciences, QU Health, Qatar University, Doha, Qatar
Rheumatic heart disease (RHD) is a heritable inflammatory condition characterized by carditis, arthritis, and systemic disease. Although remaining neglected, the last 3 years has seen some promising advances in RHD research. Whilst it is clear that RHD can be triggered by recurrent group A streptococcal infections, the mechanisms driving clinical progression are still poorly understood. This review summarizes our current understanding of the genetics implicated in this process and the genetic determinants that predispose some people to RHD. The evidence demonstrating the importance of individual cell types and cellular states in delineating causal genetic variants is discussed, highlighting phenotype/genotype correlations where possible. Genetic fine mapping and functional studies in extreme phenotypes, together with large-scale omics studies including genomics, transcriptomics, epigenomics, and metabolomics, are expected to provide new information not only on RHD but also on the mechanisms of other autoimmune diseases and facilitate future clinical translation.
Introduction
Acute rheumatic fever (ARF) is a systemic autoimmune disease that can develop after upper respiratory tract infection with group A Streptococcus (GAS). Recurrent rheumatic fever and associated sustained and abnormal inflammatory responses damage the heart valves, leading to rheumatic heart disease (RHD) (Figure 1) (3). RHD remains a common cause of acquired heart disease in young adults and children in many developing countries and, to a lesser extent, in developed countries (4, 5). Globally, GAS upper respiratory tract infections are only exceeded by HIV, tuberculosis, and malaria in terms of consequent morbidity and mortality (1). In 2015, an ~33 million RHD cases worldwide were reported, leading to over 300,000 deaths and over 10 million disability-adjusted life years (4). Despite these figures, the research community and health policymakers generally neglect RHD, which received the lowest research funding relative to global disease burden over the last 5 years (6). Moreover, RHD is rarely discussed in the context of autoimmune disease in the literature.
Figure 1. A schematic showing the possible pathogenic pathways giving rise to rheumatic heart disease (RHD) after group A streptococcal (GAS) infection. (A) RF/RHD is thought to be initiated by infection with rheumatogenic strains of GAS. A strong familial predisposition and the fact that only 60% of ARF patients develop RHD (1) indicate that the disease may only develop in those who are genetically predisposed. Genetic factors reported in RHD are mainly in immune response components including innate immunity genes. (B) Recurrent GAS infection leads to the development of autoreactive T cells and the production of cross-reactive autoantibodies. Recently, inflammasome activation has been shown to play an important role in the development of autoreactive T cells through persistent release of IL-1β (2). (C) Recruitment of autoreactive immune cells and cross-reactive autoantibodies to the valve interstices cause tissue damage. Recurrent and prolonged inflammation cause ongoing tissue damage, tissue fibrosis, and calcification. Figure was generated using Biorender.com.
ARF can present with several different symptoms a few weeks after a GAS episode. Rheumatic fever can be diagnosed based on the Jones criteria, developed in 1944, and then revised in 1992 and in 2015 by the American Heart Association (AHA) (Table 1) (7). ARF can be grouped into three main symptom constellations: musculoskeletal, cardiac, and neuropsychiatric. Fever is the most common presentation (>90% of patients), followed by arthritis (60–80% of patients). However, the most serious manifestation of ARF is carditis, which causes scarring and sometimes thickening of the valves, leading to permanent valvular damage that contributes to the morbidity and mortality of the disease. Some ARF patients have asymptomatic episodes, hence exposing the patients to more severe damage (8). Anywhere between 35 and 76% of ARF patients develop carditis (9, 10), predominantly left-sided, with the mitral valve affected in almost 100% of younger patients, followed by the aortic valve in 20–30% of cases and the tricuspid valve in 15–40% of cases (10). The right heart valves are nearly always affected in association with mitral or aortic disease and pulmonary valve is rarely affected. It remains unclear why particular valves are targeted in different individuals. Sydenham's chorea is the neuropsychiatric manifestation of ARF and is characterized by involuntary movements of the limbs and trunk (11). In ~20% of patients, chorea is the first disease manifestation, and while it may resolve over a few months in some patients it can last for several years, with mitral stenosis a common late complication (12).
Table 1. Revised Jones criteria for rheumatic fever diagnosis [adopted from (7)].
The molecular mechanisms underpinning the progression of ARF to RHD are poorly understood. An autoimmune response initiated by molecular mimicry is the most commonly proposed mechanism (13), since the GAS M protein and human proteins such as cardiac myosin show structural similarities, leading to antibody cross-reactivity and cellular immune responses to human tissue (14). Cross-reactivity may not be the mechanism of progression, since autoantibodies have been detected in healthy individuals and cardiac myosin is not expressed at the cell surface (15). Approximately 60% of ARF patients develop RHD (1), and although there is a proven association between GAS infection and RHD, the triggered autoimmune process in RHD can occur autonomously after removing the stimulus (10), suggesting that after initiation of the autoimmune response via molecular mimicry, host factors, most likely genetic factors, play an important role in disease progression in susceptible individuals. This hypothesis is supported by multiple strands of evidence. First, monozygotic twins are highly concordant for the disease compared to dizygotic twins (44 vs. 12%, respectively), although these findings date to studies reported between 1933 and 1964 (16). Second, a strong familial predisposition for ARF has been noted, and the estimated heritability based on twin studies is 0.6 (95% CI 0.41–0.81), similar to the heritability of other autoimmune genetic disorders (17). Third, GAS outbreaks (18) and endemic exposure to GAS in indigenous populations (19) show that only 3–6% of individuals develop ARF, indicating that the disease may only develop in those who are genetically predisposed (20). Finally, a number of susceptibility loci have now been reported, mainly in immune-related genes and including both human leucocyte antigen (HLA) (Table 2) and non-HLA genes (35). This article summarizes the known genetic associations with RHD and discusses their role in disease pathogenesis and their potential for risk prediction and prevention.
Table 2. Some of the HLA class II alleles reported to be associated with ARF/RHD in different populations, illustrating the complex landscape of the disease.
Search Strategy and Selection Criteria
References for this Review were identified through searches of PubMed with the search terms “rheumatic heart,” “wide genome association,” “single nucleotide polymorphism,” and “onset” from 1990 until July, 2020. Only papers published in English were reviewed. The final reference list was generated on the basis of originality and relevance to the broad scope of this Review. We would like to stress that this is not a systematic Review.
The Genetics of RHD
RHD is not a single gene disorder rather a complex multifactorial disease arising from the interaction between multiple genetic factors, GAS infection, and the resulting uncontrolled inflammatory response (Figure 1). Each of these genetic factors individually confers a small risk and explains a small fraction of disease heritability. Over the last few decades, genome-wide association studies (GWAS) have been the major tool used to identify genetic loci that predispose to different complex diseases. Interestingly, many of the identified susceptibility loci in RHD are clustered in key immunological pathways and are shared with other autoimmune diseases, showing evidence of natural positive selection due to their contribution to fighting infectious diseases during evolution (36–38). However, only three GWAS have been performed in RHD patients, with most other associations detected in candidate gene studies of small numbers of genes known to be functionally linked to RHD. Although important genes have been shown to be associated with RHD using the candidate-based approach, it has been difficult to successfully validate the results in different populations.
Screening the Genome: HLA and IGH Associations
Data from the three GWASs performed in RHD patients detected associations in HLA and immunoglobulin heavy (IGH) loci (Table 3). The Pacific Islands Rheumatic Heart Disease Genetics Network performed the first GWAS in RHD in 2017 (39), which examined 1006 indigenous people in different Oceanian countries using the 300k Illumina HumanCore platform, which is considered a low-density GWAS chip. The authors justified the use of low-density GWAS by predicting that fewer variants would be needed due to linkage disequilibrium in Oceanian populations stretching over greater distances than in any other population. A novel susceptibility signal was identified in the IGH locus; the IGHV4-61*02 allele was associated with a 1.4-fold increased disease risk (odds ratio (OR) 1.43; 95% CI 1.27–1.61; p = 4.1 x 10−9), regarded as a small to moderate effect (41). The IGH locus is a difficult region to study, as knowledge of its polymorphisms is limited and it is poorly tagged in current arrays; indeed, only 16 genotype variants were included in this study out of the entire 1,255 kilobase locus. However, the discovery of genome-wide significance of a specific gene segment is a very important step forward to understanding the immunogenetic background of RHD.
The second genome-wide association study focused on the indigenous population of Australia using the 550K Illumina Infinium HumanCoreExome platform (21). By screening 398 RHD cases, the authors identified HLA-DQA1 (rs9272622; OR = 0.897; p = 1.86 × 10−7, for protective allele C) as the strongest association. HLA-DQB1*0601 (OR = 1.07; p = 4.06 × 10−4) was also identified as a risk allele, and HLA-DQA1*0301 (OR = 0.92; p = 2.71 × 10−4) was identified as protective. Two risk haplotypes were identified [HLA-DQA1*0101-DQB1*0503 (OR 1.44) and HLA-DQA1*0103-DQB1*0601 (OR 1.27)] and one protective haplotype [HLA-DQA1*0301-DQB1*0402 (OR 0.3)]. The study also identified HLA-DRB1*0803 (OR = 1.06; p = 0.005) as a susceptibility locus.
The third GWAS study was of a combined Indian (510 cases) and Fijian (162 cases) population, which was then confirmed in a third population from the UK (150 cases) (40). The strongest signal (rs201026476) in this study was found in the 3' untranslated region of the PBX2 gene sited in HLA class III. The second signal was in a regulatory single nucleotide polymorphism (SNP) in HLA-DQB1, and the third strongest signal was in HLA class I (HLA-B). The signal in PBX2 was believed to have a regulatory element, and the authors speculated that it may influence the expression of other HLA class III genes including complement genes such as C2 and C4. These genes and complement pathways have been proposed to play an important role in the inflammatory response in RHD patients (42). Interestingly, after conditional analysis, the HLA class III signal remained associated with RHD, indicating that the effect of rs201026476 is independent of class I and class II.
The association between RHD and HLA alleles has been reported in other small candidate marker studies. Table 2 summarizes HLA class II alleles reported to be associated with ARF/RHD in different populations. Interestingly an HLA-DQA1*0104 and DQB1*0503 association with RHD was reported in Japanese RHD patients (30). To our knowledge, HLA-DRB1*08 has not previously been reported, but HLA-DRB1*07 was reported to be a risk allele in small studies of Turkish, Pakistani, and Latvian populations (24, 29, 31). Interestingly, in an Egyptian population, the DRB1*07-DQA1*02 haplotype was associated with the mitral valve lesion subgroup of the disease (27). Other reported risk alleles include DRB1*01, DRB1*13, DRB1*15, and DQB1*08 [reviewed in (43)]. Some other HLA alleles have been reported to be protective of the disease, with examples including DRB1*04-DQA1*03 in Turkish patients and DRB1*06-DQB1*06 in Latvian patients (24, 29). In the Australian indigenous population, the HLA-DQA1*0301 and HLA-DQA1*0301-DQB1*0402 haplotypes were identified as protective.
These heterogeneous results demonstrate the challenge of translating genetic findings into clinically meaningful associations. Intriguingly, the genome-wide study on the indigenous Australian population did not replicate the IGHV4-61*02 association found in the Oceanian populations. However, both populations are indigenous, and genetic heterogeneity may have affected the results. Alternatively, the variable results may represent different exposures to different GAS strains with varying interactions with protective or risk-associated HLA alleles. Additionally, epistasis between HLA alleles and non-HLA variants may explain some of the discrepancies (44, 45). In general, the HLA and IGH associations may support the molecular mimicry hypotheses as one of the mechanisms of pathogenesis in RHD, albeit in indigenous populations.
Candidate Genes
The candidate gene approach relies on forming biological hypotheses to predict the potential gene or locus for investigation. However, when the exact mechanism of the disease is unknown, this method may not identify its full genetic basis. Nevertheless, candidate gene case-control approaches have identified important genetic loci associated with different autoimmune diseases and have shed light on their molecular mechanisms. Relatively few candidate gene studies have been reported in RHD compared to other autoimmune diseases, and most associations remain invalidated in independent populations [reviewed in (46)]. In this section, we discuss some of these reported loci.
As post-infection autoimmune diseases involve a persistent inflammatory reaction, tumor necrosis factor-α (TNF-α) was one of the earliest genes to be studied in RHD, showing strong associations in different populations. A meta-analysis of seven studies including 735 RHD cases and 926 controls from seven different populations showed that the TNFA-308A allele was associated with the disease (OR 3; 95% CI 1.2–10.6; p = 0.02) (47). However, it is unknown whether this association is due to linkage disequilibrium (LD), the correlation between nearby variants, with HLA class I and II. HLA class II and class III are in high LD, so fine mapping of this region is important to identify causative variants. Interleukin (IL)10 is another immune modulatory gene that has been studied in RHD patients from different populations and other autoimmune diseases. Three promoter polymorphisms in IL10 (-1082A/G, -829C/T, and -592C/A) have been particularly well examined. The IL10-ACC promoter haplotype, associated with increased IL-10 production, has been shown to be a protective locus in RHD (OR = 0.6, 95% CIs 0.4–0.9; p = 0.01) (48). However, a meta-analysis of three different populations found that this locus shows only a non-statistically significant trend toward protection against RHD (OR = 0.70, 95% CIs 0.47–1.05; p = 0.08) (49). The angiotensin I-converting enzyme gene insertion/deletion polymorphism is another example of a locus showing different associations in different populations (50, 51). A recent meta-analysis of nine studies including 1333 RHD cases and 1212 controls from seven different populations showed no correlation with the disease or disease subtypes (52).
A few more genes have been linked to RHD in different studies, but there are currently no meta-analyses to confirm or reject these findings. Macrophage migration inhibitory factor (MIF) promoter polymorphisms -173C and -794 (5-8 CATT repeats) have been associated with RHD age of onset in a Saudi Arabian population (53). The C allele at−173 was associated with higher MIF expression in T cell lines (54). In vitro and in vivo studies found that a greater number of CATT repeats is associated with higher gene expression (55, 56). Interestingly, these promoter polymorphisms had a dual impact on the development of RHD in the studied population. The -173C allele was associated with reduction in RHD risk and later disease onset, the -794 CATT6 allele was associated with increased risk, and the -794 CATT5 allele was associated with reduced risk. The dual influence of these alleles on other autoimmune diseases has previously been reported. In systemic lupus erythematosus (SLE), higher expression MIF alleles (173C and 794 extended alleles) were associated with a lower disease risk (57). However, lower expression alleles were associated with reduced end-stage organ involvement (58). Recently, MIF was reported to have a regulatory effect on IL-1β production via the NLRP3 inflammasome (59). Mononuclear cells derived from ARF patients peripheral blood were shown to have persistent production of IL-1β after GAS infection, which was suppressed after hydroxychloroquine treatment (2). These new findings open up new avenues for research on the role of inflammasomes in RHD and perhaps targeting the pathway therapeutically (2).
IL17 is another immune-potent gene recently found to be associated with RHD. IL-17, also called IL-17A, is released by a subset of T helper cells known as Th17 cells. The IL17 promoter polymorphism (rs2275913), which has been linked to other autoimmune diseases, was associated with RHD in patients from India (OR 2.76; p = 0.021). Interestingly, this association was linked to mitral valve lesions (OR 2.74; p = 0.039) (60). The IL17 rs2275913-A allele was associated with higher expression of IL-17 (61). Wang et al. showed that IL-17 is needed to protect the host from GAS (62). On the other hand, prolonged and excessive Th17 exposure can cause valvular damage by enhancing neutrophil infiltration (63). IL-17 activates endothelial cells via the p38 MAPK pathway, which in turn induces the expression of several molecules such as VCAM and ICAM (64), which have important roles in the migration of different immune cells to the heart, particularly to the valves (65). However, the role of tissue-specific Th17 cells in heart and valvular tissue damage or repair has yet to be defined (66).
Another important pathway implicated in RHD is the lectin pathway, which includes mannose binding lectin-2 (MBL2), ficolin-1 (FCN1), and ficolin-2 (FCN2) [for a comprehensive review, see Beltrame et al. (42)]. FCN2 promoter polymorphisms (−986/−602/−4 GGA) have been linked to low protein levels and an increased risk of RHD (OR 1.6, p = 0.02), and the AGA haplotype was protective of the disease (67). FNC1 promoter polymorphisms (−1981A, 542A, −144A, and -33T) were associated with increased gene expression and were protective; however, two alleles (−1981A and −144A) increased the risk of developing valvular stenosis and mitral insufficiency (OR 0.75, p = 0.009 and OR 3.37, p = 0.027, respectively) (68). The authors proposed that alleles associated with increased production of FCN1 may help to eliminate the infection; however, prolonged inflammation under the influence of these alleles may cause tissue damage. This dual role of FCN1, similar to the dual role of MIF, highlights the complexity of the gene-environment interaction in RHD after GAS infection.
Approaches for Determining RHD Disease Variants
Large scale omics studies including genomic, transcriptomic, epigenomic, and metabolomic approaches combined with advanced computational approaches are now providing information about the mechanisms of many infectious, autoimmune, and rheumatic diseases (69). Unfortunately, very few such studies have been performed in RHD. Nevertheless, the last two decades have seen a huge amount of “big” omics data generated in many other rheumatic autoimmune diseases (70) sharing many features with RHD such as symptoms and progression. In addition, rheumatic autoimmune diseases share many disease genetic variants, a phenomenon called pleiotropy. Therefore, the results of rheumatic autoimmune disease studies can be utilized for further investigation in RHD and may lead to potential drug repositioning (71, 72).
GWAS have successfully identified important genomic regions for rheumatic autoimmune diseases, and advances in array design and statistical approaches have significantly increased the power of GWAS. For example, several studies have used fine mapping to identify likely causative variants from GWAS using a combination of genotyping arrays with specific sets of SNPs and statistical approaches to define a smaller group of SNPs that are statistically causative. This has allowed researchers to focus on this group of SNPs, known as credible sets, for further functional studies. In a recent study, Immunochip arrays were used to identify rheumatoid arthritis (RA) risk markers in African American patients (73). However, GWAS mapping precision, that is the physical distances between the top associated variant and the causative variants, is limited, as each locus identified by GWAS is co-inherited (in LD) with other genes, regulatory elements, or non-coding transcripts. Fine mapping by imputation following genotyping arrays, where all variants are statistically inferred from the reference genome, provides a cost-effective strategy to identify many disease alleles. Wu et al. (74) showed that imputation strategies are comparable to whole genome sequencing (WGS)-based approaches for common variants but not for rare variants. This was mainly due to that fact that rare variants are not adequately presented in the reference panels such as the 1000 Genome Project, HapMap, and 1000 UK Genomes. Moreover, samples representing populations from developing countries are not sufficiently represented in these panel sets, so rare causal variants for diseases common in these countries, such as RHD, may not be present. WGS-based approaches can overcome the imputation approach limitations. In a multinational study, van Rheenen et al. identified a novel rare functional variant associated with amyotrophic lateral sclerosis using WGS-based GWAS (75). Although this approach is the current method of choice for association studies in large research centers, it is rarely applied in developing countries due to its prohibitive cost. Furthermore, it has been recommended that the usual GWAS statistical significance threshold of an observed p-value < 5 x 10−8 should be increased to 5 x 10−11 (74), which requires large numbers of participants.
Whole exome sequencing (WES) has been increasingly applied to study rare variants in different autoimmune diseases due to its low sequencing cost. For example, by dividing sarcoidosis patients into groups based on HLA markers, rare variants were detected by WES in a limited sample size, which were associated with prognosis (76). HLA subgrouping helps to reduce genetic heterogeneity. Moreover, applying WES on individuals with extreme phenotypes increases the study power by increasing the number of relevant functional variants in the study, thereby increasing the probability of detecting new associations (77). This successful identification of disease-causing variants using WES is due to the fact that the low frequency variants are enriched in exonic regions (78). However, variations in regulatory elements influence gene expression in many immune cells in autoimmune diseases (79), and the majority of confirmed disease-associated loci in autoimmunity are non-coding variants (80). It has been estimated that up to 80% of the genome is responsible for gene regulation (81). Traditional WES approaches do not cover these elements, but extended exome sequencing has been used to identify these elements (82). This approach helps to reduce costs and covers all possible important regions in autoimmune conditions. However, to really understand how these regulatory elements contribute to autoimmune diseases, we need to understand how these elements regulate gene function at the cellular level (83).
Once a regulatory element has been associated with a complex disease, follow-up experiments are often conducted to understand the mechanism by which the variants influence phenotype. As an initial step, researchers carry out gene expression and DNA methylation and histone modification studies in various cell types. Genetic variants that influence gene expression level are called expression quantitative trait loci (eQTLs). eQTLs are classified according to their location (local or distant) or mode of action (cis or trans) [for a review, see (84)]. eQTLs are influenced by the environment, and some eQTLs behave differently under different conditions (85). Therefore, an important aspect of eQTL studies is to correlate the genetic locus with the correct cell type in the correct activation state in both health and disease. Although cell separation techniques may affect cellular status, they have provided relevant information regarding gene expression (known as transcript abundance) in many cases (86, 87). An interesting new approach is single-cell eQTL and the establishment of the Single-cell eQTLGen Consortium (sc-eQTLGen), which aims to map the upstream and downstream interactions of disease-related eQTLs in individual immune cell types (88). Unfortunately, in RHD, the immune cells most relevant to its pathogenesis remain unknown. In a recent study, Kim et al. reported disruption in the IL-1β granulocyte-macrophage colony-stimulating factor (GM-CSF) cytokine pathway in peripheral blood mononuclear cells from rheumatic fever patients stimulated with heat-inactivated GAS strains and studied gene expression in these cells (89). They identified hydroxychloroquine as an adjunctive therapy that might limit the dysregulation of this pathway. Further work to identify disease-causing pathways using specific cell populations might provide accurate insights into disease pathogenesis and correlations between gene expression and genetic variants.
Final Remarks and Conclusions
In the COVID-19 era, it is becoming even more important to identify pleiotropic immune mechanisms in rheumatic diseases since there may be shared mechanisms and pathways that are exploitable from the therapeutic perspective (90). The environmental trigger for autoimmunity in RHD is known; however, the inflammatory processes that occur after the infection and lead to RHD are still not well-understood. Understanding these processes at the cellular and molecular levels is crucial to deciphering RHD pathogenesis. In the last 3 years, while we have seen some advances in RHD genetics, GWAS studies have provided inconsistent results. One explanation is that these studies were mainly carried out in indigenous populations and trans-ethnic differences influence disease mechanisms. Another important factor is the widespread diversity of GAS strains, which may evoke different immune signaling pathways. Establishing which cytokine and cell signaling pathways are dysregulated in specific cell types in a specific ethnic population and the identification of the genetic variants that cause these processes are required to gain further insights into the immunological mechanisms of RHD. Exome sequencing and fine mapping will produce robust results. However, the lack of reference control data from different developing countries will challenge the identification of rare variants. Population genome project initiatives in developing countries and the Middle East, such as Qatar Genome (91), are expected to provide good reference data to help with these efforts.
Future studies should focus on understanding the function of any discovered variant through transcriptomic, epigenomics, proteomics, and metabolomics studies. Transcriptomic studies will define new disease pathways and pathogenic cell types, while epigenomic studies will link genetic promoters and enhancers with target genes. In rheumatic diseases, different cells play different roles; therefore, it is important to identify pathologically important immune cells. Single-cell technologies combined with advanced computational approaches and/or machine learning developed over the last decade have enabled integrative analysis of single-cell multi-omics data in different immunological diseases (92, 93). Future RHD studies will require collaboration between wet lab scientists, clinicians, and computational biologists to validate associations. Moreover, sequencing individuals with extreme phenotypes such as early disease onset will help to detect the strongest and most precise genetic associations. An important lesson emerging from other rheumatic diseases is that genetic associations influence disease phenotype and not only generic disease susceptibility. Therefore, detailed and accurate clinical phenotyping would also be helpful to identify new genetic factors that may influence disease progression and the underlying mechanisms. Pursuing this course will open up new avenues for early diagnosis, new therapeutic targets, and vaccine development.
Author Contributions
AA wrote the first draft of the manuscript and drafted the figure. MA-M edited and revised the manuscript and contributed figures. All authors revised the manuscript critically and approved the final version of the manuscript.
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.
References
1. Carapetis JR, Steer AC, Mulholland EK, Weber M. The global burden of group A streptococcal diseases. Lancet Infect Dis. (2005) 5:685–94. doi: 10.1016/S1473-3099(05)70267-X
2. Kuriakose T, Kanneganti T-D. Is Inflammasome a potential target of prophylaxis in rheumatic heart disease? Circulation. (2018) 138:2662–5. doi: 10.1161/CIRCULATIONAHA.118.036771
3. Marijon E, Mirabel M, Celermajer DS, Jouven X. Rheumatic heart disease. Lancet. (2012) 379:953–64. doi: 10.1016/S0140-6736(11)61171-9
4. Watkins DA, Johnson CO, Colquhoun SM, Karthikeyan G, Beaton A, Bukhman G, et al. Global, Regional, and National Burden of Rheumatic Heart Disease, 1990-2015. N Engl J Med. (2017) 377:713–22. doi: 10.1056/NEJMoa1603693
5. Pastore S, De Cunto A, Benettoni A, Berton E, Taddio A, Lepore L. The resurgence of rheumatic fever in a developed country area: the role of echocardiography. Rheumatology. (2011) 50:396–400. doi: 10.1093/rheumatology/keq290
6. Macleod CK, Bright P, Steer AC, Kim J, Mabey D, Parks T. Neglecting the neglected: the objective evidence of underfunding in rheumatic heart disease. Trans R Soc Trop Med Hyg. (2019) 113:287–90. doi: 10.1093/trstmh/trz014
7. Szczygielska I, Hernik E, Kołodziejczyk B, Gazda A, Maślińska M, Gietka P. Rheumatic fever – new diagnostic criteria. Rheumatology. (2018) 56:37–41. doi: 10.5114/reum.2018.74748
8. Reményi B, Wilson N, Steer A, Ferreira B, Kado J, Kumar K, et al. World Heart Federation criteria for echocardiographic diagnosis of rheumatic heart disease–an evidence-based guideline. Nat Rev Cardiol. (2012) 9:297–309. doi: 10.1038/nrcardio.2012.7
9. Breda L, Miulli E, Marzetti V, Chiarelli F, Marcovecchio ML. Rheumatic fever: a disease still to be kept in mind. Rheumatology. (2013) 52:953. doi: 10.1093/rheumatology/ket016
10. Sika-Paotonu D, Beaton A, Raghu A, Steer A, Carapetis J. Acute rheumatic fever rheumatic heart disease. In: Ferretti JJ, Stevens DL, Fischetti VA, editors. Streptococcus pyogenes : Basic Biology to Clinical Manifestations. Oklahoma City, OK: University of Oklahoma Health Sciences Center. Available online at: http://www.ncbi.nlm.nih.gov/books/NBK425394/ (accessed November 9, 2019)
11. Burke RJ, Chang C. Diagnostic criteria of acute rheumatic fever. Autoimmun Rev. (2014) 13:503–7. doi: 10.1016/j.autrev.2014.01.036
12. Punukollu M, Mushet N, Linney M, Hennessy C, Morton M. Neuropsychiatric manifestations of Sydenham's chorea: a systematic review. Dev Med Child Neurol. (2016) 58:16–28. doi: 10.1111/dmcn.12786
13. Guilherme L, Kalil J, Cunningham M. Molecular mimicry in the autoimmune pathogenesis of rheumatic heart disease. Autoimmunity. (2006) 39:31–9. doi: 10.1080/08916930500484674
14. Cunningham MW. Rheumatic fever, autoimmunity and molecular mimicry: the streptococcal connection. Int Rev Immunol. (2014) 33:314–29. doi: 10.3109/08830185.2014.917411
15. Bright PD, Mayosi BM, Martin WJ. An immunological perspective on rheumatic heart disease pathogenesis: more questions than answers. Heart. (2016) 102:1527–32. doi: 10.1136/heartjnl-2015-309188
16. Engel ME, Stander R, Vogel J, Adeyemo AA, Mayosi BM. Genetic susceptibility to acute rheumatic fever: a systematic review and meta-analysis of twin studies. PLoS ONE. (2011) 6:e25326. doi: 10.1371/journal.pone.0025326
17. Selmi C, Lu Q, Humble MC. Heritability versus the role of the environment in autoimmunity. J Autoimmun. (2012) 39:249–2. doi: 10.1016/j.jaut.2012.07.011
18. Martin WJ, Steer AC, Smeesters PR, Keeble J, Inouye M, Carapetis J, et al. Post-infectious group A streptococcal autoimmune syndromes and the heart. Autoimmun Rev. (2015) 14:710–25. doi: 10.1016/j.autrev.2015.04.005
19. Carapetis JR, Currie BJ, Mathews JD. Cumulative incidence of rheumatic fever in an endemic region: a guide to the susceptibility of the population? Epidemiol Infect. (2000) 124:239–44. doi: 10.1017/s0950268800003514
20. Bryant PA, Robins-Browne R, Carapetis JR, Curtis N. Some of the people, some of the time: susceptibility to acute rheumatic fever. Circulation. (2009) 119:742–53. doi: 10.1161/CIRCULATIONAHA.108.792135
21. Gray L-A, D'Antoine HA, Tong SYC, McKinnon M, Bessarab D, Brown N, et al. Genome-wide analysis of genetic risk factors for rheumatic heart disease in Aboriginal Australians provides support for pathogenic molecular mimicry. J Infect Dis. (2017) 216:1460–70. doi: 10.1093/infdis/jix497
22. Hallioglu O, Mesci L, Ozer S. DRB1, DQA1, DQB1 genes in Turkish children with rheumatic fever. Clin Exp Rheumatol. (2005) 23:117–20.
23. Kudat H, Telci G, Sozen AB, Oguz F, Akkaya V, Ozcan M, et al. The role of HLA molecules in susceptibility to chronic rheumatic heart disease. Int J Immunogenet. (2006) 33:41–4. doi: 10.1111/j.1744-313X.2006.00562.x
24. Haydardedeoglu FE, Tutkak H, Köse K, Düzgün N. Genetic susceptibility to rheumatic heart disease and streptococcal pharyngitis: association with HLA-DR alleles. Tissue Antigens. (2006) 68:293–6. doi: 10.1111/j.1399-0039.2006.00678.x
25. Gündogdu F, Islamoglu Y, Pirim I, Gurlertop Y, Dogan H, Arslan S, et al. Human leukocyte antigen (HLA) class I and II alleles in Turkish patients with rheumatic heart disease. J Heart Valve Dis. (2007) 16:293–9.
26. Karakurt C, Celiloglu C, Ozgen U, Yesilada E, Yologlu S, Zabriskie J, et al. Presence of a D8/17 B lymphocyte marker and HLA-DR subgroups in patients with rheumatic heart disease. Anadolu Kardiyol Derg. (2011) 11:314–8. doi: 10.5152/akd.2011.082
27. Guédez Y, Kotby A, El-Demellawy M, Galal A, Thomson G, Zaher S, et al. HLA class II associations with rheumatic heart disease are more evident and consistent among clinically homogeneous patients. Circulation. (1999) 99:2784–2790. doi: 10.1161/01.cir.99.21.2784
28. El-Hagrassy N, El-Chennawi F, Zaki ME-S, Fawzy H, Zaki A, Joseph N. HLA class I and class II HLA DRB profiles in Egyptian children with rheumatic valvular disease. Pediatr Cardiol. (2010) 31:650–6. doi: 10.1007/s00246-010-9663-3
29. Stanevicha V, Eglite J, Zavadska D, Sochnevs A, Shantere R, Gardovska D. HLA class II DR and DQ genotypes and haplotypes associated with rheumatic fever among a clinically homogeneous patient population of Latvian children. Arthritis Res Ther. (2007) 9:R58. doi: 10.1186/ar2216
30. Koyanagi T, Koga Y, Nishi H, Toshima H, Sasazuki T, Imaizumi T, et al. DNA typing of HLA class II genes in Japanese patients with rheumatic heart disease. J Mol Cell Cardiol. (1996) 28:1349–53. doi: 10.1006/jmcc.1996.0125
31. Rehman S, Akhtar N, Ahmad W, Ayub Q, Mehdi SQ, Mohyuddin A. Human leukocyte antigen (HLA) class II association with rheumatic heart disease in Pakistan. J Heart Valve Dis. (2007) 16:300–4.
32. Okello E, Beaton A, Mondo CK, Kruszka P, Kiwanuka N, Odoi-Adome R, et al. Rheumatic heart disease in Uganda: the association between MHC class II HLA DR alleles and disease: a case control study. BMC Cardiovasc Disord. (2014) 14:28. doi: 10.1186/1471-2261-14-28
33. Sreekanth MS, Esdan Basha SK, Arun Kumar G, Govindaraju S, Pradeep Nayar N, Pitchappan RM. Association of IL-1β +3953 C and HLA-DRB1*15 with coronary artery and rheumatic heart diseases in South India. Hum Immunol. (2016) 77:1275–9. doi: 10.1016/j.humimm.2016.08.003
34. Hernandezpacheco G. MHC class II alleles in Mexican patients with rheumatic heart disease. Int J Cardiol. (2003) 92:49–54. doi: 10.1016/S0167-5273(03)00040-8
35. Muhamed B, Parks T, Sliwa K. Genetics of rheumatic fever and rheumatic heart disease. Nat Rev Cardiol. (2020) 17:145–54. doi: 10.1038/s41569-019-0258-2
36. Cotsapas C, Voight BF, Rossin E, Lage K, Neale BM, Wallace C, et al. Pervasive sharing of genetic effects in autoimmune disease. PLoS Genet. (2011) 7:e1002254. doi: 10.1371/journal.pgen.1002254
37. Parkes M, Cortes A, van Heel DA, Brown MA. Genetic insights into common pathways and complex relationships among immune-mediated diseases. Nat Rev Genet. (2013) 14:661–73. doi: 10.1038/nrg3502
38. Raj T, Kuchroo M, Replogle JM, Raychaudhuri S, Stranger BE, De Jager PL. Common risk alleles for inflammatory diseases are targets of recent positive selection. Am J Hum Genet. (2013) 92:517–29. doi: 10.1016/j.ajhg.2013.03.001
39. Parks T, Mirabel MM, Kado J, Auckland K, Nowak J, Rautanen A, et al. Association between a common immunoglobulin heavy chain allele and rheumatic heart disease risk in Oceania. Nat Commun. (2017) 8:14946. doi: 10.1038/ncomms14946
40. Auckland K, Mittal B, Cairns BJ, Garg N, Kumar S, Mentzer AJ, et al. The human leukocyte antigen locus and rheumatic heart disease susceptibility in South Asians and Europeans. Sci Rep. (2020) 10:9004. doi: 10.1038/s41598-020-65855-8
41. Gutierrez-Arcelus M, Rich SS, Raychaudhuri S. Autoimmune diseases - connecting risk alleles with molecular traits of the immune system. Nat Rev Genet. (2016) 17:160–74. doi: 10.1038/nrg.2015.33
42. Beltrame MH, Catarino SJ, Goeldner I, Boldt ABW, de Messias-Reason IJ. The lectin pathway of complement and rheumatic heart disease. Front Pediatr. (2014) 2:148. doi: 10.3389/fped.2014.00148
43. Guilherme L, Köhler KF, Postol E, Kalil J. Genes, autoimmunity and pathogenesis of rheumatic heart disease. Ann Pediatr Cardiol. (2011) 4:13–21. doi: 10.4103/0974-2069.79617
44. Terao C, Yoshifuji H, Matsumura T, Naruse TK, Ishii T, Nakaoka Y, et al. Genetic determinants and an epistasis of LILRA3 and HLA-B*52 in Takayasu arteritis. Proc Natl Acad Sci USA. (2018) 115:13045–50. doi: 10.1073/pnas.1808850115
45. Osgood JA, Knight JC. Translating GWAS in rheumatic disease: approaches to establishing mechanism and function for genetic associations with ankylosing spondylitis. Brief Funct Genomics. (2018) 17:308–18. doi: 10.1093/bfgp/ely015
46. Muhamed B, Shaboodien G, Engel ME. Genetic variants in rheumatic fever and rheumatic heart disease. Am J Med Genet C Semin Med Genet. (2020) 184:159–77. doi: 10.1002/ajmg.c.31773
47. Zheng R-L, Zhang H, Jiang W-L. Tumor necrosis factor-alpha 308G>A polymorphism and risk of rheumatic heart disease: a meta-analysis. Sci Rep. (2014) 4:4731. doi: 10.1038/srep04731
48. Abdallah AM, Alnuzha A, Al-Mazroea AH, Eldardear AE, AlSamman AY, Almohammadi Y, et al. IL10 promoter polymorphisms are associated with rheumatic heart disease in Saudi Arabian patients. Pediatr Cardiol. (2016) 37:99–105. doi: 10.1007/s00246-015-1245-y
49. Dai W, Ye Z, Lu H, Su Q, Li H, Li L. Meta-analysis of the relationship between single nucleotide polymorphism of IL-10-1082G/A and rheumatic heart disease. Oncotarget. (2018) 9:12343–50. doi: 10.18632/oncotarget.23901
50. Al-Harbi KM, Almuzaini IS, Morsy MM, Abdelaziz NA, Al-Balawi AM, Abdallah AM. Angiotensin-converting enzyme gene insertion/deletion polymorphism in Saudi patients with rheumatic heart disease. Saudi Med J. (2015) 36:176–80. doi: 10.15537/smj.2015.2.10267
51. Gupta U, Mishra A, Rathore SS, Agarwal SK, Pande S, Garg N, et al. Association of angiotensin I-converting enzyme gene insertion/deletion polymorphism with rheumatic heart disease in Indian population and meta-analysis. Mol Cell Biochem. (2013) 382:75–82. doi: 10.1007/s11010-013-1719-2
52. Tian Y, Ge Z, Xing Y, Sun Y, Ying J. Correlation of angiotensin I-converting enzyme gene insertion/deletion polymorphism with rheumatic heart disease: a meta-analysis. Biosci Rep. (2016) 36:e00412. doi: 10.1042/BSR20160151
53. Abdallah AM, Al-Mazroea AH, Al-Harbi WN, Al-Harbi NA, Eldardear AE, Almohammadi Y, et al. Impact of MIF gene promoter variations on risk of rheumatic heart disease and its age of onset in Saudi Arabian patients. Front Immunol. (2016) 7:98. doi: 10.3389/fimmu.2016.00098
54. Donn R, Alourfi Z, De Benedetti F, Meazza C, Zeggini E, Lunt M, et al. Mutation screening of the macrophage migration inhibitory factor gene: positive association of a functional polymorphism of macrophage migration inhibitory factor with juvenile idiopathic arthritis. Arthritis Rheum. (2002) 46:2402–9. doi: 10.1002/art.10492
55. Adamali H, Armstrong ME, McLaughlin AM, Cooke G, McKone E, Costello CM, et al. Macrophage migration inhibitory factor enzymatic activity, lung inflammation, and cystic fibrosis. Am J Respir Crit Care Med. (2012) 186:162–9. doi: 10.1164/rccm.201110-1864OC
56. Kang I, Bucala R. The immunobiology of MIF: function, genetics and prospects for precision medicine. Nat Rev Rheumatol. (2019) 15:427–37. doi: 10.1038/s41584-019-0238-2
57. Lang T, Foote A, Lee JPW, Morand EF, Harris J. MIF: implications in the pathoetiology of systemic lupus erythematosus. Front Immunol. (2015) 6:577. doi: 10.3389/fimmu.2015.00577
58. Sreih A, Ezzeddine R, Leng L, LaChance A, Yu G, Mizue Y, et al. Dual effect of the macrophage migration inhibitory factor gene on the development and severity of human systemic lupus erythematosus. Arthritis Rheum. (2011) 63:3942–3951. doi: 10.1002/art.30624
59. Shin MS, Kang Y, Wahl ER, Park H-J, Lazova R, Leng L, et al. Macrophage migration inhibitory factor regulates U1 small nuclear RNP immune complex-mediated activation of the NLRP3 inflammasome. Arthritis Rheumatol. (2019) 71:109–20. doi: 10.1002/art.40672
60. Poomarimuthu M, Elango S, Solomon PR, Soundrapandian S, Mariakuttikan J. Association of IL17 and IL23R gene polymorphisms with rheumatic heart disease in South Indian population. Immunol Invest. (2018) 47:754–64. doi: 10.1080/08820139.2018.1493053
61. Espinoza JL, Takami A, Nakata K, Onizuka M, Kawase T, Akiyama H, et al. A genetic variant in the IL-17 promoter is functionally associated with acute graft-versus-host disease after unrelated bone marrow transplantation. PLoS ONE. (2011) 6:e26229. doi: 10.1371/journal.pone.0026229
62. Wang B, Dileepan T, Briscoe S, Hyland KA, Kang J, Khoruts A, et al. Induction of TGF-beta1 and TGF-beta1-dependent predominant Th17 differentiation by group A streptococcal infection. Proc Natl Acad Sci USA. (2010) 107:5937–42. doi: 10.1073/pnas.0904831107
63. Yeh C-Y, Shun C-T, Kuo Y-M, Jung C-J, Hsieh S-C, Chiu Y-L, et al. Activated human valvular interstitial cells sustain interleukin-17 production to recruit neutrophils in infective endocarditis. Infect Immun. (2015) 83:2202–12. doi: 10.1128/IAI.02965-14
64. Mai J, Nanayakkara G, Lopez-Pastrana J, Li X, Li Y-F, Wang X, et al. Interleukin-17A promotes aortic endothelial cell activation via transcriptionally and post-translationally activating p38 mitogen-activated protein kinase (MAPK) pathway. J Biol Chem. (2016) 291:4939–54. doi: 10.1074/jbc.M115.690081
65. Guilherme L, Kalil J. Rheumatic heart disease: molecules involved in valve tissue inflammation leading to the autoimmune process and Anti-S. pyogenes vaccine. Front Immunol. (2013) 4:352. doi: 10.3389/fimmu.2013.00352
66. Dominguez-Villar M, Hafler DA. Regulatory T cells in autoimmune disease. Nat Immunol. (2018) 19:665–673. doi: 10.1038/s41590-018-0120-4
67. Messias-Reason IJ, Schafranski MD, Kremsner PG, Kun JFJ. Ficolin 2 (FCN2) functional polymorphisms and the risk of rheumatic fever and rheumatic heart disease: FCN2 polymorphisms in autoimmune disease. Clin Exp Immunol. (2009) 157:395–9. doi: 10.1111/j.1365-2249.2009.03975.x
68. Catarino SJ, Andrade FA, Boldt ABW, Guilherme L, Messias-Reason IJ. Sickening or healing the heart? The association of ficolin-1 and rheumatic fever. Front Immunol. (2018) 9:3009. doi: 10.3389/fimmu.2018.03009
69. Zhernakova A, Withoff S, Wijmenga C. Clinical implications of shared genetics and pathogenesis in autoimmune diseases. Nat Rev Endocrinol. (2013) 9:646–59. doi: 10.1038/nrendo.2013.161
70. Matucci-Cerinic M, Seferović PM. Heart involvement in autoimmune rheumatic diseases: the “phantom of the opera.” Rheumatology. (2006) 45:iv1–3. doi: 10.1093/rheumatology/kel316
71. Asimit JL, Rainbow DB, Fortune MD, Grinberg NF, Wicker LS, Wallace C. Stochastic search and joint fine-mapping increases accuracy and identifies previously unreported associations in immune-mediated diseases. Nat Commun. (2019) 10:3216. doi: 10.1038/s41467-019-11271-0
72. Jia X, Shi N, Feng Y, Li Y, Tan J, Xu F, et al. Identification of 67 pleiotropic genes associated with seven autoimmune/autoinflammatory diseases using multivariate statistical analysis. Front Immunol. (2020) 11:30. doi: 10.3389/fimmu.2020.00030
73. Danila MI, Laufer VA, Reynolds RJ, Yan Q, Liu N, Gregersen PK, et al. Dense genotyping of immune-related regions identifies loci for rheumatoid arthritis risk and damage in African Americans. Mol Med. (2017) 23:177–87. doi: 10.2119/molmed.2017.00081
74. Wu Y, Zheng Z, Visscher PM, Yang J. Quantifying the mapping precision of genome-wide association studies using whole-genome sequencing data. Genome Biol. (2017) 18:86. doi: 10.1186/s13059-017-1216-0
75. PARALS Registry SLALOM Group SLAP Registry FALS Sequencing Consortium SLAGEN Consortium NNIPPS Study Group, et al. Genome-wide association analyses identify new risk variants and the genetic architecture of amyotrophic lateral sclerosis. Nat Genet. (2016) 48:1043–8. doi: 10.1038/ng.3622
76. Lahtela E, Kankainen M, Sinisalo J, Selroos O, Lokki M-L. Exome sequencing identifies susceptibility loci for sarcoidosis prognosis. Front Immunol. (2019) 10:2964. doi: 10.3389/fimmu.2019.02964
77. Peloso GM, Rader DJ, Gabriel S, Kathiresan S, Daly MJ, Neale BM. Phenotypic extremes in rare variant study designs. Eur J Hum Genet. (2016) 24:924–930. doi: 10.1038/ejhg.2015.197
78. Tennessen JA, Bigham AW, O'Connor TD, Fu W, Kenny EE, Gravel S, et al. Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. (2012) 337:64–9. doi: 10.1126/science.1219240
79. NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium, Gutierrez-Arcelus M, Baglaenko Y, Arora J, Hannes S, Luo Y, et al. Allele-specific expression changes dynamically during T cell activation in HLA and other autoimmune loci. Nat Genet. (2020) 52:247–53. doi: 10.1038/s41588-020-0579-4
80. Seldin MF. The genetics of human autoimmune disease: A perspective on progress in the field and future directions. J Autoimmun. (2015) 64:1–12. doi: 10.1016/j.jaut.2015.08.015
81. Elkon R, Agami R. Characterization of noncoding regulatory DNA in the human genome. Nat Biotechnol. (2017) 35:732–46. doi: 10.1038/nbt.3863
82. Eriksson D, Bianchi M, Landegren N, Nordin J, Dalin F, Mathioudaki A, et al. Extended exome sequencing identifies BACH2 as a novel major risk locus for Addison's disease. J Intern Med. (2016) 280:595–608. doi: 10.1111/joim.12569
83. Vockley CM, Barrera A, Reddy TE. Decoding the role of regulatory element polymorphisms in complex disease. Curr Opin Genet Dev. (2017) 43:38–45. doi: 10.1016/j.gde.2016.10.007
84. Albert FW, Kruglyak L. The role of regulatory variation in complex traits and disease. Nat Rev Genet. (2015) 16:197–212. doi: 10.1038/nrg3891
85. Idaghdour Y, Awadalla P. Exploiting gene expression variation to capture gene-environment interactions for disease. Front Genet. (2013) 3:228. doi: 10.3389/fgene.2012.00228
86. Inshaw JRJ, Cutler AJ, Burren OS, Stefana MI, Todd JA. Approaches and advances in the genetic causes of autoimmune disease and their implications. Nat Immunol. (2018) 19:674–84. doi: 10.1038/s41590-018-0129-8
87. Zeng B, Lloyd-Jones LR, Montgomery GW, Metspalu A, Esko T, Franke L, et al. Comprehensive multiple eQTL detection and its application to GWAS interpretation. Genetics. (2019) 212:905–18. doi: 10.1534/genetics.119.302091
88. van der Wijst M, de Vries D, Groot H, Trynka G, Hon C, Bonder M, et al. The single-cell eQTLGen consortium. eLife. (2020) 9:e52155. doi: 10.7554/eLife.52155
89. Kim ML, Martin WJ, Minigo G, Keeble JL, Garnham AL, Pacini G, et al. Dysregulated IL-1β-GM-CSF axis in acute rheumatic fever that is limited by hydroxychloroquine. Circulation. (2018) 138:2648–61. doi: 10.1161/CIRCULATIONAHA.118.033891
90. Ehrenfeld M, Tincani A, Andreoli L, Cattalini M, Greenbaum A, Kanduc D, et al. Covid-19 and autoimmunity. Autoimmun Rev. (2020) 19:102597. doi: 10.1016/j.autrev.2020.102597
91. Al Thani A, Fthenou E, Paparrodopoulos S, Al Marri A, Shi Z, Qafoud F, et al. Qatar biobank cohort study: study design and first results. Am J Epidemiol. (2019) 188:1420–33. doi: 10.1093/aje/kwz084
92. Cheung P, Khatri P, Utz PJ, Kuo AJ. Single-cell technologies — studying rheumatic diseases one cell at a time. Nat Rev Rheumatol. (2019) 15:340–54. doi: 10.1038/s41584-019-0220-z
Keywords: rheumatic heart, autoimmune diseases, group A streptococcus, exome sequencing, genetic association
Citation: Abdallah AM and Abu-Madi M (2021) The Genetic Control of the Rheumatic Heart: Closing the Genotype-Phenotype Gap. Front. Med. 8:611036. doi: 10.3389/fmed.2021.611036
Received: 28 September 2020; Accepted: 07 January 2021;
Published: 24 March 2021.
Edited by:
Konstantinos Thomas, University General Hospital Attikon, GreeceReviewed by:
George Lazaros, Hippokration General Hospital, GreeceChristos Koutsianas, National and Kapodistrian University of Athens, Greece
Copyright © 2021 Abdallah and Abu-Madi. 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: Marawan Abu-Madi, YWJ1bWFkaSYjeDAwMDQwO3F1LmVkdS5xYQ==