Next generation sequencing reveals novel homozygous frameshift in PUS7 and splice acceptor variants in AASS gene leading to intellectual disability,developmental delay,dysmorphic feature and microcephaly

Intellectual developmental disorder with abnormal behavior, microcephaly and short stature (IDDABS), (OMIM# 618342) is an autosomal recessive condition described as developmental delay, poor or absent speech, intellectual disability, short stature, mild to progressive microcephaly, delayed psychomotor development, hyperactivity, seizure, along with mild to swear aggressive behavior. Homozygous frameshift mutation in Pseudouridine Synthase 7, Putative; (PUS7) OMIM# 616,261 NM_019042.3 and splice acceptor variants in Alpha-Aminoadipic Semialdehyde Synthase; (AASS) OMIM# 605,113 NM_005763.3 was funded. Whole exome sequencing (WES) technique was used as tool to identify the molecular diagnostic test. Different bioinformatics analysis done for WES data and we identified two novel mutations one as frameshift mutation c.606_607delGA, p.Ser282CysfsTer9 in the PUS7 gene and splice acceptor variants c.1767–1 G > A in the AASS gene has been reported. The pattern of family segregation maintained the pathogenicity of this variation associated with abnormal behavior, intellectual developmental disorder, microcephaly along with short stature IDDABS. Further, the WES data was validated in the family having other affected individuals and healthy controls (n = 100) was done using Sanger sequencing. Finally, our results further explained the role of WES in the disease diagnosis and elucidated that the mutation in PUS7 and AASS genes may lead an important role for the development of IDDABS in Saudi family.     

RTTN Mutations Cause Primary Microcephaly and Primordial Dwarfism in Humans

Primary microcephaly is a developmental brain anomaly that results from defective proliferation of neuroprogenitors in the germinal periventricular zone. More than a dozen genes are known to be mutated in autosomal-recessive primary microcephaly in isolation or in association with a more generalized growth deficiency (microcephalic primordial dwarfism), but the genetic heterogeneity is probably more extensive. In a research protocol involving autozygome mapping and exome sequencing, we recruited a multiplex consanguineous family who is affected by severe microcephalic primordial dwarfism and tested negative on clinical exome sequencing. Two candidate autozygous intervals were identified, and the second round of exome sequencing revealed a single intronic variant therein (c.2885+8A>G [p.Ser963^∗] in RTTN exon 23). RT-PCR confirmed that this change creates a cryptic splice donor and thus causes retention of the intervening 7 bp of the intron and leads to premature truncation. On the basis of this finding, we reanalyzed the exome file of a second consanguineous family affected by a similar phenotype and identified another homozygous change in RTTN as the likely causal mutation. Combined linkage analysis of the two families confirmed that RTTN maps to the only significant linkage peak. Finally, through international collaboration, a Canadian multiplex family affected by microcephalic primordial dwarfism and biallelic mutation of RTTN was identified. Our results expand the phenotype of RTTN-related disorders, hitherto limited to polymicrogyria, to include microcephalic primordial dwarfism with a complex brain phenotype involving simplified gyration.     

A novel homozygous mutation in <Emphasis Type="Italic">SZT2</Emphasis> gene in Saudi family with developmental delay,macrocephaly and epilepsy

Epileptic encephalopathies are genetically heterogeneous disorders which leads to epilepsy and cause neurological disorders. Seizure threshold 2 (SZT2) gene located on chromosome 1p34.2 encodes protein mainly expressed predominantly in the parietal and frontal cortex and dorsal root ganglia in the brain. Previous studies in mice showed that mutation in this gene can confers low seizure threshold, enhance epileptogenesis and in human may leads to facial dysmorphism, intellectual disability, seizure and macrocephaly. Objective of this study was to find out novel gene or novel mutation related to the gene phenotype. We have identified a large consanguineous Saudi family segregating developmental delay, intellectual disability, epilepsy, high forehead and macrocephaly. Exome sequencing was performed in affected siblings of the family to study the novel mutation. Whole exome sequencing data analysis, confirmed by subsequent Sanger sequencing validation study. Our results showed a novel homozygous mutation (c.9368G>A) in a substitution of a conserved glycine residue into a glutamic acid in the exon 67 of SZT2 gene. The mutation was ruled out in 100 unrelated healthy controls. The missense variant has not yet been reported as pathogenic in literature or variant databases. In conclusion, the here detected homozygous SZT2 variant might be the causative mutation that further explain epilepsy and developmental delay in this Saudi family.     

Two Novel Mutations in LAMC2 Gene in Iranian Families Affected by Junctional Epidermolysis Bullosa

Background:Junctional epidermolysis bullosa (JEB) is an autosomal recessive skin disorder with defective adhesion of dermal- epidermal within the lamina lucida region of the basement membrane zone. The main characterization of JEB is blistering and fragile skin and mucous membrane. Laminins are noncollagenous part of basement membrane and classified as a family of extracellular matrix glycoprotein. Laminins contain three chains: Laminin α, Laminin β and Laminin γ. LAMC2 (laminin subunit gamma 2) gene encodes γ subunit of laminin and its mutation contributes to JEB. Here, we report a disease-causing nonsense mutation and a large deletion mutation in LAMC2 gene in two families affected by JEB.Methods:Whole exome sequencing (WES) was carried out on the mother of patient in family I and the patient himself in family II to detect the underlying mutations. Then, sanger sequencing was performed to confirm the identified mutations.Results:Next generation sequencing (NGS) data analysis of the first family showed a novel, nonsense mutation in LAMC2 gene (LAMC2: {"type":"entrez-nucleotide","attrs":{"text":"NM_005562","term_id":"1653960840"}}NM_005562: exon14:c.C2143T: p.R715X). The heterozygous state of the mutation was confirmed by sanger sequencing in the parents and unaffected brother. In Family II, NGS data had no coverage in the large area of LAMC2 gene. Thus, to confirm the possible deletion sanger sequencing was done and blasting of sequence showed the deleted region of 9.4 kb (exon10-17) in LAMC2 gene.Conclusion:In summary, current study reported a novel disease-causing premature termination codon (PTC) mutation in LAMC2 gene and a large deletion mutation in patients affected by JEB.Key Words: Junctional Epidermolysis Bullosa, LAMC2 gene, Novel mutation, Skin disorder     

Waardenburg syndrome type II in a Chinese pedigree caused by frameshift mutation in the SOX10 gene

Waardenburg syndrome (WS) is a congenital hereditary disease, attributed to the most common symptoms of sensorineural deafness and iris hypopigmentation. It is also known as the hearing-pigmentation deficient syndrome. Mutations on SOXl0 gene often lead to congenital deafness and has been shown to play an important role in the pathogenesis of WS. We investigated one family of five members, with four patients exhibiting the classic form of WS2, whose DNA samples were analyzed by the technique of Whole-exome sequencing (WES). From analysis of WES data, we found that both the mother and all three children in the family have a heterozygous mutation on the Sex Determining Region Y - Box 10 (SOX10) gene. The mutation was c.298_300delinsGG in exon 2 of SOX10 ({"type":"entrez-nucleotide","attrs":{"text":"NM_006941","term_id":"1732746331"}}NM_006941), which leads to a frameshift of nine nucleotides, hence the amino acids (p. S100Rfs*9) are altered and the protein translation may be terminated prematurely. Further flow cytometry confirmed significant down-regulation of SOX10 protein, which indicated the SOX10 gene mutation was responsible for the pathogenesis of WS2 patients. In addition, we speculated that some other mutated genes might be related to disease phenotype in this family, which might also participate in promoting the progression of WS2.     

Whole exome sequencing,in silico and functional studies confirm the association of the GJB2 mutation p.Cys169Tyr with deafness and suggest a role for the TMEM59 gene in the hearing process

The development of next generation sequencing techniques has facilitated the detection of mutations at an unprecedented rate. These efficient tools have been particularly beneficial for extremely heterogeneous disorders such as autosomal recessive non-syndromic hearing loss, the most common form of genetic deafness. GJB2 mutations are the most common cause of hereditary hearing loss. Amongst them the NM_004004.5: c.506G > A (p.Cys169Tyr) mutation has been associated with varying severity of hearing loss with unclear segregation patterns. In this study, we report a large consanguineous Emirati family with severe to profound hearing loss fully segregating the GJB2 missense mutation p.Cys169Tyr. Whole exome sequencing (WES), in silico, splicing and expression analyses ruled out the implication of any other variants and confirmed the implication of the p.Cys169Tyr mutation in this deafness family. We also show preliminary murine expression analysis that suggests a link between the TMEM59 gene and the hearing process. The present study improves our understanding of the molecular pathogenesis of hearing loss. It also emphasizes the significance of combining next generation sequencing approaches and segregation analyses especially in the diagnosis of disorders characterized by complex genetic heterogeneity.     

Whole exome sequencing of a Saudi family and systems biology analysis identifies CPED1 as a putative causative gene to Celiac Disease

Celiac disease (CD) is a gastrointestinal disorder whose genetic basis is not fully understood. Therefore, we studied a Saudi family with two CD affected siblings to discover the causal genetic defect. Through whole exome sequencing (WES), we identified that both siblings have inherited an extremely rare and deleterious CPED1 genetic variant (c.241 A > G; p.Thr81Ala) segregating as autosomal recessive mutation, suggesting its putative causal role in the CD. Saudi population specific minor allele frequency (MAF) analysis has confirmed its extremely rare prevalence in homozygous condition (MAF is 0.0004). The Sanger sequencing analysis confirmed the absence of this homozygous variant in 100 sporadic Saudi CD cases. Genotype-Tissue Expression (GTEx) data has revealed that CPED1 is abundantly expressed in gastrointestinal mucosa. By using a combination of systems biology approaches like protein 3D modeling, stability analysis and nucleotide sequence conservation analysis, we have further established that this variant is deleterious to the structural and functional aspects of CPED1 protein. To the best of our knowledge, this variant has not been previously reported in CD or any other gastrointestinal disease. The cell culture and animal model studies could provide further insight into the exact role of CPED1 p.Thr81Ala variant in the pathophysiology of CD. In conclusion, by using WES and systems biology analysis, present study for the first-time reports CPED1 as a potential causative gene for CD in a Saudi family with potential implications to both disease diagnosis and genetic counseling.     

Genetic analysis of hereditary multiple exostoses in Tunisian families: a novel frame-shift mutation in the <Emphasis Type="Italic">EXT1</Emphasis> gene

Hereditary multiple exostoses (HME) is an autosomal dominant orthopaedic disorder most frequently caused by mutations in the EXT1 gene. The aim of the present study is to determine the underlying molecular defect of HME in two multigenerational Tunisian families with 21 affected members and to examine the degree of intrafamilial variability. Linkage analysis was performed using three microsatellite markers encompassing the EXT1 locus and mutation screening was carried out by direct sequencing. In family 1, evidence for linkage to EXT1 was obtained on the basis of a maximum LOD score of 4.26 at θ = 0.00 with D8S1694 marker. Sequencing of the EXT1 revealed a heterozygous G > T transversion (c.1019G>T) in exon 2, leading to a missense mutation at the codon 340 (p.Arg340Leu). In family 2 we identified a novel heterozygous 1 bp deletion in the exon 1 (c.529_531delA) leading to a premature codon stop and truncated EXT1 protein expression (p.Lys177LysfsX15). This mutation was associated with the evidence of an intrafamilial clinical variability and considered to be a novel disease-causing mutation in the EXT1 gene. These findings provide additional support for the involvement of EXT1 gene in the HME disease.     

NOTCH3基因c.520T>G突变致CADASIL一家系报告及文献复习

目的:研究探讨一CADASIL家系的临床特征及基因突变情况。方法:收集同一家系中3例CADASIL患者的临床资料,并对3例患者及先证者之兄进行全外显子测序(Whole Exome Sequencing, WES)。结果:该家系中3例患者临床表现多样,女性患者均有头痛病史,先证者及先证者之姐中年起病,先证者临床表现缺乏特异性,主要表现为头昏,认知功能检查正常,心理评估示轻度焦虑抑郁状态。先证者之姐主要表现为假性球麻痹及锥体束受损,认知功能检查示重度痴呆。先证者之女自4岁起诊断为癫痫-失神发作,认知功能检查示轻度认知功能障碍。影像学显示该家系3例患者均有脑白质病变,且随着年龄增大呈进行性发展,WES显示3例患者均存在NOTCH3基因第4外显子区域杂合突变:c.520T>G,导致氨基酸改变p.Cys174Gly。结论:NOTCH3基因c.520T>G所致该家系的临床表现具有多样性,且该家系中下一代起病较早,临床表现可与父代具有较大异质性,影像学表现可在青少年时期出现,并呈现进行加重的趋势。WES显示该家系中NOTCH3基因突变为第4外显子的杂合突变,该位点突变致CADASIL为国内首次报道。     

Identification and characterization of the compound heterozygous variants of TYR gene in a northern Chinese family with Oculocutaneous albinism type 1

Oculocutaneous albinism (OCA) is a genetically heterogeneous disease and is most inherited in an autosomal recessive manner. The characteristic manifestation of OCA is due to disfunction of melanin synthesis. OCA1 is the most severe subtype of OCA and is caused by homozygous or compound heterozygous variants in tyrosinase (TYR) gene, which is the key gene for melanin synthesis. This study aimed to identify the genetic variants of a northern Chinese family with OCA1. Clinical information and peripheral blood samples were collected. PCR amplification and Sanger sequencing were used to detect the entire exons and adjacent flanking sequences of TYR gene. Functional prediction of variants was performed by various bioinformatic analyses, while the pathogenicity classification of variants was evaluated according to ACMG standards and guidelines. A missense variant NM_000372.5:c.107G > C;NP_000363.1:p.C36S was discovered in TYR gene which converted cysteine to serine. Another variant in intron, NM_000372.5:c.1037–7 T > A, also affected the function of TYR gene. We verified the pathogenicity of the intron variant with a pCAS2 mini-gene based splicing assay and found that c.1037–7 T > A led to an insertion of 5 bp upstream from the common acceptor site of exon 3, which caused a frameshift TYR:c.1037–7 T > A:p.G346Efs*11. The results showed that the compound heterozygous variants c.107G > C:p.C36S and c.1037–7 T > A:p.G346Efs*11 of TYR gene were the pathogenic variants for this OCA1 family.     

Analysis of RAS gene mutations in cytogenetically normal de novo acute myeloid leukemia patients reveals some novel alterations

Rat sarcoma gene (RAS) holds great importance in pathogenesis of acute myeloid leukemia (AML). The activated mutations in Neuroblastoma rat sarcoma viral oncogene homolog (NRAS) and Kirsten rat sarcoma viral oncogene homolog (KRAS) confers proliferative and survival signals, deliberating numerous effects on overall survival and progression free survival in AML patients. In this study thirty one (31) blood samples of adult newly diagnosed AML patients were collected to identify possible incidence of mutations through amplification of KRAS (exon 1 and 2) and NRAS gene (exon 1 and 2) using polymerase chain reaction (PCR). Amplicons were then subjected to sequencing and were analyzed through Geneious Prime 2019. Five of thirty one (16.12%) patients had altered sites in either NRAS or KRAS. The NRAS mutations were observed in three AML patients (N = 3, 9.67%). A novel missense mutation NRAS-I36R (239 T > G) representing a substitution of single nucleotide basepair found in NRAS exon 1 while exon 2 was detected with heterozygous mutation NRAS-E63X (318G > T) and insertion (A), resulting in frameshift of the amino acid sequence and insertion of two nucleotide basepairs (TA) in two of the patients. KRAS mutations (N = 2, 6.45%) were found in exon 1 whereas no mutations in KRAS exon 2 were detected in our patient cohort. Mutation in KRAS Exon 1, KRAS-D30N (280G > A) was observed in two patients and one of them also had a novel heterozygous mutation KRAS-L16N (240G > C). In addition there was no statistically significant association of mutRAS gene of AML patients with several prognostic markers including age, gender, karyotyping, CD34 positivity, cytogenetic abnormalities, total leukocyte count, white blood cell count and French-American-British (FAB) classification. However, the presence of mutRAS gene were strongly associated (p = 0.001) with increased percentage of bone marrow blasts. The prevalence of mutations in correlation with clinical and hematological parameter is useful for risk stratification in AML patients.     

Clinical whole exome sequencing revealed de novo heterozygous stop-gain and missense variants in the STXBP1 gene associated with epilepsy in Saudi families

Intellectual disability and developmental encephalopathies are mostly linked with infant epilepsy. Epileptic encephalopathy is a term that is used to define association between developmental delay and epilepsy. Mutations in the STXBP1 (Syntaxin-binding protein 1) gene have been previously reported in association with multiple severe early epileptic encephalopathies along with many neurodevelopmental disorders. Among the disorders produced due to any mutations in the STXBP1 gene is developmental and epileptic encephalopathy 4 (OMIM: 612164), is an autosomal dominant neurologic disorder categorized by the onset of tonic seizures in early infancy (usually in the first months of life). In this article, we report two Saudi families one with de novo heterozygous stop-gain mutation c.364C > T and a novel missense c. 305C > A p.Ala102Glu in exon 5 of the STXBP1 gene (OMIM: 602926) lead to development of epileptic encephalopathy 4. The variants identified in the current study broadened the genetic spectrum of STXBP1 gene related with diseases, which will help to add in the literature and benefit to the studies addressing this disease in the future.     

The unclassified variant: c.2044AD>G,p.T682A (het.) in exon 12 of the GLI3 gene in a patient with oral-facial-digital syndrome type II (Mohr syndrome) phenotype

The gene mutation for oral-facial digital syndrome type II (Mohr syndrome) is unknown. We describe a Saudi female infant with Mohr syndrome. An unclassified variant: c.2044 A>G, p.T682A in exon 12 of the GLI3 gene in a heterozygous state was identified in the infant. Mutation Taster (http://www.mutationtaster.org) considers this variant as “disease causing”. However, when the unaffected parents were tested, the father was found to have the same variant, also in a heterozygous state. Hence, the pathogenic role of this variant seems unlikely; although apparent non-penetrance remains a possibility.     

Development of simple and effective PCR based assay to detect PCCA mutation (c.425G > A) among Saudi carriers and functional study of the homozygous PCCA mutations

The aim of this study is to develop a rapid and effective method to screen for Saudi carriers of one of the most common propionic acidemia mutations (c.425G > A) and to study the functional impact of this mutation. Using allele-specific primers, we have developed a qPCR assay that clearly distinguishes heterozygotes from mutated and wild type homozygotes that overcome the dependence on labor-intensive gene sequencing. We show here that (i) qPCR rapid test has strong accuracy in detecting (c.425G > A) mutation in heterozygotes and homozygotes individuals and that the Ct-value cut-offs were estimated to be and 23.37 ± 0.04 (CV-6 %, 95 %CI-7.25) for homozygote, 25.06 ± 0.02 (CV-3.5 %, 95 %CI-7.85) for heterozygote PCCA c.425G > A mutation and 29.55 ± 0.002 (CV-11 %, 95 %CI-1.41) for PCCA wild type; (ii) the incidence of PA heterozygotes/carriers in Saudi population is about 550/100,000; (iii) skin fibroblast assays show that homozygote c.425G > A mutation induced propionyl-CoA carboxylase activity abrogation, (iv) PA patients showed an increased level of propionyl carnitine C3 in blood and 3-hydroxy propionic acid and methyl citrate in urine. Conclusion: qPCR represent an effective strategy to assess for PCCA mutation carriers in the Saudi population and we believe that will help in preventing homozygosity in the population after been implemented in pre-marriage screening program.     

Infantile-onset ascending hereditary spastic paraplegia with bulbar involvement due to the novel ALS2 mutation c.2761C > T

Recessive mutations in the alsin gene cause three clinically distinct motor neuron diseases: juvenile amyotrophic lateral sclerosis (ALS2), juvenile primary lateral sclerosis (JPLS) and infantile-onset ascending hereditary spastic paraplegia (IAHSP). A total of 23 different ALS2 mutations have been described for the three disorders so far. Most of these mutations result in a frameshift leading to a premature truncation of the alsin protein. We report the novel ALS2 truncating mutation c.2761C > T; p.R921X detected by homozygosity mapping and sequencing in two infants affected by IAHSP with bulbar involvement. The mutation c.2761C > T resides in the pleckstrin domain, a characteristic segment of guanine nucleotide exchange factors of the Rho GTPase family, which is involved in the overall neuronal development or maintenance. This study highlights the importance of using homozygosity mapping combined with candidate gene analysis to identify the underlying genetic defect as in this Saudi consanguineous family.     

A novel splicing mutation in the ABCA1 gene,causing Tangier disease and familial HDL deficiency in a large family

Tangier disease is a rare disorder of lipoprotein metabolism that presents with extremely low levels of HDL cholesterol and apoprotein A-I. It is caused by mutations in the ATP-binding cassette transporter A1 (ABCA1) gene. Clinical heterogeneity and mutational pattern of Tangier disease are poorly characterized. Moreover, also familial HDL deficiency may be caused by mutations in ABCA1 gene.ATP-binding cassette transporter A1 (ABCA1) gene mutations in a patient with Tangier disease, who presented an uncommon clinical history, and in his family were found and characterized. He was found to be compound heterozygous for two intronic mutations of ABCA1 gene, causing abnormal pre-mRNAs splicing. The novel c.1510-1G?>?A mutation was located in intron 12 and caused the activation of a cryptic splice site in exon 13, which determined the loss of 22 amino acids of exon 13 with the introduction of a premature stop codon. Five heterozygous carriers of this mutation were also found in proband's family, all presenting reduced HDL cholesterol and ApoAI (0.86?±?0.16?mmol/L and 92.2?±?10.9?mg/dL respectively), but not the typical features of Tangier disease, a phenotype compatible with the diagnosis of familial HDL deficiency. The other known mutation c.1195-27G?>?A was confirmed to cause aberrant retention of 25 nucleotides of intron 10 leading to the insertion of a stop codon after 20 amino acids of exon 11. Heterozygous carriers of this mutation also showed the clinical phenotype of familial HDL deficiency.Our study extends the catalog of pathogenic intronic mutations affecting ABCA1 pre-mRNA splicing. In a large family, a clear demonstration that the same mutations may cause Tangier disease (if in compound heterozygosis) or familial HDL deficiency (if in heterozygosis) is provided.     

Whole Exome Sequencing Identifies Recessive PKHD1 Mutations in a Chinese Twin Family with Caroli Disease

Background

Mutations in PKHD1 cause autosomal recessive Caroli disease, which is a rare congenital disorder involving cystic dilatation of the intrahepatic bile ducts. However, the mutational spectrum of PKHD1 and the phenotype-genotype correlations have not yet been fully established.

Methods

Whole exome sequencing (WES) was performed on one twin sample with Caroli disease from a Chinese family from Shandong province. Routine Sanger sequencing was used to validate the WES and to carry out segregation studies. We also described the PKHD1 mutation associated with the genotype-phenotype of this twin.

Results

A combination of WES and Sanger sequencing revealed the genetic defect to be a novel compound heterozygous genotype in PKHD1, including the missense mutation c.2507 T>C, predicted to cause a valine to alanine substitution at codon 836 (c.2507T>C, p.Val836Ala), and the nonsense mutation c.2341C>T, which is predicted to result in an arginine to stop codon at codon 781 (c.2341C>T, p.Arg781*). This compound heterozygous genotype co-segregates with the Caroli disease-affected pedigree members, but is absent in 200 normal chromosomes.

Conclusions

Our findings indicate exome sequencing can be useful in the diagnosis of Caroli disease patients and associate a compound heterozygous genotype in PKHD1 with Caroli disease, which further increases our understanding of the mutation spectrum of PKHD1 in association with Caroli disease.     

Biallelic variant in DACH1, encoding Dachshund Homolog 1, defines a novel candidate locus for recessive postaxial polydactyly type A

Polydactyly or hexadactyly is characterized by an extra digit/toe with or without a bone. Currently, variants in ten genes have been implicated in the non-syndromic form of polydactyly.DNA from a single affected individual having bilateral postaxial polydactyly was subjected to whole exome sequencing (WES), followed by Sanger sequencing. Homology modeling was performed for the identified variant and advance microscopy imaging approaches were used to reveal the localization of the DACH1 protein at the base of primary cilia.A disease-causing biallelic missense variant (c.563G > A; p.Cys188Tyr; NM_080760.5) was identified in the DACH1 gene segregating perfectly within the family. Structural analysis using homology modeling of the DACH1 protein revealed secondary structure change that might result in loss of function or influence downstream interactions. Moreover, siRNA-mediated depletion of DACH1 showed a key role of DACH1 in ciliogenesis and cilia function.This study provides the first evidence of involvement of the DACH1 gene in digits development in humans and its role in primary cilia. This signifies the importance and yet unexplored role of DACH1.     

Evidence of digenic inheritance in autoinflammation-associated genes

Familial Mediterranean fever (FMF) has traditionally been considered as a monogenic autosomal recessive disorder caused by mutations in the MEFV gene with highest incidence among Mediterranean populations. In a considerable number of patients with typical FMF, only one MEFV mutation was identified and the possibility that more than one autoinflammatory gene may be responsible for their disease was investigated. In the present study, an extensive search for possible mutations in three hereditary recurrent fever (HRF) genes was performed in 128 MEFV heterozygous Greek–Cypriots clinically diagnosed based on their phenotype with FMF-like disease from a previous study. Sequence analysis was performed for MVK, TNFRSF1A and NLRP3 genes which is also known to cause HRFs. In total, three patients were identified with heterozygous mutations and a second mutation in an autoinflammatory gene. Two patients carried a MEFV mutation and a NLRP3 mutation, and an additional third carried a MEFV mutation and a TNFRSF1A mutation. Patient 1 carried MEFV p.[Val726Ala] (NM_000243.2:c.2177T >C) and NLRP3 p.[Val198Met] (NM_001243133.1:c.592G >A) variants and patient 2 carried MEFV p.[Glu148Gln] (NM_000243.2:c.442G >C) variant which is of uncertain significance and NLRP3 p.[Arg176Trp] (NM_001243133.1:c.526C >T). Lastly, patient 3 was identified to carry MEFV p.[Met694Val] (NM_000243.2:c.2080A >G) and TNFRSF1A p.[Arg121Gln] (NM_001065.3:c.362G >A) variants. The results from this study indicate that screening of genes known to cause HRFs in patients already identified with a single MEFV mutation, can reveal quite rare but potentially causative mutational combinations at different loci. Such interaction provide further evidence for possible locus–locus interactions and phenotypes resulting from digenic inheritance.