An integrated organoid omics map extends modeling potential of kidney disease.

Kidney organoids are a promising model to study kidney disease, but their use is constrained by limited knowledge of their functional protein expression profile. Here, we define the organoid proteome and transcriptome trajectories over culture duration and upon exposure to TNFα, a cytokine stressor. Older organoids increase deposition of extracellular matrix but decrease expression of glomerular proteins. Single cell transcriptome integration reveals that most proteome changes localize to podocytes, tubular and stromal cells. TNFα treatment of organoids results in 322 differentially expressed proteins, including cytokines and complement components. Transcript expression of these 322 proteins is significantly higher in individuals with poorer clinical outcomes in proteinuric kidney disease. Key TNFα-associated protein (C3 and VCAM1) expression is increased in both human tubular and organoid kidney cell populations, highlighting the potential for organoids to advance biomarker development. By integrating kidney organoid omic layers, incorporating a disease-relevant cytokine stressor and comparing with human data, we provide crucial evidence for the functional relevance of the kidney organoid model to human kidney disease.

SGLT2 inhibitors mitigate kidney tubular metabolic and mTORC1 perturbations in youth-onset type 2 diabetes.

The molecular mechanisms of sodium-glucose cotransporter-2 (SGLT2) inhibitors (SGLT2i) remain incompletely understood. Single-cell RNA sequencing and morphometric data were collected from research kidney biopsies donated by young persons with type 2 diabetes (T2D), aged 12 to 21 years, and healthy controls (HCs). Participants with T2D were obese and had higher estimated glomerular filtration rates and mesangial and glomerular volumes than HCs. Ten T2D participants had been prescribed SGLT2i (T2Di[+]) and 6 not (T2Di[-]). Transcriptional profiles showed SGLT2 expression exclusively in the proximal tubular (PT) cluster with highest expression in T2Di(-) patients. However, transcriptional alterations with SGLT2i treatment were seen across nephron segments, particularly in the distal nephron. SGLT2i treatment was associated with suppression of transcripts in the glycolysis, gluconeogenesis, and tricarboxylic acid cycle pathways in PT, but had the opposite effect in thick ascending limb. Transcripts in the energy-sensitive mTORC1-signaling pathway returned toward HC levels in all tubular segments in T2Di(+), consistent with a diabetes mouse model treated with SGLT2i. Decreased levels of phosphorylated S6 protein in proximal and distal tubules in T2Di(+) patients confirmed changes in mTORC1 pathway activity. We propose that SGLT2i treatment benefits the kidneys by mitigating diabetes-induced metabolic perturbations via suppression of mTORC1 signaling in kidney tubules.

Precision nephrology identified tumor necrosis factor activation variability in minimal change disease and focal segmental glomerulosclerosis.

The diagnosis of nephrotic syndrome relies on clinical presentation and descriptive patterns of injury on kidney biopsies, but not specific to underlying pathobiology. Consequently, there are variable rates of progression and response to therapy within diagnoses. Here, an unbiased transcriptomic-driven approach was used to identify molecular pathways which are shared by subgroups of patients with either minimal change disease (MCD) or focal segmental glomerulosclerosis (FSGS). Kidney tissue transcriptomic profile-based clustering identified three patient subgroups with shared molecular signatures across independent, North American, European, and African cohorts. One subgroup had significantly greater disease progression (Hazard Ratio 5.2) which persisted after adjusting for diagnosis and clinical measures (Hazard Ratio 3.8). Inclusion in this subgroup was retained even when clustering was limited to those with less than 25% interstitial fibrosis. The molecular profile of this subgroup was largely consistent with tumor necrosis factor (TNF) pathway activation. Two TNF pathway urine markers were identified, tissue inhibitor of metalloproteinases-1 (TIMP-1) and monocyte chemoattractant protein-1 (MCP-1), that could be used to predict an individual's TNF pathway activation score. Kidney organoids and single-nucleus RNA-sequencing of participant kidney biopsies, validated TNF-dependent increases in pathway activation score, transcript and protein levels of TIMP-1 and MCP-1, in resident kidney cells. Thus, molecular profiling identified a subgroup of patients with either MCD or FSGS who shared kidney TNF pathway activation and poor outcomes. A clinical trial testing targeted therapies in patients selected using urinary markers of TNF pathway activation is ongoing.

A single natural RNA modification can destabilize a U•A-T-rich RNA•DNA-DNA triple helix.

Recent studies suggest noncoding RNAs interact with genomic DNA, forming RNA•DNA-DNA triple helices, as a mechanism to regulate transcription. One way cells could regulate the formation of these triple helices is through RNA modifications. With over 140 naturally occurring RNA modifications, we hypothesize that some modifications stabilize RNA•DNA-DNA triple helices while others destabilize them. Here, we focus on a pyrimidine-motif triple helix composed of canonical U•A-T and C•G-C base triples. We employed electrophoretic mobility shift assays and microscale thermophoresis to examine how 11 different RNA modifications at a single position in an RNA•DNA-DNA triple helix affect stability: 5-methylcytidine (m5C), 5-methyluridine (m5U or rT), 3-methyluridine (m3U), pseudouridine (Ψ), 4-thiouridine (s4U), N 6-methyladenosine (m6A), inosine (I), and each nucleobase with 2'-O-methylation (Nm). Compared to the unmodified U•A-T base triple, some modifications have no significant change in stability (Um•A-T), some have ∼2.5-fold decreases in stability (m5U•A-T, Ψ•A-T, and s4U•A-T), and some completely disrupt triple helix formation (m3U•A-T). To identify potential biological examples of RNA•DNA-DNA triple helices controlled by an RNA modification, we searched RMVar, a database for RNA modifications mapped at single-nucleotide resolution, for lncRNAs containing an RNA modification within a pyrimidine-rich sequence. Using electrophoretic mobility shift assays, the binding of DNA-DNA to a 22-mer segment of human lncRNA Al157886.1 was destabilized by ∼1.7-fold with the substitution of m5C at known m5C sites. Therefore, the formation and stability of cellular RNA•DNA-DNA triple helices could be influenced by RNA modifications.

Secondary Structural Model of MALAT1 Becomes Unstructured in Chronic Myeloid Leukemia and Undergoes Structural Rearrangement in Cervical Cancer.

Long noncoding RNAs (lncRNAs) influence cellular function through binding events that often depend on the lncRNA secondary structure. One such lncRNA, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), is upregulated in many cancer types and has a myriad of protein- and miRNA-binding sites. Recently, a secondary structural model of MALAT1 in noncancerous cells was proposed to form 194 hairpins and 13 pseudoknots. That study postulated that, in cancer cells, the MALAT1 structure likely varies, thereby influencing cancer progression. This work analyzes how that structural model is expected to change in K562 cells, which originated from a patient with chronic myeloid leukemia (CML), and in HeLa cells, which originated from a patient with cervical cancer. Dimethyl sulfate-sequencing (DMS-Seq) data from K562 cells and psoralen analysis of RNA interactions and structure (PARIS) data from HeLa cells were compared to the working structural model of MALAT1 in noncancerous cells to identify sites that likely undergo structural alterations. MALAT1 in K562 cells is predicted to become more unstructured, with almost 60% of examined hairpins in noncancerous cells losing at least half of their base pairings. Conversely, MALAT1 in HeLa cells is predicted to largely maintain its structure, undergoing 18 novel structural rearrangements. Moreover, 50 validated miRNA-binding sites are affected by putative secondary structural changes in both cancer types, such as miR-217 in K562 cells and miR-20a in HeLa cells. Structural changes unique to K562 cells and HeLa cells provide new mechanistic leads into how the structure of MALAT1 may mediate cancer in a cell-type specific manner.

Naturally occurring modified ribonucleosides.

The chemical identity of RNA molecules beyond the four standard ribonucleosides has fascinated scientists since pseudouridine was characterized as the "fifth" ribonucleotide in 1951. Since then, the ever-increasing number and complexity of modified ribonucleosides have been found in viruses and throughout all three domains of life. Such modifications can be as simple as methylations, hydroxylations, or thiolations, complex as ring closures, glycosylations, acylations, or aminoacylations, or unusual as the incorporation of selenium. While initially found in transfer and ribosomal RNAs, modifications also exist in messenger RNAs and noncoding RNAs. Modifications have profound cellular outcomes at various levels, such as altering RNA structure or being essential for cell survival or organism viability. The aberrant presence or absence of RNA modifications can lead to human disease, ranging from cancer to various metabolic and developmental illnesses such as Hoyeraal-Hreidarsson syndrome, Bowen-Conradi syndrome, or Williams-Beuren syndrome. In this review article, we summarize the characterization of all 143 currently known modified ribonucleosides by describing their taxonomic distributions, the enzymes that generate the modifications, and any implications in cellular processes, RNA structure, and disease. We also highlight areas of active research, such as specific RNAs that contain a particular type of modification as well as methodologies used to identify novel RNA modifications. This article is categorized under: RNA Processing > RNA Editing and Modification.

Secondary Structural Model of Human MALAT1 Reveals Multiple Structure-Function Relationships.

Human metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) is an abundant nuclear-localized long noncoding RNA (lncRNA) that has significant roles in cancer. While the interacting partners and evolutionary sequence conservation of MALAT1 have been examined, much of the structure of MALAT1 is unknown. Here, we propose a hypothetical secondary structural model for 8425 nucleotides of human MALAT1 using three experimental datasets that probed RNA structures in vitro and in various human cell lines. Our model indicates that approximately half of human MALAT1 is structured, forming 194 helices, 13 pseudoknots, five structured tetraloops, nine structured internal loops, and 13 intramolecular long-range interactions that give rise to several multiway junctions. Evolutionary conservation and covariation analyses support 153 of 194 helices in 51 mammalian MALAT1 homologs and 42 of 194 helices in 53 vertebrate MALAT1 homologs, thereby identifying an evolutionarily conserved core that likely has important functional roles in mammals and vertebrates. Data mining revealed that RNA modifications, somatic cancer-associated mutations, and single-nucleotide polymorphisms may induce structural rearrangements that sequester or expose binding sites for several cancer-associated microRNAs. Our findings reveal new mechanistic leads into the roles of MALAT1 by identifying several intriguing structure-function relationships in which the dynamic structure of MALAT1 underlies its biological functions.

Stability of an RNA•DNA-DNA triple helix depends on base triplet composition and length of the RNA third strand.

Recent studies suggest noncoding RNAs interact with genomic DNA, forming an RNA•DNA-DNA triple helix that regulates gene expression. However, base triplet composition of pyrimidine motif RNA•DNA-DNA triple helices is not well understood beyond the canonical U•A-T and C•G-C base triplets. Using native gel-shift assays, the relative stability of 16 different base triplets at a single position, Z•X-Y (where Z = C, U, A, G and X-Y = A-T, G-C, T-A, C-G), in an RNA•DNA-DNA triple helix was determined. The canonical U•A-T and C•G-C base triplets were the most stable, while three non-canonical base triplets completely disrupted triple-helix formation. We further show that our RNA•DNA-DNA triple helix can tolerate up to two consecutive non-canonical A•G-C base triplets. Additionally, the RNA third strand must be at least 19 nucleotides to form an RNA•DNA-DNA triple helix but increasing the length to 27 nucleotides does not increase stability. The relative stability of 16 different base triplets in DNA•DNA-DNA and RNA•RNA-RNA triple helices was distinctly different from those in RNA•DNA-DNA triple helices, showing that base triplet stability depends on strand composition being DNA and/or RNA. Multiple factors influence the stability of triple helices, emphasizing the importance of experimentally validating formation of computationally predicted triple helices.

Riboswitch diversity and distribution.

Riboswitches are commonly used by bacteria to detect a variety of metabolites and ions to regulate gene expression. To date, nearly 40 different classes of riboswitches have been discovered, experimentally validated, and modeled at atomic resolution in complex with their cognate ligands. The research findings produced since the first riboswitch validation reports in 2002 reveal that these noncoding RNA domains exploit many different structural features to create binding pockets that are extremely selective for their target ligands. Some riboswitch classes are very common and are present in bacteria from nearly all lineages, whereas others are exceedingly rare and appear in only a few species whose DNA has been sequenced. Presented herein are the consensus sequences, structural models, and phylogenetic distributions for all validated riboswitch classes. Based on our findings, we predict that there are potentially many thousands of distinct bacterial riboswitch classes remaining to be discovered, but that the rarity of individual undiscovered classes will make it increasingly difficult to find additional examples of this RNA-based sensory and gene control mechanism.

Singlet glycine riboswitches bind ligand as well as tandem riboswitches.

The glycine riboswitch often occurs in a tandem architecture, with two ligand-binding domains (aptamers) followed by a single expression platform. Based on previous observations, we hypothesized that "singlet" versions of the glycine riboswitch, which contain only one aptamer domain, are able to bind glycine if appropriate structural contacts are maintained. An initial alignment of 17 putative singlet riboswitches indicated that the single consensus aptamer domain is flanked by a conserved peripheral stem-loop structure. These singlets were sorted into two subtypes based on whether the active aptamer domain precedes or follows the peripheral stem-loop, and an example of each subtype of singlet riboswitch was characterized biochemically. The singlets possess glycine-binding affinities comparable to those of previously published tandem examples, and the conserved peripheral domains form A-minor interactions with the single aptamer domain that are necessary for ligand-binding activity. Analysis of sequenced genomes identified a significant number of singlet glycine riboswitches. Based on these observations, we propose an expanded model for glycine riboswitch gene control that includes singlet and tandem architectures.

Control of bacterial exoelectrogenesis by c-AMP-GMP.

Major changes in bacterial physiology including biofilm and spore formation involve signaling by the cyclic dinucleotides c-di-GMP and c-di-AMP. Recently, another second messenger dinucleotide, c-AMP-GMP, was found to control chemotaxis and colonization by Vibrio cholerae. We have identified a superregulon of genes controlled by c-AMP-GMP in numerous Deltaproteobacteria, including Geobacter species that use extracellular insoluble metal oxides as terminal electron acceptors. This exoelectrogenic process has been studied for its possible utility in energy production and bioremediation. Many genes involved in adhesion, pilin formation, and others that are important for exoelectrogenesis are controlled by members of a variant riboswitch class that selectively bind c-AMP-GMP. These RNAs constitute, to our knowledge, the first known specific receptors for c-AMP-GMP and reveal that this molecule is used by many bacteria to control specialized physiological processes.

Structural, functional, and taxonomic diversity of three preQ1 riboswitch classes.

Previously, two riboswitch classes have been identified that sense and respond to the hypermodified nucleobase called prequeuosine1 (preQ1). The enormous expansion of available genomic DNA sequence data creates new opportunities to identify additional representatives of the known riboswitch classes and to discover novel classes. We conducted bioinformatics searches on microbial genomic DNA data sets to discover numerous additional examples belonging to the two previously known riboswitch classes for preQ1 (classes preQ1-I and preQ1-II), including some structural variants that further restrict ligand specificity. Additionally, we discovered a third preQ1-binding riboswitch class (preQ1-III) that is structurally distinct from previously known classes. These findings demonstrate that numerous organisms monitor the concentrations of this modified nucleobase by exploiting one or more riboswitch classes for this widespread compound.

Mechanism and distribution of glmS ribozymes.

Among the nine classes of ribozymes that have been experimentally validated to date is the metabolite-responsive self-cleaving ribozyme called glmS. This RNA is almost exclusively located in the 5'-untranslated region of bacterial mRNAs that code for the production of GlmS proteins, which catalyze the synthesis of the aminosugar glucosamine-6-phosphate (GlcN6P). Each glmS ribozyme forms a conserved catalytic core that selectively binds GlcN6P and uses this metabolite as a cofactor to promote ribozyme self-cleavage. Metabolite-induced self-cleavage results in down-regulation of glmS gene expression, and thus the ribozyme functions as a key riboswitch component to permit feedback regulation of GlcN6P levels. Representatives of glmS ribozymes also serve as excellent experimental models to elucidate how RNAs fold to recognize small molecule ligands and promote chemical transformations.

An expanded collection and refined consensus model of glmS ribozymes.

Self-cleaving glmS ribozymes selectively bind glucosamine-6-phosphate (GlcN6P) and use this metabolite as a cofactor to promote self-cleavage by internal phosphoester transfer. Representatives of the glmS ribozyme class are found in Gram-positive bacteria where they reside in the 5' untranslated regions (UTRs) of glmS messenger RNAs that code for the essential enzyme L-glutamine:D-fructose-6-phosphate aminotransferase. By using comparative sequence analyses, we have expanded the number of glmS ribozyme representatives from 160 to 463. All but two glmS ribozymes are present in glmS mRNAs and most exhibit striking uniformity in sequence and structure, which are features that make representatives attractive targets for antibacterial drug development. However, our discovery of rare variants broadens the consensus sequence and structure model. For example, in the Deinococcus-Thermus phylum, several structural variants exist that carry additional stems within the catalytic core and changes to the architecture of core-supporting substructures. These findings reveal that glmS ribozymes have a broader phylogenetic distribution than previously known and suggest that additional rare structural variants may remain to be discovered.

Large-scale analysis of yeast filamentous growth by systematic gene disruption and overexpression.

Under certain conditions of nutrient stress, the budding yeast Saccharomyces cerevisiae initiates a striking developmental transition to a filamentous form of growth, resembling developmental transitions required for virulence in closely related pathogenic fungi. In yeast, filamentous growth involves known mitogen-activated protein kinase and protein kinase A signaling modules, but the full scope of this extensive filamentous response has not been delineated. Accordingly, we have undertaken the first systematic gene disruption and overexpression analysis of yeast filamentous growth. Standard laboratory strains of yeast are nonfilamentous; thus, we constructed a unique set of reagents in the filamentous Sigma1278b strain, encompassing 3627 integrated transposon insertion alleles and 2043 overexpression constructs. Collectively, we analyzed 4528 yeast genes with these reagents and identified 487 genes conferring mutant filamentous phenotypes upon transposon insertion and/or gene overexpression. Using a fluorescent protein reporter integrated at the MUC1 locus, we further assayed each filamentous growth mutant for aberrant protein levels of the key flocculence factor Muc1p. Our results indicate a variety of genes and pathways affecting filamentous growth. In total, this filamentous growth gene set represents a wealth of yeast biology, highlighting 84 genes of uncharacterized function and an underappreciated role for the mitochondrial retrograde signaling pathway as an inhibitor of filamentous growth.