T-cell receptor repertoire analysis of CD4-positive T cells from blood and an affected organ in an autoimmune mouse model.

One hallmark of some autoimmune diseases is the variability of symptoms among individuals. Organs affected by the disease differ between patients, posing a challenge in diagnosing the affected organs. Although numerous studies have investigated the correlation between T cell antigen receptor (TCR) repertoires and the development of infectious and immune diseases, the correlation between TCR repertoires and variations in disease symptoms among individuals remains unclear. This study aimed to investigate the correlation of TCRα and ß repertoires in blood T cells with the extent of autoimmune signs that varies among individuals. We sequenced TCRα and ß of CD4+ CD44high CD62Llow T cells in the blood and stomachs of mice deficient in autoimmune regulator (Aire) (AIRE KO), a mouse model of human autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy. Data analysis revealed that the degree of similarity in TCR sequences between the blood and stomach varied among individual AIRE KO mice and reflected the extent of T cell infiltration in the stomach. We identified a set of TCR sequences whose frequencies in blood might correlate with extent of the stomach manifestations. Our results propose a potential of using TCR repertoires not only for diagnosing disease development but also for diagnosing affected organs in autoimmune diseases.

Characterization of a Novel Capsid Assembly Modulator for the Treatment of Chronic Hepatitis B Virus Infection.

The standard of care for the treatment of chronic hepatitis B (CHB) is typically lifelong treatment with nucleos(t)ide analogs (NAs), which suppress viral replication and provide long-term clinical benefits. However, infectious virus can still be detected in patients who are virally suppressed on NA therapy, which may contribute to the failure of these agents to cure most CHB patients. Accordingly, new antiviral treatment options are being developed to enhance the suppression of hepatitis B virus (HBV) replication in combination with NAs ("antiviral intensification"). Here, we describe GS-SBA-1, a capsid assembly modulator (CAM) belonging to class CAM-E, that demonstrates potent inhibition of extracellular HBV DNA in vitro (EC50 [50% effective concentration] = 19 nM) in HBV-infected primary human hepatocytes (PHHs) as well as in vivo in an HBV-infected immunodeficient mouse model. GS-SBA-1 has comparable activities across HBV genotypes and nucleos(t)ide-resistant mutants in HBV-infected PHHs. In addition, GS-SBA-1 demonstrated in vitro additivity in combination with tenofovir alafenamide (TAF). The administration of GS-SBA-1 to PHHs at the time of infection prevents covalently closed circular DNA (cccDNA) formation and, hence, decreases HBV RNA and antigen levels (EC50 = 80 to 200 nM). Furthermore, GS-SBA-1 prevents the production of extracellular HBV RNA-containing viral particles in vitro. Collectively, these data demonstrate that GS-SBA-1 is a potent CAM that has the potential to enhance viral suppression in combination with an NA.

Detection of parechovirus-A in hospitalized children with acute lower respiratory infection in Myanmar, 2017-2018.

Parechovirus-A (PeV-A) causes emerging infection in children, and clinical presentation depends on genotype. The virus has been investigated mainly in developed countries; however, data from developing countries, especially in Asia, are sparse. This study investigated whether PeV-A circulated in children in Myanmar. This retrospective study evaluated PeV-A in nasopharyngeal samples from children aged 1 month to 12 years who were hospitalized with acute lower respiratory infection at Yankin Children Hospital, Yangon, Myanmar, during the period from May 2017 to April 2019. Real-time polymerase chain reaction (PCR) was used to detect PeV-A, and PCR-positive samples were used for genotyping and phylogenetic analysis. In total, 11/570 (1.9%) of samples were positive for PeV-A; 7 were successfully genotyped by sequencing the VP3/VP1 region, as follows: PeV-A1 (n = 4), PeV-A5 (n = 1), PeV-A6 (n = 1), and PeV-A14 (n = 1). Median age was 10.0 months (interquartile range 4.0-12.0 months), and other respiratory viruses were detected in all cases. Phylogenetic analysis showed that all detected PeV-A1 strains were in clade 1 A, which was a minor clade worldwide. Four PeV-A genotypes were detected in Myanmar. The clinical impact of PeV-A in children should be evaluated in future studies.

Neonatal Exposure to Lipopolysaccharide Promotes Neurogenesis of Subventricular Zone Progenitors in the Developing Neocortex of Ferrets.

Lipopolysaccharide (LPS) is a natural agonist of toll-like receptor 4 that serves a role in innate immunity. The current study evaluated the LPS-mediated regulation of neurogenesis in the subventricular zone (SVZ) progenitors, that is, the basal radial glia and intermediate progenitors (IPs), in ferrets. Ferret pups were subcutaneously injected with LPS (500 µg/g of body weight) on postnatal days (PDs) 6 and 7. Furthermore, 5-ethynyl-2'-deoxyuridine (EdU) and 5-bromo-2'-deoxyuridine (BrdU) were administered on PDs 5 and 7, respectively, to label the post-proliferative and proliferating cells in the inner SVZ (iSVZ) and outer SVZ (oSVZ). A significantly higher density of BrdU single-labeled proliferating cells was observed in the iSVZ of LPS-exposed ferrets than in controls but not in post-proliferative EdU single-labeled and EdU/BrdU double-labeled self-renewing cells. BrdU single-labeled cells exhibited a lower proportion of Tbr2 immunostaining in LPS-exposed ferrets (22.2%) than in controls (42.6%) and a higher proportion of Ctip2 immunostaining in LPS-exposed ferrets (22.2%) than in controls (8.6%). The present findings revealed that LPS modified the neurogenesis of SVZ progenitors. Neonatal LPS exposure facilitates the proliferation of SVZ progenitors, followed by the differentiation of Tbr2-expressing IPs into Ctip2-expressing immature neurons.

Forward-Backward Sweep Method for the System of HJB-FP Equations in Memory-Limited Partially Observable Stochastic Control.

Memory-limited partially observable stochastic control (ML-POSC) is the stochastic optimal control problem under incomplete information and memory limitation. To obtain the optimal control function of ML-POSC, a system of the forward Fokker-Planck (FP) equation and the backward Hamilton-Jacobi-Bellman (HJB) equation needs to be solved. In this work, we first show that the system of HJB-FP equations can be interpreted via Pontryagin's minimum principle on the probability density function space. Based on this interpretation, we then propose the forward-backward sweep method (FBSM) for ML-POSC. FBSM is one of the most basic algorithms for Pontryagin's minimum principle, which alternately computes the forward FP equation and the backward HJB equation in ML-POSC. Although the convergence of FBSM is generally not guaranteed in deterministic control and mean-field stochastic control, it is guaranteed in ML-POSC because the coupling of the HJB-FP equations is limited to the optimal control function in ML-POSC.

Decentralized Stochastic Control with Finite-Dimensional Memories: A Memory Limitation Approach.

Decentralized stochastic control (DSC) is a stochastic optimal control problem consisting of multiple controllers. DSC assumes that each controller is unable to accurately observe the target system and the other controllers. This setup results in two difficulties in DSC; one is that each controller has to memorize the infinite-dimensional observation history, which is not practical, because the memory of the actual controllers is limited. The other is that the reduction of infinite-dimensional sequential Bayesian estimation to finite-dimensional Kalman filter is impossible in general DSC, even for linear-quadratic-Gaussian (LQG) problems. In order to address these issues, we propose an alternative theoretical framework to DSC-memory-limited DSC (ML-DSC). ML-DSC explicitly formulates the finite-dimensional memories of the controllers. Each controller is jointly optimized to compress the infinite-dimensional observation history into the prescribed finite-dimensional memory and to determine the control based on it. Therefore, ML-DSC can be a practical formulation for actual memory-limited controllers. We demonstrate how ML-DSC works in the LQG problem. The conventional DSC cannot be solved except in the special LQG problems where the information the controllers have is independent or partially nested. We show that ML-DSC can be solved in more general LQG problems where the interaction among the controllers is not restricted.

Memory-Limited Partially Observable Stochastic Control and Its Mean-Field Control Approach.

Control problems with incomplete information and memory limitation appear in many practical situations. Although partially observable stochastic control (POSC) is a conventional theoretical framework that considers the optimal control problem with incomplete information, it cannot consider memory limitation. Furthermore, POSC cannot be solved in practice except in special cases. In order to address these issues, we propose an alternative theoretical framework, memory-limited POSC (ML-POSC). ML-POSC directly considers memory limitation as well as incomplete information, and it can be solved in practice by employing the technique of mean-field control theory. ML-POSC can generalize the linear-quadratic-Gaussian (LQG) problem to include memory limitation. Because estimation and control are not clearly separated in the LQG problem with memory limitation, the Riccati equation is modified to the partially observable Riccati equation, which improves estimation as well as control. Furthermore, we demonstrate the effectiveness of ML-POSC for a non-LQG problem by comparing it with the local LQG approximation.

Lineage EM algorithm for inferring latent states from cellular lineage trees.

SUMMARY: Phenotypic variability in a population of cells can work as the bet-hedging of the cells under an unpredictably changing environment, the typical example of which is the bacterial persistence. To understand the strategy to control such phenomena, it is indispensable to identify the phenotype of each cell and its inheritance. Although recent advancements in microfluidic technology offer us useful lineage data, they are insufficient to directly identify the phenotypes of the cells. An alternative approach is to infer the phenotype from the lineage data by latent-variable estimation. To this end, however, we must resolve the bias problem in the inference from lineage called survivorship bias. In this work, we clarify how the survivorship bias distorts statistical estimations. We then propose a latent-variable estimation algorithm without the survivorship bias from lineage trees based on an expectation-maximization (EM) algorithm, which we call lineage EM algorithm (LEM). LEM provides a statistical method to identify the traits of the cells applicable to various kinds of lineage data. AVAILABILITY AND IMPLEMENTATION: An implementation of LEM is available at https://github.com/so-nakashima/Lineage-EM-algorithm. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Connection between the Bacterial Chemotactic Network and Optimal Filtering.

The chemotactic network of Escherichia coli has been studied extensively both biophysically and information theoretically. Nevertheless, connection between these two aspects is still elusive. In this work, we report such a connection. We derive an optimal filtering dynamics under the assumption that E. coli's sensory system optimally infers the binary information whether it is swimming up or down along an exponential ligand gradient from noisy sensory signals. Then we show that a standard biochemical model of the chemotactic network is mathematically equivalent to this information-theoretically optimal dynamics. Moreover, we demonstrate that an experimentally observed nonlinear response relation can be reproduced from the optimal dynamics. These results suggest that the biochemical network of E. coli chemotaxis is designed to optimally extract the binary information along an exponential gradient in a noisy condition.

Delayed Postnatal Growth and Anterior Pituitary Development in Growth-Retarded (grt) Female Mice.

Growth-retarded (grt) mice display primary congenital hypothyroidism due to the hyporesponsiveness of their thyroid glands to thyroid-stimulating hormone (TSH). We examined somatic growth, anterior pituitary development, and hormonal profiles in female grt mice and normal ones. Although growth in grt females was suppressed 2 weeks after birth, the measured growth parameters and organ weights gradually increased and finally reached close to the normal levels. Grt mice exhibited delayed eye and vaginal openings and remained in a state of persistent diestrus thereafter, plasma estrogen levels being lower than those in normal mice. Grt mice that received normal-donor thyroids showed accelerated growth and their body weights increased up to the sham-normal levels, indicating the importance of early thyroid hormone supplementation. In the anterior pituitary, there were fewer growth hormone (GH) and prolactin (PRL) cells in grt mice than in normal mice as examined at 12 weeks after birth, but the numbers of these cells did not differ from those in normal mice after 24 weeks. Grt mice had more TSH cells than normal mice until 48 weeks. Plasma GH levels in grt mice were lower than those in normal mice at 2 weeks, but did not differ substantially after 5 weeks. Compared with normal mice, grt mice had significantly lower plasma PRL and thyroxine levels, but notably higher TSH levels until 48 weeks. These findings indicate that thyroid hormone deficiency in grt mice causes delayed development and growth, and inappropriate development of GH, PRL and TSH cells, followed by the abnormal secretion of hormones by these pituitary cells.

Forward and Backward Bellman Equations Improve the Efficiency of the EM Algorithm for DEC-POMDP.

Decentralized partially observable Markov decision process (DEC-POMDP) models sequential decision making problems by a team of agents. Since the planning of DEC-POMDP can be interpreted as the maximum likelihood estimation for the latent variable model, DEC-POMDP can be solved by the EM algorithm. However, in EM for DEC-POMDP, the forward-backward algorithm needs to be calculated up to the infinite horizon, which impairs the computational efficiency. In this paper, we propose the Bellman EM algorithm (BEM) and the modified Bellman EM algorithm (MBEM) by introducing the forward and backward Bellman equations into EM. BEM can be more efficient than EM because BEM calculates the forward and backward Bellman equations instead of the forward-backward algorithm up to the infinite horizon. However, BEM cannot always be more efficient than EM when the size of problems is large because BEM calculates an inverse matrix. We circumvent this shortcoming in MBEM by calculating the forward and backward Bellman equations without the inverse matrix. Our numerical experiments demonstrate that the convergence of MBEM is faster than that of EM.

Immunophenotyping of Nonneoplastic and Neoplastic Histiocytes in Cats and Characterization of a Novel Cell Line Derived From Feline Progressive Histiocytosis.

Histiocytic proliferative diseases are rare in cats, and their pathogenesis is poorly understood. In the present study, 25 cases of histiocytic sarcoma (HS) and 6 of feline progressive histiocytosis (FPH) were examined, and survival times were recorded in 19 cases. The immunophenotypes of tumor cells in these cases as well as of nonneoplastic feline histiocytes were characterized using formalin-fixed, paraffin-embedded tissues. An FPH cell line (AS-FPH01) and xenotransplant mouse model of FPH were also established. The median survival time of HS (150 days) was significantly shorter than that of FPH (470 days). Immunohistochemically, nonneoplastic histiocytes were immunopositive for various combinations of Iba-1, HLA-DR, E-cadherin, CD204, CD163, CD208, and MAC387. By immunohistochemistry, dermal interstitial dendritic cells (iDCs) and macrophages were CD204+/E-cadherin-, while epidermal Langerhans cells (LCs) were CD204-/E-cadherin+. Neoplastic cells of 4 FPH and 18 HS were CD204+/E-cadherin- (iDC/macrophage immunophenotype), while 2 FPH and 2 HS were CD204-/E-cadherin+ (LC immunophenotype), and 5 HS were CD204+/E-cadherin+ (LC-like cell immunophenotype). Furthermore, immunohistochemical and western blot analyses of AS-FPH01 cells derived from E-cadherin-negative FPH revealed that cultured cells were immunopositive for both CD204 and E-cadherin in vitro and in vivo. These results indicate that the neoplastic cells of feline HS and FPH were variably positive for iDC/macrophage and LC markers, and their immunophenotype changed in different microenvironments. The novel cell line established in the present study may serve as an experimental model of FPH that will enable further molecular and therapeutic studies on this disease.

Evaluation of a new apparatus combining a mobile device and hand-held roller for the detection of spinal deformities: association with the sum of the right and left angles of trunk inclination.

PURPOSE: To investigate the intra- and inter-observer reliabilities of the newly developed i-Scolioroller for scoliosis screening, and to determine the optimal i-Scolioroller measurement cutoff values for identifying adolescent scoliosis with a Cobb angle ≥ 20°. METHODS: The i-Scolioroller displays the right- and left-side maximum inclination angle (Rmax, Lmax) during the forward bending test (FBT), as well as the angle of trunk inclination (ATI, i.e., whether the Rmax or Lmax is greater). Sum-ATI is defined as the sum of Rmax and Lmax. Intra-class correlation coefficients (ICC) of the ATI and sum-ATI measurements were calculated to analyze the intra- and inter-observer reliabilities for 10 plaster torsos in FBT positions obtained from patients with idiopathic scoliosis. The optimal cutoff values for scoliosis were determined using receiver operating characteristic (ROC) analysis of i-Scolioroller measurements versus Cobb angles obtained from the upright whole-spine radiographs of 112 adolescent outpatients. RESULTS: The intra-observer ICCs for the ATI/sum-ATI for 3 observers were 0.851/0.856, 0.786/0.900, and 0.772/0.796, respectively, while the corresponding inter-observer ICCs for all participants were 0.733/0.745. On ROC analysis, an ATI of 8° was the optimal cutoff value for scoliosis (sensitivity and specificity: 79.2% and 70.0%, respectively). The optimal cutoff value for sum-ATI was 11° (sensitivity and specificity: 86.1% and 82.5%, respectively). The areas under the ROC curves were 0.859 for ATI and 0.908 for sum-ATI. CONCLUSION: The optimal cutoff values for identifying scoliosis using the i-Scolioroller were a combination of 11° for the sum-ATI and 8° for the ATI.

Adult spinal deformity and its relationship with height loss: a 34-year longitudinal cohort study.

BACKGROUND: Age-related height loss is a normal physical change that occurs in all individuals over 50 years of age. Although many epidemiological studies on height loss have been conducted worldwide, none have been long-term longitudinal epidemiological studies spanning over 30 years. This study was designed to investigate changes in adult spinal deformity and examine the relationship between adult spinal deformity and height loss. METHODS: Fifty-three local healthy subjects (32 men, 21 women) from Furano, Hokkaido, Japan, volunteered for this longitudinal cohort study. Their heights were measured in 1983 and again in 2017. Spino-pelvic parameters were compared between measurements obtained in 1983 and 2017. Individuals with height loss were then divided into two groups, those with degenerative spondylosis and those with degenerative lumbar scoliosis, and different characteristics were compared between the two groups. RESULTS: The mean age of the subjects was 44.4 (31-55) years at baseline and 78.6 (65-89) years at the final follow-up. The mean height was 157.4 cm at baseline and 153.6 cm at the final follow-up, with a mean height loss of 3.8 cm over 34.2 years. All parameters except for thoracic kyphosis were significantly different between measurements taken in 1983 and 2017 (p < 0.05). Height loss in both sexes was related to changes in pelvic parameters including pelvic incidence-lumbar lordosis (R = 0.460 p = 0.008 in men, R = 0.553 p = 0.012 in women), pelvic tilt (R = 0.374 p = 0.035 in men, R = 0.540 p = 0.014 in women), and sagittal vertical axis (R = 0.535 p = 0.002 in men, R = 0.527 p = 0.017 in women). Greater height loss was more commonly seen in women (p = 0.001) and in patients with degenerative lumbar scoliosis (p = 0.02). CONCLUSIONS: This longitudinal study revealed that height loss is more commonly observed in women and is associated with adult spinal deformity and degenerative lumbar scoliosis. Height loss is a normal physical change with aging, but excessive height loss is due to spinal kyphosis and scoliosis leading to spinal malalignment. Our findings suggest that height loss might be an early physical symptom for spinal malalignment.

Reliable target ligand detection by noise-induced receptor cluster formation.

Intracellular reactions are intrinsically stochastic. Nonetheless, cells can reliably respond to the changing environment by sensing their target molecules sensitively and specifically, even with the existence of abundant structurally-similar non-target molecules. The mechanism of how the cells can balance and achieve such different characteristics is not yet fully understood. In this work, we demonstrate that these characteristics can be attained by a ligand-induced stochastic cluster formation of receptors via the noise-induced symmetry breaking, in which the intrinsic stochasticity works to enhance sensitivity and specificity. We also show that the noise-induced cluster formation enables cells to detect the target ligand reliably by compensating the abundant non-target ligands in the environment. The proposed mechanism may lead to a deeper understanding of a biological function of the receptor clustering and provide an alternative candidate for the reliable ligand detection to the kinetic proofreading.

Discovery of the pan-genotypic hepatitis C virus NS3/4A protease inhibitor voxilaprevir (GS-9857): A component of Vosevi®.

Treatment of hepatitis C virus (HCV) infection has been historically challenging due the high viral genetic complexity wherein there are eight distinct genotypes and at least 86 viral subtypes. While HCV NS3/4A protease inhibitors are an established treatment option for genotype 1 infection, limited coverage of genotypes 2 and/or 3 combined with serum alanine transaminase (ALT) elevations for some compounds has limited the broad utility of this therapeutic class. Our discovery efforts were focused on identifying an NS3/4A protease inhibitor with pan-genotypic antiviral activity, improved coverage of resistance associated substitutions, and a decreased risk of hepatotoxicity. Towards this goal, distinct interactions with the conserved catalytic triad of the NS3/4A protease were identified that improved genotype 3 antiviral activity. We further discovered that protein adduct formation strongly correlated with clinical ALT elevation for this therapeutic class. Improving metabolic stability and decreasing protein adduct formation through structural modifications ultimately resulted in voxilaprevir. Voxilaprevir, in combination with sofosbuvir and velpatasvir, has demonstrated pan-genotypic antiviral clinical activity. Furthermore, hepatotoxicity was not observed in Phase 3 clinical trials with voxilaprevir, consistent with our design strategy. Vosevi® (sofosbuvir, velpatasvir, and voxilaprevir) is now an approved pan-genotypic treatment option for the most difficult-to-cure individuals who have previously failed direct acting antiviral therapy.

Structure-guided discovery of a novel, potent, and orally bioavailable 3,5-dimethylisoxazole aryl-benzimidazole BET bromodomain inhibitor.

The bromodomain and extra-terminal (BET) family of proteins, consisting of the bromodomains containing protein 2 (BRD2), BRD3, BRD4, and the testis-specific BRDT, are key epigenetic regulators of gene transcription and has emerged as an attractive target for anticancer therapy. Herein, we describe the discovery of a novel potent BET bromodomain inhibitor, using a systematic structure-based approach focused on improving potency, metabolic stability, and permeability. The optimized dimethylisoxazole aryl-benzimidazole inhibitor exhibited high potency towards BRD4 and related BET proteins in biochemical and cell-based assays and inhibited tumor growth in two proof-of-concept preclinical animal models.

α-Ionone, an Apocarotenoid, Induces Plant Resistance to Western Flower Thrips, Frankliniella occidentalis, Independently of Jasmonic Acid.

Apocarotenoids, such as ß-cyclocitral, α-ionone, ß-ionone, and loliolide, are derived from carotenes via chemical or enzymatic processes. Recent studies revealed that ß-cyclocitral and loliolide play an important role in various aspects of plant physiology, such as stress responses, plant growth, and herbivore resistance. However, information on the physiological role of α-ionone is limited. We herein investigated the effects of α-ionone on plant protection against herbivore attacks. The pretreatment of whole tomato (Solanum lycopersicum) plants with α-ionone vapor decreased the survival rate of western flower thrips (Frankliniella occidentalis) without exhibiting insecticidal activity. Exogenous α-ionone enhanced the expression of defense-related genes, such as basic ß-1,3-glucanase and basic chitinase genes, in tomato leaves, but not that of jasmonic acid (JA)- or loliolide-responsive genes. The pretreatment with α-ionone markedly decreased egg deposition by western flower thrips in the JA-insensitive Arabidopsis (Arabidopsis thaliana) mutant coi1-1. We also found that common cutworm (Spodoptera litura) larvae fed on α-ionone-treated tomato plants exhibited a reduction in weight. These results suggest that α-ionone induces plant resistance to western flower thrips through a different mode of action from that of JA and loliolide.

Long-term hindlimb unloading causes a preferential reduction of medullary thymic epithelial cells expressing autoimmune regulator (Aire).

Hindlimb unloading (HU) of rodents has been used as a ground-based model of spaceflight. In this study, we investigated the detailed impact of 14-day HU on the murine thymus. Thymic mass and cell number were significantly reduced after 14 days of hindlimb unloading, which was accompanied by an increment of plasma corticosterone. Although corticosterone reportedly causes selective apoptosis of CD4+CD8+ thymocytes (CD4+CD8+DPs) in mice treated with short-term HU, the reduction of thymocyte cellularity after the 14-day HU was not selective for CD4+CD8+DPs. In addition to the thymocyte reduction, the cellularity of thymic epithelial cells (TECs) was also reduced by the 14-day HU. Flow cytometric and RNA-sequencing analysis suggested that medullary TECs (mTECs) were preferentially reduced after HU. Moreover, immunohistochemical staining suggested that the 14-day HU caused a reduction of the mTECs expressing autoimmune regulator (Aire). Our data suggested that HU impacts both thymocytes and TECs. Consequently, these data imply that thymic T cell repertoire formation could be disturbed during spaceflight-like stress.

Possible involvement of thyrotropin-releasing hormone receptor 3 in the release of prolactin in the metamorphosing bullfrog larvae.

In amphibians, thyrotropin (TSH), corticotropin (ACTH) and prolactin (PRL) are regarded as the major pituitary hormones involved in metamorphosis, their releasing factors being corticotropin-releasing factor (CRF), arginine vasotocin (AVT), and thyrotropin-releasing hormone (TRH), respectively. It is also known that thyrotropes and corticotropes are equipped with CRF type-2 receptor and AVT V1b receptor, respectively. As for PRL cells, information about the type of receptor for TRH (TRHR) through which the action of TRH is mediated to induce the release of PRL is lacking. In order to fill this gap, an attempt was made to characterize the TRHR subtype existing in the PRL cells of the anterior pituitary gland of the bullfrog, Rana catesbeiana. We cloned cDNAs for three types of bullfrog TRHRs, namely TRHR1, TRHR2 and TRHR3, and confirmed that all of them are functional receptors for TRH by means of reporter gene assay. Analyses with semi-quantitative reverse transcription-PCR and in situ hybridization revealed that TRHR3 mRNA is expressed in the anterior lobe and that the signals reside mostly in the PRL cells. It was also noted that the expression levels of TRHR3 mRNA in the anterior pituitary as well as in the PRL cells of metamorphosing tadpoles elevate as metamorphosis progresses. Since the pattern of changes in TRHR3 mRNA levels in the larval pituitary is almost similar to that previously observed in the pituitary PRL mRNA and plasma PRL levels, we provide a view that TRHR3 mediates the action of TRH on the PRL cells to induce the release of PRL that is prerequisite for growth and metamorphosis in amphibians.