DUF1220-domain copy number implicated in human brain-size pathology and evolution

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.

A genetic-epigenetic interplay at 1q21.1 locus underlies CHD1L-mediated vulnerability to primary progressive multiple sclerosis

Multiple Sclerosis (MS) is a heterogeneous inflammatory and neurodegenerative disease with an unpredictable course towards progressive disability. Treating progressive MS is challenging due to limited insights into the underlying mechanisms. We examined the molecular changes associated with primary progressive MS (PPMS) using a cross-tissue (blood and post-mortem brain) and multilayered data (genetic, epigenetic, transcriptomic) from independent cohorts. In PPMS, we found hypermethylation of the 1q21.1 locus, controlled by PPMS-specific genetic variations and influencing the expression of proximal genes (CHD1L, PRKAB2) in the brain. Evidence from reporter assay and CRISPR/dCas9 experiments supports a causal link between methylation and expression and correlation network analysis further implicates these genes in PPMS brain processes. Knock-down of CHD1L in human iPSC-derived neurons and knock-out of chd1l in zebrafish led to developmental and functional deficits of neurons. Thus, several lines of evidence suggest a distinct genetic-epigenetic-transcriptional interplay in the 1q21.1 locus potentially contributing to PPMS pathogenesis.

Copy number variations in autistic children

Autism spectrum disorder (ASD) manifests as a neurodevelopmental condition marked by challenges in social communication, interaction and the performing of repetitive behaviors. The prevalence of autism increases markedly on an annual basis; however, the etiology remains incompletely understood. Cytogenetically visible chromosomal abnormalities, including copy number variations (CNVs), have been shown to contribute to the pathogenesis of ASD. More than 1% of ASD conditions can be explained based on a known genetic locus, whereas CNVs account for 5-10% of cases. However, there are no studies on the Saudi Arabian population for the detection of CNVs linked to ASD, to the best of our knowledge. Therefore, the aim of the present study was to explore the prevalence of CNVs in autistic Saudi Arabian children. Genomic DNA was extracted from the peripheral blood of 14 autistic children along with four healthy control children and then array-based comparative genomic hybridization (aCGH) was used to detect CNVs. Bioinformatics analysis of the aCGH results showed the presence of recurrent and non-recurrent deletion/duplication CNVs in several regions of the genome of autistic children. The most frequent CNVs were 1q21.2, 3p26.3, 4q13.2, 6p25.3, 6q24.2, 7p21.1, 7q34, 7q11.1, 8p23.2, 13q32.3, 14q11.1-q11.2 and 15q11.1-q11.2. In the present study, CNVs in autistic Saudi Arabian children were identified to improve the understanding of the etiology of autism and facilitate its diagnosis. Additionally, the present study identified certain possible pathogenic genes in the CNV region associated with several developmental and neurogenetic diseases.

Mutation and evolution: Conceptual possibilities

Although random mutation is central to models of evolutionary change, a lack of clarity remains regarding the conceptual possibilities for thinking about the nature and role of mutation in evolution. We distinguish several claims at the intersection of mutation, evolution, and directionality and then characterize a previously unrecognized category: complex conditioned mutation. Empirical evidence in support of this category suggests that the historically famous fluctuation test should be revisited, and new experiments should be undertaken with emerging experimental techniques to facilitate detecting mutation rates within specific loci at an ultra-high, individual base pair resolution.

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.

Evolutionary honing in and mutational replacement: how long-term directed mutational responses to specific environmental pressures are possible

Recent results have shown that the human malaria-resistant hemoglobin S mutation originates de novo more frequently in the gene and in the population where it is of adaptive significance, namely, in the hemoglobin subunit beta gene compared to the nonresistant but otherwise identical 20A[Formula: see text]T mutation in the hemoglobin subunit delta gene, and in sub-Saharan Africans, who have been subject to intense malarial pressure for many generations, compared to northern Europeans, who have not. This finding raises a fundamental challenge to the traditional notion of accidental mutation. Here, we address this finding with the replacement hypothesis, according to which preexisting genetic interactions can lead directly and mechanistically to mutations that simplify and replace them. Thus, an evolutionary process under selection can gradually hone in on interactions of importance for the currently evolving adaptations, from which large-effect mutations follow that are relevant to these adaptations. We exemplify this hypothesis using multiple types of mutation, including gene fusion mutations, gene duplication mutations, A[Formula: see text]G mutations in RNA-edited sites and transcription-associated mutations, and place it in the broader context of a system-level view of mutation origination called interaction-based evolution. Potential consequences include that similarity of mutation pressures may contribute to parallel evolution in genetically related species, that the evolution of genome organization may be driven by mutational mechanisms, that transposable element movements may also be explained by replacement, and that long-term directed mutational responses to specific environmental pressures are possible. Such mutational phenomena need to be further tested by future studies in natural and artificial settings.

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.

Identification of 1q21.1 microduplication in a family: A case report

BACKGROUND

Copy number variation (CNV) has become widely recognized in recent years due to the extensive use of gene screening in developmental disorders and epilepsy research. 1q21.1 microduplication syndrome is a rare CNV disease that can manifest as multiple congenital developmental disorders, autism spectrum disorders, congenital malformations, and congenital heart defects with genetic heterogeneity.

CASE SUMMARY

We reported a pediatric patient with 1q21.1 microduplication syndrome, and carried out a literature review to determine the correlation between 1q21.1 microduplication and its phenotypes. We summarized the patient’s medical history and clinical symptoms, and extracted genomic DNA from the patient, her parents, elder brother, and sister. The patient was an 8-mo-old girl who was hospitalized for recurrent convulsions over a 2-mo period. Whole exon sequencing and whole genome low-depth sequencing (CNV-seq) were then performed. Whole exon sequencing detected a 1.58-Mb duplication in the CHR1:145883867-147465312 region, which was located in the 1q21.1 region. Family analysis showed that the pathogenetic duplication fragment, which was also detected in her elder brother’s DNA originated from the mother.

CONCLUSION

Whole exon sequencing combined with quantitative polymerase chain reaction can provide an accurate molecular diagnosis in children with 1q21.1 microduplication syndrome, which is of great significance for genetic counseling and early intervention.

Molecular underpinnings of exceptional response in primary malignant melanoma of the esophagus to anti-PD-1 monotherapy

BACKGROUND

Accumulating data suggest that mucosal melanoma, well known for its poor response to immune checkpoint blockade (ICB) and abysmal prognosis, is a heterogeneous subtype of melanoma with distinct genomic and clinical characteristics between different anatomic locations of the primary lesions. Primary malignant melanoma of the esophagus (PMME) is a rare, highly aggressive disease with a poorer prognosis compared with that of non-esophageal mucosal melanoma (NEMM). In this study, we retrospectively analyzed the efficacy of anti-programmed death (PD)-1 in patients with PMME and explored its molecular basis.

METHODS

The response and survival of patients with PMME and NEMM under anti-PD-1 monotherapy were retrospectively analyzed. To explore the molecular mechanisms of the difference in therapeutic efficacy between PMME and NEMM, we performed genomic analysis, bulk RNA sequencing, and multiplex immunohistochemistry staining.

RESULTS

We found that PMME (n=28) responded better to anti-PD-1 treatment than NEMM (n=64), with a significantly higher objective response rate (33.3% (95% CI 14.3% to 52.3%) vs 6.6% (95% CI 0.2% to 12.9%)) and disease control rate (74.1% (95% CI 56.4% to 91.7%) vs 37.7% (95% CI 25.2% to 50.2%)). Genomic sequencing analysis revealed that the genomic aberration landscape of PMME predominated in classical cancer driver genes, with approximately half of PMME cases harboring mutations in BRAF, N/KRAS, and NF1. In contrast, most NEMM cases were triple wild-type. Transcriptome analysis revealed that, compared with NEMM, PMME displayed more significant proliferation and inflammatory features with higher expression of genes related to antigen presentation and differentiation, and a less immunosuppressive signature with lower expression of inhibitory immune checkpoints and dedifferentiation-related genes. The multiplex immunohistochemical analysis also demonstrated higher CD8⁺ T-cell infiltration in PMME than in NEMM.

CONCLUSIONS

PMME is an outlier of mucosal melanoma showing a malicious phenotype but a particularly high response rate to ICB because of its distinct molecular characteristics. Patient stratification based on anatomic origin can facilitate clinical decision-making in patients with mucosal melanoma following the verification of our results in future prospective studies.

Association Study between DUF1220 Copy Number and Severity of Social Impairment in Sex-balanced Simplex Cases of Autism

Introduction

Copy number variations (CNVs), which are genetic factors responsible for human evolution, have emerged as underlying pathogenic factors for a number of diseases such as autism spectrum disorders (ASD). DUF1220 coding sequences have been shown to be positively associated with the severity of symptoms in familial/multiplex cases of autism. However, this association has not been confirmed in simplex autism, and the potential impact of gender/sex has not been studied.

Methods

Using saliva samples taken from Iranian children with non-syndromic simplex autism, different ethnicity/race and genetic backgrounds from previous studies, we assessed the association between DUF1220 CNVs and Autism Diagnostic Interview-Revised (ADI-R) domain scores in both males and females.

Results

In the male and female combined group with autism, in line with previous reports, our findings showed that there were no significant associations between DUF1220 CNVs with either total ADI-R score, social, communication, or repetitive diagnostic scores in simplex autism cases. Interestingly, however, in sex classified groups, despite the insignificant results, our findings in girls with autism showed a negative trend between DUF1220 CNVs and severity of symptoms for the social interaction and communication domains. By contrast, in male children with autism, the results showed a positive trend.

Conclusion

It seems that association of DUF1220 CNV with the severity of symptoms in simplex children with autism may follow a sexually dimorphic pattern that needs to be re-examined in prospective studies.

Adaptive sequence divergence forged new neurodevelopmental enhancers in humans

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

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.

A rigorous in silico genomic interrogation at 1p13.3 reveals 16 autosomal dominant candidate genes in syndromic neurodevelopmental disorders

Genome-wide chromosomal microarray is extensively used to detect copy number variations (CNVs), which can diagnose microdeletion and microduplication syndromes. These small unbalanced chromosomal structural rearrangements ranging from 1 kb to 10 Mb comprise up to 15% of human mutations leading to monogenic or contiguous genomic disorders. Albeit rare, CNVs at 1p13.3 cause a variety of neurodevelopmental disorders (NDDs) including development delay (DD), intellectual disability (ID), autism, epilepsy, and craniofacial anomalies (CFA). Most of the 1p13.3 CNV cases reported in the pre-microarray era encompassed a large number of genes and lacked the demarcating genomic coordinates, hampering the discovery of positional candidate genes within the boundaries. In this study, we present four subjects with 1p13.3 microdeletions displaying DD, ID, autism, epilepsy, and CFA. In silico comparative genomic mapping with three previously reported subjects with CNVs and 22 unreported DECIPHER CNV cases has resulted in the identification of four different sub-genomic loci harboring five positional candidate genes for DD, ID, and CFA at 1p13.3. Most of these genes have pathogenic variants reported, and their interacting genes are involved in NDDs. RT-qPCR in various human tissues revealed a high expression pattern in the brain and fetal brain, supporting their functional roles in NDDs. Interrogation of variant databases and interacting protein partners led to the identification of another set of 11 potential candidate genes, which might have been dysregulated by the position effect of these CNVs at 1p13.3. Our studies define 1p13.3 as a genomic region harboring 16 NDD candidate genes and underscore the critical roles of small CNVs in in silico comparative genomic mapping for disease gene discovery. Our candidate genes will help accelerate the isolation of pathogenic heterozygous variants from exome/genome sequencing (ES/GS) databases.

cAMP Signaling in Cancer: A PKA-CREB and EPAC-Centric Approach

Cancer is one of the most common causes of death globally. Despite extensive research and considerable advances in cancer therapy, the fundamentals of the disease remain unclear. Understanding the key signaling mechanisms that cause cancer cell malignancy may help to uncover new pharmaco-targets. Cyclic adenosine monophosphate (cAMP) regulates various biological functions, including those in malignant cells. Understanding intracellular second messenger pathways is crucial for identifying downstream proteins involved in cancer growth and development. cAMP regulates cell signaling and a variety of physiological and pathological activities. There may be an impact on gene transcription from protein kinase A (PKA) as well as its downstream effectors, such as cAMP response element-binding protein (CREB). The position of CREB downstream of numerous growth signaling pathways implies its oncogenic potential in tumor cells. Tumor growth is associated with increased CREB expression and activation. PKA can be used as both an onco-drug target and a biomarker to find, identify, and stage tumors. Exploring cAMP effectors and their downstream pathways in cancer has become easier using exchange protein directly activated by cAMP (EPAC) modulators. This signaling system may inhibit or accelerate tumor growth depending on the tumor and its environment. As cAMP and its effectors are critical for cancer development, targeting them may be a useful cancer treatment strategy. Moreover, by reviewing the material from a distinct viewpoint, this review aims to give a knowledge of the impact of the cAMP signaling pathway and the related effectors on cancer incidence and development. These innovative insights seek to encourage the development of novel treatment techniques and new approaches.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

PDE4DIP in health and diseases

Cyclic-AMP (cAMP), the first second messenger to be identified, is synthesized, and is universally utilized as a second messenger, and plays important roles in integrity, and function of organs, including heart. Through its coupling with other intracellular messengers, cAMP facilitates excitation-contraction coupling, increases heart rate and conduction velocity. It is degraded by a class of enzymes called cAMP-dependent phosphodiesterase (PDE), with PDE3 and PDE4 being the predominant isoforms in the heart. This highly diverse class of enzymes degrade cAMP and through anchoring proteins generates dynamic microdomains to target specific proteins and control specific cell functions in response to various stimuli. The impaired function of the anchoring protein either by inherited genetic mutations or acquired injuries results in altered intracellular targeting, and blunted responsiveness to stimulating pathways and contributes to pathological cardiac remodeling, cardiac arrhythmias and reduced cell survival. Recent genetic studies provide compelling evidence for an association between the variants in the anchoring protein PDE4DIP and atrial fibrillation, stroke, and heart failure.

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.

De novo mutation rates at the single-mutation resolution in a human HBB gene-region associated with adaptation and genetic disease

Although it is known that the mutation rate varies across the genome, previous estimates were based on averaging across various numbers of positions. Here, we describe a method to measure the origination rates of target mutations at target base positions and apply it to a 6-bp region in the human hemoglobin subunit beta (HBB) gene and to the identical, paralogous hemoglobin subunit delta (HBD) region in sperm cells from both African and European donors. The HBB region of interest (ROI) includes the site of the hemoglobin S (HbS) mutation, which protects against malaria, is common in Africa, and has served as a classic example of adaptation by random mutation and natural selection. We found a significant correspondence between de novo mutation rates and past observations of alleles in carriers, showing that mutation rates vary substantially in a mutation-specific manner that contributes to the site frequency spectrum. We also found that the overall point mutation rate is significantly higher in Africans than in Europeans in the HBB region studied. Finally, the rate of the 20A→T mutation, called the “HbS mutation” when it appears in HBB, is significantly higher than expected from the genome-wide average for this mutation type. Nine instances were observed in the African HBB ROI, where it is of adaptive significance, representing at least three independent originations; no instances were observed elsewhere. Further studies will be needed to examine mutation rates at the single-mutation resolution across these and other loci and organisms and to uncover the molecular mechanisms responsible.

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.

Global gene expression analysis of systemic sclerosis myofibroblasts demonstrates a marked increase in the expression of multiple NBPF genes

Myofibroblasts are the key effector cells responsible for the exaggerated tissue fibrosis in Systemic Sclerosis (SSc). Despite their importance to SSc pathogenesis, the specific transcriptome of SSc myofibroblasts has not been described. The purpose of this study was to identify transcriptome differences between SSc myofibroblasts and non-myofibroblastic cells. Alpha smooth muscle actin (α-SMA) expressing myofibroblasts and α-SMA negative cells were isolated employing laser capture microdissection from dermal cell cultures from four patients with diffuse SSc of recent onset. Total mRNA was extracted from both cell populations, amplified and analyzed employing microarrays. Results for specific genes were validated by Western blots and by immunohistochemistry. Transcriptome analysis revealed 97 differentially expressed transcripts in SSc myofibroblasts compared with non-myofibroblasts. Annotation clustering of the SSc myofibroblast-specific transcripts failed to show a TGF-β signature. The most represented transcripts corresponded to several different genes from the Neuroblastoma Breakpoint Family (NBPF) of genes. NBPF genes are highly expanded in humans but are not present in murine or rat genomes. In vitro studies employing cultured SSc dermal fibroblasts and immunohistochemistry of affected SSc skin confirmed increased NBPF expression in SSc. These results indicate that SSc myofibroblasts represent a unique cell lineage expressing a specific transcriptome that includes very high levels of transcripts corresponding to numerous NBPF genes. Elevated expression of NBPF genes in SSc myofibroblasts suggests that NBPF gene products may play a role in SSc pathogenesis and may represent a novel therapeutic target.

22q11.2 Low Copy Repeats Expanded in the Human Lineage

Segmental duplications or low copy repeats (LCRs) constitute duplicated regions interspersed in the human genome, currently neglected in standard analyses due to their extreme complexity. Recent functional studies have indicated the potential of genes within LCRs in synaptogenesis, neuronal migration, and neocortical expansion in the human lineage. One of the regions with the highest proportion of duplicated sequence is the 22q11.2 locus, carrying eight LCRs (LCR22-A until LCR22-H), and rearrangements between them cause the 22q11.2 deletion syndrome. The LCR22-A block was recently reported to be hypervariable in the human population. It remains unknown whether this variability also exists in non-human primates, since research is strongly hampered by the presence of sequence gaps in the human and non-human primate reference genomes. To chart the LCR22 haplotypes and the associated inter- and intra-species variability, we de novo assembled the region in non-human primates by a combination of optical mapping techniques. A minimal and likely ancient haplotype is present in the chimpanzee, bonobo, and rhesus monkey without intra-species variation. In addition, the optical maps identified assembly errors and closed gaps in the orthologous chromosome 22 reference sequences. These findings indicate the LCR22 expansion to be unique to the human population, which might indicate involvement of the region in human evolution and adaptation. Those maps will enable LCR22-specific functional studies and investigate potential associations with the phenotypic variability in the 22q11.2 deletion syndrome.

Modelling segmental duplications in the human genome

Background

Segmental duplications (SDs) are long DNA sequences that are repeated in a genome and have high sequence identity. In contrast to repetitive elements they are often unique and only sometimes have multiple copies in a genome. There are several well-studied mechanisms responsible for segmental duplications: non-allelic homologous recombination, non-homologous end joining and replication slippage. Such duplications play an important role in evolution, however, we do not have a full understanding of the dynamic properties of the duplication process.

Results

We study segmental duplications through a graph representation where nodes represent genomic regions and edges represent duplications between them. The resulting network (the SD network) is quite complex and has distinct features which allow us to make inference on the evolution of segmantal duplications. We come up with the network growth model that explains features of the SD network thus giving us insights on dynamics of segmental duplications in the human genome. Based on our analysis of genomes of other species the network growth model seems to be applicable for multiple mammalian genomes.

Conclusions

Our analysis suggests that duplication rates of genomic loci grow linearly with the number of copies of a duplicated region. Several scenarios explaining such a preferential duplication rates were suggested.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07789-7).

Dissecting Molecular Genetic Mechanisms of 1q21.1 CNV in Neuropsychiatric Disorders

Pathogenic copy number variations (CNVs) contribute to the etiology of neurodevelopmental/neuropsychiatric disorders (NDs). Increased CNV burden has been found to be critically involved in NDs compared with controls in clinical studies. The 1q21.1 CNVs, rare and large chromosomal microduplications and microdeletions, are detected in many patients with NDs. Phenotypes of duplication and deletion appear at the two ends of the spectrum. Microdeletions are predominant in individuals with schizophrenia (SCZ) and microcephaly, whereas microduplications are predominant in individuals with autism spectrum disorder (ASD) and macrocephaly. However, its complexity hinders the discovery of molecular pathways and phenotypic networks. In this review, we summarize the recent genome-wide association studies (GWASs) that have identified candidate genes positively correlated with 1q21.1 CNVs, which are likely to contribute to abnormal phenotypes in carriers. We discuss the clinical data implicated in the 1q21.1 genetic structure that is strongly associated with neurodevelopmental dysfunctions like cognitive impairment and reduced synaptic plasticity. We further present variations reported in the phenotypic severity, genomic penetrance and inheritance.

Chromatin accessibility in neuropsychiatric disorders

Epigenetic mechanisms have recently emerged as critical regulators of brain function in health and disease. By controlling the accessibility and the expression of specific genes, these pathways can mediate transient and long-lasting changes in neuronal function in both physiological and pathological contexts. Due to the extreme complexity of the epigenetic regulatory landscape, emerging methods that directly assay chromatin accessibility are of particular interest. Here, I review recent insights gained on open and closed chromatin states in the brain, with emphasis on neuropsychiatric disorders. These advances generated an invaluable wealth of information that can help us better understand gene regulation in the brain and its impairments that contribute to the development of disease.

1q21.1 distal copy number variants are associated with cerebral and cognitive alterations in humans

Low-frequency 1q21.1 distal deletion and duplication copy number variant (CNV) carriers are predisposed to multiple neurodevelopmental disorders, including schizophrenia, autism and intellectual disability. Human carriers display a high prevalence of micro- and macrocephaly in deletion and duplication carriers, respectively. The underlying brain structural diversity remains largely unknown. We systematically called CNVs in 38 cohorts from the large-scale ENIGMA-CNV collaboration and the UK Biobank and identified 28 1q21.1 distal deletion and 22 duplication carriers and 37,088 non-carriers (48% male) derived from 15 distinct magnetic resonance imaging scanner sites. With standardized methods, we compared subcortical and cortical brain measures (all) and cognitive performance (UK Biobank only) between carrier groups also testing for mediation of brain structure on cognition. We identified positive dosage effects of copy number on intracranial volume (ICV) and total cortical surface area, with the largest effects in frontal and cingulate cortices, and negative dosage effects on caudate and hippocampal volumes. The carriers displayed distinct cognitive deficit profiles in cognitive tasks from the UK Biobank with intermediate decreases in duplication carriers and somewhat larger in deletion carriers-the latter potentially mediated by ICV or cortical surface area. These results shed light on pathobiological mechanisms of neurodevelopmental disorders, by demonstrating gene dose effect on specific brain structures and effect on cognitive function.

The psychiatric phenotypes of 1q21 distal deletion and duplication

Copy number variants are amongst the most highly penetrant risk factors for psychopathology and neurodevelopmental deficits, but little information about the detailed clinical phenotype associated with particular variants is available. We present the largest study of the microdeletion and -duplication at the distal 1q21 locus, which has been associated with schizophrenia and intellectual disability, in order to investigate the range of psychiatric phenotypes. Clinical and cognitive data from 68 deletion and 55 duplication carriers were analysed with logistic regression analysis to compare frequencies of mental disorders between carrier groups and controls, and linear mixed models to compare quantitative phenotypes. Both children and adults with copy number variants at 1q21 had high frequencies of psychopathology. In the children, neurodevelopmental disorders were most prominent (56% for deletion, 68% for duplication carriers). Adults had increased prevalence of mood (35% for deletion [OR = 6.6 (95% CI: 1.4-40.1)], 55% for duplication carriers [8.3 (1.4-55.5)]) and anxiety disorders (24% [1.8 (0.4-8.4)] and 55% [10.0 (1.9-71.2)]). The adult group, which included mainly genetically affected parents of probands, had an IQ in the normal range. These results confirm high prevalence of neurodevelopmental disorders associated with CNVs at 1q21 but also reveal high prevalence of mood and anxiety disorders in a high-functioning adult group with these CNVs. Because carriers of neurodevelopmental CNVs who show relevant psychopathology but no major cognitive impairment are not currently routinely receiving clinical genetic services widening of genetic testing in psychiatry may be considered.

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.

Disorders Associated With Diverse, Recurrent Deletions and Duplications at 1q21.1

The subchromosomal region 1q21.1 is one of the hotspots in the human genome for deletions and reciprocal duplications, owing to the existence of hundreds of segmental duplications. Recurrent deletions and duplications in this region are thought to be causative in patients with variable clinical manifestations. Based on the genomic locations, deletions and duplications at the 1q21.1 locus have been associated with distinguishable syndromes: chromosome 1q21.1 deletion syndrome, chromosome 1q21.1 duplication syndrome, and thrombocytopenia-absent radius (TAR) syndrome, which is partially due to deletions at the proximal 1q21.1 region. We report here diverse, recurrent deletions and duplications at the 1q21.1 locus in 36 patients from a cohort of 5,200 individuals. Among the 36 patients, 18 patients carry 1q21.1 deletions, nine individuals have reciprocal duplications at 1q21.1, two patients share an identical short deletion, and the remaining seven possess variable sizes of duplications at the proximal 1q21.1 region. Furthermore, we provide cytogenetic characterization and detailed clinical features for each patient. Notably, duplications at the proximal 1q21.1 region have not been associated with a defined disorder in publications. However, recurrent duplications at the proximal 1q21.1 region among the seven patients strongly suggested that the variants are likely pathogenic. The common phenotypical features of those disorders are also summarized to facilitate clinical diagnoses and genetic counseling.

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.

On the origin of schizophrenia: Testing evolutionary theories in the post-genomic era

Considering the relatively high heritability of schizophrenia and the fact that it significantly reduces the reproductive fitness of affected individuals, it is not clear how the disorder is still maintained in human populations at a disproportionally high prevalence. Many theories propose that the disorder is a result of a trade-off between costs and benefits of the evolution of exclusively human adaptations. There have also been suggestions that schizophrenia risk alleles are accompanied with increase in fitness of affected persons or their relatives in both past and current social contexts. The discoveries of novel schizophrenia-related genes and the advancements in comparative genomics (especially comparisons of the human genome and the genomes of related species, such as chimpanzees and extinct hominids) have finally made certain evolutionary theories testable. In this paper, we review the current understanding of the genetics of schizophrenia, the basic principles of evolution that complement our understanding of the subject, and the latest genetic studies that examine long-standing evolutionary theories of schizophrenia using novel methodologies and data. We find that the origin of schizophrenia is complex and likely governed by different evolutionary mechanisms that are not mutually exclusive. Furthermore, the most recent evidence implies that schizophrenia cannot be comprehended as a trait that has elevated fitness in human evolutionary lineage, but has been a mildly deleterious by-product of specific patterns of the evolution of the human brain. In other words, novel findings do not support previous hypotheses stating that schizophrenia risk genes have an evolutionary advantage.

Mutations in ANKLE2, a ZIKA Virus Target, Disrupt an Asymmetric Cell Division Pathway in Drosophila Neuroblasts to Cause Microcephaly

The apical Par complex, which contains atypical protein kinase C (aPKC), Bazooka (Par-3), and Par-6, is required for establishing polarity during asymmetric division of neuroblasts in Drosophila, and its activity depends on L(2)gl. We show that loss of Ankle2, a protein associated with microcephaly in humans and known to interact with Zika protein NS4A, reduces brain volume in flies and impacts the function of the Par complex. Reducing Ankle2 levels disrupts endoplasmic reticulum (ER) and nuclear envelope morphology, releasing the kinase Ballchen-VRK1 into the cytosol. These defects are associated with reduced phosphorylation of aPKC, disruption of Par-complex localization, and spindle alignment defects. Importantly, removal of one copy of ballchen or l(2)gl suppresses Ankle2 mutant phenotypes and restores viability and brain size. Human mutational studies implicate the above-mentioned genes in microcephaly and motor neuron disease. We suggest that NS4A, ANKLE2, VRK1, and LLGL1 define a pathway impinging on asymmetric determinants of neural stem cell division.

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.

A catalog of single nucleotide changes distinguishing modern humans from archaic hominins

Throughout the past decade, studying ancient genomes has provided unique insights into human prehistory, and differences between modern humans and other branches like Neanderthals can enrich our understanding of the molecular basis of unique modern human traits. Modern human variation and the interactions between different hominin lineages are now well studied, making it reasonable to go beyond fixed genetic changes and explore changes that are observed at high frequency in present-day humans. Here, we identify 571 genes with non-synonymous changes at high frequency. We suggest that molecular mechanisms in cell division and networks affecting cellular features of neurons were prominently modified by these changes. Complex phenotypes in brain growth trajectory and cognitive traits are likely influenced by these networks and other non-coding changes presented here. We propose that at least some of these changes contributed to uniquely human traits, and should be prioritized for experimental validation.

An eight-case 1q21 region series: novel aberrations and clinical variability with new features

BACKGROUND

Rearrangement of the 1q21 region of chromosome 1 manifests as multiple phenotypes, including microcephaly, intellectual disability, dysmorphic facial features, eye abnormalities, cardiac defects, genitourinary anomalies, autism spectrum disorder, psychiatric conditions and seizures. Herein, we describe eight patients with 1q21 deletion and duplication syndromes, and novel deletions and findings.

METHODS

Chromosomal microarray analysis was performed to identify the existence of copy number variation. Quantitative polymerase chain reaction was applied using specific primers for the control and 1q21 region of chromosome 1. Mutational analysis was performed in case 5 using direct genomic sequencing for exons 1-6 in RBM8A.

RESULTS

Copy number variation analysis identified seven deletions and one duplication of the 1q21 region in the eight patients. In addition, four variations were de novo, and two deletions are reported here for the first time. One of the cases (case 7) presents moderate intellectual disability and dysmorphic facial findings, whereas chromosomal microarray analysis showed that case 7 had an 889-kb deletion in the 1q21 proximal region (GPR89A, PDZK1, CD160, POLR3C and NBPF12).

CONCLUSION

Although the deletion in case 5 did not include the thrombocytopenia-absent radius syndrome critical region or the RBM8A gene, he had pectoral muscle hypoplasia, radius and humerus hypoplasia and short curved ribs, which are indicative of a potential thrombocytopenia-absent radius region modifier. The findings in case 7 suggest that the proximal part of the 1q21 microdeletion syndrome region might be very important for the onset of clinical manifestations. Some novel findings were observed in the presented cases, such as radius and humerus hypoplasia and brain stem hypoplasia. The presented findings expand the spectrum of 1q21 aberrations and provide evidence of genotype-phenotype correlations for this region.

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.

De novo emergence and potential function of human-specific tandem repeats in brain-related loci

Tandem repeats (TRs) are widespread in the genomes of all living organisms. In eukaryotes, they are found in both coding and noncoding regions and have potential roles in the regulation of cellular processes such as transcription, translation and in the modification of protein structure. Recent studies have highlighted TRs as a key regulator of gene expression and a potential contributor to human evolution. Thus, TRs are emerging as an important source of variation that can result in differential gene expression at intra- and inter-species levels. In this study, we performed a genome-wide survey to identify TRs that have emerged in the human lineage. We further examined these loci to explore their potential functional significance for human evolution. We identified 152 human-specific TR (HSTR) loci containing a repeat unit of more than ten bases, with most of them showing a repeat count of two. Gene set enrichment analysis showed that HSTR-associated genes were associated with biological functions in brain development and synapse function. In addition, we compared gene expression of human HSTR loci with orthologues from non-human primates (NHP) in seven different tissues. Strikingly, the expression level of HSTR-associated genes in brain tissues was significantly higher in human than in NHP. These results suggest the possibility that de novo emergence of TRs could have resulted in altered gene expression in humans within a short-time frame and contributed to the rapid evolution of human brain function.

Transcriptional fates of human-specific segmental duplications in brain

Despite the importance of duplicate genes for evolutionary adaptation, accurate gene annotation is often incomplete, incorrect, or lacking in regions of segmental duplication. We developed an approach combining long-read sequencing and hybridization capture to yield full-length transcript information and confidently distinguish between nearly identical genes/paralogs. We used biotinylated probes to enrich for full-length cDNA from duplicated regions, which were then amplified, size-fractionated, and sequenced using single-molecule, long-read sequencing technology, permitting us to distinguish between highly identical genes by virtue of multiple paralogous sequence variants. We examined 19 gene families as expressed in developing and adult human brain, selected for their high sequence identity (average >99%) and overlap with human-specific segmental duplications (SDs). We characterized the transcriptional differences between related paralogs to better understand the birth-death process of duplicate genes and particularly how the process leads to gene innovation. In 48% of the cases, we find that the expressed duplicates have changed substantially from their ancestral models due to novel sites of transcription initiation, splicing, and polyadenylation, as well as fusion transcripts that connect duplication-derived exons with neighboring genes. We detect unannotated open reading frames in genes currently annotated as pseudogenes, while relegating other duplicates to nonfunctional status. Our method significantly improves gene annotation, specifically defining full-length transcripts, isoforms, and open reading frames for new genes in highly identical SDs. The approach will be more broadly applicable to genes in structurally complex regions of other genomes where the duplication process creates novel genes important for adaptive traits.

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.

iCopyDAV: Integrated platform for copy number variations-Detection, annotation and visualization

Discovery of copy number variations (CNVs), a major category of structural variations, have dramatically changed our understanding of differences between individuals and provide an alternate paradigm for the genetic basis of human diseases. CNVs include both copy gain and copy loss events and their detection genome-wide is now possible using high-throughput, low-cost next generation sequencing (NGS) methods. However, accurate detection of CNVs from NGS data is not straightforward due to non-uniform coverage of reads resulting from various systemic biases. We have developed an integrated platform, iCopyDAV, to handle some of these issues in CNV detection in whole genome NGS data. It has a modular framework comprising five major modules: data pre-treatment, segmentation, variant calling, annotation and visualization. An important feature of iCopyDAV is the functional annotation module that enables the user to identify and prioritize CNVs encompassing various functional elements, genomic features and disease-associations. Parallelization of the segmentation algorithms makes the iCopyDAV platform even accessible on a desktop. Here we show the effect of sequencing coverage, read length, bin size, data pre-treatment and segmentation approaches on accurate detection of the complete spectrum of CNVs. Performance of iCopyDAV is evaluated on both simulated data and real data for different sequencing depths. It is an open-source integrated pipeline available at https://github.com/vogetihrsh/icopydav and as Docker’s image at http://bioinf.iiit.ac.in/icopydav/.

Hyponatremia: An Unusual Presentation in a Neonate With Chromosome 1q21.1 Deletion Syndrome

Chromosome 1q21.1 deletion syndrome is associated with a wide variety of clinical features including mild to moderate mental retardation, microcephaly, cardiac abnormalities, and cataracts. We report an unusual case of a premature neonate with persistent hyponatremia, markedly elevated plasma arginine vasopressin level (32.7 pg/mL), and clinical findings consistent with the syndrome of inappropriate antidiuretic hormone secretion (SIADH). The patient, who also had microcephaly and dextrocardia, was subsequently diagnosed with chromosome 1q21.1 deletion syndrome. Further evaluation revealed hypothalamic abnormalities, features not previously described with this syndrome. To our knowledge, this is the first report of SIADH associated with congenital hypothalamic anomalies in a neonate with chromosome 1q21.1 deletion syndrome. We also report our experience using tolvaptan, a vasopressin receptor antagonist, in this patient to effectively maintain eunatremia.

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.

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.

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.

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.