The molecular phenotype of kisspeptin neurons in the medial amygdala of female mice

Reproduction is regulated through the hypothalamic-pituitary-gonadal (HPG) axis, largely via the action of kisspeptin neurons in the hypothalamus. Importantly, Kiss1 neurons have been identified in other brain regions, including the medial amygdala (MeA). Though the MeA is implicated in regulating aspects of both reproductive physiology and behavior, as well as non-reproductive processes, the functional roles of MeA Kiss1 neurons are largely unknown. Additionally, besides their stimulation by estrogen, little is known about how MeA Kiss1 neurons are regulated. Using a RiboTag mouse model in conjunction with RNA-seq, we examined the molecular profile of MeA Kiss1 neurons to identify transcripts that are co-expressed in MeA Kiss1 neurons of female mice and whether these transcripts are modulated by estradiol (E2) treatment. RNA-seq identified >13,800 gene transcripts co-expressed in female MeA Kiss1 neurons, including genes for neuropeptides and receptors implicated in reproduction, metabolism, and other neuroendocrine functions. Of the >13,800 genes co-expressed in MeA Kiss1 neurons, only 45 genes demonstrated significantly different expression levels due to E2 treatment. Gene transcripts such as Kiss1, Gal, and Oxtr increased in response to E2 treatment, while fewer transcripts, such as Esr1 and Cyp26b1, were downregulated by E2. Dual RNAscope and immunohistochemistry was performed to validate co-expression of MeA Kiss1 with Cck and Cartpt. These results are the first to establish a profile of genes actively expressed by MeA Kiss1 neurons, including a subset of genes regulated by E2, which provides a useful foundation for future investigations into the regulation and function of MeA Kiss1 neurons.

Hypothalamic Kiss1 neurons in the AVPV and ARC have been well-studied and are known to be regulated by testosterone (T) and estradiol (E 2 ) (13-15). In the ARC, sex steroids inhibit ARC Kiss1 expression, whereas the removal of sex steroids increases ARC Kiss1 levels, suggesting that kisspeptin neurons in the ARC are involved in sex steroid negative feedback (13)(14)(15). In contrast, in the AVPV, Kiss1 levels are increased by sex steroids, particularly E 2 (13-15), suggesting that AVPV Kiss1 neurons participate in the E 2-mediated positive feedback that triggers the preovulatory LH surge in females (21,22). E 2 regulation of AVPV and ARC Kiss1 levels occurs specifically via estrogen receptor a (ERa), which is highly expressed in both of these Kiss1 populations (13, 14,18,23,24). Recent RNA-seq studies of AVPV and ARC Kiss1 neurons reported other co-expressed genes in these specific cell populations that also respond to E 2 treatment (25)(26)(27), including hormone receptors such as Pgr, Ghsr, and Npr2 (26,27).
MeA Kiss1 expression is also regulated by sex steroids. Similar to AVPV Kiss1 expression, MeA Kiss1 levels dramatically increase with T or E 2 exposure in both sexes and fall to nearly undetectable levels in the absence of sex steroids (17,18,20,28). However, the nonaromatizable androgen, DHT, has no stimulatory effect on MeA Kiss1 expression (17), which suggests that the observed upregulation of MeA Kiss1 expression by T is mediated by E 2 signaling after aromatization; this possibility is supported by the presence of high aromatase expression in the MeA region (17,(29)(30)(31). As with Kiss1 in the AVPV and ARC (14,23), data from EraKO and ERbKO mice show that the ERa receptor subtype is required for E 2 's upregulation of MeA Kiss1 levels (18). Despite these findings demonstrating that Kiss1 expression in the MeA is potently stimulated by E 2 , the functional significance of this E 2 stimulation is still currently unknown. In contrast to their E 2 -induced upregulation, MeA Kiss1 levels are inhibited by GABA signaling through GABA B R. This is evidenced by very high Kiss1 expression in the MeA, but not the AVPV or ARC, of GABA B R knockout mice of both sexes (19,20). Whether this GABA effect is direct or indirect is still unknown, though MeA Kiss1 neurons are reported to expressed GABA B R as determined with in situ hybridization (19).
At present, E 2 and GABA, acting via ERa and GABA B R respectively, are the only known regulators of MeA Kiss1 neurons. In fact, almost nothing is known about the identities of other receptors, co-transmitters and signaling factors, and transcription factors that are expressed by this specific kisspeptin neuron population, or whether other genes in MeA kisspeptin neurons are also altered by E 2 signaling. This lack of knowledge of the phenotype of MeA Kiss1 neurons has limited our understanding of possible functions of this kisspeptin population. The MeA has numerous behavioral and physiological functions, including effects on puberty, as well as reproductive physiology and behavior (32-37) and other non-reproductive behaviors (38)(39)(40)(41)(42)(43)(44). Classic studies found that lesions of the entire MeA disrupt ovarian cycles and impair E 2 positive feedback in female rodents (35)(36)(37). Conversely, acute electrical stimulation of the MeA region of E 2 -primed ovariectomized (OVX) females induced high LH secretion (34), indicating that the MeA might facilitate E 2 positive feedback of the LH surge. However, the mechanisms by which the MeA influences reproductive hormone release are still unknown, as is the identity of the specific MeA cell types responsible for these effects. Because kisspeptin can potently stimulate GnRH neurons (10,11) and the MeA projects both directly and indirectly to the POA where GnRH neurons reside (45)(46)(47)(48), it remains possible that MeA Kiss1 neurons participate in the regulation of the reproductive axis by acting directly or indirectly on GnRH neurons. Some recent optogenetic and chemogenetic mouse studies have begun to address this possibility (39,(49)(50)(51), but this is still currently not well understood. Thus, at present, the various functional roles of MeA kisspeptin neurons remain unknown, in part due to our limited understanding of the detailed phenotype of those specific kisspeptin neurons.
Recent studies have begun to detail the molecular phenotype of the two kisspeptin populations in the hypothalamus. Using the RiboTag technique coupled with RNA-seq, we recently reported the identification of >13,000 genes that are expressed in AVPV Kiss1 neurons of female mice (25). In that study, we also identified numerous AVPV Kiss1 neuron transcripts that are differentially regulated by E 2 , such as Pgr, Th, Cartpt, and Gal. Two subsequent reports from Hrabovszky and colleagues used a different approach to examine E 2 regulation of RNA transcripts in ARC and AVPV kisspeptin neurons (26,27). These various RNA profiling studies provide insight into potential mechanisms of regulation of hypothalamic Kiss1 neurons and, hence, reproductive status. However, similar RNA-seq analyses of kisspeptin neurons in the MeA have not yet been reported. In the current study, we used the RiboTag technique (52,53), which allows for the isolation of mRNAs that are actively being translated into proteins, along with RNA-seq to identify the molecular phenotype of MeA Kiss1 neurons of female mice under different E 2 conditions.

Animals
We used the RiboTag mouse model (54,55) crossed with Kiss1Cre mice to identify actively-translated gene transcripts that are co-expressed specifically in MeA Kiss1 neurons, following methods previously described (25,54,55). Briefly, RiboTag (Rpl22 HA+ ) mice have a wild-type C-terminal exon floxed on the Rpl22 gene that is followed by three copies of a hemagglutinin (HA) epitope sequence inserted prior to a stop codon (54). Cell-typespecific recombination can be induced by crossing RiboTag mice with a cell type-specific Cre recombinase-expressing mouse line, which leads to Cre-mediated recombination and expression of HA tags on ribosomes only in cells expressing Cre recombinase (54,56). Here, Rpl22 HA+ mice were crossed with Kiss1Cre mice (courtesy of Carol Elias) (57) to generate Kiss1Cre+/Rpl22 HA+/+ female mice to be used for this study. Kiss1Cre+/Rpl22 HA+/+ mice express the HAtagged ribosomes only in Kiss1-expressing cells, permitting isolation of ribosome-associated transcripts from just Kiss1 cells in specific brain regions, such as the MeA. In addition to adult (8-10 weeks old) Kiss1Cre+/Rpl22 HA+/+ females, a small cohort of adult Kiss1Cre-/ Rpl22 HA+/+ females was used as controls. These control mice do not have HA-tagged ribosomes in Kiss1-expressing cells and therefore, no Kiss1 cell-specific RNA should be isolated. Tail DNA was used to genotype mice via polymerase chain reaction (PCR) to confirm genotypes of Kiss1Cre+/Rpl22 HA+/+ and Kiss1Cre-/Rpl22 HA+/+ mice (henceforth referred to as "Kiss Ribo " or "Control" mice, respectively). Additionally, any mice with germline recombination were excluded from the study.
All mice were housed 2-3 per cage (Kiss Ribo and Control mice) or 2-5 per cage (Kiss tdTom mice) in a 12hr:12hr light:dark cycle (lights off at 18:00h), with access to food and water ad libitum. All animal procedures were approved by local IACUC committees at the University of California, San Diego (Kiss Ribo and Control mice) or Albany Medical College (Kiss tdTom mice).

Hormone treatment and tissue processing
MeA Kiss1 expression is known to be stimulated by E 2 (20). Thus, all Kiss Ribo and Control mice were ovariectomized (OVX) at 8 weeks of age, under isoflurane anesthesia, and pre-treated at this time with high dose E 2 (2 mm of 1:30 E2: cholesterol powder) via subcutaneous Silastic capsule for 4 days. This dose of E 2 in mice is known to increase MeA Kiss1 expression, as well as induce E 2 negative feedback to suppress LH levels (13, 14,17,18,20,25). This E 2 pre-treatment was used to drive sufficient Cre expression in MeA Kiss1 cells in all Kiss Ribo mice and promote a high degree of incorporation of HAtagged ribosomes in these neurons, regardless of subsequent E 2 and OVX treatment conditions (25). For consistency between genotypes, Cre-controls were similarly given the E 2 pre-treatment. All pretreatment E 2 implants were removed after 4 days and all mice were given 7 days to wash out any residual circulating E 2 . After the 1-week washout period, half of the mice were re-implanted with a new E 2 Silastic capsule (E 2 group, n = 12 Kiss Ribo ; n=4 Controls), while the remaining females received no additional E 2 treatment and served as the OVX group (n = 12 Kiss Ribo ; n=4 Controls). 5 days after receiving the second E 2 implant (or no implant), all females were euthanized between 11:00h and 14:00h. Blood was collected at this time via retroorbital bleed, and brains were collected fresh frozen on dry ice. Blood serum was assayed for LH concentrations to confirm low LH levels in the E 2 -treated group (indicating proper E 2 negative feedback) and elevated LH levels in the OVX group (indicating lack of E 2 negative feedback). Serum LH was measured via a highly sensitive mouse LH radioimmunoassay performed by the University of Virginia Ligand Assay Core (lower detection limit: 0.04 ng/mL; average reportable range: 0.04-75 ng/mL). As expected, E 2 -treated OVX females had significantly reduced mean LH levels compared to OVX females lacking E 2 (0.22 ± 0.04 ng/mL vs 3.04 ± 0.22 ng/mL, respectively, p<0.05).
Fresh frozen brains were processed for RiboTag immunoprecipitation. The brain was micro-dissected on a coronal plane and 2 consecutive 400 µm thick slices spanning the MeA region were micro-punched bilaterally using a 2 mm diameter sampling tool. To ensure sufficient yield of isolated mRNA following immunoprecipitation, MeA micro-punches from n = 4 mice were pooled for each treatment (E 2 and OVX groups). The total number of pooled samples per group were as follows: n = 3 for both OVX Kiss Ribo and E 2 Kiss Ribo , n = 1 for both OVX and E 2 Cre-controls. Prior to the RiboTag immunoprecipitation, all pooled MeA micro-punch samples were stored at -80°C in 1.7 mL Eppendorf tubes.
For the ISH/IHC co-expression experiments, female Kiss tdTom mice (n=3) were ovariectomized and received a similar E 2 Silastic capsule for 5 days, as was done for the Kiss Ribo mice. Brains were then collected in 4% paraformaldehyde, transferred to 30% sucrose 24 hours later, and stored at 4°C prior to slicing. Brains were then sectioned at 25 µm/slice, and sections containing the MeA mounted on SuperFrost Plus slides (Fisher Scientific), air dried, and stored at -80°C until the assay.

Immunoprecipitation and RNA extraction
RiboTag immunoprecipitation on pooled MeA samples from Kiss Ribo and Control females was performed following published protocols (54,56). Some modifications to the original protocol were performed to maximize isolation of ribosomes and their attached mRNA transcripts from MeA Kiss1 neurons (25). Specifically, pooled MeA samples were homogenized in a homogenization buffer solution (72% H2O, 9.6% NP-40, 9.6% 2M KCl, 3.2% 1.5M Tris -pH 7.4, 1.2% 1M MgCl2, 2% Cyclohexamide 5mg/ml, 1% Protease Inhibitors, 1% heparin 100mg/ml, 0.5% RNAsin, and 0.1% DTT) at 3% weight by volume. The samples were then centrifuged at 10K rpm for 10 minutes at 4°C. After centrifugation, 10% of the lysate was saved in a separate tube as the "Input" sample, which contains mRNA from all cell types present in the MeA micropunch, including but not limited to Kiss1 cells. To store the Input samples, 350µL of lysis buffer (from Qiagen Kit #74034) was added and briefly vortexed, and then the samples flash frozen and stored at -80°C. The remaining lysate was then used for the immunoprecipitation procedure. To precipitate Kiss1 cell ribosomes and their associated RNA transcripts from each lysate sample, 0.25µL of antibody/100µL lysate of Biolegend Purified anti-Ha.11 Epitope Tag Antibody (#910501) was added to the remaining lysate volume and incubated for 2 hours on a gentle sample rotator at 4°C. Magnetic beads (Pierce Protein A/G #88803; 25µL beads/100µL sample) were washed with homogenization buffer (described above). Following the antibody incubation, the sample was transferred to the washed magnetic beads and incubated again at 4°C for 1 hour on gentle rotation. The lysate was removed by placing the sample tubes on a magnetic tube rack and the beads washed 3 times for 10 minutes each on a gentle rotator at 4°C with a high salt buffer (53.4% H2O, 30% 2M KCl, 10% NP-40, 3.3% 1.5M Tris -pH 7.4, 1.2% 1M MgCl2, 2% Cyclohexamide 5mg/ml, and 0.05% DTT), using fresh high salt buffer for each wash. Immediately following the washes, 350µL lysis buffer (Qiagen Kit #74034) was added to each sample and vortexed for 30 seconds. Sample tubes were then secured in a magnetic tube rack and the resulting lysate (termed the "IP" sample) was separated from the magnetic beads and transferred into another tube, flash frozen, and stored at -80°C. The IP samples only contain RNA from the HA-tagged ribosomes, which in this experiment were only present in Kiss1 neurons. Thus, the IP samples from Kiss Ribo mice contained RNA specific to MeA Kiss1 neurons. To extract RNA from the Input and IP samples, the Qiagen RNeasy Plus Micro Kit (#74034) kit was used per kit instructions, and isolated RNA was stored in aliquots at -80°C until RT-PCR and RNA-sequencing.

Examination of Cck co-expression in
Kiss1 neurons using double-label in situ hybridization and immunohistochemistry A double-label ISH/IHC co-detection assay was performed to examine the co-expression of Cck, a gene known to be expressed in the MeA (59, 60), in Kiss1 neurons. This assay was performed on MeA-containing brain slices from adult Kiss tdTom female mice (1 slice/mouse, n = 3 mice total). The assay measured co-expression of Cck mRNA (cholecystokinin) in MeA Kiss1 neurons (labeled with tdtomato). To perform the co-expression assay, we used Advanced Cell Diagnostics' RNAscope ® multiplex fluorescent V2 ISH assay with IHC co-detection, with the following RNAscope ® catalog probe: Mm-Cck-C3 -Mus musculus cholecystokinin (Cck) mRNA. For IHC detection of tdtomato, Rockland rabbit anti-RFP (red fluorescent protein, #600-401-379) primary antibody was used along with Invitrogen goat anti-rabbit IgG Alexa Fluor 594 (#A-11037) secondary antibody.

RT-PCR confirmation of Kiss1 cell-specific mRNAs
To validate the efficacy of the Kiss1Cre/RiboTag method to isolate mRNA from MeA Kiss1 neurons, we performed RT-PCR on Input and IP samples for Kiss1 and Cck, which is known to be expressed in the MeA where Kiss1 cells are located (59,60) and was identified in the first experiment ( Figure 1) to be co-expressed in Kiss1 neurons. To confirm specificity, we also performed RT-PCR for Avp, a gene known to be expressed in the MeA, but not typically in the posterior MeA where Kiss1 expression is observed (61). RT-PCR was performed using the iScript cDNA Synthesis Kit (Bio-Rad #1708891), per kit instructions, to synthesize cDNA from 10ng of RNA from the Input and IP samples. RedTaq mix (Sigma #P0982), forward and reverse primers, and RNase-free water. The PCR conditions were as follows: 94 x15', (94 x 30", annealing x 30", 72 x 60") repeat 30 times, 72 x 5', 4 x 5'. The annealing temperature was 57°C for Kiss1 and Cck and 62°C for Avp.
2.6 RNA-seq of Cre+ samples and bioinformatic analysis RNA sequencing was performed by the Genomics Center at UC San Diego's Institute for Genomic Medicine, using RNA from the MeA IP samples of Kiss Ribo mice, for both E 2 and OVX (no E 2 ) treatments. Agilent High Sensitivity RNA ScreenTape System was used to determine the quality of all MeA IP RNA samples prior to RNA sequencing, and library preparation was only performed on IP samples with an RNA integrity number >7.6. The library was created using Illumina unstranded mRNA library kits with polyA enrichment. The Illumina HiSeq4000 platform (SR75 run type) was used to perform RNA sequencing. Raw RNA sequencing data quality control analysis was performed using Fastqc v0.11.8, then read trimming was performed using Trimmomatic 0.38, followed by alignment using STAR v2.6.0a, and then quantification of reads (RSEMv1.3.0) using GRCm38.p6/mm10 and Mus-musculus.GRCm38.98.gtf.
The Center for Computational Biology and Bioinformatics at UC San Diego performed all RNA sequencing analyses and statistics. All analyses were limited to known protein coding regions. To prepare for data exploration and preprocessing using edgeR Bioconductor package (62) written in R (63), integration and annotation of sample inputs was performed using per-gene-per-sample counts from both the count preparation and quality control were used with the persample RNA-seq metadata. The edgeR Bioconductor package and limma package (64) were used to explore and pre-process annotated data for determining differential gene expression specifically in transcripts produced by MeA Kiss1 neurons. In order to eliminate low-expressing genes from comparative analyses, only transcripts with a mean CPM >1 for all 3 samples in at least one hormone treatment group (OVX or E 2 ) were used for differential expression (OVX vs E 2 ) and biological KEGG pathway analyses. The voom technique (65,66), which utilizes limma and edgeR Bioconductor packages, was used to determine differential expression of any retained transcripts between the two treatment groups (E 2 versus OVX). To test for the presence of specific genes and differential gene expression in annotated functions, pathways, and diseases, the overall data and differential gene data were analyzed with WebGestalt (67). We then performed a KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis (68) using PathView (69) in both E 2 and OVX groups to examine the primary KEGG pathways observed in the overall MeA Kiss1 RNA sequencing dataset, as well as for any pathways that were represented differently between E 2 and OVX groups. The data set containing all identified gene transcripts is available in the Gene Expression Omnibus (geo) data repository (70).

Validation of gene co-expression using double-label ISH and IHC co-detection
A double-label ISH/IHC co-detection assay was performed to validate the co-expression of Cartpt, a gene found to be highly expressed in our RNA-seq dataset, in MeA Kiss1 neurons. Using brain slices from 3 female Kiss tdTom mice (1 slice/mouse), a doublelabel ISH/IHC for Cartpt mRNA and TdTomato, representing MeA Kiss1 cells, was performed using Advanced Cell Diagnostics' RNAscope ® multiplex fluorescent V2 ISH assay with IHC codetection, as described previously for the detection of Cck in Kiss1 cells. This assay used the RNAscope catalog probe for Cartpt: Mm-Cartpt-C2 -Mus musculus cocaine-and amphetamine-regulated transcript prepropeptide (Cartpt) mRNA.
The manufacturer's protocol was followed for multiplex fluorescent V2 ISH assay with IHC co-detection for fixed frozen tissue, with some modifications to preserve tissue integrity. Briefly, slides containing brain slices with the MeA region were briefly washed with Milli-Q water, baked at 60°C for 1 hour, and post-fixed in fresh 4% paraformaldehyde for 2 hours at 4°C. The tissue was dehydrated in 50%, 70%, and 100% ethanol washes, treated with hydrogen peroxide for 10 minutes, washed with Milli-Q water, and baked for an additional 30 minutes at 60°C. Following baking, the tissue was incubated at 98-102°C for 5 minutes in 1X RNAscope ® target retrieval reagent, and briefly washed with Milli-Q water followed by 1X PBS-T. 200µL primary antibody (1:1000, anti-RFP) was added to each slide and slides were incubated overnight at 4°C. On Day 2, the slides were washed 3 times in 1X PBS-T, transferred to fresh 4% paraformaldehyde for 30 minutes at room temperature, followed by four washes in 1X PBS-T. 2-4 drops of Protease Plus was added to each slide and incubated at 40°C for 30 minutes, followed by 3 brief washes in Milli-Q water. Excess liquid was removed and a probe mix (Probe Diluent + RNAscope probe) was added to each slide and hybridized to the tissue at 40°C for 2 hours. The slides were washed twice with RNAscope ® 1X Wash Buffer at room temperature, followed by 3 amplification hybridization steps (AMP 1 and 2 for 30 minutes at 40°C , and AMP 3 for 15 minutes at 40°C) with two washes with 1X wash buffer between each amplification hybridization step. After hybridization, the HRP-C2 fluorescent signal was developed for each respective probe by adding 2-4 drops of HRP-C2 to each sample and incubating for 15min at 40°C, washing twice with 1X wash buffer, incubating with 150µL diluted Opal ™ dye for 30 minutes at 40°C, and washing twice with 1X wash buffer. This was repeated again, using HRP-C3 (instead of HRP-C2) to develop the HRP-C3 fluorescent signal. The Opal ™ dyes used were: Opal 520 (Cck) and Opal 690 (Cartpt). After developing the HRP-C3 signal, the 200mL secondary antibody (Goat antirabbit), diluted 1:100 in co-detection Antibody Diluent, was incubated on each slide for 30 minutes at room temperature, followed by two washes with PBS-T. After incubating with the secondary antibody, 4 drops of DAPI were added to each slide, incubated for 30 seconds, excess DAPI removed, and the slides then immediately coverslipped using ProLong Gold Antifade Mountant. Slides were stored at 4°C in the dark.

Microscopy analyses of kisspeptin and
Cck or Cartpt co-expression Using a Zeiss LSM 880 confocal microscope, images of Cck or Cartpt mRNA with TdTomato fluorescent staining were obtained at 40X (oil) magnification for one unilateral brain slice containing the MeA for each female. For each female, the number of TdTomato (reporter for kisspeptin) cells were identified (minimum 79 cells/ mouse) and the number of Kiss1 cells that co-expressed Cck or Cartpt were counted, using the RNAscope manufacturer's criteria as a guideline (cell = ≥15 clustered dots). The mean percent of MeA Kiss1 cells that contained Cck or Carpt was calculated.

Statistical analysis
T-tests were performed to compare LH levels between the OVX and E 2 -treated females. To statistically test for differences in gene expression between OVX and E 2 -treated females, the voom technique, which uses limma and edgeR Bioconductor packages (64)(65)(66), was used. This technique using simple linear regression models and produces a modified t-statistic that is interpreted like other t-statistics, with the exception that the standard errors have been moderated across genes using a simple Bayesian model (64)(65)(66). The reported p-value for all RNA-seq statistical analyses is an Benjamini-Hochberg-adjusted p-value to account for the number of genes analyzed in the dataset (71). The statistical significance was set so the adjusted p-value < 0.05.

Validation and specificity of the Ribotag isolation of MeA Kiss1 cell transcripts
Genes that are co-expressed in MeA kisspeptin neurons are unknown, though Cck is known to be expressed in the MeA region (59, 60), suggesting it may also be expressed specifically in kisspeptin neurons in this region. As a search for a positive control gene that is expressed in MeA kisspeptin neurons, we first performed a doublelabel assay to assess if Cck is in fact present in MeA kisspeptin neurons. Indeed, we observed that most MeA Kiss1 cells of female mice also express high levels of Cck (Figure 1), with 88% of MeA Kiss1 cells coexpressing Cck mRNA. Thus, Cck expression was selected as a positive control for subsequent Ribotag pulldown validation purposes.
In order to validate the Ribotag isolation procedure, we performed RT-PCR using cDNA from both the Input samples (mRNA from all cells in the MeA micropunches, including Kiss1 cells) and IP samples (mRNA from only Kiss1 cell ribosomes) to confirm that the immunoprecipitation procedure was specific to Kiss1 cells. As expected, both the Input and IP samples of Kiss Ribo mice contained Kiss1 transcripts ( Figure 2). As also expected, Kiss1 was also found in the Input samples of Control (Cre-) mice, but was not detected in the IP samples of these Controls (Figure 2), confirming that the immuno-pulldown was specific to Kiss1 cells that had undergone Cre-mediated recombination to express the HA+ tag. Our positive control gene, Cck, was also found to be highly expressed in all the Input samples as well as in IP samples of both E 2 -treated and OVX Kiss Ribo mice, but not in the IP of Control mice (Figure 2), further validating the selectivity of the Ribotag isolation technique. Given that little is currently known about what genes are or are not expressed in MeA Kiss1 neurons, it is difficult to identify a true negative control. We selected Avp because Avp is expressed in the MeA region, but typically more in the rostral part of the MeA (61), whereas Kiss1 is expressed in the more caudal (posterior) part of the MeA. Thus, AVP and Kiss1 are in different anatomical sub-regions of the MeA and unlikely to be co-expressed in the same cells, serving as a useful negative control for our Riobotag pulldown selectivity. Indeed, here we confirmed that Avp mRNA was present, at low to moderate levels, in all Input samples but was absent in all IP samples (Figure 2), suggesting that MeA Kiss1 cells do not express Avp but other MeA cell types do.

Identification of RNA transcripts in female MeA Kiss1 cells
The goal of the current study is to identify what neuropeptides, neurotransmitter synthesis and transport factors, and receptors are made by MeA Kiss1 neurons to begin to understand the regulation and function of this Kiss1 neuronal population. RNA quality for all E 2 and OVX Kiss Ribo samples was sufficient for RNA-seq, with RNA integrity values of at least 7.6 for all 6 Kiss Ribo IP samples. As expected, the RNA quality for our Cre-control samples, which lack the HA+ tag in Kiss1 cells, was low, <5.8. Due to the low integrity of the RNA from our control samples, only the Kiss Ribo IP samples were processed for RNA-seq. The library size for each Kiss Ribo sample was ≥ 25M total reads per sample, well above the 5-25M reads per sample recommended by Illumina, and well over the 10M reads aligned reach threshold, indicating a high RNA-seq quality. RNA-seq of the Kiss Ribo samples identified approximately 13,800 different mRNA transcripts, including Kiss1 and Cck, as well as other transcripts such as Esr1, Esr2, Ar, Gal, Cartpt, Pdyn, Oxtr, and Npy2r ( Figure 3). Our first analysis identified the top 75 genes expressed by MeA Kiss1 cells, regardless of E 2 treatment. The top 3 most highly expressed genes (Sptbn1, Sptan1, and Sptbn2) are coding genes for spectrin proteins, which are involved in actin crosslinking, cell communication, and cell regulation. Other very highly expressed transcripts included several genes related to intra-cell signaling, such as Gprasp1, Atp1a3, Calm1, and Ywhag, genes important for protein synthesis and regulation, including Cpe, Hspa8, Hsp90ab1, Eef1a1, and Ubb, and genes involved with secretion and synapses, like App, Syp, and Rtn1. The 75 transcripts with the highest overall mean expression (i.e., the mean expression of all OVX and E 2 samples combined) in MeA Kiss1 cells are listed in Table 1.
Many transcripts for genes encoding neuropeptides, receptors, and proteins involved in neurotransmitter synthesis or transport were found to be highly expressed in the RNA-seq data. In addition to Gabbr1 and Gabbr2 being expressed in the RNA-seq dataset, which supports previous co-expression double-label assays with GABA B R and Kiss1 in the MeA (19), the RNAseq identified new transcripts that appear to be expressed by MeA Kiss1 cells, some of which are implicated in reproductive function or other behaviors. Aside from Kiss1, RNAseq identified other neuropeptides implicated in modulating reproduction, metabolism, and stress, such as Cck, Pdyn, Tac1, Gal, Cartpt, Trh, Sst, Npy, Agt, Vgf, and Penk (Table 2). Transcripts related to GABA synthesis and transport, like Vgat, Gad1, and Gad2, were also expressed, as were genes important for glutamate transport and synthesis, such as Glud1, Vglut1, and Vglut2 (Table 2). Table 2 provides a more detailed list of ligand (or their related enzymes) gene transcripts identified in this RNA-seq dataset. A few ligands of interest that were absent include Avp, C3, Calca, Nmu, Oxt, Ghrh, Gnrh1, and Tshb.
We examined transcripts for receptors in the RNA-seq data, which might provide insight into how MeA Kiss1 neurons are regulated, either by hormones or neural signaling. Sex steroid receptors were present, as might be predicted given estrogen's known ability to upregulate Kiss1 in the MeA. Specifically, estrogen receptor alpha (Esr1), androgen receptor (Ar), and progesterone receptor (Pgr) were each moderately expressed, while estrogen receptor beta (Esr2) had lower expression levels ( Table 3). Some additional receptors expressed in the RNA-seq data include Gabbr1, Gabbr2, Cnr1, Crhr1, Crhr2, Npy1r, Npy2r, Npy5r, Tacr1, and Thra. A more detailed list of the receptors expressed by MeA Kiss1 cells is in Table 3. Of note, the following gene transcripts are for some receptors that were absent (not expressed) in the current RNA-seq dataset: Chrna5, Ahrnb4, Pacapr1, Gpr50, Nmur2, P2rx2, Aplnr, Bdkrb2, Chrnb1, Lpar4, Npy4r, Rxfp2, Sctr, and Tbxa2r.

E 2 -mediated differential expression of MeA Kiss1 cell gene transcripts
Kiss1 gene expression in the MeA is known to be upregulated with E 2 treatment whereas Kiss1 levels in the MeA are very low or absent in gonadectomized mice lacking E 2 . Thus, we hypothesized that other gene transcripts in MeA Kiss1 cells might also change expression levels in the presence/absence of E 2 . Of the approximately 13,800 transcripts identified in the current RNA-seq data, only 45 genes had significantly different expression levels following 5-day E 2 treatment, in comparison to OVX mice lacking E 2 (Table 4; Figure 4). 10 of these genes were more highly expressed in OVX females (Figure 3, orange dots; Figure 4), while the remaining 35 genes were more highly expressed in E 2 -treated females (Figure 3, green dots; Figure 4). Table 4 provides a complete list of these differentially expressed genes, while Figure 4 is a heat map representing the expression patterns of each differentially expressed transcript for each of the 3 IP samples per treatment. As expected, Kiss1 is more highly expressed in E 2 -treated females in comparison to OVX females (Table 4, Figure 4). In addition to Kiss1, E 2 treatment upregulated transcripts encoding Galanin (Gal) and oxytocin receptor (Oxtr) as well as those for fibroblast growth factor receptor 1 (Fgfr1) and prokineticin receptor 2 (Prokr2), two genes whose loss of function results in Kallmann syndrome (72,73), and relaxin family peptide receptor 1 (Rxfp1), which is implicated in regulating sperm motility, pregnancy and birth (74, 75). Scg2 and Ecel1 (Endothelin Converting Enzyme Like 1) are important for neuropeptide release and are also higher in E 2 -treated than OVX females. E 2 treatment also resulted in greater expression in several genes linked to obesity and/or diabetes (insulin receptor substrate 2, Irs2 (76); Neuropeptide Y receptor 2, Npy2r (77); Transcription Elongation Regulator 1 Like, Tcerg1l (78)), nervous A glimma plot showing the expression of >13,800 gene transcripts identified by RNA-seq to be produced in MeA Kiss1 cells. Each dot represents a single gene transcript. Some transcripts were upregulated by E 2 treatment (green dots), whereas other transcripts were downregulated by E 2 treatment (orange dots). Grey dots represent gene transcripts whose expression did not differ between OVX and E 2 treatment groups. The y-axis represents the magnitude of difference in expression levels between OVX and E 2 samples, with greater positive or negative values indicating a greater difference in gene expression between the two treatment groups. Y-axis values close to 0 indicate little change in expression levels between OVX and E 2 treatments. The x-axis represents the average log-expression of each transcript, thus plotting the average of the mean expression levels of the OVX and E 2 treatment groups. Several genes of interest are marked by arrows to indicate their location in the glimma plot.   Inka2). Interestingly, most of the gene transcripts produced by MeA Kiss1 cells were not found to be significantly regulated by E 2 in the present study (Figure 3, grey dots; Figure 4).

Biological KEGG pathways represented by gene transcripts present in MeA Kiss1 cells
Biological pathway analysis, using KEGG pathways (68), was completed to identify potential functions of MeA Kiss1 cells. Pathway analysis examines how sets of gene transcripts cluster together to affect various biological processes. The top biological KEGG pathways with the lowest false discovery rate (pGFdr) were identified (Table 5). These pathways involve signaling pathways, such as the MAPK signaling pathway, calcium signaling pathway, and   neurotrophin signaling pathway. Other interesting pathways include pathways for amphetamine addiction and alcoholism, as well as those involved in regulating diseases, such as basal cell carcinoma and herpes simplex infection. There were no KEGG pathways identified that were significantly altered by E 2 treatment.

Validation of RNA-seq gene expression findings using double-label ISH/IHC
The RNA-seq data suggested that MeA Kiss1 cells express >13,800 gene transcripts, almost all of which have never been reported before for this specific Kiss1 population. Therefore, in addition to validating our immuno-pulldown procedure via RT-PCR (Figure 2), we performed double-label ISH/IHC on brains from female Kiss tdTom mice to confirm the co-expression of a gene identified in the RNA-seq dataset that was not previously known to be present in MeA Kiss1 neurons (Cartpt; Figure 5, Table 2). We found that Cartpt mRNA was highly expressed in MeA brain slices, including very high overlap with MeA kisspeptin cells (tdTomato fluorescence; Figure 5).

Discussion
The MeA has been implicated in modulating numerous physiological and behavioral processes (32-44), including aspects of reproductive physiology. In females, lesions to the MeA disrupt ovarian cycles and prevent the E 2 -mediated LH surge, whereas acute electrical stimulation of the MeA stimulate LH release (34-37). However, the exact cellular and molecular mechanisms and specific cell types for these MeA effects on reproduction remain unknown. Kisspeptin is able to stimulate the reproductive axis, via direct action on GnRH neurons, but very little is known about the functions, reproductive or otherwise, of Kiss1 neurons in the MeA, or how these neurons are regulated. Recent studies have begun to examine the detailed molecular profiles of AVPV and ARC Kiss1 neurons, including identifying numerous gene transcripts in these populations that are regulated by E 2 (25-27), Transcripts are listed in descending order based on total mean expression (average log CPM for E 2 and GDX mice). p <0.05 signifies significant differential expression between E 2 and OVX conditions, denoted in bold type. In this study, we used the Ribotag technique to selectively isolate actively-translated mRNAs in MeA Kiss1 neurons. We first validated the success and specificity of this technique using RT-PCR. First, we found that both the Input (all cells in the MeA micropunch) and IP samples (only mRNAs from Kiss1 cells) were positive for Kiss1 mRNA from E 2 -treated females, as we would predict. Likewise, in OVX samples, Kiss1 was detected in both Input and IP samples as expected, with lower mRNA levels (lighter bands) than for E 2 samples given that MeA Kiss1 expression is known to be stimulated by E 2 . The specificity of our technique was further confirmed by our samples from Cre-controls in which Kiss1 was present in the Input samples, as expected because they contain mRNA from all MeA cells in the micropunch including Kiss1 cells, but not in the IP samples. The lack of Kiss1 (and Cck) in the Cre-IP samples indicates that our Ribotag immuno-pulldown technique was specific to isolating mRNA from just Kiss1 cells in the micropunch. Our initial histological identification of Cck co-expression with Kiss1-TdTomato cells in the MeA (Figure 1) enabled us to use Cck as a positive control for MeA kisspeptin neurons in the RT-PCR validation step. Indeed, similar to Kiss1, we found that Cck mRNA was present in both Input and IP samples in E 2 -treated and OVX mice, with no Cck expression found in the IP samples of Cre-control mice, further supporting the specificity of the Ribotag technique. Knowing that Avp and Kiss1 are present in different regions of the MeA, we examined our Input and IP samples for the presence/absence of Avp as a negative control for kisspeptin cell transcripts. Whereas Input samples contained Avp mRNA, indicating its presence in the MeA region, no Avp was detected in any IP samples as we would predict if Avp is not expressed in Kiss1 cells. We also confirmed the co-    MAPK signaling pathway, and calcium signaling pathway, which are important processes for all cells, not just Kiss1 cells. Of interest, 2 of the 5 addiction pathways, amphetamine addiction and alcoholism, were represented in the top KEGG pathways, potentially because of the heavy involvement of dopamine, dynorphin, and glutamate in these pathways. Interestingly, many of the KEGG pathways identified for MeA Kiss1 cells were the same pathways identified for AVPV Kiss1 cells (25), which may suggest that some of these functions may be generalized to many cells in the brain and/or suggestive of undiscovered shared functions of Kiss1 neurons from several brain areas. We did not find any KEGG pathways that contained transcripts that were more represented in GDX or E 2 -treated mice. This result is not surprising given the low number of differentially expressed genes identified between the two hormone treatment groups, as discussed more below.
The mean expression analysis of gene transcripts in our dataset revealed >13,800 gene transcripts produced by MeA Kiss1 neurons, of which almost all are newly identified for this specific cell population. We identified genes for several neuropeptides found in other Kiss1 neuronal populations, such as Vgf, Penk, Tac1, Pdyn, Pnoc, and Gal, as well as several sex steroid receptors, Ar, Pgr, Esr1, and Esr2, though Esr2 levels were much lower than the other sex steroid receptors. The actions and targets of the identified neuropeptides co-expressed in MeA Kiss1 neurons remain unknown but may give clues into understanding the potential functions of these neurons. Interestingly, many of the identified neuropeptide genes, including the ones listed above, are also expressed in AVPV Kiss1 neurons (25,27), and Kiss1 is upregulated by E 2 in both the MeA and AVPV (13-15, 17, 18), which may indicate AVPV and MeA Kiss1 cells share some similar functions. The presence of Ar in MeA Kiss1 neurons is interesting because previous studies demonstrated that T or E 2 , but not DHT, can increase MeA Kiss1 expression (17). However, the MeA is a sexually dimorphic brain structure, both in anatomy and genes expressed (82)(83)(84), and this is hypothesized to be produced in part by androgen signaling. Thus, it may not be surprising that MeA Kiss1 cells (which have some sexually dimorphic qualities), express Ar, even if Kiss1 expression itself is not altered by AR signaling. Interestingly, we found that Cyp19a1, the gene encoding for aromatase, is not expressed in MeA Kiss1 cells, based on our current RNA-seq dataset, suggesting that the conversion of androgens to estrogens is occurring elsewhere. It is also possible that other gene transcripts in MeA Kiss1  (19,20). Our present RNA-seq data support this finding as transcripts for both Gababr1 and Gababr2 subunits were highly expressed. Importantly, the expression of receptor transcripts by MeA Kiss1 cells does not automatically mean that the Kiss1 gene itself is regulated via activation of these receptors, as the receptor signaling may influence other genes or processes in these cells. Surprisingly, only 45 of the 13,000+ gene transcripts produced by MeA Kiss1 neurons were found to be significantly regulated by E 2 status, with a majority of these differentially expressed transcripts increasing with E 2 treatment. This contrasts with AVPV Kiss1 cells in which 683 transcripts were altered by similar E 2 treatment (25). Regardless, it is possible that some of these E 2 -senesitive transcripts in MeA Kiss1 neurons are involved in E 2 -regulated behavior and physiology, and future research can examine the roles of transcripts like Gal and Oxtr specifically in MeA Kiss1 neurons in such processes. It is unknown at present if E 2 's effects on these 45 differentially expressed genes are due to E 2 action directly or indirectly on MeA kisspeptin cells or via which ER subtype the effects are mediated. Supporting the possibility of some direct E 2 regulation, our RNA-seq dataset indicated that MeA Kiss1 cells express Esr1 and, to a lesser degree, Esr2. However, some of the differential expression of gene transcripts may be due to indirect actions of E 2 on other afferent cells that communicate with MeA kisspeptin neurons. Because so many transcripts were not differentially expressed with E 2 treatment, it is possible that these MeA kisspeptin neurons have primary functions that are not dependent on E 2 status. Future research might focus on the genes that were most highly expressed in these neurons (Tables 1-3) regardless of hormone status.
Recent studies by our group and Göcz and colleagues (2022) and have begun to examine the gene transcripts made by kisspeptin neurons in the murine AVPV (25) and ARC (26), respectively, and evaluated the role of E 2 in regulating these gene transcripts (see Table 6 for a summary and comparison of these studies). In ARC kisspeptin neurons, 2,329 genes were found to be regulated by E 2 , with about an equal number of genes being upregulated vs downregulated by E 2 (26). When these samples were analyzed differently with low-expressing genes removed, there were still 1,583 genes expressed by ARC kisspeptin neurons that were significantly altered by E 2 treatment (27). In contrast to the ARC, kisspeptin neurons in the AVPV of the same mice only had 222 genes that were E 2 regulated, with 142 of those genes being upregulated by E 2 (27). Thus, more gene transcripts produced by ARC kisspeptin neurons appear to be responsive to E 2 than in AVPV kisspeptin neurons. Our group previously examined gene transcripts made by female AVPV kisspeptin neurons using a different mouse model, E 2 treatment paradigm, and RNA isolation techniques than used by Göcz and colleagues. In that prior study, we found a higher number of gene transcripts in AVPV kisspeptin neurons, 683 transcripts, that were regulated by E 2 (25). It is possible the differences in isolation techniques (isolation from the ribosomes vs isolation from the entire kisspeptin neuron), E 2 hormone treatment paradigms, or the parameters used to exclude low-expressing genes (summarized in Table 6) resulted in the different number of E 2 -regualted transcripts between these two AVPV studies. Despite these methodological differences, in both studies of AVPV kisspeptin neurons, more gene transcripts were upregulated by E 2 than downregulated by E 2 , and over 50 of these gene transcripts that were E 2 regulated were identified in both studies (25,27). In the current study, we used the same mice and methodology as our prior AVPV study (25), and found that only 45 of >13,800 gene transcripts of MeA kisspeptin neurons were regulated by E 2 , with Kiss1 being one of them. Thus, our present data suggest that, in comparison to AVPV and ARC cells, most gene transcripts produced by MeA kisspeptin neurons are not sensitive to E 2 and may be more strongly regulated by other signaling factors. Additional research is still needed to understand how MeA kisspeptin neurons are regulated, and the present RNA-seq dataset provides a starting point for such future research.
The current RNA-seq data greatly expands upon the limited information regarding what signaling factors and receptors are produced by MeA Kiss1 neurons. We wanted to further validate our findings by performing an IHC/ISH for on a gene (Cartpt) identified for the first time by our RNAseq to be present in MeA kisspeptin neurons. This histological assay determined that virtually all of kisspeptin cells in the MeA region express Cartpt, supporting the very high expression levels of Cartpt in the RNA-seq dataset. At present, it is technically difficult to detect sufficient Kiss1 mRNA in the MeA using RNAscope methodology owing to lower levels of Kiss1 mRNA expressed per cell in this brain region than in the hypothalamus. One limitation in the current study is that the identification of MeA kisspeptin cells was achieved using a fluorescent reporter (TdTomato), instead of looking at actual Kiss1 mRNA expression. When more sensitive methodology becomes available, future research should examine MeA Kiss1 co-expression with other transcripts of interest identified in this RNA-seq dataset, to both confirm and better understand the co-expression of neuropeptides/receptors/enzymes and active Kiss1 expression.
In summary, the current dataset is the first large-scale examination of the identities of neuropeptides, receptors, and neurotransmitter-related genes produced by MeA Kiss1 cells. Other than being regulated by E 2 and GABA, how Kiss1 in the MeA is regulated remains completely unknown. The current study identified additional receptor transcripts made by these kisspeptin neurons, which could drive future research directions on the regulation of these neurons and the Kiss1 gene itself. We also identified many neuropeptides and signaling factors not previously known to be produced by MeA Kiss1 neurons, some of which are implicated in MeA-influenced processes like reproduction and metabolism, along with transcripts important for basic cell maintenance and survival. Interestingly, while > 13,800 transcripts are made in these cells, only 45 transcripts were significantly regulated by E 2 . Whether and how this focused E 2 regulation is related to these neurons' functional roles remains to be determined. Very little is known about Kiss1 cells in the MeA and the current dataset therefore provides a valuable starting point for future research to examine the regulation and function of Kiss1 neurons in the MeA.

Data availability statement
The data presented in the study are deposited in the Gene Expression Omnibus data repository, accession number GSE224788.

Ethics statement
The animal study was reviewed and approved by IACUC committees at Albany Medical College and the University of California, San Diego.