RESUMO
Background: Autoimmune diseases (AIDs) share a common molecular etiology and often present overlapping clinical presentations. Thus, this study aims to explore the complex molecular basis of AID by whole exome sequencing and computational biology analysis. Methods: Molecular screening of the consanguineous AID family and the computational biology characterization of the potential variants were performed. The potential variants were searched against the exome data of 100 healthy individuals and 30 celiac disease patients. Result: A complex inheritance pattern of PAK2 (V43A), TAP2 (F468Y), and PLCL1 (V473I) genetic variants was observed in the three probands of the AID family. The PAK2 variant (V43A) is a novel one, but TAP2 (F468Y) and PLCL1 (V473I) variants are extremely rare in local Arab (SGHP and GME) and global (gnomAD) databases. All these variants were localized in functional domains, except for the PAK2 variant (V43A) and were predicted to alter the structural (secondary structure elements, folding, active site confirmation, stability, and solvent accessibility) and functional (gene expression) features. Therefore, it is reasonable to postulate that the dysregulation of PAK2, TAP2, and PLCL1 genes is likely to elicit autoimmune reactions by altering antigen processing and presentation, T cell receptor signaling, and immunodeficiency pathways. Conclusion: Our findings highlight the importance of exploring the alternate inheritance patterns in families presenting complex autoimmune diseases, where classical genetic models often fail to explain their molecular basis. These findings may have potential implications for developing personalized therapies for complex disease patients.