Score tests of genetic association in the presence of linkage based on the additive genetic gamma frailty model

Nuclear families with multiple affected sibs are often collected for genetic linkage analysis of complex diseases. Once linkage evidence is established, dense markers are often typed in the linked region for genetic association analysis based on linkage disequilibrium (LD). Detection of association in the presence of linkage localizes disease genes more accurately than the methods that rely on linkage alone. However, test of association due to LD in the linked region needs to account for dependency of the allele transmissions to different sibs within a family. In this paper, we define a joint model for genetic linkage and association and derive the corresponding joint survival function of age of onset for the sibs within a sibship. The joint survival function is a function of both the inheritance vector and the genotypes at the candidate marker locus. Based on this joint survival function, we derive score tests for genetic association. The proposed methods utilize the phenotype data of all the sibs and have the advantages of family-based designs which can avoid the potential spurious association caused by population admixture. In addition, the methods can account for variable age of onset or age at censoring and possible covariate effects, and therefore provide important tools for modelling disease heterogeneity. Simulation studies and application to the data sets from the 12th Genetic Analysis Workshop indicate that the proposed methods have correct type 1 error rates and increased power over other existing methods for testing allelic association.     

Family-based association analysis with tightly linked markers

Refining genomic regions which have been identified by linkage analysis to contain a disease susceptibility locus has proven to be a challenging task. Detecting association between the disease and a genetic marker can significantly narrow down the candidate region. Since an adequate sample of families is already available from the genome scan, family-based association tests may be used to search for association. The use of haplotypes consisting of tightly linked markers can be more powerful for detecting association than the use of individual markers. An extension of the transmission/disequilibrium test to allow the simultaneous analysis of more than one marker locus is complicated by ambiguity of phase in some families of the sample. The present paper shows that a recently proposed method for the analysis of nuclear families with a single affected child can be viewed as a special application of a more general principle. This observation justifies several modifications, potentially increasing the power, as well as an extension of the method to allow the analysis of general nuclear families. Finally, the problem of missing parental genotypes is discussed.     

Gene mapping by linkage and association analysis

Genetic analysis is used to map genes, including disease loci, to positions within the human genome. Linkage analysis depends on the co-segregation of a gene (locus) and a phenotype through a pedigree, while association analysis, or linkage disequilibrium mapping, depends on measuring deviation from the random occurrence of alleles in a haplotype in unrelated individuals or nuclear families. Complex computer programs may be used in both forms of analysis. In recent years most interest has focused on identifying genes involved in common, multifactorial diseases. Here I review some current and developing techniques of genetic analysis and give references to where further information can be obtained.     

Evidence for asthma susceptibility genes on chromosome 11 in an African-American population

Initial genome-wide scan data provided suggestive evidence for linkage of the asthma phenotype in African-American (AA), but not Caucasian, families to chromosome 11q markers (peak at D11S1985; LOD=2). To refine this region, mapping analysis of 91 AA families (51 multiplex families and 40 asthmatic case-parent trios) was performed with an additional 17 markers flanking the initial peak linkage marker. Multipoint analyses of the 51 multiplex families yielded significant evidence of linkage with a peak non-parametric linkage score of 4.38 at marker D11S1337 (map position 68.6 cM). Furthermore, family-based association and transmission disequilibrium tests conducted on all 91 families showed significant evidence of linkage in the presence of disequilibrium for several individual markers in this region. A putative susceptibility locus was estimated to be at map position 70.8 cM with a confidence interval spanning the linkage peak. Evidence from both linkage and association analyses suggest that this region of chromosome 11 contains one or more susceptibility genes for asthma in these AA families.     

A comparative genetic map of the turkey genome

Genetic markers (microsatellites and SNPs) were used to create and compare maps of the turkey and chicken genomes. A physical map of the chicken genome was built by comparing sequences of turkey markers with the chicken whole-genome sequence by BLAST analysis. A genetic linkage map of the turkey genome (Meleagris gallopavo) was developed by segregation analysis of genetic markers within the University of Minnesota/Nicholas Turkey Breeding Farms (UMN/NTBF) resource population. This linkage map of the turkey genome includes 314 loci arranged into 29 linkage groups. An additional 40 markers are tentatively placed within linkage groups based on two-point LOD scores and 16 markers remain unlinked. Total map distance contained within linkage groups is 2,011 cM with the longest linkage group (47 loci) measuring 413.3 cM. Average marker interval over the 29 linkage groups was 6.4 cM. All but one turkey linkage group could be aligned with the physical map of the chicken genome. The present genetic map of the turkey provides a comparative framework for future genomic studies.     

Linkage maps of microsatellite DNA markers for the Pacific oyster Crassostrea gigas

We constructed male and female consensus linkage maps for the Pacific oyster Crassostrea gigas, using a total of 102 microsatellite DNA markers typed in 11-day-old larvae from three families. We identified 11 and 12 linkage groups in the male and female consensus maps, respectively. Alignment of these separate maps, however, suggests 10 linkage groups, which agrees with the haploid chromosome number. The male linkage map comprises 88 loci and spans 616.1 cM, while the female map comprises 86 loci and spans 770.5 cM. The male and the female maps share 74 loci; 2 markers remain unlinked. The estimated coverages for the consensus linkage maps are 79% for the male and 70-75% for the female, on the basis of two estimates of genome length. Ninety-five percent of the genome is expected to lie within 16 and 21 cM of markers on the male and female maps, respectively, while 95% of simulated minimum distances to the male and female maps are within 10.1 and 13.6 cM, respectively. Females have significantly more recombination than males, across 118 pairs of linked markers in common to the parents of the three families. Significant differences in recombination and orders of markers are also evident among same-sex parents of different families as well as sibling parents of opposite sex. These observations suggest that polymorphism for chromosomal rearrangements may exist in natural populations, which could have profound implications for interpreting the evolutionary genetics of the oyster. These are the first linkage maps for a bivalve mollusc that use microsatellite DNA markers, which should enable them to be transferred to other families and to be useful for further genetic analyses such as QTL mapping.     

Analysis of porcine MUC4 gene as a candidate gene for prolificacy QTL on SSC13 in an Iberian × Meishan F2 population

Background

Canine atopic dermatitis (AD) is a common, heritable, chronic allergic skin condition prevalent in the West Highland White Terrier (WHWT). In canine AD, environmental allergens trigger an inflammatory response causing visible skin lesions and chronic pruritus that can lead to secondary bacterial and yeast infections. The disorder shares many of the clinical and histopathological characteristics of human AD and represents an animal model of this disorder that could be used to further elucidate genetic causes of human AD. Microsatellite markers genotyped in families of WHWTs affected with AD were used to perform a genome-wide linkage study in order to isolate chromosomal regions associated with the disorder.

Results

Blood samples and health questionnaires were collected from 108 WHWTs spanning three families. A linkage simulation using these 108 dogs showed high power to detect a highly penetrant mutation. Ninety WHWTs were genotyped using markers from the Minimal Screening Set 2 (MSS-2). Two hundred and fifty six markers were informative and were used for linkage analysis. Using a LOD score of 2.7 as a significance threshold, no chromosomal regions were identified with significant linkage to AD. LOD scores greater than 1.0 were located in a 56 cM region of chromosome 7.

Conclusions

The study was unable to detect any chromosomal regions significantly linked to canine AD. This could be a result of factors such as environmental modification of phenotype, incorrect assignment of phenotype, a mutation of low penetrance, or incomplete genome coverage. A genome-wide SNP association study in a larger cohort of WHWTs may prove more successful by providing higher density coverage and higher statistical power.     

Fine-mapping chromosome 20 in 230 systemic lupus erythematosus sib pair and multiplex families: evidence for genetic epistasis with chromosome 16q12

The presence of systemic lupus erythematosus (SLE) susceptibility genes on chromosome 20 is suggested by the observation of genetic linkage in several independent SLE family collections. To further localize the genetic effects, we typed 59 microsatellites in the two best regions, as defined by genome screens. Genotypes were analyzed for statistical linkage and/or association with SLE, by use of a combination of nonparametric linkage methods, family-based tests of association (transmission/disequilibrium and pedigree disequilibrium tests), and haplotype-sharing statistics (haplotype runs test), in a set of 230 SLE pedigrees. Maximal evidence for linkage to SLE was to 20p12 (LOD = 2.84) and 20q13.1 (LOD = 1.64) in the white pedigrees. Subsetting families on the basis of evidence for linkage to 16q12 significantly improved the LOD scores at both chromosome 20 locations (20p12 LOD = 5.06 and 20q13 LOD = 3.65), consistent with epistasis. We then typed 162 single-nucleotide polymorphism markers across a 1.3-Mb candidate region on 20q13.1 and identified several SNPs that demonstrated significant evidence for association. These data provide additional support for linkage and association to 20p12 and 20q13.1 in SLE and further refine the intervals of interest. These data further suggest the possibility of epistatic relationships among loci within the 20q12, 20q13, and 16q12 regions in SLE families.     

Mapping Approaches to Gene Identification in Humans

Although a number of human genes that cause disease have been traced through the defective product, most genetic defects are recognized only by phenotype. When the biochemical defect is unknown, a gene can be located only through molecular approaches based on coinheritance (genetic linkage) of the disease phenotype with a particular allele of a polymorphic DNA marker that has already been mapped to a specific chromosomal region. Linkage studies in affected families have already localized genes for several important diseases, including cystic fibrosis. Finding a genetic linkage in families in which a disease segregates requires that the human genetic map have a large number of polymorphic markers; when the map is dense enough, any disease gene can be located by linkage to a known marker. Many DNA segments with a high degree of polymorphism are being found and mapped as markers in normal reference pedigrees. Genetic linkage mapping has implications even broader than its application to prenatal diagnosis or therapeutic strategy; analyzing mutations in important genes will illuminate basic mechanisms in molecular biology and the early events that lead to cancer and other disorders.     

SLEB3 in systemic lupus erythematosus (SLE) is strongly related to SLE families ascertained through neuropsychiatric manifestations

Seizures and psychosis are neuropsychiatric (NP) manifestations of a large number of systemic lupus erythematosus (SLE) patients. Since NP manifestations were part of the SLE phenotype for some, but not all SLE affecteds, we hypothesized that those SLE patient families with NP manifestations might be more genetically homogeneous at loci important to NP-related SLE, and hence have increased power to detect linkage. We identified 23 families of European-American (EA) origin and 20 families of African-American (AA) origin, in which at least one SLE patient in each family was diagnosed with the presence of NP manifestations. A total of 318 microsatellite markers at an average marker density of 11 cM were genotyped. Uncertainty of the genetic model led us to perform the initial genome scan by a multipoint non-parametric allele sharing linkage method. Once the evidence of linkage was suggestive, we then performed parametric model-based linkage by maximizing the relevant parameters to define a parsimonious genetic model. We found the maximum multipoint parametric LOD score was 5.19 and the non-parametric linkage score (Zlr) was 3.12 ( P=9x10(-4)) for EA NP pedigrees at 4p16, previously identified as SLEB3. The segregation behavior of this linked locus suggests a dominant mode of inheritance with an almost 100% homogeneous genetic effect in these pedigrees. The results demonstrated a significant increase of LOD score to detect SLEB3 when the families were further ascertained through NP, compared with the analysis of all EA SLE families together.     

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.     

Lysinuric Protein Intolerance (LPI) Gene Maps to the Long Arm of Chromosome 14

Lysinuric protein intolerance (LPI) is an autosomal recessive disease characterized by defective transport of cationic amino acids and by hyperammonemia. Linkage analysis in 20 Finnish LPI families assigned the LPI gene locus to the proximal long arm of chromosome 14. Recombinations placed the locus between framework markers D14S72 and MYH7, a 10-cM interval in which the markers D14S742, D14S50, D14S283, and TCRA showed no recombinations with the phenotype. The phenotype was in highly significant linkage disequilibrium with markers D14S50, D14S283, and TCRA. The strongest allelic association obtained with marker TCRA, resulting in a P(excess) value of .98, suggests that the LPI gene locus lies in close proximity to this marker, probably within a distance of < 100 kb.     

Replication study supports evidence for linkage to 9p24 in obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a severe psychiatric illness that is characterized by intrusive and senseless thoughts and impulses (obsessions) and by repetitive behaviors (compulsions). Family, twin, and segregation studies support the presence of both genetic and environmental susceptibility factors, and the only published genome scan for OCD identified a candidate region on 9p24 at marker D9S288 that met criteria for suggestive significance (Hanna et al. 2002). In an attempt to replicate this finding, we genotyped 50 pedigrees with OCD, using microsatellite markers spanning the 9p24 candidate region, and analyzed the data, using parametric and nonparametric linkage analyses under both a narrow phenotype model (DSM-IV OCD definite; 41 affected sib pairs) and a broad phenotype model (DSM-IV OCD definite and probable; 50 affected sib pairs). Similar to what was described by Hanna et al. (2002), our strongest findings came with the dominant parameters and the narrow phenotype model: the parametric signal peaked at marker D9S1792 with an HLOD of 2.26 ( alpha =0.59), and the nonparametric linkage signal (NPL) peaked at marker D9S1813 with an NPL of 2.52 (P=.006). These findings are striking in that D9S1813 and D9S1792 lie within 0.5 cM (<350 kb) of the original 9p24 linkage signal at D9S288; furthermore, pedigree-based association analyses also implicated the 9p24 candidate region by identifying two markers (D9S288 and GATA62F03) with modest evidence (P=.046 and .02, respectively) for association.     

Hemifacial microsomia: progress in understanding the genetic basis of a complex malformation syndrome

Hemifacial microsomia (HFM) is a common birth defect involving first and second branchial arch derivatives. The phenotype is extremely variable. In addition to craniofacial anomalies there may be cardiac, vertebral and central nervous system defects. The majority of cases are sporadic, but there is substantial evidence for genetic involvement in this condition, including rare familial cases that exhibit autosomal dominant inheritance. As an approach towards identifying molecular pathways involved in ear and facial development, we have ascertained both familial and sporadic cases of HFM. A genome wide search for linkage in two families with features of HFM was performed to identify the disease loci. In one family data were highly suggestive of linkage to a region of approximately 10.7 cM on chromosome 14q32, with a maximum multipoint lod score of 3.00 between microsatellite markers D14S987 and D14S65. This locus harbours the Goosecoid gene, an excellent candidate for HFM based on mouse expression and phenotype data. Coding region mutations were sought in the familial cases and in 120 sporadic cases, and gross rearrangements of the gene were excluded by Southern blotting. Evidence for genetic heterogeneity is provided by the second family, in which linkage was excluded from this region.     

Genetic linkage analysis in familial breast and ovarian cancer: results from 214 families. The Breast Cancer Linkage Consortium.

Breast cancer is known to have an inherited component, consistent in some families with autosomal dominant inheritance; in such families the disease often occurs in association with ovarian cancer. Previous genetic linkage studies have established that in some such families disease occurrence is linked to markers on chromosome 17q. This paper reports the results of a collaborative linkage study involving 214 breast cancer families, including 57 breast-ovarian cancer families; this represents almost all the known families with 17q linkage data. Six markers on 17q, spanning approximately 30 cM, were typed in the families. The aims of the study were to define more precisely the localization of the disease gene, the extent of genetic heterogeneity and the characteristics of linked families and to estimate the penetrance of the 17q gene. Under the assumption of no genetic heterogeneity, the strongest linkage evidence was obtained with D17S588 (maximum LOD score [Zmax] = 21.68 at female recombination fraction [theta f] = .13) and D17S579 (Zmax = 13.02 at theta f = .16). Multipoint linkage analysis allowing for genetic heterogeneity provided evidence that the predisposing gene lies between the markers D17S588 and D17S250, an interval whose genetic length is estimated to be 8.3 cM in males and 18.0 cM in females. This position was supported over other intervals by odds of 66:1. The location of the gene with respect to D17S579 could not be determined unequivocally. Under the genetic model used in the analysis, the best estimate of the proportion of linked breast-ovarian cancer families was 1.0 (lower LOD-1 limit 0.79). In contrast, there was significant evidence of genetic heterogeneity among the families without ovarian cancer, with an estimated 45% being linked. These results suggest that a gene(s) on chromosome 17q accounts for the majority of families in which both early-onset breast cancer and ovarian cancer occur but that other genes predisposing to breast cancer exist. By examining the fit of the linkage data to different penetrance functions, the cumulative risk associated with the 17q gene was estimated to be 59% by age 50 years and 82% by age 70 years. The corresponding estimates for the breast-ovary families were 67% and 76%, and those for the families without ovarian cancer were 49% and 90%; these penetrance functions did not differ significantly from one another.     

The promises and problems of linkage analysis by using the current canine genome map

There have been major advances in the canine genome map over recent years, with an agreed karyotype for all 38 pairs of autosomes and an integrated linkage-radiation-hybrid map now available. An individual dog breed represents a closed, isolated population, and many suffer genetic diseases. Some are homologs of human conditions, and dogs represent naturally occurring, potentially useful large mammal models. Other conditions may be unique to this species, but study of these offers new insights into metabolic pathways. Linkage analysis may identify a locus linked to a particular condition. Dog pedigrees facilitate these studies by being relatively more informative than human families, although planned breeding strategies to maximize the informativeness may be required. A study of a late-onset, acquired heart disease, dilated cardiomyopathy (DCM), in Newfoundland dogs highlights some of the problems. The disease occurs naturally in pedigrees within the UK in dogs kept as pets and breeding stock. Problems encountered include confirmation of the mode of inheritance and determination of phenotype. A genome-wide linkage analysis study with over 200 markers failed to detect linkage. The study indicates that more detailed linkage maps are required, and further synteny information between dog and human gained before study of diseases affecting naturally occurring families of dogs from inbred pedigrees may have maximal chance of success in order to utilize this model.     

Contribution of the hepatic lipase gene to the atherogenic lipoprotein phenotype in familial combined hyperlipidemia

Familial combined hyperlipidemia (FCH) is a common genetic lipid disorder with a frequency of 1-2% in the population. In addition to the hypercholesterolemia and/or hypertriglyceridemia that affected individuals exhibit, small, dense LDL particles and decreased HDL-cholesterol levels are traits frequently associated with FCH. Recently, we reported that families with FCH and families enriched for coronary artery disease (CAD) share genetic determinants for the atherogenic lipoprotein phenotype (ALP), a profile presenting with small, dense LDL particles, decreased HDL-cholesterol levels, and increased triglyceride levels. Other studies in normolipidemic populations have shown that the hepatic lipase (HL) gene is linked to HDL-cholesterol levels and that a polymorphism within the HL promoter (-514C-->T) is associated with increased HDL-cholesterol levels as well as larger, more buoyant LDL particles. In the present study, we tested whether the HL gene locus also contributes to ALP in a series of Dutch FCH families using nonparametric sibpair linkage analysis and association analysis. Evidence for linkage of LDL particle size (P < 0.019), HDL-cholesterol (P < 0.003), and triglyceride levels (P < 0.026) to the HL gene locus was observed. A genome scan in a subset of these families exhibited evidence for linkage of PPD (LOD = 2.2) and HDL-cholesterol levels (LOD = 1.2) to the HL gene locus as well. The -514C-->T promoter polymorphism was significantly associated (P < 0.0001) with higher HDL-cholesterol levels in the unrelated males of this population, but not in unrelated females. No association was observed between the polymorphism and LDL particle size or triglyceride levels. Our results provide support that ALP is a multigenic trait and suggest that the relationship between small, dense LDL particles, HDL-cholesterol, and triglyceride levels in FCH families is due, in part, to common genetic factors.     

Constructing genetic linkage maps under a tetrasomic model

An international consortium has launched the whole-genome sequencing of potato, the fourth most important food crop in the world. Construction of genetic linkage maps is an inevitable step for taking advantage of the genome projects for the development of novel cultivars in the autotetraploid crop species. However, linkage analysis in autopolyploids, the kernel of linkage map construction, is theoretically challenging and methodologically unavailable in the current literature. We present here a theoretical analysis and a statistical method for tetrasomic linkage analysis with dominant and/or codominant molecular markers. The analysis reveals some essential properties of the tetrasomic model. The method accounts properly for double reduction and incomplete information of marker phenotype in regard to the corresponding phenotype in estimating the coefficients of double reduction and recombination frequency and in testing their significance by using the marker phenotype data. Computer simulation was developed to validate the analysis and the method and a case study with 201 AFLP and SSR markers scored on 228 full-sib individuals of autotetraploid potato is used to illustrate the utility of the method in map construction in autotetraploid species.     

A susceptibility locus for migraine with aura, on chromosome 4q24

Migraine is a complex neurovascular disorder with substantial evidence supporting a genetic contribution. Prior attempts to localize susceptibility loci for common forms of migraine have not produced conclusive evidence of linkage or association. To date, no genomewide screen for migraine has been published. We report results from a genomewide screen of 50 multigenerational, clinically well-defined Finnish families showing intergenerational transmission of migraine with aura (MA). The families were screened using 350 polymorphic microsatellite markers, with an average intermarker distance of 11 cM. Significant evidence of linkage was found between the MA phenotype and marker D4S1647 on 4q24. Using parametric two-point linkage analysis and assuming a dominant mode of inheritance, we found for this marker a maximum LOD score of 4.20 under locus homogeneity (P=.000006) or locus heterogeneity (P=.000011). Multipoint parametric (HLOD = 4.45; P=.0000058) and nonparametric (NPL(all) = 3.43; P=.0007) analyses support linkage in this region. Statistically significant linkage was not observed in any other chromosomal region.     

Fine mapping and SNP analysis of positional candidates at the preeclampsia susceptibility locus (PREG1) on chromosome 2

Genome scans in Icelandic, Australian and New Zealand, and Finnish families have localized putative susceptibility loci for preeclampsia/ eclampsia to chromosome 2. The locus mapped in the Australian and New Zealand study (designated PREG1) was thought to be the same locus as that identified in the Icelandic study. In both these studies, two distinct quantitative trait locus (QTL) regions were evident on chromosome 2. Here, we describe our fine mapping of the PREG1 locus and a genetic analysis of two positional candidate genes. Twenty-five additional microsatellite markers were genotyped within the 74-cM linkage region defined by the combined Icelandic and Australian and New Zealand genome scans. The overall position and shape of the localization evidence obtained using nonparametric multipoint analysis did not change from that seen previously in our 10-cM resolution genome scan; two peaks were displayed, one on chromosome 2p at marker D2S388 (107.46 cM) and the other on chromosome 2q at 151.5 cM at marker D2S2313. Using the robust two-point linkage analysis implemented in the Analyze program, all 25 markers gave positive LOD scores with significant evidence of linkage being seen at marker D2S2313 (151.5 cM), achieving a LOD score of 3.37 under a strict diagnostic model. Suggestive evidence of linkage was seen at marker D2S388 (107.46 cM) with a LOD score of 2.22 under the general diagnostic model. Two candidate genes beneath the peak on chromosome 2p were selected for further analysis using public single nucleotide polymorphisms (SNPs) within these genes. Maximum LOD scores were obtained for an SNP in TACR1 (LOD = 3.5) and for an SNP in TCF7L1 (LOD = 3.33), both achieving genome-wide significance. However, no evidence of association was seen with any of the markers tested. These data strongly support the presence of a susceptibility gene on chromosome 2p11-12 and substantiate the possibility of a second locus on chromosome 2q23.