Search

Genome-wide association studies (GWAS) have identified loci for kidney disease, but the causal variants, genes, and pathways remain unknown. Here we identify two kidney disease genes Dipeptidase 1 (DPEP1) and Charged Multivesicular Body Protein 1 A (CHMP1A) via the triangulation of kidney function GWAS, human kidney expression, and methylation quantitative trait loci. Using single-cell chromatin accessibility and genome editing, we fine map the region that controls the expression of both genes. Mouse genetic models demonstrate the causal roles of both genes in kidney disease. Cellular studies indicate that both Dpep1 and Chmp1a are important regulators of a single pathway, ferroptosis and lead to kidney disease development via altering cellular iron trafficking., Identifying causal variants and genes is an essential step in interpreting GWAS loci. Here, the authors investigate a kidney disease GWAS locus with functional genomics data, CRISPR editing and mouse experiments to identify DPEP1 and CHMP1A as putative kidney disease genes via ferroptosis.

The functional interpretation of genome-wide association studies (GWAS) is challenging due to the cell-type-dependent influences of genetic variants. Here, we generated comprehensive maps of expression quantitative trait loci (eQTLs) for 659 microdissected human kidney samples and identified cell-type-eQTLs by mapping interactions between cell type abundances and genotypes. By partitioning heritability using stratified linkage disequilibrium score regression to integrate GWAS with single-cell RNA sequencing and single-nucleus assay for transposase-accessible chromatin with high-throughput sequencing data, we prioritized proximal tubules for kidney function and endothelial cells and distal tubule segments for blood pressure pathogenesis. Bayesian colocalization analysis nominated more than 200 genes for kidney function and hypertension. Our study clarifies the mechanism of commonly used antihypertensive and renal-protective drugs and identifies drug repurposing opportunities for kidney disease.

Outbreaks of mass mortalities occurred in Macrobrachium rosenbergii farms in Gaoyou county, Jiangsu Province of China. The bacterial isolates from M. rosenbergii exhibited the same phenotypic traits and biochemical characteristics, and were identified as non-O1/O139 Vibrio cholerae according to biochemical characteristics and molecular identification. In challenge test, M. rosenbergii infected with non-O1/O139 V. cholerae GXFL1-4 developed similar pathological signs to the naturally diseased prawns, and LD 50 of the strain to M. rosenbergii was 4.5 × 10 6 CFU/mL at 96 h post-infection. Histopathological analysis revealed that hepatopancreas and intestines of diseased M. rosenbergii exhibited obvious inflammatory responses to non-O1/O139 V. cholerae infection. Detection virulence factors of the strain GXFL1-4 showed that the bacteria produced caseinase, lipase, amylase, lecithinase and hemolysin, and carried toxR , hlyA , ompW , ompU , hap , rtxA and rtxC virulence related genes, supporting the strong virulence to M. rosenbergii . Additionally, the immune related gene expression in M. rosenbergii evaluated by qRT-PCR analysis showed that HSP70 , Crustin , Lysozyme , TRL1 , ALF1 , Lectin , Peroxinectin , proPO and SOD immune related genes were significantly up-regulated at 6 and 12 h after infection with GXFL1-4. The results of our study suggested that non-O1/O139 V. cholerae was an etiological element in the mass mortalities of M. rosenbergii and this study provided preliminary insights into the diversity in the immune response of M. rosenbergii to the bacterial invasion.

Bacillus thuringiensis BMB181 (Bt BMB181), with high melanin production, is an ideal candidate for industrial scale production of light-stable insecticides. However, its melanogenic pathways remain unclear. In the present study, we demonstrated that Bt BMB181 failed to produce melanin after treatment with mesotrione, an inhibitor of 4-hydroxyphenylpyruvate dioxygenase in the homogentisic acid pathway of melanin synthesis. Heterologous expression experiments suggested that homogentisate-1,2-dioxygenase (HmgA) in Bt BMB171 functions normally, yet HmgA in Bt BMB181 had lost its activity, at least partly. Using the CRISPR-Cas9 system, the hmgA gene in Bt BMB171 was knocked out and the mutant strain gained the ability to produce melanin. Furthermore, the complemented strain reverted to the wild-type phenotype. In addition, overexpression of its own hmgA gene in Bt BMB181 also resulted in failure to produce the pigment. BLAST results indicated that the amino acid alteration (G272E) in HmgA of Bt BMB181 was caused by a single point mutation (815G→ A). The enzyme activity of purified HmgA171 was more than 10-fold higher than that of HmgA181. Finally, we determined that the mutation in hmgA was responsible for melanin accumulation in Bt BMB181. Our results provided new insights into the synthesis and regulation of melanin production in B.thuringiensis and will promote its future industrial application.

Genome-wide association studies (GWAS) for kidney function identified hundreds of risk regions; however, the causal variants, target genes, cell types, and disease mechanisms remain poorly understood. Here, we performed transcriptome-wide association studies (TWAS), summary Mendelian randomization, and MetaXcan to identify genes whose expression mediates the genotype effect on the phenotype. Our analyses identified Dachshund homolog 1 (DACH1), a cell-fate determination factor. GWAS risk variant was associated with lower DACH1 expression in human kidney tubules. Human and mouse kidney single-cell open chromatin data (snATAC-Seq) prioritized estimated glomerular filtration rate (eGFR) GWAS variants located on an intronic regulatory region in distal convoluted tubule cells. CRISPR-Cas9-mediated gene editing confirmed the role of risk variants in regulating DACH1 expression. Mice with tubule-specific Dach1 deletion developed more severe renal fibrosis both in folic acid and diabetic kidney injury models. Mice with tubule-specific Dach1 overexpression were protected from folic acid nephropathy. Single-cell RNA sequencing, chromatin immunoprecipitation, and functional analysis indicated that DACH1 controls the expression of cell cycle and myeloid chemotactic factors, contributing to macrophage infiltration and fibrosis development. In summary, integration of GWAS, TWAS, single-cell epigenome, expression analyses, gene editing, and functional validation in different mouse kidney disease models identified DACH1 as a kidney disease risk gene.

Recent genetic data can offer important insights into the roles of lipoprotein subfractions and particle sizes in preventing coronary artery disease (CAD), as previous observational studies have often reported conflicting results. In this study, we first used the LD score regression to estimate the genetic correlation of 77 subfraction traits with traditional lipid profile and identified 27 traits that may involve distinct genetic mechanisms. We then used Mendelian randomization (MR) to estimate the causal effect of these traits on the risk of CAD. In univariable MR, the concentration and content of medium high-density lipoprotein (HDL) particles showed a protective effect against coronary artery disease. The effect was not attenuated in multivariable MR that adjusted for traditional lipid profile. The multivariable MR analyses also found that small HDL particles and smaller mean HDL particle diameter may have a protective effect. We identified four genetic markers for HDL particle size and CAD.

Determining the epigenetic program that generates unique cell types in the kidney is critical for understanding cell-type heterogeneity during tissue homeostasis and injury response. Here, we profile open chromatin and gene expression in developing and adult mouse kidneys at single cell resolution. We show critical reliance of gene expression on distal regulatory elements (enhancers). We reveal key cell type-specific transcription factors and major gene-regulatory circuits for kidney cells. Dynamic chromatin and expression changes during nephron progenitor differentiation demonstrates that podocyte commitment occurs early and is associated with sustained Foxl1 expression. Renal tubule cells follow a more complex differentiation, where Hfn4a is associated with proximal and Tfap2b with distal fate. Mapping single nucleotide variants associated with human kidney disease implicates critical cell types, developmental stages, genes, and regulatory mechanisms. The single cell multi-omics atlas reveals key chromatin remodeling events and gene expression dynamics associated with kidney development., Epigenetic and transcriptional dynamics are critical for both tissue homeostasis and injury response in the kidney. Leveraging a single cell multiomics atlas of the developing mouse kidney, the authors reveal key events in chromatin regulation and gene expression dynamics during postnatal development.

Mapping single-cell RNA-seq data onto IHC images enables the annotation of cell populations at high phenotypic resolution., Highly multiplexed immunohistochemistry (mIHC) enables the staining and quantification of dozens of antigens in a tissue section with single-cell resolution. However, annotating cell populations that differ little in the profiled antigens or for which the antibody panel does not include specific markers is challenging. To overcome this obstacle, we have developed an approach for enriching mIHC images with single-cell RNA sequencing data, building upon recent experimental procedures for augmenting single-cell transcriptomes with concurrent antigen measurements. Spatially-resolved Transcriptomics via Epitope Anchoring (STvEA) performs transcriptome-guided annotation of highly multiplexed cytometry datasets. It increases the level of detail in histological analyses by enabling the systematic annotation of nuanced cell populations, spatial patterns of transcription, and interactions between cell types. We demonstrate the utility of STvEA by uncovering the architecture of poorly characterized cell types in the murine spleen using published cytometry and mIHC data of this organ.

The functional interpretation of GWAS remains challenging due to cell-type dependent influences of genetic variants.Here, we generated comprehensive maps of expression quantitative trait loci (eQTL) for 659 microdissected human kidney samples and identified cell-type eQTLs by mapping interactions between cell type abundance and genotype. Separately, we generated single cell open chromatin maps (by snATAC-seq) for human kidney samples. We highlight critical enrichment of proximal tubules in kidney function and endothelial cells and distal tubule segments in blood pressure by partitioning heritability using stratified LD-score regression to integrate GWAS with scRNA-seq and snATAC-seq data. Bayesian colocalization analysis nominated more than 200 genes for kidney function and hypertension. Our study clarifies the mechanism of the most commonly used antihypertensive and renal protective drugs and identifies drug repurposing opportunities for kidney disease.One Sentence SummaryWe define causal cell types, genes and mechanism for kidney dysfunction.

Beta-amyloid deposition is a defining feature of Alzheimer's disease (AD). How genetic risk factors, like APOE and TREM2, intersect with cellular responses to beta-amyloid in human tissues is not fully understood. Using single-nucleus RNA sequencing of postmortem human brain with varied APOE and TREM2 genotypes and neuropathology, we identified distinct microglia subpopulations, including a subpopulation of CD163-positive amyloid-responsive microglia (ARM) that are depleted in cases with APOE and TREM2 risk variants. We validated our single-nucleus RNA sequencing findings in an expanded cohort of AD cases, demonstrating that APOE and TREM2 risk variants are associated with a significant reduction in CD163-positive amyloid-responsive microglia. Our results showcase the diverse microglial response in AD and underscore how genetic risk factors influence cellular responses to underlying pathologies.

Hypoxia is a critical, but frequently overlooked problem, which commonly exists in Chinese mitten crab rearing. However, little information is available on the molecular mechanisms of the detrimental effects of hypoxia in this species. In the present study, crabs were subjected to acute hypoxia stress (DO 1.0 mg/L), followed by reoxygenation (DO 6.8 mg/L). Hepatopancreas from five groups of crabs (three to four crabs per group), including normoxia, hypoxia for one and six hours, and reoxygenation for one and 12 h, were used for transcriptome sequencing. The pooled total RNA of all samples were utilized to reconstruct a reference transcriptome with PacBio RS II sequencing, obtaining 49.19 G clean data, with a mean length of 2,180 bp. Seventeen cDNA libraries were constructed and sequenced to identify differentially expressed genes (DEGs) among the different samples (FDR 0.05 and |log2 fold change| ≥1). A total of 103 and 251 DEGs were identified when exposed to hypoxia for one and six hours, respectively. Totally 462 and 673 DEGs were identified during reoxygenation at 1 and 12 h, respectively. Among these DEGs, two transcripts with complete ORFs were identified to encode hypoxia-inducible factor 1 (Es-Hif-1α/β), which is a transcriptional activator of various genes correlated to the cellular adaptive responses to hypoxia. Es-Hif-1a/β expressions were significantly upregulated when exposed to hypoxia treatment, and no significant difference was observed for Es-Hif-1α between hypoxic treatment for 6 h and reoxygenation. The significant KEGG enrichment revealed that the DEGs under hypoxia were mainly enriched in "PPAR signaling pathway", "Gap junction" and "Phototransduction-fly". Compared with crabs in normoxia, even with 12 h of reoxygenation, the hepatopancreas transcriptome under hypoxia remained severely affected, implying its adverse effect on growth and development, or increased susceptibility to diseases.

Determining the epigenetic program that generates unique cell types in the kidney is critical for understanding cell-type heterogeneity during tissue homeostasis and injury response.Here, we profiled open chromatin and gene expression in developing and adult mouse kidneys at single cell resolution. We show critical reliance of gene expression on distal regulatory elements (enhancers). We define key cell type-specific transcription factors and major gene-regulatory circuits for kidney cells. Dynamic chromatin and expression changes during nephron progenitor differentiation demonstrated that podocyte commitment occurs early and is associated with sustainedFoxl1expression. Renal tubule cells followed a more complex differentiation, whereHfn4awas associated with proximal andTfap2bwith distal fate. Mapping single nucleotide variants associated with human kidney disease identified critical cell types, developmental stages, genes, and regulatory mechanisms.We provide a global single cell resolution view of chromatin accessibility of kidney development. The dataset is available via interactive public websites.

Bone marrow mesenchymal lineage cells are a heterogeneous cell population involved in bone homeostasis and diseases such as osteoporosis. While it is long postulated that they originate from mesenchymal stem cells, the true identity of progenitors and their in vivo bifurcated differentiation routes into osteoblasts and adipocytes remain poorly understood. Here, by employing large scale single cell transcriptome analysis, we computationally defined mesenchymal progenitors at different stages and delineated their bi-lineage differentiation paths in young, adult and aging mice. One identified subpopulation is a unique cell type that expresses adipocyte markers but contains no lipid droplets. As non-proliferative precursors for adipocytes, they exist abundantly as pericytes and stromal cells that form a ubiquitous 3D network inside the marrow cavity. Functionally they play critical roles in maintaining marrow vasculature and suppressing bone formation. Therefore, we name them marrow adipogenic lineage precursors (MALPs) and conclude that they are a newly identified component of marrow adipose tissue.

Age-related macular degeneration (AMD) preferentially affects distinct cell types and topographic regions in retina. To characterize the impact of AMD on gene expression changes across retinal cell types and regions, we generated both single-cell RNA-seq (scRNA-seq) and bulk RNA-seq data from macular and peripheral retina in postmortem human donors with and without AMD. The scRNA-seq data revealed 11 major cell types with many previously reported AMD risk genes showing substantial cell type and region specificity. Cell type proportional changes with advancing AMD stage were significant for Müller glia, rods, astrocytes, microglia and endothelium.

Bone marrow mesenchymal lineage cells are a heterogeneous cell population involved in bone homeostasis and diseases such as osteoporosis. While it is long postulated that they originate from mesenchymal stem cells (MSCs), the true identity of MSCs and their in vivo bifurcated differentiation routes into osteoblasts and adipocytes remain poorly understood. Here, by employing single cell transcriptome analysis, we identified MSCs and delineated their bi-lineage differentiation paths in young, adult and aging mice. Among several newly discovered mesenchymal subpopulations, one is a distinct population of adipose-lineage cells that we named marrow environment regulating adipose cells (MERAs). MERAs are non-proliferative, post-progenitor cells that express many mature adipocyte markers but are devoid of lipid droplets. They are abundant in the bone marrow of young mice, acting as pericytes and stromal cells that form numerous connections among themselves and with other cells inside bone, including endothelial cells. Genetic ablation of MERAs disrupts marrow vessel structure, promotes de novo bone formation. Taken together, MERAs represent a unique population of adipose lineage cells that exist only in the bone marrow with critical roles in regulating bone and vessel homeostasis.

Histology provides a unique window into the cellular and molecular architecture of tissues and is a critical component of biomedical research and clinical practice. Highly-multiplexed immunohistochemistry 1-6 (mIHC) enables the routine staining and quantification of dozens of antigens in the same tissue section with single-cell resolution. However, the amount of cell types and states that can be simultaneously identified by mIHC is limited. In contrast, cells are finely disaggregated into distinct types in single-cell transcriptomic analyses but spatial information is lost. To bridge this gap, we developed an approach for enriching mIHC histology slides with single-cell RNA-seq data, building upon recent experimental procedures for augmenting single-cell transcriptomes with concurrent antigen measurements 7, 8 . Our approach, Spatially-resolved Transcriptomics via Epitope Anchoring (STvEA), increases the level of detail in histological analyses by enabling detection of subtle cell populations, spatial patterns of transcription, and cell-to-cell interactions. It provides an improvement in throughput, resolution, and simplicity with respect to existing spatially-resolved methods for simultaneous proteomics and transcriptomics. We demonstrate the utility of STvEA by uncovering the architecture of poorly characterized cell populations in the murine spleen using published mIHC images.

Recently developed technologies for digital imaging and highly multiplexed immunohistochemistry (mIHC) are enabling complex quantitative descriptions of tissue architecture. Imaging mass cytometry, multiplexed ion beam imaging, and co-detection by indexing can be used to profile the expression of dozens of proteins in a tissue section with single-cell resolution. Despite this progress, the amount of cell types and states that can be simultaneously identified by mIHC in a tissue section is still limited. Standard methods for identifying cell populations cluster cells according to expression similarities of the profiled antigens. These clusters are then manually annotated using previous knowledge of cell population markers. However, this is generally a partial, subjective, and biased process. Clusters often differ little in their antigenic profile and the interpretation of those differences is unclear. Moreover, the design of comprehensive antibody panels that include specific markers for every cell type and state present in the tissue is usually unfeasible. To overcome these limitations and increase the level of annotation of mIHC data, we propose an approach for enriching mIHC slides with single-cell RNA-seq data. Currently available single-cell RNA-seq technologies profile the expression of thousands of genes in each cell, allowing researchers to finely classify cells based on their gene expression profile. Some of the most recent approaches, like CITE-seq, REAP-seq, and Ab-seq, allow augmenting single-cell transcriptomes with concurrent protein measurements by staining dissociated cells with oligo-tagged antibodies. These approaches can be therefore used to determine the quantitative relation between gene and antigen expression levels. Here, we build upon CITE-seq and computational methods for the integration of single-cell omics data to identify and annotate cell populations in mIHC images in an automated way based on single-cell gene expression data from the same tissue. Our method, Spatially-resolved Transcriptomics via Epitope Anchoring (STvEA), performs transcriptome-guided annotation of mIHC datasets. It increases the phenotypic resolution of histological analyses by enabling annotation of subtle cell populations, spatial patterns of transcription, and interactions between cell types. We use the murine spleen as a test system to benchmark the stability and performance of STvEA and demonstrate its utility by using it to annotate published CODEX and CyTOF datasets of this organ. Citation Format: Kiya W. Govek, Emma C. Troisi, Zhen Miao, Steven Woodhouse, Pablo G. Camara. Single-cell transcriptomic analysis of highly-multiplexed IHC images via antigen mapping [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PO-078.

Outbreaks of mass mortalities among cultured Macrobrachium nipponense occurred in a commercial hatchery during the autumn of 2017 in Jiangsu province, P. R. China, and non-O1 Vibrio cholerae was isolated and identified as causal agents of M. nipponense, with a LD50 value 4.09 × 104 CFU/mL. Detection of virulence-associated genes by PCR indicated that XL1 was positive for Mp, HlyA, RtxA, OmpU, Ace, Zot and T6SS. Furthermore, the results of extracellular enzyme analysis revealed that the strain can produce lecithinase, amylase, gelatinase and hemolysin. Histopathological analysis revealed that the hepatic tubule lumen and the gap between the hepatic tubules became larger, and the brush border disappeared in the hepatopancreas. Quantitive real-time PCR (qRT-PCR) was undertaken to measure mRNA expression levels for thirteen immune related genes in M. nipponense after non-O1 V. cholerae infection. The transcriptional analysis of these immune related genes demonstrated that the expression levels of dorsal, relish, p38, crustin1, crustin2, crustin3, hemocyanin, i-lysozyme, anti-lipopolysaccharide factors 1, anti-lipopolysaccharide factors 2, prophenoloxidase were significantly up-regulated in hemolymph of M. nipponense post-infection. These results revealed varying expression profiles and clear transcriptional activation of these immune related genes in hemolymph, which will contribute to better understand the pathogenesis and host defensive system in non-O1 V. cholerae invasion.

Plant secondary metabolites are known to not only play a key role in the adaptation of plants to their environment, but also represent an important source of active pharmaceuticals. Alpinia oxyphylla capsular fruits, made up of seeds and pericarps, are commonly used in traditional East Asian medicines. In clinical utilization of these capsular fruits, inconsistent processing approaches (i.e., hulling pericarps or not) are employed, with the potential of leading to differential pharmacological effects. Therefore, an important question arises whether the content levels of pharmacologically active chemicals between the seeds and pericarps of A. oxyphylla are comparable. Nine secondary metabolites present in A. oxyphylla capsular fruits, including flavonoids (e.g., tectochrysin, izalpinin, chrysin, apigenin-4',7-dimethylether and kaempferide), diarylheptanoids (e.g., yakuchinone A and B and oxyphyllacinol) and sesquiterpenes (e.g., nootkatone), were regarded as representative constituents with putative pharmacological activities. This work aimed to investigate the abundance of the nine constituents in the seeds and pericarps of A. oxyphylla. Thirteen batches of A. oxyphylla capsular fruits were gathered from different production regions. Accordingly, an ultra-fast high performance liquid chromatography/quadrupole tandem mass spectrometry (UFLC-MS/MS) method was developed and validated. We found that: (1) the nine secondary metabolites were differentially concentrated in seeds and fruit capsules; (2) nootkatone is predominantly distributed in the seeds; in contrast, the flavonoids and diarylheptanoids are mainly deposited in the capsules; and (3) the content levels of the nine secondary metabolites occurring in the capsules varied greatly among different production regions, although the nootkatone levels in the seeds were comparable among production regions. These results are helpful to evaluating and elucidating pharmacological activities of A. oxyphylla capsular fruits. Additionally, it may be of interest to elucidate the mechanisms involved in the distinct accumulation profiles of these secondary metabolites between seeds and pericarps.