Multivariate QTL linkage analysis suggests a QTL for platelet count on chromosome 19q

Platelet count is a highly heritable trait with genetic factors responsible for around 80% of the phenotypic variance. We measured platelet count longitudinally in 327 monozygotic and 418 dizygotic twin pairs at 12, 14 and 16 years of age. We also performed a genome-wide linkage scan of these twins and their families in an attempt to localize QTLs that influenced variation in platelet concentrations. Suggestive linkage was observed on chromosome 19q13.13-19q13.31 at 12 (LOD=2.12, P=0.0009), 14 (LOD=2.23, P=0.0007) and 16 (LOD=1.01, P=0.016) years of age and multivariate analysis of counts at all three ages increased the LOD to 2.59 (P=0.0003). A possible candidate in this region is the gene for glycoprotein VI, a receptor involved in platelet aggregation. Smaller linkage peaks were also seen at 2p, 5p, 5q, 10p and 15q. There was little evidence for linkage to the chromosomal regions containing the genes for thrombopoietin (3q27) and the thrombopoietin receptor (1q34), suggesting that polymorphisms in these genes do not contribute substantially to variation in platelet count between healthy individuals.

A genome scan for eye color in 502 twin families: Most variation is due to a QTL on chromosome 15q

We have rated eye color on a 3-point scale (1=blue/grey, 2=hazel/green, 3=brown) in 502 twin families and carried out a 5-10 cM genome scan (400-757 markers). We analyzed eye color as a threshold trait and performed multipoint sib pair linkage analysis using variance components analysis in Mx. A lod of 19.2 was found at the marker D15S1002, less than 1 cM from OCA2, which has been previously implicated in eye color variation. We estimate that 74% of variance in eye color liability is due to this QTL and a further 18% due to polygenic effects. However, a large shoulder on this peak suggests that other loci affecting eye color may be telomeric of OCA2 and inflating the QTL estimate. No other peaks reached genome-wide significance, although lods >2 were seen on 5p and 14q and lods >1 were additionally seen on chromosomes 2, 3, 6, 7, 8, 9, 17 and 18. Most of these secondary peaks were reduced or eliminated when we repeated the scan as a two locus analysis with the 15q linkage included, although this does not necessarily exclude them as false positives. We also estimated the interaction between the 15q QTL and the other marker locus but there was only minor evidence for additive x additive epistasis. Elaborating the analysis to the full two-locus model including non-additive main effects and interactions did not strengthen the evidence for epistasis. We conclude that most variation in eye color in Europeans is due to polymorphism in OCA2 but that there may be modifiers at several other loci.

Lack of association between genetic markers on chromosome 16q22-Q24 and type 1 diabetes in Russian affected families

Aim To evaluate whether the T1D susceptibility locus on chromosome 16q contributes to the genetic susceptibility to T1D in Russian patients. Method Thirteen microsatellite markers, spanning a 47-centimorgan genomic region on 16q22-q24 were evaluated for linkage to T1D in 98 Russian multiplex families. Multipoint logarithm of odds (LOD) ratio (MLS) and nonparametric LOD (NPL) values were computed for each marker, using GENEHUNTER 2.1 software. Four microsatellites (D16S422, D16S504, D16S3037, and D16S3098) and 6 biallelic markers in 2 positional candidate genes, ICSBP1 and NQO1, were additionally tested for association with T1D in 114 simplex families, using transmission disequilibrium test (TDT). Results A peak of linkage (MLS = 1.35, NPL = 0.91) was shown for marker D16S750, but this was not significant (P = 0.18). The subsequent linkage analysis in the subset of 46 multiplex families carrying a common risk HLA-DR4 haplotype increased peak MLS and NPL values to 1.77 and 1.22, respectively, but showed no significant linkage (P = 0.11) to T1D in the 16q22-q24 genomic region. TDT analysis failed to find significant association between these markers and disease, even after the conditioning for the predisposing HLA-DR4 haplotype. Conclusion Our results did not support the evidence for the susceptibility locus to T1D on chromosome 16q22-24 in the Russian family data set. The lack of association could reflect genetic heterogeneity of type 1 diabetes in diverse ethnic groups.

The combined SPE : ToxY-PAM phytotoxicity assay; application and appraisal of a novel biomonitoring tool for the aquatic environment

Mounting concerns regarding the environmental impact of herbicides has meant a growing requirement for accurate, timely information regarding herbicide residue contamination of, in particular, aquatic systems. Conventional methods of detection remain limited in terms of practicality due to high costs of operation and the specialised information that analysis provides. A new phytotoxicity bioassay was trialled for the detection of herbicide residues in filter-purified (Milli-Q) as well as natural waters. The performance of the system, which combines solid-phase extraction (SPE) with the ToxY-PAM dual-channel yield analyser (Heinz Walz GmbH), was tested alongside the traditional method of liquid chromatography-mass spectrometry (LC-MS). The assay methodology was found to be highly sensitive (LOD 0.1 ng L-1 diuron) with good reproducibility. The study showed that the assay protocol is time effective and can be employed for the aquatic screening of herbicide residues in purified as well as natural waters.

The importance of modelling heterogeneity in complex disease: application to NIMH Schizophrenia Genetics Initiative data

As for other complex diseases, linkage analyses of schizophrenia (SZ) have produced evidence for numerous chromosomal regions, with inconsistent results reported across studies. The presence of locus heterogeneity appears likely and may reduce the power of linkage analyses if homogeneity is assumed. In addition, when multiple heterogeneous datasets are pooled, intersample variation in the proportion of linked families ( a) may diminish the power of the pooled sample to detect susceptibility loci, in spite of the larger sample size obtained. We compare the significance of linkage. findings obtained using allele- sharing LOD scores ( LODexp) - which assume homogeneity - and heterogeneity LOD scores ( HLOD) in European American and African American NIMH SZ families. We also pool these two samples and evaluate the relative power of the LODexp and two different heterogeneity statistics. One of these ( HLOD- P) estimates the heterogeneity parameter a only in aggregate data, while the second ( HLOD- S) determines a separately for each sample. In separate and combined data, we show consistently improved performance of HLOD scores over LODexp. Notably, genome-wide significant evidence for linkage is obtained at chromosome 10p in the European American sample using a recessive HLOD score. When the two samples are combined, linkage at the 10p locus also achieves genome-wide significance under HLOD- S, but not HLOD- P. Using HLOD- S, improved evidence for linkage was also obtained for a previously reported region on chromosome 15q. In linkage analyses of complex disease, power may be maximised by routinely modelling locus heterogeneity within individual datasets, even when multiple datasets are combined to form larger samples.

Genomewide linkage study in 1,176 affected sister pair families identifies a significant susceptibility locus for endometriosis on chromosome 10q26

Endometriosis is a common gynecological disease that affects up to 10% of women in their reproductive years. It causes pelvic pain, severe dysmenorrhea, and subfertility. The disease is defined as the presence of tissue resembling endometrium in sites outside the uterus. Its cause remains uncertain despite 150 years of hypothesis-driven research, and thus the therapeutic options are limited. Disease predisposition is inherited as a complex genetic trait, which provides an alternative route to understanding the disease. We seek to identify susceptibility loci, using a positional-cloning approach that starts with linkage analysis to identify genomic regions likely to harbor these genes. We conducted a linkage study of 1,176 families ( 931 from an Australian group and 245 from a U. K. group), each with at least two members-mainly affected sister pairs-with surgically diagnosed disease. We have identified a region of significant linkage on chromosome 10q26 ( maximum LOD score [MLS] of 3.09; genomewide P = .047) and another region of suggestive linkage on chromosome 20p13 MLS p 2.09). Minor peaks with MLS > 1.0) were found on chromosomes 2, 6, 7, 8, 12, 14, 15, and 17. This is the first report of linkage to a major locus for endometriosis. The findings will facilitate discovery of novel positional genetic variants that influence the risk of developing this debilitating disease. Greater understanding of the aberrant cellular and molecular mechanisms involved in the etiology and pathophysiology of endometriosis should lead to better diagnostic methods and targeted treatments.

A genome scan for epidermal skin pattern in adolescent twins reveals suggestive linkage on 12p13.31

Deterioration in stratum corneum reticular patterning (skin pattern or skin wrinkling) has been associated with increased rates of solar keratoses and skin cancer. A previous analysis of data from the twin sample used in this investigation has shown that 86% of the variation in skin pattern is genetic at age 12 and 62% in an adult sample (mean age 47.5). Variation due to genetic influences is likely to be influenced by more than one locus. Here, we present results of a genome-wide linkage scan of skin pattern in adolescent twins and siblings from 428 nuclear twin families. Sib-pair linkage analysis was performed on skin pattern data collected from twins at age 12 (378 informative families) and 14 (316 families). Suggestive linkage was found at marker D12S397 (12p13.31, logarithm of the odds (lod) 1.94), when the effect of the trait locus was modelled to influence the skin pattern equally at both ages 12 and 14. In the same analysis, a peak was seen at 4q23 with a lod score of 1.55. A possible candidate for the peak at 12p13.31 is the protease inhibitor, alpha-2-macroglobulin.

Replicated linkage for eye color on 15q using comparative ratings of sibling pairs

The aim of the study was to perform a genetic linkage analysis for eye color, for comparative data. Similarity in eye color of mono- and dizygotic twins was rated by the twins' mother, their father and/or the twins themselves. For 4748 twin pairs the similarity in eye color was available on a three point scale (not at all alike-somewhat alike-completely alike), absolute eye color on individuals was not assessed. The probability that twins were alike for eye color was calculated as a weighted average of the different responses of all respondents on several different time points. The mean probability of being alike for eye color was 0.98 for MZ twins (2167 pairs), whereas the mean probability for DZ twins was 0.46 (2537 pairs), suggesting very high heritability for eye color. For 294 DZ twin pairs genome-wide marker data were available. The probability of being alike for eye color was regressed on the average amount of IBD sharing. We found a peak LOD-score of 2.9 at chromosome 15q, overlapping with the region recently implicated for absolute ratings of eye color in Australian twins [Zhu, G., Evans, D. M., Duffy, D. L., Montgomery, G. W., Medland, S. E., Gillespie, N. A., Ewen, K. R., Jewell, M., Liew, Y. W., Hayward, N. K., Sturm, R. A., Trent, J. M., and Martin, N. G. (2004). Twin Res. 7:197-210] and containing the OCA2 gene, which is the major candidate gene for eye color [Sturm, R. A. Teasdale, R. D, and Box, N. F. (2001). Gene 277:49-62]. Our results demonstrate that comparative measures on relatives can be used in genetic linkage analysis.

Genome-wide scan of IQ finds significant linkage to a quantitative trait locus on 2q

A genome-wide linkage scan of 795 microsatellite markers (761 autosomal, 34 X chromosome) was performed on Multidimensional Aptitude Battery subtests and verbal, performance and full scale scores, the WAIS-R Digit Symbol subtest, and two word-recognition tests (Schonell Graded Word Reading Test, Cambridge Contextual Reading Test) highly predictive of IQ. The sample included 361 families comprising 2-5 siblings who ranged in age from 15.7 to 22.2 years; genotype, but not phenotype, data were available for 81% of parents. A variance components analysis which controlled for age and sex effects showed significant linkage for the Cambridge reading test and performance IQ to the same region on chromosome 2, with respective LOD scores of 4.15 and 3.68. Suggestive linkage (LOD score > 2.2) for various measures was further supported on chromosomes 6, 7, 11, 14, 21 and 22. Where location of linkage peaks converged for IQ subtests within the same scale, the overall scale score provided increased evidence for linkage to that region over any individual subtest. Association studies of candidate genes, particularly those involved in neural transmission and development, will be directed to genes located under the linkage peaks identified in this study.

Genomewide Linkage Scan of 409 European-ancestry and African American Families with Schizophrenia: Suggestive Evidence of Linkage at 8p23.3-p21.2 and 11p13.1-q14.1 in the Combined Sample

We report the clinical characteristics of a schizophrenia sample of 409 pedigrees-263 of European ancestry ( EA) and 146 of African American ancestry ( AA)-together with the results of a genome scan ( with a simple tandem repeat polymorphism interval of 9 cM) and follow-up fine mapping. A family was required to have a proband with schizophrenia ( SZ) and one or more siblings of the proband with SZ or schizoaffective disorder. Linkage analyses included 403 independent full-sibling affected sibling pairs ( ASPs) ( 279 EA and 124 AA) and 100 all-possible half-sibling ASPs ( 15 EA and 85 AA). Nonparametric multipoint linkage analysis of all families detected two regions with suggestive evidence of linkage at 8p23.3-q12 and 11p11.2-q22.3 ( empirical Z likelihood-ratio score [ Z(lr)] threshold >= 2.65) and, in exploratory analyses, two other regions at 4p16.1-p15.32 in AA families and at 5p14.3-q11.2 in EA families. The most significant linkage peak was in chromosome 8p; its signal was mainly driven by the EA families. Z(lr) scores >= 2.0 in 8p were observed from 30.7 cM to 61.7 cM ( Center for Inherited Disease Research map locations). The maximum evidence in the full sample was a multipoint Z(lr) of 3.25 ( equivalent Kong-Cox LOD of 2.30) near D8S1771 ( at 52 cM); there appeared to be two peaks, both telomeric to neuregulin 1 ( NRG1). There is a paracentric inversion common in EA individuals within this region, the effect of which on the linkage evidence remains unknown in this and in other previously analyzed samples. Fine mapping of 8p did not significantly alter the significance or length of the peak. We also performed fine mapping of 4p16.3-p15.2, 5p15.2-q13.3, 10p15.3-p14, 10q25.3-q26.3, and 11p13-q23.3. The highest increase in Z(lr) scores was observed for 5p14.1-q12.1, where the maximum Z(lr) increased from 2.77 initially to 3.80 after fine mapping in the EA families.

Genome-wide linkage scan for loci influencing plasma triglycerides

Background: Plasma triglyceride concentration is known to be a significant risk factor for cardiovascular disease (CVD). Previous studies have found that the level of triglycerides is strongly influenced by genetic factors. Methods: To identify quantitative trait loci influencing triglycerides, we conducted a genome-wide linkage scan on data from 485 Australian adult dizygotic twin pairs. Prior to linkage analysis, triglyceride values were adjusted for the effects of covariates including age, sex, time since last meal, time of blood collection (CT) and time to plasma separation. Results: The heritability estimate for ln(triglyceride) adjusted for all above fixed effects was 0.49. The highest multipoint LOD score observed was 2.94 (genome-wide p=0.049) on chromosome 7 (at 65cM). This 7p region contains several candidate genes. Two other regions with suggestive multipoint LOD scores were also identified on chromosome 4 (LOD score=2.26 at 62cM) and chromosome X (LOD score=2.01 at 81cM). Conclusions: The linkage peaks found represent newly identified regions for more detailed study, in particular the significant linkage observed on chromosome 7p13. \ (c) 2006 Elsevier B.V. All rights reserved.

Linkage analyses of event-related potential slow wave phenotypes recorded in a working memory task

Working memory is an essential component of wide-ranging cognitive functions. It is a complex genetic trait probably influenced by numerous genes that individually have only a small influence. These genes may have an amplified influence on phenotypes closer to the gene action. In this study, event-related potential (ERP) phenotypes recorded during a working-memory task were collected from 656 adolescents from 299 families for whom genotypes were available. Univariate linkage analyses using the MERLIN variance-components method were conducted on slow wave phenotypes recorded at multiple sites while participants were required to remember the location of a target. Suggestive linkage (LOD > 2.2) was found on chromosomes 4, 5, 6, 10, 17, and 20. After correcting for multiple testing, suggestive linkage remained on chromosome 10. Empirical thresholds were computed for the most promising phenotypes. Those on chromosome 10 remained suggestive. A number of genes reported to regulate neural differentiation and function (i.e. NRP1, ANK3, and CHAT) were found under these linkage peaks and may influence the levels of neural activity occurring in individuals participating in a spatial working-memory task.

Direct mesh: A multiresolution approach to terrain visualization

Terrain can be approximated by a triangular mesh consisting millions of 3D points. Multiresolution triangular mesh (MTM) structures are designed to support applications that use terrain data at variable levels of detail (LOD). Typically, an MTM adopts a tree structure where a parent node represents a lower-resolution approximation of its descendants. Given a region of interest (ROI) and a LOD, the process of retrieving the required terrain data from the database is to traverse the MTM tree from the root to reach all the nodes satisfying the ROI and LOD conditions. This process, while being commonly used for multiresolution terrain visualization, is inefficient as either a large number of sequential I/O operations or fetching a large amount of extraneous data is incurred. Various spatial indexes have been proposed in the past to address this problem, however level-by-level tree traversal remains a common practice in order to obtain topological information among the retrieved terrain data. A new MTM data structure called direct mesh is proposed. We demonstrate that with direct mesh the amount of data retrieval can be substantially reduced. Comparing with existing MTM indexing methods, a significant performance improvement has been observed for real-life terrain data.

Further evidence of linkage at 7p22 with familial hyperaldosteronism type II

Primary aldosteronism (PAL) is caused by the autonomous over-production of aldosterone. Once thought rare, it is now reported to be responsible for 5–10% of hypertension. Familial hyperaldosteronism type II (FH-II), unlike familial hyperaldosteronism type I, is not glucocorticoid-remediable and not associated with the hybrid CYP11B1/CYP11B2 gene mutation. At least five times more common than FH-I, FH-II is clinically, biochemically and morphologically indistinguishable from apparently sporadic PAL, suggesting that its incidence maybe even higher. Studies performed in collaboration with C Stratakis (NIH, Bethesda) on our largest Australian FH-II family (eight affected members) demonstrated linkage at chromosome 7p22. Similar linkage at this region was also found in a South American FH-II family (DNA provided by MI New, Presbyterian Hospital, New York). Mutations in the exons and intron/exon boundaries of the PRKARIB gene (which resides at 7p22 and is closely related to PRKARIA gene mutated in Carney complex) have been excluded in our largest Australian FH-II family. Using more finely spaced markers, we have confirmed linkage at 7p22 in these 2 families, and identified a second Australian family with evidence of linkage at this locus. The combined multipoint LOD score for these 3 families is 4.87 (θ=0) with markers D7S462 and D7S2424, which exceeds the critical threshold for genome-wide significance suggested by Lander and Kruglyak (1995), providing strong support for this locus harbouring mutations responsible for FH-II. A newly identified recombination event in our largest Australian family has narrowed the region of linkage by 1.8 Mb, permitting exclusion of approximately half the genes residing in the original reported 5Mb linked locus. In addition, we have strongly excluded linkage to these key markers in two Australian families (maximum multipoint LOD scores −3.51 and −2.77), supporting the notion that FH-II may be genetically heterogeneous. In order to identify candidate genes at 7p22, more closely spaced markers will be used to refine the locus, as well as single nucleotide polymorphism analysis.

A D-optimal designed population pharmacokinetic study of itraconazole capsules and solution in adults with cystic fibrosis

Background: Oral itraconazole (ITRA) is used for the treatment of allergic bronchopulmonary aspergillosis in patients with cystic fibrosis (CF) because of its antifungal activity against Aspergillus species. ITRA has an active hydroxy-metabolite (OH-ITRA) which has similar antifungal activity. ITRA is a highly lipophilic drug which is available in two different oral formulations, a capsule and an oral solution. It is reported that the oral solution has a 60% higher relative bioavailability. The influence of altered gastric physiology associated with CF on the pharmacokinetics (PK) of ITRA and its metabolite has not been previously evaluated. Objectives: 1) To estimate the population (pop) PK parameters for ITRA and its active metabolite OH-ITRA including relative bioavailability of the parent after administration of the parent by both capsule and solution and 2) to assess the performance of the optimal design. Methods: The study was a cross-over design in which 30 patients received the capsule on the first occasion and 3 days later the solution formulation. The design was constrained to have a maximum of 4 blood samples per occasion for estimation of the popPK of both ITRA and OH-ITRA. The sampling times for the population model were optimized previously using POPT v.2.0.[1] POPT is a series of applications that run under MATLAB and provide an evaluation of the information matrix for a nonlinear mixed effects model given a particular design. In addition it can be used to optimize the design based on evaluation of the determinant of the information matrix. The model details for the design were based on prior information obtained from the literature, which suggested that ITRA may have either linear or non-linear elimination. The optimal sampling times were evaluated to provide information for both competing models for the parent and metabolite and for both capsule and solution simultaneously. Blood samples were assayed by validated HPLC.[2] PopPK modelling was performed using FOCE with interaction under NONMEM, version 5 (level 1.1; GloboMax LLC, Hanover, MD, USA). The PK of ITRA and OH‑ITRA was modelled simultaneously using ADVAN 5. Subsequently three methods were assessed for modelling concentrations less than the LOD (limit of detection). These methods (corresponding to methods 5, 6 & 4 from Beal[3], respectively) were (a) where all values less than LOD were assigned to half of LOD, (b) where the closest missing value that is less than LOD was assigned to half the LOD and all previous (if during absorption) or subsequent (if during elimination) missing samples were deleted, and (c) where the contribution of the expectation of each missing concentration to the likelihood is estimated. The LOD was 0.04 mg/L. The final model evaluation was performed via bootstrap with re-sampling and a visual predictive check. The optimal design and the sampling windows of the study were evaluated for execution errors and for agreement between the observed and predicted standard errors. Dosing regimens were simulated for the capsules and the oral solution to assess their ability to achieve ITRA target trough concentration (Cmin,ss of 0.5-2 mg/L) or a combined Cmin,ss for ITRA and OH-ITRA above 1.5mg/L. Results and Discussion: A total of 241 blood samples were collected and analysed, 94% of them were taken within the defined optimal sampling windows, of which 31% where taken within 5 min of the exact optimal times. Forty six per cent of the ITRA values and 28% of the OH-ITRA values were below LOD. The entire profile after administration of the capsule for five patients was below LOD and therefore the data from this occasion was omitted from estimation. A 2-compartment model with 1st order absorption and elimination best described ITRA PK, with 1st order metabolism of the parent to OH-ITRA. For ITRA the clearance (ClItra/F) was 31.5 L/h; apparent volumes of central and peripheral compartments were 56.7 L and 2090 L, respectively. Absorption rate constants for capsule (kacap) and solution (kasol) were 0.0315 h-1 and 0.125 h-1, respectively. Comparative bioavailability of the capsule was 0.82. There was no evidence of nonlinearity in the popPK of ITRA. No screened covariate significantly improved the fit to the data. The results of the parameter estimates from the final model were comparable between the different methods for accounting for missing data, (M4,5,6)[3] and provided similar parameter estimates. The prospective application of an optimal design was found to be successful. Due to the sampling windows, most of the samples could be collected within the daily hospital routine, but still at times that were near optimal for estimating the popPK parameters. The final model was one of the potential competing models considered in the original design. The asymptotic standard errors provided by NONMEM for the final model and empirical values from bootstrap were similar in magnitude to those predicted from the Fisher Information matrix associated with the D-optimal design. Simulations from the final model showed that the current dosing regimen of 200 mg twice daily (bd) would provide a target Cmin,ss (0.5-2 mg/L) for only 35% of patients when administered as the solution and 31% when administered as capsules. The optimal dosing schedule was 500mg bd for both formulations. The target success for this dosing regimen was 87% for the solution with an NNT=4 compared to capsules. This means, for every 4 patients treated with the solution one additional patient will achieve a target success compared to capsule but at an additional cost of AUD $220 per day. The therapeutic target however is still doubtful and potential risks of these dosing schedules need to be assessed on an individual basis. Conclusion: A model was developed which described the popPK of ITRA and its main active metabolite OH-ITRA in adult CF after administration of both capsule and solution. The relative bioavailability of ITRA from the capsule was 82% that of the solution, but considerably more variable. To incorporate missing data, using the simple Beal method 5 (using half LOD for all samples below LOD) provided comparable results to the more complex but theoretically better Beal method 4 (integration method). The optimal sparse design performed well for estimation of model parameters and provided a good fit to the data.