Genetic evidence that the human CYP2R1 enzyme is a key vitamin D 25-hydroxylase

The synthesis of bioactive vitamin D requires hydroxylation at the 1 alpha and 25 positions by cytochrome P450 enzymes in the kidney and liver, respectively. The mitochondrial enzyme CYP27B1 catalyzes 1 alpha-hydroxylation in the kidney but the identity of the hepatic 25-hydroxylase has remained unclear for >30 years. We previously identified the microsomal CYP2R1 protein as a potential candidate for the liver vitamin D 25-hydroxylase based on the enzyme's biochemical properties, conservation, and expression pattern. Here, we report a molecular analysis of a patient with low circulating levels of 25-hydroxyvitamin D and classic symptoms of vitamin D deficiency. This individual was found to be homozygous for a transition mutation in exon 2 of the CYP2R1 gene on chromosome 11p15.2. The inherited mutation caused the substitution of a proline for an evolutionarily conserved leucine at amino acid 99 in the CYP2R1 protein and eliminated vitamin D 25-hydroxylase enzyme activity. These data identify CYP2R1 as a biologically relevant vitamin D 25-hydroxylase and reveal the molecular basis of a human genetic disease, selective 25-hydroxyvitamin D deficiency.

De-orphanization of cytochrome P450 2R1: a microsomal vitamin D 25-hydroxilase

The conversion of vitamin D into an active ligand for the vitamin D receptor requires 25-hydroxylation in the liver and 1alpha-hydroxylation in the kidney. Mitochondrial and microsomal vitamin D 25-hydroxylase enzymes catalyze the first reaction. The mitochondrial activity is associated with sterol 27-hydroxylase, a cytochrome P450 (CYP27A1); however, the identity of the microsomal enzyme has remained elusive. A cDNA library prepared from hepatic mRNA of sterol 27-hydroxylase-deficient mice was screened with a ligand activation assay to identify an evolutionarily conserved microsomal cytochrome P450 (CYP2R1) with vitamin D 25-hydroxylase activity. Expression of CYP2R1 in cells led to the transcriptional activation of the vitamin D receptor when either vitamin D2 or D3 was added to the medium. Thin layer chromatography and radioimmunoassays indicated that the secosteroid product of CYP2R1 was 25-hydroxyvitamin D3. Co-expression of CYP2R1 with vitamin D 1alpha-hydroxylase (CYP27B1) elicited additive activation of vitamin D3, whereas co-expression with vitamin D 24-hydroxylase (CYP24A1) caused inactivation. CYP2R1 mRNA is abundant in the liver and testis, and present at lower levels in other tissues. The data suggest that CYP2R1 is a strong candidate for the microsomal vitamin D 25-hydroxylase.

Isolation and characterization of PDE9A, a novel human cGMP-specific phosphodiesterase

We have cloned and characterized the first human isozyme in a new family of cyclic nucleotide phosphodiesterases, PDE9A. By sequence homology in the catalytic domain, PDE9A is almost equidistant from all eight known mammalian PDE families but is most similar to PDE8A (34% amino acid identity) and least like PDE5A (28% amino acid identity). We report the cloning of human cDNA encoding a full-length protein of 593 amino acids, including a 261-amino acid region located near the C terminus that is homologous to the approximately 270-amino acid catalytic domain of other PDEs. PDE9A is expressed in all eight tissues examined as a approximately 2. 0-kilobase mRNA, with highest levels in spleen, small intestine, and brain. The full-length PDE9A was expressed in baculovirus fused to an N-terminal 9-amino acid FLAG tag. Kinetic analysis of the baculovirus-expressed enzyme shows it to be a very high affinity cGMP-specific PDE with a Km of 170 nM for cGMP and 230 microM for cAMP. The Km for cGMP makes PDE9A one of the highest affinity PDEs known. The Vmax for cGMP (4.9 nmol/min/microg recombinant enzyme) is about twice as fast as that of PDE4 for cAMP. The enzyme is about twice as active in vitro in 1-10 mM Mn2+ than in the same concentration of Mg2+ or Ca2+. PDE9A is insensitive (up to 100 microM) to a variety of PDE inhibitors including rolipram, vinpocetine, SKF-94120, dipyridamole, and 3-isobutyl-1-methyl-xanthine but is inhibited (IC50 = 35 microM) by zaprinast, a PDE5 inhibitor. PDE9A lacks a region homologous to the allosteric cGMP-binding regulatory regions found in the cGMP-binding PDEs: PDE2, PDE5, and PDE6.

Developmental fate and cellular maturity encoded in human regulatory DNA landscapes

Cellular-state information between generations of developing cells may be propagated via regulatory regions. We report consistent patterns of gain and loss of DNase I-hypersensitive sites (DHSs) as cells progress from embryonic stem cells (ESCs) to terminal fates. DHS patterns alone convey rich information about cell fate and lineage relationships distinct from information conveyed by gene expression. Developing cells share a proportion of their DHS landscapes with ESCs; that proportion decreases continuously in each cell type as differentiation progresses, providing a quantitative benchmark of developmental maturity. Developmentally stable DHSs densely encode binding sites for transcription factors involved in autoregulatory feedback circuits. In contrast to normal cells, cancer cells extensively reactivate silenced ESC DHSs and those from developmental programs external to the cell lineage from which the malignancy derives. Our results point to changes in regulatory DNA landscapes as quantitative indicators of cell-fate transitions, lineage relationships, and dysfunction.

Transcriptional Programming of Normal and Inflamed Human Epidermis at Single-Cell Resolution

Perturbations in the transcriptional programs specifying epidermal differentiation cause diverse skin pathologies ranging from impaired barrier function to inflammatory skin disease. However, the global scope and organization of this complex cellular program remain undefined. Here we report single-cell RNA sequencing profiles of 92,889 human epidermal cells from 9 normal and 3 inflamed skin samples. Transcriptomics-derived keratinocyte subpopulations reflect classic epidermal strata but also sharply compartmentalize epithelial functions such as cell-cell communication, inflammation, and WNT pathway modulation. In keratinocytes, ∼12% of assessed transcript expression varies in coordinate patterns, revealing undescribed gene expression programs governing epidermal homeostasis. We also identify molecular fingerprints of inflammatory skin states, including S100 activation in the interfollicular epidermis of normal scalp, enrichment of a CD1C+CD301A+ myeloid dendritic cell population in psoriatic epidermis, and IL1βhiCCL3hiCD14+ monocyte-derived macrophages enriched in foreskin. This compendium of RNA profiles provides a critical step toward elucidating epidermal diseases of development, differentiation, and inflammation.

Mammalian wax biosynthesis. I. Identification of two fatty acyl-Coenzyme A reductases with different substrate specificities and tissue distributions

The conversion of fatty acids to fatty alcohols is required for the synthesis of wax monoesters and ether lipids. The mammalian enzymes that synthesize fatty alcohols have not been identified. Here, an in silico approach was used to discern two putative reductase enzymes designated FAR1 and FAR2. Expression studies in intact cells showed that FAR1 and FAR2 cDNAs encoded isozymes that reduced fatty acids to fatty alcohols. Fatty acyl-CoA esters were the substrate of FAR1, and the enzyme required NADPH as a cofactor. FAR1 preferred saturated and unsaturated fatty acids of 16 or 18 carbons as substrates, whereas FAR2 preferred saturated fatty acids of 16 or 18 carbons. Confocal light microscopy indicated that FAR1 and FAR2 were localized in the peroxisome. The FAR1 mRNA was detected in many mouse tissues with the highest level found in the preputial gland, a modified sebaceous gland. The FAR2 mRNA was more restricted in distribution and most abundant in the eyelid, which contains wax-laden meibomian glands. Both FAR mRNAs were present in the brain, a tissue rich in ether lipids. The data suggest that fatty alcohol synthesis in mammals is accomplished by two fatty acyl-CoA reductase isozymes that are expressed at high levels in tissues known to synthesize wax monoesters and ether lipids.

Isolation and characterization of PDE8A, a novel human cAMP-specific phosphodiesterase

As part of our efforts to understand the regulation of intracellular cAMP and to generate new targets for pharmacological intervention, we have cloned and characterized the first isozyme in a new family of cyclic nucleotide phosphodiesterases, PDE8A. PDE8A is most similar to PDE4 (38.5% amino acid identity in the catalytic domain), but is clearly not a member of any of the seven known PDE families. We report the cloning of human cDNA encoding the C-terminal 713 amino acids of the protein, including a 283 amino acid region located near the C-terminus that is homologous to the approximately 270 amino acid catalytic domain of other PDEs. In addition, we found cDNA sequences consistent with alternative 5' mRNA splicing analogous to that seen in other PDE genes. PDE8A is expressed in a wide variety of tissues as a approximately 4.5 kb mRNA, with highest levels in testis, ovary, small intestine, and colon. The C-terminal 545 amino acids of PDE8A (the region shared among all splice variants) were expressed in baculovirus. Kinetic analysis of the baculovirus expressed enzyme shows it to be a very high affinity cAMP specific PDE with a Km of 55 nM for cAMP and 124 microM for cGMP. The Vmax (150 pmol/min/microgram recombinant enzyme) is about 10 times slower than that of PDE4. The cAMP hydrolytic activity of PDE8A is not modulated by cGMP at concentrations up to 100 microM. The enzyme requires the presence of at least 1 mM Mn2+ or Mg2+ for maximal activity in vitro, while 100 mM Ca2+ restores only about 20% activity. PDE8A is insensitive (up to 100 microM) to a variety of PDE inhibitors including rolipram, zaprinast, vinpocetine, SKF-94120, and IBMX, but is inhibited (IC50 = 9 microM) by the PDE inhibitor dipyridimole. To give PDE8A a descriptive name that distinguishes it from the other two known high affinity cAMP-specific phosphodiesterases (PDE4 and PDE7), we denote PDE8A as the high affinity cAMP-specific and IBMX-insensitive PDE.

Type 4 phosphodiesterase inhibitors have clinical and in vitro anti-inflammatory effects in atopic dermatitis

Increased cyclic AMP-phosphodiesterase activity in peripheral blood leukocytes is associated with the immune and inflammatory hyperreactivity that characterizes atopic dermatitis. Atopic phosphodiesterase has high sensitivity to a variety of enzyme inhibitors, suggesting an increased therapeutic advantage. The objective of this study was to use in vitro assays to identify a potent phosphodiesterase inhibitor and then to investigate its effectiveness in treating atopic dermatitis. Leukocyte enzyme activity was measured by radioenzyme assay, whereas prostaglandin E2 and interleukins 10 (IL-10) and 4 (IL-4) were measured in 24-h culture supernatants of mononuclear leukocytes by immunoassays. The effect of a topical phosphodiesterase inhibitor on atopic dermatitis lesional skin was assessed by double-blind, paired comparisons of active drug and placebo ointments applied to symmetrically involved sites over a 28-d period. Using in vitro, assays, we demonstrated the ability of selective high-potency phosphodiesterase inhibitors to reduce prostaglandin E2, IL-10, and IL-4 production in atopic mononuclear leukocyte cultures. We selected the Type 4 phosphodiesterase inhibitor, CP80,633, based on its inhibitory potency, for clinical testing by topical, bilateral paired comparisons in 20 patients with atopic dermatitis and demonstrated significant reductions of all inflammatory parameters. Phosphodiesterase inhibitors modulate several pathways contributing to the exaggerated immune and inflammatory responses, which characterize atopic dermatitis. This in vivo demonstration of anti-inflammatory efficacy may provide a useful alternative to the over-reliance on corticosteroid therapy in atopic disease.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Simultaneous determination of aflatoxin M1, ochratoxin A, zearalenone and α-zearalenol in milk by UHPLC-MS/MS

In this study, a sensitive and rapid method has been developed for the simultaneous determination of aflatoxin M1, ochratoxin A, zearalenone and α-zearalenol in milk by ultra high performance liquid chromatography combined with electrospray ionisation triple quadrupole tandem mass spectrometry (UHPLC-ESI-MS/MS). The milk samples were purified using Oasis HLB cartridge. The matrix effects were evaluated by determining the signal suppression-enhancement (SSE) and corrected by external matrix-matched calibration. The limits of quantity (LOQ) of the mycotoxins were in the range of 0.003-0.015μgkg(-1). The high correlation coefficients (R(2)⩾0.996) were obtained in the range of 0.01-1.00μgkg(-1) of the mycotoxins, along with good recovery (87.0-109%), repeatability (3.4-9.9%) and intra-laboratory reproducibility (4.0-9.9%) at the concentrations of 0.025, 0.1 and 0.5μgkg(-1). The detected rates of the mycotoxins were from 16.7% to 96.7% in raw milk, liquid milk and milk powder samples collected from the dairy farms and supermarkets in Beijing. The method proposed is suitable for the simultaneous determination of aflatoxin M1, ochratoxin A, zearalenone, and α-zearalenol, and could be performed for analysing the mycotoxins in milk.

Biarylcarboxylic acids and -amides: inhibition of phosphodiesterase type IV versus [3H]rolipram binding activity and their relationship to emetic behavior in the ferret

In addition to having desirable inhibitory effects on inflammation, anaphylaxis, and smooth muscle contraction, PDE-IV inhibitors also produce undesirable side effects including nausea and vomiting. In general, compounds that inhibit PDE-IV also potently displace [3H]rolipram from a high-affinity binding site in rat cortex. While this binding site has not been identified, it has been proposed to be an allosteric binding site on the PDE-IV enzyme. Preliminary studies have suggested that the emetic potency of PDE-IV inhibitors is correlated with affinity for the brain rolipram binding site rather than potency at inhibiting PDE-IV enzyme activity. Efforts to eliminate the emetic potential of PDE-IV inhibitors were directed toward developing compounds with decreased [3H]rolipram binding affinity while retaining PDE-IV potency. Thus, a novel series of 4-(3-alkoxy-4-methoxyphenyl)benzoic acids and their corresponding carboxamides were prepared and evaluated for their PDE-IV inhibitory and rolipram binding site properties. Modification of the catechol ether moiety led to phenylbutoxy and phenylpentoxy analogues that provided the desired activity profile. Specifically, 4-[3-(5-phenylpentoxy)-4-methoxyphenyl]-2-methylbenzoic acid, 18, was found to exhibit potent PDE-IV inhibitory activity (IC50 0.41 microM) and possessed 400 times weaker activity than rolipram for the [3H]rolipram binding site. In vivo, compound 18 was efficacious in the guinea pig aerosolized antigen induced airway obstruction assay (ED50 8.8 mg/kg, po) and demonstrated a significant reduction in emetic side effects (ferret, 20% emesis at 30 mg/kg, po).

Estimating absolute methylation levels at single-CpG resolution from methylation enrichment and restriction enzyme sequencing methods

Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to map complete DNA methylomes. These include whole-genome bisulfite sequencing (WGBS, MethylC-seq, or BS-seq), reduced-representation bisulfite sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. These methods yield largely comparable results but differ significantly in extent of genomic CpG coverage, resolution, quantitative accuracy, and cost, at least while using current algorithms to interrogate the data. None of these existing methods provides single-CpG resolution, comprehensive genome-wide coverage, and cost feasibility for a typical laboratory. We introduce methylCRF, a novel conditional random fields–based algorithm that integrates methylated DNA immunoprecipitation (MeDIP-seq) and methylation-sensitive restriction enzyme (MRE-seq) sequencing data to predict DNA methylation levels at single-CpG resolution. Our method is a combined computational and experimental strategy to produce DNA methylomes of all 28 million CpGs in the human genome for a fraction (<10%) of the cost of whole-genome bisulfite sequencing methods. methylCRF was benchmarked for accuracy against Infinium arrays, RRBS, WGBS sequencing, and locus-specific bisulfite sequencing performed on the same human embryonic stem cell line. methylCRF transformation of MeDIP-seq/MRE-seq was equivalent to a biological replicate of WGBS in quantification, coverage, and resolution. We used conventional bisulfite conversion, PCR, cloning, and sequencing to validate loci where our predictions do not agree with whole-genome bisulfite data, and in 11 out of 12 cases, methylCRF predictions of methylation level agree better with validated results than does whole-genome bisulfite sequencing. Therefore, methylCRF transformation of MeDIP-seq/MRE-seq data provides an accurate, inexpensive, and widely accessible strategy to create full DNA methylomes.

Factors affecting the lactoferrin concentration in bovine milk

Lactoferrin (LF) concentrations in the milk with different levels of the somatic cell count score were examined using an ELISA to determine whether milk LF concentration is influenced by parity of the cow, stage of lactation, and the somatic cell count. The study animals were 198 Chinese Holstein cows randomly chosen from more than 1,600 cows in 4 dairy farms in the Beijing area. The cows had shown no sign of mastitis for 2 mo. Daily milk production was recorded, and milk samples were taken from individual cow samples. The LF concentration varied between 31.78 and 485.63 microg/mL in milk from normal animals. Lactoferrin was significantly associated with stage of lactation (r = 0.557) and daily milk production (r = -0.472). Nevertheless, there was no significant relationship with parity. Moreover, milk LF concentration tended to be correlated with the somatic cell count score (r = 0.375). This finding suggests that milk LF may be helpful as an indicator for intramammary infection in dairy cows.

Glucocorticoids increase pulmonary beta-adrenergic receptors in fetal rabbit

beta-Adrenergic agonists stimulate surfactant release and decrease fluid in lung alveoli of fetuses. Both effects are most evident toward the end of gestation. We used [3H] dihydroalprenolol (DHA) to investigate the development of pulmonary beta-adrenergic receptors in rabbit fetuses and to study the effect of glucocorticoid treatment on the beta-receptor number. In the lung particulate preparation, DHA binding was rapid, reversible, stereoselective, and of high affinity. The order of potency for adrenergic agonists in competing for DHA binding was isoproterenol > epinephrine = norepinephrine, which is typical of interactions at a beta 1-adrenergic receptor. Using DHA, we demonstrated that the concentration of pulmonary beta-receptors increased significantly between 28 and 31 days of gestation; however, there was no change in the dissociation constant during gestation. After injecting betamethasone (0.17 mg/kg, 24 hours) into rabbits at 25 days of pregnancy, we found that the concentration of pulmonary beta-receptors increased from 44.2 +/- 6.6 fmol/mg protein in untreated fetuses to 77.9 +/- 5.6 fmol/mg protein in treated fetuses. However, this treatment did not affect the DHA binding sites in the fetal rabbit heart. Maternal treatment with the T3 analogue 3,5-dimethyl-3'-isopropyl-L-thyronine (0.5-1 mg/kg) at a dosage which increased both surfactant synthesis and release did not alter pulmonary receptor concentration. Our results indicate that the concentration of pulmonary beta-adrenergic receptors increases in the fetus at term and suggest that this increase is stimulated by endogenous glucocorticoid in fetal circulation.

Mammalian wax biosynthesis. II. Expression cloning of wax synthase cDNAs encoding a member of the acyltransferase enzyme family

Wax monoesters are synthesized by the esterification of fatty alcohols and fatty acids. A mammalian enzyme that catalyzes this reaction has not been isolated. We used expression cloning to identify cDNAs encoding a wax synthase in the mouse preputial gland. The wax synthase gene is located on the X chromosome and encodes a member of the acyltransferase family of enzymes that synthesize neutral lipids. Expression of wax synthase in cultured cells led to the formation of wax monoesters from straight chain saturated, unsaturated, and polyunsaturated fatty alcohols and acids. Polyisoprenols also were incorporated into wax monoesters by the enzyme. The wax synthase had little or no ability to synthesize cholesteryl esters, diacylglycerols, or triacylglycerols, whereas other acyltransferases, including the acyl-CoA:monoacylglycerol acyltransferase 1 and 2 enzymes and the acyl-CoA:diacylglycerol acyltransferase 1 and 2 enzymes, exhibited modest wax monoester synthesis activities. Confocal light microscopy indicated that the wax synthase was localized in membranes of the endoplasmic reticulum. Wax synthase mRNA was abundant in tissues rich in sebaceous glands such as the preputial gland and eyelid and was present at lower levels in other tissues. Coexpression of cDNAs specifying fatty acyl-CoA reductase 1 and wax synthase led to the synthesis of wax monoesters. The data suggest that wax monoester synthesis in mammals involves a two step biosynthetic pathway catalyzed by fatty acyl-CoA reductase and wax synthase enzymes.

Transcription restores DNA repair to heterochromatin, determining regional mutation rates in cancer genomes

Somatic mutations in cancer are more frequent in heterochromatic and late-replicating regions of the genome. We report that regional disparities in mutation density are virtually abolished within transcriptionally silent genomic regions of cutaneous squamous cell carcinomas (cSCCs) arising in an XPC(-/-) background. XPC(-/-) cells lack global genome nucleotide excision repair (GG-NER), thus establishing differential access of DNA repair machinery within chromatin-rich regions of the genome as the primary cause for the regional disparity. Strikingly, we find that increasing levels of transcription reduce mutation prevalence on both strands of gene bodies embedded within H3K9me3-dense regions, and only to those levels observed in H3K9me3-sparse regions, also in an XPC-dependent manner. Therefore, transcription appears to reduce mutation prevalence specifically by relieving the constraints imposed by chromatin structure on DNA repair. We model this relationship among transcription, chromatin state, and DNA repair, revealing a new, personalized determinant of cancer risk.

Classification of human chronic inflammatory skin disease based on single-cell immune profiling

Inflammatory conditions represent the largest class of chronic skin disease, but the molecular dysregulation underlying many individual cases remains unclear. Single-cell RNA sequencing (scRNA-seq) has increased precision in dissecting the complex mixture of immune and stromal cell perturbations in inflammatory skin disease states. We single-cell-profiled CD45+ immune cell transcriptomes from skin samples of 31 patients (7 atopic dermatitis, 8 psoriasis vulgaris, 2 lichen planus (LP), 1 bullous pemphigoid (BP), 6 clinical/histopathologically indeterminate rashes, and 7 healthy controls). Our data revealed active proliferative expansion of the Treg and Trm components and universal T cell exhaustion in human rashes, with a relative attenuation of antigen-presenting cells. Skin-resident memory T cells showed the greatest transcriptional dysregulation in both atopic dermatitis and psoriasis, whereas atopic dermatitis also demonstrated recurrent abnormalities in ILC and CD8+ cytotoxic lymphocytes. Transcript signatures differentiating these rash types included genes previously implicated in T helper cell (TH2)/TH17 diatheses, segregated in unbiased functional networks, and accurately identified disease class in untrained validation data sets. These gene signatures were able to classify clinicopathologically ambiguous rashes with diagnoses consistent with therapeutic response. Thus, we have defined major classes of human inflammatory skin disease at the molecular level and described a quantitative method to classify indeterminate instances of pathologic inflammation. To make this approach accessible to the scientific community, we created a proof-of-principle web interface (RashX), where scientists and clinicians can visualize their patient-level rash scRNA-seq-derived data in the context of our TH2/TH17 transcriptional framework.

The innate immune response, clinical outcomes, and ex vivo HCV antiviral efficacy of a TLR7 agonist (PF-4878691)

Hepatitis C virus (HCV) infection is an issue of global concern, and studies are ongoing to identify new therapies that are both effective and safe. PF-4878691 is a Toll-like receptor 7 (TLR7) agonist modeled so as to dissociate its antiviral activities from its inflammatory activities. In a proof-of-mechanism study in healthy volunteers who received doses of 3, 6, and 9 mg of PF-4878691 twice a week for 2 weeks, PF-4878691 induced biomarkers of the immune and interferon (IFN) responses in a dose-dependent and dose-frequency-related manner. A novel finding was induction of TLR7 expression in vivo in response to PF-4878691, leading to an amplified biomarker response. A nonresponder at the 9-mg dose had a polymorphism in the IFN-α receptor 1 subunit (Val168Leu). Two subjects who had received 9-mg doses experienced serious adverse events (SAEs), characterized by flu-like symptoms, hypotension, and lymphopenia, leading to early termination of the study. TLR7 stimulation results in a pharmacologic response at levels commensurate with predicted antiviral efficacy, but these doses are associated with SAEs, raising concerns about the therapeutic window of this class of compounds for the treatment of HCV infection.

Molecular genetics of 3beta-hydroxy-Delta5-C27-steroid oxidoreductase deficiency in 16 patients with loss of bile acid synthesis and liver disease

The 3beta-hydroxy-Delta(5)-C(27)-steroid oxidoreductase (C(27) 3beta-HSD) is a membrane-bound enzyme of the endoplasmic reticulum that catalyzes an early step in the synthesis of bile acids from cholesterol. Subjects with autosomal recessive mutations in the encoding gene, HSD3B7, on chromosome 16p11.2-12 fail to synthesize bile acids and develop a form of progressive liver disease characterized by cholestatic jaundice and malabsorption of lipids and lipid-soluble vitamins from the gastrointestinal tract. The gene encoding the human C(27) 3beta-HSD enzyme was isolated previously, and a 2-bp deletion in exon 6 of HSD3B7 was identified in a well characterized subject with this disorder. Here, we report a molecular analysis of 15 additional patients from 13 kindreds with C(27) 3beta-HSD deficiency. Twelve different mutations were identified in the HSD3B7 gene on chromosome 16p11.2-12. Ten mutations were studied in detail and shown to cause complete loss of enzyme activity and, in two cases, alterations in the size or amount of the transcribed mRNA. Mutations were inherited in homozygous form in 13 subjects from 10 families and compound heterozygous form in four subjects from three families. We conclude that a diverse spectrum of mutations in the HSD3B7 gene underlies this rare form of neonatal cholestasis.

Single-Cell Profiling Reveals Divergent, Globally Patterned Immune Responses in Murine Skin Inflammation

Inflammatory response heterogeneity has impeded high-resolution dissection of diverse immune cell populations during activation. We characterize mouse cutaneous immune cells by single-cell RNA sequencing, after inducing inflammation using imiquimod and oxazolone dermatitis models. We identify 13 CD45⁺ subpopulations, which broadly represent most functionally characterized immune cell types. Oxazolone pervasively upregulates Jak2/Stat3 expression across T cells and antigen-presenting cells (APCs). Oxazolone also induces Il4/Il13 expression in newly infiltrating basophils, and Il4ra and Ccl24, most prominently in APCs. In contrast, imiquimod broadly upregulates Il17/Il22 and Ccl4/Ccl5. A comparative analysis of single-cell inflammatory transcriptional responses reveals that APC response to oxazolone is tightly restricted by cell identity, whereas imiquimod enforces shared programs on multiple APC populations in parallel. These global molecular patterns not only contrast immune responses on a systems level but also suggest that the mechanisms of new sources of inflammation can eventually be deduced by comparison to known signatures.

•

Oxazolone pervasively upregulates Jak2/Stat3 expression across T cells and APCs

•

Il4/Il13 induction in skin by oxazolone is dominated by infiltrating basophils

•

Imiquimod broadly increases Il17/Il22 and Ccl4/Ccl5, extending to non-T cells

•

Oxazolone induces more highly compartmentalized immune cell responses than imiquimod

Comparison of muscarinic and beta adrenergic receptors between bovine peripheral lung and tracheal smooth muscles: a striking difference in the receptor concentration

This study was undertaken to compare the activity of muscarinic and beta adrenergic receptors in bovine peripheral lung to the corresponding receptor activity in tracheal smooth muscle. We used [3H] quinuclidinyl benzilate (QNB) and [3H]dihydroalprenolol (DHA) to measure muscarinic and beta receptor activity, respectively. Binding to QNB and DHA at 25 degrees C was rapid, reversible, saturable and of high affinity. The order of potency for cholinergic and adrenergic agents competing for binding was compatible with muscarinic and beta 2 adrenergic potencies. We found that the concentration of muscarinic receptor binding sites was 37-fold greater in the tracheal muscle preparation (2805 +/- 309 fmol/mg protein) than in the peripheral lung preparation (76 +/- 28 fmol/mg protein). Unlike muscarinic receptors, the lung contained 8-fold higher concentration of the beta adrenergic receptors than did the tracheal muscle (1588 +/- 417 vs. 199 +/- 42 fmol/mg protein). The dissociation constant or the agonist's inhibitory constant (Ki) for either receptor binding site, however, was not significantly different between the two tissues. Furthermore, in vitro contraction studies showed that the response of tracheal muscle strips to methacholine was markedly greater than the response of peripheral lung strips, a finding consistent with the QNB binding result. The muscle but not the peripheral lung strip exhibited a relaxing response to epinephrine. Our data indicate a striking quantitative difference in muscarinic and beta adrenergic receptors between lung tissue and tracheal muscle, and that each receptor in the lung is qualitatively similar to the corresponding receptor in the muscle.

Single-Cell Transcriptomics Reveals Spatial and Temporal Turnover of Keratinocyte Differentiation Regulators

Keratinocyte differentiation requires intricately coordinated spatiotemporal expression changes that specify epidermis structure and function. This article utilizes single-cell RNA-seq data from 22,338 human foreskin keratinocytes to reconstruct the transcriptional regulation of skin development and homeostasis genes, organizing them by differentiation stage and also into transcription factor (TF)–associated modules. We identify groups of TFs characterized by coordinate expression changes during progression from the undifferentiated basal to the differentiated state and show that these TFs also have concordant differential predicted binding enrichment in the super-enhancers previously reported to turn over between the two states. The identified TFs form a core subset of the regulators controlling gene modules essential for basal and differentiated keratinocyte functions, supporting their nomination as master coordinators of keratinocyte differentiation. Experimental depletion of the TFs ZBED2 and ETV4, both predicted to promote the basal state, induces differentiation. Furthermore, our single-cell RNA expression analysis reveals preferential expression of antioxidant genes in the basal state, suggesting keratinocytes actively suppress reactive oxygen species to maintain the undifferentiated state. Overall, our work demonstrates diverse computational methods to advance our understanding of dynamic gene regulation in development.