Variation in histone configurations correlates with gene expression across nine inbred strains of mice.

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Systems genetics uncover new loci containing functional gene candidates in Mycobacterium tuberculosis-infected Diversity Outbred mice.

Mycobacterium tuberculosis infects two billion people across the globe, and results in 8-9 million new tuberculosis (TB) cases and 1-1.5 million deaths each year. Most patients have no known genetic basis that predisposes them to disease. Here, we investigate the complex genetic basis of pulmonary TB by modelling human genetic diversity with the Diversity Outbred mouse population. When infected with M. tuberculosis, one-third develop early onset, rapidly progressive, necrotizing granulomas and succumb within 60 days. The remaining develop non-necrotizing granulomas and survive longer than 60 days. Genetic mapping using immune and inflammatory mediators; and clinical, microbiological, and granuloma correlates of disease identified five new loci on mouse chromosomes 1, 2, 4, 16; and three known loci on chromosomes 3 and 17. Further, multiple positively correlated traits shared loci on chromosomes 1, 16, and 17 and had similar patterns of allele effects, suggesting these loci contain critical genetic regulators of inflammatory responses to M. tuberculosis. To narrow the list of candidate genes, we used a machine learning strategy that integrated gene expression signatures from lungs of M. tuberculosis-infected Diversity Outbred mice with gene interaction networks to generate scores representing functional relationships. The scores were used to rank candidates for each mapped trait, resulting in 11 candidate genes: Ncf2, Fam20b, S100a8, S100a9, Itgb5, Fstl1, Zbtb20, Ddr1, Ier3, Vegfa, and Zfp318. Although all candidates have roles in infection, inflammation, cell migration, extracellular matrix remodeling, or intracellular signaling, and all contain single nucleotide polymorphisms (SNPs), SNPs in only four genes (S100a8, Itgb5, Fstl1, Zfp318) are predicted to have deleterious effects on protein functions. We performed methodological and candidate validations to (i) assess biological relevance of predicted allele effects by showing that Diversity Outbred mice carrying PWK/PhJ alleles at the H-2 locus on chromosome 17 QTL have shorter survival; (ii) confirm accuracy of predicted allele effects by quantifying S100A8 protein in inbred founder strains; and (iii) infection of C57BL/6 mice deficient for the S100a8 gene. Overall, this body of work demonstrates that systems genetics using Diversity Outbred mice can identify new (and known) QTLs and functionally relevant gene candidates that may be major regulators of complex host-pathogens interactions contributing to granuloma necrosis and acute inflammation in pulmonary TB.

Immunogenicity of the Ad26.COV2.S Vaccine for COVID-19.

Importance: Control of the global COVID-19 pandemic will require the development and deployment of safe and effective vaccines. Objective: To evaluate the immunogenicity of the Ad26.COV2.S vaccine (Janssen/Johnson & Johnson) in humans, including the kinetics, magnitude, and phenotype of SARS-CoV-2 spike-specific humoral and cellular immune responses. Design, Setting, and Participants: Twenty-five participants were enrolled from July 29, 2020, to August 7, 2020, and the follow-up for this day 71 interim analysis was completed on October 3, 2020; follow-up to assess durability will continue for 2 years. This study was conducted at a single clinical site in Boston, Massachusetts, as part of a randomized, double-blind, placebo-controlled phase 1 clinical trial of Ad26.COV2.S. Interventions: Participants were randomized to receive 1 or 2 intramuscular injections with 5 × 1010 viral particles or 1 × 1011 viral particles of Ad26.COV2.S vaccine or placebo administered on day 1 and day 57 (5 participants in each group). Main Outcomes and Measures: Humoral immune responses included binding and neutralizing antibody responses at multiple time points following immunization. Cellular immune responses included immunospot-based and intracellular cytokine staining assays to measure T-cell responses. Results: Twenty-five participants were randomized (median age, 42; age range, 22-52; 52% women, 44% male, 4% undifferentiated), and all completed the trial through the day 71 interim end point. Binding and neutralizing antibodies emerged rapidly by day 8 after initial immunization in 90% and 25% of vaccine recipients, respectively. By day 57, binding and neutralizing antibodies were detected in 100% of vaccine recipients after a single immunization. On day 71, the geometric mean titers of spike-specific binding antibodies were 2432 to 5729 and the geometric mean titers of neutralizing antibodies were 242 to 449 in the vaccinated groups. A variety of antibody subclasses, Fc receptor binding properties, and antiviral functions were induced. CD4+ and CD8+ T-cell responses were induced. Conclusion and Relevance: In this phase 1 study, a single immunization with Ad26.COV2.S induced rapid binding and neutralization antibody responses as well as cellular immune responses. Two phase 3 clinical trials are currently underway to determine the efficacy of the Ad26.COV2.S vaccine. Trial Registration: ClinicalTrials.gov Identifier: NCT04436276.

Identification of novel loci controlling inflammatory bowel disease susceptibility utilizing the genetic diversity of wild-derived mice.

Inflammatory bowel disease (IBD) is a complex disorder that imposes a growing health burden. Multiple genetic associations have been identified in IBD, but the mechanisms underlying many of these associations are poorly understood. Animal models are needed to bridge this gap, but conventional laboratory mouse strains lack the genetic diversity of human populations. To more accurately model human genetic diversity, we utilized a panel of chromosome (Chr) substitution strains, carrying chromosomes from the wild-derived and genetically divergent PWD/PhJ (PWD) strain on the commonly used C57BL/6J (B6) background, as well as their parental B6 and PWD strains. Two models of IBD were used, TNBS- and DSS-induced colitis. Compared with B6 mice, PWD mice were highly susceptible to TNBS-induced colitis, but resistant to DSS-induced colitis. Using consomic mice, we identified several PWD-derived loci that exhibited profound effects on IBD susceptibility. The most pronounced of these were loci on Chr1 and Chr2, which yielded high susceptibility in both IBD models, each acting at distinct phases of the disease. Leveraging transcriptomic data from B6 and PWD immune cells, together with a machine learning approach incorporating human IBD genetic associations, we identified lead candidate genes, including Itga4, Pip4k2a, Lcn10, Lgmn, and Gpr65.

Weak Epistasis Generally Stabilizes Phenotypes in a Mouse Intercross.

The extent and strength of epistasis is commonly unresolved in genetic studies, and observed epistasis is often difficult to interpret in terms of biological consequences or overall genetic architecture. We investigated the prevalence and consequences of epistasis by analyzing four body composition phenotypes--body weight, body fat percentage, femoral density, and femoral circumference--in a large F2 intercross of B6-lit/lit and C3.B6-lit/lit mice. We used Combined Analysis of Pleiotropy and Epistasis (CAPE) to examine interactions for the four phenotypes simultaneously, which revealed an extensive directed network of genetic loci interacting with each other, circulating IGF1, and sex to influence these phenotypes. The majority of epistatic interactions had small effects relative to additive effects of individual loci, and tended to stabilize phenotypes towards the mean of the population rather than extremes. Interactive effects of two alleles inherited from one parental strain commonly resulted in phenotypes closer to the population mean than the additive effects from the two loci, and often much closer to the mean than either single-locus model. Alternatively, combinations of alleles inherited from different parent strains contribute to more extreme phenotypes not observed in either parental strain. This class of phenotype-stabilizing interactions has effects that are close to additive and are thus difficult to detect except in very large intercrosses. Nevertheless, we found these interactions to be useful in generating hypotheses for functional relationships between genetic loci. Our findings suggest that while epistasis is often weak and unlikely to account for a large proportion of heritable variance, even small-effect genetic interactions can facilitate hypotheses of underlying biology in well-powered studies.

The detection and characterization of pleiotropy: discovery, progress, and promise.

The impact of a single genetic locus on multiple phenotypes, or pleiotropy, is an important area of research. Biological systems are dynamic complex networks, and these networks exist within and between cells. In humans, the consideration of multiple phenotypes such as physiological traits, clinical outcomes and drug response, in the context of genetic variation, can provide ways of developing a more complete understanding of the complex relationships between genetic architecture and how biological systems function in health and disease. In this article, we describe recent studies exploring the relationships between genetic loci and more than one phenotype. We also cover methodological developments incorporating pleiotropy applied to model organisms as well as humans, and discuss how stepping beyond the analysis of a single phenotype leads to a deeper understanding of complex genetic architecture.

CAPE: an R package for combined analysis of pleiotropy and epistasis.

Contemporary genetic studies are revealing the genetic complexity of many traits in humans and model organisms. Two hallmarks of this complexity are epistasis, meaning gene-gene interaction, and pleiotropy, in which one gene affects multiple phenotypes. Understanding the genetic architecture of complex traits requires addressing these phenomena, but interpreting the biological significance of epistasis and pleiotropy is often difficult. While epistasis reveals dependencies between genetic variants, it is often unclear how the activity of one variant is specifically modifying the other. Epistasis found in one phenotypic context may disappear in another context, rendering the genetic interaction ambiguous. Pleiotropy can suggest either redundant phenotype measures or gene variants that affect multiple biological processes. Here we present an R package, R/cape, which addresses these interpretation ambiguities by implementing a novel method to generate predictive and interpretable genetic networks that influence quantitative phenotypes. R/cape integrates information from multiple related phenotypes to constrain models of epistasis, thereby enhancing the detection of interactions that simultaneously describe all phenotypes. The networks inferred by R/cape are readily interpretable in terms of directed influences that indicate suppressive and enhancing effects of individual genetic variants on other variants, which in turn account for the variance in quantitative traits. We demonstrate the utility of R/cape by analyzing a mouse backcross, thereby discovering novel epistatic interactions influencing phenotypes related to obesity and diabetes. R/cape is an easy-to-use, platform-independent R package and can be applied to data from both genetic screens and a variety of segregating populations including backcrosses, intercrosses, and natural populations. The package is freely available under the GPL-3 license at http://cran.r-project.org/web/packages/cape.

Network-based analysis predicts interacting genetic modifiers from a meta-mapping study of spike-wave discharge in mice.

Absence seizures are characterized by brief lapses in awareness accompanied by a hallmark spike-and-wave discharge (SWD) electroencephalographic pattern and are common to genetic generalized epilepsies (GGEs). While numerous genes have been associated with increased risk, including some Mendelian forms with a single causal allele, most cases of GGE are idiopathic and there are many unknown genetic modifiers of GGE influencing risk and severity. In a previous meta-mapping study, crosses between transgenic C57BL/6 and C3HeB/FeJ strains, each carrying one of three SWD-causing mutations (Gabrg2tm1Spet(R43Q) , Scn8a8j or Gria4spkw1 ), demonstrated an antagonistic epistatic interaction between loci on mouse chromosomes 2 and 7 influencing SWD. These results implicate universal modifiers in the B6 background that mitigate SWD severity through a common pathway, independent of the causal mutation. In this study, we prioritized candidate modifiers in these interacting loci. Our approach integrated human genome-wide association results with gene interaction networks and mouse brain gene expression to prioritize candidate genes and pathways driving variation in SWD outcomes. We considered candidate genes that are functionally associated with human GGE risk genes and genes with evidence for coding or non-coding allele effects between the B6 and C3H backgrounds. Our analyses output a summary ranking of gene pairs, one gene from each locus, as candidates for explaining the epistatic interaction. Our top-ranking gene pairs implicate microtubule function, cytoskeletal stability and cell cycle regulation as novel hypotheses about the source of SWD variation across strain backgrounds, which could clarify underlying mechanisms driving differences in GGE severity in humans.

Probing the basis of disease heterogeneity in multiple sclerosis using genetically diverse mice.

Multiple sclerosis (MS) is a complex disease with significant heterogeneity in disease course and progression. Genetic studies have identified numerous loci associated with MS risk, but the genetic basis of disease progression remains elusive. To address this, we leveraged the Collaborative Cross (CC), a genetically diverse mouse strain panel, and experimental autoimmune encephalomyelitis (EAE). The thirty-two CC strains studied captured a wide spectrum of EAE severity, trajectory, and presentation, including severe-progressive, monophasic, relapsing remitting, and axial rotary (AR)-EAE, accompanied by distinct immunopathology. Sex differences in EAE severity were observed in six strains. Quantitative trait locus analysis revealed distinct genetic linkage patterns for different EAE phenotypes, including EAE severity and incidence of AR-EAE. Machine learning-based approaches prioritized candidate genes for loci underlying EAE severity ( Abcc4 and Gpc6 ) and AR-EAE ( Yap1 and Dync2h1 ). This work expands the EAE phenotypic repertoire and identifies novel loci controlling unique EAE phenotypes, supporting the hypothesis that heterogeneity in MS disease course is driven by genetic variation. Summary: The genetic basis of disease heterogeneity in multiple sclerosis (MS) remains elusive. We leveraged the Collaborative Cross to expand the phenotypic repertoire of the experimental autoimmune encephalomyelitis (EAE) model of MS and identify loci controlling EAE severity, trajectory, and presentation.

Functional network changes in hippocampal CA1 after status epilepticus predict spatial memory deficits in rats.

Status epilepticus (SE) is a common neurological emergency, which has been associated with subsequent cognitive impairments. Neuronal death in hippocampal CA1 is thought to be an important mechanism of these impairments. However, it is also possible that functional interactions between surviving neurons are important. In this study we recorded in vivo single-unit activity in the CA1 hippocampal region of rats while they performed a spatial memory task. From these data we constructed functional networks describing pyramidal cell interactions. To build the networks, we used maximum entropy algorithms previously applied only to in vitro data. We show that several months following SE pyramidal neurons display excessive neuronal synchrony and less neuronal reactivation during rest compared with those in healthy controls. Both effects predict rat performance in a spatial memory task. These results provide a physiological mechanism for SE-induced cognitive impairment and highlight the importance of the systems-level perspective in investigating spatial cognition.

Speed modulation of hippocampal theta frequency correlates with spatial memory performance.

Hippocampal theta rhythm is believed to play a critical role in learning and memory. In animal models of temporal lobe epilepsy (TLE), there is evidence that alterations of hippocampal theta oscillations are involved in the cognitive impairments observed in this model. However, hippocampal theta frequency and amplitude at both the local field potential (LFP) and single unit level are strongly modulated by running speed, suggesting that the integration of locomotor information into memory processes may also be critical for hippocampal processing. Here, we investigate whether hippocampal speed-theta integration influences spatial memory and whether it could account for the memory deficits observed in TLE rats. LFPs were recorded in both Control (CTR) and TLE rats as they were trained in a spatial alternation task. TLE rats required more training sessions to perform the task at CTR levels. Both theta frequency and power were significantly lower in the TLE group. In addition, speed/theta frequency correlation coefficients and regression slopes varied from session to session and were worse in TLE. Importantly, there was a strong relationship between speed/theta frequency parameters and performance. Our analyses reveal that speed/theta frequency correlation with performance cannot merely be explained by the direct influence of speed on behavior. Therefore, variations in the coordination of theta frequency with speed may participate in learning and memory processes. Impairments of this function could explain at least partially memory deficits in epilepsy.

Intralesional application of ribavirin in two American flamingos (Phoenicopterus ruber) with poxvirus infection.

We report the successful treatment of poxvirus lesions in two juvenile American flamingos (Phoenicopterus ruber) with experimental low-dose intralesional ribavirin injection. In the first flamingo, the size and location of a beak verrucosity interfered with feeding, and after multiple surgical interventions, an experimental therapy of low-dose intralesional ribavirin was implemented with close blood parameter monitoring to minimize any potential side effects due to systemic antiviral administration. The second flamingo had a poxvirus lesion on the tibiotarsus, which recurred after unsuccessful conservative medical treatment and surgical intervention and a course of intralesional ribavirin therapy was implemented. Regression of the lesions in both flamingos commenced within 3 days of ribavirin treatment resulting in complete resolution within 6 weeks of onset of ribavirin treatment.

A genetic locus complements resistance to Bordetella pertussis-induced histamine sensitization.

Histamine plays pivotal role in normal physiology and dysregulated production of histamine or signaling through histamine receptors (HRH) can promote pathology. Previously, we showed that Bordetella pertussis or pertussis toxin can induce histamine sensitization in laboratory inbred mice and is genetically controlled by Hrh1/HRH1. HRH1 allotypes differ at three amino acid residues with P263-V313-L331 and L263-M313-S331, imparting sensitization and resistance respectively. Unexpectedly, we found several wild-derived inbred strains that carry the resistant HRH1 allotype (L263-M313-S331) but exhibit histamine sensitization. This suggests the existence of a locus modifying pertussis-dependent histamine sensitization. Congenic mapping identified the location of this modifier locus on mouse chromosome 6 within a functional linkage disequilibrium domain encoding multiple loci controlling sensitization to histamine. We utilized interval-specific single-nucleotide polymorphism (SNP) based association testing across laboratory and wild-derived inbred mouse strains and functional prioritization analyses to identify candidate genes for this modifier locus. Atg7, Plxnd1, Tmcc1, Mkrn2, Il17re, Pparg, Lhfpl4, Vgll4, Rho and Syn2 are candidate genes within this modifier locus, which we named Bphse, enhancer of Bordetella pertussis induced histamine sensitization. Taken together, these results identify, using the evolutionarily significant diversity of wild-derived inbred mice, additional genetic mechanisms controlling histamine sensitization.

SLC39A2 and FSIP1 polymorphisms as potential modifiers of arsenic-related bladder cancer.

Arsenic is a carcinogen that contaminates drinking water worldwide. Accumulating evidence suggests that both exposure and genetic factors may influence susceptibility to arsenic-induced malignancies. We sought to identify novel susceptibility loci for arsenic-related bladder cancer in a US population with low to moderate drinking water levels of arsenic. We first screened a subset of bladder cancer cases using a panel of approximately 10,000 non-synonymous single nucleotide polymorphisms (SNPs). Top ranking hits on the SNP array then were considered for further analysis in our population-based case-control study (n = 832 cases and 1,191 controls). SNPs in the fibrous sheath interacting protein 1 (FSIP1) gene (rs10152640) and the solute carrier family 39, member 2 (SLC39A2) in the ZIP gene family of metal transporters (rs2234636) were detected as potential hits in the initial scan and validated in the full case-control study. The adjusted odds ratio (OR) for the FSIP1 polymorphism was 2.57 [95% confidence interval (CI) 1.13, 5.85] for heterozygote variants (AG) and 12.20 (95% CI 2.51, 59.30) for homozygote variants (GG) compared to homozygote wild types (AA) in the high arsenic group (greater than the 90th percentile), and unrelated in the low arsenic group (equal to or below the 90th percentile) (P for interaction = 0.002). For the SLC39A2 polymorphism, the adjusted ORs were 2.96 (95% CI 1.23, 7.15) and 2.91 (95% CI 1.00, 8.52) for heterozygote (TC) and homozygote (CC) variants compared to homozygote wild types (TT), respectively, and close to one in the low arsenic group (P for interaction = 0.03). Our findings suggest novel variants that may influence risk of arsenic-associated bladder cancer and those who may be at greatest risk from this widespread exposure.

Differential daytime and night-time stomatal behavior in plants from North American deserts.

Night-time stomatal conductance (g(night)) occurs in many ecosystems, but the g(night) response to environmental drivers is relatively unknown, especially in deserts. Here, we conducted a Bayesian analysis of stomatal conductance (g) (N=5013) from 16 species in the Sonoran, Chihuahuan, Mojave and Great Basin Deserts (North America). We partitioned daytime g (g(day)) and g(night) responses by describing g as a mixture of two extreme (dark vs high light) behaviors. Significant g(night) was observed across 15 species, and the g(night) and g(day) behavior differed according to species, functional type and desert. The transition between extreme behaviors was determined by light environment, with the transition behavior differing between functional types and deserts. Sonoran and Chihuahuan C(4) grasses were more sensitive to vapor pressure difference (D) at night and soil water potential (Ψ(soil)) during the day, Great Basin C(3) shrubs were highly sensitive to D and Ψ(soil) during the day, and Mojave C(3) shrubs were equally sensitive to D and Ψ(soil) during the day and night. Species were split between the exhibition of isohydric or anisohydric behavior during the day. Three species switched from anisohydric to isohydric behavior at night. Such behavior, combined with differential D, Ψ(soil) and light responses, suggests that different mechanisms underlie g(day) and g(night) regulation.

Polystyrene Degradation by Exiguobacterium sp. RIT 594: Preliminary Evidence for a Pathway Containing an Atypical Oxygenase.

The widespread use of plastics has led to their increasing presence in the environment and subsequent pollution. Some microorganisms degrade plastics in natural ecosystems and the associated metabolic pathways can be studied to understand the degradation mechanisms. Polystyrene (PS) is one of the more recalcitrant plastic polymers that is degraded by only a few bacteria. Exiguobacterium is a genus of Gram-positive poly-extremophilic bacteria known to degrade PS, thus being of biotechnological interest, but its biochemical mechanisms of degradation have not yet been elucidated. Based solely on genome annotation, we initially proposed PS degradation by Exiguobacterium sp. RIT 594 via depolymerization and epoxidation catalyzed by a ring epoxidase. However, Fourier transform infrared (FTIR) spectroscopy analysis revealed an increase of carboxyl and hydroxyl groups with biodegradation, as well as of unconjugated C-C double bonds, both consistent with dearomatization of the styrene ring. This excludes any aerobic pathways involving side chain epoxidation and/or hydroxylation. Subsequent experiments confirmed that molecular oxygen is critical to PS degradation by RIT 594 because degradation ceased under oxygen-deprived conditions. Our studies suggest that styrene breakdown by this bacterium occurs via the sequential action of two enzymes encoded in the genome: an orphan aromatic ring-cleaving dioxygenase and a hydrolase.

Shadows of complexity: what biological networks reveal about epistasis and pleiotropy.

Pleiotropy, in which one mutation causes multiple phenotypes, has traditionally been seen as a deviation from the conventional observation in which one gene affects one phenotype. Epistasis, or gene-gene interaction, has also been treated as an exception to the Mendelian one gene-one phenotype paradigm. This simplified perspective belies the pervasive complexity of biology and hinders progress toward a deeper understanding of biological systems. We assert that epistasis and pleiotropy are not isolated occurrences, but ubiquitous and inherent properties of biomolecular networks. These phenomena should not be treated as exceptions, but rather as fundamental components of genetic analyses. A systems level understanding of epistasis and pleiotropy is, therefore, critical to furthering our understanding of human genetics and its contribution to common human disease. Finally, graph theory offers an intuitive and powerful set of tools with which to study the network bases of these important genetic phenomena.

System-Level Analysis of Alzheimer's Disease Prioritizes Candidate Genes for Neurodegeneration.

Alzheimer's disease (AD) is a debilitating neurodegenerative disorder. Since the advent of the genome-wide association study (GWAS) we have come to understand much about the genes involved in AD heritability and pathophysiology. Large case-control meta-GWAS studies have increased our ability to prioritize weaker effect alleles, while the recent development of network-based functional prediction has provided a mechanism by which we can use machine learning to reprioritize GWAS hits in the functional context of relevant brain tissues like the hippocampus and amygdala. In parallel with these developments, groups like the Alzheimer's Disease Neuroimaging Initiative (ADNI) have compiled rich compendia of AD patient data including genotype and biomarker information, including derived volume measures for relevant structures like the hippocampus and the amygdala. In this study we wanted to identify genes involved in AD-related atrophy of these two structures, which are often critically impaired over the course of the disease. To do this we developed a combined score prioritization method which uses the cumulative distribution function of a gene's functional and positional score, to prioritize top genes that not only segregate with disease status, but also with hippocampal and amygdalar atrophy. Our method identified a mix of genes that had previously been identified in AD GWAS including APOE, TOMM40, and NECTIN2(PVRL2) and several others that have not been identified in AD genetic studies, but play integral roles in AD-effected functional pathways including IQSEC1, PFN1, and PAK2. Our findings support the viability of our novel combined score as a method for prioritizing region- and even cell-specific AD risk genes.

Effects of kinship correction on inflation of genetic interaction statistics in commonly used mouse populations.

It is well understood that variation in relatedness among individuals, or kinship, can lead to false genetic associations. Multiple methods have been developed to adjust for kinship while maintaining power to detect true associations. However, relatively unstudied are the effects of kinship on genetic interaction test statistics. Here, we performed a survey of kinship effects on studies of six commonly used mouse populations. We measured inflation of main effect test statistics, genetic interaction test statistics, and interaction test statistics reparametrized by the Combined Analysis of Pleiotropy and Epistasis (CAPE). We also performed linear mixed model (LMM) kinship corrections using two types of kinship matrix: an overall kinship matrix calculated from the full set of genotyped markers, and a reduced kinship matrix, which left out markers on the chromosome(s) being tested. We found that test statistic inflation varied across populations and was driven largely by linkage disequilibrium. In contrast, there was no observable inflation in the genetic interaction test statistics. CAPE statistics were inflated at a level in between that of the main effects and the interaction effects. The overall kinship matrix overcorrected the inflation of main effect statistics relative to the reduced kinship matrix. The two types of kinship matrices had similar effects on the interaction statistics and CAPE statistics, although the overall kinship matrix trended toward a more severe correction. In conclusion, we recommend using an LMM kinship correction for both main effects and genetic interactions and further recommend that the kinship matrix be calculated from a reduced set of markers in which the chromosomes being tested are omitted from the calculation. This is particularly important in populations with substantial population structure, such as recombinant inbred lines in which genomic replicates are used.