Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands

CONTEXT People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the etiology of these important psychiatric disorders. OBJECTIVES To undertake a comprehensive genome-wide linkage study of neuroticism using large study samples that have been measured multiple times and to compare the results between countries for replication and across time within countries for consistency. DESIGN Genome-wide linkage scan. SETTING Twin individuals and their family members from Australia and the Netherlands. PARTICIPANTS Nineteen thousand six hundred thirty-five sibling pairs completed self-report questionnaires for neuroticism up to 5 times over a period of up to 22 years. Five thousand sixty-nine sibling pairs were genotyped with microsatellite markers. METHODS Nonparametric linkage analyses were conducted in MERLIN-REGRESS for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time-specific measures of neuroticism from each country to investigate consistency of linkage results. RESULTS Three chromosomal regions exceeded empirically derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 Kosambi cM (cM) (Dutch study sample), 14q 103 cM (Dutch study sample), and 18q 117 cM (combined Australian and Dutch study sample), but only 14q retained significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and 1 (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice. CONCLUSIONS Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbor causal variants for neuroticism and its related psychiatric disorders and can inform prioritization of results from genome-wide association studies.

Genome-wide association studies of quantitative traits with related individuals: Little (power) lost but much to be gained

For complex disease genetics research in human populations, remarkable progress has been made in recent times with the publication of a number of genome-wide association scans (GWAS) and subsequent statistical replications. These studies have identified new genes and pathways implicated in disease, many of which were not known before. Given these early successes, more GWAS are being conducted and planned, both for disease and quantitative phenotypes. Many researchers and clinicians have DNA samples available on collections of families, including both cases and controls. Twin registries around the world have facilitated the collection of large numbers of families, with DNA and multiple quantitative phenotypes collected on twin pairs and their relatives. In the design of a new GWAS with a fixed budget for the number of chips, the question arises whether to include or exclude related individuals. It is commonly believed to be preferable to use unrelated individuals in the first stage of a GWAS because relatives are 'over-matched' for genotypes. In this study, we quantify that for GWAS of a quantitative phenotype, relative to a sample of unrelated individuals surprisingly little power is lost when using relatives. The advantages of using relatives are manifold, including the ability to perform more quality control, the choice to perform within-family tests of association that are robust to population stratification, and the ability to perform joint linkage and association analysis. Therefore, the advantages of using relatives in GWAS for quantitative traits may well outweigh the small disadvantage in terms of statistical power.

Within-family outliers: Segregating alleles or environmental effects? A linkage analysis of height from 5815 sibling pairs

Most information in linkage analysis for quantitative traits comes from pairs of relatives that are phenotypically most discordant or concordant. Confounding this, within-family outliers from non-genetic causes may create false positives and negatives. We investigated the influence of within-family outliers empirically, using one of the largest genome-wide linkage scans for height. The subjects were drawn from Australian twin cohorts consisting of 8447 individuals in 2861 families, providing a total of 5815 possible pairs of siblings in sibships. A variance component linkage analysis was performed, either including or excluding the within-family outliers. Using the entire dataset, the largest LOD scores were on chromosome 15q (LOD 2.3) and 11q (1.5). Excluding within-family outliers increased the LOD score for most regions, but the LOD score on chromosome 15 decreased from 2.3 to 1.2, suggesting that the outliers may create false negatives and false positives, although rare alleles of large effect may also be an explanation. Several regions suggestive of linkage to height were found after removing the outliers, including 1q23.1 (2.0), 3q22.1 (1.9) and 5q32 (2.3). We conclude that the investigation of the effect of within-family outliers, which is usually neglected, should be a standard quality control measure in linkage analysis for complex traits and may reduce the noise for the search of common variants of modest effect size as well as help identify rare variants of large effect and clinical significance. We suggest that the effect of within-family outliers deserves further investigation via theoretical and simulation studies.

Variants in EMX2 and PTEN do not contribute to risk of endometriosis

Endometriosis has a genetic component, and significant linkage has been found to a region on chromosome 10q. Two candidate genes, EMX2 and PTEN, implicated in both endometriosis and endometrial cancer, lie on chromosome 10q. We hypothesized that variation in EMX2 and/or PTEN could contribute to the risk of endometriosis and may account for some of the linkage signal on 10q. We genotyped single nucleotide polymorphisms (SNPs) in a case-control design to evaluate association between endometriosis and common variations in these two genes. The genotyping and statistical analysis were based on samples collected from Australian volunteers. The cases were 768 unrelated women with surgically confirmed endometriosis selected from affected sister pair (ASP) families participating in the Australian Genes behind Endometriosis Study. The controls were 768 female participants in twin studies who, based on screening questions, did not have a diagnosis of endometriosis. Genotypes of 22 SNPs in the EMX2 gene and 15 SNPs in the PTEN gene were the main outcome measures. Statistical analysis provided measures of linkage disequilibrium and association. Permutation testing showed no globally significant association between any SNPs or haplotypes and endometriosis for either gene. It is unlikely that the EMX2 or PTEN gene variants investigated contribute to risk for initiation and/or development of endometriosis.

A case of true hermaphroditism reveals an unusual mechanism of twinning

Traditionally twins are classified as dizygous or fraternal and monozygous or identical (Hall Twinning, 362, 2003 and 735-743). We report a rare case of 46,XX/46,XY twins: Twin A presented with ambiguous genitalia and Twin B was a phenotypically normal male. These twins demonstrate a third, previously unreported mechanism for twinning. The twins underwent initial investigation with 17-hydroxyprogesterone and testosterone levels, pelvic ultrasound and diagnostic laparoscopy. Cytogenetic analysis was performed on peripheral blood cells and skin fibroblasts. Histological examination and Fluorescence in situ hybridization studies on touch imprints were performed on gonadal biopsies. DNA analysis using more than 6,000 DNA markers was performed on skin fibroblast samples from the twins and on peripheral blood samples from both parents. Twin A was determined to be a true hermaphrodite and Twin B an apparently normal male. Both twins had a 46,XX/46,XY chromosome complement in peripheral lymphocytes, skin fibroblasts, and gonadal biopsies. The proportion of XX to XY cells varied between the twins and the tissues evaluated. Most significantly the twins shared 100% of maternal alleles and approximately 50% of paternal alleles in DNA analysis of skin fibroblasts. The twins are chimeric and share a single genetic contribution from their mother but have two genetic contributions from their father thus supporting the existence of a third, previously unreported type of twinning.

Personality, health and lifestyle in a questionnaire family study: A comparison between highly cooperative and less cooperative families

The effect of nonresponse on health and lifestyle measures has received extensive study, showing at most relatively modest effects. Nonresponse bias with respect to personality has been less thoroughly investigated. The present study uses data from responding individuals as a proxy for the missing data of their nonresponding family members to examine the presence of nonresponse bias for personality traits and disorders as well as health and lifestyle traits. We looked at the Big Five personality traits, borderline personality disorder (BPD) features, attention-deficit/hyperactivity disorder, Anger, and several measures of health (Body Mass Index, migraine) and lifestyle (smoking, alcohol use). In general, outcomes tend to be slightly more favorable for individuals from highly cooperative families compared to individuals from less cooperative families. The only significant difference was found for BPD features (p = .001). However, the absolute difference in mean scores is very small, less than 1 point for a scale ranging from 0 to 72. In conclusion, survey data on personality, health and lifestyle are relatively unbiased with respect to nonresponse.

On the probability of dizygotic twins being concordant for two alleles at multiple polymorphic loci

Accurate determination of same-sex twin zygosity is important for medical, scientific and personal reasons. Determination may be based upon questionnaire data, blood group, enzyme isoforms and fetal membrane examination, but assignment of zygosity must ultimately be confirmed by genotypic data. Here methods are reviewed for calculating average probabilities of correctly concluding a twin pair is monozygotic, given they share the same genotypes across all loci for commonly utilized multiplex short tandem repeat (STR) kits.

Migraine with aura and migraine without aura are not distinct entities: further evidence from a large Dutch population study

It is often debated whether migraine with aura (MA) and migraine without aura (MO) are etiologically distinct disorders. A previous study using latent class analysis (LCA) in Australian twins showed no evidence for separate subtypes of MO and MA. The aim of the present study was to replicate these results in a population of Dutch twins and their parents, siblings and partners (N = 10,144). Latent class analysis of International Headache Society (IHS)-based migraine symptoms resulted in the identification of 4 classes: a class of unaffected subjects (class 0), a mild form of nonmigrainous headache (class 1), a moderately severe type of migraine (class 2), typically without neurological symptoms or aura (8% reporting aura symptoms), and a severe type of migraine (class 3), typically with neurological symptoms, and aura symptoms in approximately half of the cases. Given the overlap of neurological symptoms and nonmutual exclusivity of aura symptoms, these results do not support the MO and MA subtypes as being etiologically distinct. The heritability in female twins of migraine based on LCA classification was estimated at .50 (95% confidence intervals [CI] .27 - .59), similar to IHS-based migraine diagnosis (h2 = .49, 95% CI .19-.57). However, using a dichotomous classification (affected-unaffected) decreased heritability for the IHS-based classification (h2 = .33, 95% CI .00-.60), but not the LCA-based classification (h2 = .51, 95% CI .23-.61). Importantly, use of the LCA-based classification increased the number of subjects classified as affected. The heritability of the screening question was similar to more detailed LCA and IHS classifications, suggesting that the screening procedure is an important determining factor in genetic studies of migraine.

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

A comparison of DNA pools constructed following whole genome amplification for two-stage SNP genotyping designs

Genotyping in DNA pools reduces the cost and the time required to complete large genotyping projects. The aim of the present study was to evaluate pooling as part of a strategy for fine mapping in regions of significant linkage. Thirty-nine single nucleotide polymorphisms (SNPs) were analyzed in two genomic DNA pools of 384 individuals each and results compared with data after typing all individuals used in the pools. There were no significant differences using data from either 2 or 8 heterozygous individuals to correct frequency estimates for unequal allelic amplification. After correction, the mean difference between estimates from the genomic pool and individual allele frequencies was .033. A major limitation of the use of DNA pools is the time and effort required to carefully adjust the concentration of each individual DNA sample before mixing aliquots. Pools were also constructed by combining DNA after Multiple Displacement Amplification (MDA). The MDA pools gave similar results to pools constructed after careful DNA quantitation (mean difference from individual genotyping .040) and MDA provides a rapid method to generate pools suitable for some applications. Pools provide a rapid and cost-effective screen to eliminate SNPs that are not polymorphic in a test population and can detect minor allele frequencies as low as 1% in the pooled samples. With current levels of accuracy, pooling is best suited to an initial screen in the SNP validation process that can provide high-throughput comparisons between cases and controls to prioritize SNPs for subsequent individual genotyping.

Genomewide significant linkage to migrainous headache on chromosome 5q21

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

Zygosity diagnosis in the absence of genotypic data: An approach using latent class analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

Genetic and environmental influences on human responses to odors

Olfaction, the sense of smell, has many important functions in humans. Human responses to odors show substantial individual variation. Olfactory receptor genes have been identified and other genes may also influence olfaction. However, the proportion of phenotypic variation in odor response due to genetic variation remains largely unknown. Little is also known about which genes modify specific responses to odors. This study aimed to elucidate genetic and environmental influences on human responses to odors. Individuals from Finnish families (n=146) and Australian (n=413), British (n=163), Danish (n=336), and Finnish (n=399) twins rated intensity and pleasantness of a set of 12 (families) or 6 (twins) odors and tried to identify the odors. In addition, the participants rated their own sense of smell and annoyance experienced with different environmental odors. The odor stimuli of a commercial smell test (The Brief Smell Identification Test; banana, chocolate, cinnamon, gasoline, lemon, onion, paint thinner, pineapple, rose, smoke, soap, and turpentine) were presented in the family study. Based on the results of the family study and a literature survey, a new set of odor stimuli (androstenone, chocolate, cinnamon, isovaleric acid, lemon, and turpentine) was designed for the twin studies. In the family sample, heritabilities of the traits were estimated and underlying genomic regions were searched using a genome-wide linkage scan. In the pooled twin sample, variation in the measured traits was decomposed into genetic and environmental components using quantitative genetic modeling. In addition, associations between nongenetic factors (e.g., sex, age, and smoking) and olfactory-related traits were explored. Suggestive evidence for a genetic linkage for pleasantness of cinnamon at a locus on chromosome 4q32.3 emerged from the family sample. High heritability for the pleasantness of cinnamon was found in the family but not the twin study. Heritability of perceived intensity of androstenone odor was determined to be ~30% in the twin sample. A strong genetic correlation between perceived intensity and pleasantness of androstenone, in the absence of any environmental correlation, indicated that only the genetic correlation explained the phenotypic correlation between the traits (r=-0.27) and that the traits were influenced by an overlapping set of genes. Self-rated olfactory function appeared to reflect the odor annoyance experienced rather than actual olfactory acuity or genetic involvement. Results from nongenetic analyses supported the speculated superiority of females' olfactory abilities, the age-related diminishing of olfactory acuity, and the influences of experience-dependent factors on odor responses. This was the first study to estimate heritabilities and perform linkage screens for individual odors. A genetic effect was detected for only a few responses to specific odors, suggesting the predominance of environmental effects in odor perceptions.

A matter of taste : genetic and environmental influences on responses to sweetness

Diet is a major player in the maintenance of health and onset of many diseases of public health importance. The food choice is known to be largely influenced by sensory preferences. However, in many cases it is unclear whether these preferences and dietary behaviors are innate or acquired. The aim of this thesis work was to study the extent to which the individual differences in dietary responses, especially in liking for sweet taste, are influenced by genetic factors. Several traits measuring the responses to sweetness and other dietary variables were applied in four studies: in British (TwinsUK) and Finnish (FinnTwin12 and FinnTwin16) twin studies and in a Finnish migraine family study. All the subjects were adults and they participated in chemosensory measurements (taste and smell tests) and filled in food behavior questionnaires. Further, it was studied, whether the correlations among the variables are mediated by genetic or environmental factors and where in the genome the genes influencing the heritable traits are located. A study of young adult Finnish twins (FinnTwin16, n=4388) revealed that around 40% of the food use is attributable to genetic factors and that the common, childhood environment does not affect the food use even shortly after moving from the parents home. Both the family study (n=146) and the twin studies (British twins, n=663) showed that around half of the variation in the liking for sweetness is inherited. The same result was obtained both by the chemosensory measurements (heritability 41-49%) and the questionnaire variables (heritability 31-54%). By contrast, the intensity perception of sweetness or the responses to saltiness were not influenced by genetic factors. Further, a locus influencing the use-frequency of sweet foods was identified on chromosome 16p. A closer examination of the relationships among the variables based on 663 British twins revealed that several genetic and environmental correlations exist among the different measures of liking for sweetness. However, these correlations were not very strong (range 0.06-0.55) implying that the instruments used measure slightly different aspects of the phenomenon. In addition, the assessment of the associations among responses to fatty foods, dieting behaviors, and body mass index in twin populations (TwinsUK n=1027 and FinnTwin12 n=299) showed that the dieting behaviors (cognitive restraint, uncontrolled eating, and emotional eating) mediate the relationship between obesity and diet. In conclusion, the work increased the understanding of the background variables of human eating behavior. Genetic effects were shown to underlie the variation of many dietary traits, such as liking for sweet taste, use of sweet foods, and dieting behaviors. However, the responses to salty taste were shown to be mainly determined by environmental factors and thus should more easily be modifiable by dietary education, exposure, and learning than sweet taste preferences. Although additional studies are needed to characterize the genetic element located on chromosome 16 that influences the use-frequency of sweet foods, the results underline the importance of inherited factors on human eating behavior.