Genome‐wide association study of word reading: Overlap with risk genes for neurodevelopmental disorders

Reading disabilities (RD) are the most common neurocognitive disorder, affecting 5% to 17% of children in North America. These children often have comorbid neurodevelopmental/psychiatric disorders, such as attention deficit/hyperactivity disorder (ADHD). The genetics of RD and their overlap with other disorders is incompletely understood. To contribute to this, we performed a genome‐wide association study (GWAS) for word reading. Then, using summary statistics from neurodevelopmental/psychiatric disorders, we computed polygenic risk scores (PRS) and used them to predict reading ability in our samples. This enabled us to test the shared aetiology between RD and other disorders. The GWAS consisted of 5.3 million single nucleotide polymorphisms (SNPs) and two samples; a family‐based sample recruited for reading difficulties in Toronto (n = 624) and a population‐based sample recruited in Philadelphia [Philadelphia Neurodevelopmental Cohort (PNC)] (n = 4430). The Toronto sample SNP‐based analysis identified suggestive SNPs (P ~ 5 × 10^?7) in the ARHGAP23 gene, which is implicated in neuronal migration/axon pathfinding. The PNC gene‐based analysis identified significant associations (P < 2.72 × 10^?6) for LINC00935 and CCNT1, located in the region of the KANSL2/CCNT1/LINC00935/SNORA2B/SNORA34/MIR4701/ADCY6 genes on chromosome 12q, with near significant SNP‐based analysis. PRS identified significant overlap between word reading and intelligence (R² = 0.18, P = 7.25 × 10^?181), word reading and educational attainment (R² = 0.07, P = 4.91 × 10^?48) and word reading and ADHD (R² = 0.02, P = 8.70 × 10^?6; threshold for significance = 7.14 × 10^?3). Overlap was also found between RD and autism spectrum disorder (ASD) as top‐ranked genes were previously implicated in autism by rare and copy number variant analyses. These findings support shared risk between word reading, cognitive measures, educational outcomes and neurodevelopmental disorders, including ASD.     

Modifier gene study of meconium ileus in cystic fibrosis: statistical considerations and gene mapping results

Cystic fibrosis (CF) is a monogenic disease due to mutations in the CFTR gene. Yet, variability in CF disease presentation is presumed to be affected by modifier genes, such as those recently demonstrated for the pulmonary aspect. Here, we conduct a modifier gene study for meconium ileus (MI), an intestinal obstruction that occurs in 16–20% of CF newborns, providing linkage and association results from large family and case–control samples. Linkage analysis of modifier traits is different than linkage analysis of primary traits on which a sample was ascertained. Here, we articulate a source of confounding unique to modifier gene studies and provide an example of how one might overcome the confounding in the context of linkage studies. Our linkage analysis provided evidence of a MI locus on chromosome 12p13.3, which was segregating in up to 80% of MI families with at least one affected offspring (HLOD = 2.9). Fine mapping of the 12p13.3 region in a large case–control sample of pancreatic insufficient Canadian CF patients with and without MI pointed to the involvement of ADIPOR2 in MI (p = 0.002). This marker was substantially out of Hardy–Weinberg equilibrium in the cases only, and provided evidence of a cohort effect. The association with rs9300298 in the ADIPOR2 gene at the 12p13.3 locus was replicated in an independent sample of CF families. A protective locus, using the phenotype of no-MI, mapped to 4q13.3 (HLOD = 3.19), with substantial heterogeneity. A candidate gene in the region, SLC4A4, provided preliminary evidence of association (p = 0.002), warranting further follow-up studies. Our linkage approach was used to direct our fine-mapping studies, which uncovered two potential modifier genes worthy of follow-up.     

Analogs of RGDVY and GRGD peptides inhibit Mycobacterium kansasii phagocytosis

Continuing our research on Mycobacteria kansasii phagocytosis inhibition, we have examined in that context three series of peptides derived from the RGDVY and GRGD sequences. It was found that the levels of the inhibitory activity depend on the amino acid composition as well as on the particular peptide sequence. Distinct inhibitory activity was found in the case of thymopentin (RKDVY), the active fragment of thymopoietin. In this case the Mycobacterium phagocytosis inhibition should be combined with general immunostimulatory activity of RKDVY peptide. Our examination of a series of GRGDV analogs with a successively prolonged oligo-Gly linker inserted into the peptide chain showed that the distance between the Arg and Asp residues required for such an activity should be about 9A.     

The genetics of cross-sectional and longitudinal body mass index

There has been a lack of consistency in detecting chromosomal loci that are linked to obesity-related traits. This may be due, in part, to the phenotype definition. Many studies use a one-time, single measurement as a phenotype while one's weight often fluctuates considerably throughout adulthood. Longitudinal data from the Framingham Heart Study were used to derive alternative phenotypes that may lead to more consistent findings. Body mass index (BMI), a measurement for obesity, is known to increase with age and then plateau or decline slightly; the decline phase may represent a threshold or survivor effect. We propose to use the weight gain phase of BMI to derive phenotypes useful for linkage analysis of obesity. Two phenotypes considered in the present study are the average of and the slope of the BMI measurements in the gain phase (gain mean and gain slope). For comparison, we also considered the average of all BMI measurements available (overall mean). Linkage analysis using the gain mean phenotype exhibited two markers with LOD scores greater than 3, with the largest score of 3.52 on chromosome 4 at ATA2A03. In contrast, no LOD scores greater than 3 were observed when overall mean was used. The gain slope produced weak evidence for linkage on chromosome 4 with a multipoint LOD score of 1.77 at GATA8A05. Our analysis shows how omitting the decline phase of BMI in the definition of obesity phenotypes can result in evidence for linkage which might have been otherwise overlooked.     

Construction of the model for the Genetic Analysis Workshop 14 simulated data: genotype-phenotype relationships,gene interaction,linkage, association,disequilibrium, and ascertainment effects for a complex phenotype

The Genetic Analysis Workshop 14 simulated dataset was designed 1) To test the ability to find genes related to a complex disease (such as alcoholism). Such a disease may be given a variety of definitions by different investigators, have associated endophenotypes that are common in the general population, and is likely to be not one disease but a heterogeneous collection of clinically similar, but genetically distinct, entities. 2) To observe the effect on genetic analysis and gene discovery of a complex set of gene x gene interactions. 3) To allow comparison of microsatellite vs. large-scale single-nucleotide polymorphism (SNP) data. 4) To allow testing of association to identify the disease gene and the effect of moderate marker x marker linkage disequilibrium. 5) To observe the effect of different ascertainment/disease definition schemes on the analysis. Data was distributed in two forms. Data distributed to participants contained about 1,000 SNPs and 400 microsatellite markers. Internet-obtainable data consisted of a finer 10,000 SNP map, which also contained data on controls. While disease characteristics and parameters were constant, four "studies" used varying ascertainment schemes based on differing beliefs about disease characteristics. One of the studies contained multiplex two- and three-generation pedigrees with at least four affected members. The simulated disease was a psychiatric condition with many associated behaviors (endophenotypes), almost all of which were genetic in origin. The underlying disease model contained four major genes and two modifier genes. The four major genes interacted with each other to produce three different phenotypes, which were themselves heterogeneous. The population parameters were calibrated so that the major genes could be discovered by linkage analysis in most datasets. The association evidence was more difficult to calibrate but was designed to find statistically significant association in 50% of datasets. We also simulated some marker x marker linkage disequilibrium around some of the genes and also in areas without disease genes. We tried two different methods to simulate the linkage disequilibrium.     

Length of the peptide chain influences the immunomodulatory activity of peptides related to p53 protein.

We have shown that the thymopoietin-like octapeptides derived from DNA-binding domain of p53 protein and of its mutated forms differ in their immunomodulatory properties. A strong increase of immunostimulative activity was observed for GMNRSPIL (mutated protein) in comparison with GMNRRPIL (wild-type of p53 protein) peptide. Here the elongated sequences of respective protein fragments were synthesized and investigated by plaque forming cells and delayed type hypersensitivity tests. The change of immunomodulatory activity toward immunosuppression was observed: NSSC(Acm)MGGMNRRPILTIITLE (1, wild-type) was inactive in both tests, and the C(Acm)MGGMNRSPILTIITLE (II) and YMC(Acm)NSSC(Acm)MGGMNRSPILTIITLE (III) (mutated p53 protein fragments) peptides produced immunosuppression in plaque forming cells as well as in delayed type hypersensitivity tests.     

Malic enzyme 2 may underlie susceptibility to adolescent-onset idiopathic generalized epilepsy

Idiopathic generalized epilepsy (IGE) is a class of genetically determined, phenotypically related epilepsy syndromes. Linkage analysis identified a chromosome 18 locus predisposing to a number of adolescent-onset IGEs. We report a single-nucleotide polymorphism (SNP) association analysis of the region around the marker locus with the high LOD score. This analysis, which used both case-control and family-based association methods, yielded strong evidence that malic enzyme 2 (ME2) is the gene predisposing to IGE. We also observed association among subgroups of IGE syndromes. An ME2-centered nine-SNP haplotype, when present homozygously, increases the risk for IGE (odds ratio 6.1; 95% confidence interval 2.9-12.7) compared with any other genotype. Both the linkage analysis and the association analysis support recessive inheritance for the locus, which is compatible with the fact that ME2 is an enzyme. ME2 is a genome-coded mitochondrial enzyme that converts malate to pyruvate and is involved in neuronal synthesis of the neurotransmitter gamma-aminobutyric acid (GABA). The results suggest that GABA synthesis disruption predisposes to common IGE and that clinical seizures are triggered when mutations at other genes, or perhaps other insults, are present.     

Unraveling the complex genetic model for cystic fibrosis: pleiotropic effects of modifier genes on early cystic fibrosis-related morbidities

The existence of pleiotropy in disorders with multi-organ involvement can suggest therapeutic targets that could ameliorate overall disease severity. Here we assessed pleiotropy of modifier genes in cystic fibrosis (CF). CF, caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, affects the lungs, liver, pancreas and intestines. However, modifier genes contribute to variable disease severity across affected organs, even in individuals with the same CFTR genotype. We sought to determine whether SLC26A9, SLC9A3 and SLC6A14, that contribute to meconium ileus in CF, are pleiotropic for other early-affecting CF co-morbidities. In the Canadian CF population, we assessed evidence for pleiotropic effects on (1) pediatric lung disease severity (n = 815), (2) age at first acquisition of Pseudomonas aeruginosa (P. aeruginosa) (n = 730), and (3) prenatal pancreatic damage measured by immunoreactive trypsinogen (n = 126). A multiple-phenotype analytic strategy assessed evidence for pleiotropy in the presence of phenotypic correlation. We required the same alleles to be associated with detrimental effects. SLC26A9 was pleiotropic for meconium ileus and pancreatic damage (p = 0.002 at rs7512462), SLC9A3 for meconium ileus and lung disease (p = 1.5 × 10^?6 at rs17563161), and SLC6A14 for meconium ileus and both lung disease and age at first P. aeruginosa infection (p = 0.0002 and p = 0.006 at rs3788766, respectively). The meconium ileus risk alleles in SLC26A9, SLC9A3 and SLC6A14 are pleiotropic, increasing risk for other early CF co-morbidities. Furthermore, co-morbidities affecting the same organ tended to associate with the same genes. The existence of pleiotropy within this single disorder suggests that complementary therapeutic strategies to augment solute transport will benefit multiple CF-associated tissues.     

Evidence for linkage of migraine in Rolandic epilepsy to known 1q23 FHM2 and novel 17q22 genetic loci

Migraine headaches are a common comorbidity in Rolandic epilepsy (RE) and familial aggregation of migraine in RE families suggests a genetic basis not mediated by seizures. We performed a genome‐wide linkage analysis of the migraine phenotype in 38 families with RE to localize potential genetic contribution, with a follow‐up in an additional 21 families at linked loci. We used two‐point and multipoint LOD (logarithm of the odds) score methods for linkage, maximized over genetic models. We found evidence of linkage to migraine at chromosome 17q12‐22 [multipoint HLOD (heterogeneity LOD) 4.40, recessive, 99% penetrance], replicated in the second dataset (HLOD 2.61), and suggestive evidence at 1q23.1‐23.2, centering over the FHM2 locus (two‐point LOD 3.00 and MP HLOD 2.52). Sanger sequencing in 14 migraine‐affected individuals found no coding mutations in the FHM2 gene ATP1A2. There was no evidence of pleiotropy for migraine and either reading or speech disorder, or the electroencephalographic endophenotype of RE when the affected definition was redefined as those with migraine or the comorbid phenotype, and pedigrees were reanalyzed for linkage. In summary, we report a novel migraine susceptibility locus at 17q12‐22, and a second locus that may contribute to migraine in the general population at 1q23.1‐23.2. Comorbid migraine in RE appears genetically influenced, but we did not obtain evidence that the identified susceptibility loci are consistent with pleiotropic effects on other comorbidities in RE. Loci identified here should be fine‐mapped in individuals from RE families with migraine, and prioritized for analysis in other types of epilepsy‐associated migraine.     

An alternative foundation for the planning and evaluation of linkage analysis. II. Implications for multiple test adjustments

The 'multiple testing problem' currently bedevils the field of genetic epidemiology. Briefly stated, this problem arises with the performance of more than one statistical test and results in an increased probability of committing at least one Type I error. The accepted/conventional way of dealing with this problem is based on the classical Neyman-Pearson statistical paradigm and involves adjusting one's error probabilities. This adjustment is, however, problematic because in the process of doing that, one is also adjusting one's measure of evidence. Investigators have actually become wary of looking at their data, for fear of having to adjust the strength of the evidence they observed at a given locus on the genome every time they conduct an additional test. In a companion paper in this issue (Strug & Hodge I), we presented an alternative statistical paradigm, the 'evidential paradigm', to be used when planning and evaluating linkage studies. The evidential paradigm uses the lod score as the measure of evidence (as opposed to a p value), and provides new, alternatively defined error probabilities (alternative to Type I and Type II error rates). We showed how this paradigm separates or decouples the two concepts of error probabilities and strength of the evidence. In the current paper we apply the evidential paradigm to the multiple testing problem - specifically, multiple testing in the context of linkage analysis. We advocate using the lod score as the sole measure of the strength of evidence; we then derive the corresponding probabilities of being misled by the data under different multiple testing scenarios. We distinguish two situations: performing multiple tests of a single hypothesis, vs. performing a single test of multiple hypotheses. For the first situation the probability of being misled remains small regardless of the number of times one tests the single hypothesis, as we show. For the second situation, we provide a rigorous argument outlining how replication samples themselves (analyzed in conjunction with the original sample) constitute appropriate adjustments for conducting multiple hypothesis tests on a data set.