A methylation clock model of mild SARS-CoV-2 infection provides insight into immune dysregulation.

DNA methylation comprises a cumulative record of lifetime exposures superimposed on genetically determined markers. Little is known about methylation dynamics in humans following an acute perturbation, such as infection. We characterized the temporal trajectory of blood epigenetic remodeling in 133 participants in a prospective study of young adults before, during, and after asymptomatic and mildly symptomatic SARS-CoV-2 infection. The differential methylation caused by asymptomatic or mildly symptomatic infections was indistinguishable. While differential gene expression largely returned to baseline levels after the virus became undetectable, some differentially methylated sites persisted for months of follow-up, with a pattern resembling autoimmune or inflammatory disease. We leveraged these responses to construct methylation-based machine learning models that distinguished samples from pre-, during-, and postinfection time periods, and quantitatively predicted the time since infection. The clinical trajectory in the young adults and in a diverse cohort with more severe outcomes was predicted by the similarity of methylation before or early after SARS-CoV-2 infection to the model-defined postinfection state. Unlike the phenomenon of trained immunity, the postacute SARS-CoV-2 epigenetic landscape we identify is antiprotective.

SARS-CoV-2 Outbreak Dynamics in an Isolated US Military Recruit Training Center With Rigorous Prevention Measures.

BACKGROUND: Marine recruits training at Parris Island experienced an unexpectedly high rate of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, despite preventive measures including a supervised, 2-week, pre-entry quarantine. We characterize SARS-CoV-2 transmission in this cohort. METHODS: Between May and November 2020, we monitored 2,469 unvaccinated, mostly male, Marine recruits prospectively during basic training. If participants tested negative for SARS-CoV-2 by quantitative polymerase chain reaction (qPCR) at the end of quarantine, they were transferred to the training site in segregated companies and underwent biweekly testing for 6 weeks. We assessed the effects of coronavirus disease 2019 (COVID-19) prevention measures on other respiratory infections with passive surveillance data, performed phylogenetic analysis, and modeled transmission dynamics and testing regimens. RESULTS: Preventive measures were associated with drastically lower rates of other respiratory illnesses. However, among the trainees, 1,107 (44.8%) tested SARS-CoV-2-positive, with either mild or no symptoms. Phylogenetic analysis of viral genomes from 580 participants revealed that all cases but one were linked to five independent introductions, each characterized by accumulation of mutations across and within companies, and similar viral isolates in individuals from the same company. Variation in company transmission rates (mean reproduction number R 0 ; 5.5 [95% confidence interval [CI], 5.0, 6.1]) could be accounted for by multiple initial cases within a company and superspreader events. Simulations indicate that frequent rapid-report testing with case isolation may minimize outbreaks. CONCLUSIONS: Transmission of wild-type SARS-CoV-2 among Marine recruits was approximately twice that seen in the community. Insights from SARS-CoV-2 outbreak dynamics and mutations spread in a remote, congregate setting may inform effective mitigation strategies.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle.

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.

Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity.

Immediate-early response genes (IEGs) are rapidly and transiently induced following an extracellular signal. Elucidating the IEG response patterns in single cells (SCs) requires assaying large numbers of timed samples at high accuracy while minimizing handling effects. To achieve this, we developed and validated RNA stabilization Buffer for Examination of Single-cell Transcriptomes (RNA-Best), a versatile single-step cell and tissue preservation protocol that stabilizes RNA in intact SCs without perturbing transcription patterns. We characterize for the first time SC heterogeneity in IEG responses to pulsatile gonadotropin-releasing hormone (GnRH) stimuli in pituitary gonadotrope cells. Our study identifies a gene-specific hierarchical pattern of all-or-none transcript induction elicited by increasing concentrations of GnRH. This quantal pattern of gene activation raises the possibility that IEG activation, when accurately resolved at the SC level, may be mediated by gene bits that behave as pure binary switches.

Modeling and high-throughput experimental data uncover the mechanisms underlying Fshb gene sensitivity to gonadotropin-releasing hormone pulse frequency.

Neuroendocrine control of reproduction by brain-secreted pulses of gonadotropin-releasing hormone (GnRH) represents a longstanding puzzle about extracellular signal decoding mechanisms. GnRH regulates the pituitary gonadotropin's follicle-stimulating hormone (FSH) and luteinizing hormone (LH), both of which are heterodimers specified by unique ß subunits (FSHß/LHß). Contrary to Lhb, Fshb gene induction has a preference for low-frequency GnRH pulses. To clarify the underlying regulatory mechanisms, we developed three biologically anchored mathematical models: 1) parallel activation of Fshb inhibitory factors (e.g. inhibin α and VGF nerve growth factor-inducible), 2) activation of a signaling component with a refractory period (e.g. G protein), and 3) inactivation of a factor needed for Fshb induction (e.g. growth differentiation factor 9). Simulations with all three models recapitulated the Fshb expression levels obtained in pituitary gonadotrope cells perifused with varying GnRH pulse frequencies. Notably, simulations altering average concentration, pulse duration, and pulse frequency revealed that the apparent frequency-dependent pattern of Fshb expression in model 1 actually resulted from variations in average GnRH concentration. In contrast, models 2 and 3 showed "true" pulse frequency sensing. To resolve which components of this GnRH signal induce Fshb, we developed a high-throughput parallel experimental system. We analyzed over 4,000 samples in experiments with varying near-physiological GnRH concentrations and pulse patterns. Whereas Egr1 and Fos genes responded only to variations in average GnRH concentration, Fshb levels were sensitive to both average concentration and true pulse frequency. These results provide a foundation for understanding the role of multiple regulatory factors in modulating Fshb gene activity.

Characterization of Gonadotrope Secretoproteome Identifies Neurosecretory Protein VGF-derived Peptide Suppression of Follicle-stimulating Hormone Gene Expression.

Reproductive function is controlled by the pulsatile release of hypothalamic gonadotropin-releasing hormone (GnRH), which regulates the expression of the gonadotropins luteinizing hormone and FSH in pituitary gonadotropes. Paradoxically, Fshb gene expression is maximally induced at lower frequency GnRH pulses, which provide a very low average concentration of GnRH stimulation. We studied the role of secreted factors in modulating gonadotropin gene expression. Inhibition of secretion specifically disrupted gonadotropin subunit gene regulation but left early gene induction intact. We characterized the gonadotrope secretoproteome and global mRNA expression at baseline and after Gαs knockdown, which has been found to increase Fshb gene expression (1). We identified 1077 secreted proteins or peptides, 19 of which showed mRNA regulation by GnRH or/and Gαs knockdown. Among several novel secreted factors implicated in Fshb gene regulation, we focused on the neurosecretory protein VGF. Vgf mRNA, whose gene has been implicated in fertility (2), exhibited high induction by GnRH and depended on Gαs In contrast with Fshb induction, Vgf induction occurred preferentially at high GnRH pulse frequency. We hypothesized that a VGF-derived peptide might regulate Fshb gene induction. siRNA knockdown or extracellular immunoneutralization of VGF augmented Fshb mRNA induction by GnRH. GnRH stimulated the secretion of the VGF-derived peptide NERP1. NERP1 caused a concentration-dependent decrease in Fshb gene induction. These findings implicate a VGF-derived peptide in selective regulation of the Fshb gene. Our results support the concept that signaling specificity from the cell membrane GnRH receptor to the nuclear Fshb gene involves integration of intracellular signaling and exosignaling regulatory motifs.

Growth differentiation factor 9 (GDF9) forms an incoherent feed-forward loop modulating follicle-stimulating hormone ß-subunit (FSHß) gene expression.

Gonadotropin-releasing hormone (GnRH) is secreted in brief pulses from the hypothalamus and regulates follicle-stimulating hormone ß-subunit (FSHß) gene expression in pituitary gonadotropes in a frequency-sensitive manner. The mechanisms underlying its preferential and paradoxical induction of FSHß by low frequency GnRH pulses are incompletely understood. Here, we identify growth differentiation factor 9 (GDF9) as a GnRH-suppressed autocrine inducer of FSHß gene expression. GDF9 gene transcription and expression were preferentially decreased by high frequency GnRH pulses. GnRH regulation of GDF9 was concentration-dependent and involved ERK and PKA. GDF9 knockdown or immunoneutralization reduced FSHß mRNA expression. Conversely, exogenous GDF9 induced FSHß expression in immortalized gonadotropes and in mouse primary pituitary cells. GDF9 exposure increased FSH secretion in rat primary pituitary cells. GDF9 induced Smad2/3 phosphorylation, which was impeded by ALK5 knockdown and by activin receptor-like kinase (ALK) receptor inhibitor SB-505124, which also suppressed FSHß expression. Smad2/3 knockdown indicated that FSHß induction by GDF9 involved Smad2 and Smad3. FSHß mRNA induction by GDF9 and GnRH was synergistic. We hypothesized that GDF9 contributes to a regulatory loop that tunes the GnRH frequency-response characteristics of the FSHß gene. To test this, we determined the effects of GDF9 knockdown on FSHß induction at different GnRH pulse frequencies using a parallel perifusion system. Reduction of GDF9 shifted the characteristic pattern of GnRH pulse frequency sensitivity. These results identify GDF9 as contributing to an incoherent feed-forward loop, comprising both intracellular and secreted components, that regulates FSHß expression in response to activation of cell surface GnRH receptors.

Low-variance RNAs identify Parkinson's disease molecular signature in blood.

The diagnosis of Parkinson's disease (PD) is usually not established until advanced neurodegeneration leads to clinically detectable symptoms. Previous blood PD transcriptome studies show low concordance, possibly resulting from the use of microarray technology, which has high measurement variation. The Leucine-rich repeat kinase 2 (LRRK2) G2019S mutation predisposes to PD. Using preclinical and clinical studies, we sought to develop a novel statistically motivated transcriptomic-based approach to identify a molecular signature in the blood of Ashkenazi Jewish PD patients, including LRRK2 mutation carriers. Using a digital gene expression platform to quantify 175 messenger RNA (mRNA) markers with low coefficients of variation (CV), we first compared whole-blood transcript levels in mouse models (1) overexpressing wild-type (WT) LRRK2, (2) overexpressing G2019S LRRK2, (3) lacking LRRK2 (knockout), and (4) and in WT controls. We then studied an Ashkenazi Jewish cohort of 34 symptomatic PD patients (both WT LRRK2 and G2019S LRRK2) and 32 asymptomatic controls. The expression profiles distinguished the four mouse groups with different genetic background. In patients, we detected significant differences in blood transcript levels both between individuals differing in LRRK2 genotype and between PD patients and controls. Discriminatory PD markers included genes associated with innate and adaptive immunity and inflammatory disease. Notably, gene expression patterns in levodopa-treated PD patients were significantly closer to those of healthy controls in a dose-dependent manner. We identify whole-blood mRNA signatures correlating with LRRK2 genotype and with PD disease state. This approach may provide insight into pathogenesis and a route to early disease detection.

Protocol for high-throughput DNA methylation profiling in rat tissues using automated reduced representation bisulfite sequencing.

Although reduced representation bisulfite sequencing (RRBS) measures DNA methylation (DNAme) across CpG-rich genomic regions with high sensitivity, the assay can be time-consuming and prone to batch effects. Here, we present a high-throughput, automated RRBS protocol starting with DNA extraction from frozen rat tissues. We describe steps for RRBS library preparation, library quality control, and sequencing. We also detail an optimized pipeline for sequencing data processing. This protocol has been applied successfully to DNAme profiling across multiple rat tissues. For complete details on the use and execution of this protocol, please refer to Nair et al.1.

Molecular adaptations in response to exercise training are associated with tissue-specific transcriptomic and epigenomic signatures.

Regular exercise has many physical and brain health benefits, yet the molecular mechanisms mediating exercise effects across tissues remain poorly understood. Here we analyzed 400 high-quality DNA methylation, ATAC-seq, and RNA-seq datasets from eight tissues from control and endurance exercise-trained (EET) rats. Integration of baseline datasets mapped the gene location dependence of epigenetic control features and identified differing regulatory landscapes in each tissue. The transcriptional responses to 8 weeks of EET showed little overlap across tissues and predominantly comprised tissue-type enriched genes. We identified sex differences in the transcriptomic and epigenomic changes induced by EET. However, the sex-biased gene responses were linked to shared signaling pathways. We found that many G protein-coupled receptor-encoding genes are regulated by EET, suggesting a role for these receptors in mediating the molecular adaptations to training across tissues. Our findings provide new insights into the mechanisms underlying EET-induced health benefits across organs.

Multi-omic identification of key transcriptional regulatory programs during endurance exercise training.

Transcription factors (TFs) play a key role in regulating gene expression and responses to stimuli. We conducted an integrated analysis of chromatin accessibility, DNA methylation, and RNA expression across eight rat tissues following endurance exercise training (EET) to map epigenomic changes to transcriptional changes and determine key TFs involved. We uncovered tissue-specific changes and TF motif enrichment across all omic layers, differentially accessible regions (DARs), differentially methylated regions (DMRs), and differentially expressed genes (DEGs). We discovered distinct routes of EET-induced regulation through either epigenomic alterations providing better access for TFs to affect target genes, or via changes in TF expression or activity enabling target gene response. We identified TF motifs enriched among correlated epigenomic and transcriptomic alterations, DEGs correlated with exercise-related phenotypic changes, and EET-induced activity changes of TFs enriched for DEGs among their gene targets. This analysis elucidates the unique transcriptional regulatory mechanisms mediating diverse organ effects of EET.

Association of survival and disease progression with chromosomal instability: a genomic exploration of colorectal cancer.

During disease progression the cells that comprise solid malignancies undergo significant changes in gene copy number and chromosome structure. Colorectal cancer provides an excellent model to study this process. To indentify and characterize chromosomal abnormalities in colorectal cancer, we performed a statistical analysis of 299 expression and 130 SNP arrays profiled at different stages of the disease, including normal tissue, adenoma, stages 1-4 adenocarcinoma, and metastasis. We identified broad (> 1/2 chromosomal arm) and focal (< 1/2 chromosomal arm) events. Broad amplifications were noted on chromosomes 7, 8q, 13q, 20, and X and broad deletions on chromosomes 4, 8p, 14q, 15q, 17p, 18, 20p, and 22q. Focal events (gains or losses) were identified in regions containing known cancer pathway genes, such as VEGFA, MYC, MET, FGF6, FGF23, LYN, MMP9, MYBL2, AURKA, UBE2C, and PTEN. Other focal events encompassed potential new candidate tumor suppressors (losses) and oncogenes (gains), including CCDC68, CSMD1, POLR1D, and PMEPA1. From the expression data, we identified genes whose expression levels reflected their copy number changes and used this relationship to impute copy number changes to samples without accompanying SNP data. This analysis provided the statistical power to show that deletions of 8p, 4p, and 15q are associated with survival and disease progression, and that samples with simultaneous deletions in 18q, 8p, 4p, and 15q have a particularly poor prognosis. Annotation analysis reveals that the oxidative phosphorylation pathway shows a strong tendency for decreased expression in the samples characterized by poor prognosis.

Optimization of the Omni-ATAC protocol to chromatin accessibility profiling in snap-frozen rat adipose and muscle tissues.

ATAC-seq is a fast and sensitive method for the epigenomic profiling of open chromatin and for mapping of transcription factor binding sites [1]. Despite the development of the Omni-ATAC protocol for the profiling of chromatin accessibility in frozen tissues [2], studies in adipose tissue have been restricted due to technical challenges including the high lipid content of adipocytes and reproducibility issues between replicates. Here, we provide a modified Omni-ATAC protocol that achieves high data reproducibility in various tissue types from rat, including adipose and muscle tissues [3].•This protocol describes a methodology that enables chromatin accessibility profiling from snap-frozen rat adipose and muscle tissues.•The technique comprises an optimized bead-based tissue homogenization process that substitutes to Dounce homogenization, reduces variability in the experimental procedure, and is adaptable to various tissue types.•In comparison with the Omni-ATAC protocol, the method described here results in improved ATAC-seq data quality that complies with ENCODE quality standards.

Multi-omics profiling of single nuclei from frozen archived postmortem human pituitary tissue.

Concomitant profiling of transcriptome and chromatin accessibility in isolated nuclei can reveal gene regulatory control mechanisms in health and disease. We report a single nucleus multi-omics analysis protocol optimized for frozen archived postmortem human pituitaries that is also effective for frozen ovine and murine pituitaries and human skeletal muscle biopsies. Its main advantages are that (1) it is not limited to fresh tissue, (2) it avoids tissue dissociation-induced transcriptional changes, and (3) it includes a novel, automated quality control pipeline. For complete details on the use and execution of this protocol, please refer to Ruf-Zamojski et al. (2021) and Zhang et al. (2022).

Single nucleus transcriptome and chromatin accessibility of postmortem human pituitaries reveal diverse stem cell regulatory mechanisms.

Despite their importance in tissue homeostasis and renewal, human pituitary stem cells (PSCs) are incompletely characterized. We describe a human single nucleus RNA-seq and ATAC-seq resource from pediatric, adult, and aged postmortem pituitaries (snpituitaryatlas.princeton.edu) and characterize cell-type-specific gene expression and chromatin accessibility programs for all major pituitary cell lineages. We identify uncommitted PSCs, committing progenitor cells, and sex differences. Pseudotime trajectory analysis indicates that early-life PSCs are distinct from the other age groups. Linear modeling of same-cell multiome data identifies regulatory domain accessibility sites and transcription factors that are significantly associated with gene expression in PSCs compared with other cell types and within PSCs. We identify distinct deterministic mechanisms that contribute to heterogeneous marker expression within PSCs. These findings characterize human stem cell lineages and reveal diverse mechanisms regulating key PSC genes and cell type identity.

Single nucleus multi-omics regulatory landscape of the murine pituitary.

To provide a multi-omics resource and investigate transcriptional regulatory mechanisms, we profile the transcriptome, chromatin accessibility, and methylation status of over 70,000 single nuclei (sn) from adult mouse pituitaries. Paired snRNAseq and snATACseq datasets from individual animals highlight a continuum between developmental epigenetically-encoded cell types and transcriptionally-determined transient cell states. Co-accessibility analysis-based identification of a putative Fshb cis-regulatory domain that overlaps the fertility-linked rs11031006 human polymorphism, followed by experimental validation illustrate the use of this resource for hypothesis generation. We also identify transcriptional and chromatin accessibility programs distinguishing each major cell type. Regulons, which are co-regulated gene sets sharing binding sites for a common transcription factor driver, recapitulate cell type clustering. We identify both cell type-specific and sex-specific regulons that are highly correlated with promoter accessibility, but not with methylation state, supporting the centrality of chromatin accessibility in shaping cell-defining transcriptional programs. The sn multi-omics atlas is accessible at snpituitaryatlas.princeton.edu.

Sedentary and Trained Older Men Have Distinct Circulating Exosomal microRNA Profiles at Baseline and in Response to Acute Exercise.

Exercise has multi-systemic benefits and attenuates the physiological impairments associated with aging. Emerging evidence suggests that circulating exosomes mediate some of the beneficial effects of exercise via the transfer of microRNAs between tissues. However, the impact of regular exercise and acute exercise on circulating exosomal microRNAs (exomiRs) in older populations remains unknown. In the present study, we analyzed circulating exomiR expression in endurance-trained elderly men (n = 5) and age-matched sedentary males (n = 5) at baseline (Pre), immediately after a forty minute bout of aerobic exercise on a cycle ergometer (Post), and three hours after this acute exercise (3hPost). Following the isolation and enrichment of exosomes from plasma, exosome-enriched preparations were characterized and exomiR levels were determined by sequencing. The effect of regular exercise on circulating exomiRs was assessed by comparing the baseline expression levels in the trained and sedentary groups. The effect of acute exercise was determined by comparing baseline and post-training expression levels in each group. Regular exercise resulted in significantly increased baseline expression of three exomiRs (miR-486-5p, miR-215-5p, miR-941) and decreased expression of one exomiR (miR-151b). Acute exercise altered circulating exomiR expression in both groups. However, exomiRs regulated by acute exercise in the trained group (7 miRNAs at Post and 8 at 3hPost) were distinct from those in the sedentary group (9 at Post and 4 at 3hPost). Pathway analysis prediction and reported target validation experiments revealed that the majority of exercise-regulated exomiRs are targeting genes that are related to IGF-1 signaling, a pathway involved in exercise-induced muscle and cardiac hypertrophy. The immediately post-acute exercise exomiR signature in the trained group correlates with activation of IGF-1 signaling, whereas in the sedentary group it is associated with inhibition of IGF-1 signaling. While further validation is needed, including measurements of IGF-1/IGF-1 signaling in blood or skeletal muscle, our results suggest that training status may counteract age-related anabolic resistance by modulating circulating exomiR profiles both at baseline and in response to acute exercise.

MDM2 gene amplification is correlated to tumor progression but not to the presence of SNP309 or TP53 mutational status in primary colorectal cancers.

Mdm2 is the main regulator of p53 and is amplified in approximately 7% of all human cancers. MDM2 gene amplification as well as expression has been correlated to an increased tumorigenic potential. We have analyzed the prevalence of MDM2 gene amplifications and SNP309 in 284 colorectal tumors using a relatively new highly sensitive PCR/ligase detection reaction method in relation to TP53 mutational status and genomic instability. We found MDM2 to be amplified in 9% of the 284 colorectal cancers analyzed and a significantly higher proportion of tumors with high MDM2 gene amplification retained a wild-type p53 gene (P = 0.058). MDM2 gene amplification was significantly correlated to advanced tumor stage. Several small-molecule MDM2 antagonists have already been identified that either physically inhibit the p53-MDM2 binding or the E3 ligase function of MDM2. Our results suggest that MDM2 is a promising target for this type of cancer therapy in a substantial subgroup of colorectal cancers.

CpG island methylator phenotype associates with low-degree chromosomal abnormalities in colorectal cancer.

PURPOSE: Aberrant promoter methylation and genomic instability occur frequently during colorectal cancer development. CpG island methylator phenotype (CIMP) has been shown to associate with microsatellite instability, and BRAF mutation and is often found in the right-side colon. Nevertheless, the relative importance of CIMP and chromosomal instability (CIN) for tumorigenesis has yet to be thoroughly investigated in sporadic colorectal cancers. EXPERIMENTAL DESIGN: We determined CIMP in 161 primary colorectal cancers and 66 matched normal mucosae using a quantitative bisulfite/PCR/ligase detection reaction (LDR)/Universal Array assay. The validity of CIMP was confirmed in a subset of 60 primary tumors using MethyLight assay and five independent markers. In parallel, CIN was analyzed in the same study cohort using Affymetrix 50K Human Mapping arrays. RESULTS: The identified CIMP-positive cancers correlate with microsatellite instability (P = 0.075) and the BRAF mutation V600E (P = 0.00005). The array-based high-resolution analysis of chromosomal aberrations indicated that the degree of aneuploidy is spread over a wide spectrum among analyzed colorectal cancers. Whether CIN was defined by copy number variations in selected microsatellite loci (criterion 1) or considered as a continuous variable (criterion 2), CIMP-positive samples showed a strong correlation with low-degree chromosomal aberrations (P = 0.075 and P = 0.012, respectively). Similar correlations were observed when CIMP was determined by MethyLight assay (P = 0.001 and P = 0.013, respectively). CONCLUSION: CIMP-positive tumors generally possess lower chromosomal aberrations, which may only be revealed using a genome-wide approach. The significant difference in the degree of chromosomal aberrations between CIMP-positive and the remainder of samples suggests that epigenetic (CIMP) and genetic (CIN) abnormalities may arise from independent molecular mechanisms of tumor progression.

Cytogenetic, Genomic, and Functional Characterization of Pituitary Gonadotrope Cell Lines.

LßT2 and αT3-1 are important, widely studied cell line models for the pituitary gonadotropes that were generated by targeted tumorigenesis in transgenic mice. LßT2 cells are more mature gonadotrope precursors than αT3-1 cells. Microsatellite authentication patterns, chromosomal characteristics, and their intercellular variation have not been reported. We performed microsatellite and cytogenetic analysis of both cell types at early passage numbers. Short tandem repeat (STR) profiling was consistent with a mixed C57BL/6J × BALB/cJ genetic background, with distinct patterns for each cell type. Spectral karyotyping in αT3-1 cells revealed cell-to-cell variation in chromosome composition and pseudodiploidy. In LßT2 cells, chromosome counting and karyotyping demonstrated pseudotriploidy and high chromosomal variation among cells. Chromosome copy number variation was confirmed by single-cell DNA sequencing. Chromosomal compositions were consistent with a male sex for αT3-1 and a female sex for LßT2 cells. Among LßT2 stocks used in multiple laboratories, we detected two genetically similar but distinguishable lines via STR authentication, LßT2a and LßT2b. The two lines differed in morphological appearance, with LßT2a having significantly smaller cell and nucleus areas. Analysis of immediate early gene and gonadotropin subunit gene expression revealed variations in basal expression and responses to continuous and pulsatile GnRH stimulation. LßT2a showed higher basal levels of Egr1, Fos, and Lhb but lower Fos induction. Fshb induction reached significance only in LßT2b cells. Our study highlights the heterogeneity in gonadotrope cell line genomes and provides reference STR authentication patterns that can be monitored to improve experimental reproducibility and facilitate comparisons of results within and across laboratories.