The effect of linkage disequilibrium on linkage analysis of incomplete pedigrees

Dense SNP maps can be highly informative for linkage studies. But when parental genotypes are missing, multipoint linkage scores can be inflated in regions with substantial marker-marker linkage disequilibrium (LD). Such regions were observed in the Affymetrix SNP genotypes for the Genetic Analysis Workshop 14 (GAW14) Collaborative Study on the Genetics of Alcoholism (COGA) dataset, providing an opportunity to test a novel simulation strategy for studying this problem. First, an inheritance vector (with or without linkage present) is simulated for each replicate, i.e., locations of recombinations and transmission of parental chromosomes are determined for each meiosis. Then, two sets of founder haplotypes are superimposed onto the inheritance vector: one set that is inferred from the actual data and which contains the pattern of LD; and one set created by randomly selecting parental alleles based on the known allele frequencies, with no correlation (LD) between markers. Applying this strategy to a map of 176 SNPs (66 Mb of chromosome 7) for 100 replicates of 116 sibling pairs, significant inflation of multipoint linkage scores was observed in regions of high LD when parental genotypes were set to missing, with no linkage present. Similar inflation was observed in analyses of the COGA data for these affected sib pairs with parental genotypes set to missing, but not after reducing the marker map until r(2 )between any pair of markers was ≤ 0.05. Additional simulation studies of affected sib pairs assuming uniform LD throughout a marker map demonstrated inflation of significance levels at r(2 )values greater than 0.05. When genotypes are available only from two affected siblings in many families in a sample, trimming SNP maps to limit r(2 )to 0–0.05 for all marker pairs will prevent inflation of linkage scores without sacrificing substantial linkage information. Simulation studies on the observed pedigree structures and map can also be used to determine the effect of LD on a particular study

Statistical corrections of linkage data suggest predominantly cis regulations of gene expression

Morley et al. (Nature 2004, 430:743–747) detected significant linkages to the expression levels of 142 genes (of 3554) at a reported threshold of genome-wide p = 0.001 (LOD ≈ 5.3), using 14 three-generation Centre d'Etude du Polymorphisme Humain pedigrees. Most of the linkages (77%) were trans, i.e., more than 5 Mb from the expressed gene. However, the analysis did not account for the expected anti-conservative effect of the skewed distribution of score- or regression-based statistics in large sibships, or for the possible variance distortion due to correlations among tests. Therefore, we re-analyzed their data, using a robust score statistic for the entire pedigrees and correcting the p-values for skewness. We found that a LOD of 5.3 had a skewness-corrected genome-wide p-value of 0.016 instead of 0.001 (a result that we confirmed using simulation), with around 50 expected false positives. We then further corrected for correlation among the (skew-corrected) p-values by using Efron's method for obtaining the empirical null distribution. Setting a threshold of FDR = 10% (Z = 6.4, LOD = 8.9), we detected linkage for the expression levels of 22 genes, 19 of which are cis. Limiting the analysis to cis regions, linkage was detected to the expression levels of 46 genes with 4.6 expected false positives (FDR = 10%)

QuickSNP: an automated web server for selection of tagSNPs

Although large-scale genetic association studies involving hundreds to thousands of SNPs have become feasible, the associated cost is substantial. Even with the increased efficiency introduced by the use of tagSNPs, researchers are often seeking ways to maximize resource utilization given a set of SNP-based gene-mapping goals. We have developed a web server named QuickSNP in order to provide cost-effective selection of SNPs, and to fill in some of the gaps in existing SNP selection tools. One useful feature of QuickSNP is the option to select only gene-centric SNPs from a chromosomal region in an automated fashion. Other useful features include automated selection of coding non-synonymous SNPs, SNP filtering based on inter-SNP distances and information regarding the availability of genotyping assays for SNPs and whether they are present on whole genome chips. The program produces user-friendly summary tables and results, and a link to a UCSC Genome Browser track illustrating the position of the selected tagSNPs in relation to genes and other genomic features. We hope the unique combination of features of this server will be useful for researchers aiming to select markers for their genotyping studies. The server is freely available and can be accessed at the URL http://bioinformoodics.jhmi.edu/quickSNP.pl

Location Preferences of Family Firms: Strategic Decision Making or “Home Sweet Home”?

Selecting a business location is among the most important strategic decisions for family firms. Yet the separate demands of the family and the business often prove difficult to balance. A comparison of location preferences in family and nonfamily firms provides insight into the family influence on strategic decision making.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/67069/2/10_1111_j_1741-6248_1992_00271_x.pd

Runs of Homozygosity Implicate Autozygosity as a Schizophrenia Risk Factor

Autozygosity occurs when two chromosomal segments that are identical from a common ancestor are inherited from each parent. This occurs at high rates in the offspring of mates who are closely related (inbreeding), but also occurs at lower levels among the offspring of distantly related mates. Here, we use runs of homozygosity in genome-wide SNP data to estimate the proportion of the autosome that exists in autozygous tracts in 9,388 cases with schizophrenia and 12,456 controls. We estimate that the odds of schizophrenia increase by ∼17% for every 1% increase in genome-wide autozygosity. This association is not due to one or a few regions, but results from many autozygous segments spread throughout the genome, and is consistent with a role for multiple recessive or partially recessive alleles in the etiology of schizophrenia. Such a bias towards recessivity suggests that alleles that increase the risk of schizophrenia have been selected against over evolutionary time

ACSL6 Is Associated with the Number of Cigarettes Smoked and Its Expression Is Altered by Chronic Nicotine Exposure

Individuals with schizophrenia tend to be heavy smokers and are at high risk for tobacco dependence. However, the nature of the comorbidity is not entirely clear. We previously reported evidence for association of schizophrenia with SNPs and SNP haplotypes in a region of chromosome 5q containing the SPEC2, PDZ-GEF2 and ACSL6 genes. In this current study, analysis of the control subjects of the Molecular Genetics of Schizophrenia (MGS) sample showed similar pattern of association with number of cigarettes smoked per day (numCIG) for the same region. To further test if this locus is associated with tobacco smoking as measured by numCIG and FTND, we conducted replication and meta-analysis in 12 independent samples (n\u3e16,000) for two markers in ACSL6 reported in our previous schizophrenia study. In the meta-analysis of the replication samples, we found that rs667437 and rs477084 were significantly associated with numCIG (p = 0.00038 and 0.00136 respectively) but not with FTND scores. We then used in vitro and in vivo techniques to test if nicotine exposure influences the expression of ACSL6 in brain. Primary cortical culture studies showed that chronic (5-day) exposure to nicotine stimulated ACSL6 mRNA expression. Fourteen days of nicotine administration via osmotic mini pump also increased ACSL6 protein levels in the prefrontal cortex and hippocampus of mice. These increases were suppressed by injection of the nicotinic receptor antagonist mecamylamine, suggesting that elevated expression ofACSL6 requires nicotinic receptor activation. These findings suggest that variations in theACSL6 gene may contribute to the quantity of cigarettes smoked. The independent associations of this locus with schizophrenia and with numCIG in non-schizophrenic subjects suggest that this locus may be a common liability to both conditions

HIV-Neutralizing Activity of Cationic Polypeptides in Cervicovaginal Secretions of Women in HIV-Serodiscordant Relationships

HIV exposed seronegative (HESN) women represent the population most in need of a prophylactic antiviral strategy. Mucosal cationic polypeptides can potentially be regulated for this purpose and we here aimed to determine their endogenous expression and HIV neutralizing activity in genital secretions of women at risk of HIV infection.Cervicovaginal secretions (CVS) of Kenyan women in HIV-serodiscordant relationships (HESN, n = 164; HIV seropositive, n = 60) and low-risk controls (n = 72) were assessed for the cationic polypeptides HNP1–3, LL-37 and SLPI by ELISA and for HIV neutralizing activity by a PBMC-based assay using an HIV primary isolate. Median levels of HNP1–3 and LL-37 in CVS were similar across study groups. Neither HSV-2 serostatus, nor presence of bacterial vaginosis, correlated with levels of HNP1–3 or LL-37 in the HESN women. However, an association with their partner's viral load was observed. High viral load (>10,000 HIV RNA copies/ml plasma) correlated with higher levels of HNP1–3 and LL-37 (p = 0.04 and 0.03, respectively). SLPI was most abundant in the low-risk group and did not correlate with male partner's viral load in the HESN women. HIV neutralizing activity was found in CVS of all study groups. In experimental studies, selective depletion of cationic polypeptides from CVS rendered the remaining CVS fraction non-neutralizing, whereas the cationic polypeptide fraction retained the activity. Furthermore, recombinant HNP1–3 and LL-37 could induce neutralizing activity when added to CVS lacking intrinsic activity.These findings show that CVS from HESN, low-risk, and HIV seropositive women contain HIV neutralizing activity. Although several innate immune proteins, including HNP1–3 and LL-37, contribute to this activity these molecules can also have inflammatory properties. This balance is influenced by hormonal and environmental factors and in the present HIV serodiscordant couple cohort study we show that a partner's viral load is associated with levels of such molecules

A Rare Functional Noncoding Variant at the GWAS-Implicated MIR137/MIR2682 Locus Might Confer Risk to Schizophrenia and Bipolar Disorder

Schizophrenia (SZ) genome-wide association studies (GWASs) have identified common risk variants in >100 susceptibility loci; however, the contribution of rare variants at these loci remains largely unexplored. One of the strongly associated loci spans MIR137 (miR137) and MIR2682 (miR2682), two microRNA genes important for neuronal function. We sequenced ∼6.9 kb MIR137/MIR2682 and upstream regulatory sequences in 2,610 SZ cases and 2,611 controls of European ancestry. We identified 133 rare variants with minor allele frequency (MAF) <0.5%. The rare variant burden in promoters and enhancers, but not insulators, was associated with SZ (p = 0.021 for MAF < 0.5%, p = 0.003 for MAF < 0.1%). A rare enhancer SNP, 1:g.98515539A>T, presented exclusively in 11 SZ cases (nominal p = 4.8 × 10−4). We further identified its risk allele T in 2 of 2,434 additional SZ cases, 11 of 4,339 bipolar (BP) cases, and 3 of 3,572 SZ/BP study controls and 1,688 population controls; yielding combined p values of 0.0007, 0.0013, and 0.0001 for SZ, BP, and SZ/BP, respectively. The risk allele T of 1:g.98515539A>T reduced enhancer activity of its flanking sequence by >50% in human neuroblastoma cells, predicting lower expression of MIR137/MIR2682. Both empirical and computational analyses showed weaker transcription factor (YY1) binding by the risk allele. Chromatin conformation capture (3C) assay further indicated that 1:g.98515539A>T influenced MIR137/MIR2682, but not the nearby DPYD or LOC729987. Our results suggest that rare noncoding risk variants are associated with SZ and BP at MIR137/MIR2682 locus, with risk alleles decreasing MIR137/MIR2682 expression

Additive genetic variation in schizophrenia risk is shared by populations of African and European descent

Previous studies have emphasized ethnically heterogeneous human leukocyte antigen (HLA) classical allele associations to rheumatoid arthritis (RA) risk. We fine-mapped RA risk alleles within the major histocompatibility complex (MHC) in 2782 seropositive RA cases and 4315 controls of Asian descent. We applied imputation to determine genotypes for eight class I and II HLA genes to Asian populations for the first time using a newly constructed pan-Asian reference panel. First, we empirically measured high imputation accuracy in Asian samples. Then we observed the most significant association in HLA-DRβ1 at amino acid position 13, located outside the classical shared epitope (Pomnibus = 6.9 × 10(-135)). The individual residues at position 13 have relative effects that are consistent with published effects in European populations (His > Phe > Arg > Tyr ≅ Gly > Ser)--but the observed effects in Asians are generally smaller. Applying stepwise conditional analysis, we identified additional independent associations at positions 57 (conditional Pomnibus = 2.2 × 10(-33)) and 74 (conditional Pomnibus = 1.1 × 10(-8)). Outside of HLA-DRβ1, we observed independent effects for amino acid polymorphisms within HLA-B (Asp9, conditional P = 3.8 × 10(-6)) and HLA-DPβ1 (Phe9, conditional P = 3.0 × 10(-5)) concordant with European populations. Our trans-ethnic HLA fine-mapping study reveals that (i) a common set of amino acid residues confer shared effects in European and Asian populations and (ii) these same effects can explain ethnically heterogeneous classical allelic associations (e.g. HLA-DRB1*09:01) due to allele frequency differences between populations. Our study illustrates the value of high-resolution imputation for fine-mapping causal variants in the MHC

Reproductive Cycle-associated Mood Symptoms In Women With Major Depression And Bipolar Disorder

Background: We sought to determine the prevalence of, and association between, reproductive cycle-associated mood symptoms in women with affective disorders. We hypothesized that symptoms would correlate with each other across a woman's reproductive life span in both major depression (MDD) and bipolar I disorder (BP). Methods: 2412 women with, MDD or BP were asked standardized questions about mood symptoms prior to menstruation, within a month of childbirth and during perimenopause. Lifetime rates for each of these symptom types were determined and an odds ratio was calculated correlating each of the types with the others. Results: Of 2524 women with mood disorders, 67.7% reported premenstrual symptoms. Of those at risk, 20.9% reported postpartum symptoms and 26.4% reported perimenopausal symptoms. The rates did not differ between women with MDD and BP but were significantly different from women who were never ill. The symptoms were significantly correlated in women with MDD with odds ratios from 1.66 to 1.82, but were not in women with BP. Limitations: This is a secondary analysis of a sample that was collected for other purposes and is based upon retrospecitve reporting. Conclusions: Reproductive cycle-associated mood symptoms were commonly reported in women with mood disorders and did not differ based on diagnosis. In MDD, but not BP, the occurrence of these symptoms was trait-like as the presence of one predicted the occurrence of the others. Further prospective study is required to clarify the determinants of this trait