Search

Neurofibromatosis type 1 (NF1) is a tumor predisposition syndrome caused by heterozygous NF1 gene mutations. Patients with NF1 present with pleiotropic somatic secondary manifestations, including development of bone pseudarthrosis after fracture. Somatic NF1 gene mutations were reproducibly identified in patient-derived pseudarthrosis specimens, suggesting a local mosaic cell population including somatic pathologic cells. The somatic cellular pathogenesis of NF1 pseudarthroses remains unclear, though defects in osteogenesis have been posited. Here, we applied time-series single-cell RNA-sequencing (scRNA-seq) to patient-matched control and pseudarthrosis-derived primary bone stromal cells (BSCs). We show that osteogenic specification to an osteoblast progenitor cell population was evident for control bone-derived cells and haploinsufficient pseudarthrosis-derived cells. Similar results were observed for somatic patient fracture-derived NF1

Proper embryonic and postnatal skeletal development require coordination of myriad complex molecular mechanisms. Disruption of these processes, through genetic mutation, contributes to variation in skeletal development. We developed a high-throughput N-ethyl-N-nitrosourea (ENU)-induced saturation mutagenesis skeletal screening approach in mice to identify genes required for proper skeletal development. Here, we report initial results from live-animal X-ray and dual-energy X-ray absorptiometry (DXA) imaging of 27,607 G3 mice from 806 pedigrees, testing the effects of 32,198 coding/splicing mutations in 13,020 genes. A total of 39.7% of all autosomal genes were severely damaged or destroyed by mutations tested twice or more in the homozygous state. Results from our study demonstrate the feasibility of in vivo mutagenesis to identify mouse models of skeletal disease. Furthermore, our study demonstrates how ENU mutagenesis provides opportunities to create and characterize putative hypomorphic mutations in developmentally essential genes. Finally, we present a viable mouse model and case report of recessive skeletal disease caused by mutations in FAM20B. Results from this study, including engineered mouse models, are made publicly available via the online Mutagenetix database. © 2021 American Society for Bone and Mineral Research (ASBMR).

Talipes equinovarus (clubfoot, TEV) is a congenital rotational foot deformity occurring in 1 per 1000 births with increased prevalence in males compared with females. The genetic etiology of isolated clubfoot (iTEV) remains unclear. Using a genome-wide association study, we identified a locus within FSTL5, encoding follistatin-like 5, significantly associated with iTEV. FSTL5 is an uncharacterized gene whose potential role in embryonic and postnatal development was previously unstudied. Utilizing multiple model systems, we found that Fstl5 was expressed during later stages of embryonic hindlimb development, and, in mice, expression was restricted to the condensing cartilage anlage destined to form the limb skeleton. In the postnatal growth plate, Fstl5 was specifically expressed in prehypertrophic chondrocytes. As Fstl5 knockout rats displayed no gross malformations, we engineered a conditional transgenic mouse line (Fstl5LSL) to overexpress Fstl5 in skeletal osteochondroprogenitors. We observed that hindlimbs were slightly shorter and that bone mineral density was reduced in adult male, but not female, Prrx1-cre;Fstl5LSL mice compared with control. No overt clubfoot-like deformity was observed in Prrx1-cre;Fstl5LSL mice, suggesting FSTL5 may function in other cell types to contribute to iTEV pathogenesis. Interrogating published mouse embryonic single-cell expression data showed that Fstl5 was expressed in cell lineage subclusters whose transcriptomes were associated with neural system development. Moreover, our results suggest that lineage-specific expression of the Fstl genes correlates with their divergent roles as modulators of transforming growth factor beta and bone morphogenetic protein signaling. Results from this study associate FSTL5 with iTEV and suggest a potential sexually dimorphic role for Fstl5 in vivo.

During genome replication, polymerase epsilon (Pol ε) acts as the major leading-strand DNA polymerase. Here we report the identification of biallelic mutations in POLE, encoding the Pol ε catalytic subunit POLE1, in 15 individuals from 12 families. Phenotypically, these individuals had clinical features closely resembling IMAGe syndrome (intrauterine growth restriction [IUGR], metaphyseal dysplasia, adrenal hypoplasia congenita, and genitourinary anomalies in males), a disorder previously associated with gain-of-function mutations in CDKN1C. POLE1-deficient individuals also exhibited distinctive facial features and variable immune dysfunction with evidence of lymphocyte deficiency. All subjects shared the same intronic variant (c.1686+32C>G) as part of a common haplotype, in combination with different loss-of-function variants in trans. The intronic variant alters splicing, and together the biallelic mutations lead to cellular deficiency of Pol ε and delayed S-phase progression. In summary, we establish POLE as a second gene in which mutations cause IMAGe syndrome. These findings add to a growing list of disorders due to mutations in DNA replication genes that manifest growth restriction alongside adrenal dysfunction and/or immunodeficiency, consolidating these as replisome phenotypes and highlighting a need for future studies to understand the tissue-specific development roles of the encoded proteins.

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093=1.7E-15), ABO (P-value_ rs687621=7.3E-10) and SOX6 (P-value_rs1455114=2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Adolescent idiopathic scoliosis (AIS) is the most common pediatric musculoskeletal disorder, affecting about 2–3% of the children worldwide. The underlying etiology of AIS has proven elusive, in part due to the lack of appropriate animal models and systems that would enable laboratory-based research. In contrast, genetic studies in patient populations have highlighted several candidate genes suggesting possible neuromuscular origins. Developing genetically defined animal models to facilitate hypothesis testing is a high priority for ongoing AIS research. Continued gene discovery efforts supported by next-generation genomic platforms will yield exciting new insights into AIS disease pathways, opening the prospect for pharmaceutical interventions.

Neurofibromatosis type 1 (NF1) is an autosomal dominant disease caused by mutations in NF1. Among the earliest manifestations is tibial pseudoarthrosis and persistent nonunion after fracture. To further understand the pathogenesis of pseudoarthrosis and the underlying bone remodeling defect, pseudoarthrosis tissue and cells cultured from surgically resected pseudoarthrosis tissue from NF1 individuals were analyzed using whole-exome and whole-transcriptome sequencing as well as genomewide microarray analysis. Genomewide analysis identified multiple genetic mechanisms resulting in somatic biallelic NF1 inactivation; no other genes with recurring somatic mutations were identified. Gene expression profiling identified dysregulated pathways associated with neurofibromin deficiency, including phosphoinositide 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling pathways. Unlike aggressive NF1-associated malignancies, tibial pseudoarthrosis tissue does not harbor a high frequency of somatic mutations in oncogenes or other tumor-suppressor genes, such as p53. However, gene expression profiling indicates that pseudoarthrosis tissue has a tumor-promoting transcriptional pattern, despite lacking tumorigenic somatic mutations. Significant overexpression of specific cancer-associated genes in pseudoarthrosis highlights a potential for receptor tyrosine kinase inhibitors to target neurofibromin-deficient pseudoarthrosis and promote proper bone remodeling and fracture healing.

Neurofibromatosis type 1 (NF1) is a common genetic disorder caused by mutations in the NF1 gene. Recalcitrant bone healing following fracture (i.e. pseudarthrosis) is one of the most problematic skeletal complications associated with NF1. The etiology of this condition is still unclear; thus, pharmacological options for clinical management are limited. Multiple studies have shown the reduced osteogenic potential of Nf1-deficient osteoprogenitors. A recent transcriptome profiling investigation revealed that EREG and EGFR, encoding epiregulin and its receptor Epidermal Growth Factor Receptor 1, respectively, were among the top over-expressed genes in cells of the NF1 pseudarthrosis site. Because EGFR stimulation is known to inhibit osteogenic differentiation, we hypothesized that increased EREG and EGFR expression in NF1-deficient skeletal progenitors may contribute to their reduced osteogenic differentiation potential. In this study, we first confirmed via single-cell mRNA sequencing that EREG over-expression was associated with NF1 second hit somatic mutations in human bone cells, whereas Transforming Growth Factor beta 1 (TGFβ1) expression was unchanged. Second, using ex-vivo recombined Nf1-deficient mouse bone marrow stromal cells (mBMSCs), we show that this molecular signature is conserved between mice and humans, and that epiregulin generated by these cells is overexpressed and active, whereas soluble TGFβ1 expression and activity are not affected. However, blocking either epiregulin function or EGFR signaling by EGFR1 or pan EGFR inhibition (using AG-1478 and Poziotinib respectively) did not correct the differentiation defect of Nf1-deficient mBMSCs, as measured by the expression of Alpl, Ibsp and alkaline phosphatase activity. These results suggest that clinically available pharmacological strategies aimed at inhibiting EGFR signaling are unlikely to have a significant benefit for the management of bone non-union in children with NF1 PA.

Macrodactyly is a discrete congenital anomaly consisting of enlargement of all tissues localized to the terminal portions of a limb, typically within a 'nerve territory'. The classic terminology for this condition is 'lipofibromatous hamartoma of nerve' or Type I macrodactyly. The peripheral nerve, itself, is enlarged both in circumference and in length. It is not related to neurofibromatosis (NF1), nor is it associated with vascular malformations, such as in the recently reported CLOVES syndrome. The specific nerve pathophysiology in this form of macrodactyly has not been well described and a genetic etiology for this specific form of enlargement is unknown. To identify the genetic cause of macrodactyly, we used whole-exome sequencing to identify somatic mutations present in the affected nerve of a single patient. We confirmed a novel mutation in PIK3CA (R115P) present in the patient's affected nerve tissue but not in blood DNA. Sequencing PIK3CA exons identified gain-of-function mutations (E542K, H1047L or H1047R) in the affected tissue of five additional unrelated patients; mutations were absent in blood DNA available from three patients. Immunocytochemistry confirmed AKT activation in cultured cells from the nerve of a macrodactyly patient. Additionally, we found that the most abnormal structure within the involved nerve in a macrodactylous digit is the perineurium, with additional secondary effects on the axon number and size. Thus, isolated congenital macrodactyly is caused by somatic activation of the PI3K/AKT cell-signaling pathway and is genetically and biochemically related to other overgrowth syndromes.

The periosteum contributes to bone repair and maintenance of cortical bone mass. In contrast to the understanding of bone development within the epiphyseal growth plate, factors that regulate periosteal osteogenesis have not been studied as intensively. Osteofibrous dysplasia (OFD) is a congenital disorder of osteogenesis and is typically sporadic and characterized by radiolucent lesions affecting the cortical bone immediately under the periosteum of the tibia and fibula. We identified germline mutations in MET, encoding a receptor tyrosine kinase, that segregate with an autosomal-dominant form of OFD in three families and a mutation in a fourth affected subject from a simplex family and with bilateral disease. Mutations identified in all families with dominant inheritance and in the one simplex subject with bilateral disease abolished the splice inclusion of exon 14 in MET transcripts, which resulted in a MET receptor (MET(Δ14)) lacking a cytoplasmic juxtamembrane domain. Splice exclusion of this domain occurs during normal embryonic development, and forced induction of this exon-exclusion event retarded osteoblastic differentiation in vitro and inhibited bone-matrix mineralization. In an additional subject with unilateral OFD, we identified a somatic MET mutation, also affecting exon 14, that substituted a tyrosine residue critical for MET receptor turnover and, as in the case of the MET(Δ14) mutations, had a stabilizing effect on the mature protein. Taken together, these data show that aberrant MET regulation via the juxtamembrane domain subverts core MET receptor functions that regulate osteogenesis within cortical diaphyseal bone.

Polymorphisms associated with genes coding for glutathione S-transferase enzymes are known to influence metabolism of different carcinogens and have been associated with incidence of various types of cancer. We have determined the GST M1 and GST T1 'null' genotype frequency in 81 patients with chronic myeloid leukaemia (CML) and 123 racially and geographically matched control individuals by multiplex polymerase chain reaction (PCR). GST M1 null genotype frequencies in CML and controls were 28.4% and 27.7%, respectively. GST T1 null genotype frequencies in CML and controls were 19.8% and 7.3%, respectively. The GST T1 null genotype frequency in CML patients is significantly different from that in controls (odds ratio (OR) 3.12, 95% confidence interval (CI) 1.3-7.45, P=0.008).

Scoliosis is a complex genetic disorder of the musculoskeletal system, characterized by three-dimensional rotation of the spine. Curvatures caused by malformed vertebrae (congenital scoliosis (CS)) are apparent at birth. Spinal curvatures with no underlying vertebral abnormality (idiopathic scoliosis (IS)) most commonly manifest during adolescence. The genetic and biological mechanisms responsible for IS remain poorly understood due largely to limited experimental models. Here we describe zygotic ptk7 (Zptk7) mutant zebrafish, deficient in a critical regulator of Wnt signalling, as the first genetically defined developmental model of IS. We identify a novel sequence variant within a single IS patient that disrupts PTK7 function, consistent with a role for dysregulated Wnt activity in disease pathogenesis. Furthermore, we demonstrate that embryonic loss-of-gene function in maternal-zygotic ptk7 mutants (MZptk7) leads to vertebral anomalies associated with CS. Our data suggest novel molecular origins of, and genetic links between, congenital and idiopathic forms of disease., Scoliosis is a complex genetic disorder characterized by spinal curvature. Here, the authors present experimental zebrafish models of idiopathic and congenital scoliosis and suggest a role for dysregulated Wnt activity in scoliosis aetiology.

Neurofibromatosis type 1 (NF1) is an autosomal dominant disease caused by mutations in NF1. Among the earliest manifestations is tibial pseudoarthrosis and persistent nonunion after fracture. To further understand the pathogenesis of pseudoarthrosis and the underlying bone remodeling defect, pseudoarthrosis tissue and cells cultured from surgically resected pseudoarthrosis tissue from NF1 individuals were analyzed using whole-exome and whole-transcriptome sequencing as well as genomewide microarray analysis. Genomewide analysis identified multiple genetic mechanisms resulting in somatic biallelic NF1 inactivation; no other genes with recurring somatic mutations were identified. Gene expression profiling identified dysregulated pathways associated with neurofibromin deficiency, including phosphoinositide 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling pathways. Unlike aggressive NF1-associated malignancies, tibial pseudoarthrosis tissue does not harbor a high frequency of somatic mutations in oncogenes or other tumor-suppressor genes, such as p53. However, gene expression profiling indicates that pseudoarthrosis tissue has a tumor-promoting transcriptional pattern, despite lacking tumorigenic somatic mutations. Significant overexpression of specific cancer-associated genes in pseudoarthrosis highlights a potential for receptor tyrosine kinase inhibitors to target neurofibromin-deficient pseudoarthrosis and promote proper bone remodeling and fracture healing.

Genetic studies of orthopaedic conditions, particularly those of childhood onset, enjoyed many successes in the early days of disease gene discovery. Some of the earliest high-profile discoveries centered on diseases with skeletal manifestations, such as neurofibromatosis1,2, achondroplasia3,4, and brachydactyly5-8. These successes were largely attributable to the availability of extended families with multiple affected members. Multigenerational families were informative in linkage mapping, also called inheritance mapping (Table I). Linkage mapping is an unbiased approach that seeks to identify genomic regions and genes that are inherited or “linked” among all affected members of the family. Linkage mapping proved particularly powerful for solving diseases that (1) have Mendelian inheritance, being caused by a single mutation in a family; (2) display high penetrance, such that essentially all family members with the mutation show the disease; and (3) are prevalent enough to permit the identification of several extended families. Linking a disease to a specific region of the genome was an important milestone for many early genetic studies, with such highlights as the mapping of the neurofibromatosis-1 gene (NF1) to chromosome 17q11.2 and the achondroplasia gene (FGFR3, which encodes fibroblast growth factor receptor 3) to 4p16.3. Once linkage mapping focused attention on a particular genomic region, the next major step in discovery was analysis of the region to identify the culprit gene harboring a disease-causing mutation. This process was labor and resource-intensive because linked regions were typically large and poorly annotated and encoded numerous genes9. A now classic, and celebrated, example of the immense effort required to find a disease mutation is the almost ten-year search for the Huntington disease gene10. Nevertheless, the combination of linkage mapping and individual gene sequencing in families, traditionally called positional cloning, was highly …

Sequencing technologies promising the “$1000 genome” have developed at a staggering pace, driven in part by the HapMap Project1 (which set out to catalog common variations throughout the genome in multiple ethnic populations) and the 1000 Genomes Project2 (a consortium to catalog rare variations from different ethnic populations worldwide). Commercial whole-genome sequencing platforms have yet to break the $1000 mark. However, sequencing throughput and accuracy are no longer limiting factors. Despite recent successes in patient-oriented whole-genome sequencing3-6, analyzing and interpreting these data sets represent a major challenge, particularly for the ∼98% of the genome that does not encode proteins. In 2009, an alternative to whole-genome sequencing was introduced7. Whole-exome sequencing uses high-throughput next-generation technologies to sequence only the portions of the genome that encode proteins (termed the exome). The goal of whole-exome sequencing is to identify sequence mutations that alter the amino acid content of a protein, potentially altering the protein’s function in a cell and leading to disease. Proof of concept for this approach was documented in a landmark paper that used whole-exome sequencing to identify the genetic cause of Freeman-Sheldon syndrome (OMIM [Online Mendelian Inheritance in Man] 193700)7. This study broke the technical and cost limitation barriers (compared with sequencing an entire genome) by sequencing only the protein-coding regions (∼1% of the genome) and showed that a study of only a few patients was sufficient to identify a disease-causing gene. Whole-exome sequencing has revolutionized disease-gene discovery and, supported by the 1000 Genomes Project2 and other large-scale sequencing studies (e.g., the NHLBI/NHGRI [National Heart, Lung, and Blood Institute/National Human Genome Research Institute] Exome Project), it promises to identify virtually all genes causing disorders with Mendelian inheritance—i.e., those disorders caused by single genes. In this review, we focus on …

Mature mammalian sperm contain a complex population of RNAs some of which might regulate spermatogenesis while others probably play a role in fertilization and early development. Due to this limited knowledge, the biological functions of sperm RNAs remain enigmatic. Here we report the first characterization of the global transcriptome of the sperm of fertile stallions. The findings improved understanding of the biological significance of sperm RNAs which in turn will allow the discovery of sperm-based biomarkers for stallion fertility. The stallion sperm transcriptome was interrogated by analyzing sperm and testes RNA on a 21,000-element equine whole-genome oligoarray and by RNA-seq. Microarray analysis revealed 6,761 transcripts in the sperm, of which 165 were sperm-enriched, and 155 were differentially expressed between the sperm and testes. Next, 70 million raw reads were generated by RNA-seq of which 50% could be aligned with the horse reference genome. A total of 19,257 sequence tags were mapped to all horse chromosomes and the mitochondrial genome. The highest density of mapped transcripts was in gene-rich ECA11, 12 and 13, and the lowest in gene-poor ECA9 and X; 7 gene transcripts originated from ECAY. Structural annotation aligned sperm transcripts with 4,504 known horse and/or human genes, rRNAs and 82 miRNAs, whereas 13,354 sequence tags remained anonymous. The data were aligned with selected equine gene models to identify additional exons and splice variants. Gene Ontology annotations showed that sperm transcripts were associated with molecular processes (chemoattractant-activated signal transduction, ion transport) and cellular components (membranes and vesicles) related to known sperm functions at fertilization, while some messenger and micro RNAs might be critical for early development. The findings suggest that the rich repertoire of coding and non-coding RNAs in stallion sperm is not a random remnant from spermatogenesis in testes but a selectively retained and functionally coherent collection of RNAs.

Studies of the Y chromosome in primates, rodents and carnivores provide compelling evidence that the male specific region of Y (MSY) contains functional genes, many of which have specialized roles in spermatogenesis and male-fertility. Little similarity, however, has been found between the gene content and sequence of MSY in different species. This hinders the discovery of species-specific male fertility genes and limits our understanding about MSY evolution in mammals. Here, a detailed MSY gene catalogue was developed for the horse – an odd-toed ungulate. Using direct cDNA selection from horse testis, and sequence analysis of Y-specific BAC clones, 37 horse MSY genes/transcripts were identified. The genes were mapped to the MSY BAC contig map, characterized for copy number, analyzed for transcriptional profiles by RT-PCR, examined for the presence of ORFs, and compared to other mammalian orthologs. We demonstrate that the horse MSY harbors 20 X-degenerate genes with known orthologs in other eutherian species. The remaining 17 genes are acquired or novel and have so far been identified only in the horse or donkey Y chromosomes. Notably, 3 transcripts were found in the heterochromatic part of the Y. We show that despite substantial differences between the sequence, gene content and organization of horse and other mammalian Y chromosomes, the functions of MSY genes are predominantly related to testis and spermatogenesis. Altogether, 10 multicopy genes with testis-specific expression were identified in the horse MSY, and considered likely candidate genes for stallion fertility. The findings establish an important foundation for the study of Y-linked genetic factors governing fertility in stallions, and improve our knowledge about the evolutionary processes that have shaped Y chromosomes in different mammalian lineages.

Sperm mRNA transcriptional profiles can be used to evaluate male fertility, yet differences between species in sperm attributes and packaging require adjustments in sperm RNA isolation protocols. The objective was to optimize RNA isolation methodology for fresh, frozen, and extended ejaculates, and epididymal sperm of stallions. Additionally, a protocol for RNA isolation from testis biopsies was established. Separation of sperm from somatic cells was critical for assuring the isolation of sperm-specific RNA. The highest purity was obtained by centrifuging ejaculates and epididymal sperm at 200 x g for 30 min through a 40% Equipure silanized silica particle solution. Sperm RNA isolation was more efficient with TRIzol reagent than with a spin-column based method; it resulted in 2 microg of total RNA per 100 x 10(6) sperm. To evaluate RNA quantity and quality, we used a NanoDrop spectrophotometer and Agilent Bioanalyzer. A protocol for reverse transcriptase PCR with equine primers for PRM2 and PTPRC genes was developed to determine sperm RNA contamination with genomic DNA or RNA from somatic cells. By these methods, hybridization- and sequencing-quality RNA was isolated from 11 samples of stallion sperm. Stallion testis biopsy with a 14 gauge 22 mm deep biopsy needle yielded approximately 12 microg of good quality total RNA, and could serve as an alternative to excision surgery for sample procurement. Compared to RNA isolation from testis, the sperm required advanced processing and RNA quality control. The described methodologies provided a foundation to establish functional genomic studies of stallion fertility.

Following the recent development of high-resolution gene maps and generation of several basic tools and resources to use them in analyzing traits that are economically important to horse owners, genome analysis in horses is witnessing a shift towards developing an ability to analyze complex traits. The likelihood of this happening in the very near future is great, mainly because of the recent availability of the whole genome sequence in the horse. The latter has triggered the development of novel tools like SNP-chip and expression arrays that will permit rapid genome-wide analysis. While these tools will be used for a range of multi-factorial disease traits, attempts are underway to develop focused tools that can target reproduction, fertility and sex determination. For this, a catalog of sex and reproduction related (SRR) genes is being developed in horses. A recently developed dense map of the horse Y chromosome will provide genes that are expressed exclusively in males and, therefore, have an impact on stallion fertility. Overall, these advances in equine genome analysis hold promise for improved diagnosis and treatment of various conditions in horses.