Loss of RNA expression and allele-specific expression associated with congenital heart disease

Increased endothelial sclerostin caused by elevated DSCAM mediates multiple trisomy 21 phenotypes

Trisomy 21 (T21), a recurrent aneuploidy occurring in 1:800 births, predisposes to congenital heart disease (CHD) and multiple extracardiac phenotypes. Despite a definitive genetic etiology, the mechanisms by which T21 perturbs development and homeostasis remain poorly understood. We compared the transcriptome of CHD tissues from 49 patients with T21 and 226 with euploid CHD (eCHD). We resolved cell lineages that misexpressed T21 transcripts by cardiac single-nucleus RNA sequencing and RNA in situ hybridization. Compared with eCHD samples, T21 samples had increased chr21 gene expression; 11-fold-greater levels (P = 1.2 × 10-8) of SOST (chr17), encoding the Wnt inhibitor sclerostin; and 1.4-fold-higher levels (P = 8.7 × 10-8) of the SOST transcriptional activator ZNF467 (chr7). Euploid and T21 cardiac endothelial cells coexpressed SOST and ZNF467; however, T21 endothelial cells expressed 6.9-fold more SOST than euploid endothelial cells (P = 2.7 × 10-27). Wnt pathway genes were downregulated in T21 endothelial cells. Expression of DSCAM, residing within the chr21 CHD critical region, correlated with SOST (P = 1.9 × 10-5) and ZNF467 (P = 2.9 × 10-4). Deletion of DSCAM from T21 endothelial cells derived from human induced pluripotent stem cells diminished sclerostin secretion. As Wnt signaling is critical for atrioventricular canal formation, bone health, and pulmonary vascular homeostasis, we concluded that T21-mediated increased sclerostin levels would inappropriately inhibit Wnt activities and promote Down syndrome phenotypes. These findings imply therapeutic potential for anti-sclerostin antibodies in T21.

In-Depth Genomic Analysis: The New Challenge in Congenital Heart Disease

The use of next-generation sequencing has provided new insights into the causes and mechanisms of congenital heart disease (CHD). Examinations of the whole exome sequence have detected detrimental gene variations modifying single or contiguous nucleotides, which are characterised as pathogenic based on statistical assessments of families and correlations with congenital heart disease, elevated expression during heart development, and reductions in harmful protein-coding mutations in the general population. Patients with CHD and extracardiac abnormalities are enriched for gene classes meeting these criteria, supporting a common set of pathways in the organogenesis of CHDs. Single-cell transcriptomics data have revealed the expression of genes associated with CHD in specific cell types, and emerging evidence suggests that genetic mutations disrupt multicellular genes essential for cardiogenesis. Metrics and units are being tracked in whole-genome sequencing studies.

PCBP1 regulates alternative splicing of AARS2 in congenital cardiomyopathy

Alanyl-transfer RNA synthetase 2 (AARS2) is a nuclear encoded mitochondrial tRNA synthetase that is responsible for charging of tRNA-Ala with alanine during mitochondrial translation. Homozygous or compound heterozygous mutations in the Aars2 gene, including those affecting its splicing, are linked to infantile cardiomyopathy in humans. However, how Aars2 regulates heart development, and the underlying molecular mechanism of heart disease remains unknown. Here, we found that poly(rC) binding protein 1 (PCBP1) interacts with the Aars2 transcript to mediate its alternative splicing and is critical for the expression and function of Aars2. Cardiomyocyte-specific deletion of Pcbp1 in mice resulted in defects in heart development that are reminiscent of human congenital cardiac defects, including noncompaction cardiomyopathy and a disruption of the cardiomyocyte maturation trajectory. Loss of Pcbp1 led to an aberrant alternative splicing and a premature termination of Aars2 in cardiomyocytes. Additionally, Aars2 mutant mice with exon-16 skipping recapitulated heart developmental defects observed in Pcbp1 mutant mice. Mechanistically, we found dysregulated gene and protein expression of the oxidative phosphorylation pathway in both Pcbp1 and Aars2 mutant hearts; these date provide further evidence that the infantile hypertrophic cardiomyopathy associated with the disorder oxidative phosphorylation defect type 8 (COXPD8) is mediated by Aars2. Our study therefore identifies Pcbp1 and Aars2 as critical regulators of heart development and provides important molecular insights into the role of disruptions in metabolism on congenital heart defects.

Risk of Congenital Heart Defects in Offspring of Affected Mothers and Fathers

Background:

Smaller studies have reported a higher offspring risk of congenital heart defects (CHDs) for mothers with CHDs than for fathers with CHDs. In a large population-based study, we investigated whether offspring risk of CHD differed for mothers and fathers with CHDs.

Methods:

All people born in Denmark, 1977 to 2011, with at least 1 registered parent, were included in our cohort (n=2 341 061). Parent-child recurrence of CHDs was evaluated using risk ratios (RRs) comparing risks of CHDs in individuals with and without a parent with a CHD, estimated using log-linear binomial regression.

Results:

The RRs for any CHD in offspring were 5.39 (95% CI, 4.88–5.96) for mothers and 3.04 (95% CI, 2.59–3.57) for fathers affected with any CHD; the ratio of RRs for mothers versus fathers was 1.82 ( P <0.0001). Recurrence RRs for the same cardiac phenotype in parent and offspring were significantly stronger for mothers than for fathers for conotruncal defects (ratio of RRs, 4.98), left ventricular outlet tract obstruction (ratio of RRs, 4.98), and ventricular septal defects (ratio of RRs, 2.51) but not for atrioventricular septal defects (ratio of RRs, 1.06). Birth rates among people with CHDs, relative to the general population, were 18% higher for women than for men, regardless of parental cardiac phenotype.

Conclusions:

Recurrence risks of CHDs were significantly greater in the offspring of affected women than in the offspring of affected men. The excess maternal recurrence risks could not be explained by the slightly higher birth rates in women with CHDs.

New developments in the biology of fibroblast growth factors

The fibroblast growth factor (FGF) family is composed of 18 secreted signaling proteins consisting of canonical FGFs and endocrine FGFs that activate four receptor tyrosine kinases (FGFRs 1-4) and four intracellular proteins (intracellular FGFs or iFGFs) that primarily function to regulate the activity of voltage-gated sodium channels and other molecules. The canonical FGFs, endocrine FGFs, and iFGFs have been reviewed extensively by us and others. In this review, we briefly summarize past reviews and then focus on new developments in the FGF field since our last review in 2015. Some of the highlights in the past 6 years include the use of optogenetic tools, viral vectors, and inducible transgenes to experimentally modulate FGF signaling, the clinical use of small molecule FGFR inhibitors, an expanded understanding of endocrine FGF signaling, functions for FGF signaling in stem cell pluripotency and differentiation, roles for FGF signaling in tissue homeostasis and regeneration, a continuing elaboration of mechanisms of FGF signaling in development, and an expanding appreciation of roles for FGF signaling in neuropsychiatric diseases. This article is categorized under: Cardiovascular Diseases > Molecular and Cellular Physiology Neurological Diseases > Molecular and Cellular Physiology Congenital Diseases > Stem Cells and Development Cancer > Stem Cells and Development.

Rbfox1 expression in amacrine cells is restricted to GABAergic and VGlut3 glycinergic cells

Rbfox1 is a multifunctional RNA binding protein that regulates alternative splicing, transcription, mRNA stability and translation. Rbfox1 is an important regulator of gene networks involved in neurogenesis and neuronal function. Disruption of Rbfox function has been associated with several neurodevelopmental and neuropsychiatric disorders. We have shown earlier that Rbfox1 is expressed in retinal ganglion and amacrine cells (ACs) and that its downregulation in adult mouse retinas leads to deficiency of depth perception. In this study, we used several markers of ACs, including GABA, choline acetyltransferase (ChAT), neuropeptide Y (NPY), glycine transporter (GlyT1) and vesicular glutamate transporter 3 (VGlut3) to identify types of ACs that express Rbfox1. Expression of Rbfox1 was observed predominantly in GABAergic ACs located in the INL and GCL. All GABAergic/cholinergic starburst ACs and virtually all NPY-positive GABAergic ACs were also Rbfox1-positive. Among glycinergic ACs, a sparse population of Rbfox1/VGlut3-positive cells was identified, indicating that Rbfox1 is expressed in a very small population of glycinergic ACs. These data contributes to our understanding about molecular differences between various types of amacrine cells and the cell-specific gene networks regulated by Rbfox1.

Recombination affects allele-specific expression of deleterious variants in human populations

How the genetic composition of a population changes through stochastic processes, such as genetic drift, in combination with deterministic processes, such as selection, is critical to understanding how phenotypes vary in space and time. Here, we show how evolutionary forces affecting selection, including recombination and effective population size, drive genomic patterns of allele-specific expression (ASE). Integrating tissue-specific genotypic and transcriptomic data from 1500 individuals from two different cohorts, we demonstrate that ASE is less often observed in regions of low recombination, and loci in high or normal recombination regions are more efficient at using ASE to underexpress harmful mutations. By tracking genetic ancestry, we discriminate between ASE variability due to past demographic effects, including subsequent bottlenecks, versus local environment. We observe that ASE is not randomly distributed along the genome and that population parameters influencing the efficacy of natural selection alter ASE levels genome wide.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

RBFOX2 is required for establishing RNA regulatory networks essential for heart development

Human genetic studies identified a strong association between loss of function mutations in RBFOX2 and hypoplastic left heart syndrome (HLHS). There are currently no Rbfox2 mouse models that recapitulate HLHS. Therefore, it is still unknown how RBFOX2 as an RNA binding protein contributes to heart development. To address this, we conditionally deleted Rbfox2 in embryonic mouse hearts and found profound defects in cardiac chamber and yolk sac vasculature formation. Importantly, our Rbfox2 conditional knockout mouse model recapitulated several molecular and phenotypic features of HLHS. To determine the molecular drivers of these cardiac defects, we performed RNA-sequencing in Rbfox2 mutant hearts and identified dysregulated alternative splicing (AS) networks that affect cell adhesion to extracellular matrix (ECM) mediated by Rho GTPases. We identified two Rho GTPase cycling genes as targets of RBFOX2. Modulating AS of these two genes using antisense oligos led to cell cycle and cell-ECM adhesion defects. Consistently, Rbfox2 mutant hearts displayed cell cycle defects and inability to undergo endocardial-mesenchymal transition, processes dependent on cell-ECM adhesion and that are seen in HLHS. Overall, our work not only revealed that loss of Rbfox2 leads to heart development defects resembling HLHS, but also identified RBFOX2-regulated AS networks that influence cell-ECM communication vital for heart development.

Genomic frontiers in congenital heart disease

The application of next-generation sequencing to study congenital heart disease (CHD) is increasingly providing new insights into the causes and mechanisms of this prevalent birth anomaly. Whole-exome sequencing analysis identifies damaging gene variants altering single or contiguous nucleotides that are assigned pathogenicity based on statistical analyses of families and cohorts with CHD, high expression in the developing heart and depletion of damaging protein-coding variants in the general population. Gene classes fulfilling these criteria are enriched in patients with CHD and extracardiac abnormalities, evidencing shared pathways in organogenesis. Developmental single-cell transcriptomic data demonstrate the expression of CHD-associated genes in particular cell lineages, and emerging insights indicate that genetic variants perturb multicellular interactions that are crucial for cardiogenesis. Whole-genome sequencing analyses extend these observations, identifying non-coding variants that influence the expression of genes associated with CHD and contribute to the estimated ~55% of unexplained cases of CHD. These approaches combined with the assessment of common and mosaic genetic variants have provided a more complete knowledge of the causes and mechanisms of CHD. Such advances provide knowledge to inform the clinical care of patients with CHD or other birth defects and deepen our understanding of the complexity of human development. In this Review, we highlight known and candidate CHD-associated human genes and discuss how the integration of advances in developmental biology research can provide new insights into the genetic contributions to CHD.

M-DATA: A statistical approach to jointly analyzing de novo mutations for multiple traits

Recent studies have demonstrated that multiple early-onset diseases have shared risk genes, based on findings from de novo mutations (DNMs). Therefore, we may leverage information from one trait to improve statistical power to identify genes for another trait. However, there are few methods that can jointly analyze DNMs from multiple traits. In this study, we develop a framework called M-DATA (Multi-trait framework for De novo mutation Association Test with Annotations) to increase the statistical power of association analysis by integrating data from multiple correlated traits and their functional annotations. Using the number of DNMs from multiple diseases, we develop a method based on an Expectation-Maximization algorithm to both infer the degree of association between two diseases as well as to estimate the gene association probability for each disease. We apply our method to a case study of jointly analyzing data from congenital heart disease (CHD) and autism. Our method was able to identify 23 genes for CHD from joint analysis, including 12 novel genes, which is substantially more than single-trait analysis, leading to novel insights into CHD disease etiology.

Genetics of congenital heart disease: a narrative review of recent advances and clinical implications

Congenital heart disease (CHD) is the most common human birth defect and remains a leading cause of mortality in childhood. Although advances in clinical management have improved the survival of children with CHD, adult survivors commonly experience cardiac and non-cardiac comorbidities, which affect quality of life and prognosis. Therefore, the elucidation of genetic etiologies of CHD not only has important clinical implications for genetic counseling of patients and families but may also impact clinical outcomes by identifying at-risk patients. Recent advancements in genetic technologies, including massively parallel sequencing, have allowed for the discovery of new genetic etiologies for CHD. Although variant prioritization and interpretation of pathogenicity remain challenges in the field of CHD genomics, advances in single-cell genomics and functional genomics using cellular and animal models of CHD have the potential to provide novel insights into the underlying mechanisms of CHD and its associated morbidities. In this review, we provide an updated summary of the established genetic contributors to CHD and discuss recent advances in our understanding of the genetic architecture of CHD along with current challenges with the interpretation of genetic variation. Furthermore, we highlight the clinical implications of genetic findings to predict and potentially improve clinical outcomes in patients with CHD.

Selection and Evaluation of a Thornless and HLB-Tolerant Bud-Sport of Pummelo Citrus With an Emphasis on Molecular Mechanisms

The selection of elite bud-sports is an important breeding approach in horticulture. We discovered and evaluated a thornless pummelo bud-sport (TL) that grew more vigorously and was more tolerant to Huanglongbing (HLB) than the thorny wild type (W). To reveal the underlying molecular mechanisms, we carried out whole-genome sequencing of W, and transcriptome comparisons of W, TL, and partially recovered thorny "mutants" (T). The results showed W, TL, and T varied in gene expression, allelic expression, and alternative splicing. Most genes/pathways with significantly altered expression in TL compared to W remained similarly altered in T. Pathway and gene ontology enrichment analysis revealed that the expression of multiple pathways, including photosynthesis and cell wall biosynthesis, was altered among the three genotypes. Remarkably, two polar auxin transporter genes, PIN7 and LAX3, were expressed at a significantly lower level in TL than in both W and T, implying alternation of polar auxin transport in TL may be responsible for the vigorous growth and thornless phenotype. Furthermore, 131 and 68 plant defense-related genes were significantly upregulated and downregulated, respectively, in TL and T compared with W. These genes may be involved in enhanced salicylic acid (SA) dependent defense and repression of defense inducing callose deposition and programmed cell death. Overall, these results indicated that the phenotype changes of the TL bud-sport were associated with tremendous transcriptome alterations, providing new clues and targets for breeding and gene editing for citrus improvement.

Nanopore sequencing reveals full-length Tropomyosin 1 isoforms and their regulation by RNA-binding proteins during rat heart development

Alternative splicing (AS) contributes to the diversity of the proteome by producing multiple isoforms from a single gene. Although short‐read RNA‐sequencing methods have been the gold standard for determining AS patterns of genes, they have a difficulty in defining full‐length mRNA isoforms assembled using different exon combinations. Tropomyosin 1 (TPM1) is an actin‐binding protein required for cytoskeletal functions in non‐muscle cells and for contraction in muscle cells. Tpm1 undergoes AS regulation to generate muscle versus non‐muscle TPM1 protein isoforms with distinct physiological functions. It is unclear which full‐length Tpm1 isoforms are produced via AS and how they are regulated during heart development. To address these, we utilized nanopore long‐read cDNA sequencing without gene‐specific PCR amplification. In rat hearts, we identified full‐length Tpm1 isoforms composed of distinct exons with specific exon linkages. We showed that Tpm1 undergoes AS transitions during embryonic heart development such that muscle‐specific exons are connected generating predominantly muscle‐specific Tpm1 isoforms in adult hearts. We found that the RNA‐binding protein RBFOX2 controls AS of rat Tpm1 exon 6a, which is important for cooperative actin binding. Furthermore, RBFOX2 regulates Tpm1 AS of exon 6a antagonistically to the RNA‐binding protein PTBP1. In sum, we defined full‐length Tpm1 isoforms with different exon combinations that are tightly regulated during cardiac development and provided insights into the regulation of Tpm1 AS by RNA‐binding proteins. Our results demonstrate that nanopore sequencing is an excellent tool to determine full‐length AS variants of muscle‐enriched genes.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

Loss of Rbfox1 Does Not Affect Survival of Retinal Ganglion Cells Injured by Optic Nerve Crush

Rbfox1 is a multifunctional RNA binding protein that regulates alternative splicing, transcription, mRNA stability and translation. Its roles in neurogenesis and neuronal functions are well established. Recent studies also implicate Rbfox1 in the regulation of gene networks that support cell survival during stress. We have earlier characterized the expression of Rbfox1 in amacrine and retinal ganglion cells (RGCs) and showed that deletion of Rbfox1 in adult animals results in depth perception deficiency. The current study investigates the effect of Rbfox1 downregulation on survival of RGCs injured by optic nerve crush (ONC). Seven days after ONC, animals sustained severe degeneration of RGC axons in the optic nerve and significant loss of RGC somas. Semi-quantitative grading of optic nerve damage in control + ONC, control + tamoxifen + ONC, and Rbfox1 ^-/- + ONC groups ranged from 4.6 to 4.8 on a scale of 1 (normal; no degenerated axons were noted) to 5 (total degeneration; all axons showed degenerated organelles, axonal content, and myelin sheath), indicating a severe degeneration. Among these three ONC groups, no statistical significance was observed when any two groups were compared. The number of RGC somas were quantitatively analyzed in superior, inferior, nasal and temporal retinal quadrants at 0.5, 1, and 1.5 mm from the center of the optic disc. The average RGC densities (cells/mm²) were: control 6,438 ± 1,203; control + ONC 2,779 ± 573; control + tamoxifen 6,163 ± 861; control + tamoxifen + ONC 2,573 ± 555; Rbfox1 ^-/- 6,437 ± 893; and Rbfox1 ^-/- + ONC 2,537 ± 526. The RGC loss in control + ONC, control + tamoxifen + ONC and Rbfox1 ^-/- + ONC was 57% (P = 1.44954E-42), 58% (P = 1.37543E-57) and 61% (P = 5.552E-59) compared to RGC numbers in the relevant uninjured groups, respectively. No statistically significant difference was observed between any two groups of uninjured animals or between any two ONC groups. Our data indicate that Rbfox1-mediated pathways have no effect on survival of RGCs injured by ONC.

Genetic Etiology of Left-Sided Obstructive Heart Lesions: A Story in Development

Congenital heart disease is the most common congenital defect observed in newborns. Within the spectrum of congenital heart disease are left‐sided obstructive lesions (LSOLs), which include hypoplastic left heart syndrome, aortic stenosis, bicuspid aortic valve, coarctation of the aorta, and interrupted aortic arch. These defects can arise in isolation or as a component of a defined syndrome; however, nonsyndromic defects are often observed in multiple family members and associated with high sibling recurrence risk. This clear evidence for a heritable basis has driven a lengthy search for disease‐causing variants that has uncovered both rare and common variants in genes that, when perturbed in cardiac development, can result in LSOLs. Despite advancements in genetic sequencing platforms and broadening use of exome sequencing, the currently accepted LSOL‐associated genes explain only 10% to 20% of patients. Further, the combinatorial effects of common and rare variants as a cause of LSOLs are emerging. In this review, we highlight the genes and variants associated with the different LSOLs and discuss the strengths and weaknesses of the present genetic associations. Furthermore, we discuss the research avenues needed to bridge the gaps in our current understanding of the genetic basis of nonsyndromic congenital heart disease.

Weighted gene co-expression network analysis identified underlying hub genes and mechanisms in the occurrence and development of viral myocarditis

Background

Myocarditis is an inflammatory myocardial disease, which may lead to heart failure and sudden death. Despite extensive research into the pathogenesis of myocarditis, effective treatments for this condition remain elusive. This study aimed to explore the potential pathogenesis and hub genes for viral myocarditis.

Methods

A weighted gene co-expression network analysis (WGCNA) was performed based on the gene expression profiles derived from mouse models at different stages of viral myocarditis (GSE35182). Functional annotation was executed within the key modules. Potential hub genes were predicted based on the intramodular connectivity (IC). Finally, potential microRNAs that regulate gene expression were predicted by miRNet analysis.

Results

Three gene co-expression modules showed the strongest correlation with the acute or chronic disease stage. A significant positive correlation was detected between the acute disease stage and the turquoise module, the genes of which were mainly enriched in antiviral response and immune-inflammatory activation. Furthermore, a significant positive correlation and a negative correlation were identified between the chronic disease stage and the brown and yellow modules, respectively. These modules were mainly associated with the cytoskeleton, phosphorylation, cellular catabolic process, and autophagy. Subsequently, we predicted the underlying hub genes and microRNAs in the three modules.

Conclusions

This study revealed the main biological processes in different stages of viral myocarditis and predicted hub genes in both the acute and chronic disease stages. Our results may be helpful for developing new therapeutic targets for viral myocarditis in future research.

ORE identifies extreme expression effects enriched for rare variants

MOTIVATION

Non-coding rare variants (RVs) may contribute to Mendelian disorders but have been challenging to study due to small sample sizes, genetic heterogeneity and uncertainty about relevant non-coding features. Previous studies identified RVs associated with expression outliers, but varying outlier definitions were employed and no comprehensive open-source software was developed.

RESULTS

We developed Outlier-RV Enrichment (ORE) to identify biologically-meaningful non-coding RVs. We implemented ORE combining whole-genome sequencing and cardiac RNAseq from congenital heart defect patients from the Pediatric Cardiac Genomics Consortium and deceased adults from Genotype-Tissue Expression. Use of rank-based outliers maximized sensitivity while a most extreme outlier approach maximized specificity. Rarer variants had stronger associations, suggesting they are under negative selective pressure and providing a basis for investigating their contribution to Mendelian disorders.

AVAILABILITY AND IMPLEMENTATION

ORE, source code, and documentation are available at https://pypi.python.org/pypi/ore under the MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Aberrant Expression of a Non-muscle RBFOX2 Isoform Triggers Cardiac Conduction Defects in Myotonic Dystrophy

Myotonic dystrophy type 1 (DM1) is a multisystemic genetic disorder caused by the CTG repeat expansion in the 3'-untranslated region of DMPK gene. Heart dysfunctions occur in ∼80% of DM1 patients and are the second leading cause of DM1-related deaths. Herein, we report that upregulation of a non-muscle splice isoform of RNA-binding protein RBFOX2 in DM1 heart tissue-due to altered splicing factor and microRNA activities-induces cardiac conduction defects in DM1 individuals. Mice engineered to express the non-muscle RBFOX240 isoform in heart via tetracycline-inducible transgenesis, or CRISPR/Cas9-mediated genome editing, reproduced DM1-related cardiac conduction delay and spontaneous episodes of arrhythmia. Further, by integrating RNA binding with cardiac transcriptome datasets from DM1 patients and mice expressing the non-muscle RBFOX2 isoform, we identified RBFOX240-driven splicing defects in voltage-gated sodium and potassium channels, which alter their electrophysiological properties. Thus, our results uncover a trans-dominant role for an aberrantly expressed RBFOX240 isoform in DM1 cardiac pathogenesis.

Quantitative Phenotyping of Xenopus Embryonic Heart Pathophysiology Using Hemoglobin Contrast Subtraction Angiography to Screen Human Cardiomyopathies

Congenital heart disease (CHD) is a significant cause of mortality in infants and adults. Currently human genomic analysis has identified a number of candidate genes in these patients. These genes span diverse categories of gene function suggesting that despite the similarity in cardiac lesion, the underlying pathophysiology may be different. In fact, patients with similar CHDs can have drastically different outcomes, including a dramatic decrease in myocardial function. To test these human candidate genes for their impact on myocardial function, we need efficient animals models of disease. For this purpose, we paired Xenopus tropicalis with our microangiography technique, hemoglobin contrast subtraction angiography (HCSA). To demonstrate the gene-teratogen-physiology relationship, we modeled human cardiomyopathy in tadpoles. First we depleted the sarcomeric protein myosin heavy chain 6 (myh6) expression using morpholino oligos. Next, we exposed developing embryos to the teratogen ethanol and in both conditions showed varying degrees of cardiac dysfunction. Our results demonstrate that HCSA can distinguish biomechanical phenotypes in the context of gene dysfunction or teratogen. This approach can be used to screen numerous candidate CHD genes or suspected teratogens for their effect on cardiac function.

Genetic regulatory variation in populations informs transcriptome analysis in rare disease

Transcriptome data can facilitate the interpretation of the effects of rare genetic variants. Here, we introduce ANEVA (analysis of expression variation) to quantify genetic variation in gene dosage from allelic expression (AE) data in a population. Application of ANEVA to the Genotype-Tissues Expression (GTEx) data showed that this variance estimate is robust and correlated with selective constraint in a gene. Using these variance estimates in a dosage outlier test (ANEVA-DOT) applied to AE data from 70 Mendelian muscular disease patients showed accuracy in detecting genes with pathogenic variants in previously resolved cases and led to one confirmed and several potential new diagnoses. Using our reference estimates from GTEx data, ANEVA-DOT can be incorporated in rare disease diagnostic pipelines to use RNA-sequencing data more effectively.

Small-protein Enrichment Assay Enables the Rapid, Unbiased Analysis of Over 100 Low Abundance Factors from Human Plasma

Unbiased and sensitive quantification of low abundance small proteins in human plasma (e.g. hormones, immune factors, metabolic regulators) remains an unmet need. These small protein factors are typically analyzed individually and using antibodies that can lack specificity. Mass spectrometry (MS)-based proteomics has the potential to address these problems, however the analysis of plasma by MS is plagued by the extremely large dynamic range of this body fluid, with protein abundances spanning at least 13 orders of magnitude. Here we describe an enrichment assay (SPEA), that greatly simplifies the plasma dynamic range problem by enriching small-proteins of 2-10 kDa, enabling the rapid, specific and sensitive quantification of >100 small-protein factors in a single untargeted LC-MS/MS acquisition. Applying this method to perform deep-proteome profiling of human plasma we identify C5ORF46 as a previously uncharacterized human plasma protein. We further demonstrate the reproducibility of our workflow for low abundance protein analysis using a stable-isotope labeled protein standard of insulin spiked into human plasma. SPEA provides the ability to study numerous important hormones in a single rapid assay, which we applied to study the intermittent fasting response and observed several unexpected changes including decreased plasma abundance of the iron homeostasis regulator hepcidin.

Single-Cell RNA-Seq Technologies and Related Computational Data Analysis

Single-cell RNA sequencing (scRNA-seq) technologies allow the dissection of gene expression at single-cell resolution, which greatly revolutionizes transcriptomic studies. A number of scRNA-seq protocols have been developed, and these methods possess their unique features with distinct advantages and disadvantages. Due to technical limitations and biological factors, scRNA-seq data are noisier and more complex than bulk RNA-seq data. The high variability of scRNA-seq data raises computational challenges in data analysis. Although an increasing number of bioinformatics methods are proposed for analyzing and interpreting scRNA-seq data, novel algorithms are required to ensure the accuracy and reproducibility of results. In this review, we provide an overview of currently available single-cell isolation protocols and scRNA-seq technologies, and discuss the methods for diverse scRNA-seq data analyses including quality control, read mapping, gene expression quantification, batch effect correction, normalization, imputation, dimensionality reduction, feature selection, cell clustering, trajectory inference, differential expression calling, alternative splicing, allelic expression, and gene regulatory network reconstruction. Further, we outline the prospective development and applications of scRNA-seq technologies.

TAMM41 is required for heart valve differentiation via regulation of PINK-PARK2 dependent mitophagy

TAMM41, located within the congenital heart diseases (CHD) sensitive region of 3p25 deletion syndrome, is a mitochondrial membrane maintenance protein critical for yeast survival, but its function in higher vertebrates remains unknown. Via in vivo zebrafish model, we found that tamm41 is highly expressed in the developing heart and deficiency of which led to heart valve abnormalities. Molecular mechanistic studies revealed that TAMM41 interacts and modulates the PINK1-PARK2 dependent mitophagy pathway, thereby implicating TAMM41 in heart valve development during zebrafish embryonic cardiogenesis. Furthermore, through screening of the congenital heart diseases (CHD) sensitive region of 3p25 deletion syndrome among 118 sporadic atrioventricular septal defect (AVSD) patients, we identified three cases carrying heterozygous pathogenic intronic variants of TAMM41. All three cases lacked normal full-length TAMM41 transcripts, most likely due to specific expression of the mutant allele. Collectively, our studies highlight essential roles for TAMM41-dependent mitophagy in development of the heart and provide novel insights into the etiology of AVSD.

A Rapid Method for Directed Gene Knockout for Screening in G0 Zebrafish

Zebrafish is a powerful model for forward genetics. Reverse genetic approaches are limited by the time required to generate stable mutant lines. We describe a system for gene knockout that consistently produces null phenotypes in G0 zebrafish. Yolk injection of sets of four CRISPR/Cas9 ribonucleoprotein complexes redundantly targeting a single gene recapitulated germline-transmitted knockout phenotypes in >90% of G0 embryos for each of 8 test genes. Early embryonic (6 hpf) and stable adult phenotypes were produced. Simultaneous multi-gene knockout was feasible but associated with toxicity in some cases. To facilitate use, we generated a lookup table of four-guide sets for 21,386 zebrafish genes and validated several. Using this resource, we targeted 50 cardiomyocyte transcriptional regulators and uncovered a role of zbtb16a in cardiac development. This system provides a platform for rapid screening of genes of interest in development, physiology, and disease models in zebrafish.

ZBTB16 and metabolic syndrome: a network perspective

Metabolic syndrome is a prevalent, complex condition. The search for genetic determinants of the syndrome is currently undergoing a paradigm enhancement by adding systems genetics approaches to association studies. We summarize the current evidence on relations between an emergent new candidate, zinc finger and BTB domain containing 16 (ZBTB16) transcription factor and the major components constituting the metabolic syndrome. Information stemming from studies on experimental models with altered Zbtb16 expression clearly shows its effect on adipogenesis, cardiac hypertrophy and fibrosis, lipid levels and insulin sensitivity. Based on current evidence, we provide a network view of relations between ZBTB16 and hallmarks of metabolic syndrome in order to elucidate the potential functional links involving the ZBTB16 node. Many of the identified genes interconnecting ZBTB16 with all or most metabolic syndrome components are linked to immune function, inflammation or oxidative stress. In summary, ZBTB16 represents a promising pleiotropic candidate node for metabolic syndrome.

Hippo Signaling Plays an Essential Role in Cell State Transitions during Cardiac Fibroblast Development

During development, progenitors progress through transition states. The cardiac epicardium contains progenitors of essential non-cardiomyocytes. The Hippo pathway, a kinase cascade that inhibits the Yap transcriptional co-factor, controls organ size in developing hearts. Here, we investigated Hippo kinases Lats1 and Lats2 in epicardial diversification. Epicardial-specific deletion of Lats1/2 was embryonic lethal, and mutant embryos had defective coronary vasculature remodeling. Single-cell RNA sequencing revealed that Lats1/2 mutant cells failed to activate fibroblast differentiation but remained in an intermediate cell state with both epicardial and fibroblast characteristics. Lats1/2 mutant cells displayed an arrested developmental trajectory with persistence of epicardial markers and expanded expression of Yap targets Dhrs3, an inhibitor of retinoic acid synthesis, and Dpp4, a protease that modulates extracellular matrix (ECM) composition. Genetic and pharmacologic manipulation revealed that Yap inhibits fibroblast differentiation, prolonging a subepicardial-like cell state, and promotes expression of matricellular factors, such as Dpp4, that define ECM characteristics.

rbFOX1/MBNL1 competition for CCUG RNA repeats binding contributes to myotonic dystrophy type 1/type 2 differences

Myotonic dystrophy type 1 and type 2 (DM1, DM2) are caused by expansions of CTG and CCTG repeats, respectively. RNAs containing expanded CUG or CCUG repeats interfere with the metabolism of other RNAs through titration of the Muscleblind-like (MBNL) RNA binding proteins. DM2 follows a more favorable clinical course than DM1, suggesting that specific modifiers may modulate DM severity. Here, we report that the rbFOX1 RNA binding protein binds to expanded CCUG RNA repeats, but not to expanded CUG RNA repeats. Interestingly, rbFOX1 competes with MBNL1 for binding to CCUG expanded repeats and overexpression of rbFOX1 partly releases MBNL1 from sequestration within CCUG RNA foci in DM2 muscle cells. Furthermore, expression of rbFOX1 corrects alternative splicing alterations and rescues muscle atrophy, climbing and flying defects caused by expression of expanded CCUG repeats in a Drosophila model of DM2.

Physiology of Cardiac Development: From Genetics to Signaling to Therapeutic Strategies

The heart is one of the first organs to form and function during embryonic development. It is comprised of multiple cell lineages, each integral for proper cardiac development, and include cardiomyocytes, endothelial cells, epicardial cells and neural crest cells. The molecular mechanisms regulating cardiac development and morphogenesis are dependent on signaling crosstalk between multiple lineages through paracrine interactions, cell-ECM interactions, and cell-cell interactions, which together, help facilitate survival, growth, proliferation, differentiation and migration of cardiac tissue. Aberrant regulation of any of these processes can induce developmental disorders and pathological phenotypes. Here, we will discuss each of these processes, the genetic factors that contribute to each step of cardiac development, as well as the current and future therapeutic targets and mechanisms of heart development and disease. Understanding the complex interactions that regulate cardiac development, proliferation and differentiation is not only vital to understanding the causes of congenital heart defects, but to also finding new therapeutics that can treat both pediatric and adult cardiac disease in the near future.

A Path to Implement Precision Child Health Cardiovascular Medicine

Congenital heart defects (CHDs) affect approximately 1% of live births and are a major source of childhood morbidity and mortality even in countries with advanced healthcare systems. Along with phenotypic heterogeneity, the underlying etiology of CHDs is multifactorial, involving genetic, epigenetic, and/or environmental contributors. Clear dissection of the underlying mechanism is a powerful step to establish individualized therapies. However, the majority of CHDs are yet to be clearly diagnosed for the underlying genetic and environmental factors, and even less with effective therapies. Although the survival rate for CHDs is steadily improving, there is still a significant unmet need for refining diagnostic precision and establishing targeted therapies to optimize life quality and to minimize future complications. In particular, proper identification of disease associated genetic variants in humans has been challenging, and this greatly impedes our ability to delineate gene–environment interactions that contribute to the pathogenesis of CHDs. Implementing a systematic multileveled approach can establish a continuum from phenotypic characterization in the clinic to molecular dissection using combined next-generation sequencing platforms and validation studies in suitable models at the bench. Key elements necessary to advance the field are: first, proper delineation of the phenotypic spectrum of CHDs; second, defining the molecular genotype/phenotype by combining whole-exome sequencing and transcriptome analysis; third, integration of phenotypic, genotypic, and molecular datasets to identify molecular network contributing to CHDs; fourth, generation of relevant disease models and multileveled experimental investigations. In order to achieve all these goals, access to high-quality biological specimens from well-defined patient cohorts is a crucial step. Therefore, establishing a CHD BioCore is an essential infrastructure and a critical step on the path toward precision child health cardiovascular medicine.

Genetics and Genomics of Congenital Heart Disease

Congenital heart disease is the most common birth defect, and because of major advances in medical and surgical management, there are now more adults living with congenital heart disease (CHD) than children. Until recently, the cause of the majority of CHD was unknown. Advances in genomic technologies have discovered the genetic causes of a significant fraction of CHD, while at the same time pointing to remarkable complexity in CHD genetics. This review will focus on the evidence for genetic causes underlying CHD and discuss data supporting both monogenic and complex genetic mechanisms underlying CHD. The discoveries from CHD genetic studies draw attention to biological pathways that simultaneously open the door to a better understanding of cardiac development and affect clinical care of patients with CHD. Finally, we address clinical genetic evaluation of patients and families affected by CHD.

Single-Cell Resolution of Temporal Gene Expression during Heart Development

Activation of complex molecular programs in specific cell lineages governs mammalian heart development, from a primordial linear tube to a four-chamber organ. To characterize lineage-specific, spatiotemporal developmental programs, we performed single-cell RNA sequencing of >1,200 murine cells isolated at seven time points spanning embryonic day 9.5 (primordial heart tube) to postnatal day 21 (mature heart). Using unbiased transcriptional data, we classified cardiomyocytes, endothelial cells, and fibroblast-enriched cells, thus identifying markers for temporal and chamber-specific developmental programs. By harnessing these datasets, we defined developmental ages of human and mouse pluripotent stem-cell-derived cardiomyocytes and characterized lineage-specific maturation defects in hearts of mice with heterozygous mutations in Nkx2.5 that cause human heart malformations. This spatiotemporal transcriptome analysis of heart development reveals lineage-specific gene programs underlying normal cardiac development and congenital heart disease.