Quadruped Gait and Regulation of Apoptotic Factors in Tibiofemoral Joints following Intra-Articular rhPRG4 Injection in Prg4 Null Mice.

Camptodactyly-arthropathy-coxa vara-pericarditis (CACP) syndrome leads to diarthrodial joint arthropathy and is caused by the absence of lubricin (proteoglycan 4-PRG4), a surface-active mucinous glycoprotein responsible for lubricating articular cartilage. In this study, mice lacking the orthologous gene Prg4 served as a model that recapitulates the destructive arthrosis that involves biofouling of cartilage by serum proteins in lieu of Prg4. This study hypothesized that Prg4-deficient mice would demonstrate a quadruped gait change and decreased markers of mitochondrial dyscrasia, following intra-articular injection of both hindlimbs with recombinant human PRG4 (rhPRG4). Prg4-/- (N = 44) mice of both sexes were injected with rhPRG4 and gait alterations were studied at post-injection day 3 and 6, before joints were harvested for immunohistochemistry for caspase-3 activation. Increased stance and propulsion was shown at 3 days post-injection in male mice. There were significantly fewer caspase-3-positive chondrocytes in tibiofemoral cartilage from rhPRG4-injected mice. The mitochondrial gene Mt-tn, and myosin heavy (Myh7) and light chains (Myl2 and Myl3), known to play a cytoskeletal stabilizing role, were significantly upregulated in both sexes (RNA-Seq) following IA rhPRG4. Chondrocyte mitochondrial dyscrasias attributable to the arthrosis in CACP may be mitigated by IA rhPRG4. In a supporting in vitro crystal microbalance experiment, molecular fouling by albumin did not block the surface activity of rhPRG4.

Future directions for high-throughput splicing assays in precision medicine.

Classification of variants of unknown significance is a challenging technical problem in clinical genetics. As up to one-third of disease-causing mutations are thought to affect pre-mRNA splicing, it is important to accurately classify splicing mutations in patient sequencing data. Several consortia and healthcare systems have conducted large-scale patient sequencing studies, which discover novel variants faster than they can be classified. Here, we compare the advantages and limitations of several high-throughput splicing assays aimed at mitigating this bottleneck, and describe a data set of ~5,000 variants that we analyzed using our Massively Parallel Splicing Assay (MaPSy). The Critical Assessment of Genome Interpretation group (CAGI) organized a challenge, in which participants submitted machine learning models to predict the splicing effects of variants in this data set. We discuss the winning submission of the challenge (MMSplice) which outperformed existing software. Finally, we highlight methods to overcome the limitations of MaPSy and similar assays, such as tissue-specific splicing, the effect of surrounding sequence context, classifying intronic variants, synthesizing large exons, and amplifying complex libraries of minigene species. Further development of these assays will greatly benefit the field of clinical genetics, which lack high-throughput methods for variant interpretation.

The effects of structure on pre-mRNA processing and stability.

Pre-mRNA molecules can form a variety of structures, and both secondary and tertiary structures have important effects on processing, function and stability of these molecules. The prediction of RNA secondary structure is a challenging problem and various algorithms that use minimum free energy, maximum expected accuracy and comparative evolutionary based methods have been developed to predict secondary structures. However, these tools are not perfect, and this remains an active area of research. The secondary structure of pre-mRNA molecules can have an enhancing or inhibitory effect on pre-mRNA splicing. An example of enhancing structure can be found in a novel class of introns in zebrafish. About 10% of zebrafish genes contain a structured intron that forms a bridging hairpin that enforces correct splice site pairing. Negative examples of splicing include local structures around splice sites that decrease splicing efficiency and potentially cause mis-splicing leading to disease. Splicing mutations are a frequent cause of hereditary disease. The transcripts of disease genes are significantly more structured around the splice sites, and point mutations that increase the local structure often cause splicing disruptions. Post-splicing, RNA secondary structure can also affect the stability of the spliced intron and regulatory RNA interference pathway intermediates, such as pre-microRNAs. Additionally, RNA secondary structure has important roles in the innate immune defense against viruses. Finally, tertiary structure can also play a large role in pre-mRNA splicing. One example is the G-quadruplex structure, which, similar to secondary structure, can either enhance or inhibit splicing through mechanisms such as creating or obscuring RNA binding protein sites.

The unusual gene architecture of polyubiquitin is created by dual-specific splice sites.

BACKGROUND: The removal of introns occurs through the splicing of a 5' splice site (5'ss) with a 3' splice site (3'ss). These two elements are recognized by distinct components of the spliceosome. However, introns in higher eukaryotes contain many matches to the 5' and 3' splice-site motifs that are presumed not to be used. RESULTS: Here, we find that many of these sites can be used. We also find occurrences of the AGGT motif that can function as either a 5'ss or a 3'ss-previously referred to as dual-specific splice sites (DSSs)-within introns. Analysis of the Sequence Read Archive reveals a 3.1-fold enrichment of DSSs relative to expectation, implying synergy between the ability to function as a 5'ss and 3'ss. Despite this suggested mechanistic advantage, DSSs are 2.7- and 4.7-fold underrepresented in annotated 5' and 3' splice sites. A curious exception is the polyubiquitin gene UBC, which contains a tandem array of DSSs that precisely delimit the boundary of each ubiquitin monomer. The resulting isoforms splice stochastically to include a variable number of ubiquitin monomers. We found no evidence of tissue-specific or feedback regulation but note the 8.4-fold enrichment of DSS-spliced introns in tandem repeat genes suggests a driving role in the evolution of genes like UBC. CONCLUSIONS: We find an excess of unannotated splice sites and the utilization of DSSs in tandem repeats supports the role of splicing in gene evolution. These findings enhance our understanding of the diverse and complex nature of the splicing process.

Application of single cell multiomics points to changes in chromatin accessibility near calcitonin receptor like receptor and a possible role for adrenomedullin in the post-shock lung.

Introduction: Acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) is a commonly occurring sequelae of traumatic injury resulting from indirect insults like hypovolemic shock and/or extrapulmonary sepsis. The high lethality rate associated with these pathologies outlines the importance of clarifying the "priming" effects seen in the post-shock lung microenvironment, which are understood to bring about a dysregulated or overt immune response when triggered by a secondary systemic infectious/septic challenge culminating in ALI. In this pilot project, we test the hypothesis that application of a single cell multiomics approach can elucidate novel phenotype specific pathways potentially contributing to shock-induced ALI/ARDS. Methods: Hypovolemic shock was induced in C57BL/6 (wild-type), PD-1, PD-L1, or VISTA gene deficient male mice, 8-12 weeks old. Wild-type sham surgeries function as negative controls. A total of 24-h post-shock rodents were sacrificed, their lungs harvested and sectioned, with pools prepared from 2 mice per background, and flash frozen on liquid nitrogen. N = 2 biological replicates (representing 4 mice total) were achieved for all treatment groups across genetic backgrounds. Samples were received by the Boas Center for Genomics and Human Genetics, where single cell multiomics libraries were prepared for RNA/ATAC sequencing. The analysis pipeline Cell Ranger ARC was implemented to attain feature linkage assessments across genes of interest. Results: Sham (pre-shock) results suggest high chromatin accessibility around calcitonin receptor like receptor (CALCRL) across cellular phenotypes with 17 and 18 feature links, exhibiting positive correlation with gene expression between biological replicates. Similarity between both sample chromatin profiles/linkage arcs is evident. Post-shock wild-type accessibility is starkly reduced across replicates where the number of feature links drops to 1 and 3, again presenting similar replicate profiles. Samples from shocked gene deficient backgrounds displayed high accessibility and similar profiles to the pre-shock lung microenvironment. Conclusion: High pre-shock availability of DNA segments and their positive correlation with CALCRL gene expression suggests an apparent regulatory capacity on transcription. Post-shock gene deficient chromatin profiles presented similar results to that of pre-shock wild-type samples, suggesting an influence on CALCRL accessibility. Key changes illustrated in the pre-ALI context of shock may allow for additional resolution of "priming" and "cellular pre-activation/pre-disposition" processes within the lung microenvironment.

Identification and mechanistic basis of non-ACE2 blocking neutralizing antibodies from COVID-19 patients with deep RNA sequencing and molecular dynamics simulations.

Variants of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continue to cause disease and impair the effectiveness of treatments. The therapeutic potential of convergent neutralizing antibodies (NAbs) from fully recovered patients has been explored in several early stages of novel drugs. Here, we identified initially elicited NAbs (Ig Heavy, Ig lambda, Ig kappa) in response to COVID-19 infection in patients admitted to the intensive care unit at a single center with deep RNA sequencing (>100 million reads) of peripheral blood as a diagnostic tool for predicting the severity of the disease and as a means to pinpoint specific compensatory NAb treatments. Clinical data were prospectively collected at multiple time points during ICU admission, and amino acid sequences for the NAb CDR3 segments were identified. Patients who survived severe COVID-19 had significantly more of a Class 3 antibody (C135) to SARS-CoV-2 compared to non-survivors (15059.4 vs. 1412.7, p = 0.016). In addition to highlighting the utility of RNA sequencing in revealing unique NAb profiles in COVID-19 patients with different outcomes, we provided a physical basis for our findings via atomistic modeling combined with molecular dynamics simulations. We established the interactions of the Class 3 NAb C135 with the SARS-CoV-2 spike protein, proposing a mechanistic basis for inhibition via multiple conformations that can effectively prevent ACE2 from binding to the spike protein, despite C135 not directly blocking the ACE2 binding motif. Overall, we demonstrate that deep RNA sequencing combined with structural modeling offers the new potential to identify and understand novel therapeutic(s) NAbs in individuals lacking certain immune responses due to their poor endogenous production. Our results suggest a possible window of opportunity for administration of such NAbs when their full sequence becomes available. A method involving rapid deep RNA sequencing of patients infected with SARS-CoV-2 or its variants at the earliest infection time could help to develop personalized treatments using the identified specific NAbs.

Deep RNA sequencing of intensive care unit patients with COVID-19.

COVID-19 has impacted millions of patients across the world. Molecular testing occurring now identifies the presence of the virus at the sampling site: nasopharynx, nares, or oral cavity. RNA sequencing has the potential to establish both the presence of the virus and define the host's response in COVID-19. Single center, prospective study of patients with COVID-19 admitted to the intensive care unit where deep RNA sequencing (> 100 million reads) of peripheral blood with computational biology analysis was done. All patients had positive SARS-CoV-2 PCR. Clinical data was prospectively collected. We enrolled fifteen patients at a single hospital. Patients were critically ill with a mortality of 47% and 67% were on a ventilator. All the patients had the SARS-CoV-2 RNA identified in the blood in addition to RNA from other viruses, bacteria, and archaea. The expression of many immune modulating genes, including PD-L1 and PD-L2, were significantly different in patients who died from COVID-19. Some proteins were influenced by alternative transcription and splicing events, as seen in HLA-C, HLA-E, NRP1 and NRP2. Entropy calculated from alternative RNA splicing and transcription start/end predicted mortality in these patients. Current upper respiratory tract testing for COVID-19 only determines if the virus is present. Deep RNA sequencing with appropriate computational biology may provide important prognostic information and point to therapeutic foci to be precisely targeted in future studies.

Deep RNA Sequencing of Intensive Care Unit Patients with COVID-19.

PURPOSE: COVID-19 has impacted millions of patients across the world. Molecular testing occurring now identifies the presence of the virus at the sampling site: nasopharynx, nares, or oral cavity. RNA sequencing has the potential to establish both the presence of the virus and define the host's response in COVID-19. METHODS: Single center, prospective study of patients with COVID-19 admitted to the intensive care unit where deep RNA sequencing (>100 million reads) of peripheral blood with computational biology analysis was done. All patients had positive SARS-CoV-2 PCR. Clinical data was prospectively collected. RESULTS: We enrolled fifteen patients at a single hospital. Patients were critically ill with a mortality of 47% and 67% were on a ventilator. All the patients had the SARS-CoV-2 RNA identified in the blood in addition to RNA from other viruses, bacteria, and archaea. The expression of many immune modulating genes, including PD-L1 and PD-L2, were significantly different in patients who died from COVID-19. Some proteins were influenced by alternative transcription and splicing events, as seen in HLA-C, HLA-E, NRP1 and NRP2. Entropy calculated from alternative RNA splicing and transcription start/end predicted mortality in these patients. CONCLUSIONS: Current upper respiratory tract testing for COVID-19 only determines if the virus is present. Deep RNA sequencing with appropriate computational biology may provide important prognostic information and point to therapeutic foci to be precisely targeted in future studies. TAKE HOME MESSAGE: Deep RNA sequencing provides a novel diagnostic tool for critically ill patients. Among ICU patients with COVID-19, RNA sequencings can identify gene expression, pathogens (including SARS-CoV-2), and can predict mortality. TWEET: Deep RNA sequencing is a novel technology that can assist in the care of critically ill COVID-19 patients & can be applied to other disease.

RNA-Binding Proteins: Splicing Factors and Disease.

Pre-mRNA splicing is mediated by interactions of the Core Spliceosome and an array of accessory RNA binding proteins with cis-sequence elements. Splicing is a major regulatory component in higher eukaryotes. Disruptions in splicing are a major contributor to human disease. One in three hereditary disease alleles are believed to cause aberrant splicing. Hereditary disease alleles can alter splicing by disrupting a splicing element, creating a toxic RNA, or affecting splicing factors. One of the challenges of medical genetics is identifying causal variants from the thousands of possibilities discovered in a clinical sequencing experiment. Here we review the basic biochemistry of splicing, the mechanisms of splicing mutations, the methods for identifying splicing mutants, and the potential of therapeutic interventions.