Rare susceptibility variants for bipolar disorder suggest a role for G protein-coupled receptors.

Bipolar disorder (BD) is a prevalent mood disorder that tends to cluster in families. Despite high heritability estimates, few genetic susceptibility factors have been identified over decades of genetic research. One possible interpretation for the shortcomings of previous studies to detect causative genes is that BD is caused by highly penetrant rare variants in many genes. We explored this hypothesis by sequencing the exomes of affected individuals from 40 well-characterized multiplex families. We identified rare variants segregating with affected status in many interesting genes, and found an enrichment of deleterious variants in G protein-coupled receptor (GPCR) family genes, which are important drug targets. Furthermore, we showed targeted downstream GPCR dysregulation for some of the variants that may contribute to disease pathology. Particularly interesting was the finding of a rare and functionally relevant nonsense mutation in the corticotropin-releasing hormone receptor 2 (CRHR2) gene that tracked with affected status in one family. By focusing on rare variants in informative families, we identified key biochemical pathways likely implicated in this complex disorder.

Stability of the human sperm DNA methylome to folic acid fortification and short-term supplementation.

STUDY QUESTION: Do short-term and long-term exposures to low-dose folic acid supplementation alter DNA methylation in sperm? SUMMARY ANSWER: No alterations in sperm DNA methylation patterns were found following the administration of low-dose folic acid supplements of 400 µg/day for 90 days (short-term exposure) or when pre-fortification of food with folic acid and post-fortification sperm samples (long-term exposure) were compared. WHAT IS KNOWN ALREADY: Excess dietary folate may be detrimental to health and DNA methylation profiles due to folate's role in one-carbon metabolism and the formation of S-adenosyl methionine, the universal methyl donor. DNA methylation patterns are established in developing male germ cells and have been suggested to be affected by high-dose (5 mg/day) folic acid supplementation. STUDY DESIGN, SIZE, DURATION: This is a control versus treatment study where genome-wide sperm DNA methylation patterns were examined prior to fortification of food (1996-1997) in men with no history of infertility at baseline and following 90-day exposure to placebo (n = 9) or supplement containing 400 µg folic acid/day (n = 10). Additionally, pre-fortification sperm DNA methylation profiles (n = 19) were compared with those of a group of post-fortification (post-2004) men (n = 8) who had been exposed for several years to dietary folic acid fortification. PARTICIPANTS/MATERIALS, SETTING, METHODS: Blood and seminal plasma folate levels were measured in participants before and following the 90-day treatment with placebo or supplement. Sperm DNA methylation was assessed using the whole-genome and genome-wide techniques, MassArray epityper, restriction landmark genomic scanning, methyl-CpG immunoprecipitation and Illumina HumanMethylation450 Bead Array. MAIN RESULTS AND THE ROLE OF CHANCE: Following treatment, supplemented individuals had significantly higher levels of blood and seminal plasma folates compared to placebo. Initial first-generation genome-wide analyses of sperm DNA methylation showed little evidence of changes when comparing pre- and post-treatment samples. With Illumina HumanMethylation450 BeadChip arrays, no significant changes were observed in individual probes following low-level supplementation; when compared with those of the post-fortification cohort, there were also few differences in methylation despite exposure to years of fortified foods. LARGE SCALE DATA: Illumina HumanMethylation450 BeadChip data from this study have been submitted to the NCBI Gene Expression Omnibus under the accession number GSE89781. LIMITATIONS, REASONS FOR CAUTION: This study was limited to the number of participants available in each cohort, in particular those who were not exposed to early (pre-1998) fortification of food with folic acid. While genome-wide DNA methylation was assessed with several techniques that targeted genic and CpG-rich regions, intergenic regions were less well interrogated. WIDER IMPLICATIONS OF THE FINDINGS: Overall, our findings provide evidence that short-term exposure to low-dose folic acid supplements of 400 µg/day, over a period of 3 months, a duration of time that might occur during infertility treatments, has no major impact on the sperm DNA methylome. STUDY FUNDING/COMPETING INTERESTS: This work was supported by a grant to J.M.T. from the Canadian Institutes of Health Research (CIHR: MOP-89944). The authors have no conflicts of interest to declare.

Opioid receptor mu 1 gene, fat intake and obesity in adolescence.

Dietary preference for fat may increase risk for obesity. It is a complex behavior regulated in part by the amygdala, a brain structure involved in reward processing and food behavior, and modulated by genetic factors. Here, we conducted a genome-wide association study (GWAS) to search for gene loci associated with dietary intake of fat, and we tested whether these loci are also associated with adiposity and amygdala volume. We studied 598 adolescents (12-18 years) recruited from the French-Canadian founder population and genotyped them with 530 011 single-nucleotide polymorphisms. Fat intake was assessed with a 24-hour food recall. Adiposity was examined with anthropometry and bioimpedance. Amygdala volume was measured by magnetic resonance imaging. GWAS identified a locus of fat intake in the µ-opioid receptor gene (OPRM1, rs2281617, P=5.2 × 10(-6)), which encodes a receptor expressed in the brain-reward system and shown previously to modulate fat preference in animals. The minor OPRM1 allele appeared to have a 'protective' effect: it was associated with lower fat intake (by 4%) and lower body-fat mass (by ∼2 kg, P=0.02). Consistent with the possible amygdala-mediated inhibition of fat preference, this allele was additionally associated with higher amygdala volume (by 69 mm(3), P=0.02) and, in the carriers of this allele, amygdala volume correlated inversely with fat intake (P=0.02). Finally, OPRM1 was associated with fat intake in an independent sample of 490 young adults. In summary, OPRM1 may modulate dietary intake of fat and hence risk for obesity, and this effect may be modulated by subtle variations in the amygdala volume.

Association of apolipoprotein E polymorphisms and dietary factors in colorectal cancer.

ApoE single nucleotide polymorphisms (SNPs) Cys112Arg (Epsilon-4), and Arg158Cys (Epsilon-2) have been implicated in cardiovascular and Alzheimer's disease, but their role in colorectal cancer (CRC) has not been extensively studied. We investigated whether ApoE polymorphisms alone or in combination with dietary factors selectively contribute to mismatch-repair (MMR) proficient (microsatellite stable/low or MSS/L) vs deficient (microsatellite unstable or MSI-H) CRCs. We carried out a case-control study with 906 CRC cases and 911 unaffected controls to examine the associations between ApoE polymorphisms and dietary factors and assessed their contribution to MSS/L and MSI-H CRCs. We used unconditional logistic regression to evaluate the associations between ApoE SNPs, tumour MSI status, and dietary factors after adjusting for age and sex. All statistical tests were two-sided. No significant differences in ApoE genotype frequencies were observed between CRC cases and unaffected controls. We observed that increased dietary intake of total fat, saturated fat, cholesterol, and red meat was significantly associated with CRC. Among non-ApoE4 carriers, 2-4 and >4 red meat servings/week were associated with developing MSS/L CRC (OR=1.51, 95% CI 1.10-2.07 and OR=1.80, 95% CI 1.30-2.48, respectively), whereas among ApoE4 allele carriers, four or more red meat servings/week were associated with MSI-H CRC (OR=4.62, 95% CI 1.20-17.77) when compared with the controls. ApoE isoforms modulate the risk of MSI-H and MSS/L CRCs among high red meat consumers.

Association of LY9 in UK and Canadian SLE families.

Systemic lupus erythematosus (SLE) is a complex disease trait of unknown aetiology. Genome-wide linkage studies in human SLE identified several linkage regions, including one at 1q23, which contains multiple susceptibility genes, including the members of the signalling lymphocyte activation molecule (SLAM) locus. In mice there is a syntenic linkage region, Sle1. The SLAM genes are functionally related cell-surface receptors, which regulate signal transduction of cells in the immune system. Family-based association study in UK and Canadian SLE families identified variants in the promoter and coding region of SLAMF7 and LY9 contributing to SLE disease susceptibility. The strongest association was from rs509749, in exon 8 of LY9 (P=0.00209). rs509749 encodes a Val/Met nonsynonymous change in amino acid 602 in the cytoplasmic domain of LY9. In the parents and affected individuals from the Canadian SLE families, the risk allele of rs509049 skews the T-cell population by increasing the number of CD8+ memory T cells, while decreasing the proportion of CD4+ naïve T cells and activated T cells. Since rs509749 lies within the consensus binding site for SAP/SH2D1a, which influences downstream signalling events from LY9, the mechanism for increased CD8+ memory T cells may include differential binding SAP/SH2D1a to the cytoplasmic domain of LY9.

A tree-based model for allele-sharing-based linkage analysis in human complex diseases.

By adapting a well-known affected-relative-pair linkage model that can incorporate covariate or sub-phenotype information [Olson, 1999: Am J Hum Genet 65:1760-1769], we have developed a recursive-partitioning (RP) algorithm (tree-based model) for identifying phenotype and covariate groupings that interact with the evidence for linkage. This strategy is designed to identify subgroups of affected relative pairs demonstrating increased evidence for linkage, where subgroups are defined by pair-level or family-level covariates. After growing a full tree, we identified optimal tree size through a form of tree pruning and chose the best covariate at each split by using bootstrap algorithms. Simulation studies showed that power to detect linkage can increase in the presence of gene-environment interactions, depending on the magnitude of the interaction. As expected, however, power can decrease by examining more covariates, despite the pruning to optimize tree size. The RP model correctly identifies tree structure in a large proportion of simulations. We applied the RP model to a dataset of families with bipolar affective disorder (BPAD) where linkage regions on chromosome 18 have been previously identified. Using the all-pairs score in Genehunter, the NPL tests showed no regions with strong linkage evidence on chromosome 18. However, using the RP model, several suggestive regions were found on chromosome 18. Two covariates appeared to influence the degree of linkage: the type II BPAD subtype and a pattern of displaying mania before or after a depressive episode. The RP model has the potential to identify previously unknown gene-environmental interactions; here we have demonstrated the practical utility and potential this new methodology holds.

Variants of the SFTPA1 and SFTPA2 genes and susceptibility to tuberculosis in Ethiopia.

Lungs are the central organ affected and targeted by Mycobacterium tuberculosis and immune processes in the lung are of critical importance in the pathogenesis of tuberculosis. A major lung defense against invading pathogens is provided by surfactant protein A, a multi-chain protein encoded by the SFTPA1 and SFTPA2 genes. Here, we investigated polymorphisms in the SFTPA1 and SFTPA2 genes for association with tuberculosis in 181 Ethiopian families comprising 226 tuberculosis cases. Four polymorphisms, SFTPA1 307A, SFTPA1 776T, SFTPA2 355C, and SFTPA2 751C, were associated with tuberculosis (P=0.00008; P=0.019, P=0.029 and P=0.042, respectively). Additional subgroup analysis in male, female and more severely affected patients provided evidence for SFTPA1/2-covariate interaction. Finally, out of five intragenic haplotypes identified in the SFTPA1 gene and nine identified in the SFTPA2 gene, 1A(3) was most significantly associated with tuberculosis susceptibility (P=0.026). These findings suggest that SFTPA1 and SFTPA2 modify the risk of tuberculosis susceptibility and that this risk is influenced by additional covariates.

The HNPCC associated MSH2*1906G-->C founder mutation probably originated between 1440 CE and 1715 CE in the Ashkenazi Jewish population.

The MSH2*1906G-->C mutation was recently shown to be a rare yet highly penetrant mutation leading to colorectal cancer. The mutation was only found among Ashkenazi Jewish individuals and lies on an extended haplotype that is common in that population. This study determined that the mutation probably arose between 11 and 22 generations ago, during the time when the Ashkenazim were living in eastern Europe.

Elevated rates of schizophrenia in a familial sample with mental illness and intellectual disability.

BACKGROUND: It is unknown whether intellectual disability (ID) is more familially related to psychotic mood disorders or schizophrenia. L. S. Penrose's large sample of families with two or more members admitted to psychiatric hospitals provided a unique opportunity to investigate the familial relationship between mild ID, schizophrenia and psychotic affective disorders. METHOD: There were 183 affected relative pairs comprising probands with mild ID (95 male, 88 female) and their first or second degree relatives with schizophrenia or psychotic affective disorder. RESULTS: There were nearly twice as many relatives with a diagnosis of schizophrenia (n = 121) as relatives with affective disorders (n = 62) among the intellectually impaired probands. This excess of schizophrenia was statistically significant, even after accounting for the increased risk of hospitalization for schizophrenia (P = 0.005), and was fairly constant across the different relative types. First-degree relatives with either mental illness were more likely to be parents (n = 77) than siblings (n = 51) or children (n = 3), but there was no excess of mother-son pairs. CONCLUSIONS: These results suggest a stronger familial relationship of ID with schizophrenia than psychotic affective disorder, and lend some support to the neurodevelopmental hypothesis of schizophrenia.

The founder mutation MSH2*1906G-->C is an important cause of hereditary nonpolyposis colorectal cancer in the Ashkenazi Jewish population.

Hereditary nonpolyposis colorectal cancer (HNPCC) is caused by mutations in the mismatch-repair genes. We report here the identification and characterization of a founder mutation in MSH2 in the Ashkenazi Jewish population. We identified a nucleotide substitution, MSH2*1906G-->C, which results in a substitution of proline for alanine at codon 636 in the MSH2 protein. This allele was identified in 15 unrelated Ashkenazi Jewish families with HNPCC, most of which meet the Amsterdam criteria. Genotype analysis of 18 polymorphic loci within and flanking MSH2 suggested a single origin for the mutation. All colorectal cancers tested showed microsatellite instability and absence of MSH2 protein, by immunohistochemical analysis. In an analysis of a population-based incident series of 686 Ashkenazi Jews from Israel who have colorectal cancer, we identified 3 (0.44%) mutation carriers. Persons with a family history of colorectal or endometrial cancer were more likely to carry the mutation than were those without such a family history (P=.042), and those with colorectal cancer who carried the mutation were, on average, younger than affected individuals who did not carry it (P=.033). The mutation was not detected in either 566 unaffected Ashkenazi Jews from Israel or 1,022 control individuals from New York. In hospital-based series, the 1906C allele was identified in 5/463 Ashkenazi Jews with colorectal cancer, in 2/197 with endometrial cancer, and in 0/83 with ovarian cancer. When families identified by family history and in case series are included, 25 apparently unrelated Ashkenazi Jewish families have been found to harbor this mutation. Although this pathogenic mutation is not frequent in the Ashkenazi Jewish population (accounting for 2%-3% of colorectal cancer in those whose age at diagnosis is <60 years), it is highly penetrant and accounts for approximately one-third of HNPCC in Ashkenazi Jewish families that fulfill the Amsterdam criteria.

Imprinting and deviation from Mendelian transmission ratios.

Deviations from a Mendelian 1:1 transmission ratio have been observed in human and mouse chromosomes. With few exceptions, the underlying mechanism of the transmission-ratio distortion remains obscure. We tested a hypothesis that grandparental-origin dependent transmission-ratio distortion is related to imprinting and possibly results from the loss of embryos which carry imprinted genes with imprinting marks that have been incorrectly reset. We analyzed transmission of alleles in four regions of the human genome that carry imprinted genes presumably critical for normal embryonic growth and development: 11p15.5 (H19, IGF2, HASH2, etc.), 11p13 (WT1), 7p11-12 (GRB10), and 6q25-q27 (IGF2R), among the offspring of 31 three-generation Centre d'Etude de polymorphism Humain (CEPH) families. Deviations from expected 1:1 ratios were found in the maternal chromosomes for regions 11p15.5, 11p13, and 6q25-27 and in the paternal chromosomes for regions 11p15 and 7p11-p12. The likelihood of the results was assessed empirically to be statistically significant (p = 0.0008), suggesting that the transmission ratios in the imprinted regions significantly deviated from 1:1. We did not find deviations from a 1:1 transmission ratio in imprinted regions that are not crucial for embryo viability (13q14 and 15q11-q13). The analysis of a larger set of 51 families for the 11p15.5 region suggests that there is heterogeneity among the families with regard to the transmission of 11p15.5 alleles. The results of this study are consistent with the hypothesis that grandparental-origin dependent transmission-ratio distortion is related to imprinting and embryo loss.

Design considerations for association studies of candidate genes in families.

In genetic epidemiologic studies, investigators often use generalized linear models to evaluate the relationships between a disease trait and covariates, such as one or more candidate genes or an environmental exposure. Recently, attention has turned to study designs that mandate the inclusion of family members in addition to a proband. Standard models for analysis assume independent observations, which is unlikely to be true for family data, and the usual standard errors for the regression parameter estimates may be too large or too small, depending on the distribution of the covariates within and between families. The consequences of familial correlation on the study efficiency can be measured by a design effect that is equivalent to the relative information in a sample of unrelated individuals compared to a sample of families with the same number of individuals. We examine design effects for studies in association, and illustrate how the design effect is influenced by the intra-familial distribution of covariate values such as would be expected for a candidate gene. Typical design effects for a candidate gene range between 1.1 and 2.4, depending on the size of the family and the amount of unexplained familial correlation. These values correspond to a modest 10% increase in the required sample size up to more than doubling the requirements. Design effect values are useful in study design to compare the efficiency of studies that sample families versus independent individuals and to determine sample size requirements that account for familial correlation.

Pedigree selection and tests of linkage in a Hutterite asthma pedigree.

We explored methods for kinship and linkage analysis in a Hutterite pedigree comprising 1,544 individuals, 72 of whom were diagnosed with asthma. Subpedigrees were selected by (a) identifying nuclear families containing asthmatics, (b) identifying couples with many asthmatic descendants in an ad hoc manner, and (c) finding the most recent common ancestors of all asthmatics. Markov chain Monte Carlo (MCMC) methods were used to estimate allele sharing in the larger subpedigrees and transmission/disequilibrium tests were performed on nuclear families. On chromosome 5q near the cytokine cluster, modest evidence for linkage to asthma was obtained. Using MCMC, we were able to evaluate the evidence for linkage in complex subpedigrees of several hundred individuals, and hence, incorporate some of the co-ancestry of this founder population.

Racial differences in survival of patients on dialysis.

BACKGROUND: Recent studies have documented racial differences in the crude mortality rates of patients on dialysis. However, proper interpretation of these findings requires adjustment for potential confounders and comorbid risk factors between the racial groups. METHODS: We examined the clinical data on 3752 Caucasian patients, 451 Southeast Asian patients, 322 South Asian patients, and 319 black patients who were treated with hemodialysis or peritoneal dialysis under a Universal Health Care system in Toronto and prospectively followed between 1981 and 1995. In all patients, a number of comorbid risk factors for survival was assessed at the start of dialysis and was reassessed with their outcome status (that is, continued dialysis, transplantation, death, or loss to follow-up) at least every six months. Cox proportional hazards analysis was used to fit multivariate models predicting patient survival. Pairwise comparisons of the relative hazards of death between the racial groups were performed after stratifying for cardiovascular disease, diabetes mellitus, and hypertension at the start of dialysis, and were adjusted for differences in other comorbid risk factors. RESULTS: The risk of death in Caucasian patients was significantly increased when compared with Southeast Asian patients, South Asian patients, and black patients [multivariate relative hazards (95% CI): 1.63 (1.36 to 1.97), 1.36 (1.07 to 1.73), 1.34 (1.07 to 1.67), respectively]. Additionally, we detected an interaction between race and cigarette smoking (P < 0. 004), suggesting that in the dialysis patients who smoked, whites had a higher mortality risk compared with non-whites. CONCLUSIONS: Differences in patient survival on dialysis exist between racial groups. However, the genetic and environmental determinants that underlie these differences are presently unknown.

Linkage of tuberculosis to chromosome 2q35 loci, including NRAMP1, in a large aboriginal Canadian family.

An epidemic of tuberculosis occurred in a community of Aboriginal Canadians during the period 1987-89. Genetic and epidemiologic data were collected on an extended family from this community, and the evidence for linkage to NRAMP1, a candidate gene for susceptibility to mycobacterial diseases, was assessed. Individuals were grouped into risk (liability) classes based on vaccination, age, previous disease, and tuberculin skin-test results. Under the assumption of a dominant mode of inheritance and a relative risk of 10, which is associated with the high-risk genotypes, a maximum LOD score of 3.81 was observed for linkage between a tuberculosis-susceptibility locus and D2S424, which is located just distal to NRAMP1, in chromosome region 2q35. Significant linkage was also observed between a tuberculosis-susceptibility locus and a haplotype of 10 NRAMP1 intragenic variants. No linkage to the major histocompatibility-complex region on chromosome 6p was observed, despite distortion of transmission from one member of the oldest couple to their affected offspring. The ability to assign individuals to risk classes was crucial to the success of this study.

The impact of transmission-ratio distortion on allele sharing in affected sibling pairs.

There have been recent reports of transmission-ratio distortion (TRD) or segregation distortion in families not selected for genetic disease. If TRD exists but is ignored, linkage studies searching for disease genes in affected relatives may be misinterpreted. We show that the identical-by-descent sharing patterns for affected sib pairs are strongly affected by TRD and, further, that the estimated statistical significance of a sib-pair linkage study may be extremely biased. However, we also show that, if TRD is suspected during the planning stage of a study, the planned sample size of the study needs to be increased by only a small amount to maintain the desired power.

Validation of linkage by sampling based on environmental exposures.

After detecting linkage in one sample, most researchers will attempt to validate this finding in another sample. Three strategies for validating a primary linkage were compared, with a focus on methods that might be appropriate in the presence of gene x environment interaction. First, a validation sample was collected and analyzed using the same ascertainment procedure and methods as the primary sample. Second, a sample of families with particular exposure patterns were ascertained subsequent to a significant test for heterogeneity due to the exposure in the primary sample. A third strategy ascertained by exposure status when exposure-defined subgroup tests were significant in the primary sample. The second strategy reduced the number of false positive linkage signals identified through exposure subgroup identification (i.e., the third strategy), but in this GAW11 data that contained no qualitative gene x environment interactions, it had poor sensitivity. Power to detect heterogeneity depends on the differences in risk between exposed and unexposed.

Analysis of affected sib pairs, with covariates--with and without constraints.

Covariate models have previously been developed as an extension to affected-sib-pair methods in which the covariate effects are jointly estimated with the degree of excess allele sharing. These models can estimate the differences in sib-pair allele sharing that are associated with measurable environment or genes. When there are no covariates, the pattern of identical-by-descent allele sharing in affected sib pairs is expected to fall within a small triangular region of the potential parameter space, under most genetic models. By restriction of the estimated allele sharing to this triangle, improved power is obtained in tests for genetic linkage. When the affected-sib-pair model is generalized to allow for covariates that affect allele sharing, however, new constraints and new methods for the application of constraints are required. Three generalized constraint methods are proposed and evaluated by use of simulated data. The results compare the power of the different methods, with and without covariates, for a single-gene model with age-dependent onset and for quantitative and qualitative gene-environment and gene-gene interaction models. Covariates can improve the power to detect linkage and can be particularly valuable when there are qualitative gene-environment interactions. In most situations, the best strategy is to assume that there is no dominance variance and to obtain constrained estimates for covariate models under this assumption.

Statistical issues in a metaregression analysis of randomized trials: impact on the dietary sodium intake and blood pressure relationship.

In a meta-analysis of randomized trials of the effects of dietary sodium interventions on blood pressure, we found substantial heterogeneity among the studies. We were interested in evaluating whether measurement error, known to be a problem for dietary sodium measures, publication bias, or confounding factors could be responsible for the heterogeneity. A measurement error correction was developed that corrects both the slope and the intercept and takes into account the sample size of each study and the number of measurements taken on an individual. The measurement error correction had a minimal effect on the estimates, although it performed well in simulated data. A smoothed scatter plot was used to assess publication bias. Metaregressions provide a convenient way to jointly assess the effects of several factors, but care must be taken to fit an appropriate model.