Perturbed liver gene zonation in a mouse model of non-alcoholic steatohepatitis.

Liver zonation characterizes the separation of metabolic pathways along the lobules and is required for optimal hepatic function. Wnt signaling is a master regulator of spatial liver zonation. A perivenous-periportal Wnt activity gradient orchestrates metabolic zonation by activating gene expression in perivenous hepatocytes, while suppressing gene expression in their periportal counterparts. However, the understanding as to the liver gene zonation and zonation regulators in diseases is limited. Non-alcoholic steatohepatitis (NASH) is a chronic liver disease characterized by fat accumulation, inflammation, and fibrosis. Here, we investigated the perturbation of liver gene zonation in a mouse NASH model by combining spatial transcriptomics, bulk RNAseq and in situ hybridization. Wnt-target genes represented a major subset of genes showing altered spatial expression in the NASH liver. The altered Wnt-target gene expression levels and zonation spatial patterns were in line with the up regulation of Wnt regulators and the augmentation of Wnt signaling. Particularly, we found that the Wnt activator Rspo3 expression was restricted to the perivenous zone in control liver but expanded to the periportal zone in NASH liver. AAV8-mediated RSPO3 overexpression in controls resulted in zonation changes, and further amplified the disturbed zonation of Wnt-target genes in NASH, similarly Rspo3 knockdown in Rspo3+/- mice resulted in zonation changes of Wnt-target genes in both chow and HFD mouse. Interestingly, there were no impacts on steatosis, inflammation, or fibrosis NASH pathology from RSPO3 overexpression nor Rspo3 knockdown. In summary, our study demonstrated the alteration of Wnt signaling in a mouse NASH model, leading to perturbed liver zonation.

Multiancestry exome sequencing reveals INHBE mutations associated with favorable fat distribution and protection from diabetes.

Body fat distribution is a major, heritable risk factor for cardiometabolic disease, independent of overall adiposity. Using exome-sequencing in 618,375 individuals (including 160,058 non-Europeans) from the UK, Sweden and Mexico, we identify 16 genes associated with fat distribution at exome-wide significance. We show 6-fold larger effect for fat-distribution associated rare coding variants compared with fine-mapped common alleles, enrichment for genes expressed in adipose tissue and causal genes for partial lipodystrophies, and evidence of sex-dimorphism. We describe an association with favorable fat distribution (p = 1.8 × 10-09), favorable metabolic profile and protection from type 2 diabetes (~28% lower odds; p = 0.004) for heterozygous protein-truncating mutations in INHBE, which encodes a circulating growth factor of the activin family, highly and specifically expressed in hepatocytes. Our results suggest that inhibin ßE is a liver-expressed negative regulator of adipose storage whose blockade may be beneficial in fat distribution-associated metabolic disease.

Germline Mutations in CIDEB and Protection against Liver Disease.

BACKGROUND: Exome sequencing in hundreds of thousands of persons may enable the identification of rare protein-coding genetic variants associated with protection from human diseases like liver cirrhosis, providing a strategy for the discovery of new therapeutic targets. METHODS: We performed a multistage exome sequencing and genetic association analysis to identify genes in which rare protein-coding variants were associated with liver phenotypes. We conducted in vitro experiments to further characterize associations. RESULTS: The multistage analysis involved 542,904 persons with available data on liver aminotransferase levels, 24,944 patients with various types of liver disease, and 490,636 controls without liver disease. We found that rare coding variants in APOB, ABCB4, SLC30A10, and TM6SF2 were associated with increased aminotransferase levels and an increased risk of liver disease. We also found that variants in CIDEB, which encodes a structural protein found in hepatic lipid droplets, had a protective effect. The burden of rare predicted loss-of-function variants plus missense variants in CIDEB (combined carrier frequency, 0.7%) was associated with decreased alanine aminotransferase levels (beta per allele, -1.24 U per liter; 95% confidence interval [CI], -1.66 to -0.83; P = 4.8×10-9) and with 33% lower odds of liver disease of any cause (odds ratio per allele, 0.67; 95% CI, 0.57 to 0.79; P = 9.9×10-7). Rare coding variants in CIDEB were associated with a decreased risk of liver disease across different underlying causes and different degrees of severity, including cirrhosis of any cause (odds ratio per allele, 0.50; 95% CI, 0.36 to 0.70). Among 3599 patients who had undergone bariatric surgery, rare coding variants in CIDEB were associated with a decreased nonalcoholic fatty liver disease activity score (beta per allele in score units, -0.98; 95% CI, -1.54 to -0.41 [scores range from 0 to 8, with higher scores indicating more severe disease]). In human hepatoma cell lines challenged with oleate, CIDEB small interfering RNA knockdown prevented the buildup of large lipid droplets. CONCLUSIONS: Rare germline mutations in CIDEB conferred substantial protection from liver disease. (Funded by Regeneron Pharmaceuticals.).

Cardiac-gated parametric images from 82 Rb PET from dynamic frames and direct 4D reconstruction.

PURPOSE: Cardiac perfusion PET data can be reconstructed as a dynamic sequence and kinetic modeling performed to quantify myocardial blood flow, or reconstructed as static gated images to quantify function. Parametric images from dynamic PET are conventionally not gated, to allow use of all events with lower noise. An alternative method for dynamic PET is to incorporate the kinetic model into the reconstruction algorithm itself, bypassing the generation of a time series of emission images and directly producing parametric images. So-called "direct reconstruction" can produce parametric images with lower noise than the conventional method because the noise distribution is more easily modeled in projection space than in image space. In this work, we develop direct reconstruction of cardiac-gated parametric images for 82 Rb PET with an extension of the Parametric Motion compensation OSEM List mode Algorithm for Resolution-recovery reconstruction for the one tissue model (PMOLAR-1T). METHODS: PMOLAR-1T was extended to accommodate model terms to account for spillover from the left and right ventricles into the myocardium. The algorithm was evaluated on a 4D simulated 82 Rb dataset, including a perfusion defect, as well as a human 82 Rb list mode acquisition. The simulated list mode was subsampled into replicates, each with counts comparable to one gate of a gated acquisition. Parametric images were produced by the indirect (separate reconstructions and modeling) and direct methods for each of eight low-count and eight normal-count replicates of the simulated data, and each of eight cardiac gates for the human data. For the direct method, two initialization schemes were tested: uniform initialization, and initialization with the filtered iteration 1 result of the indirect method. For the human dataset, event-by-event respiratory motion compensation was included. The indirect and direct methods were compared for the simulated dataset in terms of bias and coefficient of variation as a function of iteration. RESULTS: Convergence of direct reconstruction was slow with uniform initialization; lower bias was achieved in fewer iterations by initializing with the filtered indirect iteration 1 images. For most parameters and regions evaluated, the direct method achieved the same or lower absolute bias at matched iteration as the indirect method, with 23%-65% lower noise. Additionally, the direct method gave better contrast between the perfusion defect and surrounding normal tissue than the indirect method. Gated parametric images from the human dataset had comparable relative performance of indirect and direct, in terms of mean parameter values per iteration. Changes in myocardial wall thickness and blood pool size across gates were readily visible in the gated parametric images, with higher contrast between myocardium and left ventricle blood pool in parametric images than gated SUV images. CONCLUSIONS: Direct reconstruction can produce parametric images with less noise than the indirect method, opening the potential utility of gated parametric imaging for perfusion PET.

Non-Rigid Event-by-Event Continuous Respiratory Motion Compensated List-Mode Reconstruction for PET.

Respiratory motion during positron emission tomography (PET)/computed tomography (CT) imaging can cause significant image blurring and underestimation of tracer concentration for both static and dynamic studies. In this paper, with the aim to eliminate both intra-cycle and inter-cycle motions, and apply to dynamic imaging, we developed a non-rigid event-by-event (NR-EBE) respiratory motion-compensated list-mode reconstruction algorithm. The proposed method consists of two components: the first component estimates a continuous non-rigid motion field of the internal organs using the internal-external motion correlation. This continuous motion field is then incorporated into the second component, non-rigid MOLAR (NR-MOLAR) reconstruction algorithm to deform the system matrix to the reference location where the attenuation CT is acquired. The point spread function (PSF) and time-of-flight (TOF) kernels in NR-MOLAR are incorporated in the system matrix calculation, and therefore are also deformed according to motion. We first validated NR-MOLAR using a XCAT phantom with a simulated respiratory motion. NR-EBE motion-compensated image reconstruction using both the components was then validated on three human studies injected with 18F-FPDTBZ and one with 18F-fluorodeoxyglucose (FDG) tracers. The human results were compared with conventional non-rigid motion correction using discrete motion field (NR-discrete, one motion field per gate) and a previously proposed rigid EBE motion-compensated image reconstruction (R-EBE) that was designed to correct for rigid motion on a target lesion/organ. The XCAT results demonstrated that NR-MOLAR incorporating both PSF and TOF kernels effectively corrected for non-rigid motion. The 18F-FPDTBZ studies showed that NR-EBE out-performed NR-Discrete, and yielded comparable results with R-EBE on target organs while yielding superior image quality in other regions. The FDG study showed that NR-EBE clearly improved the visibility of multiple moving lesions in the liver where some of them could not be discerned in other reconstructions, in addition to improving quantification. These results show that NR-EBE motion-compensated image reconstruction appears to be a promising tool for lesion detection and quantification when imaging thoracic and abdominal regions using PET.

Direct reconstruction of parametric images for brain PET with event-by-event motion correction: evaluation in two tracers across count levels.

Parametric images for dynamic positron emission tomography (PET) are typically generated by an indirect method, i.e. reconstructing a time series of emission images, then fitting a kinetic model to each voxel time activity curve. Alternatively, 'direct reconstruction', incorporates the kinetic model into the reconstruction algorithm itself, directly producing parametric images from projection data. Direct reconstruction has been shown to achieve parametric images with lower standard error than the indirect method. Here, we present direct reconstruction for brain PET using event-by-event motion correction of list-mode data, applied to two tracers. Event-by-event motion correction was implemented for direct reconstruction in the Parametric Motion-compensation OSEM List-mode Algorithm for Resolution-recovery reconstruction. The direct implementation was tested on simulated and human datasets with tracers [11C]AFM (serotonin transporter) and [11C]UCB-J (synaptic density), which follow the 1-tissue compartment model. Rigid head motion was tracked with the Vicra system. Parametric images of K 1 and distribution volume (V T = K 1/k 2) were compared to those generated by the indirect method by regional coefficient of variation (CoV). Performance across count levels was assessed using sub-sampled datasets. For simulated and real datasets at high counts, the two methods estimated K 1 and V T with comparable accuracy. At lower count levels, the direct method was substantially more robust to outliers than the indirect method. Compared to the indirect method, direct reconstruction reduced regional K 1 CoV by 35-48% (simulated dataset), 39-43% ([11C]AFM dataset) and 30-36% ([11C]UCB-J dataset) across count levels (averaged over regions at matched iteration); V T CoV was reduced by 51-58%, 54-60% and 30-46%, respectively. Motion correction played an important role in the dataset with larger motion: correction increased regional V T by 51% on average in the [11C]UCB-J dataset. Direct reconstruction of dynamic brain PET with event-by-event motion correction is achievable and dramatically more robust to noise in V T images than the indirect method.

Quantification of myocardial blood flow with (82)Rb: Validation with (15)O-water using time-of-flight and point-spread-function modeling.

BACKGROUND: We quantified myocardial blood flow with (82)Rb PET using parameters of the generalized Renkin-Crone model estimated from (82)Rb and (15)O-water images reconstructed with time-of-flight and point spread function modeling. Previous estimates of rubidium extraction have used older-generation scanners without time-of-flight or point spread function modeling. We validated image-derived input functions with continuously collected arterial samples. METHODS: Nine healthy subjects were scanned at rest and under pharmacological stress on the Siemens Biograph mCT with (82)Rb and (15)O-water PET, undergoing arterial blood sampling with each scan. Image-derived input functions were estimated from the left ventricle cavity and corrected with tracer-specific population-based scale factors determined from arterial data. Kinetic parametric images were generated from the dynamic PET images by fitting the one-tissue compartment model to each voxel's time activity curve. Mean myocardial blood flow was determined from each subject's (15)O-water k 2 images. The parameters of the generalized Renkin-Crone model were estimated from these water-based flows and mean myocardial (82)Rb K 1 estimates. RESULTS: Image-derived input functions showed improved agreement with arterial measurements after a scale correction. The Renkin-Crone model fit (a = 0.77, b = 0.39) was similar to those previously published, though b was lower. CONCLUSIONS: We have presented parameter estimates for the generalized Renkin-Crone model of extraction for (82)Rb PET using human (82)Rb and (15)O-water PET from high-resolution images using a state-of-the-art time-of-flight-capable scanner. These results provide a state-of-the-art methodology for myocardial blood flow measurement with (82)Rb PET.

Event-by-Event Continuous Respiratory Motion Correction for Dynamic PET Imaging.

UNLABELLED: Existing respiratory motion-correction methods are applied only to static PET imaging. We have previously developed an event-by-event respiratory motion-correction method with correlations between internal organ motion and external respiratory signals (INTEX). This method is uniquely appropriate for dynamic imaging because it corrects motion for each time point. In this study, we applied INTEX to human dynamic PET studies with various tracers and investigated the impact on kinetic parameter estimation. METHODS: The use of 3 tracers-a myocardial perfusion tracer, (82)Rb (n = 7); a pancreatic ß-cell tracer, (18)F-FP(+)DTBZ (n = 4); and a tumor hypoxia tracer, (18)F-fluoromisonidazole ((18)F-FMISO) (n = 1)-was investigated in a study of 12 human subjects. Both rest and stress studies were performed for (82)Rb. The Anzai belt system was used to record respiratory motion. Three-dimensional internal organ motion in high temporal resolution was calculated by INTEX to guide event-by-event respiratory motion correction of target organs in each dynamic frame. Time-activity curves of regions of interest drawn based on end-expiration PET images were obtained. For (82)Rb studies, K1 was obtained with a 1-tissue model using a left-ventricle input function. Rest-stress myocardial blood flow (MBF) and coronary flow reserve (CFR) were determined. For (18)F-FP(+)DTBZ studies, the total volume of distribution was estimated with arterial input functions using the multilinear analysis 1 method. For the (18)F-FMISO study, the net uptake rate Ki was obtained with a 2-tissue irreversible model using a left-ventricle input function. All parameters were compared with the values derived without motion correction. RESULTS: With INTEX, K1 and MBF increased by 10% ± 12% and 15% ± 19%, respectively, for (82)Rb stress studies. CFR increased by 19% ± 21%. For studies with motion amplitudes greater than 8 mm (n = 3), K1, MBF, and CFR increased by 20% ± 12%, 30% ± 20%, and 34% ± 23%, respectively. For (82)Rb rest studies, INTEX had minimal effect on parameter estimation. The total volume of distribution of (18)F-FP(+)DTBZ and Ki of (18)F-FMISO increased by 17% ± 6% and 20%, respectively. CONCLUSION: Respiratory motion can have a substantial impact on dynamic PET in the thorax and abdomen. The INTEX method using continuous external motion data substantially changed parameters in kinetic modeling. More accurate estimation is expected with INTEX.

Regional, kinetic [(18)F]FDG PET imaging of a unilateral Parkinsonian animal model.

Positron emission tomography (PET) imaging with the glucose analog 2-deoxy-2-[(18)F]fluoro-D-glucose ([(18)F] FDG) has demonstrated clinical utility for the monitoring of brain glucose metabolism alteration in progressive neurodegenerative diseases. We examined dynamic [(18)F]FDG PET imaging and kinetic modeling of atlas-based regions to evaluate regional changes in the cerebral metabolic rate of glucose in the widely-used 6-hydroxydopamine (6-OHDA) rat model of Parkinson's disease. Following a bolus injection of 18.5 ± 1 MBq [(18)F]FDG and a 60-minute PET scan, image-derived input functions from the vena cava and left ventricle were used with three models, including Patlak graphical analysis, to estimate the influx constant and the metabolic rate in ten brain regions. We observed statistically significant changes in [(18)F]FDG uptake ipsilateral to the 6-OHDA injection in the basal ganglia, olfactory bulb, and amygdala regions; and these changes are of biological relevance to the disease. These experiments provide further validation for the use of [(18)F]FDG PET imaging in this model for drug discovery and development.

A vector for double epitope tagging with a recyclable marker.

Multimeric protein complexes play diverse and vital roles in the cell, but following the composition of these complexes under varying growth conditions can be challenging. Toward that goal, we have designed a vector that permits the double epitope tagging of a protein at its carboxy terminus. One 'universal' tag, a triple repeat of the HA1 epitope, is fused with every protein to be studied, allowing the composition and stoichiometry of the proteins in a complex to be detected with a single antibody. Each protein also can be tagged with a second epitope specific for that protein. This 'specific' tag can be used to immunoprecipitate complexes containing that protein of interest. Any epitope to which a specific antibody is available can be used for this second tag. Because there are a limited number of selection markers for cloning in yeast, the kanamycin cassette, flanked by loxP sites, was incorporated into the vector to permit marker recycling using Cre-lox recombinase. This vector was used to tag 4 proteins involved in ribosome biogenesis-Ytm1, Cic1, Brx1 and Drs1. An anti-HA1 antibody could detect all four proteins in crude lysates and yielded the relative abundance of these four proteins, of which Drs1 is reproducibly less abundant than any of the others, which may have implications for the control of ribosome biogenesis. The Ytm1 protein was also tagged with the VSV epitope and can be specifically detected using an anti-VSV antibody. This vector may prove useful for exploring other protein complexes.