Search

Fig. S1. Up-regulation of thyroid hormone synthesis-related genes in the BVE-Ctnnb1null thyroid tumors. Data are expressed as fold changes between BVE-Ctnnb1null and BVE tumors.

BRAFV600E mutation is the most frequent genetic alteration in papillary thyroid cancer (PTC). β-Catenin (Ctnnb1) is a key downstream component of canonical Wnt signaling pathway and is frequently overexpressed in PTC. BRAFV600E-driven tumors have been speculated to rely on Wnt/β-catenin signaling to sustain its growth, although many details remain to be elucidated. In this study, we investigated the role of β-catenin in BrafV600E-driven thyroid cancer in a transgenic mouse model. In BrafV600E mice with wild-type (WT) Ctnnb1 (BVE-Ctnnb1WT or BVE), overexpression of β-catenin was observed in thyroid tumors. In BrafV600E mice with Ctnnb1 knockout (BVE-Ctnnb1null), thyroid tumor growth was slowed with significant reduction in papillary architecture. This was associated with increased expression of genes involved in thyroid hormone synthesis, elevated 124iodine uptake, and serum T4. The survival of BVE-Ctnnb1null mice was increased by more than 50% during 14-month observation. Mechanistically, downregulation of MAPK, PI3K/Akt, and TGFβ pathways and loss of epithelial–mesenchymal transition (EMT) were demonstrated in the BVE-Ctnnb1null tumors. Treatment with dual β-catenin/KDM4A inhibitor PKF118–310 dramatically improved the sensitivity of BVE-Ctnnb1WT tumor cells to BRAFV600E inhibitor PLX4720, resulting in significant growth arrest and apoptosis in vitro, and tumor regression and differentiation in vivo. These findings indicate that β-catenin signaling plays an important role in thyroid cancer growth and resistance to BRAFV600E inhibitors. Simultaneously targeting both Wnt/β-catenin and MAPK signaling pathways may achieve better therapeutic outcome in BRAFV600E inhibitor-resistant and/or radioiodine-refractory thyroid cancer.

Supplementary Table 2. Differentially expressed vitamin D responsive genes in BVECyp24a1-wt and BVECyp24a1-null cells

Supplementary Table 3. Top 100 differentially expressed genes in BVECyp24a1-null cells

Supplementary Table 1. Primer sequences for qPCR

CYP24A1, the primary inactivating enzyme for vitamin D, is often overexpressed in human cancers, potentially neutralizing the antitumor effects of calcitriol, the active form of vitamin D. However, it is unclear whether CYP24A1 expression serves as a functional contributor versus only a biomarker for tumor progression. In this study, we investigated the role of CYP24A1 on malignant progression of a murine model of BrafV600E-induced papillary thyroid cancer (PTC). Mice harboring wild-type Cyp24a1 (BVECyp24a1-wt) developed PTC at 5 weeks of age. Mice harboring a homozygous deletion of Cyp24a1 (BVECyp24a1-null) exhibited a 4-fold reduction in tumor growth. Notably, we found the tumorigenic potential of BVECyp24a1-null-derived tumor cells to be nearly abolished in immunocompromised nude mice. This phenotype was associated with downregulation of the MAPK, PI3K/Akt, and TGFβ signaling pathways and a loss of epithelial–mesenchymal transition (EMT) in BVECyp24a1-null cells, associated with downregulation of genes involved in EMT, tumor invasion, and metastasis. While calcitriol treatment did not decrease cell proliferation in BVECyp24a1-null cells, it strengthened antitumor responses to the BRAFV600E inhibitor PLX4720 in both BVECyp24a1-null and BVECyp24a1-wt cells. Our findings offer direct evidence that Cyp24a1 functions as an oncogene in PTC, where its overexpression activates multiple signaling cascades to promote malignant progression and resistance to PLX4720 treatment. Cancer Res; 77(8); 2161–72. ©2017 AACR.

Abstract Background Familial transthyretin (TTR) amyloidosis (ATTR) is an autosomal dominant disease with significant phenotypic heterogeneity. Its prevalence in Saudi Arabia has not previously been investigated. An existing exome variant database of Saudi individuals, sequenced to globally investigate rare diseases in the population, was mined for TTR variants and filtered for missense mutations resulting in single amino acid changes. A total of 13,906 Saudi exomes from unrelated individuals were analyzed blindly. Results Three TTR variants known to be associated with ATTR amyloidosis were identified. Additionally, three novel TTR mutations were identified. Structural analysis of the three novel variants suggests that at least two could be amyloidogenic. The most common variant associated with amyloidosis was p.Val142Ile (allele frequency 0.001). Further investigation of these variants and their translation to clinical practice may help to diagnose, monitor, and manage patients with ATTR amyloidosis. Conclusion Multiple TTR variants potentially associated with systemic ATTR amyloidosis were identified in the Saudi population. Early diagnosis and intervention, facilitated by familial genetic testing of patients with ATTR amyloidosis, may benefit in the management of this disease. Early diagnosis could be enhanced through inclusion of ATTR variants in existing population-based screening programs.

BRAFV600E mutation is the most frequent genetic alteration in papillary thyroid cancer (PTC). β-Catenin (Ctnnb1) is a key downstream component of canonical Wnt signaling pathway and is frequently overexpressed in PTC. BRAFV600E-driven tumors have been speculated to rely on Wnt/β-catenin signaling to sustain its growth, although many details remain to be elucidated. In this study, we investigated the role of β-catenin in BrafV600E-driven thyroid cancer in a transgenic mouse model. In BrafV600E mice with wild-type (WT) Ctnnb1 (BVE-Ctnnb1WT or BVE), overexpression of β-catenin was observed in thyroid tumors. In BrafV600E mice with Ctnnb1 knockout (BVE-Ctnnb1null), thyroid tumor growth was slowed with significant reduction in papillary architecture. This was associated with increased expression of genes involved in thyroid hormone synthesis, elevated 124iodine uptake, and serum T4. The survival of BVE-Ctnnb1null mice was increased by more than 50% during 14-month observation. Mechanistically, downregulation of MAPK, PI3K/Akt, and TGFβ pathways and loss of epithelial–mesenchymal transition (EMT) were demonstrated in the BVE-Ctnnb1null tumors. Treatment with dual β-catenin/KDM4A inhibitor PKF118–310 dramatically improved the sensitivity of BVE-Ctnnb1WT tumor cells to BRAFV600E inhibitor PLX4720, resulting in significant growth arrest and apoptosis in vitro, and tumor regression and differentiation in vivo. These findings indicate that β-catenin signaling plays an important role in thyroid cancer growth and resistance to BRAFV600E inhibitors. Simultaneously targeting both Wnt/β-catenin and MAPK signaling pathways may achieve better therapeutic outcome in BRAFV600E inhibitor-resistant and/or radioiodine-refractory thyroid cancer.

Background The classification and pathogenic basis of juvenile idiopathic arthritis (JIA) are a subject of some controversy. Essentially, JIA is an exclusion diagnosis that represents a phenotypically heterogeneous group of arthritis of unknown origin. Familial aggregation of JIA supports the concept of genetic influence in the pathogenesis of JIA. Objective To present the spectrum of laccase domain-containing 1 (LACC1)-associated juvenile arthritis with clinical, biochemical, and molecular genetic data of a cohort of 43 patients, including 11 previously unpublished cases. Methods We studied 11 patients with different categories of juvenile idiopathic arthritis from 5 consanguineous families, all from Saudi Arabia, except 2 patients who were of Jordanian ethnicity. Whole-exome sequencing was used to identify the disease-causing variant of LACC1. We also reviewed the clinical spectrum and molecular genetic data of previously published cases of LACC1-associated juvenile arthritis. Results This study describes 43 (29 females, 14 males) patients from consanguineous multiplex families. Most of the included patients were of Arab origin with 86% having early onset disease. The most frequent categories were systemic (19 patients) and rheumatoid factor-negative polyarticular (19 patients). Thirty-seven (86%) had progressive erosive arthritis and 10 (23.3%) had persistent limb lymphedema. None of the patients had features of macrophage activation syndrome. Genetic analysis confirmed LACC1 variant in all patients; 22 patients had common founder mutation (LACC1: c.850T > C,p.C284R), while the others showed different LACC1 variants. All patients were treated aggressively with methotrexate and sequential biologic agents. Most of them showed a poor response to treatment. Conclusion This report expands the pathogenic variants of LACC1 and the clinical spectrum associated with this genetic subset of juvenile arthritis. The predominance of autosomal-recessive inheritance and strong genetic evidence allowed us to propose LACC1-associated juvenile arthritis as a distinct disorder.

Objective: Vitamin D–dependent rickets type 2A (VDDR2A) is a rare autosomal recessive disorder caused by mutations in the vitamin D receptor gene (VDR), leading to end-organ resistance to 1,25-dihydroxyvitamin D3 (1,25lOH]2D3). The objective of this study was to investigate VDR mutations in 11 patients from 8 Turkish-Arab families. Methods: All coding exons and intron-exon boundaries of the VDR gene were amplified by polymerase chain reaction from peripheral leukocyte DNA and sequenced. The effect of splice-site mutations on mRNA splicing was evaluated by a customized VDR mini-gene assay. Results: Homozygous VDR mutations were found in all the patients, including four novel mutations: c.473G>T (p.R158L), c.1-4A>G (IVS3-2A>G), c.755+1G>T, and c.352_356delGACAG (p.D118Sfs*7). The c.1-4A>G mutation was located in the canonical splice acceptor site and 4 base pairs from the original ATG start codon. The mutation resulted in both complete (60% of transcripts) and partial exon 4 skipping (15% of transcripts). The latter was due to activation of a cryptic splice acceptor site and did not disrupt the open reading frame. Both c.755+1G>T and c.352_356delGACAG resulted in frameshifts and a premature stop codon. Clinically, all the patients required continued treatment, except for patient IV-3, who presented with alopecia, hypocalcemia, and increased 1,25(OH)2D3 at 1.5 years of age as a result of the c.1-4A>G mutation. He stopped taking medication at 6 years of age and still maintained normal height and biochemical profile. Conclusion: We have identified four novel VDR mutations. Although canonical splice-site mutations cause premRNA splicing errors that usually lead to a severe disease phenotype, mild disease can also occur due to activation of a cryptic splice site. Abbreviations: 1,25(OH)2D3 = 1,25-dihydroxyvitamin D3 (calcitriol); 25OHD3 = 25-hydroxyvitamin D3; PCR = polymerase chain reaction; PTH = parathyroid hormone; VDDR2A = vitamin D–dependent rickets type 2A; VDR = vitamin D receptor

It is well documented that drug responses are related to Absorption, Distribution, Metabolism, and Excretion (ADME) characteristics of individual patients. Several studies have identified genetic variability in pharmacogenes, that are either directly responsible for or are associated with ADME, giving rise to individualized treatments. Our objective was to provide a comprehensive overview of pharmacogenetic variation in the Saudi population. We mined next generation sequencing (NGS) data from 11,889 unrelated Saudi nationals, to determine the presence and frequencies of known functional SNP variants in 8 clinically relevant pharmacogenes (CYP2C9, CYP2C19, CYP3A5, CYP4F2, VKORC1, DPYD, TPMT and NUDT15), recommended by the Clinical Pharmacogenetics Implementation Consortium (CPIC), and collectively identified 82 such star alleles. Functionally significant pharmacogenetic variants were prevalent especially in CYP genes (excluding CYP3A5), with 10–44.4% of variants predicted to be inactive or to have decreased activity. In CYP3A5, inactive alleles (87.5%) were the most common. Only 1.8%, 0.7% and 0.7% of NUDT15, TPMT and DPYD variants respectively, were predicted to affect gene activity. In contrast, VKORC1 was found functionally, to be highly polymorphic with 53.7% of Saudi individuals harboring variants predicted to result in decreased activity and 31.3% having variants leading to increased metabolic activity. Furthermore, among the 8 pharmacogenes studied, we detected six rare variants with an aggregated frequency of 1.1%, that among several other ethnicities, were uniquely found in Saudi population. Similarly, within our cohort, the 8 pharmacogenes yielded forty-six novel variants predicted to be deleterious. Based upon our findings, 99.2% of individuals from the Saudi population carry at least one actionable pharmacogenetic variant.

Context Hypophosphatemic rickets (HR) is a group of rare hereditary renal phosphate wasting disorders caused by mutations in PHEX, FGF23, DMP1, ENPP1, CLCN5, SLC9A3R1, SLC34A1, or SLC34A3. Objective A large kindred with 5 HR patients was recruited with dominant inheritance. The study was undertaken to investigate underlying genetic defects in HR patients. Design Patients and their family members were initially analyzed for PHEX and FGF23 mutations using polymerase chain reaction sequencing and copy number analysis. Exome sequencing was subsequently performed to identify novel candidate genes. Results PHEX and FGF23 mutations were not detected in the patients. No copy number variation was observed in the genome using CytoScan HD array analysis. Mutations in DMP1, ENPP1, CLCN5, SLC9A3R1, SLC34A1, or SLC34A3 were also not found by exome sequencing. A novel c.979–96 T>A mutation in the SGK3 gene was found to be strictly segregated in a heterozygous pattern in patients and was not present in normal family members. The mutation is located 1 bp downstream of a highly conserved adenosine branch point, resulted in exon 13 skipping and in-frame deletion of 29 amino acids, which is part of the protein kinase domain and contains a Thr-320 phosphorylation site that is required for its activation. Protein tertiary structure modelling showed significant structural change in the protein kinase domain following the deletion. Conclusions The c.979–96 T>A splice mutation in the SGK3 gene causes exon 13 skipping and deletion of 29 amino acids in the protein kinase domain. The SGK3 mutation may cause autosomal dominant HR.

Context X-linked hypophosphatemic rickets (XLH) is caused by inactivating mutations in the PHEX gene and is the most common form of hereditary rickets. The splice-site mutations account for 17% of all reported PHEX mutations. The functional consequence of these splice-site mutations has not been systemically investigated. Objective The current study was undertaken to functionally annotate previously reported 22 splice-site mutations in the PHEX gene. Methods PHEX mini-genes with different splice-site mutations were created by site-directed mutagenesis and expressed in HEK293 cells. The mRNA transcripts were analyzed by RT-PCR, cloning, and sequencing. Results These splicing mutations led to a variety of consequences, including exon skipping, intron retention, and activation of cryptic splice sites. Among 22 splice-site mutations, exon skipping was the most common event accounting for 73% (16/22). Non-canonical splice-site mutations could result in splicing errors to the same extent as canonical splice-site mutations such as c.436+3G>C, c.436+4A>C, c.436+6T>C, c.437-3C>G, c.850-3C>G, c.1080-3C>A, c.1482+5G>C, c.1586+6T>C, c.1645+5G>A, c.1645+6T>C, c.1701-16T>A, c.1768+5G>A, and c.1899+5G>A. Interestingly, non-canonical (c.436+6T>C and c.1586+6T>C) and canonical splice-site mutations (c.1769-1G>C) could generate partial splicing errors (both wild-type and mutant transcripts were detected), resulting in incomplete inactivation of PHEX gene, which may explain the mild disease phenotype reported previously, providing evidence of genotype-phenotype correlation. c.1645C>T (p.R549*) had no impact on pre-mRNA splicing although it is located next to canonical splice donor site GT. Conclusions Exon skipping is the most common outcome due to splice-site mutations. Both canonical and non-canonical splice-site mutations can result in either severe or mild RNA splicing defects, contributing to phenotype heterogeneity. Non-canonical splice-site mutations should not be overlooked in genetic screening especially those located within 50 bp from canonical splice site.

Piebaldism is a rare autosomal dominant disorder characterized by congenital patches of depigmented skin and hair that lack melanocytes (MIM172800). The KIT gene, mapped to chromosome 4q12, is one of the common causative genes in autosomal dominant piebaldism and accounts for 75% of patients with piebaldism.1 Various studies in mice found that mutations involving the murine KIT locus led to dominant white spotting and abnormalities in hematopoiesis and germ cell development. More than 47 mutations have been described that are associated with piebaldism.2, 3, 4 Here we describe a Saudi nonconsanguineous family with an autosomal dominant form of inheritance, with marked phenotypic and genotypic variability. All the affected individuals had widespread depigmented patches involving more than 70% to 80% body surface area with almost total poliosis. Although all these features presented since their birth, the father noticed heterochromia iridis as an adult. We used next-generation sequencing with 214 dermatologic gene panels (DGP)5 to identify pathogenic mutations in the KIT gene that underlies the autosomal dominant form of piebaldism in this family.

Wolf–Hirschhorn syndrome (WHS) is a genomic disorder with a recognizable dysmorphology profile caused by hemizygosity at 4p16.3. Previous attempts have failed to map the minimal critical locus to a single gene, leaving open the possibility that the core phenotypic components of the syndrome are caused by the combined haploinsufficiency of multiple genes. Clinical exome sequencing and “reverse” phenotyping. We identified two patients with de novo truncating variants in WHSC1, which maps to the WHS critical locus. The phenotype of these two individuals is consistent with WHS, which suggests that haploinsufficiency of WHSC1 is sufficient to recapitulate the core phenotype (characteristic facies, and growth and developmental delay) of this classic microdeletion syndrome. Our study expands the list of microdeletion syndromes that are solved at the single-gene level, and establishes WHSC1 as a disease gene in humans. Given the severe nature of the reported variants, the full phenotypic expression of WHSC1 may be further expanded by future reports of milder variants.

Context Congenital hypothyroidism (CH) is caused by mutations in the genes for thyroid hormone synthesis. In our previous investigation of CH patients, approximately 53% of patients had mutations in either coding exons or canonical splice sites of causative genes. Noncanonical splice-site variants in the intron were detected but their pathogenic significance was not known. Objective This work aims to evaluate noncanonical splice-site variants on pre–messenger RNA (pre-mRNA) splicing of CH-causing genes. Methods Next-generation sequencing data of 55 CH cases in 47 families were analyzed to identify rare intron variants. The effects of variants on pre-mRNA splicing were investigated by minigene RNA-splicing assay. Results Four intron variants were found in 3 patients: solute carrier family 26 member 4 (SLC26A4) c.1544+9C>T and c.1707+94C>T in one patient, and solute carrier family 5 member 5 (SLC5A5) c.970-48G>C and c.1652-97A>C in 2 other patients. The c.1707+94C>T and c.970-48G>C caused exons 15 and 16 skipping, and exon 8 skipping, respectively. The remaining variants had no effect on RNA splicing. Furthermore, we analyzed 28 previously reported noncanonical splice-site variants (4 in TG and 24 in SLC26A4). Among them, 15 variants (~ 54%) resulted in aberrant splicing and 13 variants had no effect on RNA splicing. These data were compared with 3 variant-prediction programs (FATHMM-XF, FATHMM-MKL, and CADD). Among 32 variants, FATHMM-XF, FATHMM-MKL, and CADD correctly predicted 20 (63%), 17 (53%), and 26 (81%) variants, respectively. Conclusion Two novel deep intron mutations have been identified in SLC26A4 and SLC5A5, bringing the total number of solved families with disease-causing mutations to approximately 45% in our cohort. Approximately 46% (13/28) of reported noncanonical splice-site mutations do not disrupt pre-mRNA splicing. CADD provides highest prediction accuracy of noncanonical splice-site variants.

BackgroundFamilial transthyretin (TTR) amyloidosis (ATTR) is an autosomal dominant disease with significant phenotypic heterogeneity. Its prevalence in Saudi Arabia has not previously been investigated. An existing exome variant database of Saudi individuals, sequenced to globally investigate rare diseases in the population, was mined for TTR variants and filtered for missense mutations resulting in single amino acid changes. A total of 13,906 Saudi exomes from unrelated individuals were analyzed blindly. ResultsThree TTR variants known to be associated with ATTR amyloidosis were identified. Additionally, three novel TTR mutations were identified. Structural analysis of the three novel variants suggests that at least two could be amyloidogenic. The most common variant associated with amyloidosis was p.Val142Ile (allele frequency 0.001). Further investigation of these variants and their translation to clinical practice may help to diagnose, monitor and manage patients with ATTR amyloidosis.ConclusionMultiple TTR variants potentially associated with systemic ATTR amyloidosis were identified in the Saudi population. Early diagnosis and intervention, facilitated by familial genetic testing of patients with ATTR amyloidosis, may benefit in the management of this disease. Early diagnosis could be enhanced through inclusion of ATTR variants in existing population-based screening programs.

Membrane trafficking is a complex, essential process in eukaryotic cells responsible for protein transport and processing. Deficiencies in vacuolar protein sorting (VPS) proteins, key regulators of trafficking, cause abnormal intracellular segregation of macromolecules and organelles and are linked to human disease. VPS proteins function as part of complexes such as the homotypic fusion and vacuole protein sorting (HOPS) tethering complex, composed of VPS11, VPS16, VPS18, VPS33A, VPS39 and VPS41. The HOPS-specific subunit VPS41 has been reported to promote viability of dopaminergic neurons in Parkinson’s disease but to date has not been linked to human disease. Here, we describe five unrelated families with nine affected individuals, all carrying homozygous variants in VPS41 that we show impact protein function. All affected individuals presented with a progressive neurodevelopmental disorder consisting of cognitive impairment, cerebellar atrophy/hypoplasia, motor dysfunction with ataxia and dystonia, and nystagmus. Zebrafish disease modelling supports the involvement of VPS41 dysfunction in the disorder, indicating lysosomal dysregulation throughout the brain and providing support for cerebellar and microglial abnormalities when vps41 was mutated. This provides the first example of human disease linked to the HOPS-specific subunit VPS41 and suggests the importance of HOPS complex activity for cerebellar function., Sanderson et al. describe ten individuals with biallelic variants in the gene encoding the vacuolar protein sorting protein VPS41, presenting with cerebellar dysfunction and abnormal membrane-bound vesicles. Mutant zebrafish recapitulate the cerebellar dysfunction phenotype, further implicating VPS41 as a human disease gene.

Background Hereditary Multiple Exostoses (HME), also known as Multiple Osteochondromas (MO) is a rare genetic disorder characterized by multiple benign cartilaginous bone tumors, which are caused by mutations in the genes for exostosin glycosyltransferase 1 (EXT1) and exostosin glycosyltransferase 2 (EXT2). The genetic defects have not been studied in the Saudi patients. Aim of study We investigated mutation spectrum of EXT1 and EXT2 in 22 patients from 17 unrelated families. Methods Genomic DNA was extracted from peripheral leucocytes. The coding regions and intron–exon boundaries of both EXT1 and EXT2 genes were screened for mutations by PCR-sequencing analysis. Gross deletions were analyzed by MLPA analysis. Results EXT1 mutations were detected in 6 families (35%) and 3 were novel mutations: c.739G > T (p. E247*), c.1319delG (p.R440Lfs*4), and c.1786delA (p.S596Afs*25). EXT2 mutations were detected in 7 families (41%) and 3 were novel mutations: c.541delG (p.D181Ifs*89), c.583delG (p.G195Vfs*75), and a gross deletion of approximately 10 kb including promoter and exon 1. Five patients from different families had no family history and carried de novo mutations (29%, 5/17). No EXT1 and EXT2 mutations were found in the remaining four families. In total, EXT1 and EXT2 mutations were found in 77% (13/17) of Saudi HME patients. Conclusion EXT1 and EXT2 mutations contribute significantly to the pathogenesis of HME in the Saudi population. In contrast to high mutation rate in EXT 1 (65%) and low mutation rate in EXT2 (25%) in other populations, the frequency of EXT2 mutations are much higher (41%) and comparable to that of EXT1 among Saudi patients. De novo mutations are also common and the six novel EXT1/EXT2 mutations further expands the mutation spectrum of HME.

ContextVitamin D-dependent rickets type 1A (VDDR1A) is a rare autosomal recessively inherited disorder due to loss-of-function mutations in the CYP27B1 gene. CYP27B1 encodes an enzyme of 25-hydroxyvitamin D-1α-hydroxylase for converting inactive 25-OHD to biologically active 1,25-(OH)2D.ObjectiveTo identify underlying genetic defects in patients with VDDR1A.MethodsTwelve patients from 7 Turkish and 2 Saudi families were investigated. The coding exons and intron-exon boundaries of the CYP27B1 gene were amplified by Polymerase Chain Reaction (PCR) from peripheral lymphocyte DNA. PCR products were directly sequenced. The consequences of c.590G > A mutation were analyzed by in silico and functional analysis.ResultsCYP27B1 mutations were identified in all the patients. Two novel mutations were identified in two separate families: c.171delG (family 7) and c.398_400dupAAT (family 8). The intra-exon deletion of c.171delG resulted in a frameshift and premature stop codon 20 amino acids downstream from the mutation (p.L58Cfs∗20). The intra-exon duplication of c.398_400dupAAT generated a premature stop codon at the mutation site (p.W134∗). A missense c.590G > A (p.G197D) mutation was found in a patient from family 4 and caused a defect in pre-mRNA splicing. As a result, two populations of transcripts were detected: the majority of them with intron 3 retention (83%), and the minority (17%) being properly spliced transcripts with about 16% of wild-type enzymatic activity. The remaining nine patients from six families carried a previously reported c.1319_1325dupCCCACCC (F443Pfs∗24) mutation. Clinically, all the patients need continued calcitriol treatment, which was consistent with inactivation of 25-hydroxy vitamin D1α-hydroxylase activity.ConclusionTwo novel frameshift CYP27B1 mutations were identified and predicted to inactivate 25-hydroxyvitamin D-1α-hydroxylase. The loss of enzymatic activity by c.590G > A missense mutation was mainly caused by aberrant pre-mRNA splicing.

Context Congenital hypothyroidism (CH) is the most common neonatal endocrine disorder, affecting one in 3000 to 4000 newborns. Since the introduction of a newborn screening program in 1988, more than 300 cases have been identified. The underlying genetic defects have not been systematically studied. Objective To identify the mutation spectrum of CH-causing genes. Methods Fifty-five patients from 47 families were studied by next-generation exome sequencing. Results Mutations were identified in 52.7% of patients (29 of 55) in the following 11 genes: TG, TPO, DUOX2, SLC26A4, SLC26A7, TSHB, TSHR, NKX2-1, PAX8, CDCA8, and HOXB3. Among 30 patients with thyroid dyshormonogenesis, biallelic TG mutations were found in 12 patients (40%), followed by biallelic mutations in TPO (6.7%), SLC26A7 (6.7%), and DUOX2 (3.3%). Monoallelic SLC26A4 mutations were found in two patients, one of them coexisting with two tandem biallelic deletions in SLC26A7. In 25 patients with thyroid dysgenesis, biallelic mutations in TSHR were found in six patients (24%). Biallelic mutations in TSHB, PAX 8, NKX2-1, or HOXB3 were found once in four different patients. A monoallelic CDCA8 mutation was found in one patient. Most mutations were novel, including three TG, two TSHR, and one each in DUOX2, TPO, SLC26A7, TSHB, NKX2-1, PAX8, CDCA8, and HOXB3. SLC26A7 and HOXB3 were novel genes associated with thyroid dyshormonogenesis and dysgenesis, respectively. Conclusions TG and TSHR mutations are the most common genetic defects in Saudi patients with CH. The prevalence of other disease-causing mutations is low, reflecting the consanguineous nature of the population. SLC26A7 mutations appear to be associated with thyroid dyshormonogenesis.

Genome-wide association studies have been successful in elucidating the genetic basis of colorectal cancer, but there remains unexplained variability in genetic risk. To identify new risk variants and to confirm reported associations, we conducted a genome-wide association study in 1,701 colorectal cancer cases and 14,082 cancer-free controls from the Finnish population. A total of 9,068,015 genetic variants were imputed and tested, and 30 promising variants were studied in additional 11,647 cases and 12,356 controls of European ancestry. The previously reported association between the single-nucleotide polymorphism rs992157 (2q35) and colorectal cancer was independently replicated (p=2.08x10-4; OR, 1.14; 95% CI, 1.06-1.23), and it was genome-wide significant in combined analysis (p=1.50x10-9; OR, 1.12; 95% CI, 1.08-1.16). Variants at 2q35, 6p21.2, 8q23.3, 8q24.21, 10q22.3, 10q24.2, 11q13.4, 11q23.1, 14q22.2, 15q13.3, 18q21.1, 20p12.3, and 20q13.33 were associated with colorectal cancer in the Finnish population (false discovery rate

CYP24A1, the primary inactivating enzyme for vitamin D, is often overexpressed in human cancers, potentially neutralizing the antitumor effects of calcitriol, the active form of vitamin D. However, it is unclear whether CYP24A1 expression serves as a functional contributor versus only a biomarker for tumor progression. In this study, we investigated the role of CYP24A1 on malignant progression of a murine model of BrafV600E-induced papillary thyroid cancer (PTC). Mice harboring wild-type Cyp24a1 (BVECyp24a1-wt) developed PTC at 5 weeks of age. Mice harboring a homozygous deletion of Cyp24a1 (BVECyp24a1-null) exhibited a 4-fold reduction in tumor growth. Notably, we found the tumorigenic potential of BVECyp24a1-null-derived tumor cells to be nearly abolished in immunocompromised nude mice. This phenotype was associated with downregulation of the MAPK, PI3K/Akt, and TGFβ signaling pathways and a loss of epithelial–mesenchymal transition (EMT) in BVECyp24a1-null cells, associated with downregulation of genes involved in EMT, tumor invasion, and metastasis. While calcitriol treatment did not decrease cell proliferation in BVECyp24a1-null cells, it strengthened antitumor responses to the BRAFV600E inhibitor PLX4720 in both BVECyp24a1-null and BVECyp24a1-wt cells. Our findings offer direct evidence that Cyp24a1 functions as an oncogene in PTC, where its overexpression activates multiple signaling cascades to promote malignant progression and resistance to PLX4720 treatment. Cancer Res; 77(8); 2161–72. ©2017 AACR.

Hereditary sensory autonomic neuropathy type IV (HSAN-IV) is a rare autosomal recessive disorder that usually begins in infancy and is characterized by anhidrosis, insensitivity to noxious stimuli leading to self-mutilating behavior, and intellectual disability. HSAN-IV is caused by mutations in the neurotrophic tyrosine kinase receptor type 1 gene, NTRK1, encoding the high-affinity receptor of nerve growth factor (NGF) which maps to chromosome 1q21-q22. Patients with HSAN-IV lack all NGF-dependent neurons, the primary afferents and sympathetic postganglionic neurons leading to lack of pain sensation and the presence of anhidrosis, respectively. Herein, we report nine patients from nine unrelated families with HSAN-IV due to various mutations in NTRK1, five of which are novel. These are three missense and two nonsense mutations distributed in various domains of NTRK1 involved in binding of NGF. The affected patients had variable intellectual deficits, and some had delayed diagnosis of HSAN-IV. In addition to being the first report of HSAN-IV from the Arabian Peninsula, this report expands the mutational spectrum of patients with NTRK1 mutations and provides further insights for molecular and clinical diagnosis.

Publisher's version (útgefin grein)., In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44–66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding., This paper is the work of the ROHgen consortium. We thank the Sigma T2D Consortium, whose members are detailed in Supplementary Note 3. We thank the UK Biobank Resource, approved under application 19655; we acknowledge funding from the UK Medical Research Council Human Genetics Unit and MRC Doctoral Training Programme in Precision Medicine. We also thank Neil Robertson, Wellcome Trust Centre for Human Genetics, Oxford, for use of his author details management software, Authorial. Finally, we thank all the participants, researchers and funders of ROHgen cohorts. Cohort-specific acknowledgements are in Supplementary Data 2; personal acknowledgements and disclosures are in Supplementary Note 2. We thank Rachel Edwards for administrative assistance.

Hereditary sensory and autonomic neuropathies (HSANs) are a clinically and genetically heterogeneous group of disorders involving various sensory and autonomic dysfunctions. The most common symptoms of HSANs include loss of sensations of pain and temperature that frequently lead to chronic ulcerations in the feet and hands of the patient. In this case study, we present the clinical features and genetic characteristics of two affected individuals from two unrelated Saudi families presenting mutilating sensory loss and spastic paraplegia. We employed homozygosity mapping and exome sequencing which is an efficient strategy to characterize the recessive genes, thus obtaining a rapid molecular diagnosis for genetically heterogeneous disorders like HSAN. Subsequently, a nonsense mutation (c.926 C>G; p.S309⁎) in FAM134B was identified. In addition, we confirmed that the mutant FAM134B transcripts were reduced in these patients presumably disrupting the receptors of the degradative endoplasmic reticulum pathways that facilitate the autophagy processes.

Background Voltage-gated potassium channels are highly diverse proteins representing the most complex class of voltage-gated ion channels from structural and functional perspectives. Deficiency of these channels usually results in various human disorders. Objectives To describe a novel autosomal recessive syndrome associated with KCNA4 deficiency leading to congenital cataract, abnormal striatum, intellectual disability and attention deficit hyperactivity disorder. Methods We used SNP arrays, linkage analyses, autozygosity mapping, whole-exome sequencing, RT-PCR and two-electrode voltage-clamp recording. Results We identified a missense variant (p.Arg89Gln) in KCNA4 in four patients from a consanguineous family manifesting a novel syndrome of congenital cataract, abnormal striatum, intellectual disability and attention deficit hyperactivity disorder. The variant was fully segregated with the disease and absent in 747 ethnically matched exomes. Xenopus oocytes were injected with human Kv1.4 wild-type mRNA, R89Q and WT/R89Q channels. The wild type had mean current amplitude that was significantly greater than those recorded from the cells expressing the same amount of mutant mRNA. Co-expression of the wild type and mutant mRNAs resulted in mean current amplitude that was significantly different from that of the wild type. RT-PCR indicated that KCNA4 is present in mouse brain, lens and retina. KCNA4 interacts with several molecules including synaptotagmin I, DLG1 and DLG2. The channel co-localises with cholinergic amacrine and rod bipolar cells in rats and is widely distributed in the central nervous system. Based on previous studies, the channel is highly expressed in outer retina, rod inner segments, hippocampus and concentrated in axonal membranes. Conclusion KCNA4 (Kv1.4) is implicated in a novel syndrome characterised by striatal thinning, congenital cataract and attention deficit hyperactivity disorder. Our study highlights potassium channels' role in ocular and neuronal genetics.

Article, Background: Observational studies have associated adiposity with an increased risk of colorectal cancer (CRC). However, such studies do not establish a causal relationship. To minimise bias from confounding we performed a Mendelian randomisation (MR) analysis to examine the relationship between adiposity and CRC. Methods: We used SNPs associated with adult body mass index (BMI), waist-hip ratio (WHR), childhood obesity and birth weight as instrumental variables in a MR analysis of 9254 CRC cases and 18 386 controls. Results: In the MR analysis, the odds ratios (ORs) of CRC risk per unit increase in BMI, WHR and childhood obesity were 1.23 (95% CI: 1.02–1.49, P=0.033), 1.59 (95% CI: 1.08–2.34, P=0.019) and 1.07 (95% CI: 1.03–1.13, P=0.018), respectively. There was no evidence for association between birth weight and CRC (OR=1.22, 95% CI: 0.89–1.67, P=0.22). Combining these data with a concurrent MR-based analysis for BMI and WHR with CRC risk (totalling to 18 190 cases, 27 617 controls) provided increased support, ORs for BMI and WHR were 1.26 (95% CI: 1.10–1.44, P=7.7 × 10−4) and 1.40 (95% CI: 1.14–1.72, P=1.2 × 10−3), respectively. Conclusions: These data provide further evidence for a strong causal relationship between adiposity and the risk of developing CRC highlighting the urgent need for prevention and treatment of adiposity., published version, peerReviewed

Background Autosomal recessive congenital ichthyosis (ARCI) is a rare disorder of keratinization. Infants (10–15%) born with this condition are encapsulated in hyperkeratotic membrane covering the entire body and are called “collodion babies.” So far, mutations in nine different genes have been identified as causative and implicated in the pathogenesis of the clinically and genetically heterogeneous group of ARCI disorders. Among these, TGM1 is the gene most commonly mutated in ARCI. Methods We identified 11 patients from five consanguineous but unrelated families affected by ARCI. These patients manifested thick adherent polygonal large scales all over the body. All six patients with TGM1 mutations were born with collodion membrane and had ectropion and eclabium, while none of the patients with ABCA12 mutations had these features. Molecular investigations were performed using the combined approach of homozygosity mapping and Sanger sequencing. Results Here we report two novel mutations c.397_398insAGTATGAGTA (p.Tyr136Ter); c.977-978delCT (p.Ser326Cysfs*8) in TGM1 in three different, unrelated Saudi families and one novel mutation c.6900C>A (p.Phe2300Leu) and one reported mutation c.3470C>T (p.Ser1157Leu) in the ABCA12 gene in two unrelated Saudi families with ARCI. Conclusions The identification of these homozygous variants using combined approaches of homozygosity mapping with direct sequencing are the disease causing mutations in these families. Furthermore, these findings are essential for the genetic diagnostic and prognostic workup with ARCI in Saudi patients.

The genetic susceptibility to acquiring low high density lipoprotein-cholesterol (LHDLC) levels is not completely elucidated yet. In this study, we performed a common variant association study for harboring this trait in ethnic Arabs. We employed the Affymetrix high-density Axiom Genome-Wide ASI Array (Asian population) providing a coverage of 598,000 single nucleotide variations (SNPs) to genotype 5495 individuals in a two-phase study involving discovery and validation sets of experiments. The rs1800775 [1.31 (1.22-1.42); p = 3.41E-12] in the CETP gene and rs359027 [1.26 (1.16-1.36); p = 2.55E-08] in the LMCD1 gene were significantly associated with LHDLC levels. Furthermore, rs3104435 [1.26 (1.15-1.38); p = 1.19E-06] at the MATN1 locus, rs9835344 [1.16 (1.08-1.26); p = 8.75E-06] in the CNTN6 gene, rs1559997 [1.3 (1.14-1.47); p = 9.48E-06] in the SDS gene and rs1670273 [1.2 (1.1-1.31); p = 4.81E-06] in the DMN/SYNM gene exhibited suggestive association with the disorder. Seven other variants including rs1147169 in the PLCL1 gene, rs10248618 in the DNAH11, rs476155 in the GLIS3, rs7024300 in the ABCA1, intergenic rs10836699, rs11603691 in P2RX3 and rs750134 in CORO1C gene exhibited borderline protective properties. Validation and joint meta-analysis resulted in rs1800775, rs3104435 and rs359027 retaining their predisposing properties, while rs10836699 and rs11603691 showed protective properties. Our data show several predisposing variants across the genome for LHDLC levels in ethnic Arabs.

Background Inherited cystic kidney disorders are a common cause of end-stage renal disease. Over 50 ciliopathy genes, which encode proteins that influence the structure and function of the primary cilia, are implicated in cystic kidney disease. Methods To define the phenotype and genotype of cystic kidney disease in fetuses and neonates, we correlated antenatal ultrasound examination and postnatal renal ultrasound examination with targeted exon sequencing, using a renal gene panel. A cohort of 44 families in whom antenatal renal ultrasound scanning findings in affected cases included bilateral cystic kidney disease, echogenic kidneys or enlarged kidneys was investigated. Results In this cohort, disease phenotypes were severe with 36 cases of stillbirth or perinatal death. Extra renal malformations, including encephalocele, polydactyly and heart malformations, consistent with ciliopathy phenotypes, were frequently detected. Renal gene panel testing identified causative mutations in 21 out of 34 families (62%), where patient and parental DNA was available. In the remaining 10 families, where only parental DNA was available, 7 inferred causative mutations were found. Together, mutations were found in 12 different genes with a total of 13 novel pathogenic variants, including an inferred novel variant in NEK8 . Mutations in CC2D2A were the most common cause of an antenatal cystic kidney disease and a suspected ciliopathy in our cohort. Conclusions In families with ciliopathy phenotypes, mutational analysis using a targeted renal gene panel allows a rapid molecular diagnosis and provides important information for patients, parents and their physicians.

Papillary thyroid carcinoma (PTC) accounts for >80% thyroid malignancies, and BRAF(V600E) mutation is frequently found in >40% PTC. Interleukin-12 (IL-12) is a proinflammatory heterodimeric cytokine with strong antitumor activity. It is not known whether IL-12 immunotherapy is effective against Braf(V600E)-induced PTC. In the present study, we investigated the effectiveness of IL-12 immunotherapy against Braf(V600E)-induced PTC in LSL-Braf(V600E)/TPO-Cre mice. LSL-Braf(V600E)/TPO-Cre mice were created for thyroid-specific expression of Braf(V600E) under the endogenous Braf promoter, and spontaneous PTC developed at about 5 weeks of age. The mice were subjected to two treatment regimens: (1) weekly intramuscular injection of 50 μg plasmid DNA expressing a single-chain IL-12 fusion protein (scIL-12/CMVpDNA), (2) daily intraperitoneal injection of mouse recombinant IL-12 protein (mrIL-12, 100 ng per day). The role of T cells, natural killer (NK) cells, and transforming growth factor-β (TGF-β) in IL-12-mediated antitumor effects was determined by a (51)Cr-release cytotoxicity assay. Tumor size and weight were significantly reduced by either weekly intramuscular injection of scIL-12/CMVpDNA or daily intraperitoneal injection of mrIL-12, and tumor became more localized. Survival was significantly increased when treatment started at 1 week of age as compared with that at the 6 weeks of age. Both NK and CD8(+) T cells were involved in the cytotoxicity against tumor cells and their antitumor activity was significantly reduced in tumor-bearing mice. TGF-β also inhibited the antitumor activity of NK and CD8(+) T cells. The immune suppression was completely reversed by IL-12 treatment and partially recovered by anti-TGF-β antibody. We conclude that both IL-12 gene therapy and recombinant protein therapy are effective against PTC. Given that the immune response is significantly suppressed in tumor-bearing mice and can be restored by IL-12, the current study raises a possibility of the application of IL-12 as an adjuvant therapy for thyroid cancer.

Background Pulmonary fibrosis is one of the leading indications for lung transplantation. The disease, which is of unknown aetiology, can be progressive, resulting in distortion of the extracellular matrix (ECM), inflammation, fibrosis and eventual death. Methods 13 patients born to consanguineous parents from two unrelated families presenting with interstitial lung disease were clinically investigated. Nine patients developed respiratory failure and subsequently died. Molecular genetic investigations were performed on patients' whole blood or archived tissues, and cell biological investigations were performed on patient-derived fibroblasts. Results The combination of a unique pattern of early-onset lung fibrosis (at 12–15 years old) with distinctive radiological findings, including 1) traction bronchiectasis, 2) intralobular septal thickening, 3) shrinkage of the secondary pulmonary lobules mainly around the bronchovascular bundles and 4) early type 2 respiratory failure (elevated blood carbon dioxide levels), represents a novel clinical subtype of familial pulmonary fibrosis. Molecular genetic investigation of families revealed a hypomorphic variant in S100A3 and a novel truncating mutation in S100A13, both segregating with the disease in an autosomal recessive manner. Family members that were either heterozygous carriers or wild-type normal for both variants were unaffected. Analysis of patient-derived fibroblasts demonstrated significantly reduced S100A3 and S100A13 expression. Further analysis demonstrated aberrant intracellular calcium homeostasis, mitochondrial dysregulation and differential expression of ECM components. Conclusion Our data demonstrate that digenic inheritance of mutations in S100A3 and S100A13 underlie the pathophysiology of pulmonary fibrosis associated with a significant reduction of both proteins, which suggests a calcium-dependent therapeutic approach for management of the disease., New evidence links an atypical form of pulmonary fibrosis with digenic mutations in the genes for calcium binding proteins S100A3 and S100A13. This implicates calcium homeostasis in the aetiology and pathogenesis of pulmonary fibrosis. http://bit.ly/2LyUQwb

Introduction ADCY5-related hyperkinesia encompasses a heterogeneous group of phenotypes, including paroxysmal chorea, myoclonus, and dystonia. The disease is attributed to mutations of ADCY5, which encodes an adenylate cyclase enzyme. The disease can occur in a sporadic or familial pattern. With exception of one study, all reports on familial ADCY5-related hyperkinesia were associated with an autosomal dominant inheritance. Herein, we describe a native Arabian Bedouin family with an autosomal recessive ADCY5-related disorder and expand the genotypic and phenotypic spectrum of this disorder. Methods The pedigree included 4 generations of a family with 6 affected individuals. The patients were examined clinically and radiologically. Homozygosity mapping and Whole Exome Sequencing (WES) were used to identify a variant, predicted to be pathogenic, which segregated with disease in this family. Results All patients presented with early-onset dystonia and myoclonus. The patients had delayed motor and language milestones, axial hypotonia, severe anxiety, social phobia, and isolation. One patient had dilated cardiomyopathy. WES of one affected individual revealed a novel homozygous missense mutation (c.1762G > A, p.D588N) of ADCY5, that segregated with disease in an autosomal recessive manner, and was absent in more than 1000 ethnically-matched chromosomes. The mutation replaces a highly conserved nucleotide and is predicted to be deleterious. Conclusion This study reports the second family with autosomal recessive childhood-onset ADCY5-related disorder and expands our understanding of phenotype/genotype correlations of this disorder.

Over 120 Type 2 diabetes (T2D) loci have been identified from genome-wide association studies (GWAS), mainly from Caucasian populations. Very limited knowledge is available on the Saudi Arabian population. In this study, 122 previously reported T2D-related variants from 84 loci were examined in a Saudi Arabian cohort of 1,578 individuals (659 T2D cases and 919 controls). Eleven single nucleotide polymorphisms (SNPs) corresponding to nine independent loci had a P value −4 ( = 0.05/122) were applied, none of the SNPs would have reached the significance level. Nine of the SNPs with a P value

Purpose Establishing links between Mendelian phenotypes and genes enables the proper interpretation of variants therein. Autozygome, a rich source of homozygous variants, has been successfully utilized for the high throughput identification of novel autosomal recessive disease genes. Here, we highlight the utility of the autozygome for the high throughput confirmation of previously published tentative links to diseases. Methods Autozygome and exome analysis of patients with suspected Mendelian phenotypes. All variants were classified according to the American College of Medical Genetics and Genomics guidelines. Results We highlight 30 published candidate genes (ACTL6B, ADAM22, AGTPBP1, APC, C12orf4, C3orf17 (NEPRO), CENPF, CNPY3, COL27A1, DMBX1, FUT8, GOLGA2, KIAA0556, LENG8, MCIDAS, MTMR9, MYH11, QRSL1, RUBCN, SLC25A42, SLC9A1, TBXT, TFG, THUMPD1, TRAF3IP2, UFC1, UFM1, WDR81, XRCC2, ZAK) in which we identified homozygous likely deleterious variants in patients with compatible phenotypes. We also identified homozygous likely deleterious variants in 18 published candidate genes (ABCA2, ARL6IP1, ATP8A2, CDK9, CNKSR1, DGAT1, DMXL2, GEMIN4, HCN2, HCRT, MYO9A, PARS2, PLOD3, PREPL, SCLT1, STX3, TXNRD2, WIPI2) although the associated phenotypes are sufficiently different from the original reports that they represent phenotypic expansion or potentially distinct allelic disorders. Conclusions Our results should facilitate the timely relabeling of these candidate disease genes in relevant databases to improve the yield of clinical genomic sequencing.