Family-based methods for linkage and association analysis

Statistics to prioritize rare variants in family-based sequencing studies with disease subtypes

Family-based sequencing studies are increasingly used to find rare genetic variants of high risk for disease traits with familial clustering. In some studies, families with multiple disease subtypes are collected and the exomes of affected relatives are sequenced for shared rare variants (RVs). Since different families can harbor different causal variants and each family harbors many RVs, tests to detect causal variants can have low power in this study design. Our goal is rather to prioritize shared variants for further investigation by, for example, pathway analyses or functional studies. The transmission-disequilibrium test prioritizes variants based on departures from Mendelian transmission in parent-child trios. Extending this idea to families, we propose methods to prioritize RVs shared in affected relatives with two disease subtypes, with one subtype more heritable than the other. Global approaches condition on a variant being observed in the study and assume a known probability of carrying a causal variant. In contrast, local approaches condition on a variant being observed in specific families to eliminate the carrier probability. Our simulation results indicate that global approaches are robust to misspecification of the carrier probability and prioritize more effectively than local approaches even when the carrier probability is misspecified.

Replication of previous autism-GWAS hits suggests the association between NAA1, SORCS3, and GSDME and autism in the Han Chinese population

Background

Autism is a severe neurodevelopmental disorder characterized by social interaction deficits, impairments in communication, and restricted and repetitive stereotyped behavior and activities. Family and twin studies suggested an essential role of genetic factors in the etiology of autism spectrum disorder (ASD). Also, other studies found SORCS3 and GSDME (DFNA5) might be involved in brain development and susceptible to ASD.

Methods

In this study, 17 genome-wide significant SNPs reported in previous ASD genome-wide association studies (GWAS) and 7 SNPs in strong linkage disequilibrium with known ASD GWAS hits were selected to investigate the association between these SNPs and autism in the Han Chinese population. Then, 10 tagSNPs in SORCS3 and 11 tagSNPs in GSDME were selected to analyze the association between these SNPs and autism. The selected 24 SNPs and tagSNPs were genotyped using the Agena MassARRAY SNP genotyping assay in 757 Han Chinese autism trios.

Results

Rs1484144 in NAA11 was significantly associated with autism; significance remained after the Bonferroni correction (P < 0.0022). Also, rs79879286, rs12154597, and rs12540919 near GSDME, as well as rs9787523 and rs3750261 in SORCS3, were nominally associated with autism.

Conclusion

Our study suggests that rs1484144 in NAA11 is a significant SNP for autism in the Han Chinese population, while SORCS3 and GSDME might be the susceptibility genes for autism in this population.

A pedigree-based cohort to study the genetic risk factors for cardiometabolic diseases: study design, baseline characteristics and preliminary results

Background

We initiated the Fujian Tulou Pedigree-based Cohort (FTPC) as the integration of extended pedigrees and prospective cohort to clarify the genetic and environmental risk factors of cardiometabolic diseases.

Methods

FTPC was carried out in Nanjing County, Fujian Province, China from August 2015 to December 2017 to recruit probands with the same surnames and then enroll their first-degree and more distant relatives. The participants were asked to complete questionnaire interview, physical examination, and blood collection. According to the local genealogical booklets and family registry, we reconstructed extended pedigrees to estimate the heritability of cardiometabolic traits. The follow-up of FTPC is scheduled every 5 years in the future.

Results

The baseline survey interviewed 2,727 individuals in two clans. A total of 1,563 adult subjects who completed all baseline examinations were used to reconstruct pedigrees and 452 extended pedigrees were finally identified, including one seven-generation pedigree, two five-generation pedigrees, 23 four-generation pedigrees, 186 three-generation pedigrees, and 240 two-generation pedigrees. The average age of the participants was 57.4 years, with 43.6% being males. The prevalence of hypertension, diabetes and dyslipidemia in FTPC were 49.2, 10.0, and 45.2%, respectively. Based on the pedigree structure, the heritability of systolic blood pressure, diastolic blood pressure, fast blood glucose, total cholesterol, triglyceride, high-density lipoprotein, and low-density lipoprotein was estimated at 0.379, 0.306, 0.386, 0.452, 0.568, 0.852, and 0.387, respectively.

Conclusion

As an extended pedigree cohort in China, FTPC will provide an important source to study both genetic and environmental risk factors prospectively.

Benchmarking statistical methods for analyzing parent-child dyads in genetic association studies

Genetic association studies of child health outcomes often employ family-based study designs. One of the most popular family-based designs is the case-parent trio design that considers the smallest possible nuclear family consisting of two parents and their affected child. This trio design is particularly advantageous for studying relatively rare disorders because it is less prone to type 1 error inflation due to population stratification compared to population-based study designs (e.g., case-control studies). However, obtaining genetic data from both parents is difficult, from a practical perspective, and many large studies predominantly measure genetic variants in mother-child dyads. While some statistical methods for analyzing parent-child dyad data (most commonly involving mother-child pairs) exist, it is not clear if they provide the same advantage as trio methods in protecting against population stratification, or if a specific dyad design (e.g., case-mother dyads vs. case-mother/control-mother dyads) is more advantageous. In this article, we review existing statistical methods for analyzing genome-wide marker data on dyads and perform extensive simulation experiments to benchmark their type I errors and statistical power under different scenarios. We extend our evaluation to existing methods for analyzing a combination of case-parent trios and dyads together. We apply these methods on genotyped and imputed data from multiethnic mother-child pairs only, case-parent trios only or combinations of both dyads and trios from the Gene, Environment Association Studies consortium (GENEVA), where each family was ascertained through a child affected by nonsyndromic cleft lip with or without cleft palate. Results from the GENEVA study corroborate the findings from our simulation experiments. Finally, we provide recommendations for using statistical genetic association methods for dyads.

A Novel Gene CDC27 Causes SLE and Is Associated With the Disease Activity

Background

As genetic genetic factors are important in SLE, so screening causative genes is of great significance for the prediction and early prevention in people who may develop SLE. At present, it is very difficult to screen causative genes through pedigrees. The analytical method described herein can be used to screen causative genes for SLE and other complex diseases through pedigrees.

Methods

For the first time, 24 lupus pedigrees were analyzed by combining whole exon sequencing and a variety of biological information tools including common-specific analysis, pVAAST (pedigree variant annotation, analysis and search tool), Exomiser (Combining phenotype and PPI associated analysis), and FARVAT (family based gene burden), and the causative genes of these families with lupus identified. Selected causative genes in peripheral-blood mononuclear cells (PBMCs) were evaluated by quantitative polymerase chain reaction (qPCR).

Results

Cell division cycle 27 (CDC27) was screened out by common-specific analysis and Exomiser causative gene screening. FARVAT analysis on these families detected only CDC27 at the extremely significant level (false discovery rate <0.05) by three family-based burden analyses (BURDEN, CALPHA, and SKATO). QPCR was performed to detect for CDC27 in the PBMCs of the SLE family patients, sporadic lupus patients, and healthy people. Compared with the healthy control group, CDC27 expression was low in lupus patients (familial and sporadic patients) (P<0.05) and correlated with lupus activity indicators: negatively with C-reactive protein (CRP) (P<0.05) and erythrocyte sedimentation rate (P<0.05) and positively with complement C3 and C4 (P<0.05). The CDC27 expression was upregulated in PBMCs from SLE patients with reduced lupus activity after immunotherapy (P<0.05). Based on Receiver operating characteristic (ROC) curve analysis, the sensitivity and specificity of CDC27 in diagnosing SLE were 82.30% and 94.40%.

Conclusion

The CDC27 gene, as found through WES combined with multiple analytical method may be a causative gene of lupus. CDC27 may serve as a marker for the diagnosis of SLE and is closely related to the lupus activity. We hope that the analytical method in this study will be used to screen causative genes for other diseases through small pedigrees, especially among non-close relatives.

Identification and selection of healthy spermatozoa in heterozygous carriers of the Phe508del-variant of the CFTR-gene in assisted reproduction

The pathogenic variant Phe508del of the CFTR-gene is the most frequent cause of cystic fibrosis (CF). Whereas male CF-patients are infertile due to bilateral agenesis of the efferent ducts, the fertility status of male heterozygous carriers is uncertain. We aimed at demonstrating the involvement of the CFTR-ion channel during sperm capacitation and to potentially select variant-free spermatozoa in heterozygous carriers of the CFTR-variant using flow cytometry (FC). Using FC and sorting, single cell polymerase chain reaction, immuno-fluorescent staining an experimental study was performed on nine fertile semen donors and three heterozygous infertile men carrying the Phe508del gene variant. Chemical inhibition of CFTR interfered with sperm capacitation. Most viable spermatozoa of heterozygous carriers of the Phe508del variant of the CFTR-gene show immune-fluorescent CFTR. Sperm capacitation in Phe508del carriers was similar to that in healthy semen donors. Distribution of the Phe508del allele in trio data of CF-affected families corresponded to the expected recessive inheritance pattern. Infertility in Phe508del heterozygous men is unlikely to be caused by the pathogenic variant although some contribution cannot be excluded. Normal sperm capacitation in carriers of pathogenic variants of the Phe508del-gene may in part explain the high prevalence of a potentially lethal recessive disorder.

The Complex and Diverse Genetic Architecture of Dilated Cardiomyopathy

Our insight into the diverse and complex nature of dilated cardiomyopathy (DCM) genetic architecture continues to evolve rapidly. The foundations of DCM genetics rest on marked locus and allelic heterogeneity. While DCM exhibits a Mendelian, monogenic architecture in some families, preliminary data from our studies and others suggests that at least 20% to 30% of DCM may have an oligogenic basis, meaning that multiple rare variants from different, unlinked loci, determine the DCM phenotype. It is also likely that low-frequency and common genetic variation contribute to DCM complexity, but neither has been examined within a rare variant context. Other types of genetic variation are also likely relevant for DCM, along with gene-by-environment interaction, now established for alcohol- and chemotherapy-related DCM. Collectively, this suggests that the genetic architecture of DCM is broader in scope and more complex than previously understood. All of this elevates the impact of DCM genetics research, as greater insight into the causes of DCM can lead to interventions to mitigate or even prevent it and thus avoid the morbid and mortal scourge of human heart failure.

Case-Parent Trio Studies in Cleft Lip and Palate

Cleft lip with or without cleft palate (CL/P) is one of the most common congenital malformations in humans involving various genetic and environmental risk factors. The prevalence of CL/P varies according to geographical location, ethnicity, race, gender, and socioeconomic status, affecting approximately 1 in 800 live births worldwide. Genetic studies aim to understand the mechanisms contributory to a phenotype by measuring the association between genetic variants and also between genetic variants and phenotype population. Genome-wide association studies are standard tools used to discover genetic loci related to a trait of interest. Genetic association studies are generally divided into two main design types: population-based studies and family-based studies. The epidemiological population-based studies comprise unrelated individuals that directly compare the frequency of genetic variants between (usually independent) cases and controls. The alternative to population-based studies (case–control designs) includes various family-based study designs that comprise related individuals. An example of such a study is a case–parent trio design study, which is commonly employed in genetics to identify the variants underlying complex human disease where transmission of alleles from parents to offspring is studied. This article describes the fundamentals of case–parent trio study, trio design and its significances, statistical methods, and limitations of the trio studies.

Huntingtin gene CAG repeat size affects autism risk: Family-based and case-control association study

The Huntingtin (HTT) gene contains a CAG repeat in exon 1, whose expansion beyond 39 repeats consistently leads to Huntington's disease (HD), whereas normal-to-intermediate alleles seemingly modulate brain structure, function and behavior. The role of the CAG repeat in Autism Spectrum Disorder (ASD) was investigated applying both family-based and case-control association designs, with the SCA3 repeat as a negative control. Significant overtransmission of "long" CAG alleles (≥17 repeats) to autistic children and of "short" alleles (≤16 repeats) to their unaffected siblings (all p < 10^-5 ) was observed in 612 ASD families (548 simplex and 64 multiplex). Surprisingly, both 193 population controls and 1,188 neurological non-HD controls have significantly lower frequencies of "short" CAG alleles compared to 185 unaffected siblings and higher rates of "long" alleles compared to 548 ASD patients from the same families (p < .05-.001). The SCA3 CAG repeat displays no association. "Short" HTT alleles seemingly exert a protective effect from clinically overt autism in families carrying a genetic predisposition for ASD, while "long" alleles may enhance autism risk. Differential penetrance of autism-inducing genetic/epigenetic variants may imply atypical developmental trajectories linked to HTT functions, including excitation/inhibition imbalance, cortical neurogenesis and apoptosis, neuronal migration, synapse formation, connectivity and homeostasis.

Genetic variants affecting bone mineral density and bone mineral content at multiple skeletal sites in Hispanic children

CONTEXT

Osteoporosis is a major public health burden with significant economic costs. However, the correlates of bone health in Hispanic children are understudied.

OBJECTIVE

We aimed to identify genetic variants associated with bone mineral density (BMD) and bone mineral content (BMC) at multiple skeletal sites in Hispanic children.

METHODS

We conducted a cross-sectional genome-wide linkage analysis, genome-wide and exome-wide association analysis of BMD and BMC. The Viva La Familia Study is a family-based cohort with a total of 1030 Hispanic children (4-19 years old at baseline) conducted in Houston, TX. BMD and BMC were measured by Dual-energy X-ray absorptiometry.

RESULTS

Significant heritability were observed for BMC and BMD at multiple skeletal sites ranging between 44 and 68% (P < 2.8 × 10-9). Significant evidence for linkage was found for BMD of pelvis and left leg on chromosome 7p14, lumbar spine on 20q13 and left rib on 6p21, and BMC of pelvis on chromosome 20q12 and total body on 14q22-23 (logarithm of odds score > 3). We found genome-wide significant association between BMC of right arm and rs762920 at PVALB (P = 4.6 × 10-8), and between pelvis BMD and rs7000615 at PTK2B (P = 7.4 × 10-8). Exome-wide association analysis revealed novel association of variants at MEGF10 and ABRAXAS2 with left arm and lumber spine BMC, respectively (P < 9 × 10-7).

CONCLUSIONS

We identified novel loci associated with BMC and BMD in Hispanic children, with strongest evidence for PTK2B. These findings provide better understanding of bone genetics and shed light on biological mechanisms underlying BMD and BMC variation.

Whole genome sequencing of orofacial cleft trios from the Gabriella Miller Kids First Pediatric Research Consortium identifies a new locus on chromosome 21

Orofacial clefts (OFCs) are among the most prevalent craniofacial birth defects worldwide and create a significant public health burden. The majority of OFCs are non-syndromic, and the genetic etiology of non-syndromic OFCs is only partially determined. Here, we analyze whole genome sequence (WGS) data for association with risk of OFCs in European and Colombian families selected from a multicenter family-based OFC study. This is the first large-scale WGS study of OFC in parent-offspring trios, and a part of the Gabriella Miller Kids First Pediatric Research Program created for the study of childhood cancers and structural birth defects. WGS provides deeper and more specific genetic data than using imputation on present-day single nucleotide polymorphic (SNP) marker panels. Genotypes of case-parent trios at single nucleotide variants (SNV) and short insertions and deletions (indels) spanning the entire genome were called from their sequences using human GRCh38 genome assembly, and analyzed for association using the transmission disequilibrium test. Among genome-wide significant associations, we identified a new locus on chromosome 21 in Colombian families, not previously observed in other larger OFC samples of Latin American ancestry. This locus is situated within a region known to be expressed during craniofacial development. Based on deeper investigation of this locus, we concluded that it contributed risk for OFCs exclusively in the Colombians. This study reinforces the ancestry differences seen in the genetic etiology of OFCs, and underscores the need for larger samples when studying for OFCs and other birth defects in populations with diverse ancestry.

PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate

Nonsyndromic cleft lip and palate (NSCLP) is one of the most common craniofacial anomalies in humans, affecting more than 135,000 newborns worldwide. NSCLP has a multifactorial etiology with more than 50 genes postulated to play an etiologic role. The genetic pathway comprised of Pbx-Wnt-p63-Irf6 genes was shown to control facial morphogenesis in mice and proposed as a regulatory pathway for NSCLP. Based on these findings, we investigated whether variation in PBX1, PBX2, and TP63, and their proposed interactions were associated with NSCLP. Fourteen single nucleotide variants (SNVs) in/nearby PBX1, PBX2, and TP63 were genotyped in 780 NSCLP families of nonHispanic white (NHW) and Hispanic ethnicities. Family-based association tests were performed for individual SNVs stratified by ethnicity and family history of NSCLP. Gene-gene interactions were also tested. A significant association was found for PBX2 rs3131300 and NSCLP in combined Hispanic families (p = .003) while nominal association was found for TP63 rs9332461 in multiplex Hispanic families (p = .005). Significant haplotype associations were observed for PBX2 in NHW (p = .0002) and Hispanic families (p = .003), and for TP63 in multiplex Hispanic families (.003). An independent case-control group was used to validate findings, and significant associations were found with PBX1 rs6426870 (p = .007) and TP63 rs9332461 (p = .03). Gene-gene interactions were detected between PBX1/PBX2/TP63 with IRF6 in NHW families, and between PBX1 with WNT9B in both NHW and Hispanic families (p < .0018). This study provides the first evidence for a role of PBX1 and PBX2, additional evidence for the role of TP63, and support for the proposed PBX-WNT-TP63-IRF6 regulatory pathway in the etiology of NSCLP.

Statistical Methods and Software for Substance Use and Dependence Genetic Research

BACKGROUND

Substantial substance use disorders and related health conditions emerged dur-ing the mid-20th century and continue to represent a remarkable 21st century global burden of disease. This burden is largely driven by the substance-dependence process, which is a complex process and is influenced by both genetic and environmental factors. During the past few decades, a great deal of pro-gress has been made in identifying genetic variants associated with Substance Use and Dependence (SUD) through linkage, candidate gene association, genome-wide association and sequencing studies.

METHODS

Various statistical methods and software have been employed in different types of SUD ge-netic studies, facilitating the identification of new SUD-related variants.

CONCLUSION

In this article, we review statistical methods and software that are currently available for SUD genetic studies, and discuss their strengths and limitations.

Genetic loci associated with an earlier age at onset in multiplex schizophrenia

An earlier age at onset (AAO) has been associated with greater genetic loadings in schizophrenia. This study aimed to identify modifier loci associated with an earlier AAO of schizophrenia. A genome-wide association analysis (GWAS) was conducted in 94 schizophrenia probands with the earliest AAO and 91 with the latest AAO. Candidate single nucleotide polymorphisms (SNPs) were then genotyped in the co-affected siblings and unrelated probands. Multi-SNP genetic risk scores (GRS) composed of the candidate loci were used to distinguish patients with an early or late AAO. The 14-SNP GRS could distinguish the co-affected siblings (n = 90) of the earliest probands from those (n = 91) of the latest probands. When 132 patients with an earlier AAO and 158 patients with a later AAO were included, a significant trend in the 14-SNP GRS was detected among those unrelated probands from 4 family groups with the earliest, earlier, later, and latest AAO. The overall effect of the 14 SNPs on an AAO in schizophrenia was verified using co-affected siblings of the GWAS probands and trend effect across unrelated patients. Preliminary network analysis of these loci revealed the involvement of PARK2, a gene intensively reported in Parkinson’s disease and schizophrenia research.

Birth weight interacts with a functional variant of the oxytocin receptor gene (OXTR) to predict executive functioning in children

Genetic variation in the oxytocin receptor gene (OXTR) is associated with several psychiatric conditions characterized by deficits in executive functioning (EF). A specific OXTR variant, rs2254298, has previously been associated with brain functioning in regions implicated in EF. Moreover, birth weight variation across the entire range is associated with individual differences in cortical structure and function that underlie EF. This is the first study to examine the main and interactive effect between rs2254298 and birth weight on EF in children. The sample consisted of 310 children from an ongoing longitudinal study. EF was measured at age 4.5 using observational tasks indexing working memory, cognitive flexibility, and inhibitory control. A family-based design that controlled for population admixture, stratification, and nongenomic confounds was employed. A significant genetic association between rs2254298 and EF was observed, with more copies of the major allele (G) associated with higher EF. There was also a significant interaction between rs2254298 and birth weight, such that more copies of the major allele in combination with higher birth weight predicted better EF. Findings suggest that OXTR may be associated with discrete neurocognitive abilities in childhood, and these effects may be modulated by intrauterine conditions related to fetal growth and development.

Genetic factors influencing risk to orofacial clefts: today's challenges and tomorrow's opportunities

Orofacial clefts include cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP), which combined represent the largest group of craniofacial malformations in humans with an overall prevalence of one per 1,000 live births. Each of these birth defects shows strong familial aggregation, suggesting a major genetic component to their etiology. Genetic studies of orofacial clefts extend back centuries, but it has proven difficult to define any single etiologic mechanism because many genes appear to influence risk. Both linkage and association studies have identified several genes influencing risk, but these differ across families and across populations. Genome-wide association studies have identified almost two dozen different genes achieving genome-wide significance, and there are broad classes of 'causal genes' for orofacial clefts: a few genes strongly associated with risk and possibly directly responsible for Mendelian syndromes which include orofacial clefts as a key phenotypic feature of the syndrome, and multiple genes with modest individual effects on risk but capable of disrupting normal craniofacial development under the right circumstances (which may include exposure to environmental risk factors). Genomic sequencing studies are now underway which will no doubt reveal additional genes/regions where variants (sequence and structural) can play a role in controlling risk to orofacial clefts. The real challenge to medicine and public health is twofold: to identify specific genes and other etiologic factors in families with affected members and then to devise effective interventions for these different biological mechanisms controlling risk to complex and heterogeneous birth defects such as orofacial clefts.

Family-based association analyses of imputed genotypes reveal genome-wide significant association of Alzheimer's disease with OSBPL6, PTPRG, and PDCL3

The genetic basis of Alzheimer's disease (AD) is complex and heterogeneous. Over 200 highly penetrant pathogenic variants in the genes APP, PSEN1, and PSEN2 cause a subset of early-onset familial AD. On the other hand, susceptibility to late-onset forms of AD (LOAD) is indisputably associated to the ɛ4 allele in the gene APOE, and more recently to variants in more than two-dozen additional genes identified in the large-scale genome-wide association studies (GWAS) and meta-analyses reports. Taken together however, although the heritability in AD is estimated to be as high as 80%, a large proportion of the underlying genetic factors still remain to be elucidated. In this study, we performed a systematic family-based genome-wide association and meta-analysis on close to 15 million imputed variants from three large collections of AD families (~3500 subjects from 1070 families). Using a multivariate phenotype combining affection status and onset age, meta-analysis of the association results revealed three single nucleotide polymorphisms (SNPs) that achieved genome-wide significance for association with AD risk: rs7609954 in the gene PTPRG (P-value=3.98 × 10^-8), rs1347297 in the gene OSBPL6 (P-value=4.53 × 10^-8), and rs1513625 near PDCL3 (P-value=4.28 × 10^-8). In addition, rs72953347 in OSBPL6 (P-value=6.36 × 10^-7) and two SNPs in the gene CDKAL1 showed marginally significant association with LOAD (rs10456232, P-value=4.76 × 10^-7; rs62400067, P-value=3.54 × 10^-7). In summary, family-based GWAS meta-analysis of imputed SNPs revealed novel genomic variants in (or near) PTPRG, OSBPL6, and PDCL3 that influence risk for AD with genome-wide significance.

Replication of genome-wide association study (GWAS) susceptibility loci in a Latino bipolar disorder cohort

OBJECTIVES

Recent genome-wide association studies (GWASs) have identified numerous putative genetic polymorphisms associated with bipolar disorder (BD) and/or schizophrenia (SC). We hypothesized that a portion of these polymorphisms would also be associated with BD in the Latino American population. To identify such regions, we tested previously identified genetic variants associated with BD and/or SC and ancestral haploblocks containing these single nucleotide polymorphisms (SNPs) in a sample of Latino subjects with BD.

METHODS

A total of 2254 Latino individuals were genotyped for 91 SNPs identified in previous BD and/or SC GWASs, along with selected SNPs in strong linkage disequilibrium with these markers. Family-based single marker and haplotype association testing was performed using the PBAT software package. Empirical P-values were derived from 10 000 permutations.

RESULTS

Associations of eight a priori GWAS SNPs with BD were replicated with nominal (P≤.05) levels of significance. These included SNPs within nuclear factor I A (NFIA), serologically defined colon cancer antigen 8 (SDCCAG8), lysosomal associated membrane protein 3 (LAMP3), nuclear factor kappa B subunit 1 (NFKB1), major histocompatibility complex, class I, B (HLA-B) and 5'-nucleotidase, cytosolic II (NT5C2) and SNPs within intragenic regions microRNA 6828 (MIR6828)-solute carrier family 7 member 14 (SLC7A14) and sonic hedgehog (SHH)-long intergenic non-protein coding RNA 1006 (LINC01006). Of the 76 ancestral haploblocks that were tested for associations with BD, our top associated haploblock was located in LAMP3; however, the association did not meet statistical thresholds of significance following Bonferroni correction.

CONCLUSIONS

These results indicate that some of the gene variants found to be associated with BD or SC in other populations are also associated with BD risk in Latinos. Variants in six genes and two intragenic regions were associated with BD in our Latino sample and provide additional evidence for overlap in genetic risk between SC and BD.

Oxytocin and vasopressin hormone genes in children's externalizing problems: A cognitive endophenotype approach

Externalizing problems are among the most common mental health problems of children. Research suggests that these problems are heritable, yet little is known about the specific genes involved in their pathophysiology. The current study examined a genotype-endophenotype-phenotype model of externalizing problems in 320 preschool-aged children. Markers of the oxytocin (OXT) and arginine vasopressin (AVP) hormone genes were selected as candidates owing to their known association with psychopathology in other domains. We tested whether OXT and AVP variants were related to children's externalizing problems, as well as two cognitive endophenotypes presumed to underlie these problems: theory of mind (ToM) and executive functioning (EF). Externalizing problems were assessed at age 4.5 using a previously-validated rating scale. ToM and EF were measured with age-appropriate tasks. Using a family-based association design and controlling for non-genomic confounds, support was found for an association between a two-marker OXT haplotype (rs2740210-rs2770378) and a two-marker AVP haplotype (rs1887854-rs3761249) and externalizing problems. Specific associations of these haplotypes with ToM and EF were also observed. Further, ToM and EF were shown to independently and jointly predict externalizing problems, and to partially mediate the effects of OXT and AVP on externalizing problems. This study provides the first evidence that genetic variation in OXT and AVP may contribute to individual differences in childhood externalizing problems, and that these effects may operate through emerging neurocognitive abilities in the preschool period.

The Norwegian preeclampsia family cohort study: a new resource for investigating genetic aspects and heritability of preeclampsia and related phenotypes

Background

Preeclampsia is a major pregnancy complication without curative treatment available. A Norwegian Preeclampsia Family Cohort was established to provide a new resource for genetic and molecular studies aiming to improve the understanding of the complex pathophysiology of preeclampsia.

Methods

Participants were recruited from five Norwegian hospitals after diagnoses of preeclampsia registered in the Medical birth registry of Norway were verified according to the study’s inclusion criteria. Detailed obstetric information and information on personal and family disease history focusing on cardiovascular health was collected. At attendance anthropometric measurements were registered and blood samples were drawn. The software package SPSS 19.0 for Windows was used to compute descriptive statistics such as mean and SD. P-values were computed based on t-test statistics for normally distributed variables. Nonparametrical methods (chi square) were used for categorical variables.

Results

A cohort consisting of 496 participants (355 females and 141 males) representing 137 families with increased occurrence of preeclampsia has been established, and blood samples are available for 477 participants. Descriptive analyses showed that about 60 % of the index women’s pregnancies with birth data registered were preeclamptic according to modern diagnosis criteria. We also found that about 41 % of the index women experienced more than one preeclamptic pregnancy. In addition, the descriptive analyses confirmed that preeclamptic pregnancies are more often accompanied with delivery complications.

Conclusion

The data and biological samples collected in this Norwegian Preeclampsia Family Cohort will provide an important basis for future research. Identification of preeclampsia susceptibility genes and new biomarkers may contribute to more efficient strategies to identify mothers “at risk” and contribute to development of novel preventative therapies.

Gene-environment interaction between the oxytocin receptor (OXTR) gene and parenting behaviour on children's theory of mind

Theory of mind (ToM) is the ability to interpret and understand human behaviour by representing the mental states of others. Like many human capacities, ToM is thought to develop through both complex biological and socialization mechanisms. However, no study has examined the joint effect of genetic and environmental influences on ToM. This study examined how variability in the oxytocin receptor gene (OXTR) and parenting behavior--two widely studied factors in ToM development-interacted to predict ToM in pre-school-aged children. Participants were 301 children who were part of an ongoing longitudinal birth cohort study. ToM was assessed at age 4.5 using a previously validated scale. Parenting was assessed through observations of mothers' cognitively sensitive behaviours. Using a family-based association design, it was suggestive that a particular variant (rs11131149) interacted with maternal cognitive sensitivity on children's ToM (P = 0.019). More copies of the major allele were associated with higher ToM as a function of increasing cognitive sensitivity. A sizeable 26% of the variability in ToM was accounted for by this interaction. This study provides the first empirical evidence of gene-environment interactions on ToM, supporting the notion that genetic factors may be modulated by potent environmental influences early in development.

Genetic linkage analysis in the age of whole-genome sequencing

For many years, linkage analysis was the primary tool used for the genetic mapping of Mendelian and complex traits with familial aggregation. Linkage analysis was largely supplanted by the wide adoption of genome-wide association studies (GWASs). However, with the recent increased use of whole-genome sequencing (WGS), linkage analysis is again emerging as an important and powerful analysis method for the identification of genes involved in disease aetiology, often in conjunction with WGS filtering approaches. Here, we review the principles of linkage analysis and provide practical guidelines for carrying out linkage studies using WGS data.

Methods for association analysis and meta-analysis of rare variants in families

Advances in exome sequencing and the development of exome genotyping arrays are enabling explorations of association between rare coding variants and complex traits. To ensure power for these rare variant analyses, a variety of association tests that group variants by gene or functional unit have been proposed. Here, we extend these tests to family-based studies. We develop family-based burden tests, variable frequency threshold tests and sequence kernel association tests. Through simulations, we compare the performance of different tests. We describe situations where family-based studies provide greater power than studies of unrelated individuals to detect rare variants associated with moderate to large changes in trait values. Broadly speaking, we find that when sample sizes are limited and only a modest fraction of all trait-associated variants can be identified, family samples are more powerful. Finally, we illustrate our approach by analyzing the relationship between coding variants and levels of high-density lipoprotein (HDL) cholesterol in 11,556 individuals from the HUNT and SardiNIA studies, demonstrating association for coding variants in the APOC3, CETP, LIPC, LIPG, and LPL genes and illustrating the value of family samples, meta-analysis, and gene-level tests. Our methods are implemented in freely available C++ code.

Genetics of kidney disease and related cardiometabolic phenotypes in Zuni Indians: the Zuni Kidney Project

The objective of this study is to identify genetic factors associated with chronic kidney disease (CKD) and related cardiometabolic phenotypes among participants of the Genetics of Kidney Disease in Zuni Indians study. The study was conducted as a community-based participatory research project in the Zuni Indians, a small endogamous tribe in rural New Mexico. We recruited 998 members from 28 extended multigenerational families, ascertained through probands with CKD who had at least one sibling with CKD. We used the Illumina Infinium Human1M-Duo version 3.0 BeadChips to type 1.1 million single nucleotide polymorphisms (SNPs). Prevalence estimates for CKD, hyperuricemia, diabetes, and hypertension were 24%, 30%, 17% and 34%, respectively. We found a significant (p < 1.58 × 10^-7) association for a SNP in a novel gene for serum creatinine (PTPLAD2). We replicated significant associations for genes with serum uric acid (SLC2A9), triglyceride levels (APOA1, BUD13, ZNF259), and total cholesterol (PVRL2). We found novel suggestive associations (p < 1.58 × 10^-6) for SNPs in genes with systolic (OLFML2B), and diastolic blood pressure (NFIA). We identified a series of genes associated with CKD and related cardiometabolic phenotypes among Zuni Indians, a population with a high prevalence of kidney disease. Illuminating genetic variations that modulate the risk for these disorders may ultimately provide a basis for novel preventive strategies and therapeutic interventions.

Identification of rare variants for hypertension with incorporation of linkage information

We conducted linkage analysis using the genome-wide association study data on chromosome 3, and then assessed association between hypertension and rare variants of genes located in the regions showing evidence of linkage. The rare variants were collapsed if their minor allele frequencies were less than or equal to the thresholds: 0.01, 0.03, or 0.05. In the collapsing process, they were either unweighted or weighted by the nonparametric linkage log of odds scores in 2 different schemes: exponential weighting and cumulative weighting. Logistic regression models using the generalized estimating equations approach were used to assess association between the collapsed rare variants and hypertension adjusting for age and gender. Evidence of association from the weighted and unweighted collapsing schemes with minor allele frequencies ≤0.01, after accounting for multiple testing, was found for genes DOCK3 (p = 0.0090), ARMC8 (p = 1.29E-5), KCNAB1 (p = 5.8E-4), and MYRIP (p = 5.79E-6). DOCK3 and MYRIP are newly discovered. Incorporating linkage scores as weights was found to help identify rare causal variants with a large effect size.

Lep-MAP: fast and accurate linkage map construction for large SNP datasets

MOTIVATION

Current high-throughput sequencing technologies allow cost-efficient genotyping of millions of single nucleotide polymorphisms (SNPs) for hundreds of samples. However, the tools that are currently available for constructing linkage maps are not well suited for large datasets. Linkage maps of large datasets would be helpful in de novo genome assembly by facilitating comprehensive genome validation and refinement by enabling chimeric scaffold detection, as well as in family-based linkage and association studies, quantitative trait locus mapping, analysis of genome synteny and other complex genomic data analyses.

RESULTS

We describe a novel tool, called Lepidoptera-MAP (Lep-MAP), for constructing accurate linkage maps with ultradense genome-wide SNP data. Lep-MAP is fast and memory efficient and largely automated, requiring minimal user interaction. It uses simultaneously data on multiple outbred families and can increase linkage map accuracy by taking into account achiasmatic meiosis, a special feature of Lepidoptera and some other taxa with no recombination in one sex (no recombination in females in Lepidoptera). We demonstrate that Lep-MAP outperforms other methods on real and simulated data. We construct a genome-wide linkage map of the Glanville fritillary butterfly (Melitaea cinxia) with over 40 000 SNPs. The data were generated with a novel in-house SOLiD restriction site-associated DNA tag sequencing protocol, which is described in the online supplementary material.

AVAILABILITY AND IMPLEMENTATION

Java source code under GNU general public license with the compiled classes and the datasets are available from http://sourceforge.net/users/lep-map.

Dosage transmission disequilibrium test (dTDT) for linkage and association detection

Both linkage and association studies have been successfully applied to identify disease susceptibility genes with genetic markers such as microsatellites and Single Nucleotide Polymorphisms (SNPs). As one of the traditional family-based studies, the Transmission/Disequilibrium Test (TDT) measures the over-transmission of an allele in a trio from its heterozygous parents to the affected offspring and can be potentially useful to identify genetic determinants for complex disorders. However, there is reduced information when complete trio information is unavailable. In this study, we developed a novel approach to "infer" the transmission of SNPs by combining both the linkage and association data, which uses microsatellite markers from families informative for linkage together with SNP markers from the offspring who are genotyped for both linkage and a Genome-Wide Association Study (GWAS). We generalized the traditional TDT to process these inferred dosage probabilities, which we name as the dosage-TDT (dTDT). For evaluation purpose, we developed a simulation procedure to assess its operating characteristics. We applied the dTDT to the simulated data and documented the power of the dTDT under a number of different realistic scenarios. Finally, we applied our methods to a family study of alcohol dependence (COGA) and performed individual genotyping on complete families for the top signals. One SNP (rs4903712 on chromosome 14) remained significant after correcting for multiple testing Methods developed in this study can be adapted to other platforms and will have widespread applicability in genomic research when case-control GWAS data are collected in families with existing linkage data.

Efficient simulation of epistatic interactions in case-parent trios

Statistical approaches to evaluate interactions between single nucleotide polymorphisms (SNPs) and SNP-environment interactions are of great importance in genetic association studies, as susceptibility to complex disease might be related to the interaction of multiple SNPs and/or environmental factors. With these methods under active development, algorithms to simulate genomic data sets are needed to ensure proper type I error control of newly proposed methods and to compare power with existing methods. In this paper we propose an efficient method for a haplotype-based simulation of case-parent trios when the disease risk is thought to depend on possibly higher-order epistatic interactions or gene-environment interactions with binary exposures.

Suggestive evidence for association between L-type voltage-gated calcium channel (CACNA1C) gene haplotypes and bipolar disorder in Latinos: a family-based association study

OBJECTIVES

  Through recent genome-wide association studies (GWASs), several groups have reported significant association between variants in the calcium channel, voltage-dependent, L-type, alpha 1C subunit (CACNA1C) and bipolar disorder (BP) in European and European-American cohorts. We performed a family-based association study to determine whether CACNA1C is associated with BP in the Latino population.

METHODS

  This study included 913 individuals from 215 Latino pedigrees recruited from the USA, Mexico, Guatemala, and Costa Rica. The Illumina GoldenGate Genotyping Assay was used to genotype 58 single-nucleotide polymorphisms (SNPs) that spanned a 602.9-kb region encompassing the CACNA1C gene including two SNPs (rs7297582 and rs1006737) previously shown to associate with BP. Individual SNP and haplotype association analyses were performed using Family-Based Association Test (version 2.0.3) and Haploview (version 4.2) software.

RESULTS

  An eight-locus haplotype block that included these two markers showed significant association with BP (global marker permuted p = 0.0018) in the Latino population. For individual SNPs, this sample had insufficient power (10%) to detect associations with SNPs with minor effect (odds ratio = 1.15).

CONCLUSIONS

  Although we were not able to replicate findings of association between individual CACNA1C SNPs rs7297582 and rs1006737 and BP, we were able to replicate the GWAS signal reported for CACNA1C through a haplotype analysis that encompassed these previously reported significant SNPs. These results provide additional evidence that CACNA1C is associated with BP and provides the first evidence that variations in this gene might play a role in the pathogenesis of this disorder in the Latino population.

A general semi-parametric approach to the analysis of genetic association studies in population-based designs

Background

For genetic association studies in designs of unrelated individuals, current statistical methodology typically models the phenotype of interest as a function of the genotype and assumes a known statistical model for the phenotype. In the analysis of complex phenotypes, especially in the presence of ascertainment conditions, the specification of such model assumptions is not straight-forward and is error-prone, potentially causing misleading results.

Results

In this paper, we propose an alternative approach that treats the genotype as the random variable and conditions upon the phenotype. Thereby, the validity of the approach does not depend on the correctness of assumptions about the phenotypic model. Misspecification of the phenotypic model may lead to reduced statistical power. Theoretical derivations and simulation studies demonstrate both the validity and the advantages of the approach over existing methodology. In the COPDGene study (a GWAS for Chronic Obstructive Pulmonary Disease (COPD)), we apply the approach to a secondary, quantitative phenotype, the Fagerstrom nicotine dependence score, that is correlated with COPD affection status. The software package that implements this method is available.

Conclusions

The flexibility of this approach enables the straight-forward application to quantitative phenotypes and binary traits in ascertained and unascertained samples. In addition to its robustness features, our method provides the platform for the construction of complex statistical models for longitudinal data, multivariate data, multi-marker tests, rare-variant analysis, and others.

Statistical distributions of test statistics used for quantitative trait association mapping in structured populations

Background

Spurious associations between single nucleotide polymorphisms and phenotypes are a major issue in genome-wide association studies and have led to underestimation of type 1 error rate and overestimation of the number of quantitative trait loci found. Many authors have investigated the influence of population structure on the robustness of methods by simulation. This paper is aimed at developing further the algebraic formalization of power and type 1 error rate for some of the classical statistical methods used: simple regression, two approximate methods of mixed models involving the effect of a single nucleotide polymorphism (SNP) and a random polygenic effect (GRAMMAR and FASTA) and the transmission/disequilibrium test for quantitative traits and nuclear families. Analytical formulae were derived using matrix algebra for the first and second moments of the statistical tests, assuming a true mixed model with a polygenic effect and SNP effects.

Results

The expectation and variance of the test statistics and their marginal expectations and variances according to the distribution of genotypes and estimators of variance components are given as a function of the relationship matrix and of the heritability of the polygenic effect. These formulae were used to compute type 1 error rate and power for any kind of relationship matrix between phenotyped and genotyped individuals for any level of heritability. For the regression method, type 1 error rate increased with the variability of relationships and with heritability, but decreased with the GRAMMAR method and was not affected with the FASTA and quantitative transmission/disequilibrium test methods.

Conclusions

The formulae can be easily used to provide the correct threshold of type 1 error rate and to calculate the power when designing experiments or data collection protocols. The results concerning the efficacy of each method agree with simulation results in the literature but were generalized in this work. The power of the GRAMMAR method was equal to the power of the FASTA method at the same type 1 error rate. The power of the quantitative transmission/disequilibrium test was low. In conclusion, the FASTA method, which is very close to the full mixed model, is recommended in association mapping studies.

Follow-up association study of linkage regions reveals multiple candidate genes for carotid plaque in Dominicans

OBJECTIVE

Carotid plaque is a marker of subclinical atherosclerosis with a genetic component. The aim of this follow-up fine mapping study was to identify candidate genes for carotid plaque within four linkage regions.

METHODS

We successfully genotyped 3712 single nucleotide polymorphisms (SNPs) under the four linkage regions that were previously identified in 100 extended Dominican families. Family-based association tests were performed to investigate their associations with carotid plaque. Promising SNPs were evaluated in an independent population-based subcohort (N=941, 384 Dominicans) from the Northern Manhattan Study (NOMAS).

RESULTS

In the family study, evidence for association (p<0.0005) was found regarding several genes (NAV2, EFCAB11/TDP1, AGBL1, PTPN9, LINGO1 and LOC730118), with the strongest association at rs4143999 near EFCAB11/TDP1 (p=0.00001 for carotid presence and 0.00003 for plaque area, multiple testing corrected p≤0.02). The association in AGBL1 and PTPN9 was mainly driven by the families with linkage evidence (p=0.00008-0.00001 and 0.76-0.32, respectively, in the families with and without linkage evidence). However, these associations explained only a small portion of the observed linkage. In NOMAS, replication (p<0.05 in the whole NOMAS subcohort and p<0.10 in the smaller Dominican subcohort) was found for SNPs within/near EFCAB11, NAV2, AGBL1 and other genes.

CONCLUSION

This follow-up study has identified multiple candidate genes for carotid plaque in the Dominican population. Many of these genes have been implicated in neurodegenerative and cardiovascular diseases. Further studies with in-depth re-sequencing are needed to uncover both rare and common functional variants that contribute to the susceptibility to atherosclerosis.

Design and rationale of a genetic cohort study on congenital cardiac disease: experiences from a multi-institutional platform in Quebec

BACKGROUND

Congenital cardiac disease is the most common malformation, and a substantial source of mortality and morbidity in children and young adults. A role for genetic factors is recognised for these malformations, but overall few predisposing loci have been identified. Here we report the rationale, design, and first results of a multi-institutional congenital cardiac disease cohort, assembled mainly from the French-Canadian population of the province of Quebec and centred on families with multiple affected members afflicted by cardiac malformations.

METHODS

Families were recruited into the study, phenotyped and sampled for DNA in cardiology clinics over the first 3 years of enrolment. We performed segregation analysis and linkage simulations in the subgroup of families with left ventricular outflow tract obstruction (LVOTO).

RESULTS

A total of 1603 participants from 300 families were recruited, with 169 out of 300 (56.3%) families having more than one affected member. For the LVOTO group, we estimate heritability to be 0.46-0.52 in our cohort. Simulation analysis demonstrated sufficient power to carry out linkage analyses, with an expected mean log-of-odds (LOD) score of 3.8 in 67 pedigrees with LVOTO.

CONCLUSION

We show feasibility and usefulness of a population-based biobank for genetic investigations into the causes of congenital cardiac disease. Heritability of LVOTO is high and could be accounted for by multiple loci. This platform is ideally suited for multiple analysis approaches, including linkage analysis and novel gene sequencing approaches, and will allow to establish segregation of risk alleles at family and population levels.

Identifying causal rare variants of disease through family-based analysis of Genetics Analysis Workshop 17 data set

Linkage- and association-based methods have been proposed for mapping disease-causing rare variants. Based on the family information provided in the Genetic Analysis Workshop 17 data set, we formulate a two-pronged approach that combines both methods. Using the identity-by-descent information provided for eight extended pedigrees (n = 697) and the simulated quantitative trait Q1, we explore various traditional nonparametric linkage analysis methods; the best result is obtained by assuming between-family heterogeneity and applying the Haseman-Elston regression to each pedigree separately. We discover strong signals from two genes in two different families and weaker signals for a third gene from two other families. As an exploratory approach, we apply an association test based on a modified family-based association test statistic to all rare variants (frequency < 1% or < 3%) designated as causal for Q1. Family-based association tests correctly identified causal single-nucleotide polymorphisms for four genes (KDR, VEGFA, VEGFC, and FLT1). Our results suggest that both linkage and association tests with families show promise for identifying rare variants.

Family-based designs for genome-wide association studies

Association mapping has successfully identified common SNPs associated with many diseases. However, the inability of this class of variation to account for most of the supposed heritability has led to a renewed interest in methods - primarily linkage analysis - to detect rare variants. Family designs allow for control of population stratification, investigations of questions such as parent-of-origin effects and other applications that are imperfectly or not readily addressed in case-control association studies. This article guides readers through the interface between linkage and association analysis, reviews the new methodologies and provides useful guidelines for applications. Just as effective SNP-genotyping tools helped to realize the potential of association studies, next-generation sequencing tools will benefit genetic studies by improving the power of family-based approaches.

Evaluation of methods accounting for population structure with pedigree data and continuous outcomes

Methods to account for population structure (PS) in genome-wide association studies have been well developed in samples of unrelated individuals, but when a sample is composed of families, the task of finding and accounting for PS is not as straight forward. Family-based tests that condition on parental genotypes or their sufficient statistics are immune to biases due to PS, but are known to have low power, particularly for unselected samples. Population-based approaches that use all available data are an attractive alternative, but the methods have not been evaluated for continuous outcomes when a sample has both family and PS. Therefore, we compare through simulation the performance of population-based regression models that account for family and PS with continuous outcomes using a range of family sizes and structures, including two and three generational families with admixed and discrete PS. We find that when computation time is a concern, the Dupuis et al. efficient score test performs very well. When computational time is not an issue, a linear mixed effects model adjusting for genetic principal components tends to have slightly better power than the score test and may be preferred.

Joint testing of genotype and ancestry association in admixed families

Current genome-wide association studies (GWAS) often involve populations that have experienced recent genetic admixture. Genotype data generated from these studies can be used to test for association directly, as in a non-admixed population. As an alternative, these data can be used to infer chromosomal ancestry, and thus allow for admixture mapping. We quantify the contribution of allele-based and ancestry-based association testing under a family-design, and demonstrate that the two tests can provide non-redundant information. We propose a joint testing procedure, which efficiently integrates the two sources information. The efficiencies of the allele, ancestry and combined tests are compared in the context of a GWAS. We discuss the impact of population history and provide guidelines for future design and analysis of GWAS in admixed populations.

Methods: genetic epidemiology

Given the potential benefits of gene identification in psychiatry, genetic epidemiology has become a mainstream discipline within the field. This article discusses the main tools for gene discovery. The focus is on the designs and analytic approaches for each of these methods. Because most gene discovery has now moved to genetic association studies, and most recently to genome-wide association studies, the focus is on methods for this design. Also highlighted are the current challenges of genetic epidemiology as a prelude to future approaches that may be applied to psychiatric disorders in the coming years.

Immunological and genetic aspects of asthma and allergy

Prevalence of allergy and allergic asthma are increasing worldwide. More than half of the US population has a positive skin prick test and approximately 10% are asthmatics. Many studies have been conducted to define immunological pathways underlying allergy and asthma development and to identify the main genetic determinants. In the effort to find missing pieces of the puzzle, new genomic approaches and more standardized ones, such as the candidate gene approach, have been used collectively. This article proposes an overview of the actual knowledge about immunological and genetic aspects of allergy and asthma. Special attention has been drawn to the challenges linked to genetic research in complex traits such as asthma and to the contribution of new genomic approaches.

A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4

Case-parent trios were used in a genome-wide association study of cleft lip with and without cleft palate. SNPs near two genes not previously associated with cleft lip with and without cleft palate (MAFB, most significant SNP rs13041247, with odds ratio (OR) per minor allele = 0.704, 95% CI 0.635-0.778, P = 1.44 x 10(-11); and ABCA4, most significant SNP rs560426, with OR = 1.432, 95% CI 1.292-1.587, P = 5.01 x 10(-12)) and two previously identified regions (at chromosome 8q24 and IRF6) attained genome-wide significance. Stratifying trios into European and Asian ancestry groups revealed differences in statistical significance, although estimated effect sizes remained similar. Replication studies from several populations showed confirming evidence, with families of European ancestry giving stronger evidence for markers in 8q24, whereas Asian families showed stronger evidence for association with MAFB and ABCA4. Expression studies support a role for MAFB in palatal development.

Association of MET with social and communication phenotypes in individuals with autism spectrum disorder

Autism is a complex neurodevelopmental disorder diagnosed by impairments in social interaction, communication, and behavioral flexibility. Autism is highly heritable, but it is not known whether a genetic risk factor contributes to all three core domains of the disorder or autism results from the confluence of multiple genetic risk factors for each domain. We and others reported previously association of variants in the gene encoding the MET receptor tyrosine kinase in five independent samples. We further described enriched association of the MET promoter variant rs1858830 C allele in families with co-occurring autism and gastrointestinal conditions. To test the contribution of this functional MET promoter variant to the domains of autism, we analyzed its association with quantitative scores derived from three instruments used to diagnose and describe autism phenotypes: the Autism Diagnostic Interview-Revised (ADI-R), the Autism Diagnostic Observation Schedule (ADOS), and both the parent and the teacher report forms of the Social Responsiveness Scale (SRS). In 748 individuals from 367 families, the transmission of the MET C allele from parent to child was consistently associated with both social and communication phenotypes of autism. Stratification by gastrointestinal conditions revealed a similar pattern of association with both social and communication phenotypes in 242 individuals with autism from 118 families with co-occurring gastrointestinal conditions, but a lack of association with any domain in 181 individuals from 96 families with ASD and no co-occurring gastrointestinal condition. These data indicate that the MET C allele influences at least two of the three domains of the autism triad.

Methods: genetic epidemiology

Given the potential benefits of gene identification in psychiatry, genetic epidemiology has become a mainstream discipline within the field. This article discusses the main tools for gene discovery. The focus is on the designs and analytic approaches for each of these methods. Because most gene discovery has now moved to genetic association studies, and most recently to genome-wide association studies, the focus is on methods for this design. Also highlighted are the current challenges of genetic epidemiology as a prelude to future approaches that may be applied to psychiatric disorders in the coming years.