Human lineage-specific amplification, selection, and neuronal expression of DUF1220 domains

Understanding the Emergence of Schizophrenia in the Light of Human Evolution: New Perspectives in Genetics

Schizophrenia is a frequent and disabling disease. The persistence of the disorder despite its harmful consequences represents an evolutionary paradox. Based on recent discoveries in genetics, scientists have formulated the "price-to-pay" hypothesis: schizophrenia would be intimately related to human evolution, particularly to brain development and human-specific higher cognitive functions. The objective of the present work is to question scientific literature about the relationship between schizophrenia and human evolution from a genetic point of view. In the last two decades, research investigated the association between schizophrenia and a few genetic evolutionary markers: Human accelerated regions, segmental duplications, and highly repetitive DNA such as the Olduvai domain. Other studies focused on the action of natural selection on schizophrenia-associated genetic variants, also thanks to the complete sequencing of archaic hominins' genomes (Neanderthal, Denisova). Results suggested that a connection between human evolution and schizophrenia may exist; nonetheless, much research is still needed, and it is possible that a definitive answer to the evolutionary paradox of schizophrenia will never be found.

Inferring DNA methylation in non-skeletal tissues of ancient specimens

Genome-wide premortem DNA methylation patterns can be computationally reconstructed from high-coverage DNA sequences of ancient samples. Because DNA methylation is more conserved across species than across tissues, and ancient DNA is typically extracted from bones and teeth, previous works utilizing ancient DNA methylation maps focused on studying evolutionary changes in the skeletal system. Here we suggest that DNA methylation patterns in one tissue may, under certain conditions, be informative on DNA methylation patterns in other tissues of the same individual. Using the fact that tissue-specific DNA methylation builds up during embryonic development, we identified the conditions that allow for such cross-tissue inference and devised an algorithm that carries it out. We trained the algorithm on methylation data from extant species and reached high precisions of up to 0.92 for validation datasets. We then used the algorithm on archaic humans, and identified more than 1,850 positions for which we were able to observe differential DNA methylation in prefrontal cortex neurons. These positions are linked to hundreds of genes, many of which are involved in neural functions such as structural and developmental processes. Six positions are located in the neuroblastoma breaking point family (NBPF) gene family, which probably played a role in human brain evolution. The algorithm we present here allows for the examination of epigenetic changes in tissues and cell types that are absent from the palaeontological record, and therefore provides new ways to study the evolutionary impacts of epigenetic changes.

Neuroblastoma Breakpoint Family 3mer Higher Order Repeats/Olduvai Triplet Pattern in the Complete Genome of Human and Nonhuman Primates and Relation to Cognitive Capacity

BACKGROUND/OBJECTIVES

The ~1.6 kb NBPF repeat units in neuroblastoma breakpoint family (NBPF) genes are specific to humans and are associated with cognitive capacity in higher primates. While the number of NBPF monomers/Olduvai sequences in humans is approximately 2-3 times greater than in great apes, the difference in copy number values of canonical NBPF 3mer Higher-order repeats (HORs)/Olduvai triplets between humans and great apes is substantially larger. This study aims to analyze the organization and evolutionary significance of NBPF 3mer HORs/Olduvai triplets in fully sequenced primate genomes.

METHODS

We applied the global repeat map (GRM) algorithm to identify canonical and variant NBPF 3mer HORs/Olduvai triplets in the complete genomes of humans, chimpanzees, gorillas, and orangutans. The resulting monomer arrays were analyzed using the GRMhor algorithm to generate detailed schematic representations of NBPF HOR organization.

RESULTS

The analysis reveals a distinct difference in NBPF-related patterns among these primates, particularly in the number of tandemly organized canonical 3mer HORs/Olduvai triplets: 61 tandemly organized canonical NBPF 3mer HORs/Olduvai triplets in humans, compared to 0 in chimpanzees and orangutans, and 9 in gorillas. When considering only tandemly organized 3mer HORs/Olduvai triplets with more than three copies, the numbers adjust to 36 in humans and 0 in great apes. Furthermore, the divergence between individual NBPF monomers in humans and great apes is twice as high as that observed within great apes.

CONCLUSIONS

These findings support the hypothesis that the tandem organization of NBPF 3mer HORs/Olduvai triplets plays a crucial role in enhancing cognitive capacity in humans compared to great apes, potentially providing a significant evolutionary advantage. This effect complements the impact of the increased number of individual NBPF monomers/Olduvai sequences, together contributing to a synergistic amplification effect.

Olduvai domain expression downregulates mitochondrial pathways: implications for human brain evolution and neoteny

Encoded by the NBPF gene family, Olduvai (formerly DUF1220) protein domains have undergone the largest human lineage-specific copy number expansion of any coding region in the genome. Olduvai copy number shows a linear relationship with several brain size-related measures and cortical neuron number among primates and with normal and disease-associated (micro- and macrocephaly) variation in brain size in human populations. While Olduvai domains have been shown to promote proliferation of neural stem cells, the mechanism underlying such effects has remained unclear. Here, we investigate the function of Olduvai by transcriptome and proteome analyses of cells overexpressing NBPF1, a gene encoding 7 Olduvai domains. Our results from both RNAseq and mass spectrometry approaches suggest a potential downregulation of mitochondria. In our proteomics study, a Gene Ontology (GO) enrichment analysis for the downregulated proteins revealed a striking overrepresentation of the biological process related to the mitochondrial electron transport chain (p value: 1.81e-11) and identified deregulation of the NADH dehydrogenase activity (p value: 2.43e-11) as the primary molecular function. We verify the reduction of apparent mitochondria via live-cell imaging experiments. Given these and previous Olduvai findings, we suggest that the Olduvai-mediated, dosage-dependent reduction in available energy via mitochondrial downregulation may have resulted in a developmental slowdown such that the neurogenic window among primates, and most extremely in humans, was expanded over a greater time interval, allowing for production of greater numbers of neurons and a larger brain. We further suggest that such a slowdown may extend to other developmental processes that also exhibit neotenic features.

Unraveling mechanisms of human brain evolution

Evolutionary changes in human brain structure and function have enabled our specialized cognitive abilities. How these changes have come about genetically and functionally has remained an open question. However, new methods are providing a wealth of information about the genetic, epigenetic, and transcriptomic differences that set the human brain apart. Combined with in vitro models that allow access to developing brain tissue and the cells of our closest living relatives, the puzzle pieces are now coming together to yield a much more complete picture of what is actually unique about the human brain. The challenge now will be linking these observations and making the jump from correlation to causation. However, elegant genetic manipulations are now possible and, when combined with model systems such as organoids, will uncover a mechanistic understanding of how evolutionary changes at the genetic level have led to key differences in development and function that enable human cognition.

Independent expansion, selection and hypervariability of the TBC1D3 gene family in humans

TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on chromosome 17. We find that most humans vary along two TBC1D3 clusters where human haplotypes are highly variable in copy number, differing by as many as 20 copies, and structure (structural heterozygosity 90%). We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Lastly, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL. These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.

A potential patient stratification biomarker for Parkinson´s disease based on LRRK2 kinase-mediated centrosomal alterations in peripheral blood-derived cells

Parkinson´s disease (PD) is a common neurodegenerative movement disorder and leucine-rich repeat kinase 2 (LRRK2) is a promising therapeutic target for disease intervention. However, the ability to stratify patients who will benefit from such treatment modalities based on shared etiology is critical for the success of disease-modifying therapies. Ciliary and centrosomal alterations are commonly associated with pathogenic LRRK2 kinase activity and can be detected in many cell types. We previously found centrosomal deficits in immortalized lymphocytes from G2019S-LRRK2 PD patients. Here, to investigate whether such deficits may serve as a potential blood biomarker for PD which is susceptible to LRKK2 inhibitor treatment, we characterized patient-derived cells from distinct PD cohorts. We report centrosomal alterations in peripheral cells from a subset of early-stage idiopathic PD patients which is mitigated by LRRK2 kinase inhibition, supporting a role for aberrant LRRK2 activity in idiopathic PD. Centrosomal defects are detected in R1441G-LRRK2 and G2019S-LRRK2 PD patients and in non-manifesting LRRK2 mutation carriers, indicating that they accumulate prior to a clinical PD diagnosis. They are present in immortalized cells as well as in primary lymphocytes from peripheral blood. These findings indicate that analysis of centrosomal defects as a blood-based patient stratification biomarker may help nominate idiopathic PD patients who will benefit from LRRK2-related therapeutics.

EASTR: Identifying and eliminating systematic alignment errors in multi-exon genes

Accurate alignment of transcribed RNA to reference genomes is a critical step in the analysis of gene expression, which in turn has broad applications in biomedical research and in the basic sciences. We reveal that widely used splice-aware aligners, such as STAR and HISAT2, can introduce erroneous spliced alignments between repeated sequences, leading to the inclusion of falsely spliced transcripts in RNA-seq experiments. In some cases, the 'phantom' introns resulting from these errors make their way into widely-used genome annotation databases. To address this issue, we present EASTR (Emending Alignments of Spliced Transcript Reads), a software tool that detects and removes falsely spliced alignments or transcripts from alignment and annotation files. EASTR improves the accuracy of spliced alignments across diverse species, including human, maize, and Arabidopsis thaliana, by detecting sequence similarity between intron-flanking regions. We demonstrate that applying EASTR before transcript assembly substantially reduces false positive introns, exons, and transcripts, improving the overall accuracy of assembled transcripts. Additionally, we show that EASTR's application to reference annotation databases can detect and correct likely cases of mis-annotated transcripts.

Tandem NBPF 3mer HORs (Olduvai triplets) in Neanderthal and two novel HOR tandem arrays in human chromosome 1 T2T-CHM13 assembly

It is known that the ~ 1.6 kb Neuroblastoma BreakPoint Family (NBPF) repeats are human specific and contributing to cognitive capabilities, with increasing frequency in higher order repeat 3mer HORs (Olduvai triplets). From chimpanzee to modern human there is a discontinuous jump from 0 to ~ 50 tandemly organized 3mer HORs. Here we investigate the structure of NBPF 3mer HORs in the Neanderthal genome assembly of Pääbo et al., comparing it to the results obtained for human hg38.p14 chromosome 1. Our findings reveal corresponding NBPF 3mer HOR arrays in Neanderthals with slightly different monomer structures and numbers of HOR copies compared to humans. Additionally, we compute the NBPF 3mer HOR pattern for the complete telomere-to-telomere human genome assembly (T2T-CHM13) by Miga et al., identifying two novel tandem arrays of NBPF 3mer HOR repeats with 5 and 9 NBPF 3mer HOR copies. We hypothesize that these arrays correspond to novel NBPF genes (here referred to as NBPFA1 and NBPFA2). Further improving the quality of the Neanderthal genome using T2T-CHM13 as a reference would be of great interest in determining the presence of such distant novel NBPF genes in the Neanderthal genome and enhancing our understanding of human evolution.

Comparative physiological anthropogeny: exploring molecular underpinnings of distinctly human phenotypes

Anthropogeny is a classic term encompassing transdisciplinary investigations of the origins of the human species. Comparative anthropogeny is a systematic comparison of humans and other living nonhuman hominids (so-called "great apes"), aiming to identify distinctly human features in health and disease, with the overall goal of explaining human origins. We begin with a historical perspective, briefly describing how the field progressed from the earliest evolutionary insights to the current emphasis on in-depth molecular and genomic investigations of "human-specific" biology and an increased appreciation for cultural impacts on human biology. While many such genetic differences between humans and other hominids have been revealed over the last two decades, this information remains insufficient to explain the most distinctive phenotypic traits distinguishing humans from other living hominids. Here we undertake a complementary approach of "comparative physiological anthropogeny," along the lines of the preclinical medical curriculum, i.e., beginning with anatomy and considering each physiological system and in each case considering genetic and molecular components that are relevant. What is ultimately needed is a systematic comparative approach at all levels from molecular to physiological to sociocultural, building networks of related information, drawing inferences, and generating testable hypotheses. The concluding section will touch on distinctive considerations in the study of human evolution, including the importance of gene-culture interactions.

Proteolytic activation of human-specific Olduvai domains by the furin protease

Olduvai protein domains (formerly DUF1220) show the greatest human-specific increase in copy number of any coding region in the genome and are highly correlated with human brain evolution and cognitive disease. The majority of human copies are found within four NBPF genes organized in a variable number of a tandemly arranged three-domain blocks called Olduvai triplets. Here we show that these human-specific Olduvai domains are posttranslationally processed by the furin protease, with a cleavage site occurring once at each triplet. These findings suggest that all expanded human-specific NBPF genes encode proproteins consisting of many independent Olduvai triplet proteins which are activated by furin processing. The exceptional correlation of Olduvai copy number and brain size taken together with our new furin data, indicates the ultimate target of selection was a rapid increase in dosage of autonomously functioning Olduvai triplet proteins, and that these proteins are the primary active agent underlying Olduvai's role in human brain expansion.

Keeping the balance: Trade-offs between human brain evolution, autism, and schizophrenia

The unique qualities of the human brain are a product of a complex evolutionary process. Evolution, famously described by François Jacob as a “tinkerer,” builds upon existing genetic elements by modifying and repurposing them for new functions. Genetic changes in DNA may lead to the emergence of new genes or cause altered gene expression patterns. Both gene and regulatory element mutations may lead to new functions. Yet, this process may lead to side-effects. An evolutionary trade-off occurs when an otherwise beneficial change, which is important for evolutionary success and is under strong positive selection, concurrently results in a detrimental change in another trait. Pleiotropy occurs when a gene affects multiple traits. Antagonistic pleiotropy is a phenomenon whereby a genetic variant leads to an increase in fitness at one life-stage or in a specific environment, but simultaneously decreases fitness in another respect. Therefore, it is conceivable that the molecular underpinnings of evolution of highly complex traits, including brain size or cognitive ability, under certain conditions could result in deleterious effects, which would increase the susceptibility to psychiatric or neurodevelopmental diseases. Here, we discuss possible trade-offs and antagonistic pleiotropies between evolutionary change in a gene sequence, dosage or activity and the susceptibility of individuals to autism spectrum disorders and schizophrenia. We present current knowledge about genes and alterations in gene regulatory landscapes, which have likely played a role in establishing human-specific traits and have been implicated in those diseases.

Tandemly repeated NBPF HOR copies (Olduvai triplets): Possible impact on human brain evolution

Previously it was found that the neuroblastoma breakpoint family (NBPF) gene repeat units of ∼1.6 kb have an important role in human brain evolution and function. The higher order organization of these repeat units has been discovered by both methods, the higher order repeat (HOR)-searching method and the HLS searching method. Using the HOR searching method with global repeat map algorithm, here we identified the tandemly organized NBPF HORs in the human and nonhuman primate NCBI reference genomes. We identified 50 tandemly organized canonical 3mer NBPF HOR copies (Olduvai triplets), but none in nonhuman primates chimpanzee, gorilla, orangutan, and Rhesus macaque. This discontinuous jump in tandemly organized HOR copy number is in sharp contrast to the known gradual increase in the number of Olduvai domains (NBPF monomers) from nonhuman primates to human, especially from ∼138 in chimpanzee to ∼300 in human genome. Using the same global repeat map algorithm method we have also determined the 3mer tandems of canonical 3mer HOR copies in 20 randomly chosen human genomes (10 male and 10 female). In all cases, we found the same 3mer HOR copy numbers as in the case of the reference human genome, with no mutation. On the other hand, some point mutations with respect to reference genome are found for some NBPF monomers which are not tandemly organized in canonical HORs.

Selection of Olduvai Domains during Evolution: A Role for Primate-Specific Splicing Super-Enhancer and RNA Guanine Quadruplex in Bipartite NBPF Exons

Olduvai protein domains (also known as DUF1220 or NBPF) have undergone the greatest human-specific increase in the copy number of any coding region in the genome. Their repeat number was strongly associated with the evolutionary expansion of brain volumes, neuron counts and cognitive abilities, as well as with disorders of the autistic spectrum. Nevertheless, the domain function and cellular mechanisms underlying the positive selection of Olduvai DNA sequences in higher primates remain obscure. Here, I show that the inclusion of Olduvai exon doublets in mature transcripts is facilitated by a potent splicing enhancer that was created through duplication within the first exon. The enhancer is the strongest among the NBPF transcripts and further promotes the already high splicing activity of the unexpanded first exons of the two-exon domains, safeguarding the expanded Olduvai exon doublets in the mature transcriptome. The duplication also creates a predicted RNA guanine quadruplex that may regulate the access to spliceosomal components of the super-enhancer and influence the splicing of adjacent exons. Thus, positive Olduvai selection during primate evolution is likely to result from a combination of multiple targets in gene expression pathways, including RNA splicing.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Evolution of genetic mechanisms regulating cortical neurogenesis

The size of the cerebral cortex increases dramatically across amniotes, from reptiles to great apes. This is primarily due to different numbers of neurons and glial cells produced during embryonic development. The evolutionary expansion of cortical neurogenesis was linked to changes in neural stem and progenitor cells, which acquired increased capacity of self‐amplification and neuron production. Evolution works via changes in the genome, and recent studies have identified a small number of new genes that emerged in the recent human and primate lineages, promoting cortical progenitor proliferation and increased neurogenesis. However, most of the mammalian genome corresponds to noncoding DNA that contains gene‐regulatory elements, and recent evidence precisely points at changes in expression levels of conserved genes as key in the evolution of cortical neurogenesis. Here, we provide an overview of basic cellular mechanisms involved in cortical neurogenesis across amniotes, and discuss recent progress on genetic mechanisms that may have changed during evolution, including gene expression regulation, leading to the expansion of the cerebral cortex.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

What Are the Human-Specific Aspects of Neocortex Development?

When considering what makes us human, the development of the neocortex, the seat of our higher cognitive abilities, is of central importance. Throughout this complex developmental process, neocortical stem and progenitor cells (NSPCs) exert a priming role in determining neocortical tissue fate, through a series of cellular and molecular events. In this Perspective article, we address five questions of relevance for potentially human-specific aspects of NSPCs, (i) Are there human-specific NSPC subtypes? (ii) What is the functional significance of the known temporal differences in NSPC dynamics between human and other great apes? (iii) Are there functional interactions between the human-specific genes preferentially expressed in NSPCs? (iv) Do humans amplify certain metabolic pathways for NSPC proliferation? and finally (v) Have differences evolved during human evolution, notably between modern humans and Neandertals, that affect the performance of key genes operating in NSPCs? We discuss potential implications inherent to these questions, and suggest experimental approaches on how to answer them, hoping to provide incentives to further understand key issues of human cortical development.

Solution NMR backbone assignments of disordered Olduvai protein domain CON1 employing Hα-detected experiments

Olduvai protein domains, encoded by the NBPF gene family, are responsible for the largest increase in copy number of any protein-coding region in the human genome. This has spawned various genetics studies which have linked these domains to human brain development and divergence from our primate ancestors, as well as currently relevant cognitive diseases such as schizophrenia and autism spectrum disorder (ASD). There are six separate Olduvai domains which together form the majority of the various protein products of the NBPF genes. The six domains include three conserved domains (CON1-3), and three human-lineage-specific domains (HLS1-3) which occur in triplet. Here, we present the solution nuclear magnetic resonance backbone assignments for the CON1 domain, which has been linked to the severity of ASD. The data confirm that CON1 is an intrinsically disordered protein (IDP). Additionally, we use innovative Hα-detected experiments which allow us to not only assign the Hα atoms and N atoms of proline residues, but also to assign residues where HN-experiments suffered from peak overlap or broadening.

Case Report: Neuroblastoma Breakpoint Family Genes Associate With 1q21 Copy Number Variation Disorders

Microduplications and reciprocal microdeletions of chromosome 1q21. 1 and/or 1q21.2 have been linked to variable clinical features, but the underlying pathogenic gene(s) remain unclear. Here we report that distinct microduplications were detected on chromosome 1q21.2 (GRCh37/hg19) in a mother (255 kb in size) and her newborn daughter (443 kb in size), while the same paternal locus was wild-type. Although the two microduplications largely overlap in genomic sequence (183 kb overlapping), the mother showed no clinical phenotype while the daughter presented with several features that are commonly observed on 1q21 microduplication or microdeletion patients, including developmental delay, craniofacial dysmorphism, congenital heart disease and sensorineural hearing loss. NBPF15 and NBPF16, two involved genes that are exclusively duplicated in the proband, may be the cause of the clinical manifestations. This study supports an association between NBPF genes and 1q21 copy number variation disorders.

dSPRINT: predicting DNA, RNA, ion, peptide and small molecule interaction sites within protein domains

Domains are instrumental in facilitating protein interactions with DNA, RNA, small molecules, ions and peptides. Identifying ligand-binding domains within sequences is a critical step in protein function annotation, and the ligand-binding properties of proteins are frequently analyzed based upon whether they contain one of these domains. To date, however, knowledge of whether and how protein domains interact with ligands has been limited to domains that have been observed in co-crystal structures; this leaves approximately two-thirds of human protein domain families uncharacterized with respect to whether and how they bind DNA, RNA, small molecules, ions and peptides. To fill this gap, we introduce dSPRINT, a novel ensemble machine learning method for predicting whether a domain binds DNA, RNA, small molecules, ions or peptides, along with the positions within it that participate in these types of interactions. In stringent cross-validation testing, we demonstrate that dSPRINT has an excellent performance in uncovering ligand-binding positions and domains. We also apply dSPRINT to newly characterize the molecular functions of domains of unknown function. dSPRINT's predictions can be transferred from domains to sequences, enabling predictions about the ligand-binding properties of 95% of human genes. The dSPRINT framework and its predictions for 6503 human protein domains are freely available at http://protdomain.princeton.edu/dsprint.

Human-Specific Genes, Cortical Progenitor Cells, and Microcephaly

Over the past few years, human-specific genes have received increasing attention as potential major contributors responsible for the 3-fold difference in brain size between human and chimpanzee. Accordingly, mutations affecting these genes may lead to a reduction in human brain size and therefore, may cause or contribute to microcephaly. In this review, we will concentrate, within the brain, on the cerebral cortex, the seat of our higher cognitive abilities, and focus on the human-specific gene ARHGAP11B and on the gene family comprising the three human-specific genes NOTCH2NLA, -B, and -C. These genes are thought to have significantly contributed to the expansion of the cerebral cortex during human evolution. We will summarize the evolution of these genes, as well as their expression and functional role during human cortical development, and discuss their potential relevance for microcephaly. Furthermore, we will give an overview of other human-specific genes that are expressed during fetal human cortical development. We will discuss the potential involvement of these genes in microcephaly and how these genes could be studied functionally to identify a possible role in microcephaly.

The Polygenic Nature and Complex Genetic Architecture of Specific Learning Disorder

Specific Learning Disorder (SLD) is a multifactorial, neurodevelopmental disorder which may involve persistent difficulties in reading (dyslexia), written expression and/or mathematics. Dyslexia is characterized by difficulties with speed and accuracy of word reading, deficient decoding abilities, and poor spelling. Several studies from different, but complementary, scientific disciplines have investigated possible causal/risk factors for SLD. Biological, neurological, hereditary, cognitive, linguistic-phonological, developmental and environmental factors have been incriminated. Despite worldwide agreement that SLD is highly heritable, its exact biological basis remains elusive. We herein present: (a) an update of studies that have shaped our current knowledge on the disorder’s genetic architecture; (b) a discussion on whether this genetic architecture is ‘unique’ to SLD or, alternatively, whether there is an underlying common genetic background with other neurodevelopmental disorders; and, (c) a brief discussion on whether we are at a position of generating meaningful correlations between genetic findings and anatomical data from neuroimaging studies or specific molecular/cellular pathways. We conclude with open research questions that could drive future research directions.

Evaluation of domain of unknown function 1220 (DUF1220) for detection of human genome by quantitative polymerase chain reaction: Potential use in assessing the biodistribution of transplanted therapeutic human cells

The biodistribution profile of cell-based therapy products in animal models is important for evaluation of their safety and efficacy. Because of its quantitative nature and sensitivity, the quantitative polymerase chain reaction (qPCR) is a useful method for detecting and quantifying xenogeneic cell-derived DNA in animal models, thereby allowing a biodistribution profile to be established. Although the restriction endonuclease family from Arthrobacter luteus (Alu) of repetitive elements in human genome sequences has been used to assess the biodistribution of human cells, high background signals are detected. In the present study, we evaluate the potential of domain of unknown function 1220 (DUF1220), which is a human lineage-specific, multiple-copy gene similar to Alu sequences, for such analysis. Using qPCR analysis for DUF1220, human genome could be detected against a mouse genome background at a level comparable to that of Alu sequences with no detectable background signals. Moreover, using this approach, the human genome could be distinguished from the cynomolgus monkey genome. Further investigation of the quantitative aspects of this DUF1220-based qPCR assay might prove its usefulness for biodistribution studies of human cells transplanted into animals in the future.

An integrative understanding of comparative cognition: lessons from human brain evolution

A comprehensive understanding of animal cognition requires the integration of studies on behavior, electrophysiology, neuroanatomy, development, and genomics. Although studies of comparative cognition are receiving increasing attention from organismal biologists, most current studies focus on the comparison of behaviors and anatomical structures to understand their adaptative values. However, to understand the most potentially complex cognitive program of the human brain a greater synthesis of a multitude of disciplines is needed. In this review, we start with extensive neuroanatomic comparisons between humans and other primates. One likely specialization of the human brain is the expansion of neocortex, especially in regions for high-order cognition (e.g., prefrontal cortex). We then discuss how such an expansion can be linked to heterochrony of the brain developmental program, resulting in a greater number of neurons and enhanced computational capacity. Furthermore, alteration of gene expression in the human brain has been associated with positive selection in DNA sequences of gene regulatory regions. These results not only imply that genes associated with brain development are a major factor in the evolution of cognition, but also that high-quality whole-genome sequencing and gene manipulation techniques are needed for an integrative and functional understanding of comparative cognition in non-model organisms.

Human core duplicon gene families: game changers or game players?

Illuminating the role of specific gene duplications within the human lineage can provide insights into human-specific adaptations. The so-called human core duplicon gene families have received particular attention in this respect, due to special features, such as expansion along single chromosomes, newly acquired protein domains and signatures of positive selection. Here, we summarize the data available for 10 such families and include some new analyses. A picture emerges that suggests broad functions for these protein families, possibly through modification of core cellular pathways. Still, more dedicated studies are required to elucidate the function of core-duplicons gene families and how they have shaped adaptations and evolution of humans.

The N-terminal of NBPF15 causes multiple types of aggregates and mediates phase transition

The neuroblastoma breakpoint family (NBPF) consists of 24 members that play an important role in neuroblastoma and other cancers. NBPF is an evolutionarily recent gene family that encodes several repeats of Olduvai domain and an abundant N-terminal region. The function and biochemical properties of both Olduvai domain and the N-terminal region remain enigmatic. Human NBPF15 encodes a 670 AA protein consisting of six clades of Olduvai domains. In this study, we synthesized and expressed full-length NBPF15, and purified a range of NBPF15 truncations which were analyzed using dynamic light scattering (DLS), superdex200 (S200), small-angle X-ray scattering (SAXS), far-UV circular dichroism (CD) spectroscopy, transmission electron microscope (TEM), and crystallography. We found that proteins containing both the N-terminal region and Olduvai domain are heterogeneous with multiple types of aggregates, and some of them underwent a liquid-to-solid phase transition, probably because of the entanglement within the N-terminal coiled-coil. Proteins that contain only the Olduvai domain are homogeneous extended monomers, and those with the conserved clade 1 (CON1) have manifested a tendency to crystallize. We suggest that the entanglements between the mosaic disorder-ordered segments in NBPF15 N terminus have triggered the multiple types of aggregates and phase transition of NBPF15 proteins, which could be associated with Olduvai-related cognitive dysfunction diseases.

Decoding the development of the human hippocampus

The hippocampus is an important part of the limbic system in the human brain that has essential roles in spatial navigation and the consolidation of information from short-term memory to long-term memory^1,2. Here we use single-cell RNA sequencing and assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis to illustrate the cell types, cell linage, molecular features and transcriptional regulation of the developing human hippocampus. Using the transcriptomes of 30,416 cells from the human hippocampus at gestational weeks 16-27, we identify 47 cell subtypes and their developmental trajectories. We also identify the migrating paths and cell lineages of PAX6⁺ and HOPX⁺ hippocampal progenitors, and regional markers of CA1, CA3 and dentate gyrus neurons. Multiomic data have uncovered transcriptional regulatory networks of the dentate gyrus marker PROX1. We also illustrate spatially specific gene expression in the developing human prefrontal cortex and hippocampus. The molecular features of the human hippocampus at gestational weeks 16-20 are similar to those of the mouse at postnatal days 0-5 and reveal gene expression differences between the two species. Transient expression of the primate-specific gene NBPF1 leads to a marked increase in PROX1⁺ cells in the mouse hippocampus. These data provides a blueprint for understanding human hippocampal development and a tool for investigating related diseases.

The Driver of Extreme Human-Specific Olduvai Repeat Expansion Remains Highly Active in the Human Genome

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

The somatic mutation landscape of Chinese Colorectal Cancer

Colorectal cancer (CRC) is the fifth leading cause of cancer-related death in China. The incidence of Chinese CRC has increased dramatically with the changes of dietary and lifestyle. However, the genetic landscape of Chinese colorectal cancer mutation is still poorly understood. In this study, we have performed whole exome-sequencing analysis of 63 CRC cases. We found that Chinese CRC were hypermutated, which were enriched in ECM-receptor interaction, antigen processing and presentation, and focal adhesion. Analysis with clinical characteristics indicated that the deficiency of CRC driver gene, FCGBP and NBPF1 conferred CRC development and was showed worse survival rates, which could be the novel regulators and, diagnostic and prognostic biomarkers for Chinese CRC. Taken together, the application of whole exome-sequencing unveiled previously unsuspected somatic mutation landscape in Chinese CRCs, which may expand the understanding of disease mechanisms and provide an alternative personalized treatment for Chinese CRC patients.

Solution NMR backbone assignment reveals interaction-free tumbling of human lineage-specific Olduvai protein domains

Olduvai protein domains, encoded primarily by NBPF genes, have been linked to both human brain evolution and cognitive diseases such as autism and schizophrenia. There are six primary domains that comprise the Olduvai family: three conserved domains (CON1-3) and three human lineage-specific domains (HLS1-3), which typically occur as a triplet (HLS1, HLS2 and HLS3). Herein, we present the solution NMR assignment of the backbone chemical shifts of the separate HLS1, 2 and 3 domains of NBPF15. Our data suggest that there is no change in the structure of the separate domains when compared to the full-length triplet (HLS1-HLS2-HLS3). We also demonstrate that there is no direct interaction between the three domains.

A Third Linear Association Between Olduvai (DUF1220) Copy Number and Severity of the Classic Symptoms of Inherited Autism

OBJECTIVE

The authors previously reported that the copy number of sequences encoding an Olduvai protein domain subtype (CON1) shows a linear association with the severity of social deficits and communication impairment in individuals with autism. In this study, using an improved measurement method, the authors replicated this association in an independent population.

METHOD

The authors obtained whole genome sequence (WGS) data and phenotype data on 215 individuals from the Autism Speaks MSSNG project. They derived copy number from WGS data using a modified sequence read-depth technique. A linear mixed-effects model was used to test the association between Olduvai CON1 copy number and symptom severity as measured by the Autism Diagnostic Interview-Revised. The authors then combined data from previous studies (N=524) for final analyses.

RESULTS

A significant linear association was observed between CON1 copy number and social diagnostic score (SDS) (β=0.24) and communicative diagnostic score (CDS) (β=0.23). Using the combined data, the authors present strong significant associations of CON1 dosage with SDS (β=0.18) and CDS (β=0.13). The authors also implicate Olduvai subtypes found in two genes, NBPF1 and NBPF14 (R²=6.2%). Associations were preferentially found in multiplex versus simplex families.

CONCLUSIONS

The finding of a third dose-dependent association between Olduvai sequences and autism severity, preferentially in multiplex families, provides strong evidence that this highly duplicated and underexamined protein domain family plays an important role in inherited autism.

Evolutionary Changes in Transcriptional Regulation: Insights into Human Behavior and Neurological Conditions

Understanding the biological basis for human-specific cognitive traits presents both immense challenges and unique opportunities. Although the question of what makes us human has been investigated with several different methods, the rise of comparative genomics, epigenomics, and medical genetics has provided tools to help narrow down and functionally assess the regions of the genome that seem evolutionarily relevant along the human lineage. In this review, we focus on how medical genetic cases have provided compelling functional evidence for genes and loci that appear to have interesting evolutionary signatures in humans. Furthermore, we examine a special class of noncoding regions, human accelerated regions (HARs), that have been suggested to show human-lineage-specific divergence, and how the use of clinical and population data has started to provide functional information to examine these regions. Finally, we outline methods that provide new insights into functional noncoding sequences in evolution.

Paired involvement of human-specific Olduvai domains and NOTCH2NL genes in human brain evolution

Sequences encoding Olduvai (DUF1220) protein domains show the largest human-specific increase in copy number of any coding region in the genome and have been linked to human brain evolution. Most human-specific copies of Olduvai (119/165) are encoded by three NBPF genes that are adjacent to three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Here, employing genomic, phylogenetic, and transcriptomic evidence, we show that these NOTCH2NL/NBPF gene pairs evolved jointly, as two-gene units, very recently in human evolution, and are likely co-regulated. Remarkably, while three NOTCH2NL paralogs were added, adjacent Olduvai sequences hyper-amplified, adding 119 human-specific copies. The data suggest that human-specific Olduvai domains and adjacent NOTCH2NL genes may function in a coordinated, complementary fashion to promote neurogenesis and human brain expansion in a dosage-related manner.

Genetics of human brain evolution

During the course of evolution the human brain has increased in size and complexity, ultimately these differences are the result of changes at the genetic level. Identifying and characterizing molecular evolution requires an understanding of both the genetic underpinning of the system as well as the comparative genetic tools to identify signatures of selection. This chapter aims to describe our current understanding of the genetics of human brain evolution. Primarily this is the story of the evolution of the human brain since our last common ape ancestor, but where relevant we will also discuss changes that are unique to the primate brain (compared to other mammals) or various other lineages in the evolution of humans more generally. It will focus on genetic changes that both directly affected the development and function of the brain as well as those that have indirectly influenced brain evolution through both prenatal and postnatal environment. This review is not meant to be exhaustive, but rather to begin to construct a general framework for understanding the full array of data being generated.

Mathematically universal and biologically consistent astrocytoma genotype encodes for transformation and predicts survival phenotype

DNA alterations have been observed in astrocytoma for decades. A copy-number genotype predictive of a survival phenotype was only discovered by using the generalized singular value decomposition (GSVD) formulated as a comparative spectral decomposition. Here, we use the GSVD to compare whole-genome sequencing (WGS) profiles of patient-matched astrocytoma and normal DNA. First, the GSVD uncovers a genome-wide pattern of copy-number alterations, which is bounded by patterns recently uncovered by the GSVDs of microarray-profiled patient-matched glioblastoma (GBM) and, separately, lower-grade astrocytoma and normal genomes. Like the microarray patterns, the WGS pattern is correlated with an approximately one-year median survival time. By filling in gaps in the microarray patterns, the WGS pattern reveals that this biologically consistent genotype encodes for transformation via the Notch together with the Ras and Shh pathways. Second, like the GSVDs of the microarray profiles, the GSVD of the WGS profiles separates the tumor-exclusive pattern from normal copy-number variations and experimental inconsistencies. These include the WGS technology-specific effects of guanine-cytosine content variations across the genomes that are correlated with experimental batches. Third, by identifying the biologically consistent phenotype among the WGS-profiled tumors, the GBM pattern proves to be a technology-independent predictor of survival and response to chemotherapy and radiation, statistically better than the patient's age and tumor's grade, the best other indicators, and MGMT promoter methylation and IDH1 mutation. We conclude that by using the complex structure of the data, comparative spectral decompositions underlie a mathematically universal description of the genotype-phenotype relations in cancer that other methods miss.

Recurrent homozygous deletion of DROSHA and microduplication of PDE4DIP in pineoblastoma

Pineoblastoma is a rare and highly aggressive brain cancer of childhood, histologically belonging to the spectrum of primitive neuroectodermal tumors. Patients with germline mutations in DICER1, a ribonuclease involved in microRNA processing, have increased risk of pineoblastoma, but genetic drivers of sporadic pineoblastoma remain unknown. Here, we analyzed pediatric and adult pineoblastoma samples (n = 23) using a combination of genome-wide DNA methylation profiling and whole-exome sequencing or whole-genome sequencing. Pediatric and adult pineoblastomas showed distinct methylation profiles, the latter clustering with lower-grade pineal tumors and normal pineal gland. Recurrent variants were found in genes involved in PKA- and NF-κB signaling, as well as in chromatin remodeling genes. We identified recurrent homozygous deletions of DROSHA, acting upstream of DICER1 in microRNA processing, and a novel microduplication involving chromosomal region 1q21 containing PDE4DIP (myomegalin), comprising the ancient DUF1220 protein domain. Expresion of PDE4DIP and DUF1220 proteins was present exclusively in pineoblastoma with PDE4DIP gain.

Construction, expression, and characterization of AG11-843 and AG11-1581

AG1, a member of the DUF1220 protein family, exhibits the most extreme human lineage-specific copy number expansion of any protein-coding sequence in the genome. These variations in copy number have been linked to both brain evolution among primates and brain size in humans. Unfortunately, our current understanding of the structure and function of these proteins is limited because current cloning and expression techniques fail to consistently produce recombinant protein for in vitro studies. The present work describes a method for amino acid and DNA sequence optimization and synthesis, recombinant protein expression and analysis of two AG1 fragments, AG1^1-843 and AG1^1-1581. It was first necessary to modify the nucleotide sequence, while holding the GC content at 52.9%. The genes were then sectionally synthesized by overlap PCR. The resulting segments were cloned into the pET-15 b-sumo expression vector and subsequently transformed into BL21 (DE3) cells. After inducing their expression, the AG1^1-843 and AG1^1-1581 proteins were isolated and purified. Furthermore, using dynamic light scattering and gel filtration analysis, AG1^1-843 and AG1^1-1581 were shown to be present in tetrameric and dimeric forms in solution. To our knowledge, this is the first study to synthesize and express fragments of the DUF1220 protein family for in vitro analysis. Taken together, the proven utility and versatility of this method indicate that it can be used as an effective technique to construct and express other proteins with complicated sequences, thus providing the means to study their function and structure in vitro.

Human-Specific NOTCH2NL Genes Affect Notch Signaling and Cortical Neurogenesis

Genetic changes causing brain size expansion in human evolution have remained elusive. Notch signaling is essential for radial glia stem cell proliferation and is a determinant of neuronal number in the mammalian cortex. We find that three paralogs of human-specific NOTCH2NL are highly expressed in radial glia. Functional analysis reveals that different alleles of NOTCH2NL have varying potencies to enhance Notch signaling by interacting directly with NOTCH receptors. Consistent with a role in Notch signaling, NOTCH2NL ectopic expression delays differentiation of neuronal progenitors, while deletion accelerates differentiation into cortical neurons. Furthermore, NOTCH2NL genes provide the breakpoints in 1q21.1 distal deletion/duplication syndrome, where duplications are associated with macrocephaly and autism and deletions with microcephaly and schizophrenia. Thus, the emergence of human-specific NOTCH2NL genes may have contributed to the rapid evolution of the larger human neocortex, accompanied by loss of genomic stability at the 1q21.1 locus and resulting recurrent neurodevelopmental disorders.

Genomic trade-offs: are autism and schizophrenia the steep price of the human brain?

Evolution often deals in genomic trade-offs: changes in the genome that are beneficial overall persist even though they also produce disease in a subset of individuals. Here, we explore the possibility that such trade-offs have occurred as part of the evolution of the human brain. Specifically, we provide support for the possibility that the same key genes that have been major contributors to the rapid evolutionary expansion of the human brain and its exceptional cognitive capacity also, in different combinations, are significant contributors to autism and schizophrenia. Furthermore, the model proposes that one of the primary genes behind this trade-off may not technically be "a gene" or "genes" but rather are the highly duplicated sequences that encode the Olduvai protein domain family (formerly called DUF1220). This is not an entirely new idea. Others have proposed that the same genes involved in schizophrenia were also critical to the rapid expansion of the human brain, a view that has been expressed as "the same 'genes' that drive us mad have made us human". What is new is that a "gene", or more precisely a protein domain family, has been found that may satisfy these requirements.

Human Accelerated Regions and Other Human-Specific Sequence Variations in the Context of Evolution and Their Relevance for Brain Development

The review discusses, in a format of a timeline, the studies of different types of genetic variants, present in Homo sapiens, but absent in all other primate, mammalian, or vertebrate species, tested so far. The main characteristic of these variants is that they are found in regions of high evolutionary conservation. These sequence variations include single nucleotide substitutions (called human accelerated regions), deletions, and segmental duplications. The rationale for finding such variations in the human genome is that they could be responsible for traits, specific to our species, of which the human brain is the most remarkable. As became obvious, the vast majority of human-specific single nucleotide substitutions are found in noncoding, likely regulatory regions. A number of genes, associated with these human-specific alleles, often through novel enhancer activity, were in fact shown to be implicated in human-specific development of certain brain areas, including the prefrontal cortex. Human-specific deletions may remove regulatory sequences, such as enhancers. Segmental duplications, because of their large size, create new coding sequences, like new functional paralogs. Further functional study of these variants will shed light on evolution of our species, as well as on the etiology of neurodevelopmental disorders.

Changing the name of the NBPF/DUF1220 domain to the Olduvai domain

We are jointly proposing a new name for a protein domain of approximately 65 amino acids that has been previously termed NBPF or DUF1220. Our two labs independently reported the initial studies of this domain, which is encoded almost entirely within a single gene family. The name Neuroblastoma Breakpoint Family ( NBPF) was applied to this gene family when the first identified member of the family was found to be interrupted in an individual with neuroblastoma.

Prior to this discovery, the Pfam database had termed the domain DUF1220, denoting it as one of many protein domains of unknown function. It has been Pfam’s intention to use “DUF” nomenclature to serve only as a temporary placeholder until more appropriate names are proposed based on research findings.

We believe that additional studies of this domain, primarily from our laboratories over the past 10 years, have resulted in furthering our understanding of these sequences to the point where proposing a new name for this domain is warranted. Because of considerable data linking the domain to human-specific evolution, brain expansion and cognition, we believe a name reflecting these findings would be appropriate. With this in mind, we have chosen to name the domain (and the repeat that encodes it) Olduvai. The gene family will remain as NBPF for now. The primary domain subtypes will retain their previously assigned names (e.g. CON1-3; HLS1-3), and the three-domain block that expanded dramatically in the human lineage will be termed the Olduvai triplet.

The new name refers to Olduvai Gorge, which is a site in East Africa that has been the source of major anthropological discoveries in the early-mid 1900’s. We also chose the name as a tribute to the scientists who made important contributions to the early studies of human origins and our African genesis.

Enhancing our brains: Genomic mechanisms underlying cortical evolution

Our most distinguishing higher cognitive functions are controlled by the cerebral cortex. Comparative studies detail abundant anatomical and cellular features unique to the human developing and adult neocortex. Emerging genomic studies have further defined vast differences distinguishing developing human neocortices from related primates. These human-specific changes can affect gene function and/or expression, and result from structural variations such as chromosomal deletions and duplications, or from point mutations in coding and noncoding regulatory regions. Here, we review this rapidly growing field which aims to identify and characterize genetic loci unique to the human cerebral cortex. We catalog known human-specific genomic changes distinct from other primates, including those whose function has been interrogated in animal models. We also discuss how new model systems and technologies such as single cell RNA sequencing, primate iPSCs, and gene editing, are enabling the field to gain unprecedented resolution into function of these human-specific changes. Some neurological disorders are thought to uniquely present in humans, thus reinforcing the need to comprehensively understand human-specific gene expression in the developing brain.

High resolution measurement of DUF1220 domain copy number from whole genome sequence data

Background

DUF1220 protein domains found primarily in Neuroblastoma BreakPoint Family (NBPF) genes show the greatest human lineage-specific increase in copy number of any coding region in the genome. There are 302 haploid copies of DUF1220 in hg38 (~160 of which are human-specific) and the majority of these can be divided into 6 different subtypes (referred to as clades). Copy number changes of specific DUF1220 clades have been associated in a dose-dependent manner with brain size variation (both evolutionarily and within the human population), cognitive aptitude, autism severity, and schizophrenia severity. However, no published methods can directly measure copies of DUF1220 with high accuracy and no method can distinguish between domains within a clade.

Results

Here we describe a novel method for measuring copies of DUF1220 domains and the NBPF genes in which they are found from whole genome sequence data. We have characterized the effect that various sequencing and alignment parameters and strategies have on the accuracy and precision of the method and defined the parameters that lead to optimal DUF1220 copy number measurement and resolution. We show that copy number estimates obtained using our read depth approach are highly correlated with those generated by ddPCR for three representative DUF1220 clades. By simulation, we demonstrate that our method provides sufficient resolution to analyze DUF1220 copy number variation at three levels: (1) DUF1220 clade copy number within individual genes and groups of genes (gene-specific clade groups) (2) genome wide DUF1220 clade copies and (3) gene copy number for DUF1220-encoding genes.

Conclusions

To our knowledge, this is the first method to accurately measure copies of all six DUF1220 clades and the first method to provide gene specific resolution of these clades. This allows one to discriminate among the ~300 haploid human DUF1220 copies to an extent not possible with any other method. The result is a greatly enhanced capability to analyze the role that these sequences play in human variation and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3976-z) contains supplementary material, which is available to authorized users.

Evolution of the Human Nervous System Function, Structure, and Development

The nervous system-in particular, the brain and its cognitive abilities-is among humans' most distinctive and impressive attributes. How the nervous system has changed in the human lineage and how it differs from that of closely related primates is not well understood. Here, we consider recent comparative analyses of extant species that are uncovering new evidence for evolutionary changes in the size and the number of neurons in the human nervous system, as well as the cellular and molecular reorganization of its neural circuits. We also discuss the developmental mechanisms and underlying genetic and molecular changes that generate these structural and functional differences. As relevant new information and tools materialize at an unprecedented pace, the field is now ripe for systematic and functionally relevant studies of the development and evolution of human nervous system specializations.

The birth of a human-specific neural gene by incomplete duplication and gene fusion

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.

Human adaptation and evolution by segmental duplication

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.