Genome-Wide Association Study on Three Behaviors Tested in an Open Field in Heterogeneous Stock Rats Identifies Multiple Loci Implicated in Psychiatric Disorders

Inbred rat heredity and sex affect oral oxycodone self-administration and augmented intake in long sessions: correlations with anxiety and novelty-seeking

Oxycodone abuse begins with prescription oral oxycodone, yet vulnerability factors determining abuse are largely undefined. We evaluated genetic vulnerability in a rat model of oral oxycodone self-administration (SA): increasing oxycodone concentration/session (0.025-0.1mg/ml; 1,4,16-h) followed by extinction and reinstatement. Active licks and oxycodone intake were greater in females than males during 4-h and 16-h sessions (p< 0.001). Each sex increased intake during 16-h vs 4-h sessions (p<2e-16), but a subset of strains dramatically augmented intake at 16-h (p=0.0005). Heritability (h 2) of active licks/4-h at increasing oxycodone dose ranged from 0.30-0.53. Under a progressive ratio schedule, breakpoints were strain-dependent (p<2e-16). Cued reinstatement was greater in females (p<0.001). Naive rats were assessed by elevated plus maze (EPM), open field (OF) and novel object interaction (NOI). We correlated these behaviors with 28 parameters of oxycodone SA. Anxiety-defining EPM traits were most associated with SA in both sexes, whereas more OF and NOI traits were SA-associated in males. Sex and heredity are major determinants of motivation to take and seek oxycodone; intake augments dramatically during extended access in specific strains; and pleiotropic genes affect anxiety and multiple SA parameters.

GATA binding protein 2 mediated ankyrin repeat domain containing 26 high expression in myeloid-derived cell lines

BACKGROUND

Thrombocytopenia 2, an autosomal dominant inherited disease characterized by moderate thrombocytopenia, predisposition to myeloid malignancies and normal platelet size and function, can be caused by 5'-untranslated region (UTR) point mutations in ankyrin repeat domain containing 26 (ANKRD26). Runt related transcription factor 1 (RUNX1) and friend leukemia integration 1 (FLI1) have been identified as negative regulators of ANKRD26. However, the positive regulators of ANKRD26 are still unknown.

AIM

To prove the positive regulatory effect of GATA binding protein 2 (GATA2) on ANKRD26 transcription.

METHODS

Human induced pluripotent stem cells derived from bone marrow (hiPSC-BM) and urothelium (hiPSC-U) were used to examine the ANKRD26 expression pattern in the early stage of differentiation. Then, transcriptome sequencing of these iPSCs and three public transcription factor (TF) databases (Cistrome DB, animal TFDB and ENCODE) were used to identify potential TF candidates for ANKRD26. Furthermore, overexpression and dual-luciferase reporter experiments were used to verify the regulatory effect of the candidate TFs on ANKRD26. Moreover, using the GENT2 platform, we analyzed the relationship between ANKRD26 expression and overall survival in cancer patients.

RESULTS

In hiPSC-BMs and hiPSC-Us, we found that the transcription levels of ANKRD26 varied in the absence of RUNX1 and FLI1. We sequenced hiPSC-BM and hiPSC-U and identified 68 candidate TFs for ANKRD26. Together with three public TF databases, we found that GATA2 was the only candidate gene that could positively regulate ANKRD26. Using dual-luciferase reporter experiments, we showed that GATA2 directly binds to the 5'-UTR of ANKRD26 and promotes its transcription. There are two identified binding sites of GATA2 that are located 2 kb upstream of the TSS of ANKRD26. In addition, we discovered that high ANKRD26 expression is always related to a more favorable prognosis in breast and lung cancer patients.

CONCLUSION

We first discovered that the transcription factor GATA2 plays a positive role in ANKRD26 transcription and identified its precise binding sites at the promoter region, and we revealed the importance of ANKRD26 in many tissue-derived cancers.

Genome-wide association study reveals multiple loci for analgesia and opioid consumption behaviors associated with heroin vulnerability in outbred rats

The increased prevalence of opioid use disorder (OUD) has made it imperative to disentangle the biological mechanisms contributing to individual differences in susceptibility to OUD. OUD shows strong heritability, however genetic variants contributing toward vulnerability remain poorly defined. We performed a genome-wide association study (GWAS) using over 850 male and female heterogeneous stock rats to identify genes underlying behaviors associated with OUD such as analgesia, as well as heroin-taking, refraining and seeking behaviors. By using an animal model of OUD, we were able to identify genetic variants associated with distinct OUD behaviors while maintaining a uniform environment, an experimental design not easily achieved in humans. Furthermore, we applied an animal model capturing individual variation in OUD propensity to assess if GWAS results were associated with OUD vulnerable versus resilient behavioral phenotypes. Our findings confirm the heritability of several OUD-like behaviors, including overall phenotype. We identified several genetic variants associated with basal analgesia prior to heroin experience, heroin consumption, escalation of intake, and motivation to obtain heroin. Ets2 , a regulator of microglia functional plasticity, and its eQTL PCP4 were identified for heroin consumption, and were associated with an OUD vulnerable phenotype through phenotype wide association study analysis. Furthermore, the coding variant Phd1l2 and the eQTL MMP15 for break point are both known mediators of addiction-related behaviors, and correlated with OUD vulnerability. These findings identify novel genetic markers related to individual differences in the susceptibility to OUD-relevant behaviors.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared with its predecessor. Gene annotations are now more complete, improving the mapping precision of genomic, transcriptomic, and proteomics datasets. We jointly analyzed 163 short-read whole-genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ∼20.0 million sequence variations, of which 18,700 are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genome-wide association study of delay discounting in Heterogenous Stock rats

Delay discounting refers to the behavioral tendency to devalue rewards as a function of their delay in receipt. Heightened delay discounting has been associated with substance use disorders, as well as multiple co-occurring psychopathologies. Genetic studies in humans and animal models have established that delay discounting is a heritable trait, but only a few specific genes have been associated with delay discounting. Here, we aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogenous Stock rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at variable delays. Preference switch points were calculated for each rat and both exponential and hyperbolic functions were fitted to these indifference points. Area under the curve (AUC) and the discounting parameter k of both functions were used as delay discounting measures. GWAS for AUC, exponential k, and indifference points for a short delay identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family of nucleoside sugar transporters, was the only gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression of that gene might be responsible for the association with behavior. The gene Adgrl3, which encodes a member of the latrophilin family of G-protein coupled receptors, was the only gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

A revamped rat reference genome improves the discovery of genetic diversity in laboratory rats

The seventh iteration of the reference genome assembly for Rattus norvegicus-mRatBN7.2-corrects numerous misplaced segments and reduces base-level errors by approximately 9-fold and increases contiguity by 290-fold compared to its predecessor. Gene annotations are now more complete, significantly improving the mapping precision of genomic, transcriptomic, and proteomics data sets. We jointly analyzed 163 short-read whole genome sequencing datasets representing 120 laboratory rat strains and substrains using mRatBN7.2. We defined ~20.0 million sequence variations, of which 18.7 thousand are predicted to potentially impact the function of 6,677 genes. We also generated a new rat genetic map from 1,893 heterogeneous stock rats and annotated transcription start sites and alternative polyadenylation sites. The mRatBN7.2 assembly, along with the extensive analysis of genomic variations among rat strains, enhances our understanding of the rat genome, providing researchers with an expanded resource for studies involving rats.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

Genetic background in the rat affects endocrine and metabolic outcomes of bisphenol F exposure

Environmental bisphenol compounds like bisphenol F (BPF) are endocrine-disrupting chemicals (EDCs) affecting adipose and classical endocrine systems. Genetic factors that influence EDC exposure outcomes are poorly understood and are unaccounted variables that may contribute to the large range of reported outcomes in the human population. We previously demonstrated that BPF exposure increased body growth and adiposity in male N/NIH heterogeneous stock (HS) rats, a genetically heterogeneous outbred population. We hypothesize that the founder strains of the HS rat exhibit EDC effects that were strain- and sex-dependent. Weanling littermate pairs of male and female ACI, BN, BUF, F344, M520, and WKY rats randomly received either vehicle (0.1% EtOH) or 1.125 mg BPF/l in 0.1% EtOH for 10 weeks in drinking water. Body weight and fluid intake were measured weekly, metabolic parameters were assessed, and blood and tissues were collected. BPF increased thyroid weight in ACI males, thymus and kidney weight in BUF females, adrenal weight in WKY males, and possibly increased pituitary weight in BN males. BUF females also developed a disruption in activity and metabolic rate with BPF exposure. These sex- and strain-specific exposure outcomes illustrate that HS rat founders possess diverse bisphenol-exposure risk alleles and suggest that BPF exposure may intensify inherent organ system dysfunction existing in the HS rat founders. We propose that the HS rat will be an invaluable model for dissecting gene EDC interactions on health.

The genetic architecture of behavioral canalization

Behaviors are components of fitness and contribute to adaptive evolution. Behaviors represent the interactions of an organism with its environment, yet innate behaviors display robustness in the face of environmental change, which we refer to as 'behavioral canalization'. We hypothesize that positive selection of hub genes of genetic networks stabilizes the genetic architecture for innate behaviors by reducing variation in the expression of interconnected network genes. Robustness of these stabilized networks would be protected from deleterious mutations by purifying selection or suppressing epistasis. We propose that, together with newly emerging favorable mutations, epistatically suppressed mutations can generate a reservoir of cryptic genetic variation that could give rise to decanalization when genetic backgrounds or environmental conditions change to allow behavioral adaptation.