Developmental and Injury-induced Changes in DNA Methylation in Regenerative versus Non-regenerative Regions of the Vertebrate Central Nervous System.

BACKGROUND: Because some of its CNS neurons (e.g., retinal ganglion cells after optic nerve crush (ONC)) regenerate axons throughout life, whereas others (e.g., hindbrain neurons after spinal cord injury (SCI)) lose this capacity as tadpoles metamorphose into frogs, the South African claw-toed frog, Xenopus laevis, offers unique opportunities for exploring differences between regenerative and non-regenerative responses to CNS injury within the same organism. An earlier, three-way RNA-seq study (frog ONC eye, tadpole SCI hindbrain, frog SCI hindbrain) identified genes that regulate chromatin accessibility among those that were differentially expressed in regenerative vs non-regenerative CNS [11]. The current study used whole genome bisulfite sequencing (WGBS) of DNA collected from these same animals at the peak period of axon regeneration to study the extent to which DNA methylation could potentially underlie differences in chromatin accessibility between regenerative and non-regenerative CNS. RESULTS: Consistent with the hypothesis that DNA of regenerative CNS is more accessible than that of non-regenerative CNS, DNA from both the regenerative tadpole hindbrain and frog eye was less methylated than that of the non-regenerative frog hindbrain. Also, consistent with observations of CNS injury in mammals, DNA methylation in non-regenerative frog hindbrain decreased after SCI. However, contrary to expectations that the level of DNA methylation would decrease even further with axotomy in regenerative CNS, DNA methylation in these regions instead increased with injury. Injury-induced differences in CpG methylation in regenerative CNS became especially enriched in gene promoter regions, whereas non-CpG methylation differences were more evenly distributed across promoter regions, intergenic, and intragenic regions. In non-regenerative CNS, tissue-related (i.e., regenerative vs. non-regenerative CNS) and injury-induced decreases in promoter region CpG methylation were significantly correlated with increased RNA expression, but the injury-induced, increased CpG methylation seen in regenerative CNS across promoter regions was not, suggesting it was associated with increased rather than decreased chromatin accessibility. This hypothesis received support from observations that in regenerative CNS, many genes exhibiting increased, injury-induced, promoter-associated CpG-methylation also exhibited increased RNA expression and association with histone markers for active promoters and enhancers. DNA immunoprecipitation for 5hmC in optic nerve regeneration found that the promoter-associated increases seen in CpG methylation were distinct from those exhibiting changes in 5hmC. CONCLUSIONS: Although seemingly paradoxical, the increased injury-associated DNA methylation seen in regenerative CNS has many parallels in stem cells and cancer. Thus, these axotomy-induced changes in DNA methylation in regenerative CNS provide evidence for a novel epigenetic state favoring successful over unsuccessful CNS axon regeneration. The datasets described in this study should help lay the foundations for future studies of the molecular and cellular mechanisms involved. The insights gained should, in turn, help point the way to novel therapeutic approaches for treating CNS injury in mammals.

Dually localised proteins found in both the apicoplast and mitochondrion utilize the Golgi-dependent pathway for apicoplast targeting in Toxoplasma gondii.

BACKGROUND INFORMATION: Like other apicomplexan parasites, Toxoplasma gondii harbours a four-membraned endosymbiotic organelle - the apicoplast. Apicoplast proteins are nuclear encoded and trafficked to the organelle through the endoplasmic reticulum (ER). From the ER to the apicoplast, two distinct protein trafficking pathways can be used. One such pathway is the cell's secretory pathway involving the Golgi, whereas the other is a unique Golgi-independent pathway. Using different experimental approaches, many apicoplast proteins have been shown to utilize the Golgi-independent pathway, whereas a handful of reports show that a few proteins use the Golgi-dependent pathway. This has led to an emphasis towards the unique Golgi-independent pathway when apicoplast protein trafficking is discussed in the literature. Additionally, the molecular features that drive proteins to each pathway are not known. RESULTS: In this report, we systematically test eight apicoplast proteins, using a C-terminal HDEL sequence to assess the role of the Golgi in their transport. We demonstrate that dually localised proteins of the apicoplast and mitochondrion (TgSOD2, TgTPx1/2 and TgACN/IRP) are trafficked through the Golgi, whereas proteins localised exclusively to the apicoplast are trafficked independent of the Golgi. Mutants of the dually localised proteins that localised exclusively to the apicoplast also showed trafficking through the Golgi. Phylogenetic analysis of TgSOD2, TgTPx1/2 and TgACN/IRP suggested that the evolutionary origins of TgSOD2 and TgTPx1/2 lie in the mitochondrion, whereas TgACN/IRP appears to have originated from the apicoplast. CONCLUSIONS AND SIGNIFICANCE: Collectively, with these results, for the first time, we establish that the driver of the Golgi-dependent trafficking route to the apicoplast is the dual localisation of the protein to the apicoplast and the mitochondrion.

Comparative gene expression profiling between optic nerve and spinal cord injury in Xenopus laevis reveals a core set of genes inherent in successful regeneration of vertebrate central nervous system axons.

BACKGROUND: The South African claw-toed frog, Xenopus laevis, is uniquely suited for studying differences between regenerative and non-regenerative responses to CNS injury within the same organism, because some CNS neurons (e.g., retinal ganglion cells after optic nerve crush (ONC)) regenerate axons throughout life, whereas others (e.g., hindbrain neurons after spinal cord injury (SCI)) lose this capacity as tadpoles metamorphose into frogs. Tissues from these CNS regions (frog ONC eye, tadpole SCI hindbrain, frog SCI hindbrain) were used in a three-way RNA-seq study of axotomized CNS axons to identify potential core gene expression programs for successful CNS axon regeneration. RESULTS: Despite tissue-specific changes in expression dominating the injury responses of each tissue, injury-induced changes in gene expression were nonetheless shared between the two axon-regenerative CNS regions that were not shared with the non-regenerative region. These included similar temporal patterns of gene expression and over 300 injury-responsive genes. Many of these genes and their associated cellular functions had previously been associated with injury responses of multiple tissues, both neural and non-neural, from different species, thereby demonstrating deep phylogenetically conserved commonalities between successful CNS axon regeneration and tissue regeneration in general. Further analyses implicated the KEGG adipocytokine signaling pathway, which links leptin with metabolic and gene regulatory pathways, and a novel gene regulatory network with genes regulating chromatin accessibility at its core, as important hubs in the larger network of injury response genes involved in successful CNS axon regeneration. CONCLUSIONS: This study identifies deep, phylogenetically conserved commonalities between CNS axon regeneration and other examples of successful tissue regeneration and provides new targets for studying the molecular underpinnings of successful CNS axon regeneration, as well as a guide for distinguishing pro-regenerative injury-induced changes in gene expression from detrimental ones in mammals.

The transcription factor Tfap2e/AP-2ε plays a pivotal role in maintaining the identity of basal vomeronasal sensory neurons.

The identity of individual neuronal cell types is defined and maintained by the expression of specific combinations of transcriptional regulators that control cell type-specific genetic programs. The epithelium of the vomeronasal organ of mice contains two major types of vomeronasal sensory neurons (VSNs): 1) the apical VSNs which express vomeronasal 1 receptors (V1r) and the G-protein subunit Gαi2 and; 2) the basal VSNs which express vomeronasal 2 receptors (V2r) and the G-protein subunit Gαo. Both cell types originate from a common pool of progenitors and eventually acquire apical or basal identity through largely unknown mechanisms. The transcription factor AP-2ε, encoded by the Tfap2e gene, plays a role in controlling the development of GABAergic interneurons in the main and accessory olfactory bulb (AOB), moreover AP-2ε has been previously described to be expressed in the basal VSNs. Here we show that AP-2ε is expressed in post-mitotic VSNs after they commit to the basal differentiation program. Loss of AP-2ε function resulted in reduced number of basal VSNs and in an increased number of neurons expressing markers of the apical lineage. Our work suggests that AP-2ε, which is expressed in late phases of differentiation, is not needed to initiate the apical-basal differentiation dichotomy but for maintaining the basal VSNs' identity. In AP-2ε mutants we observed a large number of cells that entered the basal program can express apical genes, our data suggest that differentiated VSNs of mice retain a notable level of plasticity.

Clinical utility of exome sequencing in individuals with large homozygous regions detected by chromosomal microarray analysis.

BACKGROUND: Chromosomal microarray analysis (CMA) is recommended as the first-tier clinical diagnostic test for individuals with developmental disabilities. In addition to detecting copy number variations, CMA platforms with single nucleotide polymorphism probes can detect large homozygous regions within the genome, which represent potential risk for recessively inherited disorders. METHODS: To determine the frequency in which pathogenic or likely pathogenic variants can be detected in these regions of homozygosity, we performed whole exome sequencing (WES) in 53 individuals where homozygosity was detected by CMA. These patients were referred to our clinical laboratory for a variety of neurodevelopmental conditions including autism spectrum disorder, developmental delay, epilepsy, intellectual disability and microcephaly. RESULTS: In 11.3% (6/53) of cases, the analysis of homozygous variants revealed pathogenic or likely pathogenic variants in GJB2, TPP1, SLC25A15, TYR, PCCB, and NDUFV2 which are implicated in a variety of diseases. The evaluation of heterozygous variants with autosomal dominant inheritance, compound heterozygotes and variants with X-linked inheritance revealed pathogenic or likely pathogenic variants in PNPLA4, CADM1, HBB, SOS1, SFTPC, OTC and ASMT in 15.1% (8/53) of cases. Two of these patients harbored both homozygous and heterozygous variants relevant to their phenotypes (TPP1 and OTC; GJB2 and ASMT). CONCLUSIONS: Our study highlights the clinical utility of WES in individuals whose CMA uncovers homozygosity. Importantly, we show that when the phenotype is complex and homozygosity levels are high, WES can identify a significant number of relevant variants that explain neurodevelopmental phenotypes, and these mutations may lie outside of the regions of homozygosity, suggesting that the appropriate follow up test is WES rather than targeted sequencing.

Chromosomal Microarray Analysis of Consecutive Individuals with Autism Spectrum Disorders Using an Ultra-High Resolution Chromosomal Microarray Optimized for Neurodevelopmental Disorders.

Copy number variants (CNVs) detected by chromosomal microarray analysis (CMA) significantly contribute to understanding the etiology of autism spectrum disorder (ASD) and other related conditions. In recognition of the value of CMA testing and its impact on medical management, CMA is in medical guidelines as a first-tier test in the evaluation of children with these disorders. As CMA becomes adopted into routine care for these patients, it becomes increasingly important to report these clinical findings. This study summarizes the results of over 4 years of CMA testing by a CLIA-certified clinical testing laboratory. Using a 2.8 million probe microarray optimized for the detection of CNVs associated with neurodevelopmental disorders, we report an overall CNV detection rate of 28.1% in 10,351 consecutive patients, which rises to nearly 33% in cases without ASD, with only developmental delay/intellectual disability (DD/ID) and/or multiple congenital anomalies (MCA). The overall detection rate for individuals with ASD is also significant at 24.4%. The detection rate and pathogenic yield of CMA vary significantly with the indications for testing, age, and gender, as well as the specialty of the ordering doctor. We note discrete differences in the most common recurrent CNVs found in individuals with or without a diagnosis of ASD.

Rare exonic deletions implicate the synaptic organizer Gephyrin (GPHN) in risk for autism, schizophrenia and seizures.

The GPHN gene codes for gephyrin, a key scaffolding protein in the neuronal postsynaptic membrane, responsible for the clustering and localization of glycine and GABA receptors at inhibitory synapses. Gephyrin has well-established functional links with several synaptic proteins that have been implicated in genetic risk for neurodevelopmental disorders such as autism spectrum disorder (ASD), schizophrenia and epilepsy including the neuroligins (NLGN2, NLGN4), the neurexins (NRXN1, NRXN2, NRXN3) and collybistin (ARHGEF9). Moreover, temporal lobe epilepsy has been linked to abnormally spliced GPHN mRNA lacking exons encoding the G-domain of the gephyrin protein, potentially arising due to cellular stress associated with epileptogenesis such as temperature and alkalosis. Here, we present clinical and genomic characterization of six unrelated subjects, with a range of neurodevelopmental diagnoses including ASD, schizophrenia or seizures, who possess rare de novo or inherited hemizygous microdeletions overlapping exons of GPHN at chromosome 14q23.3. The region of common overlap across the deletions encompasses exons 3-5, corresponding to the G-domain of the gephyrin protein. These findings, together with previous reports of homozygous GPHN mutations in connection with autosomal recessive molybdenum cofactor deficiency, will aid in clinical genetic interpretation of the GPHN mutation spectrum. Our data also add to the accumulating evidence implicating neuronal synaptic gene products as key molecular factors underlying the etiologies of a diverse range of neurodevelopmental conditions.

SHANK1 Deletions in Males with Autism Spectrum Disorder.

Recent studies have highlighted the involvement of rare (<1% frequency) copy-number variations and point mutations in the genetic etiology of autism spectrum disorder (ASD); these variants particularly affect genes involved in the neuronal synaptic complex. The SHANK gene family consists of three members (SHANK1, SHANK2, and SHANK3), which encode scaffolding proteins required for the proper formation and function of neuronal synapses. Although SHANK2 and SHANK3 mutations have been implicated in ASD and intellectual disability, the involvement of SHANK1 is unknown. Here, we assess microarray data from 1,158 Canadian and 456 European individuals with ASD to discover microdeletions at the SHANK1 locus on chromosome 19. We identify a hemizygous SHANK1 deletion that segregates in a four-generation family in which male carriers--but not female carriers--have ASD with higher functioning. A de novo SHANK1 deletion was also detected in an unrelated male individual with ASD with higher functioning, and no equivalent SHANK1 mutations were found in >15,000 controls (p = 0.009). The discovery of apparent reduced penetrance of ASD in females bearing inherited autosomal SHANK1 deletions provides a possible contributory model for the male gender bias in autism. The data are also informative for clinical-genetics interpretations of both inherited and sporadic forms of ASD involving SHANK1.

Rare deletions at the neurexin 3 locus in autism spectrum disorder.

The three members of the human neurexin gene family, neurexin 1 (NRXN1), neurexin 2 (NRXN2), and neurexin 3 (NRXN3), encode neuronal adhesion proteins that have important roles in synapse development and function. In autism spectrum disorder (ASD), as well as in other neurodevelopmental conditions, rare exonic copy-number variants and/or point mutations have been identified in the NRXN1 and NRXN2 loci. We present clinical characterization of four index cases who have been diagnosed with ASD and who possess rare inherited or de novo microdeletions at 14q24.3-31.1, a region that overlaps exons of the alpha and/or beta isoforms of NRXN3. NRXN3 deletions were found in one father with subclinical autism and in a carrier mother and father without formal ASD diagnoses, indicating issues of penetrance and expressivity at this locus. Notwithstanding these clinical complexities, this report on ASD-affected individuals who harbor NRXN3 exonic deletions advances the understanding of the genetic etiology of autism, further enabling molecular diagnoses.

Utility of hand-held echocardiography in outpatient pediatric cardiology management.

Adult patient series have shown hand-held echocardiography (echo) units (HHE) to be accurate for rapid diagnosis and triage. This is the first study to evaluate the ability of HHE to inform decision making in outpatient pediatric cardiology. New pediatric cardiology patients in outpatient clinics staffed by six pediatric cardiologists (experience 1-17 years) were prospectively enrolled if an echocardiogram (echo) was ordered during their initial visit. After history and physical examination and before a standard echo, the cardiologists performed a bedside HHE examination (GE Vscan 1.7-3.8 MHz), documented findings, and made a clinical decision. Diagnoses and decisions based on HHE were compared with final management after the standard echo. The study enrolled 101 subjects (ages 9 days to 19 years). The cardiologists considered HHE imaging adequate for decision making for 80 of the 101 subjects. For 77 of the 80 subjects with acceptable HHE imaging (68/68 normal and 9/12 abnormal standard echoes), the HHE-based primary diagnoses and decisions agreed with the final management. The sensitivity of HHE was 75 % (95 % confidence interval [CI] 43-94 %) and the positive predictive value 100 % (95 % CI 66-100 %) for pediatric heart disease. The agreement between standard echocardiography and HHE imaging was substantial (κ = 0.82). Excluding one of the least experienced cardiologists, HHE provided the basis for correct cardiac diagnoses and management for all the subjects with acceptable HHE imaging (58/58 normal and 9/9 abnormal echoes). In outpatient pediatric cardiology, HHE has potential as a tool to complement physical examination. Further investigation is needed to evaluate how value improves with clinical experience.

Anergy in CD4 memory T lymphocytes. II. Abrogation of TCR-induced formation of membrane signaling complexes.

Memory and naive CD4 T cells have unique regulatory pathways for self/non-self discrimination. A memory cell specific regulatory pathway was revealed using superantigens to trigger the TCR. Upon stimulation by bacterial superantigens, like staphylococcal enterotoxin B (SEB), TCR proximal signaling is impaired leading to clonal tolerance (anergy). In the present report, we show that memory cell anergy results from the sequestration of the protein tyrosine kinase ZAP-70 away from the TCR/CD3ζ chain. During SEB-induced signaling, ZAP-70 is excluded from both detergent-resistant membrane microdomains and the immunological synapse, thus blocking downstream signaling. We also show that the mechanism underlying memory cell anergy must involve Fyn kinase, given that the suppression of Fyn activity restores the movement of ZAP-70 to the immunological synapse, TCR proximal signaling, and cell proliferation. Thus, toleragens, including microbial toxins, may modulate memory responses by targeting the organizational structure of memory cell signaling complexes.

Evaluating the role of reference-genome phylogenetic distance on evolutionary inference.

When a high-quality genome assembly of a target species is unavailable, an option to avoid the costly de novo assembly process is a mapping-based assembly. However, mapping shotgun data to a distant relative may lead to biased or erroneous evolutionary inference. Here, we used short-read data from a mammal (beluga whale) and a bird species (rowi kiwi) to evaluate whether reference genome phylogenetic distance can impact downstream demographic (Pairwise Sequentially Markovian Coalescent) and genetic diversity (heterozygosity, runs of homozygosity) analyses. We mapped to assemblies of species of varying phylogenetic distance (from conspecific to genome-wide divergence of >7%), and de novo assemblies created using cross-species scaffolding. We show that while reference genome phylogenetic distance has an impact on demographic analyses, it is not pronounced until using a reference genome with >3% divergence from the target species. When mapping to cross-species scaffolded assemblies, we are unable to replicate the original beluga demographic results, but are able with the rowi kiwi, presumably reflecting the more fragmented nature of the beluga assemblies. We find that increased phylogenetic distance has a pronounced impact on genetic diversity estimates; heterozygosity estimates deviate incrementally with increasing phylogenetic distance. Moreover, runs of homozygosity are largely undetectable when mapping to any nonconspecific assembly. However, these biases can be reduced when mapping to a cross-species scaffolded assembly. Taken together, our results show that caution should be exercised when selecting reference genomes. Cross-species scaffolding may offer a way to avoid a costly, traditional de novo assembly, while still producing robust, evolutionary inference.

Mechanisms of cell death of thymocytes induced by polyunsaturated, monounsaturated and trans-fatty acids.

Polyunsaturated fatty acids (PUFAs) are rapidly cytotoxic to isolated murine thymocytes, and the degree of cell death has been correlated with changes in membrane fluidity, elevation of intracellular calcium concentration and generation of reactive oxygen species. We have compared the degree of cell death and increase in membrane fluidity of C-20 and C-22 omega-3 and 6 PUFAs to those induced by monounsaturated and trans-fatty acids, and find that concentrations which induce comparable increases in membrane fluidity do not cause comparable cell death. The C-18 omega-6 causes a decrease in membrane fluidity, yet is the most potent in causing cell death. Omega-6 PUFAs are more cytotoxic than omega-3 PUFAs, while monounsaturated and trans-fats show little cytotoxicity and only at much higher concentrations. Cell death is preceded by reductions of both plasma and mitochondrial membrane potential, and occurs via apoptosis. These results indicate that cell death is due to mechanisms other than changes in membrane fluidity.

Role of calcium and ROS in cell death induced by polyunsaturated fatty acids in murine thymocytes.

We investigated the mechanisms whereby omega-3 and -6 polyunsaturated fatty acids (PUFAs) cause cell death of mouse thymocytes using flow cytometry, focusing on the respective roles of intracellular calcium concentration, [Ca(2+)](i) and reactive oxygen species (ROS). We applied the C-22, 20, and 18 carbon omega-3 (DHA, EPA, ALA) and omega-6 (DTA, ARA, and LNA) fatty acids to isolated thymocytes and monitored cell death using the DNA-binding dye, propidium iodide. When applied at 20 µM concentration, omega-3 fatty acids killed thymocytes over a period of 1 h with a potency of DHA > EPA > ALA. The omega-6 PUFAs were more potent. The C18 omega-6 fatty acid, LNA, was the most potent, followed by DHA and ARA. Cell death was always accompanied by an increase in the levels of [Ca(2+)](i) and ROS. Both increases were in proportion to the potency of the PUFAs in inducing cell death. Removing extracellular calcium did not prevent the elevation in [Ca(2+)](i) nor cell death. However, the intracellular calcium chelator, BAPTA, almost totally reduced both the elevation in [Ca(2+)](i) and cell death, while vitamin E reduced the elevation in ROS and cell death. BAPTA also prevented the elevation in ROS, but vitamin E did not prevent the elevation in [Ca(2+)](i). Thapsigargin, which depletes endoplasmic reticulum calcium, blocked the elevation in [Ca(2+)](i), but CCCP, a mitochondrial calcium uptake inhibitor, did not. These results suggest that the six PUFAs we studied kill thymocytes by causing release of calcium from endoplasmic reticulum, which causes release of ROS from mitochondria which leads to cell death.

A 2cM genome-wide scan of European Holstein cattle affected by classical BSE.

BACKGROUND: Classical bovine spongiform encephalopathy (BSE) is an acquired prion disease that is invariably fatal in cattle and has been implicated as a significant human health risk. Polymorphisms that alter the prion protein of sheep or humans have been associated with variations in transmissible spongiform encephalopathy susceptibility or resistance. In contrast, there is no strong evidence that non-synonymous mutations in the bovine prion gene (PRNP) are associated with classical BSE disease susceptibility. However, two bovine PRNP insertion/deletion polymorphisms, one within the promoter region and the other in intron 1, have been associated with susceptibility to classical BSE. These associations do not explain the full extent of BSE susceptibility, and loci outside of PRNP appear to be associated with disease incidence in some cattle populations. To test for associations with BSE susceptibility, we conducted a genome wide scan using a panel of 3,072 single nucleotide polymorphism (SNP) markers on 814 animals representing cases and control Holstein cattle from the United Kingdom BSE epidemic. RESULTS: Two sets of BSE affected Holstein cattle were analyzed in this study, one set with known family relationships and the second set of paired cases with controls. The family set comprises half-sibling progeny from six sires. The progeny from four of these sires had previously been scanned with microsatellite markers. The results obtained from the current analysis of the family set yielded both some supporting and new results compared with those obtained in the earlier study. The results revealed 27 SNPs representing 18 chromosomes associated with incidence of BSE disease. These results confirm a region previously reported on chromosome 20, and identify additional regions on chromosomes 2, 14, 16, 21 and 28. This study did not identify a significant association near the PRNP in the family sample set. The only association found in the PRNP region was in the case-control sample set and this was not significant after multiple test correction. The genome scan of the case-control animals did not identify any associations that passed a stringent genome-wide significance threshold. CONCLUSIONS: Several regions of the genome are statistically associated with the incidence of classical BSE in European Holstein cattle. Further investigation of loci on chromosomes 2, 14, 16, 20, 21 and 28 will be required to uncover any biological significance underlying these marker associations.

The Temple Grandin Genome: Comprehensive Analysis in a Scientist with High-Functioning Autism.

Autism spectrum disorder (ASD) is a heterogeneous condition with a complex genetic etiology. The objective of this study is to identify the complex genetic factors that underlie the ASD phenotype and other clinical features of Professor Temple Grandin, an animal scientist and woman with high-functioning ASD. Identifying the underlying genetic cause for ASD can impact medical management, personalize services and treatment, and uncover other medical risks that are associated with the genetic diagnosis. Prof. Grandin underwent chromosomal microarray analysis, whole exome sequencing, and whole genome sequencing, as well as a comprehensive clinical and family history intake. The raw data were analyzed in order to identify possible genotype-phenotype correlations. Genetic testing identified variants in three genes (SHANK2, ALX1, and RELN) that are candidate risk factors for ASD. We identified variants in MEFV and WNT10A, reported to be disease-associated in previous studies, which are likely to contribute to some of her additional clinical features. Moreover, candidate variants in genes encoding metabolic enzymes and transporters were identified, some of which suggest potential therapies. This case report describes the genomic findings in Prof. Grandin and it serves as an example to discuss state-of-the-art clinical diagnostics for individuals with ASD, as well as the medical, logistical, and economic hurdles that are involved in clinical genetic testing for an individual on the autism spectrum.

Clinical significance of copy number variants involving KANK1 in patients with neurodevelopmental disorders.

Copy number variants (CNV)s involving KANK1 are generally classified as variants of unknown significance. Several clinical case reports suggest that the loss of KANK1 on chromosome 9p24.3 has potential impact on neurodevelopment. These case studies are inconsistent in terms of patient phenotype and suspected pattern of inheritance. Further complexities arise because these published reports utilize a variety of genetic testing platforms with varying resolution of the 9p region; this ultimately causes uncertainty about the impacted genomic coordinates and gene transcripts. Beyond these case reports, large case-control studies and publicly available databases statistically cast doubt as to whether variants of KANK1 are clinically significant. However, these large data sources are neither easily extracted nor uniformly applied to clinical interpretation. In this report we provide an updated analysis of the data on this locus and its potential clinical relevance. This is based on a review of the literature as well as 28 patients who harbor a single copy number variant involving KANK1 with or without DOCK8 (27 of whom are not published previously) identified by our clinical laboratory using an ultra-high resolution chromosomal microarray analysis. We note that 13 of 16 patients have a documented diagnosis of autism spectrum disorder (ASD) while only two, with documented perinatal complications, have a documented diagnosis of cerebral palsy (CP). A careful review of the CNVs suggests a transcript-specific effect. After evaluation of our case series and reconsideration of the literature, we propose that KANK1 aberrations do not frequently cause CP but cannot exclude that they represent a risk factor for ASD, especially when the coding region of the shorter, alternate KANK1 transcript (termed "transcript 4" in the UCSC Genome Browser) is impacted.

Critical exon indexing improves clinical interpretation of copy number variants in neurodevelopmental disorders.

OBJECTIVE: To evaluate a new tool to aid interpretation of copy number variants (CNVs) in individuals with neurodevelopmental disabilities. METHODS: Critical exon indexing (CEI) was used to identify genes with critical exons (CEGs) from clinically reported CNVs, which may contribute to neurodevelopmental disorders (NDDs). The 742 pathogenic CNVs and 1,363 variants of unknown significance (VUS) identified by chromosomal microarray analysis in 5,487 individuals with NDDs were subjected to CEI to identify CEGs. CEGs identified in a subsequent random series of VUS were evaluated for relevance to CNV interpretation. RESULTS: CEI identified a total of 2,492 unique CEGs in pathogenic CNVs and 953 in VUS compared with 259 CEGs in 6,965 CNVs from 873 controls. These differences are highly significant (p < 0.00001) whether compared as frequency, average, or normalized by CNV size. Twenty-one percent of VUS CEGs were not represented in Online Mendelian Inheritance in Man, highlighting limitations of existing resources for identifying potentially impactful genes within CNVs. CEGs were highly correlated with other indices and known pathways of relevance. Separately, 136 random VUS reports were reevaluated, and 76% of CEGs had not been commented on. In multiple cases, further investigation yielded additional relevant literature aiding interpretation. As one specific example, we discuss GTF2I as a CEG, which likely alters interpretation of several reported duplication VUS in the Williams-Beuren region. CONCLUSIONS: Application of CEI to CNVs in individuals with NDDs can identify genes of potential clinical relevance, aid laboratories in effectively searching the clinical literature, and support the clinical reporting of poorly annotated VUS.

Rare copy number variations affecting the synaptic gene DMXL2 in neurodevelopmental disorders.

BACKGROUND: Ultra-rare genetic variants, including non-recurrent copy number variations (CNVs) affecting important dosage-sensitive genes, are important contributors to the etiology of neurodevelopmental disorders (NDDs). Pairing family-based whole-genome sequencing (WGS) with detailed phenotype data can enable novel gene associations in NDDs. METHODS: We performed WGS of six members from a three-generation family, where three individuals each had a spectrum of features suggestive of a NDD. CNVs and sequence-level variants were identified and further investigated in disease and control databases. RESULTS: We identified a novel 252-kb deletion at 15q21 that overlaps the synaptic gene DMXL2 and the gene GLDN. The microdeletion segregated in NDD-affected individuals. Additional rare inherited and de novo sequence-level variants were found that may also be involved, including a missense change in GRIK5. Multiple CNVs and loss-of-function sequence variants affecting DMXL2 were discovered in additional unrelated individuals with a range of NDDs. CONCLUSIONS: Disruption of DMXL2 may predispose to NDDs including autism spectrum disorder. The robust interpretation of private variants requires a multifaceted approach that incorporates multigenerational pedigrees and genome-wide and population-scale data.

A first generation whole genome RH map of the river buffalo with comparison to domestic cattle.

BACKGROUND: The recently constructed river buffalo whole-genome radiation hybrid panel (BBURH5000) has already been used to generate preliminary radiation hybrid (RH) maps for several chromosomes, and buffalo-bovine comparative chromosome maps have been constructed. Here, we present the first-generation whole genome RH map (WG-RH) of the river buffalo generated from cattle-derived markers. The RH maps aligned to bovine genome sequence assembly Btau_4.0, providing valuable comparative mapping information for both species. RESULTS: A total of 3990 markers were typed on the BBURH5000 panel, of which 3072 were cattle derived SNPs. The remaining 918 were classified as cattle sequence tagged site (STS), including coding genes, ESTs, and microsatellites. Average retention frequency per chromosome was 27.3% calculated with 3093 scorable markers distributed in 43 linkage groups covering all autosomes (24) and the X chromosomes at a LOD >or= 8. The estimated total length of the WG-RH map is 36,933 cR5000. Fewer than 15% of the markers (472) could not be placed within any linkage group at a LOD score >or= 8. Linkage group order for each chromosome was determined by incorporation of markers previously assigned by FISH and by alignment with the bovine genome sequence assembly (Btau_4.0). CONCLUSION: We obtained radiation hybrid chromosome maps for the entire river buffalo genome based on cattle-derived markers. The alignments of our RH maps to the current bovine genome sequence assembly (Btau_4.0) indicate regions of possible rearrangements between the chromosomes of both species. The river buffalo represents an important agricultural species whose genetic improvement has lagged behind other species due to limited prior genomic characterization. We present the first-generation RH map which provides a more extensive resource for positional candidate cloning of genes associated with complex traits and also for large-scale physical mapping of the river buffalo genome.