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

Novel insights into the genetic architecture and mechanisms of host/microbiome interactions from a multi-cohort analysis of outbred laboratory rats

The intestinal microbiome influences health and disease. Its composition is affected by host genetics and environmental exposures. Understanding host genetic effects is critical but challenging in humans, due to the difficulty of detecting, mapping and interpreting them. To address this, we analysed host genetic effects in four cohorts of outbred laboratory rats exposed to distinct but controlled environments. We found that polygenic host genetic effects were consistent across environments. We identified three replicated microbiome-associated loci. One involved a sialyltransferase gene and Paraprevotella and we found a similar association, between ST6GAL1 and Paraprevotella, in a human cohort. Given Paraprevotella's known immunity-potentiating functions, this suggests ST6GAL1's effects on IgA nephropathy and COVID-19 breakthrough infections may be mediated by Paraprevotella. Moreover, we found evidence of indirect genetic effects on microbiome phenotypes, which substantially increased their total genetic variance. Finally, we identified a novel mechanism whereby indirect genetic effects can contribute to "missing heritability".

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

Oxycodone abuse frequently begins with prescription oral oxycodone, yet vulnerability factors (e.g. sex, genetics) 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-, and 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). Both sexes increased intake between 4-h and 16-h sessions (p < 2e-16), but a subset of strains augmented intake at 16-h (p = 0.0005). Heritability (h2) of active licks during 4-h sessions at increasing oxycodone dose ranged from 0.30 to 0.53. Under a progressive ratio (PR) schedule, breakpoints were strain-dependent (p < 2e-16). Cued reinstatement was greater in females (p < 0.001). Naive rats were assessed using elevated plus maze (EPM), open field (OF), and novel object interaction (NOI) tests. We correlated these behaviors with 28 parameters of oxycodone SA. Anxiety-defining EPM traits were most associated with SA in both sexes, whereas OF and NOI traits were more associated with SA 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 anxiety correlates with multiple SA parameters across strains.

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

The increased prevalence of opioid use disorder (OUD) makes it imperative to disentangle the biological mechanisms contributing to individual differences in OUD vulnerability. OUD shows strong heritability, however genetic variants contributing to vulnerability remain poorly defined. We performed a genome-wide association study using over 850 male and female heterogeneous stock (HS) rats to identify genes underlying behaviors associated with OUD such as nociception, as well as heroin-taking, extinction 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 used a novel non-linear network-based clustering approach to characterize rats based on OUD vulnerability to assess genetic variants associated with OUD susceptibility. Our findings confirm the heritability of several OUD-like behaviors, including OUD susceptibility. Additionally, several genetic variants associated with nociceptive threshold prior to heroin experience, heroin consumption, escalation of intake, and motivation to obtain heroin were identified. Tom1, a microglial component, was implicated for nociception. Several genes involved in dopaminergic signaling, neuroplasticity and substance use disorders, including Brwd1, Pcp4, Phb1l2 and Mmp15 were implicated for the heroin traits. Additionally, an OUD vulnerable phenotype was associated with genetic variants for consumption and break point, suggesting a specific genetic contribution for OUD-like traits contributing to vulnerability. Together, these findings identify novel genetic markers related to the susceptibility to OUD-relevant behaviors in HS rats.

A cost-effective, high-throughput, highly accurate genotyping method for outbred populations

Affordable sequencing and genotyping methods are essential for large-scale genome-wide association studies. While genotyping microarrays and reference panels for imputation are available for human subjects, nonhuman model systems often lack such options. Our lab previously demonstrated an efficient and cost-effective method to genotype heterogeneous stock rats using double-digest genotyping by sequencing. However, low-coverage whole-genome sequencing offers an alternative method that has several advantages. Here, we describe a cost-effective, high-throughput, high-accuracy genotyping method for N/NIH heterogeneous stock rats that can use a combination of sequencing data previously generated by double-digest genotyping by sequencing and more recently generated by low-coverage whole-genome sequencing data. Using double-digest genotyping-by-sequencing data from 5,745 heterogeneous stock rats (mean 0.21× coverage) and low-coverage whole-genome sequencing data from 8,760 heterogeneous stock rats (mean 0.27× coverage), we can impute 7.32 million biallelic single-nucleotide polymorphisms with a concordance rate > 99.76% compared to high-coverage (mean 33.26× coverage) whole-genome sequencing data for a subset of the same individuals. Our results demonstrate the feasibility of using sequencing data from double-digest genotyping by sequencing or low-coverage whole-genome sequencing for accurate genotyping and demonstrate techniques that may also be useful for other genetic studies in nonhuman subjects.

STICI: Split-Transformer with integrated convolutions for genotype imputation

Despite advances in sequencing technologies, genome-scale datasets often contain missing bases and genomic segments, hindering downstream analyses. Genotype imputation addresses this issue and has been a cornerstone pre-processing step in genetic and genomic studies. Although various methods have been widely adopted for genotype imputation, it remains challenging to impute certain genomic regions and large structural variants. Here, we present a transformer-based framework, named STICI, for accurate genotype imputation. STICI models automatically learn genome-wide patterns of linkage disequilibrium, evidenced by much higher imputation accuracy in regions with highly linked variants. Our imputation results on the human 1000 Genomes Project and non-human genomes show that STICI can achieve high imputation accuracy comparable to the state-of-the-art genotype imputation methods, with the additional capability to impute multi-allelic variants and various types of genetic variants. STICI can be trained for any collection of genomes automatically using self-supervision. Moreover, STICI shows excellent performance without needing any special presuppositions about the underlying patterns in collections of non-human genomes, pointing to adaptability and applications of STICI to impute missing genotypes in any species.

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

The increased prevalence of opioid use disorder (OUD) makes it imperative to disentangle the biological mechanisms contributing to individual differences in OUD vulnerability. OUD shows strong heritability, however genetic variants contributing toward vulnerability remain poorly defined. We performed a genome-wide association study using over 850 male and female heterogeneous stock (HS) rats to identify genes underlying behaviors associated with OUD such as nociception, as well as heroin-taking, extinction 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 used a novel non-linear network-based clustering approach to characterize rats based on OUD vulnerability to assess genetic variants associated with OUD susceptibility. Our findings confirm the heritability of several OUD-like behaviors, including OUD susceptibility. Additionally, several genetic variants associated with nociceptive threshold prior to heroin experience, heroin consumption, escalation of intake, and motivation to obtain heroin were identified. Tom1, a microglial component, was implicated for nociception. Several genes involved in dopaminergic signaling, neuroplasticity and substance use disorders, including Brwd1, Pcp4, Phb1l2 and Mmp15 were implicated for the heroin traits. Additionally, an OUD vulnerable phenotype was associated with genetic variants for consumption and break point, suggesting a specific genetic contribution for OUD-like traits contributing to vulnerability. Together, these findings identify novel genetic markers related to the susceptibility to OUD-relevant behaviors in HS rats.

Genome-wide association study of delay discounting in Heterogeneous 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 and multiple co-occurring psychopathologies. Human and animal genetic studies have established that delay discounting is heritable, but only a few associated genes have been identified. We aimed to identify novel genetic loci associated with delay discounting through a genome-wide association study (GWAS) using Heterogeneous Stock (HS) rats, a genetically diverse outbred population derived from eight inbred founder strains. We assessed delay discounting in 650 male and female HS rats using an adjusting amount procedure in which rats chose between smaller immediate sucrose rewards or a larger reward at various delays. Preference switch points were calculated 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 one indifference point identified significant loci on chromosomes 20 and 14. The gene Slc35f1, which encodes a member of the solute carrier family, was the sole gene within the chromosome 20 locus. That locus also contained an eQTL for Slc35f1, suggesting that heritable differences in the expression might be responsible for the association with behavior. Adgrl3, which encodes a latrophilin subfamily G-protein coupled receptor, was the sole gene within the chromosome 14 locus. These findings implicate novel genes in delay discounting and highlight the need for further exploration.

A Cost-effective, High-throughput, Highly Accurate Genotyping Method for Outbred Populations

Affordable sequencing and genotyping methods are essential for large scale genome-wide association studies. While genotyping microarrays and reference panels for imputation are available for human subjects, non-human model systems often lack such options. Our lab previously demonstrated an efficient and cost-effective method to genotype heterogeneous stock rats using double-digest genotyping-by-sequencing. However, low-coverage whole-genome sequencing offers an alternative method that has several advantages. Here, we describe a cost-effective, high-throughput, high-accuracy genotyping method for N/NIH heterogeneous stock rats that can use a combination of sequencing data previously generated by double-digest genotyping-by-sequencing and more recently generated by low-coverage whole-genome-sequencing data. Using double-digest genotyping-by-sequencing data from 5,745 heterogeneous stock rats (mean 0.21x coverage) and low-coverage whole-genome-sequencing data from 8,760 heterogeneous stock rats (mean 0.27x coverage), we can impute 7.32 million bi-allelic single-nucleotide polymorphisms with a concordance rate >99.76% compared to high-coverage (mean 33.26x coverage) whole-genome sequencing data for a subset of the same individuals. Our results demonstrate the feasibility of using sequencing data from double-digest genotyping-by-sequencing or low-coverage whole-genome-sequencing for accurate genotyping, and demonstrate techniques that may also be useful for other genetic studies in non-human subjects.

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.

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.

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.

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.