Search

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion–base pair sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million base pairs of sequence containing 1956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

Among autistic individuals, a subphenotype with brain enlargement disproportionate to height (autism with disproportionate megalencephaly, or ASD-DM) seen at three years of age is associated with co-occurring intellectual disability and poorer prognoses later in life. However, many of the genes contributing to ASD-DM have yet to be delineated. In this study, we aim to identify additional ASD-DM associated genes to better define the genetic etiology of this subphenotype of autism. Here, we expand the previously studied sample size of ASD-DM individuals ten-fold by including probands from the Autism Phenome Project and Simons Simplex Collection, totaling 766 autistic individuals meeting the criteria for megalencephaly or macrocephaly and revealing 153 candidate ASD-DM genes harboringde novoprotein-impacting variants. Our findings include thirteen high confidence autism genes and seven genes previously associated with DM. Five impacted genes have previously been associated with both autism and DM, includingCHD8andPTEN. By performing functional network analysis, we also narrowed in on additional candidate genes, including one previously implicated in ASD-DM (PIK3CA) as well as 184 additional genes previously implicated in ASD or DM alone. Using zebrafish as a model, we performed CRISPR gene editing to generate knockout animals for seven of candidate genes and assessed head-size and induced seizure activity differences. From this analysis, we identified significant morphological changes in zebrafish knockout models of two genes,ythdf2andryr3. While zebrafish knockout mutants model haploinsufficiency of assayed genes, we identified ade novotandem duplication impactingYTHDF2in an ASD-DM proband. Therefore, we also transiently overexpressedYTHDF2by injection ofin vitrotranscribed human mRNA to simulate the patient-identified duplication. Following this, we observed increased head and brain size in theYTHDF2overexpression zebrafish matching that of the proband, providing a functional link betweenYTHDF2mutations and DM. Though discovery of additional mutations of ASD-DM candidate genes are required in order to fully elucidate the genetics associated with this severe form of autism, our study was able to provide support forYTHDF2as a novel putative ASD-DM gene.

Social monogamy is a reproductive strategy characterized by pair living and defense of a common territory. Pair bonding, sometimes displayed by monogamous species, is an affective construct that includes preference for a specific partner, distress upon separation, and the ability of the partner to buffer against stress. Many seahorse species show a monogamous social structure in the wild, but their pair bond has not been well studied. We examined the gene expression of lined seahorses (Hippocampus erectus) during and after the process of pairing in the laboratory as well as color change (luminance), a potential form of social communication and behavioral synchrony between pair mates. When a seahorse of either sex was interacting with its pair mate, their changes in luminance ("brightness") were correlated and larger than when interacting with an opposite-sex stranger. At the conclusion of testing, subjects were euthanized, RNA was extracted from whole brains and analyzed via RNA sequencing. Changes in gene expression in paired males versus those that were unpaired included processes governing metabolic activity, hormones and cilia. Perhaps most interesting is the overlap in gene expression change induced by pairing in both male seahorses and male prairie voles, including components of hormone systems regulating reproduction. Because of our limited sample size, we consider our results and interpretations to be preliminary, and prompts for further exploration. Future studies will expand upon these findings and investigate the neuroendocrine and genetic basis of these behaviors.

Emerging evidence links genes within human-specific segmental duplications (HSDs) to traits and diseases unique to our species. Strikingly, despite being nearly identical by sequence (>98.5%), paralogous HSD genes are differentially expressed across human cell and tissue types, though the underlying mechanisms have not been examined. We compared cross-tissue mRNA levels of 75 HSD genes from 30 families between humans and chimpanzees and found expression patterns consistent with relaxed selection on or neofunctionalization of derived paralogs. In general, ancestral paralogs exhibited greatest expression conservation with chimpanzee orthologs, though exceptions suggest certain derived paralogs may retain or supplant ancestral functions. Concordantly, analysis of long-read isoform sequencing data sets from diverse human tissues and cell lines found that about half of derived paralogs exhibited globally lower expression. To understand mechanisms underlying these differences, we leveraged data from human lymphoblastoid cell lines (LCLs) and found no relationship between paralogous expression divergence and post-transcriptional regulation, sequence divergence, or copy-number variation. Considering cis-regulation, we reanalyzed ENCODE data and recovered hundreds of previously unidentified candidate CREs in HSDs. We also generated large-insert ChIP-sequencing data for active chromatin features in an LCL to better distinguish paralogous regions. Some duplicated CREs were sufficient to drive differential reporter activity, suggesting they may contribute to divergent cis-regulation of paralogous genes. This work provides evidence that cis-regulatory divergence contributes to novel expression patterns of recent gene duplicates in humans.

Compared to its predecessors, the Telomere-to-Telomere CHM13 genome adds nearly 200 million base pairs of sequence, corrects thousands of structural errors, and unlocks the most complex regions of the human genome for clinical and functional study. We show how this reference universally improves read mapping and variant calling for 3202 and 17 globally diverse samples sequenced with short and long reads, respectively. We identify hundreds of thousands of variants per sample in previously unresolved regions, showcasing the promise of the T2T-CHM13 reference for evolutionary and biomedical discovery. Simultaneously, this reference eliminates tens of thousands of spurious variants per sample, including reduction of false positives in 269 medically relevant genes by up to a factor of 12. Because of these improvements in variant discovery coupled with population and functional genomic resources, T2T-CHM13 is positioned to replace GRCh38 as the prevailing reference for human genetics.

The enhanced cognitive abilities characterizing the human species result from specialized features of neurons and circuits, but the underlying molecular mechanisms remain largely unknown. Here we report that the hominid-specific geneLRRC37Bencodes a novel receptor expressed in a subset of human cortical pyramidal neurons (CPNs). LRRC37B protein localizes at the axon initial segment (AIS), the specialized domain triggering action potentials.LRRC37Bectopic expression in mouse CPNsin vivoleads to reduced intrinsic excitability, a distinctive feature of some classes of human CPNs. At the molecular level, LRRC37B acts as a receptor for the secreted ligand FGF13A and interacts with the voltage gated sodium channel (VGSC) beta subunit SCN1B, thereby inhibiting the channel function of VGSC, specifically at the AIS. Electrophysiological recordings in adult human cortical slices reveals that endogenous expression of LRRC37B in human CPNs reduces neuronal excitability.LRRC37Bthus acts as a species-specific modifier of human cortical neuron function, with important implications for human brain evolution and diseases.

Compared to its predecessors, the Telomere-to-Telomere CHM13 genome adds nearly 200 Mbp of sequence, corrects thousands of structural errors, and unlocks the most complex regions of the human genome to clinical and functional study. Here we demonstrate how the new reference universally improves read mapping and variant calling for 3,202 and 17 globally diverse samples sequenced with short and long reads, respectively. We identify hundreds of thousands of novel variants per sample—a new frontier for evolutionary and biomedical discovery. Simultaneously, the new reference eliminates tens of thousands of spurious variants per sample, including up to 12-fold reduction of false positives in 269 medically relevant genes. The vast improvement in variant discovery coupled with population and functional genomic resources position T2T-CHM13 to replace GRCh38 as the prevailing reference for human genetics.One Sentence SummaryThe T2T-CHM13 reference genome universally improves the analysis of human genetic variation.

INTRODUCTION: To faithfully distribute genetic material to daughter cells during cell division, spindle fibers must couple to DNA by means of a structure called the kinetochore, which assembles at each chromosome’s centromere. Human centromeres are located within large arrays of tandemly repeated DNA sequences known as alpha satellite (αSat), which often span millions of base pairs on each chromosome. Arrays of αSat are frequently surrounded by other types of tandem satellite repeats, which have poorly understood functions, along with nonrepetitive sequences, including transcribed genes. Previous genome sequencing efforts have been unable to generate complete assemblies of satellite-rich regions because of their scale and repetitive nature, limiting the ability to study their organization, variation, and function. RATIONALE: Pericentromeric and centromeric (peri/centromeric) satellite DNA sequences have remained almost entirely missing from the assembled human reference genome for the past 20 years. Using a complete, telomere-to-telomere (T2T) assembly of a human genome, we developed and deployed tailored computational approaches to reveal the organization and evolutionary patterns of these satellite arrays at both large and small length scales. We also performed experiments to map precisely which αSat repeats interact with kinetochore proteins. Last, we compared peri/centromeric regions among multiple individuals to understand how these sequences vary across diverse genetic backgrounds. RESULTS: Satellite repeats constitute 6.2% of the T2T-CHM13 genome assembly, with αSat representing the single largest component (2.8% of the genome). By studying the sequence relationships of αSat repeats in detail across each centromere, we found genome-wide evidence that human centromeres evolve through “layered expansions.” Specifically, distinct repetitive variants arise within each centromeric region and expand through mechanisms that resemble successive tandem duplications, whereas older flanking sequences shrink and diverge over time. We also revealed that the most recently expanded repeats within each αSat array are more likely to interact with the inner kinetochore protein Centromere Protein A (CENP-A), which coincides with regions of reduced CpG methylation. This suggests a strong relationship between local satellite repeat expansion, kinetochore positioning, and DNA hypomethylation. Furthermore, we uncovered large and unexpected structural rearrangements that affect multiple satellite repeat types, including active centromeric αSat arrays. Last, by comparing sequence information from nearly 1600 individuals’ Xchromosomes, we observed that individuals with recent African ancestry possess the greatest genetic diversity in the region surrounding the centromere, which sometimes contains a predominantly African αSat sequence variant. CONCLUSION: The genetic and epigenetic properties of centromeres are closely interwoven through evolution. These findings raise important questions about the specific molecular mechanisms responsible for the relationship between inner kinetochore proteins, DNA hypomethylation, and layered αSat expansions. Even more questions remain about the function and evolution of non-αSat repeats. To begin answering these questions, we have produced a comprehensive encyclopedia of peri/centromeric sequences in a human genome, and we demonstrated how these regions can be studied with modern genomic tools. Our work also illuminates the rich genetic variation hidden within these formerly missing regions of the genome, which may contribute to health and disease. This unexplored variation underlines the need for more T2T human genome assemblies from genetically diverse individuals.

In 2001, Celera Genomics and the International Human Genome Sequencing Consortium published their initial drafts of the human genome, which revolutionized the field of genomics. While these drafts and the updates that followed effectively covered the euchromatic fraction of the genome, the heterochromatin and many other complex regions were left unfinished or erroneous. Addressing this remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium has finished the first truly complete 3.055 billion base pair (bp) sequence of a human genome, representing the largest improvement to the human reference genome since its initial release. The new T2T-CHM13 reference includes gapless assemblies for all 22 autosomes plus Chromosome X, corrects numerous errors, and introduces nearly 200 million bp of novel sequence containing 2,226 paralogous gene copies, 115 of which are predicted to be protein coding. The newly completed regions include all centromeric satellite arrays and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies for the first time.

Tracking and quantifying the abundance and location of cells in the developing brain is essential in neuroscience research, enabling a greater understanding of mechanisms underlying nervous system morphogenesis. Widely used experimental methods to quantify cells labeled with fluorescent markers—such as immunohistochemistry, in situ hybridization, and expression of transgenes via stable lines or transient in utero electroporations—depend upon accurate and consistent quantification of images. Current methods to quantify fluorescently-labeled cells rely on labor-intensive manual counting approaches, such as the Fiji plugin Cell Counter, which requires custom macros to enable higher-throughput analyses. Here, we present RapID Cell Counter, a semi-automated cell-counting tool with an easy-to-implement graphical user interface, that facilitates quick and consistent quantifications of cell density within user-defined boundaries that can be divided into equally-partitioned segments. Compared to the standard manual counting approach, we show that RapID matched accuracy and consistency, and only required ∼10% of user time relative to manual counting methods, when quantifying the distribution of fluorescently-labeled neurons in mouse in utero electroporation experiments. Using RapID, we recapitulated previously published work focusing on two genes, Srgap2 and Cul5, important for projection neuron migration in the neocortex and used it to quantify projection neuron displacement in a mouse knockout model of Rbx2. Moreover, RapID is capable of quantifying other cell types in the brain with complex cell morphologies, including astrocytes and dopaminergic neurons. We propose RapID as an efficient method for neuroscience researchers to process fluorescently-labeled brain images in a consistent, accurate, and high-throughput manner. SIGNIFICANCE STATEMENT Most studies in neuroscience rely on imaging to elucidate key neurodevelopmental processes including cell migration and proliferation. Many imaging techniques, including in utero electroporation and immunohistochemistry, produce multitudes of images that require accurate quantification, often via labor-intensive manual counting by multiple individuals that may delay follow-up experiments. To address this problem, we developed RapID, an efficient and semi-automated cell counting software platform, that reduces the time spent to 1/10th compared to one of the most popular quantification methods used for imaging studies today. RapID is flexible across imaging platforms and easily implemented through a graphical user interface.

Emerging evidence links genes within human-specific segmental duplications (HSDs) to traits and diseases unique to our species. Strikingly, despite being nearly identical by sequence (>98.5%), paralogous HSD genes are differentially expressed across human cell and tissue types, though the underlying mechanisms have not been examined. We compared cross-tissue mRNA levels of 75 HSD genes from 30 families between humans and chimpanzees and found expression patterns consistent with pseudo- or neofunctionalization. In general, ancestral paralogs exhibited greatest expression conservation with chimpanzee orthologs, though exceptions suggest certain derived paralogs may retain or supplant ancestral functions. Concordantly, analysis of long-read isoform sequencing datasets from diverse human tissues and cell lines found that about half of derived paralogs exhibited globally lower expression. To understand mechanisms underlying these differences, we leveraged data from human lymphoblastoid cell lines (LCLs) and found no relationship between paralogous expression divergence and post- transcriptional regulation, sequence divergence, or copy number variation. Considering cis-regulation, we reanalyzed ENCODE data and recovered hundreds of previously unidentified candidate CREs in HSDs. We also generated large-insert ChIP-sequencing data for active chromatin features in an LCL to better distinguish paralogous regions. Some duplicated CREs were sufficient to drive differential reporter activity, suggesting they may contribute to divergent cis-regulation of paralogous genes. This work provides evidence that cis-regulatory divergence contributes to novel expression patterns of recent gene duplicates in humans.

BackgroundZebrafish have practical features that make them a useful model for higher-throughput tests of gene function using CRISPR/Cas9 editing to create ‘knockout’ models. In particular, the use of G0 mosaic mutants has potential to increase throughput of functional studies significantly but may suffer from transient effects of introducing Cas9 via microinjection. Further, a large number of computational and empirical tools exist to design CRISPR assays but often produce varied predictions across methods leaving uncertainty in choosing an optimal approach for zebrafish studies.MethodsTo systematically assess accuracy of tool predictions of on- and off-target gene editing, we subjected zebrafish embryos to CRISPR/Cas9 with 50 different guide RNAs (gRNAs) targeting 14 genes. We also investigate potential confounders of G0-based CRISPR screens by screening control embryos for spurious mutations and altered gene expression.ResultsWe compared our experimental in vivo editing efficiencies in mosaic G0 embryos with those predicted by eight commonly used gRNA design tools and found large discrepancies between methods. Assessing off-target mutations (predicted in silico and in vitro) found that the majority of tested loci had low in vivo frequencies (0 mosaic CRISPR knockout fish, we generated an RNA-seq dataset of various control larvae at 5 days post fertilization. While we found no evidence of spontaneous somatic mutations of injected larvae, we did identify several hundred differentially-expressed genes with high variability between injection types. Network analyses of shared differentially-expressed genes in the ‘mock’ injected larvae implicated a number of key regulators of common metabolic pathways, and gene-ontology analysis revealed connections with response to wounding and cytoskeleton organization, highlighting a potentially lasting effect from the microinjection process that requires further investigation.ConclusionOverall, our results provide a valuable resource for the zebrafish community for the design and execution of CRISPR/Cas9 experiments.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes., High-coverage, ultra-long-read nanopore sequencing is used to create a new human genome assembly that improves on the coverage and accuracy of the current reference (GRCh38) and includes the gap-free, telomere-to-telomere sequence of the X chromosome.

In recent years zebrafish have become commonly used as a model for studying human traits and disorders. Their small size, high fecundity, and rapid development allow for more high-throughput experiments compared to other vertebrate models. Given that zebrafish share >70% gene homologs with humans and their genomes can be readily edited using highly efficient CRISPR methods, we are now able to rapidly generate mutations impacting practically any gene of interest. Unfortunately, our ability to phenotype mutant larvae has not kept pace. To address this challenge, we have developed a protocol that obtains multiple phenotypic measurements from individual zebrafish larvae in an automated and parallel fashion, including morphological features (i.e., body length, eye area, and head size) and movement/behavior. By assaying wild-type zebrafish in a variety of conditions, we determined optimal parameters that avoid significant developmental defects or physical damage; these include morphological imaging of larvae at two time points (3 days post fertilization (dpf) and 5 dpf) coupled with motion tracking of behavior at 5 dpf. As a proof-of-principle, we tested our approach on two novel CRISPR-generated mutant zebrafish lines carrying predicted null-alleles of syngap1b and slc7a5, orthologs to two human genes implicated in autism-spectrum disorder, intellectual disability, and epilepsy. Using our optimized high-throughput phenotyping protocol, we recapitulated previously published results from mouse and zebrafish models of these candidate genes. In summary, we describe a rapid parallel pipeline to characterize morphological and behavioral features of individual larvae in a robust and consistent fashion, thereby improving our ability to better identify genes important in human traits and disorders.AUTHOR SUMMARYZebrafish (Danio rerio) are a well-established model organism for the study of neurodevelopmental disorders. Due to their small size, fast reproduction, and genetic homology with humans, zebrafish have been widely used for characterizing and screening candidate genes for many disorders, including autism-spectrum disorder, intellectual disability, and epilepsy. Although several studies have described the use of high-throughput morphological and behavioral assays, few combine multiple assays in a single zebrafish larva. Here, we optimized a platform to characterize morphometric features at two developmental time points in addition to behavioral traits of zebrafish larvae. We then used this approach to characterize two autism candidate genes (SYNGAP1 and SLC7A5) in two CRISPR-generated zebrafish null mutant models we developed in house. These data recapitulate previously published results related to enhanced seizure activity, while identifying additional defects not previously reported. We propose that our phenotyping platform represents a feasible method for maximizing the use of single zebrafish larvae in the characterization of additional mutants relevant to neurodevelopmental disorders.

Recent efforts to comprehensively characterize great ape genetic diversity using short-read sequencing and single-nucleotide variants have led to important discoveries related to selection within species, demographic history, and lineage-specific traits. Structural variants (SVs), including deletions and inversions, comprise a larger proportion of genetic differences between and within species, making them an important yet understudied source of trait divergence. Here, we used a combination of long-read and -range sequencing approaches to characterize the structural variant landscape of two additional Pan troglodytes verus individuals, one of whom carries 13% admixture from Pan troglodytes troglodytes. We performed optical mapping of both individuals followed by nanopore sequencing of one individual. Filtering for larger variants (&gt, 10 kbp) and combined with genotyping of SVs using short-read data from the Great Ape Genome Project, we identified 425 deletions and 59 inversions, of which 88 and 36, respectively, were novel. Compared with gene expression in humans, we found a significant enrichment of chimpanzee genes with differential expression in lymphoblastoid cell lines and induced pluripotent stem cells, both within deletions and near inversion breakpoints. We examined chromatin-conformation maps from human and chimpanzee using these same cell types and observed alterations in genomic interactions at SV breakpoints. Finally, we focused on 56 genes impacted by SVs in &gt, 90% of chimpanzees and absent in humans and gorillas, which may contribute to chimpanzee-specific features. Sequencing a greater set of individuals from diverse subspecies will be critical to establish the complete landscape of genetic variation in chimpanzees.

Rhesus macaque is an Old World monkey that shared a common ancestor with human ∼25 Myr ago and is an important animal model for human disease studies. A deep understanding of its genetics is therefore required for both biomedical and evolutionary studies. Among structural variants, inversions represent a driving force in speciation and play an important role in disease predisposition. Here we generated a genome-wide map of inversions between human and macaque, combining single-cell strand sequencing with cytogenetics. We identified 375 total inversions between 859 bp and 92 Mbp, increasing by eightfold the number of previously reported inversions. Among these, 19 inversions flanked by segmental duplications overlap with recurrent copy number variants associated with neurocognitive disorders. Evolutionary analyses show that in 17 out of 19 cases, the Hominidae orientation of these disease-associated regions is always derived. This suggests that duplicated sequences likely played a fundamental role in generating inversions in humans and great apes, creating architectures that nowadays predispose these regions to disease-associated genetic instability. Finally, we identified 861 genes mapping at 156 inversions breakpoints, with some showing evidence of differential expression in human and macaque cell lines, thus highlighting candidates that might have contributed to the evolution of species-specific features. This study depicts the most accurate fine-scale map of inversions between human and macaque using a two-pronged integrative approach, such as single-cell strand sequencing and cytogenetics, and represents a valuable resource toward understanding of the biology and evolution of primate species.

SUMMARY Atrophy of neurons in the prefrontal cortex (PFC) plays a key role in the pathophysiology of depression and related disorders. The ability to promote both structural and functional plasticity in the PFC has been hypothesized to underlie the fast-acting antidepressant properties of the dissociative anesthetic ketamine. Here, we report that, like ketamine, serotonergic psychedelics are capable of robustly increasing neuritogenesis and/or spinogenesis both in vitro and in vivo. These changes in neuronal structure are accompanied by increased synapse number and function, as measured by fluorescence microscopy and electrophysiology. The structural changes induced by psychedelics appear to result from stimulation of the TrkB, mTOR, and 5-HT2A signaling pathways and could possibly explain the clinical effectiveness of these compounds. Our results underscore the therapeutic potential of psychedelics and, importantly, identify several lead scaffolds for medicinal chemistry efforts focused on developing plasticity-promoting compounds as safe, effective, and fast-acting treatments for depression and related disorders., In Brief Ly et al. demonstrate that psychedelic compounds such as LSD, DMT, and DOI increase dendritic arbor complexity, promote dendritic spine growth, and stimulate synapse formation. These cellular effects are similar to those produced by the fast-acting antidepressant ketamine and highlight the potential of psychedelics for treating depression and related disorders.

Background Fear of autism has led to a decline in childhood-immunization uptake and to a resurgence of preventable infectious diseases. Identifying characteristics of parents who believe in a causal role of vaccines for autism spectrum disorder (ASD) in their child may help targeting educational activities and improve adherence to the immunization schedule. Objectives To compare caregivers of children with ASD who agree or disagree that vaccines play an etiological role in autism for 1) socio-demographics characteristics and 2) developmental and clinical profiles of their children. Methods Data from 16,525 participants with ASD under age 18 were obtained from SPARK, a national research cohort started in 2016. Caregivers completed questionnaires at registration that included questions on beliefs about the etiologic role of childhood immunizations and other factors in ASD. Data were available about family socio-demographic characteristics, first symptoms of autism, developmental regression, co-occurring psychiatric disorders, seizures, and current levels of functioning. Results Participants with ASD were 80.4% male with a mean age of 8.1 years (SD = 4.1). Overall, 16.5% of caregivers endorsed immunizations as perceived causes of autism. Compared to caregivers who disagreed with vaccines as a cause for ASD, those who believed in vaccine causation came disproportionately from ethnic minority, less educated, and less wealthy backgrounds. More often their children had experienced developmental regression involving language and other skills, were diagnosed earlier, had lost skills during the second year of life, and had worse language, adaptive, and cognitive outcomes. Conclusion One in six caregivers who participate in a national research cohort believe that child immunizations could be a cause of autism in their child. Parent social background (non-White, less educated) and child developmental features (regression in second year, poorer language skills, and worse adaptive outcomes) index caregivers who are more likely to harbor these beliefs and could benefit from targeted educational activities.

Ambient temperature is a critical environmental factor for all living organisms. It was likely an important selective force as modern humans recently colonized temperate and cold Eurasian environments. Nevertheless, as of yet we have limited evidence of local adaptation to ambient temperature in populations from those environments. To shed light on this question, we exploit the fact that humans are a cosmopolitan species that inhabit territories under a wide range of temperatures. Focusing on cold perception–which is central to thermoregulation and survival in cold environments–we show evidence of recent local adaptation on TRPM8. This gene encodes for a cation channel that is, to date, the only temperature receptor known to mediate an endogenous response to moderate cold. The upstream variant rs10166942 shows extreme population differentiation, with frequencies that range from 5% in Nigeria to 88% in Finland (placing this SNP in the 0.02% tail of the FST empirical distribution). When all populations are jointly analyzed, allele frequencies correlate with latitude and temperature beyond what can be explained by shared ancestry and population substructure. Using a Bayesian approach, we infer that the allele originated and evolved neutrally in Africa, while positive selection raised its frequency to different degrees in Eurasian populations, resulting in allele frequencies that follow a latitudinal cline. We infer strong positive selection, in agreement with ancient DNA showing high frequency of the allele in Europe 3,000 to 8,000 years ago. rs10166942 is important phenotypically because its ancestral allele is protective of migraine. This debilitating disorder varies in prevalence across human populations, with highest prevalence in individuals of European descent–precisely the population with the highest frequency of rs10166942 derived allele. We thus hypothesize that local adaptation on previously neutral standing variation may have contributed to the genetic differences that exist in the prevalence of migraine among human populations today., Author summary Some human populations were likely under strong pressure to adapt biologically to cold climates during their colonization of non-African territories in the last 50,000 years. Such putative adaptations required genetic variation in genes that could mediate adaptive responses to cold. TRPM8 is potentially one such gene, being the only known receptor for the sensation of moderate cold temperature. We show that a likely regulatory genetic variant nearby TRPM8 has several signatures of positive selection raising its frequency in Eurasian populations during the last 25,000 years. While the genetic variant was and is rare in Africa, it is now common outside of Africa, with frequencies that strongly correlate with latitude and are highest in northern European populations. Interestingly, this same genetic variant has previously been strongly associated with migraine. This suggests that adaptation to cold has potentially contributed to the variation in migraine prevalence that exists among human groups today.

© 2017 Macmillan Publishers Limited, part of Springer Nature. Segmental duplications contribute to human evolution, adaptation and genomic instability but are often poorly characterized. We investigate the evolution, genetic variation and coding potential of human-specific segmental duplications (HSDs). We identify 218 HSDs based on analysis of 322 deeply sequenced archaic and contemporary hominid genomes. We sequence 550 human and nonhuman primate genomic clones to reconstruct the evolution of the largest, most complex regions with protein-coding potential (N = 80 genes from 33 gene families). We show that HSDs are non-randomly organized, associate preferentially with ancestral ape duplications termed 'core duplicons' and evolved primarily in an interspersed inverted orientation. In addition to Homo sapiens-specific gene expansions (such as TCAF1/TCAF2), we highlight ten gene families (for example, ARHGAP11B and SRGAP2C) where copy number never returns to the ancestral state, there is evidence of mRNA splicing and no common gene-disruptive mutations are observed in the general population. Such duplicates are candidates for the evolution of human-specific adaptive traits.

Obtaining high-quality sequence continuity of complex regions of recent segmental duplication remains one of the major challenges of finishing genome assemblies. In the human and mouse genomes, this was achieved by targeting large-insert clones using costly and laborious capillary-based sequencing approaches. Sanger shotgun sequencing of clone inserts, however, has now been largely abandoned, leaving most of these regions unresolved in newer genome assemblies generated primarily by next-generation sequencing hybrid approaches. Here we show that it is possible to resolve regions that are complex in a genome-wide context but simple in isolation for a fraction of the time and cost of traditional methods using long-read single molecule, real-time (SMRT) sequencing and assembly technology from Pacific Biosciences (PacBio). We sequenced and assembled BAC clones corresponding to a 1.3-Mbp complex region of chromosome 17q21.31, demonstrating 99.994% identity to Sanger assemblies of the same clones. We targeted 44 differences using Illumina sequencing and find that PacBio and Sanger assemblies share a comparable number of validated variants, albeit with different sequence context biases. Finally, we targeted a poorly assembled 766-kbp duplicated region of the chimpanzee genome and resolved the structure and organization for a fraction of the cost and time of traditional finishing approaches. Our data suggest a straightforward path for upgrading genomes to a higher quality finished state. © 2014 Huddleston et al.

© 2016 Elsevier Ltd Duplications are the primary force by which new gene functions arise and provide a substrate for large-scale structural variation. Analysis of thousands of genomes shows that humans and great apes have more genetic differences in content and structure over recent segmental duplications than any other euchromatic region. Novel human-specific duplicated genes, ARHGAP11B and SRGAP2C, have recently been described with a potential role in neocortical expansion and increased neuronal spine density. Large segmental duplications and the structural variants they promote are also frequently stratified between human populations with a subset being subjected to positive selection. The impact of recent duplications on human evolution and adaptation is only beginning to be realized as new technologies enhance their discovery and accurate genotyping.

Rare copy-number variants (CNVs) have been implicated in autism and intellectual disability. These variants are large and affect many genes but lack clear specificity toward autism as opposed to developmental-delay phenotypes. We exploited the repeat architecture of the genome to target segmental duplication-mediated rearrangement hotspots (n = 120, median size 1.78 Mbp, range 240 kbp to 13 Mbp) and smaller hotspots flanked by repetitive sequence (n = 1,247, median size 79 kbp, range 3–96 kbp) in 2,588 autistic individuals from simplex and multiplex families and in 580 controls. Our analysis identified several recurrent large hotspot events, including association with 1q21 duplications, which are more likely to be identified in individuals with autism than in those with developmental delay (p = 0.01; OR = 2.7). Within larger hotspots, we also identified smaller atypical CNVs that implicated CHD1L and ACACA for the 1q21 and 17q12 deletions, respectively. Our analysis, however, suggested no overall increase in the burden of smaller hotspots in autistic individuals as compared to controls. By focusing on gene-disruptive events, we identified recurrent CNVs, including DPP10, PLCB1, TRPM1, NRXN1, FHIT, and HYDIN, that are enriched in autism. We found that as the size of deletions increases, nonverbal IQ significantly decreases, but there is no impact on autism severity; and as the size of duplications increases, autism severity significantly increases but nonverbal IQ is not affected. The absence of an increased burden of smaller CNVs in individuals with autism and the failure of most large hotspots to refine to single genes is consistent with a model where imbalance of multiple genes contributes to a disease state.

Background Gene innovation by duplication is a fundamental evolutionary process but is difficult to study in humans due to the large size, high sequence identity, and mosaic nature of segmental duplication blocks. The human-specific gene hydrocephalus-inducing 2, HYDIN2, was generated by a 364 kbp duplication of 79 internal exons of the large ciliary gene HYDIN from chromosome 16q22.2 to chromosome 1q21.1. Because the HYDIN2 locus lacks the ancestral promoter and seven terminal exons of the progenitor gene, we sought to characterize transcription at this locus by coupling reverse transcription polymerase chain reaction and long-read sequencing. Results 5' RACE indicates a transcription start site for HYDIN2 outside of the duplication and we observe fusion transcripts spanning both the 5' and 3' breakpoints. We observe extensive splicing diversity leading to the formation of altered open reading frames (ORFs) that appear to be under relaxed selection. We show that HYDIN2 adopted a new promoter that drives an altered pattern of expression, with highest levels in neural tissues. We estimate that the HYDIN duplication occurred ~3.2 million years ago and find that it is nearly fixed (99.9%) for diploid copy number in contemporary humans. Examination of 73 chromosome 1q21 rearrangement patients reveals that HYDIN2 is deleted or duplicated in most cases. Conclusions Together, these data support a model of rapid gene innovation by fusion of incomplete segmental duplications, altered tissue expression, and potential subfunctionalization or neofunctionalization of HYDIN2 early in the evolution of the Homo lineage. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1163-9) contains supplementary material, which is available to authorized users.

Objective: Epidemiological data have suggested maternal infection and fever to be associated with increased risk of autism spectrum disorder (ASD). Animal studies show that gestational infections perturb fetal brain development and result in offspring with the core features of autism and have demon- strated that behavioral effects of maternal immune activation are dependent on genetic susceptibility. The goal of this study was to explore the impact of ASD-associated copy number variants (CNVs) and prenatal maternal infection on clinical severity of ASD within a dataset of prenatal history and complete genetic and phenotypic findings. Methods: We analyzed data from the Simons Simplex Collection sample including 1971 children with a diagnosis of ASD aged 4 to 18 years who underwent array comparative genomic hybridization screening. Information on infection and febrile episodes during pregnancy was collected through parent interview. ASD severity was clinically measured through parent-reported interview and questionnaires. Results: We found significant interactive effects between the presence of CNVs and maternal infection during pregnancy on autistic symptomatology, such that individuals with CNVs and history of maternal infection demonstrated increased rates of social communicative impairments and repetitive/restricted behaviors. In contrast, no significant interactions were found between presence of CNVs and prenatal infections on cog- nitive and adaptive functioning of individuals with ASD. Conclusions: Our findings support a gene- environment interaction model of autism impairment, in that individuals with ASD-associated CNVs are more susceptible to the effects of maternal infection and febrile episodes in pregnancy on behavioral out- comes and suggest that these effects are specific to ASD rather than to global neurodevelopment. (J Dev Behav Pediatr 36:61-67, 2015) Index terms: autism, autism spectrum disorders, fever, gene-environment, infection, pregnancy.

© 2015 Macmillan Publishers Limited. All rights reserved. The human genome is arguably the most complete mammalian reference assembly, yet more than 160 euchromatic gaps remain and aspects of its structural variation remain poorly understood ten years after its completion. To identify missing sequence and genetic variation, here we sequence and analyse a haploid human genome (CHM1) using single-molecule, real-time DNA sequencing. We close or extend 55% of the remaining interstitial gaps in the human GRCh37 reference genome - 78% of which carried long runs of degenerate short tandem repeats, often several kilobases in length, embedded within (G+C)-rich genomic regions. We resolve the complete sequence of 26,079 euchromatic structural variants at the base-pair level, including inversions, complex insertions and long tracts of tandem repeats. Most have not been previously reported, with the greatest increases in sensitivity occurring for events less than 5 kilobases in size. Compared to the human reference, we find a significant insertional bias (3:1) in regions corresponding to complex insertions and long short tandem repeats. Our results suggest a greater complexity of the human genome in the form of variation of longer and more complex repetitive DNA that can now be largely resolved with the application of this longer-read sequencing technology.

The DYX2 locus on chromosome 6p22.2 is the most replicated region of linkage to developmental dyslexia (DD). Two candidate genes within this region have recently been implicated in the disorder: KIAA0319 and DCDC2. Variants within DCDC2 have shown association with DD in a US and a German sample. However, when we genotyped these specific variants in two large, independent UK samples, we obtained only weak, inconsistent evidence for their involvement in DD. Having previously found evidence that variation in the KIAA0319 gene confers susceptibility to DD, we sought to refine this genetic association by genotyping 36 additional SNPs in the gene. Nine SNPs, predominantly clustered around the first exon, showed the most significant association with DD in one or both UK samples, including rs3212236 in the 5' flanking region (P = 0.00003) and rs761100 in intron 1 (P = 0.0004). We have thus refined the region of association with developmental dyslexia to putative regulatory sequences around the first exon of the KIAA0319 gene, supporting the presence of functional mutations that could affect gene expression. Our data also suggests a possible interaction between KIAA0319 and DCDC2, which requires further testing.

Recurrent deletions of chromosome 15q13.3 associate with intellectual disability, schizophrenia, autism and epilepsy. To gain insight into the instability of this region, we sequenced it in affected individuals, normal individuals and nonhuman primates. We discovered five structural configurations of the human chromosome 15q13.3 region ranging in size from 2 to 3 Mb. These configurations arose recently (∼0.5-0.9 million years ago) as a result of human-specific expansions of segmental duplications and two independent inversion events. All inversion breakpoints map near GOLGA8 core duplicons-a ∼14-kb primate-specific chromosome 15 repeat that became organized into larger palindromic structures. GOLGA8-flanked palindromes also demarcate the breakpoints of recurrent 15q13.3 microdeletions, the expansion of chromosome 15 segmental duplications in the human lineage and independent structural changes in apes. The significant clustering (P = 0.002) of breakpoints provides mechanistic evidence for the role of this core duplicon and its palindromic architecture in promoting the evolutionary and disease-related instability of chromosome 15.

We report an algorithm to detect structural variation and indels from 1 base pair to 1 megabase pair within exome sequence datasets. Splitread uses one-end anchored placements to cluster the mappings of subsequences of unanchored ends to identify the size, content and location of variants with good specificity and high sensitivity. The algorithm discovers indels, structural variants, de novo events and copy-number polymorphic processed pseudogenes missed by other methods.

Summary Gene duplication is an important source of phenotypic change and adaptive evolution. We leverage a haploid hydatidiform mole to identify highly identical sequences missing from the reference genome, confirming that the cortical development gene Slit-Robo Rho GTPase-activating protein 2 ( SRGAP2 ) duplicated three times exclusively in humans. We show that the promoter and first nine exons of SRGAP2 duplicated from 1q32.1 ( SRGAP2A ) to 1q21.1 ( SRGAP2B ) ∼3.4 million years ago (mya). Two larger duplications later copied SRGAP2B to chromosome 1p12 ( SRGAP2C ) and to proximal 1q21.1 ( SRGAP2D ) ∼2.4 and ∼1 mya, respectively. Sequence and expression analyses show that SRGAP2C is the most likely duplicate to encode a functional protein and is among the most fixed human-specific duplicate genes. Our data suggest a mechanism where incomplete duplication created a novel gene function—antagonizing parental SRGAP2 function—immediately "at birth" 2–3 mya, which is a time corresponding to the transition from Australopithecus to Homo and the beginning of neocortex expansion.

A remarkable characteristic of the human major histocompatibility complex (MHC) is its extreme genetic diversity, which is maintained by balancing selection. In fact, the MHC complex remains one of the best-known examples of natural selection in humans, with well-established genetic signatures and biological mechanisms for the action of selection. Here, we present genetic and functional evidence that another gene with a fundamental role in MHC class I presentation, endoplasmic reticulum aminopeptidase 2 (ERAP2), has also evolved under balancing selection and contains a variant that affects antigen presentation. Specifically, genetic analyses of six human populations revealed strong and consistent signatures of balancing selection affecting ERAP2. This selection maintains two highly differentiated haplotypes (Haplotype A and Haplotype B), with frequencies 0.44 and 0.56, respectively. We found that ERAP2 expressed from Haplotype B undergoes differential splicing and encodes a truncated protein, leading to nonsense-mediated decay of the mRNA. To investigate the consequences of ERAP2 deficiency on MHC presentation, we correlated surface MHC class I expression with ERAP2 genotypes in primary lymphocytes. Haplotype B homozygotes had lower levels of MHC class I expressed on the surface of B cells, suggesting that naturally occurring ERAP2 deficiency affects MHC presentation and immune response. Interestingly, an ERAP2 paralog, endoplasmic reticulum aminopeptidase 1 (ERAP1), also shows genetic signatures of balancing selection. Together, our findings link the genetic signatures of selection with an effect on splicing and a cellular phenotype. Although the precise selective pressure that maintains polymorphism is unknown, the demonstrated differences between the ERAP2 splice forms provide important insights into the potential mechanism for the action of selection., Author Summary It has long been known that the extremely high levels of genetic diversity present in the major histocompatibility locus (MHC) are due to balancing selection, a type of natural selection that maintains advantageous genetic diversity in populations. The MHC encodes for molecules required for a type of antigen presentation that mediates detection of infected and cancerous cells by the immune system; the genetic diversity of the MHC thus ensures an adequate response to the wide variety of pathogens that humans encounter. Here, we show that other genes involved in the same antigen-presentation pathway are also subject to balancing selection in humans. Specifically, we show that balancing selection acts to maintain two forms of the endoplasmic reticulum aminopeptidase 2 gene (ERAP2), which encodes a protein also involved in antigen presentation. Although the two ERAP2 forms are present in a similar frequency (close to 0.5), they are associated with differences with respect to the levels of MHC molecules on the cell surface of immune cells. In summary, our findings show that natural selection maintains variants of ERAP2 that affect immune surveillance; they also establish ERAP2 as one of the few examples of balancing selection in humans where the selected variant, its functional consequences, and its influence in interpersonal diversity are known.

Background Myelin protein zero (MPZ) is a critical structural component of myelin in the peripheral nervous system. The MPZ gene is regulated, in part, by the transcription factors SOX10 and EGR2. Mutations in MPZ, SOX10, and EGR2 have been implicated in demyelinating peripheral neuropathies, suggesting that components of this transcriptional network are candidates for harboring disease-causing mutations (or otherwise functional variants) that affect MPZ expression. Methodology We utilized a combination of multi-species sequence comparisons, transcription factor-binding site predictions, targeted human DNA re-sequencing, and in vitro and in vivo enhancer assays to study human non-coding MPZ variants. Principal Findings Our efforts revealed a variant within the first intron of MPZ that resides within a previously described SOX10 binding site is associated with decreased enhancer activity, and alters binding of nuclear proteins. Additionally, the genomic segment harboring this variant directs tissue-relevant reporter gene expression in zebrafish. Conclusions This is the first reported MPZ variant within a cis-acting transcriptional regulatory element. While we were unable to implicate this variant in disease onset, our data suggests that similar non-coding sequences should be screened for mutations in patients with neurological disease. Furthermore, our multi-faceted approach for examining the functional significance of non-coding variants can be readily generalized to study other loci important for myelin structure and function.

Numerous genetic association studies have implicated the KIAA0319 gene on human chromosome 6p22 in dyslexia susceptibility. The causative variant(s) remains unknown but may modulate gene expression, given that (1) a dyslexia-associated haplotype has been implicated in the reduced expression of KIAA0319, and (2) the strongest association has been found for the region spanning exon 1 of KIAA0319. Here, we test the hypothesis that variant(s) responsible for reduced KIAA0319 expression resides on the risk haplotype close to the gene's transcription start site. We identified seven single-nucleotide polymorphisms on the risk haplotype immediately upstream of KIAA0319 and determined that three of these are strongly associated with multiple reading-related traits. Using luciferase-expressing constructs containing the KIAA0319 upstream region, we characterized the minimal promoter and additional putative transcriptional regulator regions. This revealed that the minor allele of rs9461045, which shows the strongest association with dyslexia in our sample (max p-value = 0.0001), confers reduced luciferase expression in both neuronal and non-neuronal cell lines. Additionally, we found that the presence of this rs9461045 dyslexia-associated allele creates a nuclear protein-binding site, likely for the transcriptional silencer OCT-1. Knocking down OCT-1 expression in the neuronal cell line SHSY5Y using an siRNA restores KIAA0319 expression from the risk haplotype to nearly that seen from the non-risk haplotype. Our study thus pinpoints a common variant as altering the function of a dyslexia candidate gene and provides an illustrative example of the strategic approach needed to dissect the molecular basis of complex genetic traits., Author Summary Dyslexia, or reading disability, is a common disorder caused by both genetic and environmental factors. Genetic studies have implicated a number of genes as candidates for playing a role in dyslexia. We functionally characterized one such gene (KIAA0319) to identify variant(s) that might affect gene expression and contribute to the disorder. We discovered a variant residing outside of the protein-coding region of KIAA0319 that reduces expression of the gene. This variant creates a binding site for the transcription factor OCT-1. Previous studies have shown that OCT-1 binding to a specific DNA sequence upstream of a gene can reduce the expression of that gene. In this case, reduced KIAA0319 expression could lead to improper development of regions of the brain involved in reading ability. This is the first study to identify a functional variant implicated in dyslexia. More broadly, our study illustrates the steps that can be utilized for identifying mutations causing other complex genetic disorders.

Linkage disequilibrium structure of KIAA0319 and DCDC2 , two candidate susceptibility genes for developmental dyslexia

Bioelectronic interfaces that facilitate electron transfer between the electrode and a dehydrogenase enzyme have potential applications in biosensors, biocatalytic reactors, and biological fuel cells. The secondary alcohol dehydrogenase (2 degrees ADH) from Thermoanaerobacter ethanolicus is especially well suited for the development of such bioelectronic interfaces because of its thermostability and facile production and purification. However, the natural cofactor for the enzyme, beta-nicotinamide adenine dinucleotide phosphate (NADP+), is more expensive and less stable than beta-nicotinamide adenine dinucleotide (NAD+). PCR-based, site-directed mutagenesis was performed on 2 degrees ADH in an attempt to adjust the cofactor specificity toward NAD+ by mutating Tyr218 to Phe (Y218F 2 degrees ADH). This mutation increased the Km(app) for NADP+ 200-fold while decreasing the Km(app) for NAD+ 2.5-fold. The mutant enzyme was incorporated into a bioelectronic interface that established electrical communication between the enzyme, the NAD+, the electron mediator toluidine blue O (TBO), and a gold electrode. Cyclic voltammetry, impedance spectroscopy, gas chromatography, mass spectrometry, constant potential amperometry, and chronoamperometry were used to characterize the mutant and wild-type enzyme incorporated in the bioelectronic interface. The Y218F 2 degrees ADH exhibited a fourfold increase in the turnover ratio compared to the wild type in the presence of NAD+. The electrochemical and kinetic measurements support the prediction that the Rossmann fold of the enzyme binds to the phosphate moiety of the cofactor. During the 45 min of continuous operation, NAD+ was electrically recycled 6.7 x 10(4) times, suggesting that the Y218F 2 degrees ADH-modified bioelectronic interface is stable.