Deciphering the mutational signature of congenital limb malformations

Prevalence pattern, phenotypic manifestation, and descriptive genetics of congenital limb deficiencies in Pakistan

BACKGROUND

Congenital limb deficiency (CLD) is a group of very rare disorders characterized by substantial hypoplasia or the complete absence of 1 or more bones of limbs. Congenital limb deficiency has a significant physical, clinical, and psychological burden on the affected individuals and their families. This cross-sectional study aimed to describe the prevalence pattern, phenotypic manifestations, and biodemographic factors associated with CLD in a cohort assembled from the Pakistani population from the Northwestern region.

METHODS

Through a prospective cross-sectional study, 141 individuals having 166 limbs with CLD were recruited during 2017-2021.

RESULTS

There were 77 (55%) individuals with transverse defects, 61 (43%) with longitudinal defects, and 3 (2%) with Intercalary defects. Among the patients with transverse defects, 52 had terminal amputations and 25 had symbrachydactyly. Among the longitudinal defects, thumb aplasia/hypoplasia was the most common presentation (20 patients), followed by oligodactyly (18), and radial hemimelia (18). Eighty six percent had upper-limb deficiencies, 83% had unilateral deficiencies, and 92% were sporadic in nature. The parental consanguinity was observed in 33% individuals, and 79% cases had an isolated presentation which may be indicative of the substantial role of nongenetic factors in the etiology of CLD.

CONCLUSIONS

This study demonstrates marked heterogeneity in CLD subtypes in the involvement of limbs and associated variables. There is a need to establish a national registry for CLD, molecular genetic diagnosis, and multidisciplinary medical and social rehabilitation services for these individuals.

Clinical and genetic analysis in Chinese families with synpolydactyly, and cellular localization of HOXD13 with different length of polyalanine tract

Synpolydactyly (SPD) is caused by mutations in the transcription factor gene HOXD13. Such mutations include polyalanine expansion (PAE), but further study is required for the phenotypic spectrum characteristics of HOXD13 PAE. We investigated four unrelated Chinese families with significant limb malformations. Three PAEs were found in the HOXD13 polyalanine coding region: c.172_192dup (p.Ala58_Ala64dup) in Family 1, c.169_192dup (p.Ala57_Ala64dup) in Family 2, and c.183_210dup (p.Ala62_Ala70dup) in Family 3 and Family 4. Interestingly, we identified a new manifestation of preaxial polydactyly in both hands in a pediatric patient with an expansion of seven alanines, a phenotype not previously noted in SPD patients. Comparing with the wild-type cells and mutant cells with polyalanine contractions (PACs), the HOXD13 protein with a PAE of nine-alanine or more was difficult to enter the nucleus, and easy to form inclusion bodies in the cytoplasm, and with the increase of PAE, the more inclusion bodies were formed. This study not only expanded the phenotypic spectrum of SPD, but also enriched our understanding of its pathogenic mechanisms.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

Landscape of Secondary Findings in Chinese Population: A Practice of ACMG SF v3.0 List

Clinical exome sequencing (CES) has shown great utility in the diagnosis of Mendelian disorders. CES can unravel secondary findings (SFs) unrelated to the primary diagnosis but with potential health implications. The American College of Medical Genetics and Genomics (ACMG) has published a guideline for reporting secondary findings and recently updated an ACMG SF v3.0 list comprising 73 genes. Several studies have been performed to explore the prevalence of SFs. However, the data were limited in the Chinese population. In this study, we evaluated the genetic data of 2987 individuals from the Deciphering Disorders Involving Scoliosis and COmorbidities (DISCO) study group in accordance with the ACMG SF v3.0 list. The detected variants were evaluated using the ACMG classification guidelines, HGMD, and ClinVar database. Totally, 157 (157/2987, 5.3%) individuals had reportable variants within genes associated with cancer, cardiovascular, metabolic, and miscellaneous phenotypes. We identified 63 known pathogenic (KP) variants in 72 individuals (72/2987, 2.4%) and 96 expected pathogenic (EP) variants in 105 individuals (3.5%). Forty-five individuals carried SFs in v3.0 newly added genes, which accounted for 1.5% of our cohort. Our findings could contribute to existing knowledge of secondary findings in different ethnicities and indicate the necessity for clinicians to update the SFs gene list.

Identification of Novel FBN2 Variants in a Cohort of Congenital Contractural Arachnodactyly

Congenital contractural arachnodactyly (CCA) is a rare autosomal dominant disorder of connective tissue characterized by crumpled ears, arachnodactyly, camptodactyly, large joint contracture, and kyphoscoliosis. The nature course of CCA has not been well-described. We aim to decipher the genetic and phenotypic spectrum of CCA. The cohort was enrolled in Beijing Jishuitan Hospital and Peking Union Medical College Hospital, Beijing, China, based on Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study (http://www.discostudy.org/). Exome sequencing was performed on patients’ blood DNA. A recent published CCA scoring system was validated in our cohort. Seven novel variants and three previously reported FBN2 variants were identified through exome sequencing. Two variants outside of the neonatal region of FBN2 gene were found. The phenotypes were comparable between patients in our cohort and previous literature, with arachnodactyly, camptodactyly and large joints contractures found in almost all patients. All patients eligible for analysis were successfully classified into likely CCA based on the CCA scoring system. Furthermore, we found a double disease-causing heterozygous variant of FBN2 and ANKRD11 in a patient with blended phenotypes consisting of CCA and KBG syndrome. The identification of seven novel variants broadens the mutational and phenotypic spectrum of CCA and may provide implications for genetic counseling and clinical management.

PhenoApt leverages clinical expertise to prioritize candidate genes via machine learning

In recent years, exome sequencing (ES) has shown great utility in the diagnoses of Mendelian disorders. However, after rigorous filtering, a typical ES analysis still involves the interpretation of hundreds of variants, which greatly hinders the rapid identification of causative genes. Since the interpretations of ES data require comprehensive clinical analyses, taking clinical expertise into consideration can speed the molecular diagnoses of Mendelian disorders. To leverage clinical expertise to prioritize candidate genes, we developed PhenoApt, a phenotype-driven gene prioritization tool that allows users to assign a customized weight to each phenotype, via a machine-learning algorithm. Using the ability to rank causative genes in top-10 lists as an evaluation metric, baseline analysis demonstrated that PhenoApt outperformed previous phenotype-driven gene prioritization tools by a relative increase of 22.7%-140.0% in three independent, real-world, multi-center cohorts (cohort 1, n = 185; cohort 2, n = 784; and cohort 3, n = 208). Additional trials showed that, by adding weights to clinical indications, which should be explained by the causative gene, PhenoApt performance was improved by a relative increase of 37.3% in cohort 2 (n = 471) and 21.4% in cohort 3 (n = 208). Moreover, PhenoApt could assign an intrinsic weight to each phenotype based on the likelihood of its being a Mendelian trait using term frequency-inverse document frequency techniques. When clinical indications were assigned with intrinsic weights, PhenoApt performance was improved by a relative increase of 23.7% in cohort 2 and 15.5% in cohort 3. For the integration of PhenoApt into clinical practice, we developed a user-friendly website and a command-line tool.