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

Intronic position +9 and -9 are potentially splicing sites boundary from intronic variants analysis of whole exome sequencing data

Whole exome sequencing (WES) can also detect some intronic variants, which may affect splicing and gene expression, but how to use these intronic variants, and the characteristics about them has not been reported. This study aims to reveal the characteristics of intronic variant in WES data, to further improve the clinical diagnostic value of WES. A total of 269 WES data was analyzed, 688,778 raw variants were called, among these 367,469 intronic variants were in intronic regions flanking exons which was upstream/downstream region of the exon (default is 200 bps). Contrary to expectation, the number of intronic variants with quality control (QC) passed was the lowest at the +2 and -2 positions but not at the +1 and -1 positions. The plausible explanation was that the former had the worst effect on trans-splicing, whereas the latter did not completely abolish splicing. And surprisingly, the number of intronic variants that passed QC was the highest at the +9 and -9 positions, indicating a potential splicing site boundary. The proportion of variants which could not pass QC filtering (false variants) in the intronic regions flanking exons generally accord with "S"-shaped curve. At +5 and -5 positions, the number of variants predicted damaging by software was most. This was also the position at which many pathogenic variants had been reported in recent years. Our study revealed the characteristics of intronic variant in WES data for the first time, we found the +9 and -9 positions might be a potentially splicing sites boundary and +5 and -5 positions were potentially important sites affecting splicing or gene expression, the +2 and -2 positions seem more important splicing site than +1 and -1 positions, and we found variants in intronic regions flanking exons over ± 50 bps may be unreliable. This result can help researchers find more useful variants and demonstrate that WES data is valuable for intronic variants analysis.

Perturbation of endoplasmic reticulum proteostasis triggers tissue injury in the thyroid gland

Defects in endoplasmic reticulum (ER) proteostasis have been linked to diseases in multiple organ systems. Here we examined the impact of perturbation of ER proteostasis in mice bearing thyrocyte-specific knockout of either HRD1 (to disable ER-associated protein degradation [ERAD]) or ATG7 (to disable autophagy) in the absence or presence of heterozygous expression of misfolded mutant thyroglobulin (the most highly expressed thyroid gene product, synthesized in the ER). Misfolding-inducing thyroglobulin mutations are common in humans but are said to yield only autosomal-recessive disease - perhaps because misfolded thyroglobulin protein might undergo disposal by ERAD or ER macroautophagy. We find that as single defects, neither ERAD, nor autophagy, nor heterozygous thyroglobulin misfolding altered circulating thyroxine levels, and neither defective ERAD nor defective autophagy caused any gross morphological change in an otherwise WT thyroid gland. However, heterozygous expression of misfolded thyroglobulin itself triggered significant ER stress and individual thyrocyte death while maintaining integrity of the surrounding thyroid epithelium. In this context, deficiency of ERAD (but not autophagy) resulted in patchy whole-follicle death with follicular collapse and degeneration, accompanied by infiltration of bone marrow-derived macrophages. Perturbation of thyrocyte ER proteostasis is thus a risk factor for both cell death and follicular demise.

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.

Thyroid hormone synthesis continues despite biallelic thyroglobulin mutation with cell death

Complete absence of thyroid hormone is incompatible with life in vertebrates. Thyroxine is synthesized within thyroid follicles upon iodination of thyroglobulin conveyed from the endoplasmic reticulum (ER), via the Golgi complex, to the extracellular follicular lumen. In congenital hypothyroidism from biallelic thyroglobulin mutation, thyroglobulin is misfolded and cannot advance from the ER, eliminating its secretion and triggering ER stress. Nevertheless, untreated patients somehow continue to synthesize sufficient thyroxine to yield measurable serum levels that sustain life. Here, we demonstrate that TG^{W2346R/W2346R} humans, TG^cog/cog mice, and TG^rdw/rdw rats exhibited no detectable ER export of thyroglobulin, accompanied by severe thyroidal ER stress and thyroid cell death. Nevertheless, thyroxine was synthesized, and brief treatment of TG^rdw/rdw rats with antithyroid drug was lethal to the animals. When untreated, remarkably, thyroxine was synthesized on the mutant thyroglobulin protein, delivered via dead thyrocytes that decompose within the follicle lumen, where they were iodinated and cannibalized by surrounding live thyrocytes. As the animals continued to grow goiters, circulating thyroxine increased. However, when TG^rdw/rdw rats age, they cannot sustain goiter growth that provided the dying cells needed for ongoing thyroxine synthesis, resulting in profound hypothyroidism. These results establish a disease mechanism wherein dead thyrocytes support organismal survival.

Structure and genetic variants of thyroglobulin: Pathophysiological implications

Thyroglobulin (TG) plays a main role in the biosynthesis of thyroid hormones (TH), and, thus, it is involved in a wide range of vital functions throughout the life cycle of all vertebrates. Deficiency of TH production due to TG genetic variants causes congenital hypothyroidism (CH), with devastating consequences such as intellectual disability and impaired growth if untreated. To this day, 229 variations in the human TG gene have been identified while the 3D structure of TG has recently appeared. Although TG deficiency is thought to be of autosomal recessive inheritance, the introduction of massive sequencing platforms led to the identification of a variety of monoallelic TG variants (combined with mutations in other thyroid gene products) opening new questions regarding the possibility of oligogenic inheritance of the disease. In this review we discuss remarkable advances in the understanding of the TG architecture and the pathophysiology of CH associated with TG defects, providing new insights for the management of congenital disorders as well as counseling benefits for families with a history of TG abnormalities. Moreover, we summarize relevant aspects of TH synthesis within TG and offer an updated analysis of animal and cellular models of TG deficiency for pathophysiological studies of thyroid dyshormonogenesis while highlighting perspectives for new investigations. All in all, even though there has been sustained progress in understanding the role of TG in thyroid pathophysiology during the past 50 years, functional characterization of TG variants remains an important area of study for future advancement in the field.