|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Max Delbrück Centrum for Molecular Medicine, Berlin, Germany; and
2 Inserm, U698, Hôpital Bichat, Paris; University Paris 7, France
| ABSTRACT |
|---|
|
|
|---|
internal elastic lamina; abdominal aorta; elastin; ductus arteriosus; quantitative trait loci
| INTRODUCTION |
|---|
|
|
|---|
The adult BN also presents a lower elastin content and a proportionally higher collagen content in the aorta, both in the thoracic segment, which is unaffected by the above-mentioned RIEL, and in the abdominal segment, presenting numerous ruptures (31). The elastin deficit is at least in part due to a deficiency in tropoelastin synthesis in the BN rat, compared with LE and LOU, during the period of rapid growth (31, 32).
The genetic basis for the spontaneous RIEL in abdominal aorta was recently studied by another group, in an F2 population originating from a cross between BN and New Zealand genetically hypertensive (GH) rat strains (16, 11). This latter strain exhibits very low levels of aortic RIEL, similarly to the LE and LOU strains. The authors concluded from their phenotypic data that the character for RIEL was autosomal, dominant, and possibly a single gene effect (16). Subsequent linkage studies on the same population (11) identified two major quantitative trait loci (QTLs) for aortic RIEL, on chromosomes 5 and 10, exhibiting respectively semidominant and dominant modes of inheritance, with the BN alleles being associated with increased lesion number. These results are at variance with our phenotypic data in BN x LE and BN x LOU crosses because, in our case, no F1 rats and very few F2 rats presented significant numbers of RIEL. The BN, LE, and LOU strains are all normotensive, and adult male rats of these strains show no significant difference in systolic BP (2, 32).
In view of the divergence in the phenotypic data for RIEL incidence obtained by crossing BN rats with normotensive LE or LOU and with hypertensive GH rats, we have undertaken a genome-wide scan to locate QTLs for this character in a large cohort of rats originating from a BN x LOU backcross (BC). We chose to use BC rats [(BN x LOU) x BN], a population more suited to the study of RIEL than the F2 population, because of the previously observed recessive character of the phenotype in F2 rats. We included in our study the quantification of aortic elastin and collagen, the major aortic extracellular matrix (ECM) proteins, and of renal artery RIEL, a phenotype whose genetic basis has not yet been studied, and whose relation with aortic RIEL is unclear.
In addition to the above-mentioned phenotypes, the BN rat also presents another character of pathophysiological interest, namely persistence of the ductus arteriosus (PDA). PDA is quite frequently observed in the BN rat and never in the LOU or other strains (unpublished data). We thus included this phenotype in the present linkage analysis, which represents the first quantitative genetic investigation for PDA. Our data also enables us to determine whether RIEL, decreased aortic elastin-collagen ratio, and PDA are all due to some common underlying genetic variation, or whether they are independent phenotypes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Determination of Quantitative Traits
We phenotyped 41 F1 rats (20 male and 21 females) and 209 BC rats (104 males and 105 females) at 18 wk of age for RIEL, ductus arteriosus (DA) status, and thoracic aortic (TA) elastin and collagen content as described below. A further 36 BC rats were phenotyped later for fine mapping of QTL intervals and were integrated in the final linkage analysis presented here.
RIEL
RIEL were quantified on "en face" arterial preparations as previously described (2). For the AA, renal arteries, common IAs, and proximal 5 mm of the caudal artery, the total number of RIEL was recorded. Then, for AA and IAs only, each rupture was graded on a semiquantitative scale according to its size, and a score, taking into account the severity of the phenomenon, was allotted to each rat.
Scoring for PDA
Immediately after thawing the thoracic aortas to proceed to biochemical analysis, we evaluated the persistence and state of closure of the DA under a dissecting microscope (magnification x12.5) and graded it by a semiquantitative score, taking into account its various features.
Composition of the Thoracic Aorta
Biochemical analysis of the descending thoracic aorta was performed as previously described (2, 32). Elastin, collagen, and cell proteins were all expressed both as a percentage of aortic dry weight and as mg/cm of aorta.
Genetic Studies
Microsatellite markers.
We selected 175 microsatellite markers that exhibited allelic variation between the parental strains and that cover the 20 autosomal chromosomes of the rat genome. The primer pair information of these markers was obtained from RGD (http://www.rgd.mcw.edu/). Fine mapping with 18 additional markers was performed on 36 extra BC rats. All markers are listed in Suppl. Table S1. (The online version of this article contains supplemental material.)
Genotyping.
DNA extraction followed by fluorescence-based semiautomated genotyping was performed (see online supplement). The genetic map of the 193 microsatellite markers was established using MAPMAKER/EXP V3.0 (21).
Linkage analysis.
Linkage analysis was conducted using Mapmaker software (21) using all 209 BC progeny and separately for the subsets of 105 females and 104 males. Both parametric and nonparametric linkage analyses were performed using Mapmaker /QTL V3.0 (20) and V1.9 (19) and are reported as logarithm of the odds ratio (LOD) or Z-scores respectively.
Microarray analysis and validation of microarray gene expression data.
Microarray analysis was performed on AA from 6-wk-old male BN and LOU rats (n = 5 per strain) as described previously (14) using Affymetrix RAE 230A GeneChip expression arrays. This age was chosen, as RIEL have not yet started to form in the BN AA, which occurs from
7 wk onward. Quantitative real-time PCR, performed on abdominal and TA segments as described previously (14), was used to validate myosin heavy chain 10, nonmuscle (Myh10) expression differences observed between BN and LOU using microarrays. All microarray expression data has been submitted to ArrayExpress (accession number pending).
Candidate gene analysis by sequencing.
Mutation screening was performed on candidate genes by sequencing cDNA or genomic DNA (see online supplement). All sequencing data have been submitted to GenBank (accession number pending), and all single nucleotide polymorphism (SNPs) identified to db-SNP.
Production of a congenic rat line.
To prove the involvement of the chromosome 5 locus in the RIEL phenotype and to provide a tool for future studies, a line of congenic rats was produced by introgression of the segment of interest on chromosome 5 (D5rat 59-D5rat 131) from the BN donor onto the LOU recipient genome after 8 BCs and 1 intercross to obtain homozygosity as previously described (13).
For more details concerning materials and methods, see online supplement.
| RESULTS |
|---|
|
|
|---|
|
|
|
Large RIEL in AA and IAs showed suggestive linkage to chromosomes 2, 8, and 10 (Suppl. Fig. S1, E and F). The chromosome 10 QTL, between markers D10Rat96 and D10Rat27, (Z-score 3.3 at peak marker D10Rat64) overlapped with that previously described in rats (11) and is syntenic to the recently published locus in humans linked to a pathophysiological entity associating TA aneurysm and aortic dissection with PDA (18). Possible candidate genes in this QTL were the glycoprotein Sparc and three myosin heavy chain genes: Myh11, implicated in the above human syndrome (40), Myh3, and Myh10. Mutation screening was negative for these genes, but both Sparc and Myh10 were differentially expressed on the chip (Suppl. Table S4). Quantitative PCR confirmed a 2.7-fold downregulation of Myh10 in the AA of BN vs. LOU, making it a strong candidate for RIEL in AA and IAs. We thus sequenced, in addition to the exonic region of Myh10, regions 5 kb upstream and 1 kb downstream to screen for variations in the regulatory regions, but no mutation was found. Only one synonymous SNP in exon 17 was identified (Suppl. Table S4).
Renal artery.
Renal artery RIEL showed a slightly different distribution in the various rat populations from RIEL in the AA and IAs, as a few were present in male LOU and some BC rats showed more than BN rats (Suppl. Tables S2 and S3, Fig. 1B). A renal artery of a BC rat is illustrated in Fig. 2D.
Interestingly, unlike RIEL in the AA and IAs, those in the renal arteries did not map to chromosome 5, but to chromosome 3, with a peak LOD score of 4.0 (Fig. 1B, Table 1, Suppl. Fig. S1i). Renal artery RIEL were more numerous in the heterozygous than in the homozygous BC progeny (P = 0.0005), suggesting an effect of the LOU allele over this phenotype. Fbn1 was the only potential candidate gene in this QTL, but it was not significantly differentially expressed and sequencing revealed no variations leading to amino acid changes (Suppl. Table S4).
PDA
In adult rats of the BN strain it is not rare to observe a PDA (Suppl. Table S2, Fig. 1, C and D). The most affected BN rats exhibit a diverticulum, wide open on the aortic side and nearly closed on the pulmonary side (Fig. 3, AC). This observed anatomy suggests that the aorto-pulmonary shunt effect is absent or minimal, explaining the lack of overt symptoms in affected BN rats. The mildness of the phenotype implies that complete DA closure is probably polygenic and only some of the implicated genes are affected here. A continuous variation of this trait was observed within the BN strain, ranging from a grossly visible diverticulum (as in Fig. 3, AC; 18% of the total BN population), via various intermediate stages of involution, to the presence of a normal fibrous ligamentum arteriosus (14% of the total BN population), as consistently observed in the LOU strain (Fig. 3D). These observations are reflected by the significantly different scores between BN and LOU (Suppl. Table S2). BC rats showed a wide range of scores, as shown in Fig. 3, EG, with only 3.7% overall showing a macroscopically visible diverticulum as shown in Fig. 3E, and 70% with scores lying close to those observed in LOU rats (Suppl. Table S3).
|
TA Composition
LOU rats had higher TA elastin content (% of aortic dry weight) and elastin-collagen ratio than BN rats, and the F1 population showed intermediate values, as previously described (32). In the F1 and BC cohorts, females had significantly higher values than males, and the same tendency was observed in the parental strains (Suppl. Tables S2 and S3).
TA elastin content (%) (Fig. 1E) and elastin-collagen ratio (Fig. 1F) were significantly linked to chromosome 2 (Table 1; Suppl. Fig. S1, ln) although elastin-collagen linkage was observed in male progeny only and the pattern of the LOD score curve suggested two independent loci (Fig. 1F), with highest LOD scores at D2Rat205 and D2Rat51 respectively (Table 1). These two loci overlap with the QTL linked to TA % elastin. The BN allele was associated with lower values at all these three loci (Fig. 1, E and F).
TA dry weight (mg/cm) was significantly linked to chromosomes 7 and 19, showing respective LOD scores of 4.0 and 3.3 (Suppl. Table S6, Suppl. Fig. S1K). Both elastin and collagen contents, when expressed in mg/cm, also mapped to these loci (Suppl. Fig. S1, M and O), as these parameters are linked to aortic dry weight. Suggestive linkage was also observed for TA elastin content (expressed as both % and mg/cm) to chromosome 17 (Suppl. Fig. S1, M and N). Other suggestive loci were found for the various parameters expressing TA elastin and collagen contents on chromosomes 1 and 3 (Suppl. Fig. S1, KP).
The QTL on chromosome 2 for % elastin and elastin-collagen span a relatively large interval and harbor numerous known and predicted genes. We chose to sequence Kcnd3, due to a possible role of K+ on elastin synthesis. We have detected eight synonymous SNPs in this gene (Suppl. Table S4).
Results of linkage analysis for other weight-related phenotypes are summarized in Suppl. Table S6.
No correlation was found between any of the major phenotypes, i.e., RIEL, PDA, and aortic elastin content, nor did we identify any epistatic interaction between the different loci. No significant lineage effects were detected in the different subgroups (from male Y-BN and Y-LOU matings, respectively), which we tested for at all significant QTL regions.
Congenic Line LOU.BN.D5Rat59-D5Rat131
The details of the congenic strain LOU.BN.D5Rat59-D5Rat131, with the chromosomal segment of interest from BN introgressed into the LOU genetic background, are shown in Fig. 4. The number of RIEL in this congenic line compared with parental BN, LOU, F1, and BC (homozygous BN/BN at the locus of interest) are shown in Fig. 4B. The congenic rats present
30% of the BN value for this quantitative trait, confirming our linkage data and providing us with a useful tool for future studies.
|
| DISCUSSION |
|---|
|
|
|---|
The first major observation in this study is that the principal phenotypes are not correlated in the BC population and map to separate QTLs, indicating that they are each under independent genetic control and so do not represent a "syndrome."
The major QTL on chromosome 5, linked to RIEL in AAs and IAs, was responsible for a large part (45%) of the phenotype and is identical to that described by Harris et al. (11) in a BN-GH cross. However, these authors also found a significant QTL on chromosome 10, the two loci together determining only 33.5% of the phenotype. Among several suggestive loci we found for this phenotype, our chromosome 10 locus overlapped with that previously described (11) and thus can be regarded as a confirmed QTL.
The major chromosome 5 locus contained no obvious potential candidate genes, such as lysyl oxydase or a known component of elastic fibers, and no genes in this region were differentially expressed between BN and LOU in the microarray analysis. However, considering that RIEL is likely to be due to some other underlying ECM anomaly or a dysfunction of smooth muscle or endothelial cells, we chose five genes for mutation analysis. Sulphatase-FP (Sulf1) (25), matrix metalloproteinase 16 (Mmp16), and Yamaguchi Sarcoma viral (v-yes-1) oncogene homolog (Lyn), via its participation in the regulation of uPA (11), may be involved in degradation of components of the ECM, retinol dehydrogenase 10 (all-trans) (Rdh10) in arterial smooth muscle signaling (1, 6), and potassium voltage-gated channel, Shab-related family, member 2 (Kcnb2) may influence elastin synthesis (36). However, no significant mutations were found in the coding regions of any of these genes.
The chromosome 10 locus harbors several myosin heavy chain genes. Myh11, a smooth muscle form (SM-MHC) was an interesting candidate in view of its recently described implication in a human syndrome associating PDA and aneurysms (40). Moreover, a transgenic mouse strain lacking SMC-MHC gene expression presents delayed closure of the ductus arteriosus (24). Myh3 encodes for an embryonic myosin form and Myh10 for the nonmuscle myosin heavy chain IIB (NMHC IIB), which is present in all mammalian cells and is involved in physiological functions such as cell motility, morphology and cytokinesis (34, 35). In mice this gene is required for normal cardiac development (35, 37). However, sequencing these three genes did not provide evidence for their implication in RIEL. In contrast, Myh10, but not Myh11, was downregulated in BN AA compared with LOU. This finding merits further investigation by a more integrated approach to minimize the limitations of genetic mapping and differential expression analysis, as previously discussed for hypertension (39).
To validate our findings using congenic approaches we are currently establishing several congenic lines for some of the loci identified. Here we represent a congenic line for the chromosome 5 locus, in which we were able to capture part of the phenotype of interest. This represents an essential tool to further the search for genes involved in this striking vascular phenotype of spontaneous RIEL. Indeed, looking at functionally relevant genes has shown its limitations here, and using congenics and derivation of sublines will provide an unbiased view. Although the above-mentioned myosin genes merit further investigation, it is possible that some developmental or remodeling genes or an as yet unknown, minor constituent of elastic fibers may be implicated. Indeed, the steps involved in IEL formation during arterial growth and development are poorly understood, and the role, if any, of the endothelial cells in this process is unclear. Discovery of the gene(s) present in this QTL with such a major influence on RIEL may provide insight into fundamental mechanisms involved in development of the vascular elastic fiber network.
The tight correlation that we observed between RIEL in the AA, IAs, and the proximal caudal artery within the BC population is not surprising, as these are all closely related elastic arteries and the common iliac and caudal arteries take origin in the affected AA. However, renal artery RIEL were much less significantly correlated with those in AA and IAs in males and did not correlate at all in females (results not shown). Indeed the upper AA, where the renal arteries take origin, is unaffected by RIEL in the BN rat, as this phenotype occurs only in the subrenal aorta (unpublished observations). Moreover, we found a distinct QTL for renal artery RIEL on chromosome 3, which harbors the gene for fibrillin-1 (Fbn1), a potential candidate gene for elastic fiber anomalies in view of its implication in the Marfan syndrome (22). However, no mutation was found.
Surprisingly, the LOU allele was associated with increased RIEL in the renal artery. We have previously shown that renal artery RIEL, unlike those in AA and IAs, occur in all rat strains, to varying degrees (5, 26, 27), and they may be influenced by many local factors, such as blood pressure and flow and thus, indirectly, kidney size (26, 27). We did find QTLs for kidney size on chromosomes 3 and 7 (Suppl. Table S6), but the one on chromosome 3 did not localize to the same area as that for renal artery RIEL.
It thus appears that the presence of RIEL is a complex phenotype, as superposed on a genetically determined susceptibility, which may or may not be the same throughout the arterial tree, local factors undoubtedly play a role in inducing these ruptures. The fact that the BN is the only rat strain to present the AA/IA/proximal caudal artery phenotype (2), whereas all strains present RIEL in their renal arteries and in the more distal, muscular part of the caudal artery, albeit to different degrees (5, 26), is intriguing. Thus, either the two phenotypes are distinct, as suggested by the present linkage analysis, or local factors so greatly influence the phenotype that in the adult the effect of the original determining gene may be masked in the renal artery. Indeed, RIEL, once formed in rat arteries are never completely repaired since, despite rapid replacement of the cellular elements, the gap in the IEL remains throughout life. The effect of any local factors on numbers of RIEL is thus cumulative and would be proportionately greater the older the population studied.
One of the aims of this study was to investigate whether the aortic elastin deficit in the BN rat could be the cause of RIEL formation. Previous studies have indirectly suggested that this was not the case (2, 9), and here we found no correlation between these two phenotypes in the BC population in either sex, suggesting no causal relationship. Indeed, aortic RIEL correlated with none of the aortic biochemical parameters we measured in this study. Furthermore, RIEL and aortic elastin content are controlled by genes on different chromosomes. The QTLs we found on chromosome 2, for aortic elastin content and elastin-collagen ratio, overlap with those previously described (9) and confirm that the BN allele is linked to decreased elastin content and that many different genes control this phenotype. None of our QTLs for elastin content mapped to the position of Eln, confirming that the elastin gene itself controls a very minor part of the variation in aortic elastin content (32).
Identification of genes controlling ductal closure in a quantitative manner has never been attempted before. Studies in both dog (28) and human (10) suggest a genetic component, but the exact mode of transmission and the number of genes involved remain unknown. In the rat it appears to be a rare phenotype, which we only discovered because of high frequency sampling of BN thoracic aortas. Recently, another group has confirmed our observations (3). We found two significant QTLs for PDA, on chromosomes 8 and 9, with the BN allele giving the highest scores. Among genes present in these QTLs are a large number of myosin genes, which is of interest in view of the recently reported involvement of Myh11 in a human syndrome with PDA (40). However, we found no mutations with possible functional repercussions in these various myosin genes, although several synonymous SNPs were found. Three collagen genes (Col5a3 on Chr8, Col3a1, and Col5a2 on Chr9) were also present in these loci, which merited further exploration as we found a weak correlation (P = 0.047) between this phenotype and aortic collagen content in males. However, this is may be a secondary phenomenon, as sequencing of two of these genes did not show any significant mutations in the BN rat.
A recent study in an Iranian population has reported a recessive component to PDA, implicating a single locus on chromosome 12q24 (23). Sequencing of selected candidate genes in this locus (nNOS, paxillin, and PTPN11) did not provide evidence of their implication in the phenotype. In our study, one of our three suggestive loci was in the syntenic region on the rat genome (peak marker D12Rat80, LOD score 3.5). This QTL also contained a myosin gene (Mlc2), but no mutation was found. Ductal closure is a complex process involving smooth muscle cell contraction, synthesis of ECM molecules, neointima formation etc., implicating a multitude of possible genes. Further studies are thus required, but, in any event, the finding of a rat strain with incomplete ductal closure opens up new perspectives into the study of genes intervening in the control of this process, a crucial event at birth for subsequent normal cardiovascular function.
In conclusion, we have shown here that the BN rat presents several arterial phenotypes that are of pathophysiological and of fundamental interest. Here we provide considerably detailed data concerning the various phenotypes in the same population. We conclude that these phenotypes appear to be independent, controlled by different genes, and in this respect it is surprising to find them united in one strain, which, moreover, appears to have a normal lifespan. The BN rat has also been shown to present other vascular abnormalities not studied here, such as poor survival in hypertension, related to susceptibility to cerebrovascular hemorrhage (5, 7), impaired myogenic autoregulation in the kidney (38), and a defective flow-induced arterial wall remodeling (8, 15). In addition, the BN rat also presents hydronephrosis, already reported by others (29, 33), for which we have found linkage to another distinct genetic locus to be reported elsewhere. Obviously further studies will be necessary to unravel the mechanisms underlying this accumulation of subclinical abnormalities of a predominantly vascular nature in one strain, a strain that appears to be a genetic outlier compared with most other commonly used laboratory rats (4). Lastly, these data should be of particular interest to the scientific community as the BN rat has been widely used as a "reference" strain in genetic studies and was chosen to sequence the rat genome (30).
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
* L. Kota and M. Osborne-Pellegrin contributed equally to this work. ![]()
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |