Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

Hepatitis C virus promotes virion secretion through cleavage of the Rab7 adaptor protein RILP.

Hepatitis C virus (HCV) is an enveloped RNA virus that modifies intracellular trafficking processes. The mechanisms that HCV and other viruses use to modify these events are poorly understood. In this study, we observed that two different RNA viruses, HCV and Sendai, cause inhibition of ras-related protein Rab-7 (Rab7)-dependent endosome-lysosome fusion. In both cases, viral infection causes cleavage of the Rab7 adaptor protein RILP (Rab interacting lysosomal protein), which is responsible for linking Rab7 vesicles to dynein motor complexes. RILP cleavage results in the generation of a cleaved RILP fragment (cRILP) missing the N terminus of the molecule. Although RILP localizes in a perinuclear fashion, cRILP moves to the cell periphery. Both knockdown of RILP and expression of cRILP reproduced the HCV-induced trafficking defect, and restoring full-length RILP reversed the trafficking effects of virus. For the first 3 d after electroporation of HCV RNA, intracellular virus predominates over secreted virus, but the quantity of intracellular virus then rapidly declines as secreted virus dominates. The transition from the intracellular-predominant to the secretion-predominant phenotype corresponds to the time course of cRILP generation. Expressing cRILP directly prevents intracellular virus accumulation at early times without affecting net virus production. The ability of cRILP to promote virus secretion could be prevented by a kinesin inhibitor. HCV thus modifies cellular trafficking by cleaving RILP, which serves to redirect Rab7-containing vesicles to a kinesin-dependent trafficking mode promoting virion secretion. Cleavage of a Rab adaptor protein is thus a mechanism by which viruses modify trafficking patterns of infected cells.

Trans-splicing and operons in metazoans: translational control in maternally regulated development and recovery from growth arrest.

Polycistronic mRNAs transcribed from operons are resolved via the trans-splicing of a spliced-leader (SL) RNA. Trans-splicing also occurs at monocistronic transcripts. The phlyogenetically sporadic appearance of trans-splicing and operons has made the driving force(s) for their evolution in metazoans unclear. Previous work has proposed that germline expression drives operon organization in Caenorhabditis elegans, and a recent hypothesis proposes that operons provide an evolutionary advantage via the conservation of transcriptional machinery during recovery from growth arrested states. Using a modified cap analysis of gene expression protocol we mapped sites of SL trans-splicing genome-wide in the marine chordate Oikopleura dioica. Tiled microarrays revealed the expression dynamics of trans-spliced genes across development and during recovery from growth arrest. Operons did not facilitate recovery from growth arrest in O. dioica. Instead, we found that trans-spliced transcripts were predominantly maternal. We then analyzed data from C. elegans and Ciona intestinalis and found that an enrichment of trans-splicing and operon gene expression in maternal mRNA is shared between all three species, suggesting that this may be a driving force for operon evolution in metazoans. Furthermore, we found that the majority of known terminal oligopyrimidine (TOP) mRNAs are trans-spliced in O. dioica and that the SL contains a TOP-like motif. This suggests that the SL in O. dioica confers nutrient-dependent translational control to trans-spliced mRNAs via the TOR-signaling pathway. We hypothesize that SL-trans-splicing provides an evolutionary advantage in species that depend on translational control for regulating early embryogenesis, growth and oocyte production in response to nutrient levels.

Copy number variation analysis implicates the cell polarity gene glypican 5 as a human spina bifida candidate gene.

Neural tube defects (NTDs) are common birth defects of complex etiology. Family and population-based studies have confirmed a genetic component to NTDs. However, despite more than three decades of research, the genes involved in human NTDs remain largely unknown. We tested the hypothesis that rare copy number variants (CNVs), especially de novo germline CNVs, are a significant risk factor for NTDs. We used array-based comparative genomic hybridization (aCGH) to identify rare CNVs in 128 Caucasian and 61 Hispanic patients with non-syndromic lumbar-sacral myelomeningocele. We also performed aCGH analysis on the parents of affected individuals with rare CNVs where parental DNA was available (42 sets). Among the eight de novo CNVs that we identified, three generated copy number changes of entire genes. One large heterozygous deletion removed 27 genes, including PAX3, a known spina bifida-associated gene. A second CNV altered genes (PGPD8, ZC3H6) for which little is known regarding function or expression. A third heterozygous deletion removed GPC5 and part of GPC6, genes encoding glypicans. Glypicans are proteoglycans that modulate the activity of morphogens such as Sonic Hedgehog (SHH) and bone morphogenetic proteins (BMPs), both of which have been implicated in NTDs. Additionally, glypicans function in the planar cell polarity (PCP) pathway, and several PCP genes have been associated with NTDs. Here, we show that GPC5 orthologs are expressed in the neural tube, and that inhibiting their expression in frog and fish embryos results in NTDs. These results implicate GPC5 as a gene required for normal neural tube development.

The Dm-Myb Oncoprotein Contributes to Insulator Function and Stabilizes Repressive H3K27me3 PcG Domains.

Drosophila Myb (Dm-Myb) encodes a protein that plays a key role in regulation of mitotic phase genes. Here, we further refine its role in the context of a developing tissue as a potentiator of gene expression required for proper RNA polymerase II (RNA Pol II) function and efficient H3K4 methylation at promoters. In contrast to its role in gene activation, Myb is also required for repression of many genes, although no specific mechanism for this role has been proposed. We now reveal a critical role for Myb in contributing to insulator function, in part by promoting binding of insulator proteins BEAF-32 and CP190 and stabilizing H3K27me3 Polycomb-group (PcG) domains. In the absence of Myb, H3K27me3 is markedly reduced throughout the genome, leading to H3K4me3 spreading and gene derepression. Finally, Myb is enriched at boundaries that demarcate chromatin environments, including chromatin loop anchors. These results reveal functions of Myb that extend beyond transcriptional regulation.

Identification of Isthmin 1 as a Novel Clefting and Craniofacial Patterning Gene in Humans.

Orofacial clefts are one of the most common birth defects, affecting 1-2 per 1000 births, and have a complex etiology. High-resolution array-based comparative genomic hybridization has increased the ability to detect copy number variants (CNVs) that can be causative for complex diseases such as cleft lip and/or palate. Utilizing this technique on 97 nonsyndromic cleft lip and palate cases and 43 cases with cleft palate only, we identified a heterozygous deletion of Isthmin 1 in one affected case, as well as a deletion in a second case that removes putative 3' regulatory information. Isthmin 1 is a strong candidate for clefting, as it is expressed in orofacial structures derived from the first branchial arch and is also in the same "synexpression group" as fibroblast growth factor 8 and sprouty RTK signaling antagonist 1a and 2, all of which have been associated with clefting. CNVs affecting Isthmin 1 are exceedingly rare in control populations, and Isthmin 1 scores as a likely haploinsufficiency locus. Confirming its role in craniofacial development, knockdown or clustered randomly interspaced short palindromic repeats/Cas9-generated mutation of isthmin 1 in Xenopus laevis resulted in mild to severe craniofacial dysmorphologies, with several individuals presenting with median clefts. Moreover, knockdown of isthmin 1 produced decreased expression of LIM homeobox 8, itself a gene associated with clefting, in regions of the face that pattern the maxilla. Our study demonstrates a successful pipeline from CNV identification of a candidate gene to functional validation in a vertebrate model system, and reveals Isthmin 1 as both a new human clefting locus as well as a key craniofacial patterning gene.

Sex-specific chromatin landscapes in an ultra-compact chordate genome.

BACKGROUND: In multicellular organisms, epigenome dynamics are associated with transitions in the cell cycle, development, germline specification, gametogenesis and inheritance. Evolutionarily, regulatory space has increased in complex metazoans to accommodate these functions. In tunicates, the sister lineage to vertebrates, we examine epigenome adaptations to strong secondary genome compaction, sex chromosome evolution and cell cycle modes. RESULTS: Across the 70 MB Oikopleura dioica genome, we profiled 19 histone modifications, and RNA polymerase II, CTCF and p300 occupancies, to define chromatin states within two homogeneous tissues with distinct cell cycle modes: ovarian endocycling nurse nuclei and mitotically proliferating germ nuclei in testes. Nurse nuclei had active chromatin states similar to other metazoan epigenomes, with large domains of operon-associated transcription, a general lack of heterochromatin, and a possible role of Polycomb PRC2 in dosage compensation. Testis chromatin states reflected transcriptional activity linked to spermatogenesis and epigenetic marks that have been associated with establishment of transgenerational inheritance in other organisms. We also uncovered an unusual chromatin state specific to the Y-chromosome, which combined active and heterochromatic histone modifications on specific transposable elements classes, perhaps involved in regulating their activity. CONCLUSIONS: Compacted regulatory space in this tunicate genome is accompanied by reduced heterochromatin and chromatin state domain widths. Enhancers, promoters and protein-coding genes have conserved epigenomic features, with adaptations to the organization of a proportion of genes in operon units. We further identified features specific to sex chromosomes, cell cycle modes, germline identity and dosage compensation, and unusual combinations of histone PTMs with opposing consensus functions.