The H2Bub1-deposition complex is required for human and mouse cardiogenesis.

De novo variants affecting monoubiquitylation of histone H2B (H2Bub1) are enriched in human congenital heart disease. H2Bub1 is required in stem cell differentiation, cilia function, post-natal cardiomyocyte maturation and transcriptional elongation. However, how H2Bub1 affects cardiogenesis is unknown. We show that the H2Bub1-deposition complex (RNF20-RNF40-UBE2B) is required for mouse cardiogenesis and for differentiation of human iPSCs into cardiomyocytes. Mice with cardiac-specific Rnf20 deletion are embryonic lethal and have abnormal myocardium. We then analyzed H2Bub1 marks during differentiation of human iPSCs into cardiomyocytes. H2Bub1 is erased from most genes at the transition from cardiac mesoderm to cardiac progenitor cells but is preserved on a subset of long cardiac-specific genes. When H2Bub1 is reduced in iPSC-derived cardiomyocytes, long cardiac-specific genes have fewer full-length transcripts. This correlates with H2Bub1 accumulation near the center of these genes. H2Bub1 accumulation near the center of tissue-specific genes was also observed in embryonic fibroblasts and fetal osteoblasts. In summary, we show that normal H2Bub1 distribution is required for cardiogenesis and cardiomyocyte differentiation, and suggest that H2Bub1 regulates tissue-specific gene expression by increasing the amount of full-length transcripts.

Molecular Genetics and Complex Inheritance of Congenital Heart Disease.

Congenital heart disease (CHD) is the most common congenital malformation and the leading cause of mortality therein. Genetic etiologies contribute to an estimated 90% of CHD cases, but so far, a molecular diagnosis remains unsolved in up to 55% of patients. Copy number variations and aneuploidy account for ~23% of cases overall, and high-throughput genomic technologies have revealed additional types of genetic variation in CHD. The first CHD risk genotypes identified through high-throughput sequencing were de novo mutations, many of which occur in chromatin modifying genes. Murine models of cardiogenesis further support the damaging nature of chromatin modifying CHD mutations. Transmitted mutations have also been identified through sequencing of population scale CHD cohorts, and many transmitted mutations are enriched in cilia genes and Notch or VEGF pathway genes. While we have come a long way in identifying the causes of CHD, more work is required to end the diagnostic odyssey for all CHD families. Complex genetic explanations of CHD are emerging but will require increasingly sophisticated analysis strategies applied to very large CHD cohorts before they can come to fruition in providing molecular diagnoses to genetically unsolved patients. In this review, we discuss the genetic architecture of CHD and biological pathways involved in its pathogenesis.

Histone H2B monoubiquitination regulates heart development via epigenetic control of cilia motility.

Genomic analyses of patients with congenital heart disease (CHD) have identified significant contribution from mutations affecting cilia genes and chromatin remodeling genes; however, the mechanism(s) connecting chromatin remodeling to CHD is unknown. Histone H2B monoubiquitination (H2Bub1) is catalyzed by the RNF20 complex consisting of RNF20, RNF40, and UBE2B. Here, we show significant enrichment of loss-of-function mutations affecting H2Bub1 in CHD patients (enrichment 6.01, P = 1.67 × 10-03), some of whom had abnormal laterality associated with ciliary dysfunction. In Xenopus, knockdown of rnf20 and rnf40 results in abnormal heart looping, defective development of left-right (LR) asymmetry, and impaired cilia motility. Rnf20, Rnf40, and Ube2b affect LR patterning and cilia synergistically. Examination of global H2Bub1 level in Xenopus embryos shows that H2Bub1 is developmentally regulated and requires Rnf20. To examine gene-specific H2Bub1, we performed ChIP-seq of mouse ciliated and nonciliated tissues and showed tissue-specific H2Bub1 marks significantly enriched at cilia genes including the transcription factor Rfx3 Rnf20 knockdown results in decreased levels of rfx3 mRNA in Xenopus, and exogenous rfx3 can rescue the Rnf20 depletion phenotype. These data suggest that Rnf20 functions at the Rfx3 locus regulating cilia motility and cardiac situs and identify H2Bub1 as an upstream transcriptional regulator controlling tissue-specific expression of cilia genes. Our findings mechanistically link the two functional gene ontologies that have been implicated in human CHD: chromatin remodeling and cilia function.

Evaluating optimal therapy robustness by virtual expansion of a sample population, with a case study in cancer immunotherapy.

Cancer is a highly heterogeneous disease, exhibiting spatial and temporal variations that pose challenges for designing robust therapies. Here, we propose the VEPART (Virtual Expansion of Populations for Analyzing Robustness of Therapies) technique as a platform that integrates experimental data, mathematical modeling, and statistical analyses for identifying robust optimal treatment protocols. VEPART begins with time course experimental data for a sample population, and a mathematical model fit to aggregate data from that sample population. Using nonparametric statistics, the sample population is amplified and used to create a large number of virtual populations. At the final step of VEPART, robustness is assessed by identifying and analyzing the optimal therapy (perhaps restricted to a set of clinically realizable protocols) across each virtual population. As proof of concept, we have applied the VEPART method to study the robustness of treatment response in a mouse model of melanoma subject to treatment with immunostimulatory oncolytic viruses and dendritic cell vaccines. Our analysis (i) showed that every scheduling variant of the experimentally used treatment protocol is fragile (nonrobust) and (ii) discovered an alternative region of dosing space (lower oncolytic virus dose, higher dendritic cell dose) for which a robust optimal protocol exists.

Novel determinants of mammalian primary microRNA processing revealed by systematic evaluation of hairpin-containing transcripts and human genetic variation.

Mature microRNAs (miRNAs) are processed from hairpin-containing primary miRNAs (pri-miRNAs). However, rules that distinguish pri-miRNAs from other hairpin-containing transcripts in the genome are incompletely understood. By developing a computational pipeline to systematically evaluate 30 structural and sequence features of mammalian RNA hairpins, we report several new rules that are preferentially utilized in miRNA hairpins and govern efficient pri-miRNA processing. We propose that a hairpin stem length of 36 ± 3 nt is optimal for pri-miRNA processing. We identify two bulge-depleted regions on the miRNA stem, located ∼16-21 nt and ∼28-32 nt from the base of the stem, that are less tolerant of unpaired bases. We further show that the CNNC primary sequence motif selectively enhances the processing of optimal-length hairpins. We predict that a small but significant fraction of human single-nucleotide polymorphisms (SNPs) alter pri-miRNA processing, and confirm several predictions experimentally including a disease-causing mutation. Our study enhances the rules governing mammalian pri-miRNA processing and suggests a diverse impact of human genetic variation on miRNA biogenesis.

Characterization of 16 microsatellite markers for the Oreinotinus clade of Viburnum (Adoxaceae).

PREMISE OF THE STUDY: Microsatellite loci were isolated from four species of Viburnum (Adoxaceae) to study population structure and assess species boundaries among morphologically similar South American Viburnum species of the Oreinotinus clade. METHODS AND RESULTS: Using a microsatellite-enriched library and mining next-generation sequence data, 16 microsatellites were developed. Each locus was tested on two populations of V. triphyllum and one population of V. pichinchense. For nuclear loci, one to 13 alleles were recovered, expected heterozygosity ranged from 0 to 0.8975, Simpson diversity index ranged from 0.0167 to 1.000, and Shannon diversity index ranged from 0 to 2.3670 in a given population. For the mitochondrial locus, three to six alleles were recovered and unbiased haploid diversity values ranged from 0.756 to 0.853 in a given population. CONCLUSIONS: The 16 microsatellite loci developed for the Oreinotinus clade (Viburnum, Adoxaceae) will inform investigations of population structure and species boundaries within this group.