Mapping recurrent mosaic copy number variation in human neurons

When somatic cells acquire complex karyotypes, they often are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons are limited by relatively small sample sizes. Here, we develop an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We apply this approach to 2,125 frontal cortical neurons from a neurotypical human brain. SCOVAL identifies 226 CNV neurons, which include a subclass of 65 CNV neurons with highly aberrant karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we find that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contain fewer, but longer, genes.

Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Mapping the Complex Genetic Landscape of Human Neurons

When somatic cells acquire complex karyotypes, they are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons have been limited by relatively small sample sizes. Here, we developed an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We applied this approach to 2,125 frontal cortical neurons from a neurotypical human brain. This approach identified 226 CNV neurons, as well as a class of CNV neurons with complex karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we found that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contained fewer, but longer, genes.

Comprehensive identification of somatic nucleotide variants in human brain tissue

BACKGROUND

Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells.

RESULTS

Here, the Brain Somatic Mosaicism Network conducts a coordinated, multi-institutional study to examine the ability of existing methods to detect simulated somatic single-nucleotide variants (SNVs) in DNA mixing experiments, generate multiple replicates of whole-genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual, devise strategies to discover somatic SNVs, and apply various approaches to validate somatic SNVs. These efforts lead to the identification of 43 bona fide somatic SNVs that range in variant allele fractions from ~ 0.005 to ~ 0.28. Guided by these results, we devise best practices for calling mosaic SNVs from 250× whole-genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrate that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees.

CONCLUSIONS

This study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Identification and characterization of occult human-specific LINE-1 insertions using long-read sequencing technology

Long Interspersed Element-1 (LINE-1) retrotransposition contributes to inter- and intra-individual genetic variation and occasionally can lead to human genetic disorders. Various strategies have been developed to identify human-specific LINE-1 (L1Hs) insertions from short-read whole genome sequencing (WGS) data; however, they have limitations in detecting insertions in complex repetitive genomic regions. Here, we developed a computational tool (PALMER) and used it to identify 203 non-reference L1Hs insertions in the NA12878 benchmark genome. Using PacBio long-read sequencing data, we identified L1Hs insertions that were absent in previous short-read studies (90/203). Approximately 81% (73/90) of the L1Hs insertions reside within endogenous LINE-1 sequences in the reference assembly and the analysis of unique breakpoint junction sequences revealed 63% (57/90) of these L1Hs insertions could be genotyped in 1000 Genomes Project sequences. Moreover, we observed that amplification biases encountered in single-cell WGS experiments led to a wide variation in L1Hs insertion detection rates between four individual NA12878 cells; under-amplification limited detection to 32% (65/203) of insertions, whereas over-amplification increased false positive calls. In sum, these data indicate that L1Hs insertions are often missed using standard short-read sequencing approaches and long-read sequencing approaches can significantly improve the detection of L1Hs insertions present in individual genomes.

Distinct Nature of Static and Dynamic Magnetic Stripes in Cuprate Superconductors

We present detailed neutron scattering studies of the static and dynamic stripes in an optimally doped high-temperature superconductor, La_{2}CuO_{4+y}. We observe that the dynamic stripes do not disperse towards the static stripes in the limit of vanishing energy transfer. Therefore, the dynamic stripes observed in neutron scattering experiments are not the Goldstone modes associated with the broken symmetry of the simultaneously observed static stripes, and the signals originate from different domains in the sample. These observations support real-space electronic phase separation in the crystal, where the static stripes in one phase are pinned versions of the dynamic stripes in the other, having slightly different periods. Our results explain earlier observations of unusual dispersions in underdoped La_{2-x}Sr_{x}CuO_{4} (x=0.07) and La_{2-x}Ba_{x}CuO_{4} (x=0.095).

Use of spatially distributed time-integrated sediment sampling networks and distributed fine sediment modelling to inform catchment management

Under the EU Water Framework Directive, suspended sediment is omitted from environmental quality standards and compliance targets. This omission is partly explained by difficulties in assessing the complex dose-response of ecological communities. But equally, it is hindered by a lack of spatially distributed estimates of suspended sediment variability across catchments. In this paper, we demonstrate the inability of traditional, discrete sampling campaigns for assessing exposure to fine sediment. Sampling frequencies based on Environmental Quality Standard protocols, whilst reflecting typical manual sampling constraints, are unable to determine the magnitude of sediment exposure with an acceptable level of precision. Deviations from actual concentrations range between -35 and +20% based on the interquartile range of simulations. As an alternative, we assess the value of low-cost, suspended sediment sampling networks for quantifying suspended sediment transfer (SST). In this study of the 362 km² upland Esk catchment we observe that spatial patterns of sediment flux are consistent over the two year monitoring period across a network of 17 monitoring sites. This enables the key contributing sub-catchments of Butter Beck (SST: 1141 t km² yr^-1) and Glaisdale Beck (SST: 841 t km² yr^-1) to be identified. The time-integrated samplers offer a feasible alternative to traditional infrequent and discrete sampling approaches for assessing spatio-temporal changes in contamination. In conjunction with a spatially distributed diffuse pollution model (SCIMAP), time-integrated sediment sampling is an effective means of identifying critical sediment source areas in the catchment, which can better inform sediment management strategies for pollution prevention and control.

Y-Chromosome Structural Diversity in the Bonobo and Chimpanzee Lineages

The male-specific regions of primate Y-chromosomes (MSY) are enriched for multi-copy genes highly expressed in the testis. These genes are located in large repetitive sequences arranged as palindromes, inverted-, and tandem repeats termed amplicons. In humans, these genes have critical roles in male fertility and are essential for the production of sperm. The structure of human and chimpanzee amplicon sequences show remarkable difference relative to the remainder of the genome, a difference that may be the result of intense selective pressure on male fertility. Four subspecies of common chimpanzees have undergone extended periods of isolation and appear to be in the early process of subspeciation. A recent study found amplicons enriched for testis-expressed genes on the primate X-chromosome the target of hard selective sweeps, and male-fertility genes on the Y-chromosome may also be the targets of selection. However, little is understood about Y-chromosome amplicon diversity within and across chimpanzee populations. Here, we analyze nine common chimpanzee (representing three subspecies: Pan troglodytes schweinfurthii, Pan troglodytes ellioti, and Pan troglodytes verus) and two bonobo (Pan paniscus) male whole-genome sequences to assess Y ampliconic copy-number diversity across the Pan genus. We observe that the copy number of Y chromosome amplicons is variable among chimpanzees and bonobos, and identify several lineage-specific patterns, including variable copy number of azoospermia candidates RBMY and DAZ. We detect recurrent switchpoints of copy-number change along the ampliconic tracts across chimpanzee populations, which may be the result of localized genome instability or selective forces.

Resolving complex structural genomic rearrangements using a randomized approach

Complex chromosomal rearrangements are structural genomic alterations involving multiple instances of deletions, duplications, inversions, or translocations that co-occur either on the same chromosome or represent different overlapping events on homologous chromosomes. We present SVelter, an algorithm that identifies regions of the genome suspected to harbor a complex event and then resolves the structure by iteratively rearranging the local genome structure, in a randomized fashion, with each structure scored against characteristics of the observed sequencing data. SVelter is able to accurately reconstruct complex chromosomal rearrangements when compared to well-characterized genomes that have been deeply sequenced with both short and long reads.

The genomic landscape of polymorphic human nuclear mitochondrial insertions

The transfer of mitochondrial genetic material into the nuclear genomes of eukaryotes is a well-established phenomenon that has been previously limited to the study of static reference genomes. The recent advancement of high throughput sequencing has enabled an expanded exploration into the diversity of polymorphic nuclear mitochondrial insertions (NumtS) within human populations. We have developed an approach to discover and genotype novel Numt insertions using whole genome, paired-end sequencing data. We have applied this method to a thousand individuals in 20 populations from the 1000 Genomes Project and other datasets and identified 141 new sites of Numt insertions, extending our current knowledge of existing NumtS by almost 20%. We find that recent Numt insertions are derived from throughout the mitochondrial genome, including the D-loop, and have integration biases that differ in some respects from previous studies on older, fixed NumtS in the reference genome. We determined the complete inserted sequence for a subset of these events and have identified a number of nearly full-length mitochondrial genome insertions into nuclear chromosomes. We further define their age and origin of insertion and present an analysis of their potential impact to ongoing studies of mitochondrial heteroplasmy and disease.

Identification of p.A684V missense mutation in the WFS1 gene as a frequent cause of autosomal dominant optic atrophy and hearing impairment

Optic atrophy (OA) and sensorineural hearing loss (SNHL) are key abnormalities in several syndromes, including the recessively inherited Wolfram syndrome, caused by mutations in WFS1. In contrast, the association of autosomal dominant OA and SNHL without other phenotypic abnormalities is rare, and almost exclusively attributed to mutations in the Optic Atrophy-1 gene (OPA1), most commonly the p.R445H mutation. We present eight probands and their families from the US, Sweden, and UK with OA and SNHL, whom we analyzed for mutations in OPA1 and WFS1. Among these families, we found three heterozygous missense mutations in WFS1 segregating with OA and SNHL: p.A684V (six families), and two novel mutations, p.G780S and p.D797Y, all involving evolutionarily conserved amino acids and absent from 298 control chromosomes. Importantly, none of these families harbored the OPA1 p.R445H mutation. No mitochondrial DNA deletions were detected in muscle from one p.A684V patient analyzed. Finally, wolframin p.A684V mutant ectopically expressed in HEK cells showed reduced protein levels compared to wild-type wolframin, strongly indicating that the mutation is disease-causing. Our data support OA and SNHL as a phenotype caused by dominant mutations in WFS1 in these additional eight families. Importantly, our data provide the first evidence that a single, recurrent mutation in WFS1, p.A684V, may be a common cause of ADOA and SNHL, similar to the role played by the p.R445H mutation in OPA1. Our findings suggest that patients who are heterozygous for WFS1 missense mutations should be carefully clinically examined for OA and other manifestations of Wolfram syndrome.

Autosomal dominant progressive sensorineural hearing loss due to a novel mutation in the KCNQ4 gene

OBJECTIVE

To identify the genetic etiology in a family with autosomal dominant progressive sensorineural hearing loss.

DESIGN

Prospective molecular genetic research study.

SETTING

Academic genetic research laboratory.

PARTICIPANTS

Seventeen members of a family with dominant progressive nonsyndromic sensorineural hearing loss: 9 affected, 6 unaffected, and 2 spouses.

INTERVENTIONS

Clinical data from questionnaires, interviews, serial audiograms, and medical records; genetic data from genome-wide linkage analysis and candidate gene mutation analysis.

MAIN OUTCOME MEASURES

Symptoms, age at onset, serial audiometric data, and the presence or absence of a deafness-associated mutation.

RESULTS

Affected individuals in this family presented with autosomal dominant nonsyndromic high-frequency progressive sensorineural hearing loss, with age at onset ranging from 1 to 21 years. Genome-wide linkage analysis of single-nucleotide polymorphisms yielded evidence of linkage to an 18.9-Mb region on chromosome 1p34-p36, with a multipoint logarithm of odds score of 3.6. This interval contains a known deafness gene, KCNQ4, which underlies DNFA2 deafness. Sequencing of the 14 coding exons and intron-exon junctions of KCNQ4 revealed a novel heterozygous missense mutation, c.859G>C, p.Gly287Arg. The mutation disrupts the highly conserved GYG motif (glycine-tyrosine-glycine) of the phosphate-binding loop, hypothesized to be critical in maintaining pore structure and function. All 274 controls were negative for the mutation.

CONCLUSIONS

Autosomal dominant high-frequency hearing loss is genetically heterogeneous, and linkage analysis is an efficient means of identifying the etiology in larger families. Deafness in this family is caused by a novel mutation in KCNQ4.

Congenital cataracts in two siblings with Wolfram syndrome

BACKGROUND

Wolfram syndrome is characterized by optic atrophy, insulin dependent diabetes mellitus, diabetes insipidus and deafness. There are several other associated conditions reported in the literature, but congenital or early childhood cataracts are not among them.

MATERIALS AND METHODS

Observational case series with confirmatory genetic analysis.

RESULTS

A pair of siblings, followed over 17 years, who manifest congenital or early childhood cataracts, diabetes insipidus, diabetes mellitus, optic atrophy, and deafness. They are both compound heterozygotes for mutations (V415 deletion and A684V substitution) in the WFS1 gene. Their father has congenital sensorineural hearing loss and developed optic atrophy. He is heterozygous for A684V in WFS1.

CONCLUSIONS

Wolfram syndrome should be in the differential diagnosis of genetic syndromes associated with congenital and early childhood cataracts. Here, we report on a mother who is a phenotypically normal carrier of an autosomal recessive Wolfram syndrome gene, and a father who has some of the findings of the syndrome and carries a single mutation that appears to be responsible for his hearing loss and optic atrophy. Their 2 children are compound heterozygotes and manifest the full Wolfram syndrome, in addition to cataracts.

Effects of Manduca diuresin on neonates of the tobacco hornworm, Manduca sexta

Manduca diuresin (MD), a 30 amino acid peptide, was isolated and identified from the type IIb median neurosecretory cells in the brain of adult Manduca sexta (Blackburn et al., Biochem. Biophys. Res. Commun., 181, 927-932, 1991). The synthetic hormone stimulates water loss in the adult moth both in vitro using isolated Malphigian tubule preparations and in vivo using decapitated adult moths as bioassay animals (Blackburn and Ma, Arch. Insect Biochem. Physiol., 27, 3-10, 1994). In the present study, MD was shown to have adverse physiological effects on neonates of Manduca sexta when fed synthetic hormone-treated leaf discs. With a sustained diet of MD-treated discs, neonates exhibited symptoms of reduced food consumption and marginal weight increases. In addition to slowed growth and reduced developmental rates, a high percentage of these neonates failed to molt into second instar larvae and death usually followed shortly thereafter. This is a first report that a synthetic insect neuropeptide has shown detrimental effect on insect larvae by ingestion. The possibility of utilizing this peptide in plant protection is discussed.