The p.Cys1281Tyr variant in the hinge module/flap region of thyroglobulin causes intracellular transport disorder and congenital hypothyroidism.

Congenital hypothyroidism (CH) due to thyroglobulin (TG) variants causes very low serum TG levels with normal or enlarged thyroid glands, depending on the severity of the defect, and with autosomal recessive inheritance. The purpose of this study was to functionally characterize p.Cys1281Tyr variant in the TG gene in order to increase our knowledge of the molecular mechanisms associated with CH. In order to find evidence that support the hypothesis that the p.Cys1281Tyr variant would affect the TG folding were performed amino acid prediction, 3D modeling and transient expression analysis in HEK293T cells. 18 of the 21″in silico" algorithms predict a deleterious effect of the p.Cys1281Tyr variant. The full-length 3D model p.Cys1281Tyr TG showed disulfide bond cleavage between the cysteines at positions 1249 and 1281 and rearrangement of the TG structure, while transient expression analysis indicated that p.Cys1281Tyr causes retention of the protein inside the cell. Consequently, these results show that this pathogenic variant makes it impossible for TG to fulfill its function in the biosynthesis process of thyroid hormones, causing CH. In conclusion, our results confirm the pathophysiological importance of misfolding of TG as a consequence of p.Cys1281Tyr variant located in the hinge module/flap region of TG.

Curating the gnomAD database: Report of novel variants in the thyroid peroxidase gene using in silico bioinformatics algorithms and a literature review.

Thyroid peroxidase (TPO) is a membrane-bound glycoprotein located at the apical side of the thyroid follicular cells that catalyzes both iodination and coupling of iodotyrosine residues within the thyroglobulin molecule, leading to the synthesis of thyroid hormone. Variants in TPO cause congenital hypothyroidism (CH) by iodide organification defect and are commonly inherited in an autosomal recessive fashion. In the present work, we report a detailed population analysis and bioinformatic prediction of the TPO variants indexed in the Genome Aggregation Database (gnomAD) v2.1.1. The proportion of missense cysteine variants and nonsense, frameshift, and splice acceptor/donor variants were analyzed in each ethnic group (European (Non-Finnish), European (Finnish), African/African Americans, Latino/Admixed American, East Asian, South Asian, Ashkenazi Jewish, Other). The results showed a clear predominance of frameshift variants in the East Asian (82%) and European (Finnish) (75%) population, whereas the splice site variants predominate in African/African Americans (99.46%), Other (96%), Latino/Admixed American (94%), South Asian (86%), European (Non-Finnish) (56%) and Ashkenazi Jewish (56%) populations. The analysis of the distribution of the variants indexed in gnomAD v2.1.1 database revealed that most missense variants identified in the An peroxidase domain map in exon 8, followed by exons 11, 7 and 9, and finally in descending order by exons 10, 6, 12 and 5. In total, 183 novel TPO variants were described (13 missense cysteine's variants, 158 missense variants involving the An peroxidase domain and 12 splicing acceptor or donor sites variants) which were not reported in the literature and that would have deleterious effects on prediction programs. In the gnomAD v2.1.1 population, the estimated prevalence of heterozygous carriers of the potentially damaging variants was 1:77. In conclusion, we provide an updated and curated reference source of new TPO variants for application in clinical diagnosis and genetic counseling. Also, this work contributes to elucidating the molecular basis of CH associated with TPO defects.

Mutational screening of the TPO and DUOX2 genes in Argentinian children with congenital hypothyroidism due to thyroid dyshormonogenesis.

PURPOSE: Primary congenital hypothyroidism (CH) is the most common endocrine disease in children and one of the preventable causes of both cognitive and motor deficits. We present a genetic and bioinformatics investigation of rational clinical design in 17 Argentine patients suspected of CH due to thyroid dyshormonogenesis (TDH). METHODS: Next-Generation Sequencing approach was used to identify variants in Thyroid Peroxidase (TPO) and Dual Oxidase 2 (DUOX2) genes. A custom panel targeting 7 genes associated with TDH [(TPO), Iodothyrosine Deiodinase I (IYD), Solute Carrier Family 26 Member 4 (SLC26A4), Thyroglobulin (TG), DUOX2, Dual Oxidase Maturation Factor 2 (DUOXA2), Solute Carrier Family 5 Member 5 (SLC5A5)] and 4 associated with thyroid dysembryogenesis [PAX8, FOXE1, NKX2-1, Thyroid Stimulating Hormone Receptor (TSHR)] has been designed. Additionally, bioinformatic analysis and structural modeling were carried out to predict the disease-causing potential variants. RESULTS: Four novel variants have been identified, two in TPO: c.2749-2 A > C and c.2752_2753delAG, [p.Ser918Cysfs*62] and two variants in DUOX2 gene: c.425 C > G [p.Pro142Arg] and c.2695delC [p.Gln899Serfs*21]. Eighteen identified TPO, DUOX2 and IYD variants were previously described. We identified potentially pahogenic biallelic variants in TPO and DUOX2 in 7 and 2 patients, respectively. We also detected a potentially pathogenic monoallelic variant in TPO and DUOX2 in 7 and 1 patients respectively. CONCLUSIONS: 22 variants have been identified associated with TDH. All described novel mutations occur in domains important for protein structure and function, predicting the TDH phenotype.

Curating the gnomAD database: Report of novel variants in the thyrogobulin gene using in silico bioinformatics algorithms.

Thyroglobulin (TG) is a large glycosylated protein of 2767 amino acids, secreted by the thyrocytes into the follicular lumen. It plays an essential role in the process of thyroid hormone synthesis. TG gene variants lead to permanent congenital hypothyroidism. In the present work, we report a detailed population and bioinformatic prediction analyses of the TG variants indexed in the Genome Aggregation Database (gnomAD). The results showed a clear predominance of nonsense variants in the European (Finnish), European (Non-Finnish) and Ashkenazi Jewish ethnic groups, whereas the splice site variants predominate in South Asian and African/African-American populations. In total, 282 novel TG variants were described (47 missense involving the wild-type cysteine residues, 177 missense located in the ChEL domain and 58 splice site variants) which were not reported in the literature and that would have deleterious effects in prediction programs. In the gnomAD population, the estimated prevalence of heterozygous carriers of the potentially damaging variants was 1:320. In conclusion, we provide an updated and curated reference source for the diagnosis of thyroid disease, mainly to congenital hypothyroidism due to TG deficiency. The identification and characterization of TG variants is undoubtedly a valuable approach to study the TG structure/function relations and an important tool for clinical diagnosis and genetic counseling.

A novel mutation in intron 11 donor splice site, responsible of a rare genotype in thyroglobulin gene by altering the pre-mRNA splincing process. Cell expression and bioinformatic analysis.

Thyroglobulin (TG) is a homodimeric glycoprotein synthesized by the thyroid gland. To date, two hundred twenty-seven variations of the TG gene have been identified in humans. Thyroid dyshormonogenesis due to TG gene mutations have an estimated incidence of approximately 1 in 100,000 newborns. The clinical spectrum ranges from euthyroid to mild or severe hypothyroidism. The purpose of the present study was to identify and characterize new variants in the TG gene. We report an Argentine patient with congenital hypothyroidism, enlarged thyroid gland and low levels of serum TG. Sequencing of DNA, expression of chimeric minigenes as well as bioinformatics analysis were performed. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a c.3001+5G > A, whereas the paternal mutation consists of p.Arg296*. Minigen analysis of the variant c.3001+5A performed in HeLa, CV1 and Hek293T cell lines, showed a total lack of transcript expression. So, in order to validate that the loss of expression was caused by such variation, site-directed mutagenesis was performed on the mutated clone, which previously had a pSPL3 vector change, to give rise to a wild-type clone c.3001+5G, endorsing that the mutation c.3001+5G > A is the cause of the total lack of expression. In conclusion, we demonstrate that the c.3001+5G > A mutation causes a rare genotype, altering the splicing of the pre-mRNA. This work contributes to elucidating the molecular bases of TG defects associated with congenital hypothyroidism and expands our knowledge in relation to the pathologic roles of the position 5 in the donor splice site.

p.L571P in the linker domain of rat thyroglobulin causes intracellular retention.

Thyroglobulin (TG), a large glycosylated protein secreted by thyrocytes into the thyroid follicular lumen, plays an essential role in thyroid hormone biosynthesis. Rattus norvegicus TG (rTG) is encoded by a large single copy gene, 186-kb long, located on chromosome 7 composed of 48 exons encoding a 8461-kb mRNA. Although the TG gene displays sequence variability, many missense mutations do not impose any adverse effect on the TG protein, whereas other nucleotide substitutions may affect its TG stability and/or TG intracellular trafficking. In order to gain a further understanding of the protein domains regulating its intracellular fate, we cloned a full-length cDNA from rTG into the pcDNA6/V5-His B expression vector. However, transient expression of the cDNA in HEK293T cells showed that the encoded protein was not a wild-type molecule, as it was unable to be secreted in the culture supernatant. Sequencing analyses revealed three random mutations, which accidentally emerged during the course of cloning: c.1712T>C [p.L571P] in the linker domain (amino acid positions 360 to 604), c.2027A>G [p.Q676R] in TG type 1-6 repeat and c.2720A>G [p.Q907R] in the TG type 1-7 repeat. Expression of cDNAs encoding a combination of two mutations [p.Q676R-p.Q907R], [p.L571P-p.Q907R] or [p.L571P-p.Q676R] indicated that any TG bearing the p.L571P substitution was trapped intracellularly. Indeed, we expressed the single point mutant p.L571P and confirmed that this point mutation was sufficient to cause intracellular retention of mutant TG in HEK293T cells. Endo H analysis showed that the p.L571P mutant is completely sensitive to the enzyme, whereas the will-type TG acquires full N-glycan modifications in Golgi apparatus. This data suggest that the p.L571P mutant contains the mannose-type N-glycan, that was added at the first stage of glycosylation. Complex-type N-glycan formation in the Golgi apparatus does not occur, consistent with defective endoplasmic reticulum exit of the mutant TG. Moreover, predictive analysis of the 3D linker domain showed that the p.L571P mutation would result in a significant protein conformational change. In conclusion, our studies identified a novel amino acid residue within the linker domain of TG associated with its conformational maturation and intracellular trafficking.