Galpha(12/13) have been implicated in numerous cellular processes, however, their roles in vertebrate gastrulation are largely unknown. Here, we show that during zebrafish gastrulation, suppression of both Galpha(12) and Galpha(13) signaling by overexpressing dominant negative proteins and application of antisense morpholino-modified oligonucleotide translation interference disrupted convergence and extension without changing embryonic patterning. Analyses of mesodermal cell behaviors revealed that Galpha(12/13) are required for cell elongation and efficient dorsalward migration during convergence independent of noncanonical Wnt signaling. Furthermore, Galpha(12/13) function cell-autonomously to mediate mediolateral cell elongation underlying intercalation during notochord extension, likely acting in parallel to noncanonical Wnt signaling. These findings provide the first evidence that Galpha(12) and Galpha(13) have overlapping and essential roles in distinct cell behaviors that drive vertebrate gastrulation.

Nemaline myopathy (NM) is the most common congenital myopathy, characterized by extreme weakness of the respiratory, limb, and facial muscles. Pathogenic variants in Tropomyosin 2 (TPM2), which encodes a skeletal muscle-specific actin binding protein essential for sarcomere function, cause a spectrum of musculoskeletal disorders that include NM as well as cap myopathy, congenital fiber type disproportion, and distal arthrogryposis (DA). The in vivo pathomechanisms underlying TPM2-related disorders are unknown, so we expressed a series of dominant, pathogenic TPM2 variants in Drosophila embryos and found 4 variants significantly affected muscle development and muscle function. Transient overexpression of the 4 variants also disrupted the morphogenesis of mouse myotubes in vitro and negatively affected zebrafish muscle development in vivo. We used transient overexpression assays in zebrafish to characterize 2 potentially novel TPM2 variants and 1 recurring variant that we identified in patients with DA (V129A, E139K, A155T, respectively) and found these variants caused musculoskeletal defects similar to those of known pathogenic variants. The consistency of musculoskeletal phenotypes in our assays correlated with the severity of clinical phenotypes observed in our patients with DA, suggesting disrupted myogenesis is a potentially novel pathomechanism of TPM2 disorders and that our myogenic assays can predict the clinical severity of TPM2 variants.

Paroxysmal kinesigenic dyskinesia with infantile convulsions (PKD/IC) is an episodic movement disorder with autosomal-dominant inheritance and high penetrance, but the causative genetic mutation is unknown. We have now identified four truncating mutations involving the gene PRRT2 in the vast majority (24/25) of well-characterized families with PKD/IC. PRRT2 truncating mutations were also detected in 28 of 78 additional families. PRRT2 encodes a proline-rich transmembrane protein of unknown function that has been reported to interact with the t-SNARE, SNAP25. PRRT2 localizes to axons but not to dendritic processes in primary neuronal culture, and mutants associated with PKD/IC lead to dramatically reduced PRRT2 levels, leading ultimately to neuronal hyperexcitability that manifests in vivo as PKD/IC.

Embryonic axis formation in vertebrates is initiated by the establishment of the dorsal Nieuwkoop blastula organizer, marked by the nuclear accumulation of maternal β-catenin, a transcriptional effector of canonical Wnt signaling. Known regulators of axis specification include the canonical Wnt pathway components that positively or negatively affect β-catenin. An involvement of G-protein coupled receptors (GPCRs) was hypothesized from experiments implicating G proteins and intracellular calcium in axis formation, but such GPCRs have not been identified. Mobilization of intracellular Ca(2+) stores generates Ca(2+) transients in the superficial blastomeres of zebrafish blastulae when the nuclear accumulation of maternal β-catenin marks the formation of the Nieuwkoop organizer. Moreover, intracellular Ca(2+) downstream of non-canonical Wnt ligands was proposed to inhibit β-catenin and axis formation, but mechanisms remain unclear. Here we report a novel function of Ccr7 GPCR and its chemokine ligand Ccl19.1, previously implicated in chemotaxis and other responses of dendritic cells in mammals, as negative regulators of β-catenin and axis formation in zebrafish. We show that interference with the maternally and ubiquitously expressed zebrafish Ccr7 or Ccl19.1 expands the blastula organizer and the dorsoanterior tissues at the expense of the ventroposterior ones. Conversely, Ccr7 or Ccl19.1 overexpression limits axis formation. Epistatic analyses demonstrate that Ccr7 acts downstream of Ccl19.1 ligand and upstream of β-catenin transcriptional targets. Moreover, Ccl19/Ccr7 signaling reduces the level and nuclear accumulation of maternal β-catenin and its axis-inducing activity and can also inhibit the Gsk3β -insensitive form of β-catenin. Mutational and pharmacologic experiments reveal that Ccr7 functions during axis formation as a GPCR to inhibit β-catenin, likely by promoting Ca(2+) transients throughout the blastula. Our study delineates a novel negative, Gsk3β-independent control mechanism of β-catenin and implicates Ccr7 as a long-hypothesized GPCR regulating vertebrate axis formation.

Juvenile myoclonic epilepsy (JME), accounting for approximately 25% of idiopathic generalized epilepsies, is genetically heterogeneous. Mutations in the alpha-1 subunit of the GABAA receptor (GABRA1) and EFHC1 genes have been reported in a few families with autosomal dominant (AD) JME. We have investigated the contribution of these two genes to familial JME in our cohort of 54 JME Caucasian families. Syndromic classification of JME was based on previously published criteria. We considered kindreds with at least one affected first-degree relative and the evidence of a vertical transmission as definite AD JME, and families with at least one affected second-degree relative as probable AD JME. We included 33 families meeting criteria for definitive AD JME and 21 that were classified as probable AD JME. None of these families were considered informative enough to analyze candidate loci for JME using linkage analysis. We have systematically screened coding exons of these two genes using temperature gradient capillary electrophoresis. Every heteroduplex with an abnormal mobility was sequenced. No disease-causing mutations in the GABRA1 gene were identified. Analysis of EFHC1 gene found one putative disease-causing mutation R221H that was previously reported as a tandem mutation. Several synonymous and non-synonymous coding polymorphisms were identified but the allelic frequency did not differ between controls and affected individuals. Our data suggests that the majority of familial AD JME is not caused by mutations in the GABRA1 and EFHC1 genes.

Adolescents with low von Willebrand factor (VWF) levels and heavy menstrual bleeding (HMB) experience significant morbidity. There is a need to better characterize these patients genetically and improve our understanding of the pathophysiology of bleeding. We performed whole-exome sequencing on 86 postmenarchal patients diagnosed with low VWF levels (30-50 IU/dL) and HMB and compared them with 660 in-house controls. We compared the number of rare stop-gain/stop-loss and rare ClinVar "pathogenic" variants between cases and controls, as well as performed gene burden and gene-set burden analyses. We found an enrichment in cases of rare stop-gain/stop-loss variants in genes involved in bleeding disorders and an enrichment of rare ClinVar "pathogenic" variants in genes involved in anemias. The 2 most significant genes in the gene burden analysis, CFB and DNASE2, are associated with atypical hemolytic uremia and severe anemia, respectively. VWF also surpassed exome-wide significance in the gene burden analysis (P = 7.31 × 10-6). Gene-set burden analysis revealed an enrichment of rare nonsynonymous variants in cases in several hematologically relevant pathways. Further, common variants in FERMT2, a gene involved in the regulation of hemostasis and angiogenesis, surpassed genome-wide significance. We demonstrate that adolescents with HMB and low VWF have an excess of rare nonsynonymous and pathogenic variants in genes involved in bleeding disorders and anemia. Variants of variable penetrance in these genes may contribute to the spectrum of phenotypes observed in patients with HMB and could partially explain the bleeding phenotype. By identifying patients with HMB who possess these variants, we may be able to improve risk stratification and patient outcomes.

Non-obstructive azoospermia (NOA) is the most severe form of male infertility and typically incurable. Defining the genetic basis of NOA has proven challenging, and the most advanced classification of NOA subforms is not based on genetics, but simple description of testis histology. In this study, we exome-sequenced over 1000 clinically diagnosed NOA cases and identified a plausible recessive Mendelian cause in 20%. We find further support for 21 genes in a 2-stage burden test with 2072 cases and 11,587 fertile controls. The disrupted genes are primarily on the autosomes, enriched for undescribed human "knockouts", and, for the most part, have yet to be linked to a Mendelian trait. Integration with single-cell RNA sequencing data shows that azoospermia genes can be grouped into molecular subforms with synchronized expression patterns, and analogs of these subforms exist in mice. This analysis framework identifies groups of genes with known roles in spermatogenesis but also reveals unrecognized subforms, such as a set of genes expressed across mitotic divisions of differentiating spermatogonia. Our findings highlight NOA as an understudied Mendelian disorder and provide a conceptual structure for organizing the complex genetics of male infertility, which may provide a rational basis for disease classification.

Chiari I malformation (CM1), the displacement of the cerebellum through the foramen magnum into the spinal canal, is one of the most common pediatric neurological conditions. Individuals with CM1 can present with neurological symptoms, including severe headaches and sensory or motor deficits, often as a consequence of brainstem compression or syringomyelia (SM). We conducted whole-exome sequencing (WES) on 668 CM1 probands and 232 family members and performed gene-burden and de novo enrichment analyses. A significant enrichment of rare and de novo non-synonymous variants in chromodomain (CHD) genes was observed among individuals with CM1 (combined p = 2.4 × 10-10), including 3 de novo loss-of-function variants in CHD8 (LOF enrichment p = 1.9 × 10-10) and a significant burden of rare transmitted variants in CHD3 (p = 1.8 × 10-6). Overall, individuals with CM1 were found to have significantly increased head circumference (p = 2.6 × 10-9), with many harboring CHD rare variants having macrocephaly. Finally, haploinsufficiency for chd8 in zebrafish led to macrocephaly and posterior hindbrain displacement reminiscent of CM1. These results implicate chromodomain genes and excessive brain growth in CM1 pathogenesis.

Genetic factors predictive of severe adolescent idiopathic scoliosis (AIS) are largely unknown. To identify genetic variation associated with severe AIS, we performed an exome-wide association study of 457 severe AIS cases and 987 controls. We find a missense SNP in SLC39A8 (p.Ala391Thr, rs13107325) associated with severe AIS (P = 1.60 × 10-7, OR = 2.01, CI = 1.54-2.62). This pleiotropic SNP was previously associated with BMI, blood pressure, cholesterol, and blood manganese level. We replicate the association in a second cohort (841 cases and 1095 controls) resulting in a combined P = 7.02 × 10-14, OR = 1.94, CI = 1.63-2.34. Clinically, the minor allele of rs13107325 is associated with greater spinal curvature, decreased height, increased BMI and lower plasma manganese in our AIS cohort. Functional studies demonstrate reduced manganese influx mediated by the SLC39A8 p.Ala391Thr variant and vertebral abnormalities, impaired growth, and decreased motor activity in slc39a8 mutant zebrafish. Our results suggest the possibility that scoliosis may be amenable to dietary intervention.

Despite declining sequencing costs, few methods are available for cost-effective single-nucleotide polymorphism (SNP), insertion/deletion (INDEL) and copy number variation (CNV) discovery in a single assay. Commercially available methods require a high investment to a specific region and are only cost-effective for large samples. Here, we introduce a novel, flexible approach for multiplexed targeted sequencing and CNV analysis of large genomic regions called multiplexed direct genomic selection (MDiGS). MDiGS combines biotinylated bacterial artificial chromosome (BAC) capture and multiplexed pooled capture for SNP/INDEL and CNV detection of 96 multiplexed samples on a single MiSeq run. MDiGS is advantageous over other methods for CNV detection because pooled sample capture and hybridization to large contiguous BAC baits reduces sample and probe hybridization variability inherent in other methods. We performed MDiGS capture for three chromosomal regions consisting of ∼ 550 kb of coding and non-coding sequence with DNA from 253 patients with congenital lower limb disorders. PITX1 nonsense and HOXC11 S191F missense mutations were identified that segregate in clubfoot families. Using a novel pooled-capture reference strategy, we identified recurrent chromosome chr17q23.1q23.2 duplications and small HOXC 5' cluster deletions (51 kb and 12 kb). Given the current interest in coding and non-coding variants in human disease, MDiGS fulfills a niche for comprehensive and low-cost evaluation of CNVs, coding, and non-coding variants across candidate regions of interest.

Idiopathic scoliosis occurs in 3% of individuals and has an unknown etiology. The objective of this study was to identify rare variants that contribute to the etiology of idiopathic scoliosis by using exome sequencing in a multigenerational family with idiopathic scoliosis. Exome sequencing was completed for three members of this multigenerational family with idiopathic scoliosis, resulting in the identification of a variant in the HSPG2 gene as a potential contributor to the phenotype. The HSPG2 gene was sequenced in a separate cohort of 100 unrelated individuals affected with idiopathic scoliosis and also was examined in an independent idiopathic scoliosis population. The exome sequencing and subsequent bioinformatics filtering resulted in 16 potentially damaging and rare coding variants. One of these variants, p.Asn786Ser, is located in the HSPG2 gene. The variant p.Asn786Ser also is overrepresented in a larger cohort of idiopathic scoliosis cases compared with a control population (P = 0.024). Furthermore, we identified additional rare HSPG2 variants that are predicted to be damaging in two independent cohorts of individuals with idiopathic scoliosis. The HSPG2 gene encodes for a ubiquitous multifunctional protein within the extracellular matrix in which loss of function mutation are known to result in a musculoskeletal phenotype in both mouse and humans. Based on these results, we conclude that rare variants in the HSPG2 gene potentially contribute to the idiopathic scoliosis phenotype in a subset of patients with idiopathic scoliosis. Further studies must be completed to confirm the effect of the HSPG2 gene on the idiopathic scoliosis phenotype.

Although autism spectrum disorders (ASD) are among the most heritable of all neuropsychiatric syndromes, most affected children are born to unaffected parents. Recently, we reported an average increase of 3-5% over general population risk of ASD among offspring of adults who have first-degree relatives with ASD in a large epidemiologic family sample. A next essential step is to investigate whether there are measurable characteristics of individual parents placing them at higher or lower recurrence risk, as this information could allow more personalized genetic counseling.

Transmission of SARS-CoV-2 in schools primarily for typically developing children is rare. However, less is known about transmission in schools for children with intellectual and developmental disabilities (IDD), who are often unable to mask or maintain social distancing. The objectives of this study were to determine SARS-CoV-2 positivity and in-school transmission rates using weekly screening tests for school staff and students and describe the concurrent deployment of mitigation strategies in six schools for children with IDD.

Cell proliferation has generally been considered dispensable for anteroposterior extension of embryonic axis during vertebrate gastrulation. Signal transducer and activator of transcription 3 (Stat3), a conserved controller of cell proliferation, survival and regeneration, is associated with human scoliosis, cancer and Hyper IgE Syndrome. Zebrafish Stat3 was proposed to govern convergence and extension gastrulation movements in part by promoting Wnt/Planar Cell Polarity (PCP) signaling, a conserved regulator of mediolaterally polarized cell behaviors. Here, using zebrafish stat3 null mutants and pharmacological tools, we demonstrate that cell proliferation contributes to anteroposterior embryonic axis extension. Zebrafish embryos lacking maternal and zygotic Stat3 expression exhibit normal convergence movements and planar cell polarity signaling, but transient axis elongation defect due to insufficient number of cells resulting largely from reduced cell proliferation and increased apoptosis. Pharmacologic inhibition of cell proliferation during gastrulation phenocopied axis elongation defects. Stat3 regulates cell proliferation and axis extension in part via upregulation of Cdc25a expression during oogenesis. Accordingly, restoring Cdc25a expression in stat3 mutants partially suppressed cell proliferation and gastrulation defects. During later development, stat3 mutant zebrafish exhibit stunted growth, scoliosis, excessive inflammation, and fail to thrive, affording a genetic tool to study Stat3 function in vertebrate development, regeneration, and disease.

A recurrent de novo mutation in the transcriptional corepressor CTBP1 is associated with neurodevelopmental disabilities in children (Beck et al., 2016, 2019; Sommerville et al., 2017). All reported patients harbor a single recurrent de novo heterozygous missense mutation (p.R342W) within the cofactor recruitment domain of CtBP1. To investigate the transcriptional activity of the pathogenic CTBP1 mutant allele in physiologically relevant human cell models, we generated induced pluripotent stem cells (iPSC) from the dermal fibroblasts derived from patients and normal donors. The transcriptional profiles of the iPSC-derived "early" neurons were determined by RNA-sequencing. Comparison of the RNA-seq data of the neurons from patients and normal donors revealed down regulation of gene networks involved in neurodevelopment, synaptic adhesion and anti-viral (interferon) response. Consistent with the altered gene expression patterns, the patient-derived neurons exhibited morphological and electrophysiological abnormalities, and susceptibility to viral infection. Taken together, our studies using iPSC-derived neuron models provide novel insights into the pathological activities of the CTBP1 p.R342W allele.

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

The goal of this study is to demonstrate the utility of a growth assay to quantify the functional impact of single nucleotide variants (SNVs) in SLC2A1, the gene responsible for Glut1DS.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

In cancer, missense mutations in the DNA-binding domain of TP53 are common. They abrogate canonical p53 activity and frequently confer gain-of-oncogenic function (GOF) through localization of transcriptionally active mutant p53 to non-canonical genes. We found that several recurring p53 mutations exhibit a sex difference in frequency in patients with glioblastoma (GBM). In vitro and in vivo analysis of three mutations, p53R172H, p53Y202C, and p53Y217C revealed unique interactions between cellular sex and p53 GOF mutations that determined each mutation's ability to transform male versus female primary mouse astrocytes. These phenotypic differences were correlated with sex- and p53 mutation- specific patterns of genomic localization to the transcriptional start sites of upregulated genes belonging to core cancer pathways. The promoter regions of these genes exhibited a sex difference in enrichment for different transcription factor DNA-binding motifs. Together, our data establish a novel mechanism for sex specific mutant p53 GOF activity in GBM with implications for all cancer.

Congenital clubfoot is a common birth defect that affects at least 0.1% of all births. Nearly 25% cases are familial and the remaining are sporadic in inheritance. Copy number variants (CNVs) involving transcriptional regulators of limb development, including PITX1 and TBX4, have previously been shown to cause familial clubfoot, but much of the heritability remains unexplained.