Leah C. Solberg Woods

Fine-mapping a locus for glucose tolerance using heterogeneous stock rats

Heterogeneous stock (HS) animals provide the ability to map quantitative trait loci at high resolution [<5 Megabase (Mb)] in a relatively short time period. In the current study, we hypothesized that the HS rat colony would be useful for fine-mapping a region on rat chromosome 1 that has previously been implicated in glucose regulation. We administered a glucose tolerance test to 515 HS rats and genotyped these animals with 69 microsatellite markers, spaced an average distance of <1 Mb apart, on a 67 Mb region of rat chromosome 1. Using regression modeling of inferred haplotypes based on a hidden Markov model reconstruction and mixed model analysis in which we accounted for the complex family structure of the HS, we identified one sharp peak within this region. Using positional bootstrapping, we determined the most likely location of this locus is from 205.04 to 207.48 Mb. This work demonstrates the utility of HS rats for fine-mapping complex traits and emphasizes the importance of taking into account family structure when using highly recombinant populations.

Fine-mapping diabetes-related traits, including insulin resistance, in heterogeneous stock rats

Type 2 diabetes (T2D) is a disease of relative insulin deficiency resulting from both insulin resistance and beta cell failure. We have previously used heterogeneous stock (HS) rats to fine-map a locus for glucose tolerance. We show here that glucose intolerance in the founder strains of the HS colony is mediated by different mechanisms: insulin resistance in WKY and an insulin secretion defect in ACI, and we demonstrate a high degree of variability for measures of insulin resistance and insulin secretion in HS rats. As such, our goal was to use HS rats to fine-map several diabetes-related traits within a region on rat chromosome 1. We measured blood glucose and plasma insulin levels after a glucose tolerance test in 782 male HS rats. Using 97 SSLP markers, we genotyped a 68 Mb region on rat chromosome 1 previously implicated in glucose and insulin regulation. We used linkage disequilibrium mapping by mixed model regression with inferred descent to identify a region from 198.85 to 205.9 that contains one or more quantitative trait loci (QTL) for fasting insulin and a measure of insulin resistance, the quantitative insulin sensitivity check index. This region also encompasses loci identified for fasting glucose and Insulin_AUC (area under the curve). A separate <3 Mb QTL was identified for body weight. Using a novel penalized regression method we then estimated effects of alternative haplotype pairings under each locus. These studies highlight the utility of HS rats for fine-mapping genetic loci involved in the underlying causes of T2D.

Identification of a Novel Gene for Diabetic Traits in Rats, Mice, and Humans

The genetic basis of type 2 diabetes remains incompletely defined despite the use of multiple genetic strategies. Multiparental populations such as heterogeneous stocks (HS) facilitate gene discovery by allowing fine mapping to only a few megabases, significantly decreasing the number of potential candidate genes compared to traditional mapping strategies. In the present work, we employed expression and sequence analysis in HS rats (Rattus norvegicus) to identify Tpcn2 as a likely causal gene underlying a 3.1-Mb locus for glucose and insulin levels. Global gene expression analysis on liver identified Tpcn2 as the only gene in the region that is differentially expressed between HS rats with glucose intolerance and those with normal glucose regulation. Tpcn2 also maps as a cis-regulating expression QTL and is negatively correlated with fasting glucose levels. We used founder sequence to identify variants within this region and assessed association between 18 variants and diabetic traits by conducting a mixed-model analysis, accounting for the complex family structure of the HS. We found that two variants were significantly associated with fasting glucose levels, including a nonsynonymous coding variant within Tpcn2. Studies in Tpcn2 knockout mice demonstrated a significant decrease in fasting glucose levels and insulin response to a glucose challenge relative to those in wild-type mice. Finally, we identified variants within Tpcn2 that are associated with fasting insulin in humans. These studies indicate that Tpcn2 is a likely causal gene that may play a role in human diabetes and demonstrate the utility of multiparental populations for positionally cloning genes within complex loci.

Individual variation in the propensity to attribute incentive salience to a food cue: influence of sex.

There is considerable individual variation in the propensity of animals to attribute incentive salience to discrete reward cues, but to date most of this research has been conducted in male rats. The purpose of this study was to determine whether sex influences the propensity to attribute incentive salience to a food cue, using rats from two different outbred strains (Sprague-Dawley [SD] and Heterogeneous Stock [HS]). The motivational value of a food cue was assessed in two ways: (i) by the ability of the cue to elicit approach toward it and (ii) by its ability to act as a conditioned reinforcer. We found that female SD rats acquired Pavlovian conditioned approach behavior slightly faster than males, but no sex difference was detected in HS rats, and neither strain showed a sex difference in asymptotic performance of approach behavior. Moreover, female approach behavior did not differ across estrous cycle. Compared to males, females made more active responses during the test for conditioned reinforcement, although they made more inactive responses as well. We conclude that although there are small sex differences in performance on these tasks, these are probably not due to a notable sex difference in the propensity to attribute incentive salience to a food cue.

QTL mapping in outbred populations: successes and challenges

Quantitative trait locus (QTL) mapping in animal populations has been a successful strategy for identifying genomic regions that play a role in complex diseases and traits. When conducted in an F2 intercross or backcross population, the resulting QTL is frequently large, often encompassing 30 Mb or more and containing hundreds of genes. To narrow the locus and identify candidate genes, additional strategies are needed. Congenic strains have proven useful but work less well when there are multiple tightly linked loci, frequently resulting in loss of phenotype. As an alternative, we discuss the use of highly recombinant outbred models for directly fine-mapping QTL to only a few megabases. We discuss the use of several currently available models such as the advanced intercross (AI), heterogeneous stocks (HS), the diversity outbred (DO), and commercially available outbred stocks (CO). Once a QTL has been fine-mapped, founder sequence and expression QTL mapping can be used to identify candidate genes. In this regard, the large number of alleles found in outbred stocks can be leveraged to identify causative genes and variants. We end this review by discussing some important statistical considerations when analyzing outbred populations. Fine-resolution mapping in outbred models, coupled with full genome sequence, has already led to the identification of several underlying causative genes for many complex traits and diseases. These resources will likely lead to additional successes in the coming years.

Heterogeneous stock rats: a new model to study the genetics of renal phenotypes

Chronic kidney disease is a growing medical concern, with an estimated 25.6 million people in the United States exhibiting some degree of kidney injury and/or decline in kidney function. Animal models provide great insight into the study of the genetics of complex diseases. In particular, heterogeneous stock (HS) rats represent a unique genetic resource enabling rapid fine-mapping of complex traits. However, they have not been explored as a model to study renal phenotypes. To evaluate the usefulness of HS rats in the genetics of renal traits, a time course evaluation (weeks 8-40) was performed for several renal phenotypes. As expected, a large degree of variation was seen for most renal traits. By week 24, three (of 40) rats exhibited marked proteinuria that increased gradually until week 40 and ranged from 33.7 to 80.2 mg/24 h. Detailed histological analysis confirmed renal damage in these rats. In addition, several rats consistently exhibited significant hematuria (5/41). Interestingly, these rats were not the same rats that exhibited proteinuria, indicating that susceptibility to different types of kidney injury is likely segregating within the HS population. One HS rat exhibited unilateral renal agenesis (URA), which was accompanied by a significant degree of proteinuria and glomerular and tubulointerstitial injury. The parents of this HS rat were identified and bred further. Additional offspring of this pair were observed to exhibit URA at frequency between 40% and 60%. In summary, these novel data demonstrate that HS rats exhibit variation in proteinuria and other kidney-related traits, confirming that the model harbors susceptibility alleles for kidney injury and providing the basis for further genetic studies.

Distinct genetic regulation of progression of diabetes and renal disease in the Goto-Kakizaki rat.

Goto-Kakizaki (GK) rats develop early-onset type 2 diabetes (T2D) symptoms, with signs of diabetic nephropathy becoming apparent with aging. To determine whether T2D and renal disease share similar genetic architecture, we ran a quantitative trait locus (QTL) analysis in the F2 progeny of a GK x Brown Norway (BN) rat cross. Further, to determine whether genetic components change over time, we ran the QTL analysis on phenotypes collected longitudinally, at 3, 6, 9 and 12 mo, from the same animals. We confirmed three chromosomal regions that are linked to early diabetes phenotypes (chromosomes 1, 5, and 10) and a single region involved in the late progression of the disorder (chromosome 4). A single region was identified for the onset of the renal phenotype proteinuria (chromosome 5). This region overlaps the diabetic QTL, although it is not certain whether similar genes are involved in both phenotypes. A second QTL linked to the progression of the renal phenotype was found on chromosome 7. Linkage for triglyceride and cholesterol levels were also identified (chromosomes 7 and 8, respectively). These results demonstrate that, in general, different genetic components control diabetic and renal phenotypes in a diabetic nephropathy model. Furthermore, these results demonstrate that, over time, different genetic components are involved in progression of disease from those that were involved in disease onset. This observation would suggest that clinical studies collecting participants over a wide age distribution may be diluting genetic effects and reducing power to detect true effects.

Genetic Fine-Mapping and Identification of Candidate Genes and Variants for Adiposity Traits in Outbred Rats

Obesity is a major risk factor for multiple diseases and is in part heritable, yet the majority of causative genetic variants that drive excessive adiposity remain unknown. Here, outbred heterogeneous stock (HS) rats were used in controlled environmental conditions to fine‐map novel genetic modifiers of adiposity.

Spontaneous one-kidney rats are more susceptible to develop hypertension by DOCA-NaCl and subsequent kidney injury compared with uninephrectomized rats.

There is little clinical data of how hypertension may influence individuals with nephron deficiency in the context of being born with a single kidney. We recently developed a new rat model (the heterogeneous stock-derived model of unilateral renal agenesis rat) that is born with a single kidney and exhibits progressive kidney injury and decline in kidney function with age. We hypothesized that DOCA-salt would induce a greater increase in blood pressure and therefore accelerate the progression of kidney injury in rats born with a solitary kidney compared with rats that have undergone unilateral nephrectomy. Time course evaluation of blood pressure, kidney injury, and renal hemodynamics was performed in the following six groups of animals from weeks 13 to 18: 1) DOCA-treated rats with a solitary kidney (DOCA+S group), 2) placebo-treated rats with a solitary kidney, 3) DOCA-treated control rats with two kidneys (DOCA+C group), 4) placebo-treated control rats with two kidneys, 5) DOCA-treated rats with two kidneys that underwent uninephrectomy (DOCA+UNX8 group), and 6) placebo-treated rats with two kidneys that underwent uninephrectomy. DOCA+S rats demonstrated a significant rise (P < 0.05) in blood pressure (192 ± 4 mmHg), proteinuria (205 ± 31 mg/24 h), and a decline in glomerular filtration rate (600 ± 42 μl·min(-1)·g kidney weight(-1)) relative to the DOCA+UNX8 (173 ± 3 mmHg, 76 ± 26 mg/24 h, and 963 ± 36 μl·min(-1)·g kidney weight(-1)) and DOCA+C (154 ± 2 mmHg, 7 ± 1 mg/24 h, and 1,484 ± 121 μl·min(-1)·g kidney weight(-1)) groups. Placebo-treated groups showed no significant change among the three groups. An assessment of renal injury markers via real-time PCR/Western blot analysis and histological analysis was concordant with the measured physiological parameters. In summary, congenital solitary kidney rats are highly susceptible to the induction of hypertension compared with uninephrectomized rats, suggesting that low nephron endowment is an important driver of elevated blood pressure, hastening nephron injury through the transmission of elevated systemic blood pressure and thereby accelerating decline in kidney function.

Region-Based Convolutional Neural Nets for Localization of Glomeruli in Trichrome-Stained Whole Kidney Sections

Background Histologic examination of fixed renal tissue is widely used to assess morphology and the progression of disease. Commonly reported metrics include glomerular number and injury. However, characterization of renal histology is a time-consuming and user-dependent process. To accelerate and improve the process, we have developed a glomerular localization pipeline for trichrome-stained kidney sections using a machine learning image classification algorithm.Methods We prepared 4-μm slices of kidneys from rats of various genetic backgrounds that were subjected to different experimental protocols and mounted the slices on glass slides. All sections used in this analysis were trichrome stained and imaged in bright field at a minimum resolution of 0.92 μm per pixel. The training and test datasets for the algorithm comprised 74 and 13 whole renal sections, respectively, totaling over 28,000 glomeruli manually localized. Additionally, because this localizer will be ultimately used for automated assessment of glomerular injury, we assessed bias of the localizer for preferentially identifying healthy or damaged glomeruli.Results Localizer performance achieved an average precision and recall of 96.94% and 96.79%, respectively, on whole kidney sections without evidence of bias for or against glomerular injury or the need for manual preprocessing.Conclusions This study presents a novel and robust application of convolutional neural nets for the localization of glomeruli in healthy and damaged trichrome-stained whole-renal section mounts and lays the groundwork for automated glomerular injury scoring.