Divergent origins and concerted expansion of two segmental duplications on chromosome 16

SVA retrotransposons and a low copy repeat in humans and great apes: a mobile connection

Segmental duplications (SDs) constitute a considerable fraction of primate genomes. They contribute to genetic variation and provide raw material for evolution. Groups of SDs are characterized by the presence of shared core duplicons. One of these core duplicons, low copy repeat (lcr)16a, has been shown to be particularly active in the propagation of interspersed SDs in primates. The underlying mechanisms are, however, only partially understood. Alu short interspersed elements (SINEs) are frequently found at breakpoints and have been implicated in the expansion of SDs.

Detailed analysis of lcr16a-containing SDs shows that the hominid-specific SVA (SINE-R-VNTR-Alu) retrotransposon is an integral component of the core duplicon in Asian and African great apes. In orang-utan, it provides breakpoints and contributes to both interchromosomal and intrachromosomal lcr16a mobility by inter-element recombination. Furthermore, the data suggest that in hominines (human, chimpanzee, gorilla) SVA recombination-mediated integration of a circular intermediate is the founding event of a lineage-specific lcr16a expansion. One of the hominine lcr16a copies displays large flanking direct repeats, a structural feature shared by other SDs in the human genome.

Taken together, the results obtained extend the range of SVAs’ contribution to genome evolution from RNA-mediated transduction to DNA-based recombination. In addition, they provide further support for a role of circular intermediates in SD mobilization.

Submicroscopic aberrations of chromosome 16 in prenatal diagnosis

Background

Nearly 9.89% of chromosome 16 consists of segmental duplications, which makes it prone to non-homologous recombination. The present study aimed to investigate the incidence and perinatal characteristics of submicroscopic chromosome 16 aberrations in prenatal diagnosis.

Results

A total of 2,414 consecutive fetuses that underwent prenatal chromosomal microarray analysis (CMA) between January 2016 and December 2018 were reviewed. Submicroscopic anomalies of chromosome 16 accounted for 11.1% (15/134) of all submicroscopic anomalies detected in fetuses with normal karyotype, which was larger than the percentage of anomalies in any other chromosome. The 15 submicroscopic anomalies of chromosome 16 were identified in 14 cases; 12 of them had ultrasound abnormalities. They were classified as pathogenic (N = 7), and variants of uncertain significance (N = 8). Seven fetuses with variants of uncertain significance were ended in live-born, and the remaining were end in pregnancy termination.

Conclusion

Submicroscopic aberrations of chromosome 16 are frequent findings in prenatal diagnosis, which emphasize the challenge of genetic counseling and the value of CMA. Prenatal diagnosis should lead to long-term monitoring of children with such chromosomal abnormalities for better understanding of the phenotype of chromosome 16 microdeletion and microduplication syndromes.

16p13 microduplication without CREBBP involvement: Moving toward a phenotype delineation

The short arm of chromosome 16 is one of the less stable regions of our genome, as over 10% of the euchromatic region of 16p is composed of highly complex low copy repeats that are known to be predisposed to rearrangements mediated by non-allelic homologous recombination. The 16p13.3p13.13 molecular region has been defined as the 16p duplication hotspot, and duplications of chromosome 16p13 have recently been confirmed to cause a recognizable syndrome, with CREBBP being the main phenotype-causing gene. To date, only one case report is present in the literature with a 16p13 duplication without CREBBP involvement; we describe here a second analogous case with a not previously reported 16p13.2p13.13 microduplication. This paper allows us to better delineate the clinical features of 16p13 microduplications that do not encompass CREBBP and, concurrently, to narrow the molecular region responsible for congenital heart defects in 16p duplications as well as to propose GRIN2A as a candidate gene for epilepsy.

A recurrent 16p12.1 microdeletion supports a two-hit model for severe developmental delay

We report the identification of a recurrent 520-kbp 16p12.1 microdeletion significantly associated with childhood developmental delay. The microdeletion was detected in 20/11,873 cases vs. 2/8,540 controls (p=0.0009, OR=7.2) and replicated in a second series of 22/9,254 cases vs. 6/6,299 controls (p=0.028, OR=2.5). Most deletions were inherited with carrier parents likely to manifest neuropsychiatric phenotypes (p=0.037, OR=6). Probands were more likely to carry an additional large CNV when compared to matched controls (10/42 cases, p=5.7×10^-5, OR=6.65). Clinical features of cases with two mutations were distinct from and/or more severe than clinical features of patients carrying only the co-occurring mutation. Our data suggest a two-hit model in which the 16p12.1 microdeletion both predisposes to neuropsychiatric phenotypes as a single event and exacerbates neurodevelopmental phenotypes in association with other large deletions or duplications. Analysis of other microdeletions with variable expressivity suggests that this two-hit model may be more generally applicable to neuropsychiatric disease.

Primate-specific spliced PMCHL RNAs are non-protein coding in human and macaque tissues

Background

Brain-expressed genes that were created in primate lineage represent obvious candidates to investigate molecular mechanisms that contributed to neural reorganization and emergence of new behavioural functions in Homo sapiens. PMCHL1 arose from retroposition of a pro-melanin-concentrating hormone (PMCH) antisense mRNA on the ancestral human chromosome 5p14 when platyrrhines and catarrhines diverged. Mutations before divergence of hylobatidae led to creation of new exons and finally PMCHL1 duplicated in an ancestor of hominids to generate PMCHL2 at the human chromosome 5q13. A complex pattern of spliced and unspliced PMCHL RNAs were found in human brain and testis.

Results

Several novel spliced PMCHL transcripts have been characterized in human testis and fetal brain, identifying an additional exon and novel splice sites. Sequencing of PMCHL genes in several non-human primates allowed to carry out phylogenetic analyses revealing that the initial retroposition event took place within an intron of the brain cadherin (CDH12) gene, soon after platyrrhine/catarrhine divergence, i.e. 30–35 Mya, and was concomitant with the insertion of an AluSg element. Sequence analysis of the spliced PMCHL transcripts identified only short ORFs of less than 300 bp, with low (VMCH-p8 and protein variants) or no evolutionary conservation. Western blot analyses of human and macaque tissues expressing PMCHL RNA failed to reveal any protein corresponding to VMCH-p8 and protein variants encoded by spliced transcripts.

Conclusion

Our present results improve our knowledge of the gene structure and the evolutionary history of the primate-specific chimeric PMCHL genes. These genes produce multiple spliced transcripts, bearing short, non-conserved and apparently non-translated ORFs that may function as mRNA-like non-coding RNAs.

Hominoid chromosomal rearrangements on 17q map to complex regions of segmental duplication

BACKGROUND

Chromosomal rearrangements, such as translocations and inversions, are recurrent phenomena during evolution, and both of them are involved in reproductive isolation and speciation. To better understand the molecular basis of chromosome rearrangements and their part in karyotype evolution, we have investigated the history of human chromosome 17 by comparative fluorescence in situ hybridization (FISH) and sequence analysis.

RESULTS

Human bacterial artificial chromosome/p1 artificial chromosome probes spanning the length of chromosome 17 were used in FISH experiments on great apes, Old World monkeys and New World monkeys to study the evolutionary history of this chromosome. We observed that the macaque marker order represents the ancestral organization. Human, chimpanzee and gorilla homologous chromosomes differ by a paracentric inversion that occurred specifically in the Homo sapiens/Pan troglodytes/Gorilla gorilla ancestor. Detailed analyses of the paracentric inversion revealed that the breakpoints mapped to two regions syntenic to human 17q12/21 and 17q23, both rich in segmental duplications.

CONCLUSION

Sequence analyses of the human and macaque organization suggest that the duplication events occurred in the catarrhine ancestor with the duplication blocks continuing to duplicate or undergo gene conversion during evolution of the hominoid lineage. We propose that the presence of these duplicons has mediated the inversion in the H. sapiens/P. troglodytes/G. gorilla ancestor. Recently, the same duplication blocks have been shown to be polymorphic in the human population and to be involved in triggering microdeletion and duplication in human. These results further support a model where genomic architecture has a direct role in both rearrangement involved in karyotype evolution and genomic instability in human.

Ancestral reconstruction of segmental duplications reveals punctuated cores of human genome evolution

Human segmental duplications are hotspots for nonallelic homologous recombination leading to genomic disorders, copy-number polymorphisms and gene and transcript innovations. The complex structure and history of these regions have precluded a global evolutionary analysis. Combining a modified A-Bruijn graph algorithm with comparative genome sequence data, we identify the origin of 4,692 ancestral duplication loci and use these to cluster 437 complex duplication blocks into 24 distinct groups. The sequence-divergence data between ancestral-derivative pairs and a comparison with the chimpanzee and macaque genome support a 'punctuated' model of evolution. Our analysis reveals that human segmental duplications are frequently organized around 'core' duplicons, which are enriched for transcripts and, in some cases, encode primate-specific genes undergoing positive selection. We hypothesize that the rapid expansion and fixation of some intrachromosomal segmental duplications during great-ape evolution has been due to the selective advantage conferred by these genes and transcripts embedded within these core duplications.

Who's afraid of Homo sapiens?

Understanding how humans differ from other animals, as well as how we are like them, requires comparative investigations. For the purpose of documenting the distinctive features of humans, the most informative research involves comparing humans to our closest relatives–the chimpanzees and other great apes. Psychology and anthropology have maintained a tradition of empirical comparative research on human specializations of cognition. The neurosciences, by contrast, have been dominated by the model-animal research paradigm, which presupposes the commonality of "basic" features of brain organization across species and discourages serious treatment of species differences. As a result, the neurosciences have made little progress in understanding human brain specializations. Recent developments in neuroimaging, genomics, and other non-invasive techniques make it possible to directly compare humans and nonhuman species at levels of organization that were previously inaccessible, offering the hope of gaining a better understanding of the species-specific features of the human brain. This hope will be dashed, however, if chimpanzees and other great ape species become unavailable for even non-invasive research.

Recurrent duplication-driven transposition of DNA during hominoid evolution

The underlying mechanism by which the interspersed pattern of human segmental duplications has evolved is unknown. Based on a comparative analysis of primate genomes, we show that a particular segmental duplication (LCR16a) has been the source locus for the formation of the majority of intrachromosomal duplications blocks on human chromosome 16. We provide evidence that this particular segment has been active independently in each great ape and human lineage at different points during evolution. Euchromatic sequence that flanks sites of LCR16a integration are frequently lineage-specific duplications. This process has mobilized duplication blocks (15-200 kb in size) to new genomic locations in each species. Breakpoint analysis of lineage-specific insertions suggests coordinated deletion of repeat-rich DNA at the target site, in some cases deleting genes in that species. Our data support a model of duplication where the probability that a segment of DNA becomes duplicated is determined by its proximity to core duplicons, such as LCR16a.

Evolutionary implications of pericentromeric gene expression in humans

Human pericentromeric sequences are enriched for recent sequence duplications. The continual creation and shuffling of these duplications can create novel intron-exon structures and it has been suggested that these regions have a function as gene nurseries. However, these sequences are also rich in satellite repeats which can repress transcription, and analyses of chromosomes 10 and 21 have suggested that they are transcript poor. Here, we investigate the relationship between pericentromeric duplication and transcription by analyzing the in silico transcriptional profiles within the proximal 1.5 Mb of genomic sequence on all human chromosome arms in relation to duplication status. We identify an approximately 5x excess of transcripts specific to cancer and/or testis in pericentromeric duplications compared to surrounding single copy sequence, with the expression of >50% of all transcripts in duplications being restricted to these tissues. We also identify an approximately 5x excess of transcripts in duplications which contain large quantities of interspersed repeats. These results indicate that the transcriptional profiles of duplicated and single copy sequences within pericentromeric DNA are distinct, suggesting that pericentromeric instability is unlikely to represent a common route for gene creation but may have a disproportionate effect upon genes whose function is restricted to the germ line.

Pooled genomic indexing of rhesus macaque

Pooled genomic indexing (PGI) is a method for mapping collections of bacterial artificial chromosome (BAC) clones between species by using a combination of clone pooling and DNA sequencing. PGI has been used to map a total of 3858 BAC clones covering approximately 24% of the rhesus macaque (Macaca mulatta) genome onto 4178 homologous loci in the human genome. A number of intrachromosomal rearrangements were detected by mapping multiple segments within the individual rhesus BACs onto multiple disjoined loci in the human genome. Transversal pooling designs involving shuffled BAC arrays were employed for robust mapping even with modest DNA sequence read coverage. A further innovation, short-tag pooled genomic indexing (ST-PGI), was also introduced to further improve the economy of mapping by sequencing multiple, short, mapable tags within a single sequencing reaction.

Phylogenomic approaches to common problems encountered in the analysis of low copy repeats: the sulfotransferase 1A gene family example

Background

Blocks of duplicated genomic DNA sequence longer than 1000 base pairs are known as low copy repeats (LCRs). Identified by their sequence similarity, LCRs are abundant in the human genome, and are interesting because they may represent recent adaptive events, or potential future adaptive opportunities within the human lineage. Sequence analysis tools are needed, however, to decide whether these interpretations are likely, whether a particular set of LCRs represents nearly neutral drift creating junk DNA, or whether the appearance of LCRs reflects assembly error. Here we investigate an LCR family containing the sulfotransferase (SULT) 1A genes involved in drug metabolism, cancer, hormone regulation, and neurotransmitter biology as a first step for defining the problems that those tools must manage.

Results

Sequence analysis here identified a fourth sulfotransferase gene, which may be transcriptionally active, located on human chromosome 16. Four regions of genomic sequence containing the four human SULT1A paralogs defined a new LCR family. The stem hominoid SULT1A progenitor locus was identified by comparative genomics involving complete human and rodent genomes, and a draft chimpanzee genome. SULT1A expansion in hominoid genomes was followed by positive selection acting on specific protein sites. This episode of adaptive evolution appears to be responsible for the dopamine sulfonation function of some SULT enzymes. Each of the conclusions that this bioinformatic analysis generated using data that has uncertain reliability (such as that from the chimpanzee genome sequencing project) has been confirmed experimentally or by a "finished" chromosome 16 assembly, both of which were published after the submission of this manuscript.

Conclusion

SULT1A genes expanded from one to four copies in hominoids during intra-chromosomal LCR duplications, including (apparently) one after the divergence of chimpanzees and humans. Thus, LCRs may provide a means for amplifying genes (and other genetic elements) that are adaptively useful. Being located on and among LCRs, however, could make the human SULT1A genes susceptible to further duplications or deletions resulting in 'genomic diseases' for some individuals. Pharmacogenomic studies of SULT1Asingle nucleotide polymorphisms, therefore, should also consider examining SULT1A copy number variability when searching for genotype-phenotype associations. The latest duplication is, however, only a substantiated hypothesis; an alternative explanation, disfavored by the majority of evidence, is that the duplication is an artifact of incorrect genome assembly.

The sequence and analysis of duplication-rich human chromosome 16

Human chromosome 16 features one of the highest levels of segmentally duplicated sequence among the human autosomes. We report here the 78,884,754 base pairs of finished chromosome 16 sequence, representing over 99.9% of its euchromatin. Manual annotation revealed 880 protein-coding genes confirmed by 1,670 aligned transcripts, 19 transfer RNA genes, 341 pseudogenes and three RNA pseudogenes. These genes include metallothionein, cadherin and iroquois gene families, as well as the disease genes for polycystic kidney disease and acute myelomonocytic leukaemia. Several large-scale structural polymorphisms spanning hundreds of kilobase pairs were identified and result in gene content differences among humans. Whereas the segmental duplications of chromosome 16 are enriched in the relatively gene-poor pericentromere of the p arm, some are involved in recent gene duplication and conversion events that are likely to have had an impact on the evolution of primates and human disease susceptibility.

Fluorescence in situ hybridization to chromosomes as a tool to understand human and primate genome evolution

For the last 15 years molecular cytogenetic techniques have been extensively used to study primate evolution. Molecular probes were helpful to distinguish mammalian chromosomes and chromosome segments on the basis of their DNA content rather than solely on morphological features such as banding patterns. Various landmark rearrangements have been identified for most of the nodes in primate phylogeny while chromosome banding still provides helpful reference maps. Fluorescence in situ hybridization (FISH) techniques were used with probes of different complexity including chromosome painting probes, probes derived from chromosome sub-regions and in the size of a single gene. Since more recently, in silico techniques have been applied to trace down evolutionarily derived chromosome rearrangements by searching the human and mouse genome sequence databases. More detailed breakpoint analyses of chromosome rearrangements that occurred during higher primate evolution also gave some insights into the molecular changes in chromosome rearrangements that occurred in evolution. Hardly any "fusion genes" as known from chromosome rearrangements in cancer cells or dramatic "position effects" of genes transferred to new sites in primate genomes have been reported yet. Most breakpoint regions have been identified within gene poor areas rich in repetitive elements and/or low copy repeats (segmental duplications). The progress in various molecular and molecular-cytogenetic approaches including the recently launched chimpanzee genome project suggests that these new tools will have a significant impact on the further understanding of human genome evolution.

Comparative analysis of gene-expression patterns in human and African great ape cultured fibroblasts

Although much is known about genetic variation in human and African great ape (chimpanzee, bonobo, and gorilla) genomes, substantially less is known about variation in gene-expression profiles within and among these species. This information is necessary for defining transcriptional regulatory networks that contribute to complex phenotypes unique to humans or the African great apes. We took a systematic approach to this problem by investigating gene-expression profiles in well-defined cell populations from humans, bonobos, and gorillas. By comparing these profiles from 18 human and 21 African great ape primary fibroblast cell lines, we found that gene-expression patterns could predict the species, but not the age, of the fibroblast donor. Several differentially expressed genes among human and African great ape fibroblasts involved the extracellular matrix, metabolic pathways, signal transduction, stress responses, as well as inherited overgrowth and neurological disorders. These gene-expression patterns could represent molecular adaptations that influenced the development of species-specific traits in humans and the African great apes.

The Tre2 (USP6) oncogene is a hominoid-specific gene

Gene duplication and domain accretion are thought to be the major mechanisms for the emergence of novel genes during evolution. Such events are thought to have occurred at early stages in the vertebrate lineage, but genomic sequencing has recently revealed extensive amplification events during the evolution of higher primates. We report here that the Tre2 (USP6) oncogene is derived from the chimeric fusion of two genes, USP32 (NY-REN-60), and TBC1D3. USP32 is an ancient, highly conserved gene, whereas TBC1D3 is derived from a recent segmental duplication, which is absent in most other mammals and shows rapid amplification and dispersal through the primate lineage. Remarkably, the chimeric gene Tre2 exists only in the hominoid lineage of primates. This hominoid-specific oncogene arose as recently as 21-33 million years ago, after proliferation of the TBC1D3 segmental duplication in the primate lineage. In contrast to the broad expression pattern of USP32 and TBC1D3, expression of Tre2 is testis-specific, a pattern proposed for novel genes implicated in the emergence of reproductive barriers. The sudden emergence of chimeric proteins, such as that encoded by Tre2, may have contributed to hominoid speciation.

Refinement of a chimpanzee pericentric inversion breakpoint to a segmental duplication cluster

BACKGROUND

Pericentric inversions are the most common euchromatic chromosomal differences among humans and the great apes. The human and chimpanzee karyotype differs by nine such events, in addition to several constitutive heterochromatic increases and one chromosomal fusion event. Reproductive isolation and subsequent speciation are thought to be the potential result of pericentric inversions, as reproductive boundaries form as a result of hybrid sterility.

RESULTS

Here we employed a comparative fluorescence in situ hybridization approach, using probes selected from a combination of physical mapping, genomic sequence, and segmental duplication analyses to narrow the breakpoint interval of a pericentric inversion in chimpanzee involving the orthologous human 15q11-q13 region. We have refined the inversion breakpoint of this chimpanzee-specific rearrangement to a 600 kilobase (kb) interval of the human genome consisting of entirely duplicated material. Detailed analysis of the underlying sequence indicated that this region comprises multiple segmental duplications, including a previously characterized duplication of the alpha7 neuronal nicotinic acetylcholine receptor subunit gene (CHRNA7) in 15q13.3 and several Golgin-linked-to-PML, or LCR15, duplications.

CONCLUSIONS

We conclude that, on the basis of experimental data excluding the CHRNA7 duplicon as the site of inversion, and sequence analysis of regional duplications, the most likely rearrangement site is within a GLP/LCR15 duplicon. This study further exemplifies the genomic plasticity due to the presence of segmental duplications and highlights their importance for a complete understanding of genome evolution.