Power and Robustness of Linkage Tests for Quantitative Traits in General Pedigrees

There are numerous statistical methods for quantitative trait linkage analysis in human studies. An ideal such method would have high power to detect genetic loci contributing to the trait, would be robust to non-normality in the phenotype distribution, would be appropriate for general pedigrees, would allow the incorporation of environmental covariates, and would be appropriate in the presence of selective sampling. We recently described a general framework for quantitative trait linkage analysis, based on generalized estimating equations, for which many current methods are special cases. This procedure is appropriate for general pedigrees and easily accommodates environmental covariates. In this paper, we use computer simulations to investigate the power robustness of a variety of linkage test statistics built upon our general framework. We also propose two novel test statistics that take account of higher moments of the phenotype distribution, in order to accommodate non-normality. These new linkage tests are shown to have high power and to be robust to non-normality. While we have not yet examined the performance of our procedures in the context of selective sampling via computer simulations, the proposed tests satisfy all of the other qualities of an ideal quantitative trait linkage analysis method.

Genetic analysis of tumor progression susceptibility in the mouse skin model

In the field of chemical carcinogenesis the use of animal models has proved to be a useful tool in dissecting the multistage process of tumor formation. In this regard the outbred SENCAR mouse has been the strain of choice in the analysis of skin carcinogenesis given its high sensitivity to the chemically induced acquisition of premalignant lesions, papillomas, and the later progression of these lesions into squamous cell carcinomas (SCC).^ The derivation of an inbred strain from the SENCAR stock called SSIN, that in spite of a high sensitivity to the development of papillomas lack the ability to transform these premalignant lesions into SCC, suggested that tumor promotion and progression were under the genetic control of different sets of genes.^ In the present study the nature of susceptibility to tumor progression was investigated. Analysis of F1 hybrids between the outbred SENCAR and SSIN mice suggested that there is at least one dominant gene responsible for susceptibility to tumor progression.^ Later development of another inbred strain from the outbred SENCAR stock, that had sensitivity to both tumor promotion and progression, allowed the formulation of a more accurate genetic model. Using this newly derived line, SENCAR B/Pt. and SSIN it was determined that there is one dominant tumor progression susceptibility gene. Linkage analysis showed that this gene maps to mouse chromosome 14 and it was possible to narrow the region to a 16 cM interval.^ In order to better characterize the nature of the progression susceptibility differences between these two strains, their proliferative pattern was investigated. It was found that SENCAR B/Pt, have an enlarged proliferative compartment with overexpression of cyclin D1, p16 and p21. Further studies showed an aberrant overexpression of TGF-$\beta$ in the susceptible strain, an increase in apoptosis, p53 protein accumulation and early loss of connexin 26. These results taken together suggest that papillomas in the SENCAR B/Pt. mice have higher proliferation and may have an increase in genomic instability, these two factors would contribute to a higher sensitivity to tumor progression. ^

Linkage and association of candidate obesity genes in Mexican-Americans from Starr County, Texas

Obesity and related chronic diseases represent a tremendous public health burden among Mexican Americans, a young and rapidly-expanding population. This study investigated the impact of variation within eight candidate obesity genes, which include leptin (LEP), leptin receptor (LEPR), neuropeptide Y (NPY), NPYY1 receptor (NPYY1), glucagon-like peptide-1 (GLP-1), GLP-1 receptor (GLP1R), beta-3 adrenergic receptor (β3AR), and uncoupling protein (UCP1), on variation in human obesity status and/or quantitative traits related to obesity in Mexican Americans from Starr County, Texas. The Trp64Arg polymorphism within β3AR was typed in 820 random individuals and 240 pedigrees (N = 2,044). The Arg allele frequency was significantly greater in obese versus non-obese individuals (0.20 versus 0. 15, respectively). In addition, within the random sample, the Arg allele was associated with significantly greater body weight (p = 0.031) and body mass index (BMI, p = 0.008) than the Trp allele. In the family sample, the Trp64Arg locus was also linked to percent fat (p = 0.045) but not to body weight or BMI. No linkage between obesity, diabetes, hypertension, or gallbladder disease and the Trp64Arg mutation was observed in families using affected sib pair linkage analysis or the transmission disequilibrium test. Microsatellite markers proximate to the remaining seven genes were typed in 302 individuals from 59 families. Sib pair linkage analysis provided evidence for linkage between obesity and NPY within affected sibling pairs (p = 0.042; n = 170 pairs). NPY was also linked to weight (p = 0.020), abdominal circumference (p = 0.031), hip circumference (p = 0.012), DBP (p ≤ 0.005), and a composite measure of body mass/fat (p ≤ 0.048) in all sibling pairs (n = 545 pairs). Additionally, LEP was linked to waist/hip ratio (p ≤ 0.009), total cholesterol (p ≤ 0.030), and HDL cholesterol (p ≤ 0.026), and LEPR was linked to fasting blood glucose (p ≤ 0.018) and DBP (p ≤ 0.003). Subsequent to the linkage analyses, the NPY gene was sequenced and eight variant sites identified. Two variant sites (-880I/D and 69I/D) were typed in a random sample of 914 individuals. The 880I/D variant was significantly associated with waist/hip ratio (p = 0.035) in the entire sample (N = 914) and with BMI (p = 0. 031), abdominal circumference (p = 0.044), and waist/hip ratio (p = 0.041) in a non-obese subsample (BW < 30 kg/m2, n = 594). The 69I/D variant was a rare mutation observed in only one pedigree and was not associated with obesity or body size/mass within this pedigree. Results of this study indicate that variation at or near β3AR, LEP, LEPR, and NPY may exert effects which increase obesity susceptibility and influence obesity-related measures in this population. ^

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

Genome scan of human systemic lupus erythematosus: Evidence for linkage on chromosome 1q in African-American pedigrees

Systemic lupus erythematosus (SLE) is an autoimmune disorder characterized by production of autoantibodies against intracellular antigens including DNA, ribosomal P, Ro (SS-A), La (SS-B), and the spliceosome. Etiology is suspected to involve genetic and environmental factors. Evidence of genetic involvement includes: associations with HLA-DR3, HLA-DR2, Fcγ receptors (FcγR) IIA and IIIA, and hereditary complement component deficiencies, as well as familial aggregation, monozygotic twin concordance >20%, λs > 10, purported linkage at 1q41–42, and inbred mouse strains that consistently develop lupus. We have completed a genome scan in 94 extended multiplex pedigrees by using model-based linkage analysis. Potential [log10 of the odds for linkage (lod) > 2.0] SLE loci have been identified at chromosomes 1q41, 1q23, and 11q14–23 in African-Americans; 14q11, 4p15, 11q25, 2q32, 19q13, 6q26–27, and 12p12–11 in European-Americans; and 1q23, 13q32, 20q13, and 1q31 in all pedigrees combined. An effect for the FcγRIIA candidate polymorphism) at 1q23 (lod = 3.37 in African-Americans) is syntenic with linkage in a murine model of lupus. Sib-pair and multipoint nonparametric analyses also support linkage (P < 0.05) at nine loci detected by using two-point lod score analysis (lod > 2.0). Our results are consistent with the presumed complexity of genetic susceptibility to SLE and illustrate racial origin is likely to influence the specific nature of these genetic effects.

Cloning and functional analysis of two gibberellin 3β-hydroxylase genes that are differently expressed during the growth of rice

We have cloned two gibberellin (GA) 3β-hydroxylase genes, OsGA3ox1 and OsGA3ox2, from rice by screening a genomic library with a DNA fragment obtained by PCR using degenerate primers. We have used full-scan GC-MS and Kovats retention indices to show function for the two encoded recombinant fusion proteins. Both proteins show 3β-hydroxylase activity for the steps GA20 to GA1, GA5 to GA3, GA44 to GA38, and GA9 to GA4. In addition, indirect evidence suggests that the OsGA3ox1 protein also has 2,3-desaturase activity, which catalyzes the steps GA9 to 2,3-dehydro-GA9 and GA20 to GA5 (2,3-dehydro GA20), and 2β-hydroxylase activity, which catalyzes the steps GA1 to GA8 and GA4 to GA34. Molecular and linkage analysis maps the OsGA3ox1 gene to the distal end of the short arm of chromosome 5; the OsGA3ox2 gene maps to the distal end of the short arm of chromosome 1 that corresponds to the D18 locus. The association of the OsGA3ox2 gene with the d18 locus is confirmed by sequence and complementation analysis of three d18 alleles. Complementation of the d18-AD allele with the OxGA3ox2 gene results in transgenic plants with a normal phenotype. Although both genes show transient expression, the highest level for OsGA3ox1 is from unopened flower. The highest level for OsGA3ox2 is from elongating leaves.

Genetic analysis of type 1 diabetes using whole genome approaches.

Whole genome linkage analysis of type 1 diabetes using affected sib pair families and semi-automated genotyping and data capture procedures has shown how type 1 diabetes is inherited. A major proportion of clustering of the disease in families can be accounted for by sharing of alleles at susceptibility loci in the major histocompatibility complex on chromosome 6 (IDDM1) and at a minimum of 11 other loci on nine chromosomes. Primary etiological components of IDDM1, the HLA-DQB1 and -DRB1 class II immune response genes, and of IDDM2, the minisatellite repeat sequence in the 5' regulatory region of the insulin gene on chromosome 11p15, have been identified. Identification of the other loci will involve linkage disequilibrium mapping and sequencing of candidate genes in regions of linkage.

Analysis of Plasminogen Genetic Variants in Multiple Sclerosis Patients.

Multiple sclerosis (MS) is a prevalent neurological disease of complex etiology. Here, we describe the characterization of a multi-incident MS family that nominated a rare missense variant (p.G420D) in plasminogen (PLG) as a putative genetic risk factor for MS. Genotyping of PLG p.G420D (rs139071351) in 2160 MS patients, and 886 controls from Canada, identified 10 additional probands, two sporadic patients and one control with the variant. Segregation in families harboring the rs139071351 variant, identified p.G420D in 26 out of 30 family members diagnosed with MS, 14 unaffected parents, and 12 out of 30 family members not diagnosed with disease. Despite considerably reduced penetrance, linkage analysis supports cosegregation of PLG p.G420D and disease. Genotyping of PLG p.G420D in 14446 patients, and 8797 controls from Canada, France, Spain, Germany, Belgium, and Austria failed to identify significant association with disease (P = 0.117), despite an overall higher prevalence in patients (OR = 1.32; 95% CI = 0.93-1.87). To assess whether additional rare variants have an effect on MS risk, we sequenced PLG in 293 probands, and genotyped all rare variants in cases and controls. This analysis identified nine rare missense variants, and although three of them were exclusively observed in MS patients, segregation does not support pathogenicity. PLG is a plausible biological candidate for MS owing to its involvement in immune system response, blood-brain barrier permeability, and myelin degradation. Moreover, components of its activation cascade have been shown to present increased activity or expression in MS patients compared to controls; further studies are needed to clarify whether PLG is involved in MS susceptibility.

The second generation of the International Equine Gene Mapping Workshop half-sibling linkage map

A low-density, male-based linkage map was constructed as one of the objectives of the International Equine Gene Mapping Workshop. Here we report the second generation map based on testing 503 half-sibling offspring from 13 sire families for 344 informative markers using the crimap program. The multipoint linkage analysis localized 310 markers (90%) with 257 markers being linearly ordered. The map included 34 linkage groups representing all 31 autosomes and spanning 2262 cM with an average interval between loci of 10.1 cM. This map is a milestone in that it is the first map with linkage groups assigned to each of the 31 automosomes and a single linkage group to all but three chromosomes.

Identification of a major locus conferring resistance to powdery mildew (Erysiphe polygoni DC) in mungbean (Vigna radiata L. Wilczek) by QTL analysis

A major locus conferring resistance to the causal organism of powdery mildew, Erysiphe polygoni DC,, in mungbean (Vigna radiata L. Wilczek) was identified using QTL analysis with a population of 147 recombinant inbred individuals. The population was derived from a cross between 'Berken', a highly susceptible variety, and ATF 3640, a highly resistant line. To test for response to powdery mildew, F-7 and F-8 lines were inoculated by dispersing decaying mungbean leaves with residual conidia of E. polygoni amongst the young plants to create an artificial epidemic and assayed in a glasshouse facility. To generate a linkage map, 322 RFLP clones were tested against the two parents and 51 of these were selected to screen the mapping population. The 51 probes generated 52 mapped loci, which were used to construct a linkage map spanning 350 cM of the mungbean genome over 10 linkage groups. Using these markers, a single locus was identified that explained up to a maximum of 86% of the total variation in the resistance response to the pathogen.

Genome scan meta-analysis of schizophrenia and bipolar disorder, part II: Schizophrenia

Schizophrenia is a common disorder with high heritability and a 10-fold increase in risk to siblings of probands. Replication has been inconsistent for reports of significant genetic linkage. To assess evidence for linkage across studies, rank-based genome scan meta-analysis (GSMA) was applied to data from 20 schizophrenia genome scans. Each marker for each scan was assigned to 1 of 120 30-cM bins, with the bins ranked by linkage scores (1 = most significant) and the ranks averaged across studies (R-avg) and then weighted for sample size (rootN[affected cases]). A permutation test was used to compute the probability of observing, by chance, each bin's average rank (P-AvgRnk) or of observing it for a bin with the same place (first, second, etc.) in the order of average ranks in each permutation (P-ord). The GSMA produced significant genomewide evidence for linkage on chromosome 2q (P-AvgRnk

Cosegregation analysis of natriuretic peptide genes and blood pressure in the spontaneously hypertensive rat

1. The natriuretic peptide precursor A (Nppa) and B (Nppb) genes are candidate genes for hypertension and cardiac hypertrophy in the spontaneously hypertensive rat (SHR). The purpose of the present study was to determine the role of the Nppa and Nppb genes in the development of hypertension in the SHR. 2. A cohort (n = 162) of F2 segregating intercross animals was established between strains of hypertensive SHR and normotensive Wistar-Kyoto rats. Blood pressure and heart weight were measured in each rat at 12-16 weeks of age. Rats were genotyped using 11 informative microsatellite markers, distributed in the vicinity of the Nppa marker on rat chromosome 5 including an Nppb marker. The phenotype values were compared with genotype using the computer package MAP-MAKER 3.0 (Whitehead Institute, Boston, MA, USA) to determine whether there was a link between the genetic variants of the natriuretic peptide family and blood pressure or cardiac hypertrophy. 3. A strong correlation was observed between the Nppa marker and blood pressure. A quantitative trait locus (QTL) for blood pressure on chromosome 5 was identified between the Nppa locus and the D5Mgh15 marker, less than 2 cM from the Nppa locus. The linkage score for the blood pressure QTL on chromosome 5 was 3.8 and the QTL accounted for 43% of the total variance of systolic blood pressure, 54% of diastolic blood pressure and 59% of mean blood pressure. No association was found between the Nppb gene and blood pressure. This is the first report of linkage between the Nppa marker and blood pressure in the rat. There was no correlation between the Nppa or Nppb genes or other markers in this region and either heart weight or left ventricular weight in F2 rats. 4. These findings suggest the existence of a blood pressure-dependent Nppa marker variant or a gene close to Nppa predisposing to spontaneous hypertension in the rat. It provides a strong foundation for further detailed genetic studies in congenic strains, which may help to narrow down the location of this gene and lead to positional cloning.

Analysis of Melanoma Onset: Assessing Familial Aggregation by Using Estimating Equations and Fitting Variance Components via Bayesian Random Effects Models

We investigate whether relative contributions of genetic and shared environmental factors are associated with an increased risk in melanoma. Data from the Queensland Familial Melanoma Project comprising 15,907 subjects arising from 1912 families were analyzed to estimate the additive genetic, common and unique environmental contributions to variation in the age at onset of melanoma. Two complementary approaches for analyzing correlated time-to-onset family data were considered: the generalized estimating equations (GEE) method in which one can estimate relationship-specific dependence simultaneously with regression coefficients that describe the average population response to changing covariates; and a subject-specific Bayesian mixed model in which heterogeneity in regression parameters is explicitly modeled and the different components of variation may be estimated directly. The proportional hazards and Weibull models were utilized, as both produce natural frameworks for estimating relative risks while adjusting for simultaneous effects of other covariates. A simple Markov Chain Monte Carlo method for covariate imputation of missing data was used and the actual implementation of the Bayesian model was based on Gibbs sampling using the free ware package BUGS. In addition, we also used a Bayesian model to investigate the relative contribution of genetic and environmental effects on the expression of naevi and freckles, which are known risk factors for melanoma.

Multicenter linkage study of schizophrenia loci on chromosome 22q

The hypothesis of the existence of one or more schizophrenia susceptibility loci on chromosome 22q is supported by reports of genetic linkage and association, meta-analyses of linkage, and the observation of elevated risk for psychosis in people with velocardiofacial syndrome, caused by 22q11 microdeletions. We tested this hypothesis by evaluating 10 microsatellite markers spanning 22q in a multicenter sample of 779 pedigrees. We also incorporated age at onset and sex into the analysis as covariates. No significant evidence for linkage to schizophrenia or for linkage associated with earlier age at onset, gender, or heterogeneity across sites was observed. We interpret these findings to mean that the population-wide effects of putative 22q schizophrenia susceptibility loci are too weak to detect with linkage analysis even in large samples.

Genomewide significant linkage to migrainous Hheadache on chromosome 5q21

Familial typical migraine is a common, complex disorder that shows strong familial aggregation. Using latent-class analysis (LCA), we identified subgroups of people with migraine/severe headache in a community sample of 12,245 Australian twins (60% female), drawn from two cohorts of individuals aged 23-90 years who completed an interview based on International Headache Society criteria. We report results from genomewide linkage analyses involving 756 twin families containing a total of 790 independent sib pairs ( 130 affected concordant, 324 discordant, and 336 unaffected concordant for LCA-derived migraine). Quantitative-trait linkage analysis produced evidence of significant linkage on chromosome 5q21 and suggestive linkage on chromosomes 8, 10, and 13. In addition, we replicated previously reported typical-migraine susceptibility loci on chromosomes 6p12.2-p21.1 and 1q21-q23, the latter being within 3 cM of the rare autosomal dominant familial hemiplegic migraine gene (ATP1A2), a finding which potentially implicates ATP1A2 in familial typical migraine for the first time. Linkage analyses of individual migraine symptoms for our six most interesting chromosomes provide tantalizing hints of the phenotypic and genetic complexity of migraine. Specifically, the chromosome 1 locus is most associated with phonophobia; the chromosome 5 peak is predominantly associated with pulsating headache; the chromosome 6 locus is associated with activity-prohibiting headache and photophobia; the chromosome 8 locus is associated with nausea/vomiting and moderate/severe headache; the chromosome 10 peak is most associated with phonophobia and photophobia; and the chromosome 13 peak is completely due to association with photophobia. These results will prove to be invaluable in the design and analysis of future linkage and linkage disequilibrium studies of migraine.