Fuchs endothelial corneal dystrophy (FECD) is the most common late-onset, vision-threatening corneal dystrophy in the United States, affecting about 4% of the population. Advanced FECD involves a thickening of the cornea from stromal edema and changes in Descemet membrane. To understand the relationship between FECD and central corneal thickness (CCT), we characterized common genetic variation in COL8A2 and TCF4, genes previously implicated in CCT and/or FECD. Other genes previously associated with FECD (PITX2, ZEB1, SLC4A11), and genes only known to affect CCT (COL5A1, FOXO1, AVGR8, ZNF469) were also interrogated. FECD probands, relatives and controls were recruited from 32 clinical sites; a total of 532 cases and 204 controls were genotyped and tested for association of FECD case/control status, a 7-step FECD severity scale and CCT, adjusting for age and sex. Association of FECD grade with TCF4 was highly significant (OR  = 6.01 at rs613872; p = 4.8×10−25), and remained significant when adjusted for changes in CCT (OR  = 4.84; p = 2.2×10−16). Association of CCT with TCF4 was also significant (p = 6.1×10−7), but was abolished with adjustment for FECD grade (p = 0.92). After adjusting for FECD grade, markers in other genes examined were modestly associated (p ∼ 0.001) with FECD and/or CCT. Thus, common variants in TCF4 appear to influence FECD directly, and CCT secondarily via FECD. Additionally, changes in corneal thickness due to the effect of other loci may modify disease severity, age-at-onset, or other biomechanical characteristics.

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 1 896 cases and 1 866 controls and identified two novel genetic protective factors for AMD. In addition to the CFH (p=2.3×10−64) and ARMS2 (p=1.2×10−60) loci, we observed a protective effect at rs429608, an intronic SNP in SKIV2L (p=5.3×10−15), a gene near the C2/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−4), 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.

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.

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 1417 single-nucleotide polymorphisms (SNPs) that passed quality control. None of the candidate gene or fine mapping SNPs were found significantly associated with TB susceptibility (P > 0.10). When we restricted the analysis to HIV-negative individuals, two 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, OR = 1.5493, 95% CI [1.039, 2.320]).

Complement factor H shows very strong association with Age-related Macular Degeneration (AMD), and recent data suggest that multiple causal variants are associated with disease. To refine the location of the disease associated variants, we characterized in detail the structural variation at CFH and its paralogs, including two copy number polymorphisms (CNP), CNP147 and CNP148, and several rare deletions and duplications. Examination of 34 AMD-enriched extended families (N = 293) and AMD cases (White N = 4210 Indian = 134; Malay = 140) and controls (White N = 3229; Indian = 117; Malay = 2390) demonstrated that deletion CNP148 was protective against AMD, independent of SNPs at CFH. Regression analysis of seven common haplotypes showed three haplotypes, H1, H6 and H7, as conferring risk for AMD development. Being the most common haplotype H1 confers the greatest risk by increasing the odds of AMD by 2.75-fold (95% CI = [2.51, 3.01]; p = 8.31×10−109); Caucasian (H6) and Indian-specific (H7) recombinant haplotypes increase the odds of AMD by 1.85-fold (p = 3.52×10−9) and by 15.57-fold (P = 0.007), respectively. We identified a 32-kb region downstream of Y402H (rs1061170), shared by all three risk haplotypes, suggesting that this region may be critical for AMD development. Further analysis showed that two SNPs within the 32 kb block, rs1329428 and rs203687, optimally explain disease association. rs1329428 resides in 20 kb unique sequence block, but rs203687 resides in a 12 kb block that is 89% similar to a noncoding region contained in ΔCNP148. We conclude that causal variation in this region potentially encompasses both regulatory effects at single markers and copy number.

Familial combined hyperlipidemia (FCHL) is a complex trait leading to cardiovascular disease (CVD) risk. Elevated levels and size of apolipoprotein B (apoB) and low-density lipoprotein (LDL) are associated with FCHL, which is genetically heterogeneous and is likely caused by rare variants. We carried out a linkage-based genome scan of four large FCHL pedigrees for apoB level that is independent of LDL: apoB level that is adjusted for LDL level and size. Follow-up included SNP genotyping in the region with the strongest evidence of linkage. Several regions with the evidence of linkage in individual pedigrees support the rare variant model. Evidence of linkage was strongest on chromosome 4q, with multipoint analysis in one pedigree giving LOD = 3.1 with a parametric model, and a log Bayes Factor = 1.5 from a Bayesian oligogenic approach. Of the 293 SNPs spanning the implicated region on 4q, rs6829588 completely explained the evidence of linkage. This SNP accounted for 39% of the apoB phenotypic variance, with heterozygotes for this SNP having a trait value that was ~30% higher than that of the high-frequency homozygote, thus identifying and considerably refining a strong candidate region. These results illustrate the advantage of using large pedigrees in the search for rare variants: reduced genetic heterogeneity within single pedigrees coupled with the large number of individuals segregating otherwise-rare single variants leads to high power to implicate such variants.

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.

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. Nonparametric linkage analysis, using an average of 4,404 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 chromosomes 4q21.1 (LOD=3.13; genome-wide p=0.04) in AA. In addition, a total of eleven 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.

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, (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.

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. 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.

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 LODmax 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.

PURPOSE

Diabetic retinopathy (DR) and diabetic nephropathy (DN) are serious microvascular complications of diabetes mellitus. Correlations between severity of DR and DN and computed heritability estimates for DR were determined in a large, multiethnic sample of diabetic families. The hypothesis was that (1) the severity of DR correlates with the presence and severity of nephropathy in individuals with diabetes mellitus, and (2) the severity of DR is under significant familial influence in members of multiplex diabetic families.

METHODS

The Family Investigation of Nephropathy and Diabetes (FIND) was designed to evaluate the genetic basis of DN in American Indians, European Americans, African Americans, and Mexican Americans. FIND enrolled probands with advanced DN, along with their diabetic siblings who were concordant and discordant for nephropathy. These diabetic family members were invited to participate in the FIND-Eye study to determine whether inherited factors underlie susceptibility to DR and its severity. FIND-Eye participants underwent eye examinations and had fundus photographs taken. The severity of DR was graded by using the Early Treatment Diabetic Retinopathy Study Classification (ETDRS). Sib–sib correlations were calculated with the SAGE 5.0 program FCOR, to estimate heritability of retinopathy severity.

RESULTS

This report summarizes the results for the first 2368 diabetic subjects from 767 families enrolled in FIND-Eye; nearly 50% were Mexican American, the largest single ethnicity within FIND. The overall prevalence of DR was high; 33.4% had proliferative DR; 7.5%, 22.8%, and 9.5% had severe, moderate, and mild nonproliferative DR, respectively; 26.6% had no DR. The severity of DR was significantly associated with severity of DN, both by phenotypic category and by increasing serum creatinine concentration (χ2 = 658.14, df = 20; P < 0.0001). The sib–sib correlation for DR severity was 0.1358 in the total sample and 0.1224 when limited to the Mexican-American sample. Broad sense heritabilities for DR were 27% overall and 24% in Mexican-American families. The polygenic heritability of liability for proliferative DR approximated 25% in this FIND-Eye sample.

CONCLUSIONS

These data confirm that the severity of DR parallels the presence and severity of nephropathy in individuals with diabetes mellitus. The severity of DR in members of multiplex diabetic families appears to have a significant familial connection.

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 followup 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 of 53.2 (chromosome-wide p < 0.004), a VC LOD score of 2.29, and a parametric LOD score of 2.94, based on a quantitative-trait model from MCMC segregation analysis. A weaker signal for SWE on chromosome 2q occurred in the same location as a significant linkage peak seen previously in a scan for phonological decoding. MCMC oligogenic segregation analysis identified three models of transmission for WID, which could be assigned to two distinct linkage peaks on chromosomes 12 and 15. Taken together, these results indicate a locus for efficiency and accuracy of RWR on chromosome 13, and a complex model for inheritance of RWR accuracy with loci on chromosomes 12 and 15.

Clustering of related haplotypes in haplotype-based association mapping has the potential to improve power by reducing the degrees of freedom without sacrificing important information about the underlying genetic structure. We have modified a generalized linear model approach for association analysis by incorporating a density-based clustering algorithm to reduce the number of coefficients in the model. Using the GAW 15 Problem 3 simulated data, we show that our novel method can substantially enhance power to detect association with the binary rheumatoid arthritis (RA) phenotype at the HLA-DRB1 locus on chromosome 6. In contrast, clustering did not appreciably improve performance at locus D, perhaps a consequence of a rare susceptibility allele and of the overwhelming effect of HLA-DRB1/locus C, 5 cM distal. Optimization of parameters governing the clustering algorithm identified a set of parameters that delivered nearly ideal performance in a variety of situations. The cluster-based score test was valid over a wide range of haplotype diversity, and was robust to severe departures from Hardy-Weinberg equilibrium encountered near HLA-DRB1 in RA case-control samples.

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.

Editing of mitochondrial mRNAs in kinetoplastid protozoa occurs by a series of enzymatic steps that insert and delete uridylates (U's) as specified by guide RNAs (gRNAs). The characteristics of the 3" exonuclease activity that removes the U's following cleavage during deletion editing were determined by using an in vitro precleaved deletion assay that is based on ATPase subunit 6 pre-mRNA and gA6[14] gRNA. The exonuclease in partially purified editing complexes is specific for U's. The specificity occurs in the absence of gRNA, but its activity is enhanced by the presence of gRNA. The 3" pre-mRNA fragment enhances the specificity, but not the efficiency, of U removal. The activity is sensitive to the 5" phosphate of the 3" fragment, which is not required for U removal. The ability of the 3" U's to base pair with purines in the gRNA protects them from removal, suggesting that the U-specific 3" exonuclease (exoUase) is specific for U's which are not base paired. ExoUase is stereospecific and cannot remove (Rp)α-thio-U. The specificity of the exoUase activity thus contributes to the precision of RNA editing.

RNA editing in kinetoplastid mitochondria occurs by a series of enzymatic steps that is catalyzed by a macromolecular complex. Four novel proteins and their corresponding genes were identified by mass spectrometric analysis of purified editing complexes from Trypanosoma brucei. These four proteins, TbMP81, TbMP63, TbMP42, and TbMP18, contain conserved sequences to various degrees. All four proteins have sequence similarity in the C terminus; TbMP18 has considerable sequence similarity to the C-terminal region of TbMP42, and TbMP81, TbMP63, and TbMP42 contain zinc finger motif(s). Monoclonal antibodies that are specific for TbMP63 and TbMP42 immunoprecipitate in vitro RNA editing activities. The proteins are present in the immunoprecipitates and sediment at 20S along with the in vitro editing, and RNA editing ligases TbMP52 and TbMP48. Recombinant TbMP63 and TbMP52 coimmunoprecipitate. These results indicate that these four proteins are components of the RNA editing complex and that TbMP63 and TbMP52 can interact.

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 ∼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.

RNA editing in Trypanosoma brucei inserts and deletes uridylates (U's) in mitochondrial pre-mRNAs under the direction of guide RNAs (gRNAs). We report here the development of a novel in vitro precleaved editing assay and its use to study the gRNA specificity of the U addition and RNA ligation steps in insertion RNA editing. The 5′ fragment of substrate RNA accumulated with the number of added U's specified by gRNA, and U addition products with more than the specified number of U's were rare. U addition up to the number specified occurred in the absence of ligation, but accumulation of U addition products was slowed. The 5′ fragments with the correct number of added U's were preferentially ligated, apparently by adenylylated RNA ligase since exogenously added ATP was not required and since ligation was eliminated by treatment with pyrophosphate. gRNA-specified U addition was apparent in the absence of ligation when the pre-mRNA immediately upstream of the editing site was single stranded and more so when it was base paired with gRNA. These results suggest that both the U addition and RNA ligation steps contributed to the precision of RNA editing.