Genetic mapping of a new heart rate QTL on chromosome 8 of spontaneously hypertensive rats

Abstract Background Tachycardia is commonly observed in hypertensive patients, predominantly mediated by regulatory mechanisms integrated within the autonomic nervous system. The genetic loci and genes associated with increased heart rate in hypertension, however, have not yet been identified. Methods An F2 intercross of Spontaneously Hypertensive Rats (SHR) × Brown Norway (BN) linkage analysis of quantitative trait loci mapping was utilized to identify candidate genes associated with an increased heart rate in arterial hypertension. Results Basal heart rate in SHR was higher compared to that of normotensive BN rats (365 ± 3 vs. 314 ± 6 bpm, p < 0.05 for SHR and BN, respectively). A total genome scan identified one quantitative trait locus in a 6.78 cM interval on rat chromosome 8 (8q22–q24) that was responsible for elevated heart rate. This interval contained 241 genes, of which 65 are known genes. Conclusion Our data suggest that an influential genetic region located on the rat chromosome 8 contributes to the regulation of heart rate. Candidate genes that have previously been associated with tachycardia and/or hypertension were found within this QTL, strengthening our hypothesis that these genes are, potentially, associated with the increase in heart rate in a hypertension rat model.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Population genomic tests of models of adaptive radiation in Lake Victoria Region cichlid fish

Adaptive radiation is usually thought to be associated with speciation, but the evolution of intraspecific polymorphisms without speciation is also possible. The radiation of cichlid fish in Lake Victoria (LV) is perhaps the most impressive example of a recent rapid adaptive radiation, with 600+ very young species. Key questions about its origin remain poorly characterized, such as the importance of speciation versus polymorphism, whether species persist on evolutionary time scales, and if speciation happens more commonly in small isolated or in large connected populations. We used 320 individuals from 105 putative species from Lakes Victoria, Edward, Kivu, Albert, Nabugabo and Saka, in a radiation-wide amplified fragment length polymorphism (AFLP) genome scan to address some of these questions. We demonstrate pervasive signatures of speciation supporting the classical model of adaptive radiation associated with speciation. A positive relationship between the age of lakes and the average genomic differentiation of their species, and a significant fraction of molecular variance explained by above-species level taxonomy suggest the persistence of species on evolutionary time scales, with radiation through sequential speciation rather than a single starburst. Finally the large gene diversity retained from colonization to individual species in every radiation suggests large effective population sizes and makes speciation in small geographical isolates unlikely.

IgE, IgGa, IgGb and IgG(T) serum antibody levels in offspring of two sires affected with equine recurrent airway obstruction

Equine recurrent airway obstruction (RAO) is a chronic lower airway disease of the horse caused by hypersensitivity reactions to inhaled stable dust, including mould spores such as Aspergillus fumigatus. The goals of this study were to investigate whether total serum IgE levels and allergen-specific IgE and IgG subclasses are influenced by genetic factors and/or RAO and whether quantitative trait loci (QTL) could be identified for these parameters. The offspring of two RAO-affected sires (S1: n=56 and S2: n=65) were grouped by stallion and disease status, and total serum IgE levels and specific IgE, IgGa, IgGb and IgG(T) levels against recombinant Aspergillus fumigatus 7 (rAspf7) were measured by ELISA. A panel of 315 microsatellite markers covering the 31 equine autosomes were used to genotype the stallions and their offspring. A whole-genome scan using half-sib regression interval mapping was performed for each of the IgG and IgE subclasses. There was no significant effect of disease status or sire on total IgE levels, but there was a significant effect of gender and age. rAspf7-specific IgGa levels were significantly higher in RAO-affected than in healthy horses. The offspring of S1 had significantly higher rAspf7-specific IgGa and IgE levels than those of S2. Five QTLs were significant chromosome-wide (P<0.01). QTLs for rAspf7-specific IgGa and IgE were identified on ECA 1, for rAspf7-specific IgGa and IgGb on ECA 24 and for rAspf7 IgGa on ECA 26. These results provide evidence for effects of disease status and genetics on allergen-specific IgGa and IgE.

Sequential Quantitative Trait Locus Mapping in Experimental Crosses

The etiology of complex diseases is heterogeneous. The presence of risk alleles in one or more genetic loci affects the function of a variety of intermediate biological pathways, resulting in the overt expression of disease. Hence, there is an increasing focus on identifying the genetic basis of disease by sytematically studying phenotypic traits pertaining to the underlying biological functions. In this paper we focus on identifying genetic loci linked to quantitative phenotypic traits in experimental crosses. Such genetic mapping methods often use a one stage design by genotyping all the markers of interest on the available subjects. A genome scan based on single locus or multi-locus models is used to identify the putative loci. Since the number of quantitative trait loci (QTLs) is very likely to be small relative to the number of markers genotyped, a one-stage selective genotyping approach is commonly used to reduce the genotyping burden, whereby markers are genotyped solely on individuals with extreme trait values. This approach is powerful in the presence of a single quantitative trait locus (QTL) but may result in substantial loss of information in the presence of multiple QTLs. Here we investigate the efficiency of sequential two stage designs to identify QTLs in experimental populations. Our investigations for backcross and F2 crosses suggest that genotyping all the markers on 60% of the subjects in Stage 1 and genotyping the chromosomes significant at 20% level using additional subjects in Stage 2 and testing using all the subjects provides an efficient approach to identify the QTLs and utilizes only 70% of the genotyping burden relative to a one stage design, regardless of the heritability and genotyping density. Complex traits are a consequence of multiple QTLs conferring main effects as well as epistatic interactions. We propose a two-stage analytic approach where a single-locus genome scan is conducted in Stage 1 to identify promising chromosomes, and interactions are examined using the loci on these chromosomes in Stage 2. We examine settings under which the two-stage analytic approach provides sufficient power to detect the putative QTLs.

Mapping quantitative trait loci for osteochondrosis in fetlock and hock joints and palmar/plantar osseus fragments in fetlock joints of South German Coldblood horses

The aim of this study was to identify quantitative trait loci (QTL) for osteochondrosis (OC) and palmar/plantar osseous fragments (POF) in fetlock joints in a whole-genome scan of 219 South German Coldblood horses. Symptoms of OC and POF were checked by radiography in 117 South German Coldblood horses at a mean age of 17 months. The radiographic examination comprised the fetlock and hock joints of all limbs. The genome scan included 157 polymorphic microsatellite markers. All microsatellite markers were equally spaced over the 31 autosomes and the X chromosome, with an average distance of 17.7 cM and a mean polymorphism information content (PIC) of 63%. Sixteen chromosomes harbouring putative QTL regions were further investigated by genotyping the animals with 93 additional markers. QTL that had chromosome-wide significance by non-parametric Z-means and LOD scores were found on 10 chromosomes. This included seven QTL for fetlock OC and one QTL on ECA18 associated with hock OC and fetlock OC. Significant QTL for POF in fetlock joints were located on equine chromosomes 1, 4, 8, 12 and 18. This genome scan is an important step towards the identification of genes responsible for OC in horses.

The locus for bovine dilated cardiomyopathy maps to chromosome 18

Bovine dilated cardiomyopathy (BDCMP) is a severe and terminal disease of the heart muscle observed in Holstein-Friesian cattle over the last 30 years. There is strong evidence for an autosomal recessive mode of inheritance for BDCMP. The objective of this study was to genetically map BDCMP, with the ultimate goal of identifying the causative mutation. A whole-genome scan using 199 microsatellite markers and one SNP revealed an assignment of BDCMP to BTA18. Fine-mapping on BTA18 refined the candidate region to the MSBDCMP06-BMS2785 interval. The interval containing the BDCMP locus was confirmed by multipoint linkage analysis using the software loki. The interval is about 6.7 Mb on the bovine genome sequence (Btau 3.1). The corresponding region of HSA19 is very gene-rich and contains roughly 200 genes. Although telomeric of the marker interval, TNNI3 is a possible positional and a functional candidate for BDCMP given its involvement in a human form of dilated cardiomyopathy. Sequence analysis of TNNI3 in cattle revealed no mutation in the coding sequence, but there was a G-to-A transition in intron 6 (AJ842179:c.378+315G>A). The analysis of this SNP using the study's BDCMP pedigree did not conclusively exclude TNNI3 as a candidate gene for BDCMP. Considering the high density of genes on the homologous region of HSA19, further refinement of the interval on BTA18 containing the BDCMP locus is needed.

Heterosis for biomass-related traits in Arabidopsis investigated by quantitative trait loci analysis of the triple testcross design with recombinant inbred lines

Arabidopsis thaliana has emerged as a leading model species in plant genetics and functional genomics including research on the genetic causes of heterosis. We applied a triple testcross (TTC) design and a novel biometrical approach to identify and characterize quantitative trait loci (QTL) for heterosis of five biomass-related traits by (i) estimating the number, genomic positions, and genetic effects of heterotic QTL, (ii) characterizing their mode of gene action, and (iii) testing for presence of epistatic effects by a genomewide scan and marker x marker interactions. In total, 234 recombinant inbred lines (RILs) of Arabidopsis hybrid C24 x Col-0 were crossed to both parental lines and their F1 and analyzed with 110 single-nucleotide polymorphism (SNP) markers. QTL analyses were conducted using linear transformations Z1, Z2, and Z3 calculated from the adjusted entry means of TTC progenies. With Z1, we detected 12 QTL displaying augmented additive effects. With Z2, we mapped six QTL for augmented dominance effects. A one-dimensional genome scan with Z3 revealed two genomic regions with significantly negative dominance x additive epistatic effects. Two-way analyses of variance between marker pairs revealed nine digenic epistatic interactions: six reflecting dominance x dominance effects with variable sign and three reflecting additive x additive effects with positive sign. We conclude that heterosis for biomass-related traits in Arabidopsis has a polygenic basis with overdominance and/or epistasis being presumably the main types of gene action.

Speciation with gene-flow in eastern North American red oaks (Quercus section Lobatae)

As foundational species, oaks (Quercus : Fagaceae) support the activities of both humans and wildlife. However, many oaks in North America are declining, a crisis exacerbated by the previous disappearance of other hard mast-producing trees. In addition, the economic demands placed on this drought-tolerant group may intensify if climate change extirpates other, relatively mesophytic species. Genetic tools can help address these management challenges. To this end, we developed a suite of 27 microsatellite markers, of which 22 are derived from expressed sequence tags (ESTs). Many of these markers bear significant homology to known genes and may be able to directly assay functional genetic variation. Markers obtained from enriched microsatellite libraries, on the other hand, are typically located in heterochromatic regions and should reflect demographic processes. Considered jointly, genic and genomic microsatellites can elucidate patterns of gene-flow and natural selection, which are fundamental to both an organism's evolutionary ecology and conservation biology. To this end, we employed the developed markers in an FST-based genome scan to detect the signature of divergent selection among the red oaks (Quercus section Lobatae). Three candidate genes with putative roles in stress responses demonstrated patterns of diversity consistent with adaptation to heterogeneous selective pressures. These genes may be important in both local genetic adaptation within species and divergence among them. Next, we used an isolation-with-migration model to quantify levels of gene-flow among four red oaks species during speciation. Both speciation in allopatry and speciation with gene-flow were found to be major drivers of red oak biodiversity. Loci playing a key role in speciation are also likely to be ecologically important within species

Arachnomelia in Brown Swiss cattle maps to chromosome 5

Arachnomelia in Brown Swiss cattle is a monogenic autosomal recessive inherited congenital disorder of the skeletal system giving affected calves a spidery look (OMIA ID 000059). Over a period of 20 years 15 cases were sampled in the Swiss and Italian Brown cattle population. Pedigree data revealed that all affected individuals trace back to a single acknowledged carrier founder sire. A genome scan using 240 microsatellites spanning the 29 bovine autosomes showed homozygosity at three adjacent microsatellite markers on bovine Chr 5 in all cases. Linkage analysis confirmed the localization of the arachnomelia mutation in the region of the marker ETH10. Fine-mapping and haplotype analysis using a total of 34 markers in this region refined the critical region of the arachnomelia locus to a 7.19-Mb interval on bovine Chr 5. The disease-associated IBD haplotype was shared by 36 proven carrier animals and allows marker-assisted selection. As the corresponding human and mouse chromosome segments do not contain any clear functional candidate genes for this disorder, the mutation causing arachnomelia in the Brown Swiss cattle might help to identify an unknown gene in bone development.

Genetic mapping of the belt pattern in Brown Swiss cattle to BTA3

The white belt pattern of Brown Swiss cattle is characterized by a lack of melanocytes in a stretch of skin around the midsection. This pattern is of variable width and sometimes the belt does not fully circle the body. To identify the gene responsible for this colour variation, we performed linkage mapping of the belted locus using six segregating half-sib families including 104 informative meioses for the belted character. The pedigree confirmed a monogenic autosomal dominant inheritance of the belted phenotype in Brown Swiss cattle. We performed a genome scan using 186 microsatellite markers in a subset of 88 animals of the six families. Linkage with the belt phenotype was detected at the telomeric region of BTA3. Fine-mapping and haplotype analysis using 19 additional markers in this region refined the critical region of the belted locus to a 922-kb interval on BTA3. As the corresponding human and mouse chromosome segments contain no obvious candidate gene for this coat colour trait, the mutation causing the belt pattern in the Brown Swiss cattle might help to identify an unknown gene influencing skin pigmentation.

Gene Discovery in Nonsyndromic Cleft Lip with or without Cleft Palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP), a common, complex orofacial birth defect that affects approximately 4,000 newborns each year in the United States, is caused by both genetic and environmental factors. Orofacial clefts affect the mouth and nose, causing severe deformity of the face, which require medical, dental and speech therapies. Despite having substantial genetic liability, less than 25% of the genetic contribute to NSCLP has been identified. The studies described in this thesis were performed to identify genes that contribute to NSCLP and to demonstrate the role of these genes in normal craniofacial development. Using genome scan and candidate gene approaches, novel associations with NSCLP were identified. These include MYH9 (7 SNPs, 0.009≤p<0.05), Wnt3A (4 SNPs, 0.001≤p≤0.005), Wnt11 (2 SNPs, 0.001≤p≤0.01) and CRISPLD2 (4 SNPs, 0.001≤p<0.05). The most interesting findings were for CRISPLD2. This gene is expressed in the fused mouse palate at E17.5. In zebrafish, crispld2 localized to the craniofacial region by one day post fertilization. Morpholino knockdown of crispld2 resulted in a lower survival rates and altered neural crest cell (NCC) clustering. Because NCCs form the tissues that populate the craniofacies, this NCC abnormality resulted in cartilage abnormalities of the jaw including fewer ceratobranchial cartilages forming the lower jaw (three pairs compared to five) and broader craniofacies compared to wild-type zebrafish. These findings suggest that the CRISPLD2 gene plays an important role in normal craniofacial development and perturbation of this gene in humans contributes to orofacial clefting. Overall, these results are important because they contribute to our understanding of normal craniofacial development and orofacial clefting etiology, information that can be used to develop better methods to diagnose, counsel and potentially treat NSCLP patients.

The genetic basis of thoracic aortic aneurysms and dissections: Genetic heterogeneity and mapping of TAAD1 and TAAD2 loci

Thoracic aortic aneurysms leading to aortic dissections (TAAD) are a major cause of morbidity and mortality in the United States. TAAD is a complication of some known genetic disorders, such as Marfan syndrome and Turner syndrome, but the majority of familial cases are not due to a known genetic syndrome. Previous studies by our group have established that nonsyndromic, familial TAAD is inherited in an autosomal dominant manner with decreased penetrance and variable expression. Using one large family with multiple members with TAAD for the genome wide scan, a major locus for familial TAAD was mapped to 5q13–14 (TAAD1). Nine out of 15 families studied were linked to this locus, establishing that TAAD1 was a major locus, and that there was genetic heterogeneity for the condition. Mapping of TAAD2 locus was accomplished using a single large family with multiple members with TAAD not linked to known loci of aneurysm formation. This established a second novel locus for familial TAAD on 3p24–25 (LOD score of 4.3), termed the TAAD2 locus. Two putative loci with suggestive LOD scores were mapped on 4q and 12q through a genome scan carried out using three families. TAAD phenotype in 12 families did not segregate with known loci, indicating further genetic heterogeneity. An STS-tagged BAC based contig was constructed for 7.8Mb and 25Mb critical interval of TAAD1 and TAAD2 respectively and characterized to identify the defective gene. The hypothesis that the defective genes responsible for the TAAD1 and TAAD2 encoded extracellular matrix (ECM) proteins, the major components of the elastic fiber system in the aortic media was tested. Four genes encoding ECM proteins, versican, thrombospondin-3, CRTL1, on TAAD1 and FBLN2 at TAAD2 were sequenced, but no disease-causing mutations were identified. Studies to identify the defective gene are initiated through the positional candidate gene approach using combination of bioinformatics and expression studies. The identification of the TAAD susceptibility genes will allow for presymptomatic diagnosis of individuals at risk for this life threatening disease. The identification of the molecular defects that contribute to TAAD will also further our understanding of the proteins that provide structural integrity to the aortic wall. ^

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

A human model for multigenic inheritance: Phenotypic expression in Hirschsprung disease requires both the RET gene and a new 9q31 locus

Reduced penetrance in genetic disorders may be either dependent or independent of the genetic background of gene carriers. Hirschsprung disease (HSCR) demonstrates a complex pattern of inheritance with ≈50% of familial cases being heterozygous for mutations in the receptor tyrosine kinase RET. Even when identified, the penetrance of RET mutations is only 50–70%, gender-dependent, and varies with the extent of aganglionosis. We searched for additional susceptibility genes which, in conjunction with RET, lead to phenotypic expression by studying 12 multiplex HSCR families. Haplotype analysis and extensive mutation screening demonstrated three types of families: six families harboring severe RET mutations (group I); and the six remaining families, five of which are RET-linked families with no sequence alterations and one RET-unlinked family (group II). Although the presence of RET mutations in group I families is sufficient to explain HSCR inheritance, a genome scan reveals a new susceptibility locus on 9q31 exclusively in group II families. As such, the gene at 9q31 is a modifier of HSCR penetrance. These observations imply that identification of new susceptibility factors in a complex disease may depend on classification of families by mutational type at known susceptibility genes.