Introduction

Rare monogenic disorders are often diagnosed late, primarily due to missing medical awareness and knowledge. Morbidity and mortality are commonly worsened by a protracted diagnosis. Although causal therapies are not available for most single-gene disorders, prevention of complications and avoidance of ineffective and/or contraindicated medical interventions can dramatically improve individual outcomes.1

Despite tremendous progress in genomic medicine over the last decade, the underlying genetic defect has only been identified for approximately one-third to one-half of the 7000–15 000 presumed monogenic diseases.2, 3 Purported reasons for the large proportion of unresolved monogenic diseases include complicated genetic architecture, technical limitations, variants of unclear pathogenicity and/or occurring de novo, and incomplete or incorrect phenotyping. An additional explanation that has so far received little attention, is that a proportion of the unresolved monogenic diseases in fact may represent ‘blended phenotypes’4, 5 that arise from the presence of molecular causes of more than one monogenic disease in an individual. This may be especially relevant for offspring from consanguineous matings.6 Here, we describe the first case of a patient in whom a combination of three recessive single-gene disorders gives rise to a clinical picture of congenital hypothyroidism, hypomagnesemia and hypercholesterolemia.

Materials and methods

Patients and materials

Patients were recruited by physician-initiated referral. The study was conducted in accordance with the Declaration of Helsinki protocols and approved by the institutional ethics review board of Freiburg University Hospital, Germany. Written informed consent for molecular studies was obtained from the affected individuals and/or their legal guardians in accordance with current German law (GenDG). Control samples and primary cells were collected from ancestry-, sex- and age-matched healthy individuals under the same criteria and regulations. DNA and RNA were obtained from peripheral-blood leukocytes by standard extraction procedures.7

Whole-exome sequencing data analysis

Exome sequences were enriched from DNA using the Agilent SureSelect Exome Enrichment kit V4 (Agilent Technologies, Waldbronn, Germany). Sequencing was performed on the 5500xl SOLiD System (Life Technologies, Darmstadt, Germany); 75 bp single-end reads were produced. Sequence alignment and variant calling were performed using the targeted re-sequencing workflow in LifeScope v2.5.1 (Life Technologies). The reads were aligned to the human hg19 reference genome. Variant calling stringency was set to high. Called variants were annotated and filtered using GeneTalk (http://www.gene-talk.de/).8 For the purpose of variant filtering, HomozygosityMapper (http://www.homozygositymapper.org/)9 was used to detect regions of homozygosity shared by the two patients but not their mother. In addition, regions of homozygosity were detected individually in the two patients using HomozygosityMapper. When filtering for autosomal-recessive candidate variants, we retained stop, splice, frameshift and non-synonymous variants with prevalence <1% in the 1000 Genomes Project and Exome Sequencing Project. Filtered candidate variants mapping into regions of homozygosity were prioritized, but all other filtered candidate variants were also examined for potentially causal variants. Details on additional bioinformatics analyses, including those of publicly available RNA-seq data from thyroid tissue to characterize the effects of the variant underlying hypothyroidism, are available in the Supplementary Material.

Variant confirmation and molecular studies

For the genomic Sanger and cDNA analysis, exons including intron–exon boundaries of TRPM6 (NM_017662.4; exon numbering according to NG_017036.1; variant affects splice donor site adjacent to exon 20), ABCG5 (NM_022436.2; exon numbering according to NG_008883.1; variant affects exon 10), and TG (NM_003235.4; exon numbering according to NG_015832.1; variant affects splice donor site adjacent to exon 5) were PCR-amplified using standard protocols as described previously;7 details are given in the Supplementary Material. The three variants and phenotypes have been submitted to ClinVar database (accession numbers SCV000257394 (TRPM6), SCV000257395 (ABCG5), SCV257396 (TG)).

Measurement of blood sterols

Cholesterol, its metabolite 5α-cholestanol, the plant sterols campesterol, stigmasterol and sitosterol as well as their 5α-saturated metabolites 5α-campestanol and 5α-sitostanol were determined by gas chromatography-flame ionization detection using 5a-cholestane as internal standard.10 Structure identity was verified by gas chromatography-mass spectrometry in Scan Modus comparing the spectra with those from authentic highly pure compounds.

Results

The index patients are a 19-year-old boy (individual II-01) and his 14-year-old sister (individual II-02) from a family of Turkish origin (compare with pedigree in Figure 2). Their parents are first-degree cousins; both paternal and maternal grandparents were also reportedly related. Both siblings were diagnosed with congenital primary hypothyroidism by TSH newborn screening within the first 3 days of life and given thyroid hormone daily since. In postnatal examination and follow-up, hypothyroidism was classified as a non-goiterous thyroid developmental disorder of unknown etiology. Laboratory parameters and sonography on regular follow-up visits indicated thyroid dysgenesis and secretory thyroid hormone insufficiency with an orthotopic, but small and dysplastic thyroid (Figure 1a,Supplementary Figure S1A). Normal audiometry in both children made Pendred Syndrome (MIM #274600) unlikely; no further investigations to classify primary hypothyroidism were undertaken. In the following years, the children had no developmental problems under adapted thyroxine substitution.

Figure 1
figure 1

(a) Representative transverse thyroid sonogram of the girl (individual II-02) at the age of 13 years. Volume of both lobes was reduced and echogenicity increased, but no cystic lesions were found. Measurement of the thyroid gland was performed using a Philips ultrasound machine with a 7.5 MHz transducer. For each thyroid lobe, the maximum perpendicular anteroposterior (depth) and mediolateral (width) diameters were measured on a transverse image of the largest diameter. The greatest craniocaudal (length) diameter of each lobe was measured on a longitudinal image. (b) Time course of serum magnesium levels in both index patients. Dashed lines/squares represent individual II-01 and solid lines/circles individual II-02. Arrows indicate the time range when repeated seizures occurred in the girl but do not represent the number of episodes. Reference range for serum magnesium: 0.75–1.05 mmol/l (gray area). (c) Gas chromatographic (GC) separation combined with flame ionization detection (FID) of the trimethylsilyl ethers of cholesterol, 5α-cholestanol, campesterol, and sitosterol from a serum sample of the patient (PII-02) shortly before the start of ezetimibe treatment, in comparison with extracted sterols from an healthy, age- and gender-matched Control subject. (d) Treatment with 10 mg ezetemibe for 13 months significantly reduced cholesterol as well as plant sterols (values based on one measurement at each time point).

However, at the age of 2 years, the girl was brought to the emergency room with generalized seizures. Ictal EEG showed generalized spike or polyspikes-and-wave complexes (Supplementary Figure S1B). The fit responded well to emergent benzodiazepine treatment. Later, clinical presentation and EEG changes were consistent with a diagnosis of myoclonic-astatic epilepsy (MAE). MAE, also known as Doose Syndrome (MIM #607208), is an epilepsy syndrome of early childhood usually characterized by difficult to control generalized seizures. Laboratory results at that time were unremarkable and she was started on anticonvulsant levetiracetam. Although no additional generalized seizure occurred over several months, she continued to suffer from several myoclonic events per week. On a subsequent blood draw, severe hypomagnesemia with <0.2 mmol/l was unraveled (Figure 1b, Supplementary Table S1). Serum magnesium levels had not been checked in any previous investigation. The other electrolytes, pH, and bicarbonate levels had always been normal. Anticonvulsants were discontinued and switched to oral magnesium substitution therapy, which proved effective to cure MAE instantly. Furthermore, EEG changes completely normalized even though blood magnesium levels remained subnormal (Figure 1c). Testing in the brother revealed low magnesium levels as well, although he had never been symptomatic (Figure 1c, Table 1). No changes were found in his EEG and ECG investigations. Electrolytes were normal in both parents. Further work-up of hypomagnesemia excluded Gitelman Syndrome (MIM #263800), and both siblings received a tentative diagnosis of defective intestinal resorption/loss of magnesium.

Table 1 Laboratory results of the two patients

On a routine blood draw, hypercholesterolemia was diagnosed in the girl (Table 1). As concomitant TSH was suppressed and fT3 and fT4 values were normal, hypercholesterolemia was not secondary to hypothyroidism. In spite of a low-normal BMI and no dietary particularities, both her LDL-cholesterol and total cholesterol levels were increased, while triglycerides were normal (type IIA, Table 1). No clinical signs of systemic cholesterol deposition evolved, and evaluation of intima-media thickness by high-resolution B-mode echo-Doppler sonography demonstrated normal findings for her age at several occasions. An attempt to treat hyperlipidemia by cholestyramine at the age of 9 years lowered cholesterol levels initially but high cholesterol levels were observed on several occasions thereafter. Notwithstanding, under thyroid hormone replacement therapy and oral substitution of magnesium, both patients developed well and attended school without problems.

Given that both siblings from this consanguineous marriage were affected, we suspected a potential genetic origin. As no known monogenic disorder matched our patients’ multi-faceted phenotype, we assumed an underlying recessive disease of yet unknown origin. Therefore, whole-exome sequencing (WES) of the two index patients, their parents and paternal grandparents was undertaken to search for causal variants.

High-quality WES data were obtained from the two patients, their mother, and the paternal grandparents. Detailed summary statistics are provided in Supplementary Table S1. Subsequently, variant annotation and filtering gave rise to a short list of 73 candidate variants. Although the attempts to pinpoint a single variant that could explain the observed clinical phenotype were unsuccessful, we identified separate variants in three genes known to cause single-gene disorders that collectively explained the clinical phenotype.

A splice variant in TRPM6 causes hypomagnesemia

Regions of homozygosity comprised 5.7% of the boy’s genome (166437.6 kb, hg19) and 7.2% of the girl’s genome (207969.7 kb). Focusing on the 95476.1 kb regions of homozygosity shared by the children but not their mother (Supplementary Figure S2), we identified the variant c.2667+1G>A in the donor splice site adjacent to exon 20 of TRPM6, for which both patients were homozygous (Supplementary Table S2). TRPM6 encodes a magnesium transport protein expressed predominantly in intestinal and renal epithelia, where it mediates magnesium absorption. Mutations in TRPM6 are known to cause autosomal-recessive hypomagnesemia with secondary hypocalcemia (MIM #602014).11 Indeed, the identified variant c.2667+1G>A has previously been reported in several patients of Turkish origin.12 Sanger sequencing in our patients and other family members confirmed heterozygosity of both parents and co-segregation of variant and disease (Figure 2). Hypocalcemia was not regularly observed in either of our patients, but secondary hypocalcemia in patients with TRPM6 mutations is reportedly variable.12

Figure 2
figure 2

Pedigree of the nuclear family (middle panel), coloring of the filled symbol segments indicates manifestation of phenotypes (green: hypomagnesemia; orange: sitosterolemia; blue: thyroid dyshormonogenesis). Confirmatory genomic Sanger sequencing of TRPM6, ABCG5, and TG mutations in the boy (individual II-01, left panel) and the girl (individual II-02, right panel).

A stop variant in ABCG5 causes sitosterolemia

Among the remaining candidate variants, we identified the c. 1336C>T/p.(Arg446*) variant, causing a premature stop codon in the ABCG5 gene (Supplementary Table S2). ABCG5 encodes an ATP-binding cassette (ABC) half-transporter expressed in liver, colon, and intestine, which limits intestinal absorption and promotes biliary excretion of sterols. ABCG5 mutations have been identified as a cause of autosomal-recessive sitosterolemia (MIM #210250).13 The p.(Arg446*) variant has been reported previously in a consanguineous family with sitosterolemia from Iran.13 Although clinical manifestations of sitosterolemia can be similar to those of familial hypercholesterolemia, the disease can be diagnosed by elevated serum concentrations of the plant sterols, sitosterol and campesterol. We, therefore, measured these sterols in the girl and found them significantly increased, confirming the diagnosis of sitosterolemia (Figure 1c, Table 1). Increased total cholesterol and slightly increased plant sterol concentrations were detected in the heterozygous father. The variant was confirmed by Sanger sequencing and segregated with the phenotype (Figure 2). The girl was started on 10 mg ezetimibe, which after 13 months led to a significant reduction of cholesterol, 5α-cholestanol, plant sterols and plant stanols (Figure 1d).

Hypothyroidism caused by a splice variant in TG

After loosening initial filtering criteria for exon/intron junctions from 2 to 5 bp, we identified a variant at a highly conserved +5 splice donor site in intron 5 (IVS5) of the TG gene (c.638+5G>A), encoding thyroglobulin (Supplementary Table S2). The variant was predicted to reduce splicing by a number of splice site prediction programs (Supplementary Figure S3A, Supplementary Table S3). Using a recently published method,14 which predicts how strongly a genetic variant affects RNA splicing with very high accuracy, the variant was predicted to lead to a 36.8% reduction of transcripts containing exon 5 in comparison with the wild-type sequence. For comparison, the corresponding estimate for the known splice variant c.2667+1G>A in TRPM6 was 12.5%. Recessive TG mutations are a known cause of congenital thyroid dyshormonogenesis (MIM #274700).15 TG is a large gene producing many alternatively spliced transcripts.16 Mutations affecting splice sites make up a significant portion of all reported pathogenic TG variants. In particular, missplicing of TG pre-mRNA due to another variant in the consensus splice donor site of IVS5 is known to induce congenital goiter, hypothyroidism and subsequent malignant transformation.17

We used RT-PCR analysis of patient-derived PBMCs to directly assess the consequences of the novel c.638+5G>A on TG mRNA. As predicted, the +5 variant resulted in a transcript missing exon 5 (Figures 3a and b, and Supplementary figure S3B). To exclude that this TG mRNA was a normal transcript splice variant, we performed in silico bioinformatics analyses mapping available RNA-seq data sets to TG exons (Materials and methods section, Supplementary Figure S3C). No natural TG mRNA missing exon 5 was detected either in fetal and adult thyroid tissue or in white blood cells. Different from the known IVS5+1 variant,17 semi-quantitative RT-PCR analysis of the intron IVS5+5 variant indicated that – while overall TG mRNA levels were reduced (Figure 3a, Supplementary Figure S3B) – not all transcripts were mis-spliced in our patients, as evidenced by a faint band of the size of wild-type TG. In keeping with a less penetrant effect of the IVS5+5 variant compared with the known IVS5+1 variant, thyroglobulin in patients’ sera was detectable by indirect ELISA (Table 1) and represents a plausible explanation of the mild phenotype of non-goiterous thyroid dyshormonogenesis.

Figure 3
figure 3

(a) PCR detection of exon 5 skipping in TG transcripts from the girl and healthy controls; cDNA from different leukocyte fractions was amplified with primers annealing in exon 4 (forward) and exon 6 (reverse, blue arrows in b). Expected product sizes are 317 bp in wild-type TG (blue arrow in a) and 157 bp when exon 5 is spliced out (red arrow in a). Note a faint signal at 317 bp in patient 1, indicating residual normal splicing of TG. (b) Sanger sequencing of exon-spanning cDNA amplicons shown in b. The aberrant splicing results in a frameshift at amino acid 161 (underlined). In the TG reading frame, skipping of exon 5 gives rise to a premature stop codon five codons downstream within exon 6. The predicted protein from this transcript contains 165 residues, including 161 residues of normal TG sequence. CTL, control; GRAN, granulocytes; MNC, monocytes. A full color version of this figure is available at the European Journal of Human Genetics journal online.

Together, all molecular data confirm that the clinical presentation of this family can be explained by three separate known single-gene disorders.

Discussion

Here, we report the first instance of a family in whom three separate single-gene disorders, familial hypomagnesemia, sitosterolemia, and thyroid dyshormonogenesis, give rise to a complex and previously unobserved clinical picture. Against our expectations, the molecular cause in offspring from a consanguineous marriage did not reside in a single gene transmitted in an autosomal-recessive fashion. Instead, the blended phenotype resulted from three separate molecular causes, and could only be resolved by their combined identification. These findings have several implications.

First, our findings strengthen the evidence for blended phenotypes4, 5 arising from more than one underlying monogenic disease. In our cases, these diseases were independent and the causative variants are located on different chromosomes, which distinguishes this case from blended phenotypes arising from contiguous gene deletion syndromes. Studies of clinical WES and whole-genome sequencing have previously reported on patients with phenotypes arising from two, but not three, independent single-gene disorders.4, 5, 18, 19, 20 We propose that a sizeable proportion of the presumed monogenic diseases of yet unknown origin may in fact be due to clustering of several single-gene disorders. One previous large-scale clinical WES study reported that 5% of patients with a molecular diagnosis had two underlying single-gene disorders.5 Another study in patients from consanguineous families reported a slightly higher rate (one out of 18),20 which may be an underestimate given the identified variants only partially explained the clinical phenotype in two families. It needs to be noted that both estimates were influenced by indication for WES as well as previous knowledge about the identified variants.

In our patients, mutation identification was facilitated by homozygosity mapping as well as previous knowledge about TRPM6, ABCG5, and TG. This task will be more challenging in the setting of de novo mutations, variants of unknown significance or incomplete penetrance, or when each of the causal variants follows a different mode of inheritance. The delineation of a blended phenotype will further be complicated when variants map into genes of yet unknown function. Moreover, even well-characterized variants may provoke unpredictable consequences in a complex genetic background. Variants may, for instance, interact because they directly or indirectly act on a common molecular pathway. In addition, disease genes may have similar functions during embryonic development, and converging developmental effects may also modify the clinical phenotype. Discrete features may thus be ameliorated, aggravated, or completely transformed to a novel phenotype, which – together with epigenetic factors – should contribute to the variable genotype–phenotype relationship commonly observed for single-gene disorders. Together, these issues complicate the estimation of the proportion of presumed monogenic diseases of yet unknown origin that in truth arise from a combination of several monogenic diseases.

Second, our study is the first report of three concomitant independent monogenic diseases in one person, pressing the question for the likelihood of this observation. The prevalence of each of these diseases is exceedingly low, ranging from an estimated 1:100 000 to 1:1 000 000. Therefore, in unrelated individuals and assuming independence of the harmful alleles, the chance of all three of them co-occurring is less than 1 in a quadrillion (1015). In our case, the patients are offspring from a consanguineous mating and the probability of co-occurrence of the diseases are therefore higher. The number of recessive single-gene diseases observed in offspring of such matings should be influenced by several factors: the number of variants for which the parents are heterozygous, the degree of parental relatedness and reproductive fitness as well as the penetrance and lethality of the resulting recessive diseases. The clustering of single-gene disorders will be facilitated if the underlying variants are hypomorphic or variably penetrant as was the case for TG variant we observed. Yet, the probability of observing three concomitant monogenic recessive diseases in one patient should still be very low, given the probability that a randomly sampled allele at a specific locus is identical-by-descent for first-degree cousins is 6.25% and considering the factors listed above.

Third, our findings have both practical and theoretical implications for clinical practice. Although both patients continue to receive thyroid hormone replacement therapy and oral magnesium, the lipid-lowering treatment for the girl was switched from cholestyramine to ezetimibe, a specific therapy for sitosterolemia, and put on a diet low in vegetable fat. Ezetimibe blocks the uptake of sterols in the intestine and thereby lowers cholesterol and plant sterol levels. The identification of the underlying ABCG5 variant now allows for monitoring treatment efficacy by measuring blood sterol levels. In addition, the patient will be closely followed for the potential development of cardiovascular disease, although reportedly not all patients with phytosterolemia develop cardiovascular disease.21

Advanced diagnostic procedures to investigate exceedingly rare metabolic disorders such as those of sterol absorption are readily available, but their utilization requires a concrete clinical suspicion. As effective prevention and therapy depend on correct classification of pleiotropic phenotypes, circumstantial evidence for the extraordinary is critical. WES has greatly facilitated and accelerated diagnosis in rare disorders. Its cost-effectiveness and utility in clinical diagnosis are increasingly recognized.22, 23 Even without clinical clues and a working hypothesis, unbiased comprehensive genetic analysis may deliver pathogenic variants. Promoted by success and falling costs, however, widespread use of clinical WES may interfere with clinical judgment. Although WES was instrumental in our case as well, it exemplifies a so far un(der)-recognized pitfall: the successful identification of causative variants in one gene may not be sufficient to resolve a clinical phenotype completely, necessitating the search beyond variants that can only account for part of the clinical picture. As the prevalence of blended phenotypes4, 5 from multiple monogenic disorders is currently unknown, systematic analysis and regular re-analysis of larger WES data sets and, if applicable, integration with results from complementary diagnostic procedures such as karyotyping, is needed to gauge the significance of such ‘complex-simplex’ diseases in daily practice.