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.

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.

Genetic improvement of lucerne for anthracnose (Colletotrichum trifolii) resistance

Anthracnose, caused by Colletotrichum trifolii, is one of the most serious diseases influencing lucerne persistence and productivity in eastern Australia. The disease is largely controlled by plant resistance; however, new pathotypes of C. trifolii have developed in Australia, seriously limiting the productive life of susceptible cultivars. This paper describes an incompletely recessive and quantitatively inherited resistance to C. trifolii identified in a clone (W116) from cv. Sequel. S-1, F-1, F-2 and backcross populations of W116 and D (highly susceptible clone) were studied for their reaction to C. trifolii race 1. Resistance was found to be quantitatively inherited, and quantitative trait loci associated with resistance and susceptibility were identified in a backcross population (D x W116) x D using random amplified polymorphic DNA and amplified fragment length polymorphic markers. A multi-locus region on linkage group 4 was found to contribute significantly to the resistance phenotype. The application of DNA markers to allow exploitation of this quantitatively inherited resistance in lucerne breeding is discussed.

Inheritance of proportionate dwarfism in Angus cattle

Objective To determine the mode of inheritance of congenital proportionate dwarfism in Angus and Angus crossbred cattle, initially detected in two commercial beef herds in northern New South Wales. Design Matings of normal carrier sires to unrelated cows of diverse breeds, and of one carrier sire to his unaffected daughters. An unrelated Piedmontese bull was also mated to unaffected daughters of the carrier sires. Procedure Two carrier Angus bulls and nine unaffected daughters, all of whom were completely indistinguishable from normal animals, were purchased for controlled breeding studies under known nutritional and disease conditions. Affected and carrier individuals were examined for the presence of obvious chromosomal abnormalities. Results Angus dwarfism has been successfully reproduced under controlled experimental conditions over successive years using unrelated dams and is undoubtedly heritable. The high frequency of occurrence of affected individuals (23/61 = 0.38 +/- .06) among the progeny of matings of the Angus sires to unrelated females of diverse breeding is not compatible with recessive inheritance, because of the negligible frequency of proportionate dwarfism in the breeds of the dams. Both paternal and maternal transmission of the defect was demonstrated, so that imprinting in the strict sense of a gene that is only expressed when received from the male parent appears not to be involved. Tested individuals showed no evidence of gross chromosomal abnormality. Dominant autosomal inheritance with incomplete penetrance was indicated by the lack of expression of the defective gene in the two Angus sires and in three unaffected daughters who produced dwarf calves from matings to the Piedmontese bull. Conclusions The mode of inheritance is that of a single autosomal dominant gene with a penetrance coefficient of 0.75 +/- 0.12, estimated from the observed incidence of 23/61 affected offspring of the two carrier Angus bulls mated to unrelated dams. Simple genetic models involving either (i) an unstable mutant which changes at high frequency to the expressed dominant dwarfing allele during gametogenesis, or (ii) a dominant allele with penetrance determined by an unlinked modifying locus, are shown to be compatible with the experimental data. Both models indicate that penetrance of the dwarfing gene may possibly be higher in matings involving carrier daughters of the two Angus bulls.

Genorne sequence of the thermostable Newcastle disease virus (strain I-2) reveals a possible phenotypic locus

The complete genome sequence of the Australian 1-2 heat-tolerant Newcastle disease virus (NDV) vaccine (master seed stocks) was determined and compared to the sequence of the parent virus from which it had been derived after exposure of the parent stock at 56 degrees C for 30 min. Nucleotide changes were observed at a number of positions with synonymous mutations being greater than those observed for non-synonymous mutations. Sequence data for the HN gene of a parental culture of V4 and two heat-tolerant variants of V4 were obtained. These were compared with the data for the 1-2 viruses and with published sequences for parental V4 and for a number of ND vaccine strains. Sequence analyses did not reveal the ARG 303 deletion in the HN protein, previously claimed to be responsible for the thermostable phenotype. No consistent changes were detected that would indicate involvement of the HN protein in heat resistance. The majority of alterations were observed in the L protein of the virus and it is proposed that these alterations were responsible for the heat-tolerant phenotype of the 1-2 NDV vaccine. (c) 2005 Elsevier B.V. All rights reserved.

Do predators aggregate in response to pest density in agroecosystems? Assessing within-field spatial patterns

1. The spatial heterogeneity of predator populations is an important component of ecological theories pertaining to predator-prey dynamics. Most studies within agricultural fields show spatial correlation (positive or negative) between mean predator numbers and prey abundance across a whole field over time but generally ignore the within-field spatial dimension. We used explicit spatial mapping to determine if generalist predators aggregated within a soybean field, the size of these aggregations and if predator aggregation was associated with pest aggregation, plant damage and predation rate. 2. The study was conducted at Gatton in the Lockyer Valley, 90 km west of Brisbane, Australia. Intensive sampling grids were used to investigate within-field spatial patterns. The first row of each grid was located in a lucerne field (10 m from interface) and the remaining rows were in an adjacent soybean field. At each point on the grid the abundance of foliage-dwelling and ground-dwelling pests and predators was measured, predation rates [using sentinel Helicoverpa armigera (Hubner) egg cards] and plant damage were estimated. Eight grids were sampled across two summer cropping seasons (2000/01, 2001/02). 3. Predators exhibited strong spatial patterning with regions of high and low abundance and activity within what are considered to be uniform soybean fields. Ground-dwelling and foliage-dwelling predators were often aggregated in patches approximately 40 m across. 4. Lycosidae (wolf spiders) displayed aggregation and were consistently more abundant within the lucerne, with a decreasing trap catch with distance from the lucrene/soybean interface. This trend was consistent between subsequent grids in a single field and between fields. 5. The large amount of spatial variability in within-field arthropod abundance (pests and predators) and activity (egg predation and plant damage) indicates that whole field averages were misleading. This result has serious implications for sampling of arthropod abundance and pest management decision-making based on scouting data. 6. There was a great deal of temporal change in the significant spatial patterns observed within a field at each sampling time point during a single season. Predator and pest aggregations observed in these fields were generally not stable for the entire season. 7. Predator aggregation did not correlate consistently with pest aggregation, plant damage or predation rate. Spatial patterns in predator abundance were not associated consistently with any single parameter measured. The most consistent positive association was between foliage-dwelling predators and pests (significant in four of seven grids). Inferring associations between predators and prey based on an intensive one-off sampling grid is difficult, due to the temporal variability in the abundance of each group. 8. Synthesis and applications. This study demonstrated that generalist predator populations are rarely distributed randomly and field edges and adjacent crops can have an influence on within-field predator abundance. This must be considered when estimating arthropod (pest and predator) abundance from a set of samples taken at random locations within a field.

Trait physiology and crop modelling to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. We explore whether a crop growth and development modelling framework can link phenotype complexity to underlying genetic systems in a way that strengthens molecular breeding strategies. We use gene-to-phenotype simulation studies on sorghum to consider the value to marker-assisted selection of intrinsically stable QTLs that might be generated by physiological dissection of complex traits. The consequences on grain yield of genetic variation in four key adaptive traits – phenology, osmotic adjustment, transpiration efficiency, and staygreen – were simulated for a diverse set of environments by placing the known extent of genetic variation in the context of the physiological determinants framework of a crop growth and development model. It was assumed that the three to five genes associated with each trait, had two alleles per locus acting in an additive manner. The effects on average simulated yield, generated by differing combinations of positive alleles for the traits incorporated, varied with environment type. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages with gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies. We simulated a marker-assisted selection (MAS) breeding strategy based on the analyses of gene effects. When marker scores were allocated based on the contribution of gene effects to yield in a single environment, there was a wide divergence in rate of yield gain over all environments with breeding cycle depending on the environment chosen for the QTL analysis. It was suggested that knowledge resulting from trait physiology and modelling would overcome this dependency by identifying stable QTLs. The improved predictive power would increase the utility of the QTLs in MAS. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate QTLs.

Teaching principles of genetics through SNP analysis

An understanding of inheritance requires comprehension of genetic processes at all levels, from molecules to populations. Frequently genetics courses are separated into molecular and organismal genetics and students may fail to see the relationships between them. This is particularly true with human genetics, because of the difficulties in designing experimental approaches which are consistent with ethical restrictions, student abilities and background knowledge, and available time and materials. During 2005 we used analysis of single nucleotide polymorphisms (SNPs) in two genetic regions to enhance student learning and provide a practical experience in human genetics. Students scanned databases to discover SNPs in a gene of interest, used software to design PCR primers and a restriction enzyme based assay for the alleles, and carried out an analysis of the SNP on anonymous individual and family DNAs. The project occupied eight to ten hours per week for one semester, with some time spent in the laboratory and some spent in database searching, reading and writing the report. In completing their projects, students acquired a knowledge of Mendel’s first law (through looking at inheritance patterns), Mendel’s second law and the exceptions (the concepts of linkage and linkage disequilibrium), DNA structure (primer design and restriction enzyme analysis) and function (SNPs in coding and non-coding regions), population genetics and the statistical analysis of allele frequencies, genomics, bioinformatics and the ethical issues associated with the use of human samples. They also developed skills in presentation of results by publication and conference participation. Deficiencies in their understanding (for example of inheritance patterns, gene structure, statistical approaches and report writing) were detected and guidance given during the project. SNP analysis was found to be a powerful approach to enhance and integrate student understanding of genetic concepts.