Prdm13 mediates the balance of inhibitory and excitatory neurons in somatosensory circuits

Epistasis analysis for quantitative traits by functional regression model

The critical barrier in interaction analysis for rare variants is that most traditional statistical methods for testing interactions were originally designed for testing the interaction between common variants and are difficult to apply to rare variants because of their prohibitive computational time and poor ability. The great challenges for successful detection of interactions with next-generation sequencing (NGS) data are (1) lack of methods for interaction analysis with rare variants, (2) severe multiple testing, and (3) time-consuming computations. To meet these challenges, we shift the paradigm of interaction analysis between two loci to interaction analysis between two sets of loci or genomic regions and collectively test interactions between all possible pairs of SNPs within two genomic regions. In other words, we take a genome region as a basic unit of interaction analysis and use high-dimensional data reduction and functional data analysis techniques to develop a novel functional regression model to collectively test interactions between all possible pairs of single nucleotide polymorphisms (SNPs) within two genome regions. By intensive simulations, we demonstrate that the functional regression models for interaction analysis of the quantitative trait have the correct type 1 error rates and a much better ability to detect interactions than the current pairwise interaction analysis. The proposed method was applied to exome sequence data from the NHLBI’s Exome Sequencing Project (ESP) and CHARGE-S study. We discovered 27 pairs of genes showing significant interactions after applying the Bonferroni correction (P-values < 4.58 × 10⁻¹⁰) in the ESP, and 11 were replicated in the CHARGE-S study.

prdm12b specifies the p1 progenitor domain and reveals a role for V1 interneurons in swim movements

Proper functioning of the vertebrate central nervous system requires the precise positioning of many neuronal cell types. This positioning is established during early embryogenesis when gene regulatory networks pattern the neural tube along its anteroposterior and dorsoventral axes. Dorsoventral patterning of the embryonic neural tube gives rise to multiple progenitor cell domains that go on to differentiate unique classes of neurons and glia. While the genetic program is reasonably well understood for some lineages, such as ventrally derived motor neurons and glia, other lineages are much less characterized. Here we show that prdm12b, a member of the PR domain containing-family of transcriptional regulators, is expressed in the p1 progenitor domain of the zebrafish neural tube in response to Sonic Hedgehog signaling. We find that disruption of prdm12b function leads to dorsal expansion of nkx6.1 expression and loss of p1-derived eng1b-expressing V1 interneurons, while the adjacent p0 and p2 domains are unaffected. We also demonstrate that prdm12b-deficient fish exhibit an abnormal touch-evoked escape response with excessive body contractions and a prolonged response time, as well as an inability to coordinate swimming movements, thereby revealing a functional role for V1 interneurons in locomotor circuits. We conclude that prdm12b is required for V1 interneuron specification and that these neurons control swimming movements in zebrafish.

The Prdm13 histone methyltransferase encoding gene is a Ptf1a-Rbpj downstream target that suppresses glutamatergic and promotes GABAergic neuronal fate in the dorsal neural tube

The basic helix-loop-helix (bHLH) transcriptional activator Ptf1a determines inhibitory GABAergic over excitatory glutamatergic neuronal cell fate in progenitors of the vertebrate dorsal spinal cord, cerebellum and retina. In an in situ hybridization expression survey of PR domain containing genes encoding putative chromatin-remodeling zinc finger transcription factors in Xenopus embryos, we identified Prdm13 as a histone methyltransferase belonging to the Ptf1a synexpression group. Gain and loss of Ptf1a function analyses in both frog and mice indicates that Prdm13 is positively regulated by Ptf1a and likely constitutes a direct transcriptional target. We also showed that this regulation requires the formation of the Ptf1a-Rbp-j complex. Prdm13 knockdown in Xenopus embryos and in Ptf1a overexpressing ectodermal explants lead to an upregulation of Tlx3/Hox11L2, which specifies a glutamatergic lineage and a reduction of the GABAergic neuronal marker Pax2. It also leads to an upregulation of Prdm13 transcription, suggesting an autonegative regulation. Conversely, in animal caps, Prdm13 blocks the ability of the bHLH factor Neurog2 to activate Tlx3. Additional gain of function experiments in the chick neural tube confirm that Prdm13 suppresses Tlx3(+)/glutamatergic and induces Pax2(+)/GABAergic neuronal fate. Thus, Prdm13 is a novel crucial component of the Ptf1a regulatory pathway that, by modulating the transcriptional activity of bHLH factors such as Neurog2, controls the balance between GABAergic and glutamatergic neuronal fate in the dorsal and caudal part of the vertebrate neural tube.