Linkage screen for BMD phenotypes in male and female COP and DA rat strains

Chromosomal Substitution Strategies to Localize Genomic Regions Related to Complex Traits

Chromosomal substitution strategies provide a powerful tool to anonymously reveal the relationship between DNA sequence variants and a normal or disease phenotype of interest. Even in this age of CRISPR-Cas9 genome engineering, the knockdown or overexpression of a gene provides relevant information to our understanding of complex disease only when a close association of an allelic variant with the phenotype has first been established. Limitations of genetic linkage approaches led to the development of more efficient breeding strategies to substitute chromosomal segments from one animal strain into the genetic background of a different strain, enabling a direct comparison of the phenotypes of the strains with variant(s) that differ only at a defined locus. This substitution can be a whole chromosome (consomic), a part of a chromosome (congenic), or as small as only a single or several alleles (subcongenics). In contrast to complete knockout of a specific candidate gene of interest, which simply studies the effects of complete elimination of the gene, the substitution of naturally occurring variants can provide special insights into the functional actions of wild-type alleles. Strategies for production of these inbred strains are reviewed, and a number of examples are used to illustrate the utility of these model systems. Consomic/congenic strains provide a number of experimental advantages in the study of functions of genes and their variants, which are emphasized in this article, such as replication of experimental studies; determination of temporal relationships throughout a life; rigorously controlled experiments in which relations between genotype and phenotype can be tested with the confounding effects of heterogeneous genetic backgrounds, both targeted and multilayered; and "omic" studies performed at many levels of functionality, from molecules to organelles, cells to organs, and organs to organismal behavior across the life span. The application of chromosomal substitution strategies and development of consomic/congenic rat and mouse strains have greatly expanded our knowledge of genomic variants and their phenotypic relationship to physiological functions and to complex diseases such as hypertension and cancer. © 2020 American Physiological Society. Compr Physiol 10:365-388, 2020.

Uncovering the liver's role in immunity through RNA co-expression networks

Gene co-expression analysis has proven to be a powerful tool for ascertaining the organization of gene products into networks that are important for organ function. An organ, such as the liver, engages in a multitude of functions important for the survival of humans, rats, and other animals; these liver functions include energy metabolism, metabolism of xenobiotics, immune system function, and hormonal homeostasis. With the availability of organ-specific transcriptomes, we can now examine the role of RNA transcripts (both protein-coding and non-coding) in these functions. A systems genetic approach for identifying and characterizing liver gene networks within a recombinant inbred panel of rats was used to identify genetically regulated transcriptional networks (modules). For these modules, biological consensus was found between functional enrichment analysis and publicly available phenotypic quantitative trait loci (QTL). In particular, the biological function of two liver modules could be linked to immune response. The eigengene QTLs for these co-expression modules were located at genomic regions coincident with highly significant phenotypic QTLs; these phenotypes were related to rheumatoid arthritis, food preference, and basal corticosterone levels in rats. Our analysis illustrates that genetically and biologically driven RNA-based networks, such as the ones identified as part of this research, provide insight into the genetic influences on organ functions. These networks can pinpoint phenotypes that manifest through the interaction of many organs/tissues and can identify unannotated or under-annotated RNA transcripts that play a role in these phenotypes.

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 role of statistical epistasis in the genetic architecture of complex traits has been of great interest to the genetics community since Fisher introduced the concept in 1918. However, assessing epistasis in human and model organism populations has been impeded by limited statistical power. To mitigate this limitation, we analyzed bone and body composition traits in an unusually large mouse intercross population of over 2000 mice, paired with a recently-developed computational approach that leverages information to detect interactions across multiple phenotypes. We discovered a large network of highly significant genetic interactions between variants that influence complex body composition traits. Although epistasis was abundant, the interaction network was dominated by epistasis that stabilizes phenotypes by reducing phenotypic deviation from the parent strains. Nevertheless, the observed network provides an overview of genetic architecture and specific hypotheses of how QTL combine to affect phenotypes. These findings suggest that epistatic effects are generally of lesser magnitude than main QTL effects, and therefore are unlikely to account for major components of variance, but also reinforce genetic interaction analysis as a potent tool for dissecting the biology of complex traits.

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis.

Genes influencing spinal bone mineral density in inbred F344, LEW, COP, and DA rats

Previously, we identified the regions of chromosomes 10q12-q31 and 15p16-q21 harbor quantitative trait loci (QTLs) for lumbar volumetric bone mineral density (vBMD) in female F2 rats derived from Fischer 344 (F344) x Lewis (LEW) and Copenhagen 2331 (COP) x Dark Agouti (DA) crosses. The purpose of this study is to identify the candidate genes within these QTL regions contributing to the variation in lumbar vBMD. RNA was extracted from bone tissue of F344, LEW, COP, and DA rats. Microarray analysis was performed using Affymetrix Rat Genome 230 2.0 Arrays. Genes differentially expressed among the rat strains were then ranked based on the strength of the correlation with lumbar vBMD in F2 animals derived from these rats. Quantitative PCR (qPCR) analysis was performed to confirm the prioritized candidate genes. A total of 285 genes were differentially expressed among all strains of rats with a false discovery rate less than 10%. Among these genes, 18 candidate genes were prioritized based on their strong correlation (r (2) > 0.90) with lumbar vBMD. Of these, 14 genes (Akap1, Asgr2, Esd, Fam101b, Irf1, Lcp1, Ltc4s, Mdp-1, Pdhb, Plxdc1, Rabep1, Rhot1, Slc2a4, Xpo4) were confirmed by qPCR. We identified several novel candidate genes influencing spinal vBMD in rats.

Bone mineral density variation in men is influenced by sex-specific and non sex-specific quantitative trait loci

INTRODUCTION

A major predictor of age-related osteoporotic fracture is peak areal bone mineral density (aBMD) which is a highly heritable trait. However, few linkage and association studies have been performed in men to identify the genes contributing to normal variation in aBMD. The aim of this study was to perform a genome wide scan in healthy men to identify quantitative trait loci (QTL) that were significantly linked to aBMD and to test whether any of these might be sex-specific.

METHODS

aBMD at the spine and hip were measured in 515 pairs of brothers, aged 18-61 (405 white pairs, 110 black pairs). Linkage analysis in the brother sample was compared with results in a previously published sample of 774 sister pairs to identify sex-specific quantitative trait loci (QTL).

RESULTS

A genome wide scan identified significant QTL (LOD>3.6) for aBMD on chromosomes 4q21 (hip), 7q34 (spine), 14q32 (hip), 19p13 (hip), 21q21 (hip), and 22q13 (hip). Analysis suggested that the QTL on chromosomes 7q34, 14q32, and 21q21 were male-specific whereas the others were not sex-specific.

CONCLUSIONS

This study demonstrates that six QTL were significantly linked with aBMD in men. One was linked to the spine and five were linked to the hip. When compared to published data in women from the same geographical region, the QTL on chromosomes 7, 14 and 21 were male-specific. The occurrence of sex-specific genes in humans for aBMD has important implications for the pathogenesis and treatment of osteoporosis.

Epistasis between QTLs for bone density variation in Copenhagen x dark agouti F2 rats

The variation in several of the risk factors for osteoporotic fracture, including bone mineral density (BMD), has been shown to be strongly influenced by genetic differences. However, the genetic architecture of BMD is complex in both humans and in model organisms. We previously reported quantitative trait locus (QTL) results for BMD from a genome screen of 828 F2 progeny of Copenhagen and dark agouti rats. These progeny also provide an excellent opportunity to search for epistatic effects, or interaction between genetic loci, that contribute to fracture risk. Microsatellite marker data from a 20-cM genome screen was analyzed along with weight-adjusted bone density (DXA and pQCT) phenotypic data using the R/qtl software package. Genotype and phenotype data were permuted to determine genome-wide significance thresholds for the full model and epistasis (interaction) LOD scores corresponding to an alpha level of 0.01. A novel locus on chromosome 15 and a previously reported chromosome 14 QTL demonstrated a strong epistatic effect on BMD at the femur by DXA (LOD = 5.4). Two novel QTLs on chromosomes 2 and 12 were found to interact to affect total BMD at the femur midshaft by pQCT (LOD = 5.0). These results provide new information regarding the mode of action of previously identified QTL in the rat, as well as identifying novel loci that act in combination with known QTL or with other novel loci to contribute to BMD variation.