Common genetic risk factors in ASD and ADHD co-occurring families.

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are two major neurodevelopmental disorders that frequently co-occur. However, the genetic mechanism of the co-occurrence remains unclear. The New Jersey Language and Autism Genetics Study (NJLAGS) collected more than 100 families with at least one member affected by ASD. NJLAGS families show a high prevalence of ADHD and provide a good opportunity to study shared genetic risk factors for ASD and ADHD. The linkage study of the NJLAGS families revealed regions on chromosomes 12 and 17 that are significantly associated with ADHD. Using whole-genome sequencing data on 272 samples from 73 NJLAGS families, we identified potential risk genes for ASD and ADHD. Within the linkage regions, we identified 36 genes that are associated with ADHD using a pedigree-based gene prioritization approach. KDM6B (Lysine Demethylase 6B) is the highest-ranking gene, which is a known risk gene for neurodevelopmental disorders, including ASD and ADHD. At the whole-genome level, we identified 207 candidate genes from the analysis of both small variants and structure variants, including both known and novel genes. Using enrichment and protein-protein interaction network analyses, we identified gene ontology terms and pathways enriched for ASD and ADHD candidate genes, such as cilia function and cation channel activity. Candidate genes and pathways identified in our study improve the understanding of the genetic etiology of ASD and ADHD and will lead to new diagnostic or therapeutic interventions for ASD and ADHD in the future.

Structural Variations Contribute to the Genetic Etiology of Autism Spectrum Disorder and Language Impairments.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by restrictive interests and/or repetitive behaviors and deficits in social interaction and communication. ASD is a multifactorial disease with a complex polygenic genetic architecture. Its genetic contributing factors are not yet fully understood, especially large structural variations (SVs). In this study, we aimed to assess the contribution of SVs, including copy number variants (CNVs), insertions, deletions, duplications, and mobile element insertions, to ASD and related language impairments in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Within the cohort, ~77% of the families contain SVs that followed expected segregation or de novo patterns and passed our filtering criteria. These SVs affected 344 brain-expressed genes and can potentially contribute to the genetic etiology of the disorders. Gene Ontology and protein-protein interaction network analysis suggested several clusters of genes in different functional categories, such as neuronal development and histone modification machinery. Genes and biological processes identified in this study contribute to the understanding of ASD and related neurodevelopment disorders.

MiR-29b is associated with perinatal inflammation in extremely preterm infants.

BACKGROUND: Inflammation is strongly associated with premature birth and neonatal morbidities. Increases in infant haptoglobin, haptoglobin-related protein (Hp&HpRP), and interleukin-6 (IL-6) levels are indicators of intra-amniotic inflammation (IAI) and have been linked to poor neonatal outcomes. Inflammation causes epigenetic changes, specifically suppression of miR-29 expression. The current study sought to determine whether miR-29b levels in cord blood or neonatal venous blood are associated with IAI, identified by elevated IL-6 and Hp, and subsequent clinical morbidities in the infant. METHODS: We tested 92 cord blood samples from premature newborns and 18 venous blood samples at 36 weeks corrected gestational age. MiR-29b, Hp&HpRP, and IL-6 were measured by polymerase chain reaction and enzyme-linked immunosorbent assay, respectively. RESULTS: Decreased levels of miR-29b were observed in infants exposed to IAI with elevated Hp&HpRP and IL-6 levels and in infants delivered by spontaneous preterm birth. Lower miR-29 levels were also observed in women diagnosed with histological chorioamnionitis or funisitis and in infants with cerebral palsy. Higher levels of miR-29 were measured in infants small for gestational age and in venous samples from older infants. CONCLUSIONS: MiR-29 may be an additional biomarker of IAI and a potential therapeutic target for treating poor newborn outcomes resulting from antenatal exposure to IAI. IMPACT: Decreases in miR-29b are associated with intrauterine inflammation. Hp&HpRP increases are associated with decreased miR-29b. MiR-29b may be an additional biomarker for neonatal outcomes and a potential therapeutic target for intrauterine inflammation.

Behavioral and Molecular Genetics of Reading-Related AM and FM Detection Thresholds.

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h 2  = 0.20) and FM (h 2  = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.

An integrative analysis of regional gene expression profiles in the human brain.

Studies of the brain's transcriptome have become prominent in recent years, resulting in an accumulation of datasets with somewhat distinct attributes. These datasets, which are often analyzed only in isolation, also are often collected with divergent goals, which are reflected in their sampling properties. While many researchers have been interested in sampling gene expression in one or a few brain areas in a large number of subjects, recent efforts from the Allen Institute for Brain Sciences and others have focused instead on dense neuroanatomical sampling, necessarily limiting the number of individual donor brains studied. The purpose of the present work is to develop methods that draw on the complementary strengths of these two types of datasets for study of the human brain, and to characterize the anatomical specificity of gene expression profiles and gene co-expression networks derived from human brains using different specific technologies. The approach is applied using two publicly accessible datasets: (1) the high anatomical resolution Allen Human Brain Atlas (AHBA, Hawrylycz et al., 2012) and (2) a relatively large sample size, but comparatively coarse neuroanatomical dataset described previously by Gibbs et al. (2010). We found a relatively high degree of correspondence in differentially expressed genes and regional gene expression profiles across the two datasets. Gene co-expression networks defined in individual brain regions were less congruent, but also showed modest anatomical specificity. Using gene modules derived from the Gibbs dataset and from curated gene lists, we demonstrated varying degrees of anatomical specificity based on two classes of methods, one focused on network modularity and the other focused on enrichment of expression levels. Two approaches to assessing the statistical significance of a gene set's modularity in a given brain region were studied, which provide complementary information about the anatomical specificity of a gene network of interest. Overall, the present work demonstrates the feasibility of cross-dataset analysis of human brain microarray studies, and offers a new approach to annotating gene lists in a neuroanatomical context.

Molecular Genetic Evidence for Shared Etiology of Autism and Prodigy.

Child prodigies are rare individuals with an exceptional working memory and unique attentional skills that may facilitate the attainment of professional skill levels at an age well before what is observed in the general population. Some characteristics of prodigy have been observed to be quantitatively similar to those observed in autism spectrum disorder (ASD), suggesting possible shared etiology, though objectively validated prodigies are so rare that evidence has been sparse. We performed a family-based genome-wide linkage analysis on 5 nuclear and extended families to search for genetic loci that influence the presence of both prodigy and ASD, assuming that the two traits have the same genetic etiology in the analysis model in order to find shared loci. A shared locus on chromosome 1p31-q21 reached genome-wide significance with two extended family-based linkage methods consisting of the Bayesian PPL method and the LOD score maximized over the trait parameters (i.e., MOD), yielding a simulation-based empirical significance of p = 0.000742 and p = 0.000133, respectively. Within linkage regions, we performed association analysis and assessed if copy number variants could account for the linkage signal. No evidence of specificity for either the prodigy or the ASD trait was observed. This finding suggests that a locus on chromosome 1 increases the likelihood of both prodigy and autism in these families.

Addressing the unmet need for visualizing conditional random fields in biological data.

BACKGROUND: The biological world is replete with phenomena that appear to be ideally modeled and analyzed by one archetypal statistical framework - the Graphical Probabilistic Model (GPM). The structure of GPMs is a uniquely good match for biological problems that range from aligning sequences to modeling the genome-to-phenome relationship. The fundamental questions that GPMs address involve making decisions based on a complex web of interacting factors. Unfortunately, while GPMs ideally fit many questions in biology, they are not an easy solution to apply. Building a GPM is not a simple task for an end user. Moreover, applying GPMs is also impeded by the insidious fact that the "complex web of interacting factors" inherent to a problem might be easy to define and also intractable to compute upon. DISCUSSION: We propose that the visualization sciences can contribute to many domains of the bio-sciences, by developing tools to address archetypal representation and user interaction issues in GPMs, and in particular a variety of GPM called a Conditional Random Field(CRF). CRFs bring additional power, and additional complexity, because the CRF dependency network can be conditioned on the query data. CONCLUSIONS: In this manuscript we examine the shared features of several biological problems that are amenable to modeling with CRFs, highlight the challenges that existing visualization and visual analytics paradigms induce for these data, and document an experimental solution called StickWRLD which, while leaving room for improvement, has been successfully applied in several biological research projects. Software and tutorials are available at http://www.stickwrld.org/.

Who is afraid of math? Two sources of genetic variance for mathematical anxiety.

BACKGROUND: Emerging work suggests that academic achievement may be influenced by the management of affect as well as through efficient information processing of task demands. In particular, mathematical anxiety has attracted recent attention because of its damaging psychological effects and potential associations with mathematical problem solving and achievement. This study investigated the genetic and environmental factors contributing to the observed differences in the anxiety people feel when confronted with mathematical tasks. In addition, the genetic and environmental mechanisms that link mathematical anxiety with math cognition and general anxiety were also explored. METHODS: Univariate and multivariate quantitative genetic models were conducted in a sample of 514 12-year-old twin siblings. RESULTS: Genetic factors accounted for roughly 40% of the variation in mathematical anxiety, with the remaining being accounted for by child-specific environmental factors. Multivariate genetic analyses suggested that mathematical anxiety was influenced by the genetic and nonfamilial environmental risk factors associated with general anxiety and additional independent genetic influences associated with math-based problem solving. CONCLUSIONS: The development of mathematical anxiety may involve not only exposure to negative experiences with mathematics, but also likely involves genetic risks related to both anxiety and math cognition. These results suggest that integrating cognitive and affective domains may be particularly important for mathematics and may extend to other areas of academic achievement.

Gene-environment interplay in child psychology and psychiatry--challenges and ways forward.

This special issue in the Journal of Child Psychology and Psychiatry presents several invited articles examining gene-environment interplay in child development and psychopathology. Models of gene-environment interplay have been exhaustively discussed in the literature, including an important contribution by Rutter, Moffitt and Caspi (2006) published in this journal.

Gene × smoking interactions on human brain gene expression: finding common mechanisms in adolescents and adults.

BACKGROUND: Numerous studies have examined gene × environment interactions (G × E) in cognitive and behavioral domains. However, these studies have been limited in that they have not been able to directly assess differential patterns of gene expression in the human brain. Here, we assessed G × E interactions using two publically available datasets to assess if DNA variation is associated with post-mortem brain gene expression changes based on smoking behavior, a biobehavioral construct that is part of a complex system of genetic and environmental influences. METHODS: We conducted an expression quantitative trait locus (eQTL) study on two independent human brain gene expression datasets assessing G × E for selected psychiatric genes and smoking status. We employed linear regression to model the significance of the Gene × Smoking interaction term, followed by meta-analysis across datasets. RESULTS: Overall, we observed that the effect of DNA variation on gene expression is moderated by smoking status. Expression of 16 genes was significantly associated with single nucleotide polymorphisms that demonstrated G × E effects. The strongest finding (p = 1.9 × 10⁻¹¹) was neurexin 3-alpha (NRXN3), a synaptic cell-cell adhesion molecule involved in maintenance of neural connections (such as the maintenance of smoking behavior). Other significant G × E associations include four glutamate genes. CONCLUSIONS: This is one of the first studies to demonstrate G × E effects within the human brain. In particular, this study implicated NRXN3 in the maintenance of smoking. The effect of smoking on NRXN3 expression and downstream behavior is different based upon SNP genotype, indicating that DNA profiles based on SNPs could be useful in understanding the effects of smoking behaviors. These results suggest that better measurement of psychiatric conditions, and the environment in post-mortem brain studies may yield an important avenue for understanding the biological mechanisms of G × E interactions in psychiatry.

Evaluation of a bayesian model integration-based method for censored data.

OBJECTIVE: Non-random missing data can adversely affect family-based linkage detection through loss of power and possible introduction of bias depending on how censoring is modeled. We examined the statistical properties of a previously proposed quantitative trait threshold (QTT) model developed for when censored data can be reasonably inferred to be beyond an unknown threshold. METHODS: The QTT model is a Bayesian model integration approach implemented in the PPL framework that requires neither specification of the threshold nor imputation of the missing data. This model was evaluated under a range of simulated data sets and compared to other methods with missing data imputed. RESULTS: Across the simulated conditions, the addition of a threshold parameter did not change the PPL's properties relative to quantitative trait analysis on non-censored data except for a slight reduction in the average PPL as a reflection of the lowered information content due to censoring. This remained the case for non-normally distributed data and extreme sampling of pedigrees. CONCLUSIONS: Overall, the QTT model showed the smallest loss of linkage information relative to alternative approaches and therefore provides a unique analysis tool that obviates the need for ad hoc imputation of censored data in gene mapping studies.

EARLY CONSIDERATIONS OF GENETICS IN APHASIA REHABILITATION: A NARRATIVE REVIEW.

Background: Early investigations linking language and genetics were focused on the evolution of human communication in populations with developmental speech and language disorders. Recently, studies suggest that genes may also modulate recovery from post-stroke aphasia. Aims: Our goal is to review current literature related to the influence of genetics on post-stroke recovery, and the implications for aphasia rehabilitation. We describe candidate genes implicated by empirical findings and address additional clinical considerations. Main Contribution: We describe existing evidence and mechanisms supporting future investigations into how genetic factors may modulate aphasia recovery and propose that two candidate genes, brain derived neurotrophic factor (BDNF) and apolipoprotein E (APOE), may be important considerations for future research assessing response to aphasia treatment. Evidence suggests that BDNF is important for learning, memory, and neuroplasticity. APOE influences cognitive functioning and memory in older individuals and has also been implicated in neural repair. Moreover, recent data suggest an interaction between specific alleles of the BDNF and APOE genes in influencing episodic memory. Conclusions: Genetic influences on recovery from aphasia have been largely unexplored in the literature despite evidence that genetic factors influence behaviour and recovery from brain injury. As researchers continue to explore prognostic factors that may influence response to aphasia treatment, it is time for genetic factors to be considered as a source of variability. As the field moves in the direction of personalized medicine, eventually allied health professionals may utilize genetic profiles to inform treatment decisions and education for patients and care partners.

Invasive or More Direct Measurements Can Provide an Objective Early-Stopping Ceiling for Training Deep Neural Networks on Non-invasive or Less-Direct Biomedical Data.

Early stopping is an extremely common tool to minimize overfitting, which would otherwise be a cause of poor generalization of the model to novel data. However, early stopping is a heuristic that, while effective, primarily relies on ad hoc parameters and metrics. Optimizing when to stop remains a challenge. In this paper, we suggest that for some biomedical applications, a natural dichotomy of invasive/non-invasive measurements, or more generally proximal vs distal measurements of a biological system can be exploited to provide objective advice on early stopping. We discuss the conditions where invasive measurements of a biological process should provide better predictions than non-invasive measurements, or at best offer parity. Hence, if data from an invasive measurement are available locally, or from the literature, that information can be leveraged to know with high certainty whether a model of non-invasive data is overfitted. We present paired invasive/non-invasive cardiac and coronary artery measurements from two mouse strains, one of which spontaneously develops type 2 diabetes, posed as a classification problem. Examination of the various stopping rules shows that generalization is reduced with more training epochs and commonly applied stopping rules give widely different generalization error estimates. The use of an empirically derived training ceiling is demonstrated to be helpful as added information to leverage early stopping in order to reduce overfitting.

The Rosetta Phenotype Harmonization Method Facilitates Finding a Relationship Quantitative Trait Locus for a Complex Cognitive Trait.

Genetics researchers increasingly combine data across many sources to increase power and to conduct analyses that cross multiple individual studies. However, there is often a lack of alignment on outcome measures when the same constructs are examined across studies. This inhibits comparison across individual studies and may impact the findings from meta-analysis. Using a well-characterized genotypic (brain-derived neurotrophic factor: BDNF) and phenotypic constructs (working memory and reading comprehension), we employ an approach called Rosetta, which allows for the simultaneous examination of primary studies that employ related but incompletely overlapping data. We examined four studies of BDNF, working memory, and reading comprehension with a combined sample size of 1711 participants. Although the correlation between working memory and reading comprehension over all participants was high, as expected (ρ = 0.45), the correlation between working memory and reading comprehension was attenuated in the BDNF Met/Met genotype group (ρ = 0.18, n.s.) but not in the Val/Val (ρ = 0.44) or Val/Met (ρ = 0.41) groups. These findings indicate that Met/Met carriers may be a unique and robustly defined subgroup in terms of memory and reading comprehension. This study demonstrates the utility of the Rosetta method when examining complex phenotypes across multiple studies, including psychiatric genetic studies, as shown here, and also for the mega-analysis of cohorts generally.

The MPRINT Hub Data, Model, Knowledge and Research Coordination Center: Bridging the gap in maternal-pediatric therapeutics research through data integration and pharmacometrics.

Maternal and pediatric populations have historically been considered "therapeutic orphans" due to their limited inclusion in clinical trials. Physiologic changes during pregnancy and lactation and growth and maturation of children alter pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. Precision therapy in these populations requires knowledge of these effects. Efforts to enhance maternal and pediatric participation in clinical studies have increased over the past few decades. However, studies supporting precision therapeutics in these populations are often small and, in isolation, may have limited impact. Integration of data from various studies, for example through physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling or bioinformatics approaches, can augment the value of data from these studies, and help identify gaps in understanding. To catalyze research in maternal and pediatric precision therapeutics, the Obstetric and Pediatric Pharmacology and Therapeutics Branch of the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) established the Maternal and Pediatric Precision in Therapeutics (MPRINT) Hub. Herein, we provide an overview of the status of maternal-pediatric therapeutics research and introduce the Indiana University-Ohio State University MPRINT Hub Data, Model, Knowledge and Research Coordination Center (DMKRCC), which aims to facilitate research in maternal and pediatric precision therapeutics through the integration and assessment of existing knowledge, supporting pharmacometrics and clinical trials design, development of new real-world evidence resources, educational initiatives, and building collaborations among public and private partners, including other NICHD-funded networks. By fostering use of existing data and resources, the DMKRCC will identify critical gaps in knowledge and support efforts to overcome these gaps to enhance maternal-pediatric precision therapeutics.

An eQTL biological data visualization challenge and approaches from the visualization community.

In 2011, the IEEE VisWeek conferences inaugurated a symposium on Biological Data Visualization. Like other domain-oriented Vis symposia, this symposium's purpose was to explore the unique characteristics and requirements of visualization within the domain, and to enhance both the Visualization and Bio/Life-Sciences communities by pushing Biological data sets and domain understanding into the Visualization community, and well-informed Visualization solutions back to the Biological community. Amongst several other activities, the BioVis symposium created a data analysis and visualization contest. Unlike many contests in other venues, where the purpose is primarily to allow entrants to demonstrate tour-de-force programming skills on sample problems with known solutions, the BioVis contest was intended to whet the participants' appetites for a tremendously challenging biological domain, and simultaneously produce viable tools for a biological grand challenge domain with no extant solutions. For this purpose expression Quantitative Trait Locus (eQTL) data analysis was selected. In the BioVis 2011 contest, we provided contestants with a synthetic eQTL data set containing real biological variation, as well as a spiked-in gene expression interaction network influenced by single nucleotide polymorphism (SNP) DNA variation and a hypothetical disease model. Contestants were asked to elucidate the pattern of SNPs and interactions that predicted an individual's disease state. 9 teams competed in the contest using a mixture of methods, some analytical and others through visual exploratory methods. Independent panels of visualization and biological experts judged entries. Awards were given for each panel's favorite entry, and an overall best entry agreed upon by both panels. Three special mention awards were given for particularly innovative and useful aspects of those entries. And further recognition was given to entries that correctly answered a bonus question about how a proposed "gene therapy" change to a SNP might change an individual's disease status, which served as a calibration for each approaches' applicability to a typical domain question. In the future, BioVis will continue the data analysis and visualization contest, maintaining the philosophy of providing new challenging questions in open-ended and dramatically underserved Bio/Life Sciences domains.

Improvement of automated analysis of coronary Doppler echocardiograms.

Coronary artery disease is the leading cause of heart disease, and while it can be assessed through transthoracic Doppler echocardiography (TTDE) by observing changes in coronary flow, manual analysis of TTDE is time consuming and subject to bias. In a previous study, a program was created to automatically analyze coronary flow patterns by parsing Doppler videos into a single continuous image, binarizing and separating the image into cardiac cycles, and extracting data values from each of these cycles. The program significantly reduced variability and time to complete TTDE analysis, but some obstacles such as interfering noise and varying video sizes left room to increase the program's accuracy. The goal of this current study was to refine the existing automation algorithm and heuristics by (1) moving the program to a Python environment, (2) increasing the program's ability to handle challenging cases and video variations, and (3) removing unrepresentative cardiac cycles from the final data set. With this improved analysis, examiners can use the automatic program to easily and accurately identify the early signs of serious heart diseases.

MicroRNA and MicroRNA-Target Variants Associated with Autism Spectrum Disorder and Related Disorders.

Autism spectrum disorder (ASD) is a childhood neurodevelopmental disorder with a complex and heterogeneous genetic etiology. MicroRNA (miRNA), a class of small non-coding RNAs, could regulate ASD risk genes post-transcriptionally and affect broad molecular pathways related to ASD and associated disorders. Using whole-genome sequencing, we analyzed 272 samples in 73 families in the New Jersey Language and Autism Genetics Study (NJLAGS) cohort. Families with at least one ASD patient were recruited and were further assessed for language impairment, reading impairment, and other associated phenotypes. A total of 5104 miRNA variants and 1,181,148 3' untranslated region (3' UTR) variants were identified in the dataset. After applying several filtering criteria, including population allele frequency, brain expression, miRNA functional regions, and inheritance patterns, we identified high-confidence variants in five brain-expressed miRNAs (targeting 326 genes) and 3' UTR miRNA target regions of 152 genes. Some genes, such as SCP2 and UCGC, were identified in multiple families. Using Gene Ontology overrepresentation analysis and protein-protein interaction network analysis, we identified clusters of genes and pathways that are important for neurodevelopment. The miRNAs and miRNA target genes identified in this study are potentially involved in neurodevelopmental disorders and should be considered for further functional studies.

Genetic covariation underlying reading, language and related measures in a sample selected for specific language impairment.

Specific language impairment is a developmental language disorder characterized by failure to develop language normally in the absence of a specific cause. Previous twin studies have documented the heritability of reading and language measures as well as the genetic correlation between those measures. This paper presents results from an alternative to the classical twin designs by estimating heritability from extended pedigrees. These pedigrees were previously studied as part of series of molecular genetic studies of specific language impairment where the strongest genetic findings were with reading phenotypes rather than language despite selecting pedigrees based on language impairments. To explore the relationship between reading and language in these pedigrees, variance components estimates of heritability of reading and language measures were conducted showing general agreement with the twin literature, as were genetics correlations between reading and language. Phonological short-term memory, phonological awareness and auditory processing were evaluated as candidate mediators of the reading-language genetic correlations. Only phonological awareness showed significant genetic correlations with all reading measures and several language measures while phonological short-term memory and auditory processing did not.

Increasing genotype-phenotype model determinism: application to bivariate reading/language traits and epistatic interactions in language-impaired families.

While advances in network and pathway analysis have flourished in the era of genome-wide association analysis, understanding the genetic mechanism of individual loci on phenotypes is still readily accomplished using genetic modeling approaches. Here, we demonstrate two novel genotype-phenotype models implemented in a flexible genetic modeling platform. The examples come from analysis of families with specific language impairment (SLI), a failure to develop normal language without explanatory factors such as low IQ or inadequate environment. In previous genome-wide studies, we observed strong evidence for linkage to 13q21 with a reading phenotype in language-impaired families. First, we elucidate the genetic architecture of reading impairment and quantitative language variation in our samples using a bivariate analysis of reading impairment in affected individuals jointly with language quantitative phenotypes in unaffected individuals. This analysis largely recapitulates the baseline analysis using the categorical trait data (posterior probability of linkage (PPL) = 80%), indicating that our reading impairment phenotype captured poor readers who also have low language ability. Second, we performed epistasis analysis using a functional coding variant in the brain-derived neurotrophic factor (BDNF) gene previously associated with reduced performance on working memory tasks. Modeling epistasis doubled the evidence on 13q21 and raised the PPL to 99.9%, indicating that BDNF and 13q21 susceptibility alleles are jointly part of the genetic architecture of SLI. These analyses provide possible mechanistic insights for further cognitive neuroscience studies based on the models developed herein.