Isolation and high-throughput sequencing of two closely linked epistatic hypertension susceptibility loci with a panel of bicongenic strains.

Interactions or epistasis between genetic factors may contribute to "missing heritability." While linkage analyses detect epistasis, defining the limits of the interacting segments poses a significant challenge especially when the interactions are between loci in close proximity. The goal of the present study was to isolate two such epistatic blood pressure (BP) loci on rat chromosome 5. A panel of S.LEW bicongenic strains along with the corresponding monocongenic strains was constructed. BP of each set comprising of one bicongenic and two corresponding monocongenic strains were determined along with the parental Salt-sensitive (S) strain. Epistasis was observed in one out of four sets of congenic strains, wherein systolic blood pressures (SBP) of the two monocongenic strains S.LEW(5)x6Bx9x5a and S.LEW(5)x6Bx9x5b were comparable to that of S, but the SBP of the bicongenic strain S.LEW(5)x6Bx9x5 (157 ± 4.3 mmHg) was significantly lower than that of S (196 ± 6.8 mmHg, P < 0.001). A two-way ANOVA indicated significant interactions between the LEW alleles at the two loci. The interacting loci were 2.02 Mb apart and located within genomic segments spanning 7.77 and 4.18 Mb containing 7,360 and 2,753 candidate variants, respectively. The current study demonstrates definitive evidence for epistasis and provides genetic tools for further dissection of the isolated epistatic BP loci.

Construction of two novel reciprocal conplastic rat strains and characterization of cardiac mitochondria.

Because of the lack of appropriate animal models, the potentially causal contributions of inherited mitochondrial genomic factors to complex traits are less well studied compared with inherited nuclear genomic factors. We previously detected variations between the mitochondrial DNA (mtDNA) of the Dahl salt-sensitive (S) rat and the spontaneously hypertensive rat (SHR). Specifically, multiple variations were detected in mitochondrial genes coding for subunits of proteins essential for electron transport, in mitochondrial reactive oxygen species production, and within the D-loop region. To evaluate the effects of these mtDNA variations in the absence of the corresponding nuclear genomic factors as confounding variables, novel reciprocal strains of S and SHR were constructed and characterized. When compared with that of the S rat, the heart tissue from the S.SHR(mt) conplastic strain wherein the mtDNA of the S rat was substituted with that of the SHR had a significant increase in mtDNA copy number and decrease in mitochondrial reactive oxygen species production. A corresponding increase in aerobic treadmill running capacity and a significant increase in survival that was not related to changes in blood pressure were observed in the S.SHR(mt) rats compared with the S rat. The reciprocal SHR.S(mt) rats did not differ from the SHR in any phenotype tested, suggesting lower penetrance of the S mtDNA on the nuclear genomic background of the SHR. These novel conplastic strains serve as invaluable tools to further dissect the relationship between heart function, aerobic fitness, cardiovascular disease progression, and mortality.

Positional identification of variants of Adamts16 linked to inherited hypertension.

A previously reported blood pressure (BP) quantitative trait locus on rat Chromosome 1 was isolated in a short congenic segment spanning 804.6 kb. The 804.6 kb region contained only two genes, LOC306664 and LOC306665. LOC306664 is predicted to translate into A Disintegrin-like and Metalloproteinase with Thrombospondin Motifs-16 (Adamts16). LOC306665 is a novel gene. All predicted exons of both LOC306664 and LOC306665 were sequenced. Non-synonymous variants were identified in only one of these genes, LOC306664. These variants were naturally existing polymorphisms among inbred, outbred and wild rats. The full-length rat transcript of Adamts16 was detected in multiple tissues. Similar to ADAMTS16 in humans, expression of Adamts16 was prominent in the kidney. Renal transcriptome analysis suggested that a network of genes related to BP was differential between congenic and S rats. These genes were also differentially expressed between kidney cell lines with or without knock-down of Adamts16. Adamts16 is conserved between rats and humans. It is a candidate gene within the homologous region on human Chromosome 5, which is linked to systolic and diastolic BP in the Quebec Family Study. Multiple variants, including an Ala to Pro variant in codon 90 (rs2086310) of human ADAMTS16, were associated with human resting systolic BP (SBP). Replication study in GenNet confirmed the association of two variants of ADAMTS16 with SBP, including rs2086310. Overall, our report represents a high resolution positional cloning and translational study for Adamts16 as a candidate gene controlling BP.

Mitochondrial polymorphisms in rat genetic models of hypertension.

Hypertension is a complex trait that has been studied extensively for genetic contributions of the nuclear genome. We examined mitochondrial genomes of the hypertensive strains: the Dahl Salt-Sensitive (S) rat, the Spontaneously Hypertensive Rat (SHR), and the Albino Surgery (AS) rat, and the relatively normotensive strains: the Dahl Salt-Resistant (R) rat, the Milan Normotensive Strain (MNS), and the Lewis rat (LEW). These strains were used previously for linkage analysis for blood pressure (BP) in our laboratory. The results provide evidence to suggest that variations in the mitochondrial genome do not account for observed differences in blood pressure between the S and R rats. However, variants were detected among the mitochondrial genomes of the various hypertensive strains, S, SHR, and AS, and also among the normotensive strains R, MNS, and LEW. A total of 115, 114, 106, 106, and 16 variations in mtDNA were observed between the comparisons S versus LEW, S versus MNS, S versus SHR, S versus AS, and SHR versus AS, respectively. Among the 13 genes coding for proteins of the electron transport chain, 8 genes had nonsynonymous variations between S, LEW, MNS, SHR, and AS. The lack of any sequence variants between the mitochondrial genomes of S and R rats provides conclusive evidence that divergence in blood pressure between these two inbred strains is exclusively programmed through their nuclear genomes. The variations detected among the various hypertensive strains provides the basis to construct conplastic strains and further evaluate the effects of these variants on hypertension and associated phenotypes.

Blood pressure and proteinuria effects of multiple quantitative trait loci on rat chromosome 9 that differentiate the spontaneously hypertensive rat from the Dahl salt-sensitive rat.

A blood pressure (BP) quantitative trait locus (QTL) was previously located within 117 kb on rat chromosome 9 (RNO9) using hypertensive Dahl salt-sensitive and normotensive Dahl salt-resistant rats. An independent study between two hypertensive rat strains, the Dahl salt-sensitive rat and the spontaneously hypertensive rat (SHR), also detected a QTL encompassing this 117 kb region. Dahl salt-sensitive alleles in both of these studies were associated with increased BP. To map SHR alleles that decrease BP in the Dahl salt-sensitive rat, a panel of eight congenic strains introgressing SHR alleles onto the Dahl salt-sensitive genetic background were constructed and characterized. S.SHR(9)x3B, S.SHR(9)x3A and S.SHR(9)x2B, the congenic regions of which span a portion or all of the 1 logarithm of odds (LOD) interval identified by linkage analysis, did not significantly alter BP. However, S.SHR(9), S.SHR(9)x4A, S.SHR(9)x7A, S.SHR(9)x8A and S.SHR(9)x10A, the introgressed segments of which extend distal to the 1 LOD interval, significantly reduced BP. The shortest genomic segment, BP QTL1, to which this BP-lowering effect can be traced is the differential segment of S.SHR(9)x4A and S.SHR(9)x2B, to which an urinary protein excretion QTL also maps. However, the introgressed segment of S.SHR(9)x10A, located outside of this QTL1 region, represented a second BP QTL (BP QTL2) having no detectable effects on urinary protein excretion. In summary, the data suggest that there are multiple RNO9 alleles of the SHR that lower BP of the Dahl salt-sensitive rat with or without detectable effects on urinary protein excretion and that only one of these BP QTLs, QTL1, overlaps with the 117 kb BP QTL region identified using Dahl salt-sensitive and Dahl salt-resistant rats.

Increased Expression of Rififylin in A < 330 Kb Congenic Strain is Linked to Impaired Endosomal Recycling in Proximal Tubules.

Cell surface proteins are internalized into the cell through endocytosis and either degraded within lysosomes or recycled back to the plasma membrane. While perturbations in endosomal internalization are known to modulate renal function, it is not known whether similar alterations in recycling affect renal function. Rififylin is a known regulator of endocytic recycling with E3 ubiquitin protein ligase activity. In this study, using two genetically similar strains, the Dahl Salt-sensitive rat and an S.LEW congenic strain, which had allelic variants within a < 330 kb segment containing rififylin, we tested the hypothesis that alterations in endosomal recycling affect renal function. The congenic strain had 1.59-fold higher renal expression of rififylin. Transcriptome analysis indicated that components of both endocytosis and recycling were upregulated in the congenic strain. Transcription of Atp1a1 and cell surface content of the protein product of Atp1a1, the alpha subunit of Na(+)K(+)ATPase were increased in the proximal tubules from the congenic strain. Because rififylin does not directly regulate endocytosis and it is also a differentially expressed gene within the congenic segment, we reasoned that the observed alterations in the transcriptome of the congenic strain constitute a feedback response to the primary functional alteration of recycling caused by rififylin. To test this, recycling of transferrin was studied in isolated proximal tubules. Recycling was significantly delayed within isolated proximal tubules of the congenic strain, which also had a higher level of polyubiquitinated proteins and proteinuria compared with S. These data provide evidence to suggest that delayed endosomal recycling caused by excess of rififylin indirectly affects endocytosis, enhances intracellular protein polyubiquitination and contributes to proteinuria.

Refined mapping of blood pressure quantitative trait loci using congenic strains developed from two genetically hypertensive rat models.

Previously linkage and substitution mapping were conducted between the Dahl Salt-sensitive (S) rat and the Spontaneously Hypertensive Rat (SHR) to address the hypothesis that genetic contributions to blood pressure (BP) in two genetically hypertensive rat strains are different. Among the BP quantitative trait loci (QTLs) detected, two are located on chromosome 9 within large genomic segments. The goal of the current study was to develop new iterations of congenic substrains, to further resolve both of these BP QTLs on chromosome 9 as independent congenic segments. A total of 10 new congenic substrains were developed and characterized. The newly developed congenic substrains S.SHR(9)x8Ax11A and S.SHR(9)x10Ax1, with introgressed segments of 2.05 and 6.14 Mb, represented the shortest genomic segments. Both of these congenic substrains, S.SHR(9)x8Ax11A and S.SHR(9)x10Ax1 lowered BP of the S rat by 56 mm Hg (P<0.001) and 15 mm Hg (P<0.039), respectively. The BP measurements were corroborated by radiotelemetry. Urinary protein excretion was significantly lowered by SHR alleles within S.SHR(9)x10Ax1 but not by S.SHR(9)x8Ax11A. The shorter of the two congenic segments, 2.05 Mb was further characterized and found to contain a single differentially expressed protein-coding gene, Tomoregulin-2 (Tmeff2). The protein expression of Tmeff2 was higher in the S rat compared with S.SHR(9)x8Ax11A, which also had lower cardiac hypertrophy as measured by echocardiography. Tmeff2 is known to be upregulated in patients from multiple cohorts with cardiac hypertrophy. Taken together, Tmeff2 can be prioritized as a candidate gene for hypertension and associated cardiac hypertrophy in both rats and in humans.

Augmented rififylin is a risk factor linked to aberrant cardiomyocyte function, short-QT interval and hypertension.

Using congenic strains of the Dahl salt-sensitive (S) rat introgressed with genomic segments from the normotensive Lewis rat, a blood pressure quantitative trait locus was previously mapped within 104 kb on chromosome 10. The goal of the current study was to conduct extensive phenotypic studies and to further fine-map this locus. At 14 weeks of age, the blood pressure of the congenic rats fed a low-salt diet was significantly higher by 47 mm Hg (P<0.001) compared with that of the S rat. A time-course study showed that the blood pressure effect was significant from very young ages of 50 to 52 days (13 mm Hg; P<0.01). The congenic strain implanted with electrocardiography transmitters demonstrated shorter-QT intervals and increased heart rate compared with S rats (P<0.01). The average survival of the congenic strain was shorter (134 days) compared with the S rat (175 days; P<0.0007). The critical region was narrowed to <42.5 kb containing 171 variants and a single gene, rififylin. Both the mRNA and protein levels of rififylin were significantly higher in the hearts of the congenic strain. Overexpression of rififylin is known to delay endocytic recycling. Endocytic recycling of fluorescently labeled holotransferrin from cardiomyocytes of the congenic strain was slower than that of S rats (P<0.01). Frequency of cardiomyocyte beats in the congenic strain (62±9 bpm) was significantly higher than that of the S rat (24±6 bpm; P<0.001). Taken together, our study provides evidence to suggest that early perturbations in endocytic recycling caused by the overexpression of Rffl is a novel physiological mechanism potentially underlying the development of hypertension.

Epistatic genetic determinants of blood pressure and mortality in a salt-sensitive hypertension model.

Although genetic determinants protecting against the development of elevated blood pressure (BP) are well investigated, less is known regarding their impact on longevity. We concomitantly assessed genomic regions of rat chromosomes 3 and 7 (RNO3 and RNO7) carrying genetic determinants of BP without known epistasis, for their independent and combinatorial effects on BP and the presence of genetic determinants of survival using Dahl salt-sensitive (S) strains carrying congenic segments from Dahl salt-resistant (R) rats. Although congenic and bicongenic S.R strains carried independent BP quantitative trait loci within the RNO3 and RNO7 congenic regions, only the RNO3 allele(s) independently affected survival. The bicongenic S.R strain showed epistasis between R-rat RNO3 and RNO7 alleles for BP under salt-loading conditions, with less-than-additive effects observed on a 2% NaCl diet and greater-than-additive effects observed after prolonged feeding on a 4% NaCl diet. These RNO3 and RNO7 congenic region alleles had more-than-additive effects on survival. Increased survival of bicongenic compared with RNO3 congenic rats was attributable, in part, to maintaining lower BP despite chronic exposure to an increased dietary salt (4% NaCl) intake, with both strains showing delays in reaching highest BP. R-rat RNO3 alleles were also associated with superior systolic function, with the S.R bicongenic strain showing epistasis between R-rat RNO3 and RNO7 alleles leading to compensatory hypertrophy. Whether these alleles affect survival by additional actions within other BP-regulating tissues/organs remains unexplored. This is the first report of simultaneous detection of independent and epistatic loci dictating, in part, longevity in a hypertensive rat strain.

Congenic mapping of a blood pressure QTL region on rat chromosome 10 using the Dahl salt-sensitive rat with introgressed alleles from the Milan normotensive strain.

Multiple blood pressure (BP) quantitative trait loci (QTLs) are reported on rat chromosome 10 (RNO10). Of these, QTLs detected by contrasting the genome of the hypertensive Dahl salt-sensitive (S) rat with two different relatively normotensive strains, Lewis (LEW) and the Milan normotensive strain (MNS), are reported. Because the deduced QTL regions of both S vs. LEW and S vs. MNS comparisons are within large genomic segments encompassing more than 2 cM, there was a need to further localize these QTLs and determine whether the QTLs are unique to specific strain comparisons. Previously, the S.MNS QTL1 was mapped to less than 2.6 cM as a differential segment between two congenic strains. In this study, multiple congenic strains spanning the projected interval were studied. The BP effect of each strain was interpreted as the net effect of alleles introgressed within that congenic strain. The results suggest that the MNS alleles within the previously proposed differential segment (D10Rat27-D10Rat24) do not independently lower BP of the S rat. However, another congenic strain, S.MNS(10) x 9, containing introgressed MNS alleles that are outside of the previously proposed differential segment is of interest because (1) it demonstrated a BP-lowering effect, (2) it is contained within a single congenic strain and is not based on the observed effect of a differential segment, and, more importantly, (3) it overlaps with the previously identified S.LEW BP QTL region. Identification of the same QTL affecting BP in multiple rat strains will provide further support for the QTL's involvement and importance in human essential hypertension.

Closely linked non-additive blood pressure quantitative trait loci.

There is enough evidence through linkage and substitution mapping to indicate that rat chromosome 1 harbors multiple blood pressure (BP) quantitative trait loci (QTLs). Of these, BP QTL1b was previously reported from our laboratory using congenic strains derived by introgressing normotensive alleles from the LEW rat onto the genetic background of the hypertensive Dahl salt-sensitive (S) rat. The region spanned by QTL1b is quite large (20.92 Mb), thus requiring further mapping with improved resolution so as to facilitate systematic identification of the underlying genetic determinant(s). Using congenic strains containing the LEW rat chromosomal segments on the Dahl salt-sensitive (S) rat background, further iterations of congenic substrains were constructed and characterized. Collective data obtained from this new iteration of congenic substrains provided evidence for further fragmentation of QTL1b with improved resolution. At least two separate genetic determinants of blood pressure underlie QTL1b. These are within 7.40 Mb and 7.31 Mb and are known as the QTL1b1 region and the QTL1b2 region, respectively. A genetic interaction was detected between the two BP QTLs. Interestingly, five of the previously reported differentially expressed genes located within the newly mapped QTL1b1 region remained differentially expressed. The congenic strain S.LEW(D1Mco36-D1Mco101), which harbors the QTL1b1 region alone but not the QTL1b2 region, serves as a genetic tool for further dissection of the QTL1b1 region and validation of Nr2f2 as a positional candidate gene. Overall, this study represents an intermediary yet obligatory progression towards the identification of genetic elements controlling BP.

Fine-mapping and comprehensive transcript analysis reveals nonsynonymous variants within a novel 1.17 Mb blood pressure QTL region on rat chromosome 10.

The presence of blood pressure (BP) quantitative trait loci (QTL) on rat chromosome 10 has been clearly demonstrated by linkage analysis and substitution mapping. Using congenic strains containing the LEW rat chromosomal segments on the Dahl salt-sensitive (S) rat background, further iterations of congenic substrains were constructed and characterized to fine-map a chromosome 10 region (QTL1) linked to blood pressure. Comparison of seven congenic substrains refined QTL1 to a 1.17 Mb segment flanked by D10Mco88 and D10Mco89, which are located at 71,513,116 and 72,684,774 bp, respectively. The newly defined QTL1, containing 18 genes, is captured in its entirety within a single congenic substrain. A thorough transcript analysis revealed that 3 of these 18 genes, Ccl5, Ddx52, and RGD1559577, had nonsynonymous allelic variations between the S rat and the LEW rat. None of the detected transcripts within the newly defined QTL1 are implicated directly in BP control in humans or model organisms. Therefore, the present work defines a novel blood pressure QTL with three potential quantitative trait nucleotides.

Congenic interval mapping of RNO10 reveals a complex cluster of closely-linked genetic determinants of blood pressure.

Genetic dissection of the rat genome for identifying alleles that cause abnormalities in blood pressure (BP) resulted in the mapping of a significant number of BP quantitative trait loci (QTLs). In this study we mapped at least one such BP QTL on rat chromosome 10 (RNO10) as being within the introgressed segment of a S.LEW congenic strain S.LEWx12x2x3x8 spanning 1.34 Mb from 70,725,437 bp to 72,063,232 bp. BP of 3 congenic strains that span shorter segments of this region was additionally examined. Results obtained indicate that LEW alleles that comprise a 375-kb introgressed segment of the congenic strain S.LEWx12x2x3x5 (70,725,437 bp to 71,100,513 bp) increase BP. The magnitude of change in BP exhibited by the 2 strains, S.LEWx12x2x3x8 and S.LEWx12x2x3x5, is the net phenotypic effect of the underlying genetic determinants of BP. In this respect, the current data are superior to previous QTL localization of BP QTL1, which was hypothesized based on differential congenic segments. Genetic dissection using these 2 congenic strains as tools is expected to facilitate further dissection of the regions. Meanwhile, differential congenic segments were used to predict and thereby prioritize regions for candidate gene analysis. Using this approach, 2 distinct regions of 401 kb and 409 kb within S.LEWx12x2x3x8 and a 104 kb region within S.LEWx12x2x3x5 were prioritized for further consideration. Because all of these genetic elements are located within a 1.06-Mb region of RNO10, our study has revealed a remarkable insight into a genomic module comprising very closely-linked, opposing genetic determinants of BP.

Cross-talk of expression quantitative trait loci within 2 interacting blood pressure quantitative trait loci.

Genetic dissection of the S rat genome has provided strong evidence for the presence of 2 interacting blood pressure quantitative trait loci (QTLs), termed QTL1 and QTL2, on rat chromosome 5. However, the identities of the underlying interacting genetic factors remain unknown. Further experiments targeted to identify the interacting genetic factors by the substitution mapping approach alone are difficult because of the interdependency of natural recombinations to occur at the 2 QTLs. We hypothesized that the interacting genetic factors underlying these 2 QTLs may interact at the level of gene transcription and thereby represent expression QTLs or eQTLs. To detect these interacting expression QTLs, a custom QTL chip containing the annotated genes within QTL1 and QTL2 was developed and used to conduct a transcriptional profiling study of S and 2 congenic strains that retain either 1 or both of the QTLs. The results uncovered an interaction between 2 transcription factor genes, Dmrta2 and Nfia. Furthermore, the "biological signature" elicited by these 2 transcription factors was differential between the congenic strain that retained Lewis alleles at both QTL1 and QTL2 compared with the congenic strain that retained Lewis alleles at QTL1 alone. A network of transcription factors potentially affecting blood pressure could be traced, lending support to our hypothesis.