Association of rare missense variants in the second intracellular loop of NaV1.7 sodium channels with familial autism.

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder often accompanied by intellectual disability, language impairment and medical co-morbidities. The heritability of autism is high and multiple genes have been implicated as causal. However, most of these genes have been identified in de novo cases. To further the understanding of familial autism, we performed whole-exome sequencing on five families in which second- and third-degree relatives were affected. By focusing on novel and protein-altering variants, we identified a small set of candidate genes. Among these, a novel private missense C1143F variant in the second intracellular loop of the voltage-gated sodium channel NaV1.7, encoded by the SCN9A gene, was identified in one family. Through electrophysiological analysis, we show that NaV1.7C1143F exhibits partial loss-of-function effects, resulting in slower recovery from inactivation and decreased excitability in cultured cortical neurons. Furthermore, for the same intracellular loop of NaV1.7, we found an excess of rare variants in a case-control variant-burden study. Functional analysis of one of these variants, M932L/V991L, also demonstrated reduced firing in cortical neurons. However, although this variant is rare in Caucasians, it is frequent in Latino population, suggesting that genetic background can alter its effects on phenotype. Although the involvement of the SCN1A and SCN2A genes encoding NaV1.1 and NaV1.2 channels in de novo ASD has previously been demonstrated, our study indicates the involvement of inherited SCN9A variants and partial loss-of-function of NaV1.7 channels in the etiology of rare familial ASD.

Identification of rare variants from exome sequence in a large pedigree with autism.

We carried out analyses with the goal of identifying rare variants in exome sequence data that contribute to disease risk for a complex trait. We analyzed a large, 47-member, multigenerational pedigree with 11 cases of autism spectrum disorder, using genotypes from 3 technologies representing increasing resolution: a multiallelic linkage marker panel, a dense diallelic marker panel, and variants from exome sequencing. Genome-scan marker genotypes were available on most subjects, and exome sequence data was available on 5 subjects. We used genome-scan linkage analysis to identify and prioritize the chromosome 22 region of interest, and to select subjects for exome sequencing. Inheritance vectors (IVs) generated by Markov chain Monte Carlo analysis of multilocus marker data were the foundation of most analyses. Genotype imputation used IVs to determine which sequence variants reside on the haplotype that co-segregates with the autism diagnosis. Together with a rare-allele frequency filter, we identified only one rare variant on the risk haplotype, illustrating the potential of this approach to prioritize variants. The associated gene, MYH9, is biologically unlikely, and we speculate that for this complex trait, the key variants may lie outside the exome.

Family-based genome scan for age at onset of late-onset Alzheimer's disease in whole exome sequencing data.

Alzheimer's disease (AD) is a common and complex neurodegenerative disease. Age at onset (AAO) of AD is an important component phenotype with a genetic basis, and identification of genes in which variation affects AAO would contribute to identification of factors that affect timing of onset. Increase in AAO through prevention or therapeutic measures would have enormous benefits by delaying AD and its associated morbidities. In this paper, we performed a family-based genome-wide association study for AAO of late-onset AD in whole exome sequence data generated in multigenerational families with multiple AD cases. We conducted single marker and gene-based burden tests for common and rare variants, respectively. We combined association analyses with variance component linkage analysis, and with reference to prior studies, in order to enhance evidence of the identified genes. For variants and genes implicated by the association study, we performed a gene-set enrichment analysis to identify potential novel pathways associated with AAO of AD. We found statistically significant association with AAO for three genes (WRN, NTN4 and LAMC3) with common associated variants, and for four genes (SLC8A3, SLC19A3, MADD and LRRK2) with multiple rare-associated variants that have a plausible biological function related to AD. The genes we have identified are in pathways that are strong candidates for involvement in the development of AD pathology and may lead to a better understanding of AD pathogenesis.

Growth hormone insufficiency associated with haploinsufficiency at 18q23.

Growth hormone insufficiency is a common cause of growth failure in children with the 18q- syndrome. Individuals with this syndrome have a deletion as large as 36 Mb from the long arm of chromosome 18. We have evaluated 33 children with this syndrome for growth hormone production and have identified a region of approximately 2 Mb, which is deleted in every growth hormone insufficient patient. Two genes contained in this region, myelin basic protein, and the galanin receptor, are candidate genes for the growth hormone insufficiency phenotype.

Preferential loss of the paternal alleles in the 18q- syndrome.

Individuals with the 18q- syndrome have variable deletions from the long arm of chromosome 18. They also exhibit a highly variable phenotype. To correlate genotype with phenotype accurately, extensive molecular and phenotypic analyses are needed on each affected individual. As a part of this analysis, we have determined the parental origin of the deleted chromosome in 34 individuals with the 18q- syndrome. We have found that 85% of the de novo deletions are paternal in origin. The percentage of fathers of individuals with paternally derived deletions who were > 30 years old was (not significantly) greater than that of the general population. The mothers of individuals with maternally derived deletions were near an average age for childbearing compared to the general population. Individuals with maternally derived terminal deletions had breakpoints as varied as those with paternally derived deletions. These results are consistent with the hypothesis that the reduced incidence of maternally derived deletions is not due to reduced viability, since individuals with large maternally derived deletions of chromosome 18q were found. We hypothesize that the prevalence of paternally derived deletions is due to an increased frequency of chromosome breakage in male germ cells. These results are consistent with results observed in other segmental aneusomies in which there is a high incidence of paternally derived deletions.

A genome scan in multigenerational families with dyslexia: Identification of a novel locus on chromosome 2q that contributes to phonological decoding efficiency.

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

Identification of cryptic rearrangements in patients with 18q- deletion syndrome.

The majority of patients with 18q- syndrome appear cytogenetically to have a terminal deletion of the long arm of chromosome 18. These 18q- patients are diagnosed by use of standard cytogenetic banding techniques, which have resolution insufficient for precise genotyping. In our effort to obtain a thorough genotype, we have analyzed the DNA from 35 patients who originally were diagnosed as having de novo terminal deletions of chromosome 18. Molecular analysis was performed with polymorphic markers throughout the 18q- region. Cytogenetic FISH was performed with two human 18q telomeric probes, a chromosome 18-specific alpha-satellite probe, and whole chromosome 18-specific paint. Of 35 patients previously reported to have terminal deletions of 18q, we found that 5 (14%) have more-complex cryptic rearrangements and that 3 (9%) retain the most distal portion of 18q, consistent with an interstitial rather than a terminal deletion. These findings indicate that a standard karyotype can lead to insufficient characterization in 18q- syndrome. This has important ramifications for phenotype mapping of this syndrome, as well as for proper prognosis.

Evidence for a novel osteosarcoma tumor-suppressor gene in the chromosome 18 region genetically linked with Paget disease of bone.

Paget disease of bone, or "osteitis deformans," is a bone disorder characterized by rapid bone remodeling resulting in abnormal bone formation. It is the second most common metabolic bone disease after osteoporosis, affecting 3%-5% of subjects aged >40 years. Recent evidence suggests that predisposition to Paget disease may have a genetic component. Genetic linkage analysis of families with multigenerational Paget disease shows linkage to a region of chromosome 18q near the polymorphic locus D18S42. Approximately 1% of Paget patients develop osteosarcoma, which represents an increase in risk that is several thousandfold over that of the general population. Osteosarcoma in Paget patients is the underlying basis for a significant fraction of osteosarcomas occurring after age 60 years. Our analysis of tumor-specific loss of constitutional heterozygosity (LOH) in 96 sporadic osteosarcomas has identified a putative tumor-suppressor locus that maps to chromosome 18q. We have localized this tumor-suppressor locus between D18S60 and D18S42, a region tightly linked to familial Paget disease. Analysis of osteosarcomas from patients with Paget disease revealed that these tumors also undergo LOH in this region. These findings suggest that the association between Paget disease and osteosarcoma is the result of a single gene or two tightly linked genes on chromosome 18.

Genetic pattern of prostate cancer progression.

Genetic alterations in primary prostate cancer (CaP) have been extensively studied, yet little is known about the genetic mechanisms underlying progression of primary CaP to metastatic prostate cancer. As a result, it is not possible to distinguish clinically indolent localized disease from potentially life-threatening tumors with high metastatic potential. To address this question, we collected tissue from 34 autopsy-derived metastases, samples rarely analyzed in previous studies. These were compared to a separate set of 17 prostatectomy specimens containing 22 foci of CaP associated with 49 examples of high-grade prostatic intraepithelial neoplasia (PIN), a histological precursor of CaP. We compared the loss of heterozygosity (LOH) profiles of high-grade PIN, primary CaP and metastases by analyzing 33 microsatellite markers previously found to have high frequencies of LOH in primary CaP. These markers were on chromosomes 5q, 6q, 7q, 8p, 9p, 10q, 11p, 13q, 16q, 17, 18q and 21q. In addition, markers on chromosomes 4p, 11q, 14q and 20q with no reported LOH in primary CaP were analyzed to determine the frequency of background LOH. In PIN lesions, the rate of LOH was significant only at D5S806 (20%) and D16S422 (29%). In addition, different PIN lesions within the same prostate gland were genetically diverse, indicating divergent evolution of synchronous neoplastic precursor lesions. LOH frequency was progressively higher in primary CaP and metastatic lesions. In primary CaP, significant losses occurred at the 8p, 10q, 11p, 16q, 17p, 18q and 21q loci (range 17-43%). Distinct patterns of LOH frequencies were observed in primary CaP compared with metastases. Although some loci (D16S422, D17S960, D21S156) showed similar frequencies of LOH in primary CaP and metastatic CaP, most other loci showed up to 7-fold metastasis-related increases. The metastatic samples revealed previously unrecognized prostate cancer LOH at D5S806, D6S262, D9S157, D13S133 and D13S227. These significant stage-specific differences in LOH frequency specify genetic loci that may play key roles in CaP progression and could represent clinically useful biomarkers for CaP aggressiveness.