Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus.

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Repeated genetic adaptation to altitude in two tropical butterflies.

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to altitude in two tropical butterflies, Heliconius erato and H. melpomene, which have repeatedly and independently adapted to montane habitats on either side of the Andes. We sequenced 518 whole genomes from altitudinal transects and found many regions differentiated between highland (~ 1200 m) and lowland (~ 200 m) populations. We show repeated genetic differentiation across replicate populations within species, including allopatric comparisons. In contrast, there is little molecular parallelism between the two species. By sampling five close relatives, we find that a large proportion of divergent regions identified within species have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a role for both standing genetic variation and gene flow from independently adapted species in promoting parallel local adaptation to the environment.

Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol.

Highly multiplexed approaches have become common in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals using combinatorially barcoded adapters. These strategies, however, can potentially misassigned reads to incorrect samples. Here, we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to 100 multiplexed samples per sequencing lane (639 samples in total). We created two types of sequencing libraries: four libraries of type A, where PCRs were run on individual samples before multiplexing, and three libraries of type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) than type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in type A and B libraries, respectively. We also show that increasing the number of mismatches allowed for barcode rescue to above 2 dramatically increases the number of recovered chimeric reads. We provide recommendations for developing highly multiplexed RAD-seq protocols and analysing the resulting data to minimize the generation of chimeric sequences, allowing their quantification and a finer control on the number of PCR cycles necessary to generate enough input DNA for library preparation.

Genetic studies of human-chimpanzee divergence using stem cell fusions.

Complete genome sequencing has identified millions of DNA changes that differ between humans and chimpanzees. Although a subset of these changes likely underlies important phenotypic differences between humans and chimpanzees, it is currently difficult to distinguish causal from incidental changes and to map specific phenotypes to particular genome locations. To facilitate further genetic study of human-chimpanzee divergence, we have generated human and chimpanzee autotetraploids and allotetraploids by fusing induced pluripotent stem cells (iPSCs) of each species. The resulting tetraploid iPSCs can be stably maintained and retain the ability to differentiate along ectoderm, mesoderm, and endoderm lineages. RNA sequencing identifies thousands of genes whose expression differs between humans and chimpanzees when assessed in single-species diploid or autotetraploid iPSCs. Analysis of gene expression patterns in interspecific allotetraploid iPSCs shows that human-chimpanzee expression differences arise from substantial contributions of both cis-acting changes linked to the genes themselves and trans-acting changes elsewhere in the genome. To enable further genetic mapping of species differences, we tested chemical treatments for stimulating genome-wide mitotic recombination between human and chimpanzee chromosomes, and CRISPR methods for inducing species-specific changes on particular chromosomes in allotetraploid cells. We successfully generated derivative cells with nested deletions or interspecific recombination on the X chromosome. These studies confirm an important role for the X chromosome in trans regulation of expression differences between species and illustrate the potential of this system for more detailed cis and trans mapping of the molecular basis of human and chimpanzee evolution.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species.

Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present "haplotagging," a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Rapid genotype imputation from sequence with reference panels.

Inexpensive genotyping methods are essential to modern genomics. Here we present QUILT, which performs diploid genotype imputation using low-coverage whole-genome sequence data. QUILT employs Gibbs sampling to partition reads into maternal and paternal sets, facilitating rapid haploid imputation using large reference panels. We show this partitioning to be accurate over many megabases, enabling highly accurate imputation close to theoretical limits and outperforming existing methods. Moreover, QUILT can impute accurately using diverse technologies, including long reads from Oxford Nanopore Technologies, and a new form of low-cost barcoded Illumina sequencing called haplotagging, with the latter showing improved accuracy at low coverages. Relative to DNA genotyping microarrays, QUILT offers improved accuracy at reduced cost, particularly for diverse populations that are traditionally underserved in modern genomic analyses, with accuracy nearly doubling at rare SNPs. Finally, QUILT can accurately impute (four-digit) human leukocyte antigen types, the first such method from low-coverage sequence data.

Calcium-Prolactin Secretion Coupling in Rat Pituitary Lactotrophs Is Controlled by PI4-Kinase Alpha.

The role of calcium, but not of other intracellular signaling molecules, in the release of pituitary hormones by exocytosis is well established. Here, we analyzed the contribution of phosphatidylinositol kinases (PIKs) to calcium-driven prolactin (PRL) release in pituitary lactotrophs: PI4Ks - which control PI4P production, PIP5Ks - which synthesize PI(4, 5)P2 by phosphorylating the D-5 position of the inositol ring of PI4P, and PI3KCs - which phosphorylate PI(4, 5)P2 to generate PI(3, 4, 5)P3. We used common and PIK-specific inhibitors to evaluate the strength of calcium-secretion coupling in rat lactotrophs. Gene expression was analyzed by single-cell RNA sequencing and qRT-PCR analysis; intracellular and released hormones were assessed by radioimmunoassay and ELISA; and single-cell calcium signaling was recorded by Fura 2 imaging. Single-cell RNA sequencing revealed the expression of Pi4ka, Pi4kb, Pi4k2a, Pi4k2b, Pip5k1a, Pip5k1c, and Pik3ca, as well as Pikfyve and Pip4k2c, in lactotrophs. Wortmannin, a PI3K and PI4K inhibitor, but not LY294002, a PI3K inhibitor, blocked spontaneous action potential driven PRL release with a half-time of ~20 min when applied in 10 µM concentration, leading to accumulation of intracellular PRL content. Wortmannin also inhibited increase in PRL release by high potassium, the calcium channel agonist Bay K8644, and calcium mobilizing thyrotropin-releasing hormone without affecting accompanying calcium signaling. GSK-A1, a specific inhibitor of PI4KA, also inhibited calcium-driven PRL secretion without affecting calcium signaling and Prl expression. In contrast, PIK93, a specific inhibitor of PI4KB, and ISA2011B and UNC3230, specific inhibitors of PIP5K1A and PIP5K1C, respectively, did not affect PRL release. These experiments revealed a key role of PI4KA in calcium-secretion coupling in pituitary lactotrophs downstream of voltage-gated and PI(4, 5)P2-dependent calcium signaling.

Population structure of Apodemus flavicollis and comparison to Apodemus sylvaticus in northern Poland based on RAD-seq.

BACKGROUND: Mice of the genus Apodemus are one the most common mammals in the Palaearctic region. Despite their broad range and long history of ecological observations, there are no whole-genome data available for Apodemus, hindering our ability to further exploit the genus in evolutionary and ecological genomics context. RESULTS: Here we present results from the double-digest restriction site-associated DNA sequencing (ddRAD-seq) on 72 individuals of A. flavicollis and 10 A. sylvaticus from four populations, sampled across 500 km distance in northern Poland. Our data present clear genetic divergence of the two species, with average p-distance, based on 21377 common loci, of 1.51% and a mutation rate of 0.0011 - 0.0019 substitutions per site per million years. We provide a catalogue of 117 highly divergent loci that enable genetic differentiation of the two species in Poland and to a large degree of 20 unrelated samples from several European countries and Tunisia. We also show evidence of admixture between the three A. flavicollis populations but demonstrate that they have negligible average population structure, with largest pairwise FST<0.086. CONCLUSION: Our study demonstrates the feasibility of ddRAD-seq in Apodemus and provides the first insights into the population genomics of the species.

A major locus controls a biologically active pheromone component in Heliconius melpomene.

Understanding the production, response, and genetics of signals used in mate choice can inform our understanding of the evolution of both intraspecific mate choice and reproductive isolation. Sex pheromones are important for courtship and mate choice in many insects, but we know relatively little of their role in butterflies. The butterfly Heliconius melpomene uses a complex blend of wing androconial compounds during courtship. Electroantennography in H. melpomene and its close relative Heliconius cydno showed that responses to androconial extracts were not species specific. Females of both species responded equally strongly to extracts of both species, suggesting conservation of peripheral nervous system elements across the two species. Individual blend components provoked little to no response, with the exception of octadecanal, a major component of the H. melpomene blend. Supplementing octadecanal on the wings of octadecanal-rich H. melpomene males led to an increase in the time until mating, demonstrating the bioactivity of octadecanal in Heliconius. Using quantitative trait locus (QTL) mapping, we identified a single locus on chromosome 20 responsible for 41% of the parental species' difference in octadecanal production. This QTL does not overlap with any of the major wing color or mate choice loci, nor does it overlap with known regions of elevated or reduced FST . A set of 16 candidate fatty acid biosynthesis genes lies underneath the QTL. Pheromones in Heliconius carry information relevant for mate choice and are under simple genetic control, suggesting they could be important during speciation.

An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice.

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

Characterization of a non-sexual population of Strongyloides stercoralis with hybrid 18S rDNA haplotypes in Guangxi, Southern China.

Strongyloidiasis is a much-neglected but sometimes fatal soil born helminthiasis. The causing agent, the small intestinal parasitic nematode Strongyloides stercoralis can reproduce sexually through the indirect/heterogonic life cycle, or asexually through the auto-infective or the direct/homogonic life cycles. Usually, among the progeny of the parasitic females both, parthenogenetic parasitic (females only) and sexual free-living (females and males) individuals, are present simultaneously. We isolated S. stercoralis from people living in a village with a high incidence of parasitic helminths, in particular liver flukes (Clonorchis sinensis) and hookworms, in the southern Chinese province Guangxi. We determined nuclear and mitochondrial DNA sequences of individual S. stercoralis isolated from this village and from close by hospitals and we compared these S. stercoralis among themselves and with selected published S. stercoralis from other geographic locations. For comparison, we also analyzed the hookworms present in the same location. We found that, compared to earlier studies of S. stercoralis populations in South East Asia, all S. stercoralis sampled in our study area were very closely related, suggesting a recent common source of infection for all patients. In contrast, the hookworms from the same location, while all belonging to the species Necator americanus, showed rather extensive genetic diversity even within host individuals. Different from earlier studies conducted in other geographic locations, almost all S. stercoralis in this study appeared heterozygous for different sequence variants of the 18S rDNA hypervariable regions (HVR) I and IV. In contrast to earlier investigations, except for three males, all S. stercoralis we isolated in this study were infective larvae, suggesting that the sampled population reproduces predominantly, if not exclusively through the clonal life cycles. Consistently, whole genome sequencing of individual worms revealed higher heterozygosity than reported earlier for likely sexual populations of S. stercoralis. Elevated heterozygosity is frequently associated with asexual clonal reproduction.

Distinct Expression Patterns of Osteopontin and Dentin Matrix Protein 1 Genes in Pituitary Gonadotrophs.

Cell-matrix interactions play important roles in pituitary development, physiology, and pathogenesis. In other tissues, a family of non-collagenous proteins, termed SIBLINGs, are known to contribute to cell-matrix interactions. Anterior pituitary gland expresses two SIBLING genes, Dmp1 (dentin matrix protein-1) and Spp1 (secreted phosphoprotein-1) encoding DMP1 and osteopontin proteins, respectively, but their expression pattern and roles in pituitary functions have not been clarified. Here we provide novel evidence supporting the conclusion that Spp1/osteopontin, like Dmp1/DMP1, are expressed in gonadotrophs in a sex- and age-specific manner. Other anterior pituitary cell types do not express these genes. In contrast to Dmp1, Spp1 expression is higher in males; in females, the expression reaches the peak during the diestrus phase of estrous cycle. In further contrast to Dmp1 and marker genes for gonadotrophs, the expression of Spp1 is not regulated by gonadotropin-releasing hormone in vivo and in vitro. However, Spp1 expression increases progressively after pituitary cell dispersion in both female and male cultures. We may speculate that gonadotrophs signal to other pituitary cell types about changes in the structure of pituitary cell-matrix network by osteopontin, a function consistent with the role of this secretory protein in postnatal tissue remodeling, extracellular matrix reorganization after injury, and tumorigenesis.

Genetic mapping of species differences via in vitro crosses in mouse embryonic stem cells.

Discovering the genetic changes underlying species differences is a central goal in evolutionary genetics. However, hybrid crosses between species in mammals often suffer from hybrid sterility, greatly complicating genetic mapping of trait variation across species. Here, we describe a simple, robust, and transgene-free technique to generate "in vitro crosses" in hybrid mouse embryonic stem (ES) cells by inducing random mitotic cross-overs with the drug ML216, which inhibits the DNA helicase Bloom syndrome (BLM). Starting with an interspecific F1 hybrid ES cell line between the Mus musculus laboratory mouse and Mus spretus (∼1.5 million years of divergence), we mapped the genetic basis of drug resistance to the antimetabolite tioguanine to a single region containing hypoxanthine-guanine phosphoribosyltransferase (Hprt) in as few as 21 d through "flow mapping" by coupling in vitro crosses with fluorescence-activated cell sorting (FACS). We also show how our platform can enable direct study of developmental variation by rederiving embryos with contribution from the recombinant ES cell lines. We demonstrate how in vitro crosses can overcome major bottlenecks in mouse complex trait genetics and address fundamental questions in evolutionary biology that are otherwise intractable through traditional breeding due to high cost, small litter sizes, and/or hybrid sterility. In doing so, we describe an experimental platform toward studying evolutionary systems biology in mouse and potentially in human and other mammals, including cross-species hybrids.

The relationship between basal and regulated Gnrhr expression in rodent pituitary gonadotrophs.

Hypothalamic GnRH together with gonadal steroids and activins/inhibin regulate its receptor gene (Gnrhr) expression in vivo, which leads to crucial changes in GnRHR numbers on the plasma membrane. This is accompanied by alterations in the gonadotroph sensitivity and responsiveness during physiologically relevant situations. Here we investigated basal and GnRH-regulated Gnrhr expression in rodent pituitary gonadotrophs in vitro. In pituitary cells from adult animals cultured in the absence of GnRH and steroid hormones, the Gnrhr expression was progressively reduced but not completely abolished. The basal Gnrhr expression was also operative in LßT2 immortalized gonadotrophs never exposed to GnRH. In both cell types, basal transcription was sufficient for the expression of functional GnRHRs. Continuous application of GnRH transiently elevated the Gnrhr expression in cultured pituitary cells followed by a sustained fall without affecting basal transcription. Both basal and regulated Gnrhr transcriptions were dependent on the protein kinase C signaling pathway. The GnRH-regulated Gnrhr expression was not operative in embryonal pituitary and LßT2 cells and was established neonatally, the sex-specific response patterns were formed at the juvenile-peripubertal stage and there was a strong correlation between basal and regulated gene expression during development. Thus, the age-dependent basal and regulated Gnrhr transcription could account for the initial blockade and subsequent activation of the reproductive system during development.

Spontaneous and CRH-Induced Excitability and Calcium Signaling in Mice Corticotrophs Involves Sodium, Calcium, and Cation-Conducting Channels.

Transgenic mice expressing the tdimer2(12) form of Discosoma red fluorescent protein under control of the proopiomelanocortin gene's regulatory elements are a useful model for studying corticotrophs. Using these mice, we studied the ion channels and mechanisms controlling corticotroph excitability. Corticotrophs were either quiescent or electrically active, with a 22-mV difference in the resting membrane potential (RMP) between the 2 groups. In quiescent cells, CRH depolarized the membrane, leading to initial single spiking and sustained bursting; in active cells, CRH further facilitated or inhibited electrical activity and calcium spiking, depending on the initial activity pattern and CRH concentration. The stimulatory but not inhibitory action of CRH on electrical activity was mimicked by cAMP independently of the presence or absence of arachidonic acid. Removal of bath sodium silenced spiking and hyperpolarized the majority of cells; in contrast, the removal of bath calcium did not affect RMP but reduced CRH-induced depolarization, which abolished bursting electrical activity and decreased the spiking frequency but not the amplitude of single spikes. Corticotrophs with inhibited voltage-gated sodium channels fired calcium-dependent action potentials, whereas cells with inhibited L-type calcium channels fired sodium-dependent spikes; blockade of both channels abolished spiking without affecting the RMP. These results indicate that the background voltage-insensitive sodium conductance influences RMP, the CRH-depolarization current is driven by a cationic conductance, and the interplay between voltage-gated sodium and calcium channels plays a critical role in determining the status and pattern of electrical activity and calcium signaling.

Cell Type-Specific Sexual Dimorphism in Rat Pituitary Gene Expression During Maturation.

The most obvious functional differences between mammalian males and females are related to the control of reproductive physiology and include patterns of GnRH and gonadotropin release, the timing of puberty, sexual and social behavior, and the regulation of food intake and body weight. Using the rat as the best-studied mammalian model for maturation, we examined the expression of major anterior pituitary genes in five secretory cell types of developing males and females. Corticotrophs show comparable Pomc profiles in both sexes, with the highest expression occurring during the infantile period. Somatotrophs and lactotrophs also exhibit no difference in Gh1 and Prl profiles during embryonic to juvenile age but show the amplification of Prl expression in females and Gh1 expression in males during peripubertal and postpubertal ages. Gonadotrophs exhibit highly synchronized Lhb, Fshb, Cga, and Gnrhr expression in both sexes, but the peak of expression occurs during the infantile period in females and at the end of the juvenile period in males. Thyrotrophs also show different developmental Tshb profiles, which are synchronized with the expression of gonadotroph genes in males but not in females. These results indicate the lack of influence of sex on Pomc expression and the presence of two patterns of sexual dimorphism in the expression of other pituitary genes: a time shift in the peak expression during postnatal development, most likely reflecting the perinatal sex-specific brain differentiation, and modulation of the amplitude of expression during late development, which is secondary to the establishment of the hypothalamic-pituitary-gonadal and -thyroid axes.

Paliperidone and aripiprazole differentially affect the strength of calcium-secretion coupling in female pituitary lactotrophs.

Hyperprolactinemia is a common adverse in vivo effect of antipsychotic medications that are used in the treatment of patients with schizophrenia. Here, we compared the effects of two atypical antipsychotics, paliperidone and aripiprazole, on cAMP/calcium signaling and prolactin release in female rat pituitary lactotrophs in vitro. Dopamine inhibited spontaneous cAMP/calcium signaling and prolactin release. In the presence of dopamine, paliperidone rescued cAMP/calcium signaling and prolactin release in a concentration-dependent manner, whereas aripiprazole was only partially effective. In the absence of dopamine, paliperidone stimulated cAMP/calcium signaling and prolactin release, whereas aripiprazole inhibited signaling and secretion more potently but less effectively than dopamine. Forskolin-stimulated cAMP production was facilitated by paliperidone and inhibited by aripiprazole, although the latter was not as effective as dopamine. None of the compounds affected prolactin transcript activity, intracellular prolactin accumulation, or growth hormone secretion. These data indicate that paliperidone has dual hyperprolactinemic actions in lactotrophs i) by preserving the coupling of spontaneous electrical activity and prolactin secretion in the presence of dopamine and ii) by inhibiting intrinsic dopamine receptor activity in the absence of dopamine, leading to enhanced calcium signaling and secretion. In contrast, aripiprazole acts on prolactin secretion by attenuating, but not abolishing, calcium-secretion coupling.

Loss of basal and TRH-stimulated Tshb expression in dispersed pituitary cells.

This study addresses the in vivo and in vitro expression pattern of three genes that are operative in the thyrotroph subpopulation of anterior pituitary cells: glycoprotein α-chain (Cga), thyroid-stimulating hormone ß-chain (Tshb), and TRH receptor (Trhr). In vivo, the expression of Cga and Tshb was robust, whereas the expression of Trhr was low. In cultured pituitary cells, there was a progressive decline in the expression of Cga, Tshb, and Trhr. The expression of Tshb could not be reversed via pulsatile or continuous TRH application in variable concentrations and treatment duration or by the removal of thyroid and steroid hormones from the sera. In parallel, the expression of CGA and TSHB proteins declined progressively in pituitary cells from both sexes. The lack of the effect of TRH on Tshb expression was not related to the age of pituitary cultures and the presence of functional TRH receptors. In cultured pituitary fragments, there was also a rapid decline in expression of these genes, but TRH was able to induce transient Tshb expression. In vivo, thyrotrophs were often in close proximity to each other and to somatotroph and folliculostellate cell networks and especially to the lactotroph cell network; such an organization pattern was lost in vitro. These observations suggest that the lack of influence of anterior pituitary architecture and/or intrapituitary factors probably accounts for the loss of basal and TRH-stimulated Tshb expression in dispersed pituitary cells.

Female-specific induction of rat pituitary dentin matrix protein-1 by GnRH.

Hypothalamic GnRH is the primary regulator of reproduction in vertebrates, acting via the G protein-coupled GnRH receptor (GnRHR) in pituitary gonadotrophs to control synthesis and release of gonadotropins. To identify elements of the GnRHR-coupled gene network, GnRH was applied in a pulsatile manner for 6 hours to a mixed population of perifused pituitary cells from cycling females, mRNA was extracted, and RNA sequencing analysis was performed. This revealed 83 candidate-regulated genes, including a large number coding for secreted proteins. Most notably, GnRH induces a greater than 600-fold increase in expression of dentin matrix protein-1 (Dmp1), one of five members of the small integrin-binding ligand N-linked glycoprotein gene family. The Dmp1 response is mediated by the GnRHR, not elicited by other hypothalamic releasing factors, and is approximately 20-fold smaller in adult male pituitary cells. The sex-dependent Dmp1 response is established during the peripubertal period and independent of the developmental pattern of Gnrhr expression. In vitro, GnRH-induced expression of this gene is coupled with release of DMP1 in extracellular medium through the regulated secretory pathway. In vivo, pituitary Dmp1 expression in identified gonadotrophs is elevated after ovulation. Cell signaling studies revealed that the GnRH induction of Dmp1 is mediated by the protein kinase C signaling pathway and reflects opposing roles of ERK1/2 and p38 MAPK; in addition, the response is facilitated by progesterone. These results establish that DMP1 is a novel secretory protein of female rat gonadotrophs, the synthesis and release of which are controlled by the hypothalamus through the GnRHR signaling pathway. This advance raises intriguing questions about the intrapituitary and downstream effects of this new player in GnRH signaling.

The role of cyclic nucleotides in pituitary lactotroph functions.

Lactotrophs are one of the five secretory anterior pituitary cell types specialized to synthesize and release prolactin. In vitro, these cells fire action potentials (APs) spontaneously and the accompanied Ca(2+) transients are of sufficient amplitude to keep the exocytotic pathway, the transcription of prolactin gene, and de novo hormone synthesis continuously active. Basal cyclic nucleotide production is also substantial in cultured cells but not critical for the APs secretion/transcription coupling in lactotrophs. However, elevated intracellular cAMP levels enhance the excitability of lactotrophs by stimulating the depolarizing non-selective cationic hyperpolarization-activated cyclic nucleotide-regulated and background channels, whereas cGMP inhibits it by activating Ca(2+)-controlled K(+) channels. Elevated cAMP also modulates prolactin release downstream of Ca(2+) influx by changing the kinetic of secretory pores: stimulate at low and inhibit at high concentrations. Induction of prolactin gene and lactotroph proliferation is also stimulated by elevated cAMP through protein kinase A. Together, these observations suggest that in lactotrophs cAMP exhibits complex regulatory effects on voltage-gated Ca(2+) influx and Ca(2+)-dependent cellular processes.