A powerful genome-wide feasible approach to detect parent-of-origin effects in studies of quantitative traits

There is currently a lot of interest in the role of genomic imprinting in mammalian development. Many human diseases, such as cancer, obesity, diabetes and behavioral traits, may be related to imprinted genes. When searching for genes related to complex disorders, the power of genome-wide association analysis can be improved by introducing parent-of-origin effects into the analyses. For quantitative traits, family-based TDT analysis has successfully implemented such an approach. Although attractive for several reasons, TDT-based tests are known to be less powerful than methods based on measured genotype approaches. In this study, we describe a fast, powerful method for detecting parent-of-origin effects in studies of quantitative traits using a measured genotype framework. First, for each locus studied, we estimate the probabilities of an allele's parental origin using multipoint haplotype reconstruction. Next, we introduce the parental origin of these alleles as a covariate in regression models during the second step of GRAMMAR, a fast approximation to the measured genotype approach. We show that, compared with a TDT-based analysis, our method has a higher power to detect a locus exhibiting a parent-of-origin effect. Moreover, our method is applicable to a wider range of data, including pedigree structures that are not very informative for TDT. The method gives no false positives in the absence of parent-of-origin effects, under both additive and dominant models. As this method is an extension of the rapid GRAMMAR analysis, it is fast enough to be suitable for genome-wide association scans.

Genetic determinants of circulating sphingolipid concentrations in European populations

Sphingolipids have essential roles as structural components of cell membranes and in cell signalling, and disruption of their metabolism causes several diseases, with diverse neurological, psychiatric, and metabolic consequences. Increasingly, variants within a few of the genes that encode enzymes involved in sphingolipid metabolism are being associated with complex disease phenotypes. Direct experimental evidence supports a role of specific sphingolipid species in several common complex chronic disease processes including atherosclerotic plaque formation, myocardial infarction (MI), cardiomyopathy, pancreatic beta-cell failure, insulin resistance, and type 2 diabetes mellitus. Therefore, sphingolipids represent novel and important intermediate phenotypes for genetic analysis, yet little is known about the major genetic variants that influence their circulating levels in the general population. We performed a genome-wide association study (GWAS) between 318,237 single-nucleotide polymorphisms (SNPs) and levels of circulating sphingomyelin (SM), dihydrosphingomyelin (Dih-SM), ceramide (Cer), and glucosylceramide (GluCer) single lipid species (33 traits); and 43 matched metabolite ratios measured in 4,400 subjects from five diverse European populations. Associated variants (32) in five genomic regions were identified with genome-wide significant corrected p-values ranging down to 9.08x10(-66). The strongest associations were observed in or near 7 genes functionally involved in ceramide biosynthesis and trafficking: SPTLC3, LASS4, SGPP1, ATP10D, and FADS1-3. Variants in 3 loci (ATP10D, FADS3, and SPTLC3) associate with MI in a series of three German MI studies. An additional 70 variants across 23 candidate genes involved in sphingolipid-metabolizing pathways also demonstrate association (p = 10(-4) or less). Circulating concentrations of several key components in sphingolipid metabolism are thus under strong genetic control, and variants in these loci can be tested for a role in the development of common cardiovascular, metabolic, neurological, and psychiatric diseases.

PedStr software for cutting large pedigrees for haplotyping, IBD computation and multipoint linkage analysis

We propose an automatic heuristic algorithm for splitting large pedigrees into fragments of no more than a user-specified bit size. The algorithm specifically aims to split large pedigrees where many close relatives are genotyped and to produce a set of sub-pedigrees for haplotype reconstruction, IBD computation or multipoint linkage analysis with the help of the Lander-Green-Kruglyak algorithm. We demonstrate that a set of overlapping pedigree fragments constructed with the help of our algorithm allows fast and effective haplotype reconstruction and detection of an allele's parental origin. Moreover, we compared pedigree fragments constructed with the help of our algorithm and existing programs PedCut and Jenti for multipoint linkage analysis. Our algorithm demonstrated significantly higher linkage power than the algorithm of Jenti and significantly shorter running time than the algorithm of PedCut. The software package PedStr implementing our algorithms is available at http://mga.bionet.nsc.ru/soft/index.html.

Predicting human height by Victorian and genomic methods

In the Victorian era, Sir Francis Galton showed that 'when dealing with the transmission of stature from parents to children, the average height of the two parents, ... is all we need care to know about them' (1886). One hundred and twenty-two years after Galton's work was published, 54 loci showing strong statistical evidence for association to human height were described, providing us with potential genomic means of human height prediction. In a population-based study of 5748 people, we find that a 54-loci genomic profile explained 4-6% of the sex- and age-adjusted height variance, and had limited ability to discriminate tall/short people, as characterized by the area under the receiver-operating characteristic curve (AUC). In a family-based study of 550 people, with both parents having height measurements, we find that the Galtonian mid-parental prediction method explained 40% of the sex- and age-adjusted height variance, and showed high discriminative accuracy. We have also explored how much variance a genomic profile should explain to reach certain AUC values. For highly heritable traits such as height, we conclude that in applications in which parental phenotypic information is available (eg, medicine), the Victorian Galton's method will long stay unsurpassed, in terms of both discriminative accuracy and costs. For less heritable traits, and in situations in which parental information is not available (eg, forensics), genomic methods may provide an alternative, given that the variants determining an essential proportion of the trait's variation can be identified.

An approach for cutting large and complex pedigrees for linkage analysis

Utilizing large pedigrees in linkage analysis is a computationally challenging task. The pedigree size limits applicability of the Lander-Green-Kruglyak algorithm for linkage analysis. A common solution is to split large pedigrees into smaller computable subunits. We present a pedigree-splitting method that, within a user supplied bit-size limit, identifies subpedigrees having the maximal number of subjects of interest (eg patients) who share a common ancestor. We compare our method with the maximum clique partitioning method using a large and complex human pedigree consisting of 50 patients with Alzheimer's disease ascertained from genetically isolated Dutch population. We show that under a bit-size limit our method can assign more patients to subpedigrees than the clique partitioning method, particularly when splitting deep pedigrees where the subjects of interest are scattered in recent generations and are relatively distantly related via multiple genealogic connections. Our pedigree-splitting algorithm and associated software can facilitate genome-wide linkage scans searching for rare mutations in large pedigrees coming from genetically isolated populations. The software package PedCut implementing our approach is available at http://mga.bionet.nsc.ru/soft/index.html.

Breaking loops in large complex pedigrees

For pedigrees with multiple loops, exact likelihoods could not be computed in an acceptable time frame and thus, approximate methods are used. Some of these methods are based on breaking loops and approximations of complex pedigree likelihoods using the exact likelihood of the corresponding zero-loop pedigree. Due to ignoring loops, this method results in a loss of genetic information and a decrease in the power to detect linkage. To minimize this loss, an optimal set of loop breakers has to be selected. In this paper, we present a graph theory based algorithm for automatic selection of an optimal set of loop breakers. We propose using a total relationship between measured pedigree members as a proxy to power. To minimize the loss of genetic information, we suggest selection of such breakers whose duplication in a pedigree would be accompanied by a minimal loss of total relationship between measured pedigree members. We show that our algorithm compares favorably with other existing loop-breaker selection algorithms in terms of conservation of genetic information, statistical power and CPU time of subsequent linkage analysis. We implemented our method in a software package LOOP_EDGE, which is available at http://mga.bionet.nsc.ru/nlru/.

A genomewide screen for late-onset Alzheimer disease in a genetically isolated Dutch population

Alzheimer disease (AD) is the most common cause of dementia. We conducted a genome screen of 103 patients with late-onset AD who were ascertained as part of the Genetic Research in Isolated Populations (GRIP) program that is conducted in a recently isolated population from the southwestern area of The Netherlands. All patients and their 170 closely related relatives were genotyped using 402 microsatellite markers. Extensive genealogy information was collected, which resulted in an extremely large and complex pedigree of 4,645 members. The pedigree was split into 35 subpedigrees, to reduce the computational burden of linkage analysis. Simulations aiming to evaluate the effect of pedigree splitting on false-positive probabilities showed that a LOD score of 3.64 corresponds to 5% genomewide type I error. Multipoint analysis revealed four significant and one suggestive linkage peaks. The strongest evidence of linkage was found for chromosome 1q21 (heterogeneity LOD [HLOD]=5.20 at marker D1S498). Approximately 30 cM upstream of this locus, we found another peak at 1q25 (HLOD=4.0 at marker D1S218). These two loci are in a previously established linkage region. We also confirmed the AD locus at 10q22-24 (HLOD=4.15 at marker D10S185). There was significant evidence of linkage of AD to chromosome 3q22-24 (HLOD=4.44 at marker D3S1569). For chromosome 11q24-25, there was suggestive evidence of linkage (HLOD=3.29 at marker D11S1320). We next tested for association between cognitive function and 4,173 single-nucleotide polymorphisms in the linked regions in an independent sample consisting of 197 individuals from the GRIP region. After adjusting for multiple testing, we were able to detect significant associations for cognitive function in four of five AD-linked regions, including the new region on chromosome 3q22-24 and regions 1q25, 10q22-24, and 11q25. With use of cognitive function as an endophenotype of AD, our study indicates the that the RGSL2, RALGPS2, and C1orf49 genes are the potential disease-causing genes at 1q25. Our analysis of chromosome 10q22-24 points to the HTR7, MPHOSPH1, and CYP2C cluster. This is the first genomewide screen that showed significant linkage to chromosome 3q23 markers. For this region, our analysis identified the NMNAT3 and CLSTN2 genes. Our findings confirm linkage to chromosome 11q25. We were unable to confirm SORL1; instead, our analysis points to the OPCML and HNT genes.

Optimal peeling order for pedigrees with incomplete genotypic information

The likelihood approach is common in linkage analysis of large extended pedigrees. Various peeling procedures, based on the conditional independence of separate parts of a pedigree, are typically used for likelihood calculations. A peeling order may significantly affect the complexity of such calculations, particularly for pedigrees with loops or when many pedigrees members have unknown genotypes. Several algorithms have been proposed to address this problem for pedigrees with loops. However, the problem has not been solved for pedigrees without loops until now. In this paper, we suggest a new graph theoretic algorithm for optimal selection of peeling order in zero-loop pedigrees with incomplete genotypic information. It is especially useful when multiple likelihood calculation is needed, for example, when genetic parameters are estimated or linkage with multiple marker loci is tested. The algorithm can be easily introduced into the existing software packages for linkage analysis based on the Elston-Stewart algorithm for likelihood calculation. The algorithm was implemented in a software package PedPeel, which is freely available at http://mga.bionet.nsc.ru/nlru/.

Heritability of blood pressure traits and the genetic contribution to blood pressure variance explained by four blood-pressure-related genes

OBJECTIVE

To study the heritability of four blood pressure traits and the proportion of variance explained by four blood-pressure-related genes.

METHODS

All participants are members of an extended pedigree from a Dutch genetically isolated population. Heritability and genetic correlations of systolic blood pressure, diastolic blood pressure, mean arterial pressure and pulse pressure were assessed using a variance components approach (SOLAR). Polymorphisms of the alpha-adducin (ADD1), angiotensinogen (AGT), angiotensin II type 1 receptor (AT1R) and G protein beta3 (GNB3) genes were typed.

RESULTS

Heritability estimates were significant for all four blood pressure traits, ranging between 0.24 and 0.37. Genetic correlations between systolic blood pressure, diastolic blood pressure and mean arterial pressure were high (0.93-0.98), and those between pulse pressure and diastolic blood pressure were low (0.05). The ADD1 polymorphism explained 0.3% of the variance of pulse pressure (P = 0.07), and the polymorphism of GNB3 explained 0.4% of the variance of systolic blood pressure (P = 0.02), 0.2% of mean arterial pressure (P = 0.05) and 0.3% of pulse pressure (P = 0.06).

CONCLUSION

Genetic factors contribute to a substantial proportion of blood pressure variance. In this study, the effect of polymorphisms of ADD1, AGT, AT1R and GNB3 explained a very small proportion of the heritability of blood pressure traits. As new genes associated with blood pressure are localized in the future, their effect on blood pressure variance should be calculated.

Inheritance of litter size at birth in farmed arctic foxes (Alopex lagopus, Canidae, Carnivora)

Natural populations of the arctic fox (Alopex lagopus, Canidae, Carnivora) differ drastically in their reproductive strategy. Coastal foxes, which depend on stable food resources, produce litters of moderate size. Inland foxes feed on small rodents, whose populations are characterized by cycling fluctuation. In the years with low food supply, inland fox populations have a very low rate of reproduction. In the years with high food supply, they undergo a population explosion. To gain insight into the genetic basis of the reproductive strategy of this species, we performed complex segregation analysis of the litter size in the extended pedigree of the farmed arctic foxes involving 20,665 interrelated animals. Complex segregation analysis was performed using a mixed model assuming that the trait was under control of a major gene and a large number of additive genetic and random factors. To check the significance of any major gene effect, we used Elston-Stewart transmission probability test. Our analysis demonstrated that the inheritance of this trait can be described within the frameworks of a major gene model with recessive control of low litter size. This model was also supported by the pattern of its familial segregation and by comparison of the distributions observed in the population and that expected under our model. We suggest that a system of balanced polymorphism for litter size in the farmed population might have been established in natural populations of arctic foxes as a result of adaptation to the drastic fluctuations in prey availability.

Solution for underflow problem in linkage and segregation analysis

Finding genes for complex traits is one of the major challenges of modern human genetics. Current developments of molecular techniques facilitated use of large pedigrees and marker sets of thousands of single-nucleotide polymorphisms (SNPs). However, one of the problems occurring in statistical analysis of such large data sets is that the likelihood is very low and underflow may easily occur. In this work we describe a method permitting to avoid underflow during computation of a likelihood function, using different algorithms. Our method makes practically possible analysis of thousands of individuals and thousands of SNPs. The method is easy to implement without major change of the code of existing programs. It also helps to reduce the amount of computer memory used in analysis without noticeable alteration of the program running time. The algorithm was implemented in the software packages for segregation and linkage analysis, which are available from http://mga.bionet.nsc.ru/.

Inheritance of White Head Spotting in Natural Populations of South American Water Rat (Nectomys squamipes Rodentia: Sigmodontinae)

Specimens with white head spots are present at low frequency in the natural populations of South American water rat (Nectomys squamipes) and absent in the sibling species Nectomys rattus. We analyzed the pattern of inheritance of the phenotype using complex segregation analysis of pedigrees of a captive-bred population of N. squamipes. We found that the inheritance of the white head spot in this species can be described within the framework of the major gene recessive model with incomplete penetrance of genotypes.

miLD and booLD programs for calculation and analysis of corrected linkage disequilibrium

We describe software which calculates a set of linkage disequilibrium statistics, including multiallelic D' corrected by the bootstrap and permutation. The software also provides a tool for maximum likelihood and least squares estimation and testing of a set of hierarchical hypotheses formulated within the framework of the Malecot model of the decay of linkage disequilibrium with distance. Additionally, the bootstrap approach is used for estimation of the model parameter's confidence intervals and for hypothesis testing. The programs are available from http://www.geneticepi.com/Research/software/software.html

Inheritance of litter size at birth in the Brazilian grass mouse (Akodon cursor, Sigmodontinae, Rodentia)

By means of complex segregation analysis we studied the inheritance of litter size in two large pedigrees of captive-bred colonies of the Brazilian grass mouse Akodon cursor. Genetic analysis has revealed a highly significant influence of genetic factors on the variation of litter size (heritability, h2, was estimated as 0.44). The inheritance followed the classical polygene model: neither the major-gene model nor the polygene with unequal contribution model described the data significantly better.

Segregation analysis of Scheuermann disease in ninety families from Siberia

Scheuermann disease [OMIM number 181440] is the most common cause of structural kyphosis in adolescence. Segregation analysis using a model with gender effects was applied to 90 pedigrees from Barnaul (West Siberia, Russia) ascertained through a proband with Scheuermann disease. The transmission probability model was used to detect major gene effect. A significant contribution of a major gene to the control of the pathology was established. Inheritance of the disease can be described within the framework of a dominant major gene diallele model. According to this model, Scheuermann disease should never occur in the absence of the mutant allele. All male carriers of the mutant allele develop the disease, while only a half of female carriers manifest it. We found a high frequency of idiopathic scoliosis in the families with Scheuermann disease (0.08 vs. 0.01-0.02 in general population). We also observed a succession of idiopathic scoliosis and Scheuermann disease in consecutive generations. The familial aggregation of these two spinal pathologies in the present sample may indicate a genetic unity of Scheuermann disease and idiopathic scoliosis.

Sample size required for predefined linkage decision quality

A method for estimating the sample size required to attain a predefined linkage decision quality (type I and type II errors) is proposed using the linkage test power estimate developed by Ginsburg et al. [(1996) Genet Epidemiol 13:355-366]. The method is applicable for samples of arbitrarily structured pedigrees collected via proband. Comparison of different ascertainment schemes and pedigree structures by their consequent minimal sample size was performed. For recessive and dominant inheritance with complete penetrance, the relative ranks of the ascertainment schemes are invariant regardless of the true recombination fraction value and the trait and marker gene frequencies, which enables one to point out the better scheme. The feasibility of evaluating a sampling strategy by the cost of pedigree collection is also considered, and comparison between these two methods of sample planning is performed.

Inheritance of pendulum movements in rats

The stereotyped hyperkinesis referred to as pendulum movements (PM) may be found in up to 50% of the animals in stocks of Wistar rats. The mode of inheritance of predisposition to PM was studied by two methods: 1) a classical Mendelian analysis of hybrids of the strains PM+ and PM- bred from Wistar stock for enhancement and absence of PM, respectively, and 2) a segregation analysis of pedigree data from the archive records of breeding the cataleptic GC strain. The two methods gave the same result: the inheritance of predisposition to PM can be explained by a major gene model with an incomplete penetrance of heterozygous genotype.

Segregation analysis of idiopathic scoliosis: demonstration of a major gene effect

Segregation analysis using a model with age and gender effects was applied to 101 pedigrees ascertained through a proband with idiopathic scoliosis. The transmission probability model was used to detect major gene effect. When we analyzed the pedigrees where affected status was assigned to persons with a Cobb's angle of more than 5 degrees we did not detect a significant major gene effect. However, when the affected status was assigned to persons with pronounced forms of disease only (a curve of at least 11 degrees) a significant contribution of a major causal gene could be established and inheritance could be described according to a dominant major gene diallele model, assuming incomplete sex and age dependent penetrance of genotypes. According to this model, the pronounced forms of idiopathic scoliosis should never occur in the absence of the mutant allele. This indicates that only the carriers of the mutant allele develop pronounced forms of the disease. At the same time, only a fraction of the carriers of the mutant gene should manifest the disease (30% of males and 50% of females).

Inheritance of male hybrid sterility in the house musk shrew (Suncus murinus, Insectivora, Soricidae)

Two geographic races of the house musk shrew (Suncus murinus) were crossed and intercrossed in the laboratory. Many cases of male sterility were detected among the hybrids. Segregation analysis of the pedigree data showed that the inheritance of male sterility in interracial hybrids of S. murinus can be described within the framework of monogene polyallele model with sterility of a single allele combination. This model is similar if not identical to that proposed by Dobzhansky and Muller.

Chromosomal segregation and fertility in Robertsonian chromosomal heterozygotes of the house musk shrew (Suncus murinus, Insectivora, Soricidae)

Crucial to our understanding of chromosomal variation and evolution in mammals are detailed studies of chromosomal heterozygotes, with analyses of chromosomal segregation and chromosome-derived infertility. We studied segregation and fertility in hybrids between karyotypic races of the house musk shrew Suncus murinus. These individuals were heterozygous for up to five Robertsonian fusions (Rbs) and an insertion of heterochromatin in an autosome. All variant chromosomes showed Mendelian segregation and all Rbs segregated independently of each other in the progeny of double heterozygotes. Litter size in single and even multiple Rb heterozygotes was no smaller than that in the less fertile parental strain. The effects of genetic background were more important in determining litter size than Rb heterozygosity for the shrews that we examined.

Inheritance of litter size at birth in the house musk shrew (Suncus murinus, Insectivora, Soricidae)

In this research we estimated the contribution of a major-gene effect to the control of litter size in hybrids between two local populations of the house musk shrew (Suncus murinus). Segregation analysis was performed on the basis of a mixed polygene and major-gene model. The model presumes that two parental populations may differ from each other in gene frequencies and in the values of polygenic effects but not in the major-gene contribution of the trait. Moreover, the peculiarity of the trait--litter size--is taken into account. This trait is not an individual attribute. It characterizes the parental couple and may depend on the genotypes of both parents. Results of segregation analysis of a large hybrid pedigree of Suncus murinus indicate that the parental populations differ in the allele frequency of the major gene (one population is homozygous, while the other contains the two alleles in approximately equal proportions) and in the values of average polygenic effects. Both major-gene and polygenic components are necessary for the correct description of litter size inheritance in interracial hybrids of S murinus, inasmuch as the exclusion of either of them leads to a significant drop in likelihood. The Elston-Stewart criterion also confirms the Mendelian inheritance of the major gene.

Segregation analysis of animal pedigree data from inter-population crosses

The general method of segregation analysis of pedigree data has been developed and widely used in human genetics. We modified this method to examine pedigree data coming from inter-population crosses. These kinds of pedigrees are common in laboratory and farm animal breeding. This paper describes a rationale for the method and illustrates its application to the study of inheritance of litter size and of male sterility in hybrid stock of the house musk shrew (Suncus murinus) derived from crosses of two geographically isolated populations.

On estimation of linkage test power

Two methods of estimating linkage test power are proposed. The first method is based on the maximum likelihood (ML) estimate of the recombination fraction and is intended for use with a likelihood ratio test (LRT) in the form of a chi 2 or lod score. The power is estimated through a noncentral chi 2 distribution with a specially chosen noncentrality parameter. The second method uses the LRT constructed for a simple alternative hypothesis regarding the recombination fraction value. The approximate distribution of this test and a method of estimating its power is proposed. Using simulated pedigree data, the power estimates for these two methods were shown to be satisfactory. Comparisons among these two methods and the computer simulation approach of Boehnke [1986] are performed.