Alu-Alu mediated intragenic duplications in IFT81 and MATN3 are associated with skeletal dysplasias

Adverse developmental impacts in progeny of zebrafish exposed to the agricultural herbicide atrazine during embryogenesis

Atrazine (ATZ) is an herbicide commonly used on crops in the Midwestern US and other select global regions. The US Environmental Protection Agency ATZ regulatory limit is 3 parts per billion (ppb; µg/L), but this limit is often exceeded. ATZ has a long half-life, is a common contaminant of drinking water sources, and is indicated as an endocrine disrupting chemical in multiple species. The zebrafish was used to test the hypothesis that an embryonic parental ATZ exposure alters protein levels leading to modifications in morphology and behavior in developing progeny. Zebrafish embryos (F1) were collected from adults (F0) exposed to 0, 0.3, 3, or 30 ppb ATZ during embryogenesis. Differential proteomics, morphology, and behavior assays were completed with offspring aged 120 or 144 h with no additional chemical treatment. Proteomic analysis identified differential expression of proteins associated with neurological development and disease; and organ and organismal morphology, development, and injury, specifically the skeletomuscular system. Head length and ratio of head length to total length was significantly increased in the F1 of 0.3 and 30 ppb ATZ groups (p < 0.05). Based on molecular pathway alterations, further craniofacial morphology assessment found decreased distance for cartilaginous structures, decreased surface area and distance between saccular otoliths, and a more posteriorly positioned notochord (p < 0.05), indicating delayed ossification and skeletal growth. The visual motor response assay showed hyperactivity in progeny of the 30 ppb treatment group for distance moved and of the 0.3 and 30 ppb treatment groups for time spent moving (p < 0.05). Due to the changes in saccular otoliths, an acoustic startle assay was completed and showed decreased response in the 0.3 and 30 ppb treatments (p < 0.05). These findings suggest that a single embryonic parental exposure alters cellular pathways in their progeny that lead to perturbations in craniofacial development and behavior.

Compound heterozygous IFT81 variations in a skeletal ciliopathy patient cause Bardet-Biedl syndrome-like ciliary defects

Owing to their crucial roles in development and homeostasis, defects in cilia cause ciliopathies with diverse clinical manifestations. The intraflagellar transport (IFT) machinery, containing the IFT-A and IFT-B complexes, mediates not only the intraciliary bidirectional trafficking but also import and export of ciliary proteins together with the kinesin-2 and dynein-2 motor complexes. The BBSome, containing eight subunits encoded by causative genes of Bardet-Biedl syndrome (BBS), connects the IFT machinery to ciliary membrane proteins to mediate their export from cilia. Although mutations in subunits of the IFT-A and dynein-2 complexes cause skeletal ciliopathies, mutations in some IFT-B subunits are also known to cause skeletal ciliopathies. We here show that compound heterozygous variations of an IFT-B subunit, IFT81, found in a patient with skeletal ciliopathy cause defects in its interactions with other IFT-B subunits, and in ciliogenesis and ciliary protein trafficking when one of the two variants was expressed in IFT81-knockout (KO) cells. Notably, we found that IFT81-KO cells expressing IFT81(Δ490-519), which lacks the binding site for the IFT25-IFT27 dimer, causes ciliary defects reminiscent of those found in BBS cells and those in IFT74-KO cells expressing a BBS variant of IFT74, which forms a heterodimer with IFT81. In addition, IFT81-KO cells expressing IFT81(Δ490-519) in combination with the other variant, IFT81 (L645*), which mimics the cellular conditions of the above skeletal ciliopathy patient, demonstrated essentially the same phenotype as those expressing only IFT81(Δ490-519). Thus, our data indicate that BBS-like defects can be caused by skeletal ciliopathy variants of IFT81.

Integration of Selection Signatures and Protein Interactions Reveals NR6A1, PAPPA2, and PIK3C2B as the Promising Candidate Genes Underlying the Characteristics of Licha Black Pig

Licha black (LI) pig has the specific characteristics of larger body length and appropriate fat deposition among Chinese indigenous pigs. Body length is one of the external traits that affect production performance, and fat deposition influences meat quality. However, the genetic characteristics of LI pigs have not yet been systematically uncovered. Here, the genomic information from 891 individuals of LI pigs, commercial pigs, and other Chinese indigenous pigs was used to analyze the breed characteristics of the LI pig with runs of homozygosity, haplotype, and FST selection signatures. The results showed the growth traits-related genes (i.e., NR6A1 and PAPPA2) and the fatness traits-related gene (i.e., PIK3C2B) were the promising candidate genes that closely related to the characteristics of LI pigs. In addition, the protein–protein interaction network revealed the potential interactions between the promising candidate genes and the FASN gene. The RNA expression data from FarmGTEx indicated that the RNA expression levels of NR6A1, PAPPA2, PIK3C2B, and FASN were highly correlated in the ileum. This study provides valuable molecular insights into the mechanisms that affect pig body length and fat deposition, which can be used in the further breeding process to improve meat quality and commercial profitability.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Zebrafish Models for Human Skeletal Disorders

In 2019, the Nosology Committee of the International Skeletal Dysplasia Society provided an updated version of the Nosology and Classification of Genetic Skeletal Disorders. This is a reference list of recognized diseases in humans and their causal genes published to help clinician diagnosis and scientific research advances. Complementary to mammalian models, zebrafish has emerged as an interesting species to evaluate chemical treatments against these human skeletal disorders. Due to its versatility and the low cost of experiments, more than 80 models are currently available. In this article, we review the state-of-art of this “aquarium to bedside” approach describing the models according to the list provided by the Nosology Committee. With this, we intend to stimulate research in the appropriate direction to efficiently meet the actual needs of clinicians under the scope of the Nosology Committee.

High diagnostic yield in skeletal ciliopathies using massively parallel genome sequencing, structural variant screening and RNA analyses

Skeletal ciliopathies are a heterogenous group of disorders with overlapping clinical and radiographic features including bone dysplasia and internal abnormalities. To date, pathogenic variants in at least 30 genes, coding for different structural cilia proteins, are reported to cause skeletal ciliopathies. Here, we summarize genetic and phenotypic features of 34 affected individuals from 29 families with skeletal ciliopathies. Molecular diagnostic testing was performed using massively parallel sequencing (MPS) in combination with copy number variant (CNV) analyses and in silico filtering for variants in known skeletal ciliopathy genes. We identified biallelic disease-causing variants in seven genes: DYNC2H1, KIAA0753, WDR19, C2CD3, TTC21B, EVC, and EVC2. Four variants located in non-canonical splice sites of DYNC2H1, EVC, and KIAA0753 led to aberrant splicing that was shown by sequencing of cDNA. Furthermore, CNV analyses showed an intragenic deletion of DYNC2H1 in one individual and a 6.7 Mb de novo deletion on chromosome 1q24q25 in another. In five unsolved cases, MPS was performed in family setting. In one proband we identified a de novo variant in PRKACA and in another we found a homozygous intragenic deletion of IFT74, removing the first coding exon and leading to expression of a shorter message predicted to result in loss of 40 amino acids at the N-terminus. These findings establish IFT74 as a new skeletal ciliopathy gene. In conclusion, combined single nucleotide variant, CNV and cDNA analyses lead to a high yield of genetic diagnoses (90%) in a cohort of patients with skeletal ciliopathies.

Diagnostic Approach to Macrocephaly in Children

Macrocephaly affects up to 5% of the pediatric population and is defined as an abnormally large head with an occipitofrontal circumference (OFC) >2 standard deviations (SD) above the mean for a given age and sex. Taking into account that about 2-3% of the healthy population has an OFC between 2 and 3 SD, macrocephaly is considered as "clinically relevant" when OFC is above 3 SD. This implies the urgent need for a diagnostic workflow to use in the clinical setting to dissect the several causes of increased OFC, from the benign form of familial macrocephaly and the Benign enlargement of subarachnoid spaces (BESS) to many pathological conditions, including genetic disorders. Moreover, macrocephaly should be differentiated by megalencephaly (MEG), which refers exclusively to brain overgrowth, exceeding twice the SD (3SD-"clinically relevant" megalencephaly). While macrocephaly can be isolated and benign or may be the first indication of an underlying congenital, genetic, or acquired disorder, megalencephaly is most likely due to a genetic cause. Apart from the head size evaluation, a detailed family and personal history, neuroimaging, and a careful clinical evaluation are crucial to reach the correct diagnosis. In this review, we seek to underline the clinical aspects of macrocephaly and megalencephaly, emphasizing the main differential diagnosis with a major focus on common genetic disorders. We thus provide a clinico-radiological algorithm to guide pediatricians in the assessment of children with macrocephaly.

Expanding the phenotypic spectrum of IFT81: Associated ciliopathy syndrome

Short-rib polydactyly syndromes are a heterogeneous group of disorders characterized by narrow thorax with short ribs, polydactyly and often other visceral and skeletal malformations. To date there have only been six reported patients with homozygous and compound heterozygous variants in IFT81, causing a short-rib thoracic dysplasia, with, or without, polydactyly (SRTD19: OMIM 617895). IFT81 is a protein integral to the core of the intraflagellar transport complex B (IFT-B), which is involved in anterograde transport in the cilium. We describe the case of a male infant with compound heterozygous variants in IFT81, who presented with short long bones, a narrow thorax, polydactyly, and multiple malformations. Three novel clinical features are reported including complete situs inversus, micropenis, and rectal atresia, which have not previously been associated with variants in IFT81. We reviewed the literature and identified the most consistent clinical features associated with this rare ciliopathy syndrome. We postulate that dolichocephaly and sagittal craniosynostosis may be associated with this condition, and provide a clue to considering IFT81 as the causative gene when deciphering complex ciliopathies.

Multiple epiphyseal dysplasia and related disorders: Molecular genetics, disease mechanisms, and therapeutic avenues

For the vast majority of the 6000 known rare disease the pathogenic mechanisms are poorly defined and there is little treatment, leading to poor quality of life and high healthcare costs. Genetic skeletal diseases (skeletal dysplasias) are archetypal examples of rare diseases that are chronically debilitating, often life-threatening and for which no treatments are currently available. There are more than 450 unique phenotypes that, although individually rare, have an overall prevalence of at least 1 per 4000 children. Multiple epiphyseal dysplasia (MED) is a clinically and genetically heterogeneous disorder characterized by disproportionate short stature, joint pain, and early-onset osteoarthritis. MED is caused by mutations in the genes encoding important cartilage extracellular matrix proteins, enzymes, and transporter proteins. Recently, through the use of various cell and mouse models, disease mechanisms underlying this diverse phenotypic spectrum are starting to be elucidated. For example, ER stress induced as a consequence of retained misfolded mutant proteins has emerged as a unifying disease mechanisms for several forms of MED in particular and skeletal dysplasia in general. Moreover, targeting ER stress through drug repurposing has become an attractive therapeutic avenue.

The essential role of intraflagellar transport protein IFT81 in male mice spermiogenesis and fertility

Intraflagellar transport (IFT) is an evolutionarily conserved mechanism that is indispensable for the formation and maintenance of cilia and flagella; however, the implications and functions of IFT81 remain unknown. In this study, we disrupted IFT81 expression in male germ cells starting from the spermatocyte stage. As a result, homozygous mutant males were completely infertile and displayed abnormal sperm parameters. In addition to oligozoospermia, spermatozoa presented dysmorphic and nonfunctional flagella. Histological examination of testes from homozygous mutant mice revealed abnormal spermiogenesis associated with sloughing of germ cells and the presence of numerous multinucleated giant germ cells (symblasts) in the lumen of seminiferous tubules and epididymis. Moreover, only few elongated spermatids and spermatozoa were seen in analyzed cross sections. Transmission electron microscopy showed a complete disorganization of the axoneme and para-axonemal structures such as the mitochondrial sheath, fibrous sheath, and outer dense fibers. In addition, numerous vesicles that contain unassembled microtubules were observed within developing spermatids. Acrosome structure analysis showed normal appearance, thus excluding a crucial role of IFT81 in acrosome biogenesis. These observations showed that IFT81 is an important member of the IFT process during spermatogenesis and that its absence is associated with abnormal flagellum formation leading to male infertility. The expression levels of several IFT components in testes, including IFT20, IFT25, IFT27, IFT57, IFT74, and IFT88, but not IFT140, were significantly reduced in homozygous mutant mice. Overall, our study demonstrates that IFT81 plays an essential role during spermatogenesis by modulating the assembly and elongation of the sperm flagella.

From cytogenetics to cytogenomics: whole-genome sequencing as a first-line test comprehensively captures the diverse spectrum of disease-causing genetic variation underlying intellectual disability

Background

Since different types of genetic variants, from single nucleotide variants (SNVs) to large chromosomal rearrangements, underlie intellectual disability, we evaluated the use of whole-genome sequencing (WGS) rather than chromosomal microarray analysis (CMA) as a first-line genetic diagnostic test.

Methods

We analyzed three cohorts with short-read WGS: (i) a retrospective cohort with validated copy number variants (CNVs) (cohort 1, n = 68), (ii) individuals referred for monogenic multi-gene panels (cohort 2, n = 156), and (iii) 100 prospective, consecutive cases referred to our center for CMA (cohort 3). Bioinformatic tools developed include FindSV, SVDB, Rhocall, Rhoviz, and vcf2cytosure.

Results

First, we validated our structural variant (SV)-calling pipeline on cohort 1, consisting of three trisomies and 79 deletions and duplications with a median size of 850 kb (min 500 bp, max 155 Mb). All variants were detected. Second, we utilized the same pipeline in cohort 2 and analyzed with monogenic WGS panels, increasing the diagnostic yield to 8%. Next, cohort 3 was analyzed by both CMA and WGS. The WGS data was processed for large (> 10 kb) SVs genome-wide and for exonic SVs and SNVs in a panel of 887 genes linked to intellectual disability as well as genes matched to patient-specific Human Phenotype Ontology (HPO) phenotypes. This yielded a total of 25 pathogenic variants (SNVs or SVs), of which 12 were detected by CMA as well. We also applied short tandem repeat (STR) expansion detection and discovered one pathologic expansion in ATXN7. Finally, a case of Prader-Willi syndrome with uniparental disomy (UPD) was validated in the WGS data.

Important positional information was obtained in all cohorts. Remarkably, 7% of the analyzed cases harbored complex structural variants, as exemplified by a ring chromosome and two duplications found to be an insertional translocation and part of a cryptic unbalanced translocation, respectively.

Conclusion

The overall diagnostic rate of 27% was more than doubled compared to clinical microarray (12%). Using WGS, we detected a wide range of SVs with high accuracy. Since the WGS data also allowed for analysis of SNVs, UPD, and STRs, it represents a powerful comprehensive genetic test in a clinical diagnostic laboratory setting.

High Fidelity of Mouse Models Mimicking Human Genetic Skeletal Disorders

The 2019 International Skeletal Dysplasia Society nosology update lists 441 genes for which mutations result in rare human skeletal disorders. These genes code for enzymes (33%), scaffolding proteins (18%), signal transduction proteins (16%), transcription factors (14%), cilia proteins (8%), extracellular matrix proteins (5%), and membrane transporters (4%). Skeletal disorders include aggrecanopathies, channelopathies, ciliopathies, cohesinopathies, laminopathies, linkeropathies, lysosomal storage diseases, protein-folding and RNA splicing defects, and ribosomopathies. With the goal of evaluating the ability of mouse models to mimic these human genetic skeletal disorders, a PubMed literature search identified 260 genes for which mutant mice were examined for skeletal phenotypes. These mouse models included spontaneous and ENU-induced mutants, global and conditional gene knockouts, and transgenic mice with gene over-expression or specific base-pair substitutions. The human X-linked gene ARSE and small nuclear RNA U4ATAC, a component of the minor spliceosome, do not have mouse homologs. Mouse skeletal phenotypes mimicking human skeletal disorders were observed in 249 of the 260 genes (96%) for which comparisons are possible. A supplemental table in spreadsheet format provides PubMed weblinks to representative publications of mutant mouse skeletal phenotypes. Mutations in 11 mouse genes (Ccn6, Cyp2r1, Flna, Galns, Gna13, Lemd3, Manba, Mnx1, Nsd1, Plod1, Smarcal1) do not result in similar skeletal phenotypes observed with mutations of the homologous human genes. These discrepancies can result from failure of mouse models to mimic the exact human gene mutations. There are no obvious commonalities among these 11 genes. Body BMD and/or radiologic dysmorphology phenotypes were successfully identified for 28 genes by the International Mouse Phenotyping Consortium (IMPC). Forward genetics using ENU mouse mutagenesis successfully identified 37 nosology gene phenotypes. Since many human genetic disorders involve hypomorphic, gain-of-function, dominant-negative and intronic mutations, future studies will undoubtedly utilize CRISPR/Cas9 technology to examine transgenic mice having genes modified to exactly mimic variant human sequences. Mutant mice will increasingly be employed for drug development studies designed to treat human genetic skeletal disorders.

SIGNIFICANCE

Great progress is being made identifying mutant genes responsible for human rare genetic skeletal disorders and mouse models for genes affecting bone mass, architecture, mineralization and strength. This review organizes data for 441 human genetic bone disorders with regard to heredity, gene function, molecular pathways, and fidelity of relevant mouse models to mimic the human skeletal disorders. PubMed weblinks to citations of 249 successful mouse models are provided.