Right Ventricular Sarcomere Contractile Depression and the Role of Thick Filament Activation in Human Heart Failure With Pulmonary Hypertension.

BACKGROUND: Right ventricular (RV) contractile dysfunction commonly occurs and worsens outcomes in patients with heart failure with reduced ejection fraction and pulmonary hypertension (HFrEF-PH). However, such dysfunction often goes undetected by standard clinical RV indices, raising concerns that they may not reflect aspects of underlying myocyte dysfunction. We thus sought to characterize RV myocyte contractile depression in HFrEF-PH, identify those components reflected by clinical RV indices, and uncover underlying biophysical mechanisms. METHODS: Resting, calcium-, and load-dependent mechanics were prospectively studied in permeabilized RV cardiomyocytes isolated from explanted hearts from 23 patients with HFrEF-PH undergoing cardiac transplantation and 9 organ donor controls. RESULTS: Unsupervised machine learning using myocyte mechanical data with the highest variance yielded 2 HFrEF-PH subgroups that in turn mapped to patients with decompensated or compensated clinical RV function. This correspondence was driven by reduced calcium-activated isometric tension in decompensated clinical RV function, whereas surprisingly, many other major myocyte contractile measures including peak power and myocyte active stiffness were similarly depressed in both groups. Similar results were obtained when subgroups were first defined by clinical indices, and then myocyte mechanical properties in each group compared. To test the role of thick filament defects, myofibrillar structure was assessed by x-ray diffraction of muscle fibers. This revealed more myosin heads associated with the thick filament backbone in decompensated clinical RV function, but not compensated clinical RV function, as compared with controls. This corresponded to reduced myosin ATP turnover in decompensated clinical RV function myocytes, indicating less myosin in a crossbridge-ready disordered-relaxed (DRX) state. Altering DRX proportion (%DRX) affected peak calcium-activated tension in the patient groups differently, depending on their basal %DRX, highlighting potential roles for precision-guided therapeutics. Last, increasing myocyte preload (sarcomere length) increased %DRX 1.5-fold in controls but only 1.2-fold in both HFrEF-PH groups, revealing a novel mechanism for reduced myocyte active stiffness and by extension Frank-Starling reserve in human heart failure. CONCLUSIONS: Although there are many RV myocyte contractile deficits in HFrEF-PH, commonly used clinical indices only detect reduced isometric calcium-stimulated force, which is related to deficits in basal and recruitable %DRX myosin. Our results support use of therapies to increase %DRX and enhance length-dependent recruitment of DRX myosin heads in such patients.

Patterns of Unwanted Biological and Technical Expression Variation Among 49 Human Tissues.

Tissue gene expression studies are impacted by biological and technical sources of variation, which can be broadly classified into wanted and unwanted variation. The latter, if not addressed, results in misleading biological conclusions. Methods have been proposed to reduce unwanted variation, such as normalization and batch correction. A more accurate understanding of all causes of variation could significantly improve the ability of these methods to remove unwanted variation while retaining variation corresponding to the biological question of interest. We used 17,282 samples from 49 human tissues in the Genotype-Tissue Expression data set (v8) to investigate patterns and causes of expression variation. Transcript expression was transformed to z-scores, and only the most variable 2% of transcripts were evaluated and clustered based on coexpression patterns. Clustered gene sets were assigned to different biological or technical causes based on histologic appearances and metadata elements. We identified 522 variable transcript clusters (median: 11 per tissue) among the samples. Of these, 63% were confidently explained, 16% were likely explained, 7% were low confidence explanations, and 14% had no clear cause. Histologic analysis annotated 46 clusters. Other common causes of variability included sex, sequencing contamination, immunoglobulin diversity, and compositional tissue differences. Less common biological causes included death interval (Hardy score), disease status, and age. Technical causes included blood draw timing and harvesting differences. Many of the causes of variation in bulk tissue expression were identifiable in the Tabula Sapiens data set of single-cell expression. This is among the largest explorations of the underlying sources of tissue expression variation. It uncovered expected and unexpected causes of variable gene expression and demonstrated the utility of matched histologic specimens. It further demonstrated the value of acquiring meaningful tissue harvesting metadata elements to use for improved normalization, batch correction, and analysis of both bulk and single-cell RNA-seq data.

Bile Acid Receptor Agonist Reverses Transforming Growth Factor-ß1-Mediated Fibrogenesis in Human Induced Pluripotent Stem Cells-Derived Kidney Organoids.

Chronic kidney disease progresses through the replacement of functional tissue compartments with fibrosis, a maladaptive repair process. Shifting kidney repair toward a physiologically intact architecture, rather than fibrosis, is key to blocking chronic kidney disease progression. Much research into the mechanisms of fibrosis is performed in rodent models with less attention to the human genetic context. Recently, human induced pluripotent stem cell (iPSC)-derived organoids have shown promise in overcoming the limitation. In this study, we developed a fibrosis model that uses human iPSC-based 3-dimensional renal organoids, in which exogenous transforming growth factor-ß1 (TGF-ß1) induced the production of extracellular matrix. TGF-ß1-treated organoids showed tubulocentric collagen 1α1 production by regulating downstream transcriptional regulators, Farnesoid X receptor, phosphorylated mothers against decapentaplegic homolog 3 (p-SMAD3), and transcriptional coactivator with PDZ-binding motif (TAZ). Increased nuclear TAZ expression was confirmed in the tubular epithelium in human kidney biopsies with tubular injury and early fibrosis. A dual bile acid receptor agonist (INT-767) increased Farnesoid X receptor and reduced p-SMAD3 and TAZ, attenuating TGF-ß1-induced fibrosis in kidney organoids. Finally, we show that TAZ interacted with TEA-domain transcription factors and p-SMAD3 with TAZ and TEA-domain transcription factor 4 coregulating collagen 1α1 gene transcription. In summary, we establish a novel, readily manipulable fibrogenesis model and posit a role for bile acid receptor agonism early in renal parenchymal fibrosis.

Expression Microdissection for the Analysis of miRNA in a Single-Cell Type.

Cell-specific microRNA (miRNA) expression estimates are important in characterizing the localization of miRNA signaling within tissues. Much of these data are obtained from cultured cells, a process known to significantly alter miRNA expression levels. Thus, our knowledge of in vivo cell miRNA expression estimates is poor. We previously demonstrated expression microdissection-miRNA-sequencing (xMD-miRNA-seq) to acquire in vivo estimates, directly from formalin-fixed tissues, albeit with a limited yield. In this study, we optimized each step of the xMD process, including tissue retrieval, tissue transfer, film preparation, and RNA isolation, to increase RNA yields and ultimately show strong enrichment for in vivo miRNA expression by qPCR array. These method improvements, such as the development of a noncrosslinked ethylene vinyl acetate membrane, resulted in a 23- to 45-fold increase in miRNA yield, depending on the cell type. By qPCR, miR-200a increased by 14-fold in xMD-derived small intestine epithelial cells, with a concurrent 336-fold reduction in miR-143 relative to the matched nondissected duodenal tissue. xMD is now an optimized method to obtain robust in vivo miRNA expression estimates from cells. xMD will allow formalin-fixed tissues from surgical pathology archives to make theragnostic biomarker discoveries.

Direct observation of the evolution of cell-type-specific microRNA expression signatures supports the hematopoietic origin model of endothelial cells.

The evolution of specialized cell-types is a long-standing interest of biologists, but given the deep time-scales very difficult to reconstruct or observe. microRNAs have been linked to the evolution of cellular complexity and may inform on specialization. The endothelium is a vertebrate-specific specialization of the circulatory system that enabled a critical new level of vasoregulation. The evolutionary origin of these endothelial cells is unclear. We hypothesized that Mir-126, an endothelial cell-specific microRNA may be informative. We here reconstruct the evolutionary history of Mir-126. Mir-126 likely appeared in the last common ancestor of vertebrates and tunicates, which was a species without an endothelium, within an intron of the evolutionary much older EGF Like Domain Multiple (Egfl) locus. Mir-126 has a complex evolutionary history due to duplications and losses of both the host gene and the microRNA. Taking advantage of the strong evolutionary conservation of the microRNA among Olfactores, and using RNA in situ hybridization, we localized Mir-126 in the tunicate Ciona robusta. We found exclusive expression of the mature Mir-126 in granular amebocytes, supporting a long-proposed scenario that endothelial cells arose from hemoblasts, a type of proto-endothelial amoebocyte found throughout invertebrates. This observed change of expression of Mir-126 from proto-endothelial amoebocytes in the tunicate to endothelial cells in vertebrates is the first direct observation of the evolution of a cell-type in relation to microRNA expression indicating that microRNAs can be a prerequisite of cell-type evolution.

The effect of tissue composition on gene co-expression.

Variable cellular composition of tissue samples represents a significant challenge for the interpretation of genomic profiling studies. Substantial effort has been devoted to modeling and adjusting for compositional differences when estimating differential expression between sample types. However, relatively little attention has been given to the effect of tissue composition on co-expression estimates. In this study, we illustrate the effect of variable cell-type composition on correlation-based network estimation and provide a mathematical decomposition of the tissue-level correlation. We show that a class of deconvolution methods developed to separate tumor and stromal signatures can be applied to two component cell-type mixtures. In simulated and real data, we identify conditions in which a deconvolution approach would be beneficial. Our results suggest that uncorrelated cell-type-specific markers are ideally suited to deconvolute both the expression and co-expression patterns of an individual cell type. We provide a Shiny application for users to interactively explore the effect of cell-type composition on correlation-based co-expression estimation for any cell types of interest.

A Comparison of Tissue Dissection Techniques for Diagnostic, Prognostic, and Theragnostic Analysis of Human Disease.

Histopathology has historically been the critical technique for the diagnosis and treatment of human disease. Today, genomics, transcriptomics, and proteomics from specific cells, rather than bulk tissue, have become key to understanding underlying disease mechanisms and rendering useful diagnostic information. Extraction of desired analytes, i.e., nucleic acids or proteins, from easily accessible formalin-fixed paraffin-embedded tissues allows for clinically relevant activities, such as sequencing biomarker mutations or typing amyloidogenic proteins. Genetic profiling has become routine for cancers as varied as non-small cell lung cancer and prostatic carcinoma. The five main tissue dissection techniques that have been developed thus far include: bulk scraping, manual macrodissection, manual microdissection, laser-capture microdissection, and expression microdissection. In this review, we discuss the importance of tissue dissection in clinical practice and research, the basic methods, applications, as well as some advantages and disadvantages for each modality.

Mycobacterium tuberculosis Infection Drives Mitochondria-Biased Dysregulation of Host Transfer RNA-Derived Fragments.

BACKGROUND: Mycobacterium tuberculosis (Mtb), the bacterium that causes tuberculosis, causes 10 million infections and 1.5 million deaths per year worldwide. The success of Mtb as a human pathogen is directly related to its ability to suppress host responses, which are critical for clearing intracellular pathogens. Emerging evidence suggests that key response pathways may be regulated by a novel class of small noncoding RNA, called transfer RNA (tRNA)-derived fragments (tRFs). tRFs can complex with Argonaute proteins to target and degrade messenger RNA targets, similarly to micro RNAs, but have thus far been overlooked in the context of bacterial infections. METHODS: We generated a novel miRge2.0-based tRF-analysis tool, tRFcluster, and used it to analyze independently generated and publicly available RNA-sequencing datasets to assess tRF dysregulation in host cells following infection with Mtb and other intracellular bacterial pathogens. RESULTS: We found that Mtb and Listeria monocytogenes drive dramatic tRF dysregulation, whereas other bacterial pathogens do not. Interestingly, Mtb infection uniquely increased the expression of mitochondria-derived tRFs rather than genomic-derived tRFs, suggesting an association with mitochondrial damage in Mtb infection. CONCLUSIONS: tRFs are dysregulated in some, but not all, bacterial infections. Biased dysregulation of mitochondria-derived tRFs in Mtb infection suggests a link between mitochondrial distress and tRF production.

Spatial Proteomic Approach to Characterize Skeletal Muscle Myofibers.

Skeletal muscle myofibers have differential protein expression resulting in functionally distinct slow- and fast-twitch types. While certain protein classes are well-characterized, the depth of all proteins involved in this process is unknown. We utilized the Human Protein Atlas (HPA) and the HPASubC tool to classify mosaic expression patterns of staining across 49,600 unique tissue microarray (TMA) images using a visual proteomic approach. We identified 2164 proteins with potential mosaic expression, of which 1605 were categorized as "likely" or "real." This list included both well-known fiber-type-specific and novel proteins. A comparison of the 1605 mosaic proteins with a mass spectrometry (MS)-derived proteomic dataset of single human muscle fibers led to the assignment of 111 proteins to fiber types. We additionally used a multiplexed immunohistochemistry approach, a multiplexed RNA-ISH approach, and STRING v11 to further assign or suggest fiber types of newly characterized mosaic proteins. This visual proteomic analysis of mature skeletal muscle myofibers greatly expands the known repertoire of twitch-type-specific proteins.

Big Strides in Cellular MicroRNA Expression.

A lack of knowledge of the cellular origin of miRNAs has greatly confounded functional and biomarkers studies. Recently, three studies characterized miRNA expression patterns across >78 human cell types. These combined data expand our knowledge of miRNA expression localization and confirm that many miRNAs show cell type-specific expression patterns.

Human cardiac cis-regulatory elements, their cognate transcription factors, and regulatory DNA sequence variants.

Cis-regulatory elements (CRE), short DNA sequences through which transcription factors (TFs) exert regulatory control on gene expression, are postulated to be the major sites of causal sequence variation underlying the genetics of complex traits and diseases. We present integrative analyses, combining high-throughput genomic and epigenomic data with sequence-based computations, to identify the causal transcriptional components in a given tissue. We use data on adult human hearts to demonstrate that (1) sequence-based predictions detect numerous, active, tissue-specific CREs missed by experimental observations, (2) learned sequence features identify the cognate TFs, (3) CRE variants are specifically associated with cardiac gene expression, and (4) a significant fraction of the heritability of exemplar cardiac traits (QT interval, blood pressure, pulse rate) is attributable to these variants. This general systems approach can thus identify candidate causal variants and the components of gene regulatory networks (GRN) to enable understanding of the mechanisms of complex disorders on a tissue- or cell-type basis.

Discovery of noncanonical translation initiation sites through mass spectrometric analysis of protein N termini.

Translation initiation generally occurs at AUG codons in eukaryotes, although it has been shown that non-AUG or noncanonical translation initiation can also occur. However, the evidence for noncanonical translation initiation sites (TISs) is largely indirect and based on ribosome profiling (Ribo-seq) studies. Here, using a strategy specifically designed to enrich N termini of proteins, we demonstrate that many human proteins are translated at noncanonical TISs. The large majority of TISs that mapped to 5' untranslated regions were noncanonical and led to N-terminal extension of annotated proteins or translation of upstream small open reading frames (uORF). It has been controversial whether the amino acid corresponding to the start codon is incorporated at the TIS or methionine is still incorporated. We found that methionine was incorporated at almost all noncanonical TISs identified in this study. Comparison of the TISs determined through mass spectrometry with ribosome profiling data revealed that about two-thirds of the novel annotations were indeed supported by the available ribosome profiling data. Sequence conservation across species and a higher abundance of noncanonical TISs than canonical ones in some cases suggests that the noncanonical TISs can have biological functions. Overall, this study provides evidence of protein translation initiation at noncanonical TISs and argues that further studies are required for elucidation of functional implications of such noncanonical translation initiation.

Unification of miRNA and isomiR research: the mirGFF3 format and the mirtop API.

MOTIVATION: MicroRNAs (miRNAs) are small RNA molecules (∼22 nucleotide long) involved in post-transcriptional gene regulation. Advances in high-throughput sequencing technologies led to the discovery of isomiRs, which are miRNA sequence variants. While many miRNA-seq analysis tools exist, the diversity of output formats hinders accurate comparisons between tools and precludes data sharing and the development of common downstream analysis methods. RESULTS: To overcome this situation, we present here a community-based project, miRNA Transcriptomic Open Project (miRTOP) working towards the optimization of miRNA analyses. The aim of miRTOP is to promote the development of downstream isomiR analysis tools that are compatible with existing detection and quantification tools. Based on the existing GFF3 format, we first created a new standard format, mirGFF3, for the output of miRNA/isomiR detection and quantification results from small RNA-seq data. Additionally, we developed a command line Python tool, mirtop, to create and manage the mirGFF3 format. Currently, mirtop can convert into mirGFF3 the outputs of commonly used pipelines, such as seqbuster, isomiR-SEA, sRNAbench, Prost! as well as BAM files. Some tools have also incorporated the mirGFF3 format directly into their code, such as, miRge2.0, IsoMIRmap and OptimiR. Its open architecture enables any tool or pipeline to output or convert results into mirGFF3. Collectively, this isomiR categorization system, along with the accompanying mirGFF3 and mirtop API, provide a comprehensive solution for the standardization of miRNA and isomiR annotation, enabling data sharing, reporting, comparative analyses and benchmarking, while promoting the development of common miRNA methods focusing on downstream steps of miRNA detection, annotation and quantification. AVAILABILITY AND IMPLEMENTATION: https://github.com/miRTop/mirGFF3/ and https://github.com/miRTop/mirtop. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Left ventricular fibro-fatty replacement in arrhythmogenic right ventricular dysplasia/cardiomyopathy: prevalence, patterns, and association with arrhythmias.

BACKGROUND: Left ventricular (LV) fibrofatty infiltration in arrhythmogenic right ventricular (RV) dysplasia/cardiomyopathy (ARVD/C) has been reported, however, detailed cardiovascular magnetic resonance (CMR) characteristics and association with outcomes are uncertain. We aim to describe LV findings on CMR in ARVD/C patients and their relationship with arrhythmic outcomes. METHODS: CMR of 73 subjects with ARVD/C according to the 2010 Task Force Criteria (TFC) were analyzed for LV involvement, defined as ≥ 1 of the following features: LV wall motion abnormality, LV late gadolinium enhancement (LGE), LV fat infiltration, or LV ejection fraction (LVEF) < 50%. Ventricular volumes and function, regional wall motion abnormalities, and the presence of ventricular fat or fibrosis were recorded. Findings on CMR were correlated with arrhythmic outcomes. RESULTS: Of the 73 subjects, 50.7% had CMR evidence for LV involvement. Proband status and advanced RV dysfunction were independently associated with LV abnormalities. The most common pattern of LV involvement was focal fatty infiltration in the sub-epicardium of the apicolateral LV with a "bite-like" pattern. LGE in the LV was found in the same distribution and most often had a linear appearance. LV involvement was more common with non-PKP2 genetic mutation variants, regardless of proband status. Only RV structural disease on CMR (HR 3.47, 95% CI 1.13-10.70) and prior arrhythmia (HR 2.85, 95% CI 1.33-6.10) were independently associated with arrhythmic events. CONCLUSION: Among patients with 2010 TFC for ARVD/C, CMR evidence for LV abnormalities are seen in half of patients and typically manifest as fibrofatty infiltration in the subepicardium of the apicolateral wall and are not associated with arrhythmic outcomes.

Siglec-1 Macrophages and the Contribution of IFN to the Development of Autoimmune Congenital Heart Block.

Given that diseases associated with anti-SSA/Ro autoantibodies, such as systemic lupus erythematosus and Sjögren syndrome, are linked with an upregulation of IFN and type I IFN-stimulated genes, including sialic acid-binding Ig-like lectin 1 (Siglec-1), a receptor on monocytes/macrophages, recent attention has focused on a potential role for IFN and IFN-stimulated genes in the pathogenesis of congenital heart block (CHB). Accordingly, three approaches were leveraged to address the association of IFN, IFN-stimulated genes, and the phenotype of macrophages in affected fetal cardiac tissue: 1) cultured healthy human macrophages transfected with hY3, an anti-SSA/Ro-associated ssRNA, 2) RNA isolated from freshly sorted human leukocytes/macrophages after Langendorff perfusion of three fetal hearts dying with CHB and three healthy gestational age-matched hearts, and 3) autopsy tissue from three additional human CHB hearts and one healthy heart. TLR ligation of macrophages with hY3 led to the upregulation of a panel of IFN transcripts, including SIGLEC1, a result corroborated using quantitative PCR. Using independent and agnostic bioinformatics approaches, CD45+CD11c+ and CD45+CD11c- human leukocytes flow sorted from the CHB hearts highly expressed type I IFN response genes inclusive of SIGLEC1. Furthermore, Siglec-1 expression was identified in the septal region of several affected fetal hearts. These data now provide a link between IFN, IFN-stimulated genes, and the inflammatory and possibly fibrosing components of CHB, positioning Siglec-1-positive macrophages as integral to the process.

Aortitis masquerading as intramural hematoma: When to observe, when to operate? A case report.

Radiologic evidence of aortic disease is not always consistent with the diagnosis. With a lack of accompanying symptoms or with an atypical presentation, diagnosis, and management of aortic pathology rely greatly on imaging techniques. We report the case of a 58-year-old female who presented with incidental radiographic findings consistent with a type A aortic intramural hematoma and a vague left-sided chest discomfort. After follow-up, imaging was consistent with disease progression and hematoma expansion; the affected segment was resected and pathology reported lymphoplasmacytic aortitis as the underlying etiology of the imaging findings rather than an intramural hematoma. The patient lacked symptoms or serology consistent with the rheumatologic disease, and the postoperative course was uneventful. The management of a suspected ascending intramural hematoma is controversial, especially when the patient presents with atypical signs and symptoms. Features of disease progression may warrant urgent surgical intervention.

Toward the human cellular microRNAome.

MicroRNAs are short RNAs that serve as regulators of gene expression and are essential components of normal development as well as modulators of disease. MicroRNAs generally act cell-autonomously, and thus their localization to specific cell types is needed to guide our understanding of microRNA activity. Current tissue-level data have caused considerable confusion, and comprehensive cell-level data do not yet exist. Here, we establish the landscape of human cell-specific microRNA expression. This project evaluated 8 billion small RNA-seq reads from 46 primary cell types, 42 cancer or immortalized cell lines, and 26 tissues. It identified both specific and ubiquitous patterns of expression that strongly correlate with adjacent superenhancer activity. Analysis of unaligned RNA reads uncovered 207 unknown minor strand (passenger) microRNAs of known microRNA loci and 495 novel putative microRNA loci. Although cancer cell lines generally recapitulated the expression patterns of matched primary cells, their isomiR sequence families exhibited increased disorder, suggesting DROSHA- and DICER1-dependent microRNA processing variability. Cell-specific patterns of microRNA expression were used to de-convolute variable cellular composition of colon and adipose tissue samples, highlighting one use of these cell-specific microRNA expression data. Characterization of cellular microRNA expression across a wide variety of cell types provides a new understanding of this critical regulatory RNA species.

Integrated Transcriptomic and Proteomic Analysis of Primary Human Umbilical Vein Endothelial Cells.

Understanding the molecular profile of every human cell type is essential for understanding its role in normal physiology and disease. Technological advancements in DNA sequencing, mass spectrometry, and computational methods allow us to carry out multiomics analyses although such approaches are not routine yet. Human umbilical vein endothelial cells (HUVECs) are a widely used model system to study pathological and physiological processes associated with the cardiovascular system. In this study, next-generation sequencing and high-resolution mass spectrometry to profile the transcriptome and proteome of primary HUVECs is employed. Analysis of 145 million paired-end reads from next-generation sequencing confirmed expression of 12 186 protein-coding genes (FPKM ≥0.1), 439 novel long non-coding RNAs, and revealed 6089 novel isoforms that were not annotated in GENCODE. Proteomics analysis identifies 6477 proteins including confirmation of N-termini for 1091 proteins, isoforms for 149 proteins, and 1034 phosphosites. A database search to specifically identify other post-translational modifications provide evidence for a number of modification sites on 117 proteins which include ubiquitylation, lysine acetylation, and mono-, di- and tri-methylation events. Evidence for 11 "missing proteins," which are proteins for which there was insufficient or no protein level evidence, is provided. Peptides supporting missing protein and novel events are validated by comparison of MS/MS fragmentation patterns with synthetic peptides. Finally, 245 variant peptides derived from 207 expressed proteins in addition to alternate translational start sites for seven proteins and evidence for novel proteoforms for five proteins resulting from alternative splicing are identified. Overall, it is believed that the integrated approach employed in this study is widely applicable to study any primary cell type for deeper molecular characterization.