Search

Background Human genetic susceptibility for tuberculosis (TB) has been demonstrated by several studies, but few have examined multiple innate and adaptive immunity genes comprehensively, age-specific effects, and/or resistance to Mycobacterium tuberculosis (Mtb) infection (RSTR). We hypothesized that RSTR, defined by a persistently negative tuberculin skin test, may have different genetic influences than Mtb disease. Methods We examined 29 candidate genes in pathways that mediate immune responses to Mtb in subjects in a household contact study in Kampala, Uganda. We genotyped 546 haplotype-tagging single nucleotide polymorphisms (SNPs) in 835 individuals from 481 families; 28.7% had TB, 10.5% were RSTR, and the remaining 60.8% had latent Mtb infection. Results Among our most significant findings were SNPs in TICAM2 (p=3.6×10−6) and IL1B (p=4.3×10−5) associated with TB. Multiple SNPs in IL4 and TOLLIP were associated with TB (p

OBJECTIVES Speech sound disorder (SSD) is one of the most common communication disorders, with a prevalence rate of 16% at 3 years of age, and an estimated 3.8% of children still presenting speech difficulties at 6 years of age. Several studies have identified promising associations between communication disorders and genes in brain and neuronal pathways; however, there have been few studies focusing on SSD and its associated endophenotypes. On the basis of the hypothesis that neuronal genes may influence endophenotypes common to communication disorders, we focused on three genes related to brain and central nervous system functioning: the dopamine D2 receptor (DRD2) gene, the arginine-vasopressin receptor 1a (AVPR1A) gene, and the microcephaly-associated protein gene (ASPM). METHODS We examined the association of these genes with key endophenotypes of SSD - phonological memory measured through multisyllabic and nonword repetition, vocabulary measured using the Expressive One Word Picture Vocabulary Test and Peabody Picture Vocabulary Test, and reading decoding measured using the Woodcock Reading Mastery Tests Revised - as well as with the clinical phenotype of SSD. We genotyped tag single nucleotide polymorphisms in these genes and examined 498 individuals from 180 families. RESULTS These data show that several single nucleotide polymorphisms in all three genes were associated with phonological memory, vocabulary, and reading decoding, with P less than 0.05. Notably, associations in AVPR1A (rs11832266) were significant after multiple testing correction. Gene-level tests showed that DRD2 was associated with vocabulary, ASPM with vocabulary and reading decoding, and AVPR1A with all three endophenotypes. CONCLUSION Endophenotypes common to SSD, language impairment, and reading disability are all associated with these neuronal pathway genes.

Neurodegenerative diseases affecting the macula constitute a major cause of incurable vision loss and exhibit considerable clinical and genetic heterogeneity, from early-onset monogenic disease to multifactorial late-onset age-related macular degeneration (AMD). As part of our continued efforts to define genetic causes of macular degeneration, we performed whole exome sequencing in four individuals of a two-generation family with autosomal dominant maculopathy and identified a rare variant p.Glu1144Lys in Fibrillin 2 (FBN2), a glycoprotein of the elastin-rich extracellular matrix (ECM). Sanger sequencing validated the segregation of this variant in the complete pedigree, including two additional affected and one unaffected individual. Sequencing of 192 maculopathy patients revealed additional rare variants, predicted to disrupt FBN2 function. We then undertook additional studies to explore the relationship of FBN2 to macular disease. We show that FBN2 localizes to Bruch's membrane and its expression appears to be reduced in aging and AMD eyes, prompting us to examine its relationship with AMD. We detect suggestive association of a common FBN2 non-synonymous variant, rs154001 (p.Val965Ile) with AMD in 10 337 cases and 11 174 controls (OR = 1.10; P-value = 3.79 × 10(-5)). Thus, it appears that rare and common variants in a single gene--FBN2--can contribute to Mendelian and complex forms of macular degeneration. Our studies provide genetic evidence for a key role of elastin microfibers and Bruch's membrane in maintaining blood-retina homeostasis and establish the importance of studying orphan diseases for understanding more common clinical phenotypes.

One in three people has been infected with Mycobacterium tuberculosis (MTB), and the risk for MTB infection in HIV-infected individuals is even higher. We hypothesized that HIV-positive individuals living in tuberculosis-endemic regions who do not get infected by Mycobacterium tuberculosis are genetically resistant. Using an “experiment of nature” design that proved successful in our previous work, we performed a genome-wide association study of tuberculin skin test positivity using 469 HIV-positive patients from prospective study cohorts of tuberculosis from Tanzania and Uganda to identify genetic loci associated with MTB infection in the context of HIV-infection. Among these individuals, 244 tested were tuberculin skin test (TST) positive either at enrollment or during the >8 year follow up, while 225 were not. We identified a genome-wide significant association between a dominant model of rs877356 and binary TST status in the combined cohort (Odds ratio = 0.2671, p = 1.22x10-8). Association was replicated with similar significance when examining TST induration as a continuous trait. The variant lies in the 5q31.1 region, 57kb downstream from IL9. Two-locus analyses of association of variants near rs877356 showed a haplotype comprised of rs877356 and an IL9 missense variant, rs2069885, had the most significant association (p = 1.59x10-12). We also replicated previously linked loci on chromosomes 2, 5, and 11. IL9 is a cytokine produced by mast cells and TH2 cells during inflammatory responses, providing a possible link between airway inflammation and protection from MTB infection. Our results indicate that studying uninfected, HIV-positive participants with extensive exposure increases the power to detect associations in complex infectious disease., Author summary Approximately one-third of the world’s population has been exposed to Mycobacterium tuberculosis, the bacterium that causes tuberculosis. A small number of those infected develop active disease; however, there is a substantial portion of exposed people who do not even show evidence of an immunological response. These people who appear to resist infection, as measured by a negative tuberculin skin test, represent a subpopulation from which we can learn about resistance. We used a genome-wide approach to study the genetic basis of this resistance in unique cohorts of hypervulnerable, HIV-positive individuals from Uganda and Tanzania, in which exposure was virtually assured. We identified one locus that was highly significantly associated and conferred more than 70% protection from infection. The most significant variant, rs8773656, was near IL9 and SLC25A48, and a haplotype including this variant and a missense mutation in IL9 was even more significantly associated with negative skin tests. Although it is impossible based solely on our data to determine the causal variant or genes, IL9 is an attractive candidate as its product has previously been associated with bronchial hyperresponsiveness, thereby providing a possible link between inflammation and protection from Mycobacterium tuberculosis infection.

Type 2 diabetes (T2D) affects Latinos at twice the rate seen in populations of European descent. We recently identified a risk haplotype spanning SLC16A11 that explains ∼20% of the increased T2D prevalence in Mexico. Here, through genetic fine-mapping, we define a set of tightly linked variants likely to contain the causal allele(s). We show that variants on the T2D-associated haplotype have two distinct effects: (1) decreasing SLC16A11 expression in liver and (2) disrupting a key interaction with basigin, thereby reducing cell-surface localization. Both independent mechanisms reduce SLC16A11 function and suggest SLC16A11 is the causal gene at this locus. To gain insight into how SLC16A11 disruption impacts T2D risk, we demonstrate that SLC16A11 is a proton-coupled monocarboxylate transporter and that genetic perturbation of SLC16A11 induces changes in fatty acid and lipid metabolism that are associated with increased T2D risk. Our findings suggest that increasing SLC16A11 function could be therapeutically beneficial for T2D. Video Abstract [Figure presented] Keywords: type 2 diabetes (T2D); genetics; disease mechanism; SLC16A11; MCT11; solute carrier (SLC); monocarboxylates; fatty acid metabolism; lipid metabolism; precision medicine

Genome-wide variation data with millions of genetic markers have become commonplace. However, the potential for interpretation and application of these data for clinical assessment of outcomes of interest, and prediction of disease risk, is currently not fully realized. Many common complex diseases now have numerous, well-established risk loci and likely harbor many genetic determinants with effects too small to be detected at genome-wide levels of statistical significance. A simple and intuitive approach for converting genetic data to a predictive measure of disease susceptibility is to aggregate the effects of these loci into a single measure, the genetic risk score. Here, we describe some common methods and software packages for calculating genetic risk scores and polygenic risk scores, with focus on studies of common complex diseases. We review the basic information needed, as well as important considerations for constructing genetic risk scores, including specific requirements for phenotypic and genetic data, and limitations in their application. © 2019 by John Wiley & Sons, Inc.

IMPORTANCE: Mutations in the myocilin (MYOC) gene are the most common molecularly defined cause of primary open-angle glaucoma that typically occurs in patients with high intraocular pressures (IOP). One MYOC mutation, p.Gln368Ter, has been associated with as many as 1.6% of primary open-angle glaucoma cases that had a mean maximum recorded IOP of 30 mm Hg. However, to our knowledge, the role of the p.Gln368Ter mutation in patients with normal-tension glaucoma (NTG) with an IOP of 21 mm Hg or lower has not been investigated. OBJECTIVE: To evaluate the role of the p.Gln368Ter MYOC mutation in patients with NTG. DESIGN, SETTING, AND PARTICIPANTS: In this case-control study of the prevalence of the p.Gln368Ter mutation in patients with NTG, cohort 1 was composed of 772 patients with NTG and 2152 controls from the United States (Iowa, Minnesota, and New York) and England and cohort 2 was composed of 561 patients with NTG and 2606 controls from the Massachusetts Eye and Ear Infirmary and the NEIGHBORHOOD consortium. Genotyping was conducted using real-time polymerase chain reaction that was confirmed with Sanger sequencing, the imputation of genome-wide association study data, or an analysis of whole-exome sequence data. Data analysis occurred between April 2007 and April 2018. MAIN OUTCOMES AND MEASURES: Comparison of the frequency of the p.Gln368Ter MYOC mutation between NTG cases and controls with the Fisher exact test. RESULTS: Of 6091 total participants, 3346 (54.9%) were women and 5799 (95.2%) were white. We detected the p.Gln368Ter mutation in 7 of 772 patients with NTG (0.91%) and 7 of 2152 controls (0.33%) in cohort 1 (P = .03). In cohort 2, we detected the p.Gln368Ter mutation in 4 of 561 patients with NTG (0.71%) and 10 of 2606 controls (0.38%; P = .15). When the cohorts were analyzed as a group, the p.Gln368Ter mutation was associated with NTG (odds ratio, 2.3; 95% CI, 0.98-5.3; P = .04). CONCLUSIONS AND RELEVANCE: In cohorts 1 and 2, the p.Gln368Ter mutation in MYOC was found in patients with IOPs that were 21 mm Hg or lower (NTG), although at a frequency that is lower than previously detected in patients with higher IOP. These data suggest that the p.Gln368Ter mutation may be associated with glaucoma in patients with normal IOPs as well as in patients with IOPs that are greater than 21 mm Hg.

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P

Vitamin D deficiency {defined by the blood concentration of 25-hydroxyvitamin D [25(OH)D]} has been associated with many adverse health outcomes. Genetic and nongenetic factors account for variation in 25(OH)D, but the role of interactions between these factors is unknown. To assess this, we examined 1204 women of European descent from the Carotenoids in Age-Related Eye Disease Study, an ancillary study of the Women’s Health Initiative Observational Study. Twenty-nine single nucleotide polymorphisms (SNPs) in 4 genes, GC, CYP2R1, DHCR7, and CYP24A1, from recent meta-analyses of 25(OH)D genome-wide association studies were genotyped. Associations between these SNPs and 25(OH)D were tested using generalized linear regression under an additive genetic model adjusted for age, blood draw month, and ancestry. Results were stratified by season of blood draw and, separately, vitamin D intake for the 6 SNPs showing a significant association with 25(OH)D at the P < 0.01 level. Two nonsynonymous SNPs in GC and 4 SNPs in CYP2R1 were strongly associated with 25(OH)D in individuals whose blood was drawn in summer (P ≤ 0.002) but not winter months and, independently, in individuals with vitamin D intakes ≥400 (P ≤ 0.004) but not

Author version made available in accordance with the publisher's policy., Refractive error is the most common eye disorder worldwide and is a prominent cause of blindness. Myopia affects over 30% of Western populations and up to 80% of Asians. The CREAM consortium conducted genome-wide meta-analyses, including 37,382 individuals from 27 studies of European ancestry and 8,376 from 5 Asian cohorts. We identified 16 new loci for refractive error in individuals of European ancestry, of which 8 were shared with Asians. Combined analysis identified 8 additional associated loci. The new loci include candidate genes with functions in neurotransmission (GRIA4), ion transport (KCNQ5), retinoic acid metabolism (RDH5), extracellular matrix remodeling (LAMA2 and BMP2) and eye development (SIX6 and PRSS56). We also confirmed previously reported associations with GJD2 and RASGRF1. Risk score analysis using associated SNPs showed a tenfold increased risk of myopia for individuals carrying the highest genetic load. Our results, based on a large meta-analysis across independent multiancestry studies, considerably advance understanding of the mechanisms involved in refractive error and myopia., Australian National Health and Medical Research Council.

The Illumina HumanExome BeadChip and other exome-based genotyping arrays offer inexpensive genotyping of some 240,000 mostly nonsynonymous coding variants across the human genome. The HumanExome chip, with its highly non-uniform distribution of markers and emphasis on rare coding variants, presents some unique challenges for quality control (QC) and data cleaning. Here, we describe QC procedures for HumanExome data, with examples of challenges specific to exome arrays from our experience cleaning a data set of ∼7,500 samples from the NEIGHBORHOOD Consortium. We focus on standard procedures for QC of genome-wide array data including genotype calling, sex verification, sample identity verification, relationship checking, and population structure that are complicated by the HumanExome panel's enrichment in rare, exonic variation. © 2016 by John Wiley & Sons, Inc.

International audience; Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. To identify new susceptibility loci, we performed meta-analysis on genome-wide association study (GWAS) results from eight independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significantly associated SNPs in two Australian studies (1,252 cases and 2,592 controls), three European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of the top SNPs identified three new associated loci: rs35934224[T] in TXNRD2 (odds ratio (OR) = 0.78, P = 4.05 × 10(-11)) encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] in ATXN2 (OR = 1.17, P = 8.73 × 10(-10)); and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76 × 10(-10)). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest new targets for preventative therapies.

Immunosuppression resulting from HIV infection increases the risk of progression to active tuberculosis (TB) both in individuals newly exposed to Mycobacterium tuberculosis (MTB) and in those with latent infections. We hypothesized that HIV-positive individuals who do not develop TB, despite living in areas where it is hyperendemic, provide a model of natural resistance. We performed a genome-wide association study of TB resistance by using 581 HIV-positive Ugandans and Tanzanians enrolled in prospective cohort studies of TB; 267 of these individuals developed active TB, and 314 did not. A common variant, rs4921437 at 5q33.3, was significantly associated with TB (odds ratio = 0.37, p = 2.11 × 10(-8)). This variant lies within a genomic region that includes IL12B and is embedded in an H3K27Ac histone mark. The locus also displays consistent patterns of linkage disequilibrium across African populations and has signals of strong selection in populations from equatorial Africa. Along with prior studies demonstrating that therapy with IL-12 (the cytokine encoded in part by IL12B, associated with longer survival following MTB infection in mice deficient in CD4 T cells), our results suggest that this pathway might be an excellent target for the development of new modalities for treating TB, especially for HIV-positive individuals. Our results also indicate that studying extreme disease resistance in the face of extensive exposure can increase the power to detect associations in complex infectious disease.

Author(s): Liu, Yutao; Bailey, Jessica Cooke; Helwa, Inas; Dismuke, W Michael; Cai, Jingwen; Drewry, Michelle; Brilliant, Murray H; Budenz, Donald L; Christen, William G; Chasman, Daniel I; Fingert, John H; Gaasterland, Douglas; Gaasterland, Terry; Gordon, Mae O; Igo, Robert P; Kang, Jae H; Kass, Michael A; Kraft, Peter; Lee, Richard K; Lichter, Paul; Moroi, Sayoko E; Realini, Anthony; Richards, Julia E; Ritch, Robert; Schuman, Joel S; Scott, William K; Singh, Kuldev; Sit, Arthur J; Song, Yeunjoo E; Vollrath, Douglas; Weinreb, Robert; Medeiros, Felipe; Wollstein, Gadi; Zack, Donald J; Zhang, Kang; Pericak-Vance, Margaret A; Gonzalez, Pedro; Stamer, W Daniel; Kuchtey, John; Kuchtey, Rachel W; Allingham, R Rand; Hauser, Michael A; Pasquale, Louis R; Haines, Jonathan L; Wiggs, Janey L | Abstract: PurposeNoncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG).MethodsUsing the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR-182 expression in AH between five HTG cases and five controls.ResultsOnly rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] = 1.23, 95% confidence interval [CI]: 1.11-1.42, P = 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR = 1.26, 95% CI: 1.08-1.47, P = 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P = 0.03) without controlling for medication treatment.ConclusionsOur integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression.

Background The presence of population structure in a sample may confound the search for important genetic loci associated with disease. Our four samples in the Family Investigation of Nephropathy and Diabetes (FIND), European Americans, Mexican Americans, African Americans, and American Indians are part of a genome- wide association study in which population structure might be particularly important. We therefore decided to study in detail one component of this, individual genetic ancestry (IGA). From SNPs present on the Affymetrix 6.0 Human SNP array, we identified 3 sets of ancestry informative markers (AIMs), each maximized for the information in one the three contrasts among ancestral populations: Europeans (HAPMAP, CEU), Africans (HAPMAP, YRI and LWK), and Native Americans (full heritage Pima Indians). We estimate IGA and present an algorithm for their standard errors, compare IGA to principal components, emphasize the importance of balancing information in the ancestry informative markers (AIMs), and test the association of IGA with diabetic nephropathy in the combined sample. Results A fixed parental allele maximum likelihood algorithm was applied to the FIND to estimate IGA in four samples: 869 American Indians; 1385 African Americans; 1451 Mexican Americans; and 826 European Americans. When the information in the AIMs is unbalanced, the estimates are incorrect with large error. Individual genetic admixture is highly correlated with principle components for capturing population structure. It takes ~700 SNPs to reduce the average standard error of individual admixture below 0.01. When the samples are combined, the resulting population structure creates associations between IGA and diabetic nephropathy. Conclusions The identified set of AIMs, which include American Indian parental allele frequencies, may be particularly useful for estimating genetic admixture in populations from the Americas. Failure to balance information in maximum likelihood, poly-ancestry models creates biased estimates of individual admixture with large error. This also occurs when estimating IGA using the Bayesian clustering method as implemented in the program STRUCTURE. Odds ratios for the associations of IGA with disease are consistent with what is known about the incidence and prevalence of diabetic nephropathy in these populations. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2654-x) contains supplementary material, which is available to authorized users.

Age-related macular degeneration (AMD) is the leading cause of blindness in the elderly in the developed world. We conducted a genome-wide association study in a series of families enriched for AMD and completed a meta-analysis of this new data with results from reanalysis of an existing study of a late-stage case-control cohort. We tested the top findings for replication in 1896 cases and 1866 controls and identified two novel genetic protective factors for AMD. In addition to the complement factor H (CFH) (P=2.3 × 10⁻⁶⁴) and age-related maculopathy susceptibility 2 (ARMS2) (P=1.2 × 10⁻⁶⁰) loci, we observed a protective effect at rs429608, an intronic SNP in SKIV2L (P=5.3 × 10⁻¹⁵), a gene near the complement component 2 (C2)/complement factor B (BF) locus, that indicates the protective effect may be mediated by variants other than the C2/BF variants previously studied. Haplotype analysis at this locus identified three protective haplotypes defined by the rs429608 protective allele. We also identified a new potentially protective effect at rs2679798 in MYRIP (P=2.9 × 10⁻⁴), a gene involved in retinal pigment epithelium melanosome trafficking. Interestingly, MYRIP was initially identified in the family-based scan and was confirmed in the case-control set. From these efforts, we report the identification of two novel protective factors for AMD and confirm the previously known associations at CFH, ARMS2 and C3.

PURPOSE: To identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders. METHODS: One hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN. RESULTS: Genome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; β = −0.57, 95% confidence interval CI: −0.78 to −0.36; P = 1.7 × 10−7). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (β = −3.94, 95% CI: −5.23 to −2.66; P = 1.7 × 10−9). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT. CONCLUSIONS: The results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.

Background Previous studies have shown that in addition to environmental influences, type 2 diabetes mellitus (T2DM) has a strong genetic component. The goal of the current study is to identify regions of linkage for T2DM in ethnically diverse populations. Methods Phenotypic and genotypic data were obtained from African American (AA; total number of individuals [N] = 1004), American Indian (AI; N = 883), European American (EA; N = 537), and Mexican American (MA; N = 1634) individuals from the Family Investigation of Nephropathy and Diabetes. Non-parametric linkage analysis, using an average of 4404 SNPs, was performed in relative pairs affected with T2DM in each ethnic group. In addition, family-based tests were performed to detect association with T2DM. Results Statistically significant evidence for linkage was observed on chromosome 4q21.1 (LOD = 3.13; genome-wide p = 0.04) in AA. In addition, a total of 11 regions showed suggestive evidence for linkage (estimated at LOD > 1.71), with the highest LOD scores on chromosomes 12q21.31 (LOD = 2.02) and 22q12.3 (LOD = 2.38) in AA, 2p11.1 (LOD = 2.23) in AI, 6p12.3 (LOD = 2.77) in EA, and 13q21.1 (LOD = . 2.24) in MA. While no region overlapped across all ethnic groups, at least five loci showing LOD > 1.71 have been identified in previously published studies. Conclusions The results from this study provide evidence for the presence of genes affecting T2DM on chromosomes 4q, 12q, and 22q in AA; 6p in EA; 2p in AI; and 13q in MA. The strong evidence for linkage on chromosome 4q in AA provides important information given the paucity of diabetes genetic studies in this population. Copyright © 2009 John Wiley & Sons, Ltd.

Haplotypes of closely linked single-nucleotide polymorphisms (SNPs) potentially offer greater power than individual SNPs to detect association between genetic variants and disease. We present a novel approach for association mapping in which density-based clustering of haplotypes reduces the dimensionality of the general linear model (GLM)-based score test of association implemented in the HaploStats software (Schaid et al. [2002] Am. J. Hum. Genet. 70:425-434). A flexible haplotype similarity score, a generalization of previously used measures, forms the basis, for grouping haplotypes of probable recent common ancestry. All haplotypes within a cluster are assigned the same regression coefficient within the GLM, and evidence for association is assessed with a score statistic. The approach is applicable to both binary and continuous trait data, and does not require prior phase information. Results of simulation studies demonstrated that clustering enhanced the power of the score test to detect association, under a variety of conditions, while preserving valid Type-I error. Improvement in performance was most dramatic in the presence of extreme haplotype diversity, while a slight improvement was observed even at low diversity. Our method also offers, for binary traits, a slight advantage in power over a similar approach based on an evolutionary model (Tzeng et al. [2006] Am. J. Hum. Genet. 78:231-242).

To understand the genetic architecture of dyslexia and identify the locations of genes involved, we performed linkage analyses in multigenerational families using a phonological memory phenotype--Nonword Repetition (NWR). A genome scan was first performed on 438 people from 51 families (DS-1) and linkage was assessed using variance components (VC), Bayesian oligogenic (BO), and parametric analyses. For replication, the genome scan and analyses were repeated on 693 people from 93 families (DS-2). For the combined set (DS-C), analyses were performed with all three methods in the regions that were identified in both samples. In DS-1, regions on chromosomes 4p, 6q, 12p, 17q, and 22q exceeded our initial threshold for linkage, with 17q providing a parametric LOD score of 3.2. Analysis with DS-2 confirmed the locations on chromosomes 4p and 12p. The strongest VC and BO signals in both samples were on chromosome 4p in DS-C, with a parametric multipoint LOD(max) of 2.36 for the 4p locus. Our linkage analyses of NWR in dyslexia provide suggestive and reproducible evidence for linkage to 4p12 and 12p in both samples, and significant evidence for linkage to 17q in one of the samples. These results warrant further studies of phonological memory and chromosomal regions identified here in other datasets.

Variance-components (VC) linkage analysis is a powerful model-free method for assessing linkage, but the distribution of VC logarithm of the odds ratio (LOD) scores may deviate substantially from the assumed asymptotic distribution. Typically, the null distribution of the VC-LOD score and other linkage statistics has been estimated by generating new genotype data independently of the trait data, and computing a linkage statistic for many such marker-simulated data sets. However, marker simulation is susceptible to errors in the assumed marker and map model and is computationally intensive. Here, we describe a method for generating posterior distributions of linkage statistics through simulation of trait data based on the original sample and on results from an initial scan using a Bayesian Markov-chain Monte Carlo (MCMC) approach for oligogenic segregation analysis. We use samples of oligogenic trait models taken from the posterior distribution to generate new samples of trait data, which were paired with the original marker data for analysis. Empirical P-values obtained from trait and marker simulation were similar when derived for several strong linkage signals from published linkage scans, and for analysis of data with a known, simulated, trait model. Furthermore, trait simulation produces the expected null distribution of VC-LOD scores and is computationally fast when marker identity-by-descent estimates from the original data could be reused. These results suggest that trait simulation gives valid estimates of statistical significance of linkage signals. Finally, these results also demonstrate the feasibility of obtaining empirical significance levels for evaluating Bayesian oligogenic linkage signals with either marker or trait simulation. Genet. Epidemiol. 2008. © 2007 Wiley–Liss, Inc.

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD., Author Summary Type 2 diabetes is the most common cause of severe kidney disease worldwide and diabetic kidney disease (DKD) associates with premature death. Individuals of non-European ancestry have the highest burden of type 2 DKD; hence understanding the causes of DKD remains critical to reducing health disparities. Family studies demonstrate that genes regulate the onset and progression of DKD; however, identifying these genes has proven to be challenging. The Family Investigation of Diabetes and Nephropathy consortium (FIND) recruited a large multi-ethnic collection of individuals with type 2 diabetes with and without kidney disease in order to detect genes associated with DKD. FIND discovered and replicated a DKD-associated genetic locus on human chromosome 6q25.2 (rs955333) between the SCAF8 and CNKSR genes. Findings were supported by significantly different expression of genes in this region from kidney tissue of subjects with, versus without DKD. The present findings identify a novel kidney disease susceptibility locus in individuals with type 2 diabetes which is consistent across subjects of differing ancestries. In addition, FIND results provide a rich catalogue of genetic variation in DKD patients for future research on the genetic architecture regulating this common and devastating disease.

Dyslexia is a common learning disability exhibited as a delay in acquiring reading skills despite adequate intelligence and instruction. Reading single real words (real-word reading, RWR) is especially impaired in many dyslexics. We performed a genome scan, using variance components (VC) linkage analysis and Bayesian Markov chain Monte Carlo (MCMC) joint segregation and linkage analysis, for three quantitative measures of RWR in 108 multigenerational families, with follow up of the strongest signals with parametric LOD score analyses. We used single-word reading efficiency (SWE) to assess speed and accuracy of RWR, and word identification (WID) to assess accuracy alone. Adjusting SWE for WID provided a third measure of RWR efficiency. All three methods of analysis identified a strong linkage signal for SWE on chromosome 13q. Based on multipoint analysis with 13 markers we obtained a MCMC intensity ratio (IR) of 53.2 (chromosome-wide P

Background/Aims: Complex traits pose a particular challenge to standard methods for segregation analysis (SA), and for such traits it is difficult to assess the ability of complex SA (CSA) to approximate the true mode of inheritance. Here we use an oligogenic Bayesian Markov chain Monte Carlo method for SA (OSA) to verify results from a single-locus likelihood-based CSA for data on a quantitative measure of reading ability. Methods: We compared the profile likelihood from CSA, maximized over the trait allele frequency, to the posterior distribution of genotype effects from OSA to explore differences in the overall parameter estimates from SA on the original phenotype data and the same data Winsorized to reduce the potential influence of three outlying data points. Results: Bayesian OSA revealed two modes of inheritance, one of which coincided with the QTL model from CSA. Winsorizing abolished the model originally estimated by CSA; both CSA and OSA identified only the second OSA model. Conclusion: Differences between the results from the two methods alerted us to the presence of influential data points, and identified the QTL model best supported by the data. Thus, the Bayesian OSA proved a valuable tool for assessing and verifying inheritance models from CSA.

Dyslexia is a common and complex developmental disorder manifested by unexpected difficulty in learning to read. Multiple different measures are used for diagnosis, and may reflect different biological pathways related to the disorder. Impaired phonological decoding (translation of written words without meaning cues into spoken words) is thought to be a core deficit. We present a genome scan of two continuous measures of phonological decoding ability: phonemic decoding efficiency (PDE) and word attack (WA). PDE measures both accuracy and speed of phonological decoding, whereas WA measures accuracy alone. Multipoint variance component linkage analyses (VC) and Markov chain Monte-Carlo (MCMC) multipoint joint linkage and segregation analyses were performed on 108 families. A strong signal was observed on chromosome 2 for PDE using both VC (LOD=2.65) and MCMC methods (intensity ratio (IR)=32.1). The IR is an estimate of the ratio of the posterior to prior probability of linkage in MCMC analysis. The chromosome 2 signal was not seen for WA. More detailed mapping with additional markers provided statistically significant evidence for linkage of PDE to chromosome 2, with VC-LOD=3.0 and IR=59.6 at D2S1399. Parametric analyses of PDE, using a model obtained by complex segregation analysis, provided a multipoint maximum LOD=2.89. The consistency of results from three analytic approaches provides strong evidence for a locus on chromosome 2 that influences speed but not accuracy of phonological decoding.

Dyslexia is a common, complex disorder, which is thought to have a genetic component. There have been numerous reports of linkage to several regions of the genome for dyslexia and continuous dyslexia-related phenotypes. We attempted to confirm linkage of continuous measures of (1) accuracy and efficiency of phonological decoding; and (2) accuracy of single word reading (WID) to regions on chromosomes 2p, 6p, 15q, and 18p, using 111 families with a total of 898 members. We used both single-marker and multipoint variance components linkage analysis and Markov Chain Monte Carlo (MCMC) joint segregation and linkage analysis for initial inspection of these regions. Positive results were followed with traditional parametric lod score analysis using a model estimated by MCMC segregation analysis. No positive linkage signals were found on chromosomes 2p, 6p, or 18p. Evidence of linkage of WID to chromosome 15q was found with both methods of analysis. The maximum single-marker parametric lod score of 2.34 was obtained at a distance of 3 cM from D15S143. Multipoint analyses localized the putative susceptibility gene to the interval between markers GATA50C03 and D15S143, which falls between a region implicated in a recent genome screen for attention-deficit/hyperactivity disorder, and DYX1C1, a candidate gene for dyslexia. This apparent multiplicity of linkage signals in the region for developmental disorders may be the result of errors in map and/or model specification obscuring the pleiotropic effect of a single gene on different phenotypes, or it may reflect the presence of multiple genes. This article contains supplementary material, which may be viewed at the American Journal of Medical Genetics website at http://www.interscience.wiley.com/jpages/0148-7299:1/suppmat/index.html. © 2004 Wiley-Liss, Inc.

RNA editing produces mature trypanosome mitochondrial mRNAs by uridylate (U) insertion and deletion. In insertion editing, Us are added to the pre-mRNA by a 3' terminal uridylyl transferase (TUTase) activity. We report the identification of a TUTase activity that copurifies with in vitro editing and is catalyzed by the integral editosome protein TbMP57. TbMP57 catalyzes the addition of primarily a single U to single-stranded (ss) RNA and adds the number of Us specified by a guide RNA to insertion editing-like substrates. TbMP57 is distinct from a previously identified TUTase that adds many Us to ssRNA and which we find is neither a stable editosome component nor does it add Us to editing-like substrates. Recombinant TbMP57 specifically interacts with the editosome protein TbMP81, and this interaction enhances the TUTase activity. These results suggest that TbMP57 catalyzes U addition to pre-mRNA during editing.

RNA editing in kinetoplastids is a type of post-transcriptional processing that changes mitochondrial mRNA sequences by the addition or deletion of uridines. Multiple enzymatic activities, such as endoribonuclease and RNA ligase, are associated with this process and exist in a multienzyme complex. Endonuclease activities from Trypanosoma brucei mitochondrial extracts were fractionated by sequential ion exchange and gel filtration chromatography. The RNA editing specific endonuclease activity co-fractionated with in vitro editing while another endonuclease activity with a different substrate specificity, and the majority of mtRNase P activity fractionated away from the editing activity. The pH, salt, temperature, and Mg 2+ optima of all three endonucleases were determined. All three activities are sensitive to high temperature and protease digestion. In addition, treatment with micrococcal nuclease resulted in partial disruption of the editing complex and decreased pre-cleaved in vitro insertion editing activity, suggesting that both RNA(s) and protein(s) are necessary in the intact functional complex.

RNA editing inserts and deletes uridylates (U's) in kinetoplastid mitochondrial pre-mRNAs by a series of enzymatic steps. Small guide RNAs (gRNAs) specify the edited sequence. Editing, though sometimes extensive, is precise. The effects of mutating pre-mRNA and gRNA sequences in, around, and upstream of the editing site on the specificity and efficiency of in vitro insertion editing were examined. U's could be added opposite guiding pyrimidines, but guiding purines, particularly A's, were required for efficient ligation. A base pair between mRNA and gRNA immediately upstream of the editing site was not required for insertion editing, although it greatly enhanced its efficiency and accuracy. In addition, a gRNA/mRNA duplex upstream of the editing site enhanced insertion editing when it was close to the editing site, but prevented cleavage, and hence editing, when immediately adjacent to the editing site. Thus, several aspects of mRNA-gRNA interaction, as well as gRNA base pairing with added U's, optimize editing efficiency, although they are not required for insertion editing.

The underlying ethos of dbGaP is that access to these data by secondary data analysts facilitates advancement of science. NIH has required that genome-wide association study data be deposited in the Database of Genotypes and Phenotypes (dbGaP) since 2003. In 2013, a proposed updated policy extended this requirement to next-generation sequencing data. However, recent literature and anecdotal reports suggest lingering logistical and ethical concerns about subject identifiability, informed consent, publication embargo enforcement, and difficulty in accessing dbGaP data. We surveyed the International Genetic Epidemiology Society (IGES) membership about their experiences. One hundred and seventy five (175) individuals completed the survey, a response rate of 27%. Of respondents who received data from dbGaP (43%), only 32% perceived the application process as easy but most (75%) received data within five months. Remaining challenges include difficulty in identifying an institutional signing official and an overlong application process. Only 24% of respondents had contributed data to dbGaP. Of these, 31% reported local IRB restrictions on data release; an additional 15% had to reconsent study participants before depositing data. The majority of respondents (56%) disagreed that the publication embargo period was sufficient. In response, we recommend longer embargo periods and use of varied data-sharing models rather than a one-size-fits-all approach.

Purpose We investigated whether mitochondrial DNA (mtDNA) variants affect the susceptibility of Fuchs endothelial corneal dystrophy (FECD). Methods Ten mtDNA variants defining European haplogroups were genotyped in a discovery dataset consisting of 530 cases and 498 controls of European descent from the Duke FECD cohort. Association tests for mtDNA markers and haplogroups were performed using logistic regression models with adjustment of age and sex. Subset analyses included controlling for additional effects of either the TCF4 SNP rs613872 or cigarette smoking. Our replication dataset was derived from the genome-wide association study (GWAS) of the FECD Genetics Consortium, where genotypes for three of 10 mtDNA markers were available. Replication analyses were performed to compare non-Duke cases to all GWAS controls (GWAS1, N = 3200), and to non-Duke controls (GWAS2, N = 3043). Results The variant A10398G was significantly associated with FECD (odds ratio [OR] = 0.72; 95% confidence interval [CI] = [0.53, 0.98]; P = 0.034), and remains significant after adjusting for smoking status (min P = 0.012). This variant was replicated in GWAS1 (P = 0.019) and GWAS2 (P = 0.036). Haplogroup I was significantly associated with FECD (OR = 0.46; 95% CI = [0.22, 0.97]; P = 0.041) and remains significant after adjusting for the effect of smoking (min P = 0.008) or rs613872 (P = 0.034). Conclusions The 10398G allele and Haplogroup I appear to confer significant protective effects for FECD. The effect of A10398G and Haplogroup I to FECD is likely independent of the known TCF4 variant. More data are needed to decipher the interaction between smoking and mtDNA haplogroups.

In a North American, HIV-positive, highly active antiretroviral therapy (HAART)-treated, adherent cohort of self-identified white and black patients, we previously observed that chemokine (C-C motif) receptor 5 (CCR5) -2459GA genotype had a strong association with time to achieve virologic success (TVLS) in black but not in white patients.Using 128 genome-wide ancestry informative markers, we performed a quantitative assessment of ancestry in these patients (n = 310) to determine (1) whether CCR5 -2459GA genotype is still associated with TVLS of HAART when ancestry, not self-identified race, is considered and (2) whether this association is influenced by varying African ancestry.We found that the interaction between CCR5 -2459GA genotype and African ancestry (≤ 0.125 vs. ≥ 0.425 and0.71 vs. ≥ 0.71) was significantly associated with TVLS (GG compared with AA, P = 0.044 and 0.018, respectively). Furthermore, the association between CCR5 -2459GA genotype and TVLS was stronger in patients with African ancestry ≥ 0.71 than in patients with African ancestry ≥ 0.452, in both Kaplan-Meier (log-rank P = 0.039 and 0.057, respectively, for AA, GA, and GG) and Cox proportional hazards regression (relative hazard for GG compared with AA 2.59 [95% confidence interval: 1.27 to 5.22; P = 0.01] and 2.26 [95% confidence interval: 1.18 to 4.32; P = 0.01], respectively) analyses.We observed that the association between CCR5 -2459GA genotype and TVLS of HAART increased with stronger African ancestry. Understanding the genomic mechanisms by which African ancestry influences this association is critical and requires further studies.

RNA editing in Trypanosoma brucei produces mature mRNAs by posttranscriptional insertion and deletion of uridylates (Us) by a series of catalytic steps, which include endoribonucleolytic cleavage, 3' terminal addition or removal of Us, and RNA ligation. Preedited mRNA (pre-mRNA) and guide RNA (gRNA) that are mutated at or near the editing site (ES) were used to examine the effects on the specificity of in vitro editing. Sequences that are not predicted to form a gRNA/pre-mRNA base pair immediately 5' to the ES still supported accurate editing. Substitution of a non-U nucleotide at various positions within a stretch of Us that are normally removed from the ES resulted in deletion of only the Us that were 3' to the substituted nucleotide. Overall, ES selection by the endoribonuclease, the specificity of the 3' exoribonuclease for Us, and ligation appear to act in concert to ensure the production of accurately edited RNA.

RNA editing in trypanosomes occurs by a series of enzymatic steps that are catalyzed by a macromolecular complex. The TbMP52 protein is shown to be a component of this complex, to have RNA ligase activity, and to be one of two adenylatable proteins in the complex. Regulated repression ofTbMP52blocks editing, which shows that it is a functional component of the editing complex. This repression is lethal in bloodforms of the parasite, indicating that editing is essential in the mammalian stage of the life cycle. The editing complex, which is present in all kinetoplastid parasites, may thus be a chemotherapeutic target.

RNA editing in kinetoplastid mitochondria inserts and deletes uridylates at multiple sites in pre-mRNAs as directed by guide RNAs. This occurs by a series of steps that are catalyzed by endoribonuclease, 3'-terminal uridylyl transferase, 3'-exouridylylase, and RNA ligase activities. A multiprotein complex that contains these activities and catalyzes deletion editing in vitro was enriched from Trypanosoma brucei mitochondria by sequential ion-exchange and gel filtration chromatography, followed by glycerol gradient sedimentation. The complex size is approximately 1,600 kDa, and the purified fraction contains 20 major polypeptides. A monoclonal antibody that was generated against the enriched complex reacts with an approximately 49-kDa protein and specifically immunoprecipitates in vitro deletion RNA editing activity. The protein recognized by the antibody was identified by mass spectrometry, and the corresponding gene, designated TbMP52, was cloned. Recombinant TbMP52 reacts with the monoclonal antibody. Another novel protein, TbMP48, which is similar to TbMP52, and its gene were also identified in the enriched complex. These results suggest that TbMP52 and TbMP48 are components of the RNA editing complex.

Objective: Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR. Methods: Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula. Results: The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4×10-5) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5×10-4) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5×10-4) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome. Conclusion: The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN. © 2013 Thameem et al.

Glutathione synthesis, a vital cellular process, depends on L-cystine uptake by the amino acid transporter, System x-C. Here we show that a second transporter, System X-AG, is required for normal System x-C activity and glutathione maintenance by employing somatic cell mutants of CHO-K1. Uptake by System x-C in two X-AG-null mutants is significantly lower than that of CHO-K1, either under control conditions or after prolonged treatment with an electrophile. In addition, levels of glutathione in control and treated mutant cells are less than half those of wild-type CHO-K1 or of a pseudorevertant. The significance of this reduction was tested by chemical challenge: mutants are twofold more sensitive than wild type to reactive oxygen species generated by phenylbenzoquinone and to damage produced by the anticancer drug, cisplatin. These results suggest that System X-AG provides a significant portion of the glutamate used to energize the uptake of cystine required for the synthesis of glutathione.

Refractive errors are common eye disorders of public health importance worldwide. Ocular axial length (AL) is the major determinant of refraction and thus of myopia and hyperopia. We conducted a meta-analysis of genome-wide association studies for AL, combining 12,531 Europeans and 8,216 Asians. We identified eight genome-wide significant loci for AL (RSPO1, C3orf26, LAMA2, GJD2, ZNRF3, CD55, MIP, and ALPPL2) and confirmed one previously reported AL locus (ZC3H11B). Of the nine loci, five (LAMA2, GJD2, CD55, ALPPL2, and ZC3H11B) were associated with refraction in 18 independent cohorts (n = 23,591). Differential gene expression was observed for these loci in minus-lens-induced myopia mouse experiments and human ocular tissues. Two of the AL genes, RSPO1 and ZNRF3, are involved in Wnt signaling, a pathway playing a major role in the regulation of eyeball size. This study provides evidence of shared genes between AL and refraction, but importantly also suggests that these traits may have unique pathways.

Two new Chinese hamster ovary cell (CHO-K1) mutants lacking amino acid transport System X-AG activity were isolated by [3H]aspartate suicide selection. These null mutants, Dd-B6 and Dd-B7, were analyzed by somatic cell hybridization, along with previously described partial-function mutants, Ed-A1 and Ed-B8. With respect to System X-AG activity, all four mutations fell into a single complementation group. By quantitative assay, the mutations in Ed-A1 and Ed-B8 behaved as simple recessives in fusions with wild type cells, while those in Dd-B6 and Dd-B7 were codominant. We have discovered that Ed-A1 and Ed-B8 are highly permeable to small neutral molecules. This high permeability phenotype was dominant to wild-type. Northern, Southern, and Western analyses indicated that System X-AG in CHO is not closely related to any of the three well characterized glutamate transporters represented by GLT-1, EAACI or GLAST.

Myopia is the most common human eye disorder and it results from complex genetic and environmental causes. The rapidly increasing prevalence of myopia poses a major public health challenge. Here, the CREAM consortium performs a joint meta-analysis to test single-nucleotide polymorphism (SNP) main effects and SNP × education interaction effects on refractive error in 40,036 adults from 25 studies of European ancestry and 10,315 adults from 9 studies of Asian ancestry. In European ancestry individuals, we identify six novel loci (FAM150B-ACP1, LINC00340, FBN1, DIS3L-MAP2K1, ARID2-SNAT1 and SLC14A2) associated with refractive error. In Asian populations, three genome-wide significant loci AREG, GABRR1 and PDE10A also exhibit strong interactions with education (P, This report by the Consortium for Refractive Error and Myopia uses gene-environment-wide interaction study (GEWIS) to identify genetic loci that affect environmental influence in myopia development, and identifies ethnic specific genetic loci that attribute to eye refractive errors.

The aim of this chapter is to introduce the reader to commonly used software packages and illustrate their input requirements, analysis options, strengths, and limitations. We focus on packages that perform more than one function and include a program for quality control, linkage, and association analyses. Additional inclusion criteria were (1) programs that are free to academic users and (2) currently supported, maintained, and developed. Using those criteria, we chose to review three programs: Statistical Analysis for Genetic Epidemiology (S.A.G.E.), PLINK, and Merlin. We will describe the required input format and analysis options. We will not go into detail about every possible program in the packages, but we will give an overview of the packages requirements and capabilities.

The aim of this chapter is to introduce the reader to commonly used software packages and illustrate their input requirements, analysis options, strengths, and limitations. We focus on packages that perform more than one function and include a program for quality control, linkage, and association analyses. Additional inclusion criteria were (1) programs that are free to academic users and (2) currently supported, maintained, and developed. Using those criteria, we chose to review three programs: Statistical Analysis for Genetic Epidemiology (S.A.G.E.), PLINK, and Merlin. We will describe the required input format and analysis options. We will not go into detail about every possible program in the packages, but we will give an overview of the packages requirements and capabilities.

For both model-free and model-based linkage analysis the S.A.G.E. (Statistical Analysis for Genetic Epidemiology) program package has some unique capabilities in analyzing both continuous traits and binary traits with variable age of onset. Here we highlight model-based linkage analysis of a quantitative trait (plasma dopamine β hydroxylase) that is known to be largely determined by monogenic inheritance, using a prior segregation analysis to produce the best fitting model for the trait. For a binary trait with variable age of onset (schizophrenia), we illustrate how using age of onset information to obtain a quantitative susceptibility trait leads to more statistically significant linkage signals, suggesting better power.

To describe the methods for family and case-control recruitment for a multicenter genetic and associated heritability analyses of Fuchs endothelial corneal dystrophy (FECD).Twenty-nine enrolling sites with 62 trained investigators and coordinators gathered individual and family information, graded the phenotype, and collected blood and/or saliva for genetic analysis on all individuals with and without FECD. The degree of FECD was assessed in a 0 to 6 semiquantitative scale using standardized clinical methods with pathological verification of FECD on at least 1 member of each family. Central corneal thickness was measured by ultrasonic pachymetry.Three hundred twenty-two families with 330 affected sibling pairs with FECD were enrolled and included a total of 650 sibling pairs of all disease grades. Using the entire 7-step FECD grading scale or a dichotomous definition of severe disease, heritability was assessed in families via sib-sib correlations. Both binary indicators of severe disease and semiquantitative measures of disease severity were significantly heritable, with heritability estimates of 30% for severe disease, 37% to 39% for FECD score, and 47% for central corneal thickness.Genetic risk factors have a strong role in the severity of the FECD phenotype and corneal thickness. Genotyping this cohort with high-density genetic markers followed by appropriate statistical analyses should lead to novel loci for disease susceptibility.

Background: Diabetic nephropathy (DN) is a leading cause of mortality and morbidity in patients with type 1 and type 2 diabetes. The multicenter FIND consortium aims to identify genes for DN and its associated quantitative traits, e.g. the urine albumin:creatinine ratio (ACR). Herein, the results of whole-genome linkage analysis and a sparse association scan for ACR and a dichotomous DN phenotype are reported in diabetic individuals. Methods: A genomewide scan comprising more than 5,500 autosomal single nucleotide polymorphism markers (average spacing of 0.6 cM) was performed on 1,235 nuclear and extended pedigrees (3,972 diabetic participants) ascertained for DN from African-American (AA), American-Indian (AI), European-American (EA) and Mexican-American (MA) populations. Results: Strong evidence for linkage to DN was detected on chromosome 6p (p = 8.0 × 10–5, LOD = 3.09) in EA families as well as suggestive evidence for linkage to chromosome 7p in AI families. Regions on chromosomes 3p in AA, 7q in EA, 16q in AA and 22q in MA displayed suggestive evidence of linkage for urine ACR. The linkage peak on chromosome 22q overlaps the MYH9/APOL1 gene region, previously implicated in AA diabetic and nondiabetic nephropathies. Conclusion: These results strengthen the evidence for previously identified genomic regions and implicate several novel loci potentially involved in the pathogenesis of DN.

Background Abdominal aortic aneurysm (AAA) is a complex disorder with multiple genetic risk factors. Using affected relative pair linkage analysis, we previously identified an AAA susceptibility locus on chromosome 19q13. This locus has been designated as the AAA1 susceptibility locus in the Online Mendelian Inheritance in Man (OMIM) database. Methods Nine candidate genes were selected from the AAA1 locus based on their function, as well as mRNA expression levels in the aorta. A sample of 394 cases and 419 controls was genotyped for 41 SNPs located in or around the selected nine candidate genes using the Illumina GoldenGate platform. Single marker and haplotype analyses were performed. Three genes (CEBPG, PEPD and CD22) were selected for DNA sequencing based on the association study results, and exonic regions were analyzed. Immunohistochemical staining of aortic tissue sections from AAA and control individuals was carried out for the CD22 and PEPD proteins with specific antibodies. Results Several SNPs were nominally associated with AAA (p < 0.05). The SNPs with most significant p-values were located near the CCAAT enhancer binding protein (CEBPG), peptidase D (PEPD), and CD22. Haplotype analysis found a nominally associated 5-SNP haplotype in the CEBPG/PEPD locus, as well as a nominally associated 2-SNP haplotype in the CD22 locus. DNA sequencing of the coding regions revealed no variation in CEBPG. Seven sequence variants were identified in PEPD, including three not present in the NCBI SNP (dbSNP) database. Sequencing of all 14 exons of CD22 identified 20 sequence variants, five of which were in the coding region and six were in the 3'-untranslated region. Five variants were not present in dbSNP. Immunohistochemical staining for CD22 revealed protein expression in lymphocytes present in the aneurysmal aortic wall only and no detectable expression in control aorta. PEPD protein was expressed in fibroblasts and myofibroblasts in the media-adventitia border in both aneurysmal and non-aneurysmal tissue samples. Conclusions Association testing of the functional positional candidate genes on the AAA1 locus on chromosome 19q13 demonstrated nominal association in three genes. PEPD and CD22 were considered the most promising candidate genes for altering AAA risk, based on gene function, association evidence, gene expression, and protein expression.

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), causes 9 million new cases worldwide and 2 million deaths annually. Genetic linkage and association analyses have suggested several chromosomal regions and candidate genes involved in TB susceptibility. This study examines the association of TB disease susceptibility with a selection of biologically relevant genes on regions on chromosomes 7 (IL6 and CARD11) and 20 (CTSZ and MC3R) and fine mapping of the chromosome 7p22–p21 region identified through our genome scan. We analyzed 565 individuals from Kampala, Uganda, who were previously included in our genome-wide linkage scan. Association analyses were conducted for 1,417 single-nucleotide polymorphisms (SNP) that passed quality control. None of the candidate gene or fine mapping SNPs was significantly associated with TB susceptibility (p > 0.10). When we restricted the analysis to HIV-negative individuals, 2 SNPs on chromosome 7 were significantly associated with TB susceptibility (p < 0.05). Haplotype analyses identified a significant risk haplotype in cathepsin X (CTSZ; p = 0.0281, odds ratio = 1.5493, 95% confidence interval [1.039, 2.320]).

Intermediate fine mapping has received considerable attention recently, with the goal of providing statistically precise and valid chromosomal regions for fine mapping following initial identification of broad regions that are linked to a disease. The following classes of methods have been proposed and compared in the literature: (1) LOD-support intervals, (2) generalized estimating equations, (3) bootstrap, and (4) confidence set inference framework. These methods provide confidence intervals either with coverage levels deviating from the nominal confidence levels or that are not fully efficient. Here, we propose a novel Bayesian method for constructing such intervals using affected sibling pair data. The susceptibility gene location is treated as a parameter in this method, with a uniform prior. A Metropolis-Hastings algorithm is implemented to sample from the posterior distribution and highest posterior density intervals of the disease gene locations are constructed. Correct coverage levels are maintained by our method. Both simulation studies and an application to a rheumatoid arthritis dataset demonstrate the improved efficiency of the Bayesian intervals compared with existing methods.

As alluded to in the chapter “Linkage Analysis of Qualitative Traits”, neither the Elston–Steward algorithm nor the Lander–Green approach is amenable to genetic data from large complex pedigrees and a large number of markers. In such cases, Monte Carlo estimation methods provide a viable alternative to the exact solutions. Two types of Monte Carlo methods have been developed for linkage analysis, haplotype inference, and other kinds of genetic analysis. They are Markov chain Monte Carlo (MCMC) methods and Monte Carlo methods that are based on independent samples. Approaches based on Markov chain Monte Carlo methods are more widely applicable; there is practically no limit on the size or complexity of the pedigrees, nor on the number of markers to be considered simultaneously. In this chapter, we will review the basic principles of MCMC methods for multipoint linkage analysis with extended pedigrees. Both simulations and application to data from the Framingham study will be used to compare and contrast three MCMC software packages: LOKI, MORGAN, and SIMWALK.