- Open Access
Compound heterozygous variants in DYNC2H1 in a foetus with type III short rib-polydactyly syndrome and situs inversus totalis
BMC Medical Genomics volume 15, Article number: 55 (2022)
Short-rib thoracic dysplasia 3 with or without polydactyly (SRTD3, OMIM: 613091) is an autosomal recessive disorder. SRTD3 presents clinically with a narrow thorax, short ribs, shortened tubular bones, and acetabular roof abnormalities. Clinical signs of SRTD3 vary among individuals. Pathogenic variants of DYNC2H1 (OMIM: 603297) have been reported to cause SRTD3.
We performed a detailed clinical prenatal sonographic characterization of a foetus with SRTD3. Trio whole-exome sequencing was used to identify causative variants in the family. The identified variants in the families were validated by Sanger sequencing and mass spectrometry. Multiple computational tools were used to predict the harmfulness of the two variants. A minigene splicing assay was carried out to evaluate the impact of the splice-site variant.
We evaluated prenatal sonographic images of the foetus with SRTD3, including abnormal rib curvature, narrow thorax, bilateral hypoplastic lungs, bilateral polydactyly, syndactyly, and foetal visceral situs inversus with mirror-image dextrocardia. We revealed novel compound variants of DYNC2H1 (NM_001377.3:c.11483T > G (p.Ile3828Arg) and c.2106 + 3A > T). Various statistical methods predicted that the variants would cause harmful effects on genes or gene products. The minigene assay findings suggested that c.2106 + 3A > T caused the skipping over exon 14, producing an exon 14 loss in the protein.
This study identified a foetus with SRTD3 with situs inversus totalis with mirror-image dextrocardia in a Chinese family, revealing two novel compound heterozygous dynein cytoplasmic 2 heavy chain 1 (DYNC2H1) variants, expanding the phenotypic spectrum of SRTD3. The minigene study of c.2106 + 3A > T was predicted to cause an inframe exclusion of exon 14, which was predicted to have important molecular functions. Our findings strongly supported the use of WES in prenatal diagnosis and helped to understand the correlation of genotype and phenotypes of DYNC2H1. The specific sonographic findings and the molecular diagnosis helped add experience to further our expertise in prenatal counselling for SRTD3.
Short-rib thoracic dysplasia 3 with or without polydactyly (SRTD3, OMIM: 613091) is an autosomal recessive disorder. SRTD3 presents clinically with a narrow thorax, short ribs, shortened tubular bones, and acetabular roof abnormalities . Clinical signs of SRTD3 vary among individuals. Nonskeletal manifestations also include cleft lip and palate and some major organs, such as the brain, eyes, heart, kidneys, liver, pancreas, intestinal abnormalities, and genitalia . Some foetuses die in the neonatal period due to severe restriction of the thorax secondary to respiratory insufficiency, while others survive. There are many SRTD subtypes. The four lethal SRTDs are SRTD types I-IV. SRTD types I-IV are mainly characterized by short limb deformities, thoracic stenosis, often combined with polydactyly and multiple anomalies of other organs. Each SRTD subtype has its typical features and a common overlapping phenotype . The prenatal phenotype differs and has variable expressivity. Therefore, it is a challenge to provide a timely and accurate diagnosis . Sonographic examinations provide clues, and genetic tests help to give a molecular diagnosis.
Patient and sonographic examinations
A 22-year-old woman was referred to our department for routine sonographic examination at 22.4 gestational weeks (GW). Systematic prenatal sonographic examinations, including two- and three-dimensional ultrasound and echocardiography, were obtained to a final diagnosis by at least two chief physicians. GE Volusion E10 was used to perform the sonographic examination with a 1–5 MHz probe. These sonographic scans were all performed using a GE Volusion E10 machine with a 1–5 MHz sector transducer.
Trio whole-exome sequencing (WES) of the family was performed at BGI company, Wuhan, China. In this method, genomic DNA was collected from blood from the parents and umbilical cord blood from the foetus. First, DNA was broken into 100- to 500-bp fragments. The 200- to 300-bp fragments were collected by magnetic beads (VAHTS™ DNA Clean Beads, Vazyme). Base “A” was added at the 3' overhangs, which could pair the “T” base with a special adapter. After elution and postcapture amplification, the library was enriched for 16 h at 65 °C by array hybridization (KAPA Hyper Exome, Roche). The magnitude of enrichment of the product was estimated. The qualified products were pooled, quantified, and prepared for circularization. Finally, the MGISEQ-2000 sequencing platform was used for variant detection. The average sequencing depth of the target region was ≥ 180X, and the percentage of loci with an average depth of > 20X in the target region was > 95%. After the sequencing was completed, raw reads from the machine were evaluated for sequencing quality, and low-quality and joint-contaminated reads were removed for further analysis .
Sequenced fragments were aligned with the UCSC hg19 human reference genome by BWA (Burrows Wheeler Aligner) to remove duplicates . Base quality score recalibration, SNVs (single-nucleotide variants), InDels (insertions and deletions), and genotypes were detected using GATK software (GATK version number 4.1.9) . Variation annotation and screening were conducted based on the clinical information of subjects, population databases, disease databases, and biological information prediction tools, including NCBI dbSNP, HapMap, the 1000Genomes human dataset, and a database of 100 Chinese healthy adults. ExomeDepth was used for copy number variation (CNV) detection at the exon level. The variants were named according to the Human Genome Variation Society (HGVS) nomenclature. To predict the CNV, we used the in-house database (https://phoenix.bgi.com/autocnv/). To predict the effects of missense variants, we used dbNSFP, which contains seven well-established in silico prediction programs (Scale-Invariant Feature Transform (SIFT), Polyphen2, LRT, Mutation Taster, and PhyloP) . To predict the effect of splice variants, we used prediction tools (SpliceAI, dbscSNV_ADA and dbscSNV_RF). To predict the conservation of variants, we used tools such as PhyloP Vertebrates, PhyloP Placental Mammals and GERP++_RS. Variation pathogenicity was classified according to the American College of Medical Genetics and Genomics (ACMG), the Association for Molecular Pathology (AMP) Sequence Variation Interpretation Guidelines, the Clingen Sequence Variation Interpretation Working Group and the Association for Clinical Genomic Science (ACGS), and other detailed interpretations of these guidelines .
To verify the results of gene chip capture and high throughput sequencing, the two suspected variants were validated by Sanger sequence and mass spectrometry. The DYNC2H1 variant c.2106 + 3A > T (NM_001377.3) was validated by Sanger sequencing. Primers were designed upstream and downstream of the segment. Polymerase chain reaction (PCR) amplification and Sanger sequencing were performed using the forward primer TTTGCCAGACTGCCTTATCA and the reverse primer AGACTGAGGAGGGGAGGCTA. The DYNC2H1 c.11483T > G (NM_001377.3) validation was performed by using the MassARRAY system (Agena Bioscience company) according to a previously described method . The MassARRAY System uses matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectrometry for the precise detection of DNA molecules. Genetic variants are distinguished by analysis of their mass. The probe sequence was CTACAGTCAAGTCTGAAGA.
Minigene splicing assay
A genomic segment encompassing the variant c.2106 + 3A > T sequence along with flanking intronic sequences, including exon 13, part of intron 13–14, exon 14, intron 14–15 and exon 15, was PCR-amplified from patient genomic DNA. The segments containing the wild-type (WT, c.2106 + 3A) and mutant-type (MT, c.2106 + 3T) sequences were separately cloned into the minigene vector pMini-CopGFP (Hitrobio.tech, China). The following primers were used to construct the vectors using the ClonExpress® II One Step Cloning Kit (Vazyme Biotech Co., Ltd): DYNC2H1-AF, AAGCTTGGTACCGAGCTCGGATCCGTGGCACATTTTTATAATTCTATTGATC; DYNC2H1-AR, AAGAAGGCAGACAGCCATGGGGTTTCACCATGTTGG; DYNC2H1-BF, ATGGCTGTCTGCCTTCTTCATTCTTCCTTTCCG; DYNC2H1-BR, TTAAACGGGCCCTCTAGACTCGAGCTGTGCTTCTACAGTTGCTAAGCCAGTT. The correct wild-type and mutant Minigene plasmids were transfected into HEK293T cell line for 48 h. Total RNA was extracted, and cDNA synthesis was performed (R212-01, Vazyme Biotech Co., Ltd.). Real-time PCR was performed using the primers F-GGCTAACTAGAGAACCCACTGCTTA and R-CTGTGCTTCTACAGTTGCTAAGC, and the products were sequenced.
Patient and sonographic findings
A 22-year-old woman was referred to our department for routine sonographic examination at an ultrasonic gestational age of 20.6 weeks and 22.4 weeks from the last menstrual period (LMP). This was the second pregnancy of the woman. She had an early pregnancy loss in her first pregnancy. The couple was healthy and non-consanguineous. The sonographic examination detected multiple anomalies. The measurement value (− 2 standard deviation (SD)/+ 2 SD) of the biparietal diameter (BPD) was 5.7 cm (4.9 cm/6.1 cm); the head circumference (HC) was 21.7 cm (18.3 cm/21.8 cm); the abdominal circumference (AC) was 16.7 cm (15.7 cm/19.6 cm), the femur length (FL) was 3.3 cm (3.3 cm/4.1 cm); and the humerus length (HL) was 2.5 cm (3.21 cm/4.05 cm), respectively. (Fig. 1) Both the HL and FL values were below – 2 SD. The foetus showed abnormal rib curvature, a narrow thorax, bilateral hypoplastic lungs, bilateral polydactyly, a hallux valgus angle on the foot and syndactyly. The foetus was situs inversus totalis with mirror-image dextrocardia. The liver and gallbladder inferior vena cava was on the left side of the body, the stomach and aorta were on the right, and the spleen was on the right, but the spleen was not easy to visualize in the foetus. The heart presented with mirror-image dextrocardia. The bilateral humerus was curved (Fig. 1).
An equal echo of approximately 0.4*0.3 cm in size was seen on the lateral aspect of the little fingers on both sides. A strong osseous echo could be seen inside. The ribs were short and curved. The thoracoabdominal junction was markedly depressed, and the thorax was markedly reduced (3.4/5.18 = 0.66). The thoracoabdominal transverse diameter ratio was 0.66. The size of the left lung was 2.7 * 1.1 cm, and that of the right lung was 1.6 * 1.0 cm. The two foetal kidney collecting systems were separated by approximately 0.4 cm. Echoes of 6 toes were seen on the left foot. The first and second toes of the right foot were widely spaced, and the fourth and fifth toes appeared to have syndactyly. The tibial length was 2.7 cm (− 2.8 SD); the ulnar length was 2.6 cm (− 3.2 SD). The final US diagnosis was short-rib thoracic dysplasia (narrow thorax, broken ribs with curvature, short limb deformities, polydactyly, hallux valgus angle), situs inversus totalis, and mirror-image dextrocardia. The parents chose to terminate the pregnancy due to multiple anomalies.
Two heterozygous variants of the DYNC2H1 gene associated with the subject's phenotype in SRTD3 were detected, NM_001377.3: c.11483T > G (p.Ile3828Arg) and c.2106 + 3A > T. We did not detect rare or likely pathogenic homozygous or heterozygous variants in other STRD-causing genes in family trio WES. No pathogenic chromosomal CNV variants above 1 M and LOH variants above 5 M associated with the phenotypes of the proband were detected through the analysis of the high-throughput data. c.11483T > G (p.Ile3828Arg) was detected in the father and the foetus in family trio WES. The c.11483T > G variant was located in exon 78 of the DYNC2H1 gene and caused the alteration of isoleucine to arginine at amino acid position 3828. c.11483T > G had frequencies of 0.000013 in GnomAD and 0.00019 in GnomAD-EAS (East Asian population). This variant was not found in the control population in the ESP, 1000, and EXAC databases. Various statistical methods, including conservative and evolutionary predictions, predicted that the variant would cause harmful effects on genes or gene products. Computer software tools, including SIFT, Mutation Taster, and Condel, predicted that predicted the missense variant was pathogenic. All three tools (PhyloP Vertebrates, PhyloP Placental Mammals, GERP++) showed that c.11483T > G was conserved (Table 1). c.2106 + 3A > T was detected in the mother and the foetus in family trio WES. No frequency was found in the ESP, 1000, EXAC, and GnomAD databases. All three tools predicting splicing site effects (SpliceAI, dbscSNV_RF, dbscSNV_ADA) raised the possibility that c.2106 + 3A > T caused loss of function for authentic splice sites (Table 2). The SpliceAI showed that the probability that position 11:103,004,436 (= 103,004,439–3) was used as a splice donor was decreased by 0.72. Mass spectrometry confirmed that the foetus carried c.11483 > G, which was inherited from the father. Sanger sequencing confirmed that the foetus carried c.2106 + 3A > T, which was inherited from the mother (Fig. 2).
Splicing study of the DYNC2H1 c.2106 + 3A > T variant
A minigene expression assay revealed the synthesis of a messenger RNA lacking exon 14, which generated a predicted in-frame deleted protein. The agarose gel electrophoresis result is shown in Fig. 3A. The predicted PCR amplification of the WT sequence length was 419 bp, while the MT sequence length was 266 bp, lacking the 153-bp exon 14. Sanger sequencing of the RT-PCR products confirmed that MT caused mRNA skipping over exon 14 (Fig. 3B). We have drawn a schematic diagram to illustrate aberrant splicing, as shown in Fig. 3C.
Many genes have been reported to be related to SRTD3, such as DYNC2H1, DYNC2LI1, WDR19, WDR34, WDR35, WDR60, IFT140, TTC21B, TCTEX1D2, IFT80, IFT172, EVC1, EVC2, KIAA0586, CEP120, INTU, ICK, NEK1, and C21ORF2 . Among the above genes, DYNC2H1 is the most common gene associated with SRTD3 disease. DYNC2H1 encodes a dynein protein that participates in ciliary intraflagellar transport . Mutations in DYNC2H1 will lead to the dysfunction of primary cilia, causing a heterogeneous spectrum of conditions such as skeletal dysplasias . DYNC2H1 is the pathogenic gene of SRTD3 disease. SRTD3 disease is clinically and genetically heterogeneous. There is significant variability in the course and severity of the thoracic phenotype, both between siblings with the same DYNC2H1 allele and between individuals with different alleles, suggesting that the DYNC2H1 phenotype may be influenced by modified alleles and nongenetic or epigenetic factors . More details about how each variant causes its phenotype await further study of the mechanisms of gene function.
SRTD3 generally presents prenatally with extremely shortened long bones and ribs, round metaphyseal ends and lateral spikes, a small and narrow thorax, and preaxial and postaxial polydactyly, sometimes presenting with bowing femora, genitals and multiple malformations . Among all the cases of reported SRTD, Okamoto et al. once reported a Japanese patient carrying a c.5682_5683 delAA and c.9010C > T mutations who showed an extremely severe phenotype, including rib and long bone shortness, severe thoracic hypoplasia, polydactyly, ventriculomegaly, macrocephaly, polyhydramnios, visceral heterotaxia, a congenital heart defect and hypoplasia of the iliac bones . Baujat et al. reported a baby in Turkey aged 9 months who was diagnosed with SRTD3 with situs inversus . As evidenced by this study of SRTD3 with situs inversus totalis and mirror-image dextrocardia at the prenatal stage, we added one more detailed piece of the situs inversus phenotype with SRTD3 in a Chinese foetus.
The DYNC2H1 compound heterozygous variants of the autosomal recessive inheritance pattern were consistent with SRTD3 disease. The foetal manifestation accorded with the SRTD3. No rare or likely pathogenic homozygous or heterozygous variant in other STRD-causing genes was found in family trio WES. We performed a minigene splicing study and confirmed that the variant caused a skipping over exon 14. Though no pathogenic variants were reported in this region previously, exon 14 was predicted to participate in many GO processes such as molecular_function, protein binding, catalytic activity and so on by NetGO 2.0 . According to the ACMG/AMP guidelines, c.11483T > G (p.Ile3828Arg) was classified as Variants of Uncertain Significance with the criteria PM2 and PP3 and c.2106 + 3A > T was classified as Likely Pathogenic with the criteria PM2, PP3 and PS3. On all these counts, the compound heterozygous variants in DYNC2H1 are highly likely the cause of the SRTD3 with polydactyly in this patient.
Here we presented a foetus with abnormal rib curvature, narrow thorax, bilateral hypoplastic lungs, bilateral polydactyly, syndactyly and foetal visceral situs inversus with mirror-image dextrocardia presentations caused by compound heterozygous variants in DYNC2H1. These findings expanded the spectrum of DYNC2H1 variants which was more diverse and predicted important molecular functions of exon 14 of DYNC2H1. The specific sonographic findings and the molecular diagnosis helped add experience to further our expertise in prenatal counselling for SRTD3. Our findings strongly supported the use of WES in prenatal diagnosis and further helped to understand the correlation of genotype and phenotypes of DYNC2H1 variants. Next-generation sequencing was an effective method of prenatal diagnosis of a genetic skeletal disorder.
Availability of data and materials
The novel variants revealed during the study were submitted to ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar/) under Accession Numbers SCV001806774 and SCV001806775.
Dagoneau N, Goulet M, Genevieve D, Sznajer Y, Martinovic J, Smithson S, Huber C, Baujat G, Flori E, Tecco L, et al. DYNC2H1 mutations cause asphyxiating thoracic dystrophy and short rib-polydactyly syndrome, type III. Am J Hum Genet. 2009;84(5):706–11.
Chen LS, Shi SJ, Zou PS, Ma M, Chen XH, Cao DH. Identification of novel DYNC2H1 mutations associated with short rib-polydactyly syndrome type III using next-generation panel sequencing. Genet Mol Res. 2016;15(2):8134.
Zhang W, Taylor SP, Ennis HA, Forlenza KN, Duran I, Li B, Sanchez JAO, Nevarez L, Nickerson DA, Bamshad M, et al. Expanding the genetic architecture and phenotypic spectrum in the skeletal ciliopathies. Hum Mutat. 2018;39(1):152–66.
Mei L, Huang Y, Pan Q, Su W, Quan Y, Liang D, Wu L. Targeted next-generation sequencing identifies novel compound heterozygous mutations of DYNC2H1 in a fetus with short rib-polydactyly syndrome, type III. Clin Chim Acta. 2015;447:47–51.
Wei X, Ju X, Yi X, Zhu Q, Qu N, Liu T, Chen Y, Jiang H, Yang G, Zhen R, et al. Identification of sequence variants in genetic disease-causing genes using targeted next-generation sequencing. PLOS ONE. 2011;6(12):e29500.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
San Lucas FA, Wang G, Scheet P, Peng B. Integrated annotation and analysis of genetic variants from next-generation sequencing studies with variant tools. Bioinformatics. 2012;28(3):421–2.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.
Ellis JA, Ong B. The MassARRAY((R)) system for targeted SNP genotyping. Methods Mol Biol. 2017;1492:77–94.
Badiner N, Taylor SP, Forlenza K, Lachman RS, University of Washington Center for Mendelian G, Bamshad M, Nickerson D, Cohn DH, Krakow D: Mutations in DYNC2H1, the cytoplasmic dynein 2, heavy chain 1 motor protein gene, cause short-rib polydactyly type I, Saldino-Noonan type. Clin Genet 2017; 92(2):158–165
Merrill AE, Merriman B, Farrington-Rock C, Camacho N, Sebald ET, Funari VA, Schibler MJ, Firestein MH, Cohn ZA, Priore MA, et al. Ciliary abnormalities due to defects in the retrograde transport protein DYNC2H1 in short-rib polydactyly syndrome. Am J Hum Genet. 2009;84(4):542–9.
Yi P, Li WJ, Dong MQ, Ou G. Dynein-driven retrograde intraflagellar transport is Triphasic in C. elegans Sensory Cilia. Curr Biol. 2017;27(10):1448–61.
Schmidts M, Arts HH, Bongers EM, Yap Z, Oud MM, Antony D, Duijkers L, Emes RD, Stalker J, Yntema JB, et al. Exome sequencing identifies DYNC2H1 mutations as a common cause of asphyxiating thoracic dystrophy (Jeune syndrome) without major polydactyly, renal or retinal involvement. J Med Genet. 2013;50(5):309–23.
Deng L, Cheung SW, Schmitt ES, Xiong S, Yuan M, Chen Z, Chen L, Sun L. Targeted gene panel sequencing prenatally detects two novel mutations of DYNC2H1 in a fetus with increased biparietal diameter and polyhydramnios. Birth Defects Res. 2018;110(4):364–71.
Okamoto T, Nagaya K, Kawata Y, Asai H, Tsuchida E, Nohara F, Okajima K, Azuma H. Novel compound heterozygous mutations in DYNC2H1 in a patient with severe short-rib polydactyly syndrome type III phenotype. Congenit Anom (Kyoto). 2015;55(3):155–7.
Baujat G, Huber C, El Hokayem J, Caumes R, Do Ngoc Thanh C, David A, Delezoide AL, Dieux-Coeslier A, Estournet B, Francannet C, et al. Asphyxiating thoracic dysplasia: clinical and molecular review of 39 families. J Med Genet. 2013;50(2):91–8.
Yao S, You R, Wang S, Xiong Y, Huang X, Zhu S. NetGO 20: improving large-scale protein function prediction with massive sequence, text, domain, family and network information. Nucleic Acids Res. 2021;49(W1):W469–75.
We thank the family participating in our study and the support of the programme of the Maternal and Child Health Hospital of Hubei Province.
This study was supported by the Hubei Province Natural Science Foundation under Grant Nos. 2020BCB002, 2020CFB164 and 2017CFB789 and the Hubei Province health and family planning scientific research project under Grant No. WJ2018H0132.
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the Maternal and Child Health Hospital of Hubei Province, and informed consent was obtained from the parents. Our research have been performed in accordance with the Declaration of Helsinki.
Consent for publication
The authors have declared no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cheng, C., Li, X., Zhao, S. et al. Compound heterozygous variants in DYNC2H1 in a foetus with type III short rib-polydactyly syndrome and situs inversus totalis. BMC Med Genomics 15, 55 (2022). https://doi.org/10.1186/s12920-022-01205-z