Pleiotropy has long been thought to be a common phenomenon in the human genome; however, until recently appropriate data was unavailable to test this hypothesis. Prior studies focused on assessing the prevalence of pleiotropy in only small subsets of phenotypes (≤ 53 phenotypes), without a truly comprehensive assessment of pleiotropy in the human genome. In this study, we determined the prevalence of pleiotropy, using the entire GWAS catalog (1094 disease phenotypes, 14,459 genes), as well as investigate the relationship between the degree of pleiotropy and the average effect size for each associating gene. The number of associating phenotypes per gene ranged from 1 to 53, with 44% of genes reported in the GWAS catalog associating with more than one phenotype. The proportion of genes shown to be pleiotropic has continued to increase as more studies are added to the catalog. We also found the degree of pleiotropy scales positively with a gene's average effect size (r = 0.04, p value = 0.0003) and negatively with the variance of effect sizes in genes with a given number of associating phenotypes (r = − 0.590, p value = 0.0019). Based on this and prior work, it is becoming evident that pleiotropy is a common, if not ubiquitous, phenomenon. These results have implications in understanding disease etiologies, potentially common biology underlying even disparate diseases, and in elucidating the genotype-phenotype map.

In Duchenne muscular dystrophy, dystrophin loss leads to chronic muscle damage, dysregulation of repair, fibro-fatty replacement, and weakness. We develop methodology to efficiently isolate individual nuclei from minute quantities of frozen skeletal muscle, allowing single nuclei sequencing of irreplaceable archival samples and from very small samples. We apply this method to identify cell and gene expression dynamics within human DMD and mdx mouse muscle, characterizing effects of dystrophin rescue by exon skipping therapy at single nuclei resolution. DMD exon 23 skipping events are directly observed and increased in myonuclei from treated mice. We describe partial rescue of type IIa and IIx myofibers, expansion of an MDSC-like myeloid population, recovery of repair/remodeling M2-macrophage, and repression of inflammatory POSTN1 + fibroblasts in response to exon skipping and partial dystrophin restoration. Use of this method enables exploration of cellular and transcriptomic mechanisms of dystrophin loss and repair within an intact muscle environment. Our initial findings will scaffold our future work to more directly examine muscular dystrophies and putative recovery pathways.

While the American horseshoe crab, Limulus polyphemus, has robust circadian and circatidal rhythms, virtually nothing is known about the molecular basis of these rhythms in this species or any other chelicerate. In this study, next generation sequencing was used to assemble transcriptomic reads and then putative homologs of known core and accessory circadian genes were identified in these databases. Homologous transcripts were discovered for one circadian clock input gene, five core genes, 22 accessory genes, and two possible output pathways. Alignments and functional domain analyses showed generally high conservation between the putative L. polyphemus clock genes and homologs from Drosophila melanogaster and Daphnia pulex. The presence of both cry1 and cry2 in the L. polyphemus transcriptome would classify its system as an "ancestral", type 2 clock system. In addition, a novel duplication of CYCLE, and a novel triplication of PERIOD were found. Investigations are currently underway to determine if any of these "circadian" genes also participate in the molecular processes that drive the Limulus circatidal clock.

Pleiotropy has been claimed to constrain gene evolution but specific mechanisms and extent of these constraints have been difficult to demonstrate. The expansion of molecular data makes it possible to investigate these pleiotropic effects. Few classes of genes have been characterized as intensely as human transcription factors (TFs). We therefore analyzed the evolutionary rates of full TF proteins, along with their DNA binding domains and protein-protein interacting domains (PID) in light of the degree of pleiotropy, measured by the number of TF–TF interactions, or the number of DNA-binding targets. Data were extracted from the ENCODE Chip-Seq dataset, the String v 9.2 database, and the NHGRI GWAS catalog. Evolutionary rates of proteins and domains were calculated using the PAML CodeML package. Our analysis shows that the numbers of TF-TF interactions and DNA binding targets associated with constrained gene evolution; however, the constraint caused by the number of DNA binding targets was restricted to the DNA binding domains, whereas the number of TF-TF interactions constrained the full protein and did so more strongly. Additionally, we found a positive correlation between the number of protein–PIDs and the evolutionary rates of the protein–PIDs. These findings show that not only does pleiotropy associate with constrained protein evolution but the constraint differs by domain function. Finally, we show that GWAS associated TF genes are more highly pleiotropic. The GWAS data illustrates that mutations in highly pleiotropic genes are more likely to be associated with disease phenotypes.

ObjectivePremature aging and short telomere lengths of fetal tissues are associated with spontaneous preterm labor (PTL) and preterm premature rupture of membranes (pPROM). Maintenance of telomere length is performed by the enzyme telomerase. Human telomerase reverse transcriptase (hTERT) is a subunit of telomerase, and its dysfunction affects telomere shortening. This study assessed whether maternal or fetal genetic variations in the hTERT gene are associated with PTL or pPROM.MethodsA case (PTL or pPROM) control (term birth) genetic association study was conducted in 654 non-Hispanic white mothers (438 term, 162 PTL, 54 pPROM) and 502 non-Hispanic white newborns (346 term, 116 PTB, 40 pPROM). Maternal and fetal DNA samples were genotyped for 23 single nucleotide polymorphisms (SNPs) within the hTERT gene. Allele frequencies were compared between cases and controls, stratified by PTL and pPROM. Maternal and fetal data were analyzed separately.ResultsAllelic differences in one SNP of hTERT (rs2853690) were significantly associated with both PTL (adjusted OR 2.24, 95%CI 1.64–3.06, p = 2.32e-05) and with pPROM (adjusted OR 7.54, 95%CI 3.96–14.33, p = 2.39e-07) in maternal DNA. There was no significant association between the hTERT SNPs analyzed and PTL or pPROM in the fetal samples.ConclusionhTERT polymorphisms in fetal DNA do not associate with PTL or pPROM risk; however, maternal genetic variations in hTERT may play a contributory role in risk of PTL and PPROM.