Estimating heritability using family-pooled phenotypic and genotypic data: a simulation study applied to aquaculture

Estimating heritability based on individual phenotypic and genotypic measurements can be expensive and labour-intensive in commercial aquaculture breeding. Here, the feasibility of estimating heritability using within-family means of phenotypes and allelic frequencies was investigated. Different numbers of full-sib families and family sizes across ten generations with phenotypic and genotypic information on 10 K SNPs were analysed in ten replicates. Three scenarios, representing differing numbers of pools per family (one, two and five) were considered. The results showed that using one pool per family did not reliably estimate the heritability of family means. Using simulation parameters appropriate for aquaculture, at least 200 families of 60 progeny per family divided equally in two pools per family was required to estimate the heritability of family means effectively. Although application of five pools generated more within- and between- family relationships, it reduced the number of individuals per pool and increased within-family residual variation, hence, decreased the heritability of family means. Moreover, increasing the size of pools resulted in increasing the heritability of family means towards one. In addition, heritability of family mean estimates were higher than family heritabilities obtained from Falconer’s formula due to lower intraclass correlation estimate compared to the coefficient of relationship.Subject terms: Genome evolution     

A common genetic variant in the neurexin superfamily member CNTNAP2 increases familial risk of autism

Autism is a childhood neuropsychiatric disorder that, despite exhibiting high heritability, has largely eluded efforts to identify specific genetic variants underlying its etiology. We performed a two-stage genetic study in which genome-wide linkage and family-based association mapping was followed up by association and replication studies in an independent sample. We identified a common polymorphism in contactin-associated protein-like 2 (CNTNAP2), a member of the neurexin superfamily, that is significantly associated with autism susceptibility. Importantly, the genetic variant displays a parent-of-origin and gender effect recapitulating the inheritance of autism.     

No evidence for significant linkage between bipolar affective disorder and chromosome 18 pericentromeric markers in a large series of multiplex extended pedigrees.

Bipolar affective disorder (BP) is a major neuropsychiatric disorder with high heritability and complex inheritance. Previously reported linkage between BP and DNA markers in the pericentromeric region of chromosome 18, with a parent-of-origin effect (linkage was present in pedigrees with paternal transmission and absent in pedigrees with exclusive maternal inheritance), has been a focus of interest in human genetics. We reexamined the evidence in one of the largest samples reported to date (1,013 genotyped individuals in 53 unilineal multiplex pedigrees), using 10 highly polymorphic markers and a range of parametric and nonparametric analyses. There was no evidence for significant linkage between BP and chromosome 18 pericentromeric markers in the sample as a whole, nor was there evidence for significant parent-of-origin effect (pedigrees with paternal transmission were not differentially linked to the implicated chromosomal region). Two-point LOD scores and single-locus sib-pair results gave some support for suggestive linkage, but this was not substantiated by multilocus analysis, and the results were further tempered by multiple test effects. We conclude that there is no compelling evidence for linkage between BP and chromosome 18 pericentromeric markers in this sample.     

Imprinting detection by extending a regression-based QTL analysis method

We present an extension of a regression-based quantitative-trait linkage analysis method to incorporate parent-of-origin effects. We separately regressed total, paternal, and maternal IBD sharing on traits’ squared sums and differences. We also developed a test for imprinting that indicates whether there is any difference between the paternal and maternal regression coefficients. Since this method treats the identity-by-descent information as the dependent variable that is conditioned on the trait, it can be readily applied to data from complex ascertainment processes. We performed a simulation study to examine the performance of the method. We found that when using empirical critical values, the method shows identical or higher power compared to existing methods for evaluation of parent-of-origin effect in linkage analysis of quantitative traits. Missing parental genotypes increase the type I error rate of the linkage test and decrease the power of the imprinting test. When the major gene has a low heritability, the power of the method decreases considerably, but the statistical tests still perform well. We also applied a permutation algorithm, which ensures the appropriate type I error rate for the test for imprinting. The method was applied to a data from a study of 6 body size related measures and 23 loci on chromosome 7 for 255 nuclear families. Multipoint identities-by-descent (IBD) were obtained using a modification of the SIMWALK 2 program. A parent-of-origin effect consistent with maternal imprinting was suggested at 99.67–111.26 Mb for body mass index, bioelectrical impedance analysis, waist circumference, and leptin concentration. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. An erratum to this article can be found at     

Within-family marker-assisted selection for aquaculture species

A within-family marker-assisted selection scheme was designed for typical aquaculture breeding schemes, where most traits are recorded on sibs of the candidates. Here, sibs of candidates were tested for the trait and genotyped to establish genetic marker effects on the trait. BLUP breeding values were calculated, including information of the markers (MAS) or not (NONMAS). These breeding values were identical for all family members in the NONMAS schemes, but differed between family members in the MAS schemes, making within-family selection possible. MAS had up to twice the total genetic gain of the corresponding NONMAS scheme. MAS was somewhat less effective when heritability increased from 0.06 to 0.12 or when the frequency of the positive allele was < 0.5. The relative efficiency of MAS was higher for schemes with more candidates, because of larger fullsib family sizes. MAS was also more efficient when male:female mating ratio changed from 1:1 to 1:5 or when the QTL explained more of the total genetic variation. Four instead of two markers linked to the QTL increased genetic gain somewhat. There was no significant difference in polygenic genetic gain between MAS and NONMAS for most schemes. The rates of inbreeding were lower for MAS than NON-MAS schemes, because fewer full-sibs were selected by MAS.     

Likelihood Formulation of Parent-of-Origin Effects on Segregation Analysis, Including Ascertainment

We developed a likelihood-based method for testing for parent-of-origin effect in complex diseases. The likelihood formulations model parent-of-origin effect and allow for incorporation of ascertainment, as well as differential male and female ascertainment probabilities. The results based on simulated data indicated that the estimates of parental effect (either maternal or paternal) were biased when ascertainment was ignored or when the wrong ascertainment model was used. The exception was single ascertainment, in which we proved that ignoring ascertainment does not bias the estimation of parental effect, in a simple parent-of-origin model. These results underscore the importance of considering ascertainment models when testing for parent-of-origin effect in complex diseases.     

‘Heritability’ of dispersal propensity in a patchy population

Although dispersal is often considered to be a plastic, condition-dependent trait with low heritability, growing evidence supports medium to high levels of dispersal heritability. Obtaining unbiased estimates of dispersal heritability in natural populations nevertheless remains crucial to understand the evolution of dispersal strategies and their population consequences. Here we show that dispersal propensity (i.e. the probability of dispersal between habitat patches) displays a significant heritability in the collared flycatcher Ficedula albicollis, as estimated by within-family resemblance when accounting for environmental factors. Offspring of dispersing mothers or fathers had a higher propensity to disperse to a new habitat patch themselves. The effect of parental dispersal status was additional to that of local habitat quality, as measured by local breeding population size and success, confirming previous results about condition-dependent dispersal in this population. The estimated levels of heritability varied between 0.30±0.07 and 0.47±0.10, depending on parent–offspring comparisons made and correcting for a significant assortative mating with respect to dispersal status. Siblings also displayed a significant resemblance in dispersal propensity. These results suggest that variation in between-patch natal dispersal in the collared flycatcher is partly genetically determined, and we discuss ways to quantify this genetic basis and its implications.     

No Evidence for Enrichment in Schizophrenia for Common Allelic Associations at Imprinted Loci

Most genetic studies assume that the function of a genetic variant is independent of the parent from which it is inherited, but this is not always true. The best known example of parent-of-origin effects arises with respect to alleles at imprinted loci. In classical imprinting, characteristically, either the maternal or paternal copy is expressed, but not both. Only alleles present in one of the parental copies of the gene, the expressed copy, is likely to contribute to disease. It has been postulated that imprinting is important in central nervous system development, and that consequently, imprinted loci may be involved in schizophrenia. If this is true, allowing for parent-of-origin effects might be important in genetic studies of schizophrenia. Here, we use genome-wide association data from one of the world’s largest samples (N = 695) of parent schizophrenia-offspring trios to test for parent-of-origin effects. To maximise power, we restricted our analyses to test two main hypotheses. If imprinting plays a disproportionate role in schizophrenia susceptibility, we postulated a) that alleles showing robust evidence for association to schizophrenia from previous genome-wide association studies should be enriched for parent-of-origin effects and b) that genes at loci imprinted in humans or mice should be enriched both for genome-wide significant associations, and in our sample, for parent-of-origin effects. Neither prediction was supported in the present study. We have shown, that it is unlikely that parent-of-origin effects or imprinting play particularly important roles in schizophrenia, although our findings do not exclude such effects at specific loci nor do they exclude such effects among rare alleles.     

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within-family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.     

Testing the imprinted brain: parent-of-origin effects on empathy and systemizing

Genomic imprinting is a violation of Mendel's laws that enables selection to act on genes, depending on parent of origin, but, even more controversially, on the sex of the offspring. This study tested whether there are parent-of-origin effects on the heritability of empathy and systemizing in the general population as part of a larger question concerning the role of imprinted genes in the evolution of human cognition and behaviour. The measures tested were the Empathy and Systemizing Quotients as proxies for the related terms mentalistic and mechanistic cognition in the imprinted brain theory.To test genomic imprinting hypotheses, correlations in behavioural scores between pairs of full, maternal and paternal siblings were compared. Where scores are influenced by imprinted genes, the actual correlations between pairs of siblings will differ from those expected following classical Mendelian inheritance in a predictable way depending on what kind of imprinting is influencing the trait. These theoretical predictions were used to test the fit of the data against Mendelian and imprinting models using structural equation modeling. The imprinted brain theory proposes a trade-off between maternally influenced mentalistic cognition and paternally influenced mechanistic cognition. However, the results of this study support a model of contrasting maternal and paternal influences on strong and weak empathizing and a maternal influence on systemizing. Although the sample size was insufficient to comprehensively analyse sex-limitation models, there is some evidence that heritability of systemizing is stronger in females than in males.     

Detection of parent-of-origin effects based on complete and incomplete nuclear families with multiple affected children

Parent-of-origin effects are important in studying genetic traits. More than 1% of all mammalian genes are believed to show parent-of-origin effects. Some statistical methods may be ineffective or fail to detect linkage or association for a gene with parent-of-origin effects. Based on case-parents trios, the parental-asymmetry test (PAT) is simple and powerful in detecting parent-of-origin effects. However, it is common in practice to collect nuclear families with both parents as well as nuclear families with only one parent. In this paper, when only one parent is available for each family with an arbitrary number of affected children, we firstly develop a new test statistic 1-PAT to test for parent-of-origin effects in the presence of association between an allele at the marker locus under study and a disease gene. Then we extend the PAT to accommodate complete nuclear families each with one or more affected children. Combining families with both parents and families with only one parent, the C-PAT is proposed to detect parent-of-origin effects. The validity of the test statistics is verified by simulation in various scenarios of parameter values. A power study shows that using the additional information from incomplete nuclear families in the analysis greatly improves the power of the tests, compared to that based on only complete nuclear families. Also, utilizing all affected children in each family, the proposed tests have a higher power than when only one affected child from each family is selected. Additional power comparison also demonstrates that the C-PAT is more powerful than a number of other tests for detecting parent-of-origin effects.     

A Bayesian method for simultaneously detecting Mendelian and imprinted quantitative trait loci in experimental crosses of outbred species

Genomic imprinting is interpreted as a phenomenon, in which some genes inherited from one parent are not completely expressed due to modification of the genome caused during gametogenesis. Subsequently, the expression level of an allele at the imprinted gene is changed dependent on the parental origin, which is referred to as the parent-of-origin effect. In livestock, some QTL for reproductive performance and meat productivity have been reported to be imprinted. So far, methods detecting imprinted QTL have been proposed on the basis of interval mapping, where only a single QTL was tested at a time. In this study, we developed a Bayesian method for simultaneously mapping multiple QTL, allowing the inference about expression modes of QTL in an outbred F2 family. The inference about whether a QTL is Mendelian or imprinted was made using Markov chain Monte Carlo estimation by comparing the goodness-of-fits between models, assuming the presence and the absence of parent-of-origin effect at a QTL. We showed by the analyses of simulated data sets that the Bayesian method can effectively detect both Mendelian QTL and imprinted QTL.     

Linkage analysis of alcoholism-related electrophysiological phenotypes: genome scans with microsatellites compared to single-nucleotide polymorphisms

P300 amplitude is an electrophysiological quantitative trait that is correlated with both alcoholism and smoking status. Using the Collaborative Study on the Genetics of Alcoholism data, we performed model-free linkage analysis to investigate the relationship between alcoholism, P300 amplitude, and habitual smoking. We also analyzed the effect of parent-of-origin on alcoholism, and utilized both microsatellites (MS) markers and single-nucleotide polymorphisms (SNPs). We found significant evidence of linkage for alcoholism to chromosome 10; inclusion of P300 amplitude as a covariate provided additional evidence of linkage to chromosome 12. This same region on chromosome 12 showed some evidence for a parent-of-origin effect. We found evidence of linkage for the P300 phenotype to chromosome 7 in non-smokers, and to chromosome 17 in alcoholics. The effects of alcoholism and habitual smoking on P300 amplitude appear to have separate genetic determinants. Overall, there were few differences between MS and SNP genome scans. The use of covariates and parent-of-origin effects allowed detection of linkage not seen otherwise.     

Possible genomic imprinting of three human obesity-related genetic loci

To detect potentially imprinted, obesity-related genetic loci, we performed genomewide parent-of-origin linkage analyses under an allele-sharing model for discrete traits and under a family regression model for obesity-related quantitative traits, using a European American sample of 1,297 individuals from 260 families, with 391 microsatellite markers. We also used two smaller, independent samples for replication (a sample of 370 German individuals from 89 families and a sample of 277 African American individuals from 52 families). For discrete-trait analysis, we found evidence for a maternal effect in chromosome region 10p12 across the three samples, with LOD scores of 5.69 (single-point) and 4.52 (multipoint) for the pooled sample. For quantitative-trait analysis, we found the strongest evidence for a maternal effect (single-point LOD of 2.85; multipoint LOD of 4.01 for body mass index [BMI] and 3.69 for waist circumference) in region 12q24 and for a paternal effect (single-point LOD of 4.79; multipoint LOD of 3.72 for BMI) in region 13q32, in the European American sample. The results suggest that parent-of-origin effects, perhaps including genomic imprinting, may play a role in human obesity.     

Path analysis of palmar ridge counts.

The methods for path analysis of family resemblance (Rao et al., '74) are employed to test hypotheses concerning the inheritance of a-b, b-c and c-d palmar ridge counts using the correlation data of Pateria ('74). Homogeneity chi-square tests of the various familial correlations provide no evidence for sex-linkage of either kind, and also suggest that maternal effects are absent. The path coefficient model employed here involves heritability (additive) and common sibling environment. Variance components show that both heritability and common environment are significant, and account for most of the variation at each of the three ridge count area; b-c has the highest heritability, significantly higher than that for a-b or c-d.     

Phylogeny, evolutionary history, and biogeography of Oriental-Australian rear-fanged water snakes (Colubroidea: Homalopsidae) inferred from mitochondrial and nuclear DNA sequences

Homalopsid snakes are widely distributed throughout Southeast Asia and form the ecologically dominant component of the herpetofauna over much of their range. Although they are considered well differentiated from other colubrid lineages, several aspects of their radiation including within-family relationships, temporal patterns of species diversification, and biogeographic history remain under studied. We analyzed sequence data from four genes (three mitochondrial and one nuclear) for 22 species of the Homalopsidae to generate the most comprehensive phylogeny of the family to date. We also estimated divergence times within the family using a model of independent but log-normally distributed rates of evolution in conjunction with two external fossil calibrations. Using this chronogram, we inferred historical patterns of species diversification within the family. Finally, we used previously published sequence data for 172 snake species to test for the monophyly of the Homalopsidae. Phylogenetic analysis reveals strong support for homalopsid monophyly with an estimate age of the crown group of approximately 22 MYA. The family comprises three major clades which all originated 18-20 MY. Lineage through time plots reveal that homalopsids experienced a significantly higher rate of effective cladogenesis in their early history, consistent with a hypothesis of adaptive radiation. We discuss several Miocene and Pliocene paleogeographic factors that might underlie observed patterns of temporal diversification and biogeography.     

Quantitative trait loci for foliar terpenes in a global eucalypt species

Terpenes are a diverse group of plant secondary metabolites that mediate a plethora of ecological interactions in many plant species. Despite increasing research into the genetic control of important adaptive traits in some plant species, the genetic control of terpenes in forest tree species is still relatively poorly studied. In this study, we use quantitative genetic and quantitative trait loci (QTL) analysis to investigate the genetic control of foliar terpenes in an ecologically and commercially important eucalypt species, Eucalyptus globulus. We show a moderate to high within-family broad-sense heritability and significant genetic basis to the variation in 14 of the 16 terpenes assayed. This is the first report of QTL for terpenes in this species. Eleven QTL influenced the terpenes overall. One QTL on linkage group 6 affected six of the seven different sesquiterpenes assayed (plus one monoterpene), which, in combination with highly significant correlations between these compounds, argues that their variation is influenced by a QTL with pleiotropic effect early in the biosynthetic pathway. We examine the homology of these QTL to those found in a closely related eucalypt, Eucalyptus nitens, and provide evidence that both common and unique QTL influence terpene levels.     

Kin selection and evolution of infectious disease resistance

Discoveries of mutations conferring resistance to infectious diseases have led to increased interest in the evolutionary dynamics of disease resistance. Several recent papers have estimated the historical strength of selection for mutations conferring disease resistance. These studies are based on simple population genetic models that do not take account of factors such as spatial and family structure. Such factors may have a substantial impact on the strength of natural selection through inclusive fitness effects. That is, people have a strong tendency to live with relatives and therefore have a high probability of transmitting infectious diseases to them. Thus, an allele that protects an individual against disease infection also protects that individual's family members. Because some of these family members are likely to also be carrying the allele, selection for that allele is magnified by family structure. In this paper, I use mathematical modeling techniques to explore the impact of such kin selection on the strength of selection for infectious disease resistance alleles. I show that if the resistance allele has the same proportional effect on both within- and between-family transmission, then the impact of kin selection is relatively minor. Selection coefficients are increased by 5-35%, with a greater benefit for weaker alleles. The reason is that an individual with a strong resistance allele does not need much protection from infection by family members and thus does not benefit much from their alleles. The effect of kin selection can be dramatic, however, if the resistance allele has a larger effect on between-family transmission than within-family transmission (which can occur if between-family infection rates are much smaller than within-family rates), increasing selection coefficients by as much as two- to threefold. These results show conditions when it is important to consider family structure in estimates of the strength of selection for infectious disease resistance alleles.     

maternally expressed gene1 Is a novel maize endosperm transfer cell-specific gene with a maternal parent-of-origin pattern of expression

Growth of the maize (Zea mays) endosperm is tightly regulated by maternal zygotic and sporophytic genes, some of which are subject to a parent-of-origin effect. We report here a novel gene, maternally expressed gene1 (meg1), which shows a maternal parent-of-origin expression pattern during early stages of endosperm development but biallelic expression at later stages. Interestingly, a stable reporter fusion containing the meg1 promoter exhibits a similar pattern of expression. meg1 is exclusively expressed in the basal transfer region of the endosperm. Further, we show that the putatively processed MEG1 protein is glycosylated and subsequently localized to the labyrinthine ingrowths of the transfer cell walls. Hence, the discovery of a parent-of-origin gene expressed solely in the basal transfer region opens the door to epigenetic mechanisms operating in the endosperm to regulate certain aspects of nutrient trafficking from the maternal tissue into the developing seed.