Serial segmental duplications during primate evolution result in complex human genome architecture

PhaseDancer: a novel targeted assembler of segmental duplications unravels the complexity of the human chromosome 2 fusion going from 48 to 46 chromosomes in hominin evolution

Resolving complex genomic regions rich in segmental duplications (SDs) is challenging due to the high error rate of long-read sequencing. Here, we describe a targeted approach with a novel genome assembler PhaseDancer that extends SD-rich regions of interest iteratively. We validate its robustness and efficiency using a golden-standard set of human BAC clones and in silico-generated SDs with predefined evolutionary scenarios. PhaseDancer enables extension of the incomplete complex SD-rich subtelomeric regions of Great Ape chromosomes orthologous to the human chromosome 2 (HSA2) fusion site, informing a model of HSA2 formation and unravelling the evolution of human and Great Ape genomes.

Genome-wide detection of CNV regions between Anqing six-end-white and Duroc pigs

BACKGROUND

Anqing six-end-white pig is a native breed in Anhui Province. The pigs have the disadvantages of a slow growth rate, low proportion of lean meat, and thick back fat, but feature the advantages of strong stress resistance and excellent meat quality. Duroc pig is an introduced pig breed with a fast growth rate and high proportion of lean meat. With the latter breed featuring superior growth characteristics but inferior meat quality traits, the underlying molecular mechanism that causes these phenotypic differences between Chinese and foreign pigs is still unclear.

RESULTS

In this study, copy number variation (CNV) detection was performed using the re-sequencing data of Anqing Six-end-white pigs and Duroc pigs, A total of 65,701 CNVs were obtained. After merging the CNVs with overlapping genomic positions, 881 CNV regions (CNVRs) were obtained. Based on the obtained CNVR information combined with their positions on the 18 chromosomes, a whole-genome map of the pig CNVs was drawn. GO analysis of the genes in the CNVRs showed that they were primarily involved in the cellular processes of proliferation, differentiation, and adhesion, and primarily involved in the biological processes of fat metabolism, reproductive traits, and immune processes.

CONCLUSION

The difference analysis of the CNVs between the Chinese and foreign pig breeds showed that the CNV of the Anqing six-end-white pig genome was higher than that of the introduced pig breed Duroc. Six genes related to fat metabolism, reproductive performance, and stress resistance were found in genome-wide CNVRs (DPF3, LEPR, MAP2K6, PPARA, TRAF6, NLRP4).

Formation of a Family of Long Intergenic Noncoding RNA Genes with an Embedded Translocation Breakpoint Motif in Human Chromosomal Low Copy Repeats of 22q11.2-Some Surprises and Questions

A family of long intergenic noncoding RNA (lincRNA) genes, FAM230 is formed via gene sequence duplication, specifically in human chromosomal low copy repeats (LCR) or segmental duplications. This is the first group of lincRNA genes known to be formed by segmental duplications and is consistent with current views of evolution and the creation of new genes via DNA low copy repeats. It appears to be an efficient way to form multiple lincRNA genes. But as these genes are in a critical chromosomal region with respect to the incidence of abnormal translocations and resulting genetic abnormalities, the 22q11.2 region, and also carry a translocation breakpoint motif, several intriguing questions arise concerning the presence and function of the translocation breakpoint sequence in RNA genes situated in LCR22s.

A local duplication of the Melanocortin receptor 1 locus in Astyanax

In this study, we report evidence of a novel duplication of Melanocortin receptor 1 (Mc1r) in the cavefish genome. This locus was discovered following the observation of excessive allelic diversity in a ∼820 bp fragment of Mc1r amplified via degenerate PCR from a natural population of Astyanax aeneus fish from Guerrero, Mexico. The cavefish genome reveals the presence of two closely related Mc1r open reading frames separated by a 1.46 kb intergenic region. One open reading frame corresponds to the previously reported Mc1r receptor, and the other open reading frame (duplicate copy) is 975 bp in length, encoding a receptor of 325 amino acids. Sequence similarity analyses position both copies in the syntenic region of the single Mc1r locus in 16 representative craniate genomes spanning bony fish (including Astyanax) to mammals, suggesting we discovered tandem duplicates of this important gene. The two Mc1r copies share ∼89% sequence similarity and, within Astyanax, are more similar to one another compared to other melanocortin family members. Future studies will inform the precise functional significance of the duplicated Mc1r locus and if this novel copy number variant may have adaptive significance for the Astyanax lineage.

Genetic Approaches to Understanding Psychiatric Disease

Human genetic studies have been the driving force in bringing to light the underlying biology of psychiatric conditions. As these studies fill in the gaps in our knowledge of the mechanisms at play, we will be better equipped to design therapies in rational and targeted ways, or repurpose existing therapies in previously unanticipated ways. This review is intended for those unfamiliar with psychiatric genetics as a field and provides a primer on different modes of genetic variation, the technologies currently used to probe them, and concepts that provide context for interpreting the gene-phenotype relationship. Like other subfields in human genetics, psychiatric genetics is moving from microarray technology to sequencing-based approaches as barriers of cost and expertise are removed, and the ramifications of this transition are discussed here. A summary is then given of recent genetic discoveries in a number of neuropsychiatric conditions, with particular emphasis on neurodevelopmental conditions. The general impact of genetics on drug development has been to underscore the extensive etiological heterogeneity in seemingly cohesive diagnostic categories. Consequently, the path forward is not in therapies hoping to reach large swaths of patients sharing a clinically defined diagnosis, but rather in targeting patients belonging to specific "biotypes" defined through a combination of objective, quantifiable data, including genotype.

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

BACKGROUND

Gene innovation by duplication is a fundamental evolutionary process but is difficult to study in humans due to the large size, high sequence identity, and mosaic nature of segmental duplication blocks. The human-specific gene hydrocephalus-inducing 2, HYDIN2, was generated by a 364 kbp duplication of 79 internal exons of the large ciliary gene HYDIN from chromosome 16q22.2 to chromosome 1q21.1. Because the HYDIN2 locus lacks the ancestral promoter and seven terminal exons of the progenitor gene, we sought to characterize transcription at this locus by coupling reverse transcription polymerase chain reaction and long-read sequencing.

RESULTS

5' RACE indicates a transcription start site for HYDIN2 outside of the duplication and we observe fusion transcripts spanning both the 5' and 3' breakpoints. We observe extensive splicing diversity leading to the formation of altered open reading frames (ORFs) that appear to be under relaxed selection. We show that HYDIN2 adopted a new promoter that drives an altered pattern of expression, with highest levels in neural tissues. We estimate that the HYDIN duplication occurred ~3.2 million years ago and find that it is nearly fixed (99.9%) for diploid copy number in contemporary humans. Examination of 73 chromosome 1q21 rearrangement patients reveals that HYDIN2 is deleted or duplicated in most cases.

CONCLUSIONS

Together, these data support a model of rapid gene innovation by fusion of incomplete segmental duplications, altered tissue expression, and potential subfunctionalization or neofunctionalization of HYDIN2 early in the evolution of the Homo lineage.

Structure, function, and genetics of lipoprotein (a)

Lipoprotein (a) [Lp(a)] has attracted the interest of researchers and physicians due to its intriguing properties, including an intragenic multiallelic copy number variation in the LPA gene and the strong association with coronary heart disease (CHD). This review summarizes present knowledge of the structure, function, and genetics of Lp(a) with emphasis on the molecular and population genetics of the Lp(a)/LPA trait, as well as aspects of genetic epidemiology. It highlights the role of genetics in establishing Lp(a) as a risk factor for CHD, but also discusses uncertainties, controversies, and lack of knowledge on several aspects of the genetic Lp(a) trait, not least its function.

Comparative Genomic Analyses of the Human NPHP1 Locus Reveal Complex Genomic Architecture and Its Regional Evolution in Primates

Many loci in the human genome harbor complex genomic structures that can result in susceptibility to genomic rearrangements leading to various genomic disorders. Nephronophthisis 1 (NPHP1, MIM# 256100) is an autosomal recessive disorder that can be caused by defects of NPHP1; the gene maps within the human 2q13 region where low copy repeats (LCRs) are abundant. Loss of function of NPHP1 is responsible for approximately 85% of the NPHP1 cases—about 80% of such individuals carry a large recurrent homozygous NPHP1 deletion that occurs via nonallelic homologous recombination (NAHR) between two flanking directly oriented ~45 kb LCRs. Published data revealed a non-pathogenic inversion polymorphism involving the NPHP1 gene flanked by two inverted ~358 kb LCRs. Using optical mapping and array-comparative genomic hybridization, we identified three potential novel structural variant (SV) haplotypes at the NPHP1 locus that may protect a haploid genome from the NPHP1 deletion. Inter-species comparative genomic analyses among primate genomes revealed massive genomic changes during evolution. The aggregated data suggest that dynamic genomic rearrangements occurred historically within the NPHP1 locus and generated SV haplotypes observed in the human population today, which may confer differential susceptibility to genomic instability and the NPHP1 deletion within a personal genome. Our study documents diverse SV haplotypes at a complex LCR-laden human genomic region. Comparative analyses provide a model for how this complex region arose during primate evolution, and studies among humans suggest that intra-species polymorphism may potentially modulate an individual’s susceptibility to acquiring disease-associated alleles.

Genomic instability due to the intrinsic sequence architecture of the genome, such as low copy repeats (LCRs), is a major contributor to de novo mutations that can occur in the process of human genome evolution. LCRs can mediate genomic rearrangements associated with genomic disorders by acting as substrates for nonallelic homologous recombination. Juvenile-onset nephronophthisis 1 is the most frequent genetic cause of renal failure in children. An LCR-mediated, homozygous common recurrent deletion encompassing NPHP1 is found in the majority of affected subjects, while heterozygous deletion representing the nephronophthisis 1 recessive carrier state is frequently observed amongst world populations. Interestingly, the human NPHP1 locus is located proximal to the head-to-head fusion site of two ancestral chromosomes that occurred in the great apes, which resulted in a reduction of chromosome number from 48 in nonhuman primates to the current 46 in humans. In this study, we characterized and provided evidence for the diverse genomic architecture at the NPHP1 locus and potential structural variant haplotypes in the human population. Furthermore, our analyses of primate genomes shed light on the massive changes of genomic architecture at the human NPHP1 locus and delineated a model for the emergence of the LCRs during primate evolution.

Gene structure variation in segmental duplication block C of human chromosome 7q 11.23 during primate evolution

Segmental duplication, or low-copy repeat (LCR) event, occurs during primate evolution and is an important source of genomic diversity, including gain or loss of gene function. The human chromosome 7q 11.23 is related to the William-Beuren syndrome and contains large region-specific LCRs composed of blocks A, B, and C that have different copy numbers in humans and different primates. We analyzed the structure of POM121, NSUN5, FKBP6, and TRIM50 genes in the LCRs of block C. Based on computational analysis, POM121B created by a segmental duplication acquired a new exonic region, whereas NSUN5B (NSUN5C) showed structural variation by integration of HERV-K LTR after duplication from the original NSUN5 gene. The TRIM50 gene originally consists of seven exons, whereas the duplicated TRIM73 and TRIM74 genes present five exons because of homologous recombination-mediated deletion. In addition, independent duplication events of the FKBP6 gene generated two pseudogenes at different genomic locations. In summary, these clustered genes are created by segmental duplication, indicating that they show dynamic evolutionary events, leading to structure variation in the primate genome.

Nonrecurrent 17p11.2p12 Rearrangement Events that Result in Two Concomitant Genomic Disorders: The PMP22-RAI1 Contiguous Gene Duplication Syndrome

The genomic duplication associated with Potocki-Lupski syndrome (PTLS) maps in close proximity to the duplication associated with Charcot-Marie-Tooth disease type 1A (CMT1A). PTLS is characterized by hypotonia, failure to thrive, reduced body weight, intellectual disability, and autistic features. CMT1A is a common autosomal dominant distal symmetric peripheral polyneuropathy. The key dosage-sensitive genes RAI1 and PMP22 are respectively associated with PTLS and CMT1A. Recurrent duplications accounting for the majority of subjects with these conditions are mediated by nonallelic homologous recombination between distinct low-copy repeat (LCR) substrates. The LCRs flanking a contiguous genomic interval encompassing both RAI1 and PMP22 do not share extensive homology; thus, duplications encompassing both loci are rare and potentially generated by a different mutational mechanism. We characterized genomic rearrangements that simultaneously duplicate PMP22 and RAI1, including nine potential complex genomic rearrangements, in 23 subjects by high-resolution array comparative genomic hybridization and breakpoint junction sequencing. Insertions and microhomologies were found at the breakpoint junctions, suggesting potential replicative mechanisms for rearrangement formation. At the breakpoint junctions of these nonrecurrent rearrangements, enrichment of repetitive DNA sequences was observed, indicating that they might predispose to genomic instability and rearrangement. Clinical evaluation revealed blended PTLS and CMT1A phenotypes with a potential earlier onset of neuropathy. Moreover, additional clinical findings might be observed due to the extra duplicated material included in the rearrangements. Our genomic analysis suggests replicative mechanisms as a predominant mechanism underlying PMP22-RAI1 contiguous gene duplications and provides further evidence supporting the role of complex genomic architecture in genomic instability.

Copy Number Variation at the APOL1 Locus

Two coding variants in the APOL1 gene (G1 and G2) explain most of the high rate of kidney disease in African Americans. APOL1-associated kidney disease risk inheritance follows an autosomal recessive pattern: The relative risk of kidney disease associated with inheritance of two high-risk variants is 7–30 fold, depending on the specific kidney phenotype. We wished to determine if the variability in phenotype might in part reflect structural differences in APOL1 gene. We analyzed sequence coverage from 1000 Genomes Project Phase 3 samples as well as exome sequencing data from African American kidney disease cases for copy number variation. 8 samples sequenced in the 1000 Genomes Project showed increased coverage over a ~100kb region that includes APOL2, APOL1 and part of MYH9, suggesting the presence of APOL1 copy number greater than 2. We reasoned that such duplications should be enriched in apparent G1 heterozygotes with kidney disease. Using a PCR-based assay, we observed the presence of this duplication in additional samples from apparent G0G1 or G0G2 individuals. The frequency of this APOL1 duplication was compared among cases (n = 123) and controls (n = 255) with apparent G0G1 heterozygosity. The presence of APOL1 duplication was observed in 4.06% of cases and 0.78% controls, preliminary evidence that this APOL1 duplication may alter susceptibility to kidney disease (p = 0.03). Taqman-based copy number assays confirmed the presence of 3 APOL1 copies in individuals positive for this specific duplication by PCR assay, but also identified a small number of individuals with additional APOL1 copies of presumably different structure. These observations motivate further studies to better assess the contribution of APOL1 copy number on kidney disease risk and on APOL1 function. Investigators and clinicians genotyping APOL1 should also consider whether the particular genotyping platform used is subject to technical errors when more than two copies of APOL1 are present.

5'RUNX1-3'USP42 chimeric gene in acute myeloid leukemia can occur through an insertion mechanism rather than translocation and may be mediated by genomic segmental duplications

Background

The runt-related transcription factor 1 (RUNX1) gene is a transcription factor that acts as a master regulator of hematopoiesis and represents one of the most frequent targets of chromosomal rearrangements in human leukemias. The t(7;21)(p22;q22) rearrangement generating a 5′RUNX1-3′USP42 fusion transcript has been reported in two cases of pediatric acute myeloid leukemia (AML) and further in eight adult cases of myeloid neoplasms. We describe the first case of adult AML with a 5′RUNX1-3′USP42 fusion gene generated by an insertion event instead of chromosomal translocation.

Methods

Conventional and molecular cytogenetic analyses allowed the precise characterization of the chromosomal rearrangement and breakpoints identification. Gene expression analysis was performed by quantitative real-time PCR experiments, whereas bioinformatic studies were carried out for revealing structural genomic characteristics of breakpoint regions.

Results

We identified an adult AML case bearing a ins(21;7)(q22;p15p22) generating a 5′RUNX1-3′USP42 fusion gene on der(21) chromosome and causing USP42 gene over-expression. Bioinformatic analysis of the genomic regions involved in ins(21;7)/t(7;21) showed the presence of interchromosomal segmental duplications (SDs) next to the USP42 and RUNX1 genes, that may underlie a non-allelic homologous recombination between chromosome 7 and 21 in AML.

Conclusions

We report the first case of a 5′RUNX1-3′USP42 chimeric gene generated by a chromosomal cryptic insertion in an adult AML patient. Our data revealed that there may be a pivotal role for SDs in this very rare but recurrent chromosomal rearrangement.

DNA clustering and genome complexity

Early global measures of genome complexity (power spectra, the analysis of fluctuations in DNA walks or compositional segmentation) uncovered a high degree of complexity in eukaryotic genome sequences. The main evolutionary mechanisms leading to increases in genome complexity (i.e. gene duplication and transposon proliferation) can all potentially produce increases in DNA clustering. To quantify such clustering and provide a genome-wide description of the formed clusters, we developed GenomeCluster, an algorithm able to detect clusters of whatever genome element identified by chromosome coordinates. We obtained a detailed description of clusters for ten categories of human genome elements, including functional (genes, exons, introns), regulatory (CpG islands, TFBSs, enhancers), variant (SNPs) and repeat (Alus, LINE1) elements, as well as DNase hypersensitivity sites. For each category, we located their clusters in the human genome, then quantifying cluster length and composition, and estimated the clustering level as the proportion of clustered genome elements. In average, we found a 27% of elements in clusters, although a considerable variation occurs among different categories. Genes form the lowest number of clusters, but these are the longest ones, both in bp and the average number of components, while the shortest clusters are formed by SNPs. Functional and regulatory elements (genes, CpG islands, TFBSs, enhancers) show the highest clustering level, as compared to DNase sites, repeats (Alus, LINE1) or SNPs. Many of the genome elements we analyzed are known to be composed of clusters of low-level entities. In addition, we found here that the clusters generated by GenomeCluster can be in turn clustered into high-level super-clusters. The observation of 'clusters-within-clusters' parallels the 'domains within domains' phenomenon previously detected through global statistical methods in eukaryotic sequences, and reveals a complex human genome landscape dominated by hierarchical clustering.

Characterizing regions in the human genome unmappable by next-generation-sequencing at the read length of 1000 bases

Repetitive and redundant regions of a genome are particularly problematic for mapping sequencing reads. In the present paper, we compile a list of the unmappable regions in the human genome based on the following definition: hypothetical reads with length 1 kb which cannot be uniquely mapped with zero-mismatch alignment for the described regions, considering both the forward and reverse strand. The respective collection of unmappable regions covers 0.77% of the sequence of human autosomes and 8.25% of the sex chromosomes in the reference genome GRCh37/hg19 (overall 1.23%). Not surprisingly, our unmappable regions overlap greatly with segmental duplication, transposable elements, and structural variants. About 99.8% of bases in our unmappable regions are part of either segmental duplication or transposable elements and 98.3% overlap structural variant annotations. Notably, some of these regions overlap units with important biological functions, including 4% of protein-coding genes. In contrast, these regions have zero intersection with the ultraconserved elements, very low overlap with microRNAs, tRNAs, pseudogenes, CpG islands, tandem repeats, microsatellites, sensitive non-coding regions, and the mapping blacklist regions from the ENCODE project.

Distribution of segmental duplications in the context of higher order chromatin organisation of human chromosome 7

BACKGROUND

Segmental duplications (SDs) are not evenly distributed along chromosomes. The reasons for this biased susceptibility to SD insertion are poorly understood. Accumulation of SDs is associated with increased genomic instability, which can lead to structural variants and genomic disorders such as the Williams-Beuren syndrome. Despite these adverse effects, SDs have become fixed in the human genome. Focusing on chromosome 7, which is particularly rich in interstitial SDs, we have investigated the distribution of SDs in the context of evolution and the three dimensional organisation of the chromosome in order to gain insights into the mutual relationship of SDs and chromatin topology.

RESULTS

Intrachromosomal SDs preferentially accumulate in those segments of chromosome 7 that are homologous to marmoset chromosome 2. Although this formerly compact segment has been re-distributed to three different sites during primate evolution, we can show by means of public data on long distance chromatin interactions that these three intervals, and consequently the paralogous SDs mapping to them, have retained their spatial proximity in the nucleus. Focusing on SD clusters implicated in the aetiology of the Williams-Beuren syndrome locus we demonstrate by cross-species comparison that these SDs have inserted at the borders of a topological domain and that they flank regions with distinct DNA conformation.

CONCLUSIONS

Our study suggests a link of nuclear architecture and the propagation of SDs across chromosome 7, either by promoting regional SD insertion or by contributing to the establishment of higher order chromatin organisation themselves. The latter could compensate for the high risk of structural rearrangements and thus may have contributed to their evolutionary fixation in the human genome.

Maternal complex chromosomal rearrangement leads to TCF12 microdeletion in a patient presenting with coronal craniosynostosis and intellectual disability

We report on a young child with intellectual disability and unilateral coronal craniosynostosis leading to craniofacial malformations. Standard karyotype showed an apparently balanced translocation between chromosomes 2 and 15 [t(2;15)(q21;q21.3)], inherited from his mother. Interestingly, array-CGH 180K showed a 3.64 Mb de novo deletion on chromosome 15 in the region 15q21.3q22.2, close to the chromosome 15 translocation breakpoints. This deletion leads to haploinsufficiency of TCF12 gene that can explain the coronal craniosynostosis described in the patient. Additional FISH analyses showed a complex balanced maternal chromosomal rearrangement combining the reciprocal translocation t(2;15)(q21;q21.3), and an insertion of the 15q22.1 segment into the telomeric region of the translocated 15q fragment. The genomic imbalance in the patient is likely caused by a crossing-over that occurs in the recombination loop formed during the maternal meiosis resulting in the deletion of the inserted fragment. This original case of a genomic microdeletion of TCF12 exemplifies the importance of array-CGH in the clinical investigation of apparently balanced rearrangements but also the importance of FISH analysis to identify the chromosomal mechanism causing the genomic imbalance.

Fine-scale signatures of molecular evolution reconcile models of indel-associated mutation

Genomic structural alterations that vary within species, known as large copy number variants, represent an unanticipated and abundant source of genetic diversity that associates with variation in gene expression and susceptibility to disease. Even short insertions and deletions (indels) can exert important effects on genomes by locally increasing the mutation rate, with multiple mechanisms proposed to account for this pattern. To better understand how indels promote genome evolution, we demonstrate that the single nucleotide mutation rate is elevated in the vicinity of indels, with a resolution of tens of base pairs, for the two closely related nematode species Caenorhabditis remanei and C. sp. 23. In addition to indels being clustered with single nucleotide polymorphisms and fixed differences, we also show that transversion mutations are enriched in sequences that flank indels and that many indels associate with sequence repeats. These observations are compatible with a model that reconciles previously proposed mechanisms of indel-associated mutagenesis, implicating repeat sequences as a common driver of indel errors, which then recruit error-prone polymerases during DNA repair, resulting in a locally elevated single nucleotide mutation rate. The striking influence of indel variants on the molecular evolution of flanking sequences strengthens the emerging general view that mutations can induce further mutations.

Association of the BRCA1 promoter polymorphism rs11655505 with the risk of familial breast and/or ovarian cancer

Germline mutations in the BRCA1 tumor suppressor gene predispose affected individuals to breast cancer; however, incomplete cancer penetrance and the presence of phenocopies in BRCA1 families also indicate genetic and environmental modifiers of breast cancer risk. In this study, we have tested the single nucleotide polymorphism rs1655505 of the BRCA1 promoter, as candidate for the modifier of breast cancer risk. The polymorphic variants were genotyped in BRCA1-negative (729), familial breast and/or ovarian cancer cases (FBOC), including cases with a reported maternal history (154), nonfamilal (sporadic) cases (600), hereditary breast/ovarian cases with BRCA1 mutations (190) and population controls (1,590) from Central Poland. An association with the risk of FBOC was observed for the minor (T) allele and (TT) genotype (T: p = 0.006, OR = 1.40, 95% CI = 1.10-1.79; TT: p = 0.001, OR = 2.23, 95% CI = 1.37-3.62) in female cases with a reported maternal history, specifically in women with the onset of disease after 50 years of age (T: p = 0.004, OR = 1.77, 95% CI = 1.20-2.62; TT: p = 0.001, OR = 3.7, 95% CI = 1.62-8.46). The presented evidence suggests a need to conduct larger studies on the association between genetic variations at the BRCA1 promoter and the breast cancer risk, according to maternal/paternal lineage.

Genomic hypomethylation in the human germline associates with selective structural mutability in the human genome

The hotspots of structural polymorphisms and structural mutability in the human genome remain to be explained mechanistically. We examine associations of structural mutability with germline DNA methylation and with non-allelic homologous recombination (NAHR) mediated by low-copy repeats (LCRs). Combined evidence from four human sperm methylome maps, human genome evolution, structural polymorphisms in the human population, and previous genomic and disease studies consistently points to a strong association of germline hypomethylation and genomic instability. Specifically, methylation deserts, the ∼1% fraction of the human genome with the lowest methylation in the germline, show a tenfold enrichment for structural rearrangements that occurred in the human genome since the branching of chimpanzee and are highly enriched for fast-evolving loci that regulate tissue-specific gene expression. Analysis of copy number variants (CNVs) from 400 human samples identified using a custom-designed array comparative genomic hybridization (aCGH) chip, combined with publicly available structural variation data, indicates that association of structural mutability with germline hypomethylation is comparable in magnitude to the association of structural mutability with LCR–mediated NAHR. Moreover, rare CNVs occurring in the genomes of individuals diagnosed with schizophrenia, bipolar disorder, and developmental delay and de novo CNVs occurring in those diagnosed with autism are significantly more concentrated within hypomethylated regions. These findings suggest a new connection between the epigenome, selective mutability, evolution, and human disease.

The human genome contains many loci with high incidence of structural mutations, including insertions and deletions of chromosomal segments. This excessive mutability has accelerated evolution and contributed to human disease but has yet to be explained. Segments of DNA repeated in low-copy numbers (LCRs) have been previously implicated in promoting structural mutability in specific disease-associated loci. Lack of methylation (hypomethylation) of genomic DNA has been previously associated with high structural mutability in gibbons and in human cancer cells, but the association with structural mutability in the human germline has not been explored prior to this study. Our analyses confirm the role of LCRs in promoting structural mutability on the genome scale but also reveal a surprisingly strong association of genomic instability with hypomethylation. Specifically, evolutionary analyses reveal that methylation deserts, the ∼1% fraction of the human genome with the lowest methylation in human sperm, harbor a tenfold higher number of structural mutations than genome-wide average. Moreover, the structural mutations in individuals diagnosed with schizophrenia, bipolar disorder, developmental delay, and autism are significantly more concentrated within hypomethylated regions. Our findings suggest a new connection between methylation of genomic DNA, selective structural mutability, evolution, and human disease.

What have studies of genomic disorders taught us about our genome?

The elucidation of genomic disorders began with molecular technologies that enabled detection of genomic changes which were (a) smaller than those resolved by traditional cytogenetics (less than 5 Mb) and (b) larger than what could be determined by conventional gel electrophoresis. Methods such as pulsed field gel electrophoresis (PFGE) and fluorescent in situ hybridization (FISH) could resolve such changes but were limited to locus-specific studies. The study of genomic disorders has rapidly advanced with the development of array-based techniques. These enabled examination of the entire human genome at a higher level of resolution, thus allowing elucidation of the basis of many new disorders, mechanisms that result in genomic changes that can result in copy number variation (CNV), and most importantly, a deeper understanding of the characteristics, features, and plasticity of our genome. In this chapter, we focus on the structural and architectural features of the genome, which can potentially result in genomic instability, delineate how mechanisms, such as NAHR, NHEJ, and FoSTeS/MMBIR lead to disease-causing rearrangements, and briefly describe the relationship between the leading methods presently used in studying genomic disorders. We end with a discussion on our new understanding about our genome including: the contribution of new mutation CNV to disease, the abundance of mosaicism, the extent of subtelomeric rearrangements, the frequency of de novo rearrangements associated with sporadic birth defects, the occurrence of balanced and unbalanced translocations, the increasing discovery of insertional translocations, the exploration of complex rearrangements and exonic CNVs. In the postgenomic era, our understanding of the genome has advanced very rapidly as the level of technical resolution has become higher. This leads to a greater understanding of the effects of rearrangements present both in healthy subjects and individuals with clinically relevant phenotypes.

Frequency of nonallelic homologous recombination is correlated with length of homology: evidence that ectopic synapsis precedes ectopic crossing-over

Genomic disorders constitute a class of diseases that are associated with DNA rearrangements resulting from region-specific genome instability, that is, genome architecture incites genome instability. Nonallelic homologous recombination (NAHR) or crossing-over in meiosis between sequences that are not in allelic positions (i.e., paralogous sequences) can result in recurrent deletions or duplications causing genomic disorders. Previous studies of NAHR have focused on description of the phenomenon, but it remains unclear how NAHR occurs during meiosis and what factors determine its frequency. Here we assembled two patient cohorts with reciprocal genomic disorders; deletion associated Smith-Magenis syndrome and duplication associated Potocki-Lupski syndrome. By assessing the full spectrum of rearrangement types from the two cohorts, we find that complex rearrangements (those with more than one breakpoint) are more prevalent in copy-number gains (17.7%) than in copy-number losses (2.3%); an observation that supports a role for replicative mechanisms in complex rearrangement formation. Interestingly, for NAHR-mediated recurrent rearrangements, we show that crossover frequency is positively associated with the flanking low-copy repeat (LCR) length and inversely influenced by the inter-LCR distance. To explain this, we propose that the probability of ectopic chromosome synapsis increases with increased LCR length, and that ectopic synapsis is a necessary precursor to ectopic crossing-over.

Copy number variations at the Prader-Willi syndrome region on chromosome 15 and associations with obesity in whites

Obesity is a serious health problem with strong genetic determination. Copy number variation (CNV) is a common type of genomic variant associated with some complex human diseases. However, it is not clear how CNVs contribute to the etiology of obesity. In this study, we examined 1,000 unrelated US whites to search for CNVs that may predispose to obesity. We focused our analyses on the Prader-Willi syndrome (PWS) critical region (chromosome 15q11-q13), because the PWS region is a hotspot for CNV generation and obesity is one of the major clinical manifestations for chromosome abnormalities at this region. We constructed a map containing 39 CNVs at the PWS critical region with CNV occurrence rates higher than 1%. Among them, three CNVs were significantly associated with body fat mass (P < 0.05), with a higher copy number (CN) associated with an increase of 5.08-9.77 kg in body fat mass. These three CNVs are close to two known PWS genes, NDN (necdin homolog) and C15orf2 (chromosome 15 open reading frame 2), and partially overlap with another obesity gene PWRN1 (Prader-Willi region nonprotein-coding RNA 1). Interestingly, our recently published whole genome association scan study using the same sample by examining single-nucleotide polymorphisms (SNPs) did not find any significant associations at these CNV regions, suggesting the importance of examining both CNVs and SNPs for better understanding of genetic basis of obesity. Further studies are warranted to validate these CNVs and their importance to obesity.

Rad51 inhibits translocation formation by non-conservative homologous recombination in Saccharomyces cerevisiae

Chromosomal translocations are a primary biological response to ionizing radiation (IR) exposure, and are likely to result from the inappropriate repair of the DNA double-strand breaks (DSBs) that are created. An abundance of repetitive sequences in eukaryotic genomes provides ample opportunity for such breaks to be repaired by homologous recombination (HR) between non-allelic repeats. Interestingly, in the budding yeast, Saccharomyces cerevisiae the central strand exchange protein, Rad51 that is required for DSB repair by gene conversion between unlinked repeats that conserves genomic structure also suppresses translocation formation by several HR mechanisms. In particular, Rad51 suppresses translocation formation by single-strand annealing (SSA), perhaps the most efficient mechanism for translocation formation by HR in both yeast and mammalian cells. Further, the enhanced translocation formation that emerges in the absence of Rad51 displays a distinct pattern of genetic control, suggesting that this occurs by a separate mechanism. Since hypomorphic mutations in RAD51 in mammalian cells also reduce DSB repair by conservative gene conversion and stimulate non-conservative repair by SSA, this mechanism may also operate in humans and, perhaps contribute to the genome instability that propels the development of cancer.

Mechanisms for nonrecurrent genomic rearrangements associated with CMT1A or HNPP: rare CNVs as a cause for missing heritability

Genomic rearrangements involving the peripheral myelin protein gene (PMP22) in human chromosome 17p12 are associated with neuropathy: duplications cause Charcot-Marie-Tooth disease type 1A (CMT1A), whereas deletions lead to hereditary neuropathy with liability to pressure palsies (HNPP). Our previous studies showed that >99% of these rearrangements are recurrent and mediated by nonallelic homologous recombination (NAHR). Rare copy number variations (CNVs) generated by nonrecurrent rearrangements also exist in 17p12, but their underlying mechanisms are not well understood. We investigated 21 subjects with rare CNVs associated with CMT1A or HNPP by oligonucleotide-based comparative genomic hybridization microarrays and breakpoint sequence analyses, and we identified 17 unique CNVs, including two genomic deletions, ten genomic duplications, two complex rearrangements, and three small exonic deletions. Each of these CNVs includes either the entire PMP22 gene, or exon(s) only, or ultraconserved potential regulatory sequences upstream of PMP22, further supporting the contention that PMP22 is the critical gene mediating the neuropathy phenotypes associated with 17p12 rearrangements. Breakpoint sequence analysis reveals that, different from the predominant NAHR mechanism in recurrent rearrangement, various molecular mechanisms, including nonhomologous end joining, Alu-Alu-mediated recombination, and replication-based mechanisms (e.g., FoSTeS and/or MMBIR), can generate nonrecurrent 17p12 rearrangements associated with neuropathy. We document a multitude of ways in which gene function can be altered by CNVs. Given the characteristics, including small size, structural complexity, and location outside of coding regions, of selected rare CNVs, their identification remains a challenge for genome analysis. Rare CNVs may potentially represent an important portion of "missing heritability" for human diseases.

Identification of uncommon recurrent Potocki-Lupski syndrome-associated duplications and the distribution of rearrangement types and mechanisms in PTLS

Nonallelic homologous recombination (NAHR) can mediate recurrent rearrangements in the human genome and cause genomic disorders. Smith-Magenis syndrome (SMS) and Potocki-Lupski syndrome (PTLS) are genomic disorders associated with a 3.7 Mb deletion and its reciprocal duplication in 17p11.2, respectively. In addition to these common recurrent rearrangements, an uncommon recurrent 5 Mb SMS-associated deletion has been identified. However, its reciprocal duplication predicted by the NAHR mechanism had not been identified. Here we report the molecular assays on 74 subjects with PTLS-associated duplications, 35 of whom are newly investigated. By both oligonucleotide-based comparative genomic hybridization and recombination hot spot analyses, we identified two cases of the predicted 5 Mb uncommon recurrent PTLS-associated duplication. Interestingly, the crossovers occur in proximity to a recently delineated allelic homologous recombination (AHR) hot spot-associated sequence motif, further documenting the common hot spot features shared between NAHR and AHR. An additional eight subjects with nonrecurrent PTLS duplications were identified. The smallest region of overlap (SRO) for all of the 74 PTLS duplications examined is narrowed to a 125 kb interval containing only RAI1, a gene recently further implicated in autism. Sequence complexities consistent with DNA replication-based mechanisms were identified in four of eight (50%) newly identified nonrecurrent PTLS duplications. Our findings of the uncommon recurrent PTLS-associated duplication at a relative prevalence reflecting the de novo mutation rate and the distribution of 17p11.2 duplication types in PTLS reveal insights into both the contributions of new mutations and the different underlying mechanisms that generate genomic rearrangements causing genomic disorders.

Common chromosomal fragile sites (CFS) may be involved in normal and traumatic cognitive stress memory consolidation and altered nervous system immunity

Previous reports of specific patterns of increased fragility at common chromosomal fragile sites (CFS) found in association with certain neurobehavioural disorders did not attract attention at the time due to a shift towards molecular approaches to delineate neuropsychiatric disorder candidate genes. Links with miRNA, altered methylation and the origin of copy number variation indicate that CFS region characteristics may be part of chromatinomic mechanisms that are increasingly linked with neuroplasticity and memory. Current reports of large-scale double-stranded DNA breaks in differentiating neurons and evidence of ongoing DNA demethylation of specific gene promoters in adult hippocampus may shed new light on the dynamic epigenetic changes that are increasingly appreciated as contributing to long-term memory consolidation. The expression of immune recombination activating genes in key stress-induced memory regions suggests the adoption by the brain of this ancient pattern recognition and memory system to establish a structural basis for long-term memory through controlled chromosomal breakage at highly specific genomic regions. It is furthermore considered that these mechanisms for management of epigenetic information related to stress memory could be linked, in some instances, with the transfer of the somatically acquired information to the germline. Here, rearranged sequences can be subjected to further selection and possible eventual retrotranscription to become part of the more stable coding machinery if proven to be crucial for survival and reproduction. While linkage of cognitive memory with stress and fear circuitry and memory establishment through structural DNA modification is proposed as a normal process, inappropriate activation of immune-like genomic rearrangement processes through traumatic stress memory may have the potential to lead to undesirable activation of neuro-inflammatory processes. These theories could have a significant impact on the interpretation of risks posed by heredity and the environment and the search for neuropsychiatric candidate genes.

Evolution in health and medicine Sackler colloquium: Genomic disorders: a window into human gene and genome evolution

Gene duplications alter the genetic constitution of organisms and can be a driving force of molecular evolution in humans and the great apes. In this context, the study of genomic disorders has uncovered the essential role played by the genomic architecture, especially low copy repeats (LCRs) or segmental duplications (SDs). In fact, regardless of the mechanism, LCRs can mediate or stimulate rearrangements, inciting genomic instability and generating dynamic and unstable regions prone to rapid molecular evolution. In humans, copy-number variation (CNV) has been implicated in common traits such as neuropathy, hypertension, color blindness, infertility, and behavioral traits including autism and schizophrenia, as well as disease susceptibility to HIV, lupus nephritis, and psoriasis among many other clinical phenotypes. The same mechanisms implicated in the origin of genomic disorders may also play a role in the emergence of segmental duplications and the evolution of new genes by means of genomic and gene duplication and triplication, exon shuffling, exon accretion, and fusion/fission events.

Sequencing primate genomes: what have we learned?

We summarize the progress in whole-genome sequencing and analyses of primate genomes. These emerging genome datasets have broadened our understanding of primate genome evolution revealing unexpected and complex patterns of evolutionary change. This includes the characterization of genome structural variation, episodic changes in the repeat landscape, differences in gene expression, new models regarding speciation, and the ephemeral nature of the recombination landscape. The functional characterization of genomic differences important in primate speciation and adaptation remains a significant challenge. Limited access to biological materials, the lack of detailed phenotypic data and the endangered status of many critical primate species have significantly attenuated research into the genetic basis of primate evolution. Next-generation sequencing technologies promise to greatly expand the number of available primate genome sequences; however, such draft genome sequences will likely miss critical genetic differences within complex genomic regions unless dedicated efforts are put forward to understand the full spectrum of genetic variation.

Footprints of X-to-Y gene conversion in recent human evolution

Different X-homologous regions of the male-specific portion of the human Y chromosome (MSY) are characterized by a different content of putative single nucleotide polymorphisms (SNPs), as reported in public databases. The possible role of X-to-Y nonallelic gene conversion in contributing to these differences remains poorly understood. We explored this issue by analyzing sequence variation in three regions of the MSY characterized by a different degree of X-Y similarity and a different density of putative SNPs: the PCDH11Y gene in the X-transposed (X-Y identity 99%, high putative SNP content); the TBL1Y gene in the X-degenerate (X-Y identity 86-88%, low putative SNP content); and VCY genes-containing region in the P8 palindrome (X-Y identity 95%, low putative SNP content). Present findings do not provide any evidence for gene conversion in the PCDH11Y and TBL1Y genes; they also strongly suggest that most putative SNPs of the PCDH11Y gene (and possibly the entire X-transposed region) are most likely X-Y paralogous sequence variants, which have been entered in the databases as SNPs. On the other hand, clear evidence for the VCY genes in the P8 palindrome having acted as an acceptor of X-to-Y gene conversion was obtained. A rate of 1.8 x 10(-7) X-to-Y conversions/bp/year was estimated for these genes. These findings indicate that in the VCY region of the MSY, X-to-Y gene conversion can be highly effective to increase the level of diversity among human Y chromosomes and suggest an additional explanation for the ability of the Y chromosome to retard degradation during evolution. Present data are expected to pave the way for future investigations on the role of nonallelic gene conversion in double-strand break repair and the maintenance of Y chromosome integrity.

The evolution of human segmental duplications and the core duplicon hypothesis

Duplicated sequences are important sources of genetic instability and in the evolution of new gene function within species. Hominids have a preponderance of intrachromosomal duplications organized in an interspersed fashion, as opposed to tandem duplications, which are common in other mammalian genomes such as mouse, dog, and cow. Multiple lines of evidence, including sequence divergence, comparative primate genomes, and fluorescence in situ hybridization (FISH) analyses, point to an excess of segmental duplications in the common ancestor of humans and African great apes. We find that much of the interspersed human duplication architecture within chromosomes is focused around common sequence elements referred to as "core duplicons." These cores correspond to the expansion of gene families, some of which show signatures of positive selection and lack orthologs present in other mammalian species. This genomic architecture predisposes apes and humans not only to extensive genetic diversity, but also to large-scale structural diversity mediated by nonallelic homologous recombination. In humans, many de novo large-scale genomic changes mediated by these duplications are associated with neuropsychiatric and neurodevelopmental disease. We propose that the disadvantage of a high rate of new mutations is offset by the selective advantage of newly minted genes within the cores.

Mechanisms of change in gene copy number

Deletions and duplications of chromosomal segments (copy number variants, CNVs) are a major source of variation between individual humans and are an underlying factor in human evolution and in many diseases, including mental illness, developmental disorders and cancer. CNVs form at a faster rate than other types of mutation, and seem to do so by similar mechanisms in bacteria, yeast and humans. Here we review current models of the mechanisms that cause copy number variation. Non-homologous end-joining mechanisms are well known, but recent models focus on perturbation of DNA replication and replication of non-contiguous DNA segments. For example, cellular stress might induce repair of broken replication forks to switch from high-fidelity homologous recombination to non-homologous repair, thus promoting copy number change.

Mechanisms of change in gene copy number

Deletions and duplications of chromosomal segments (copy number variants, CNVs) are a major source of variation between individual humans and are an underlying factor in human evolution and in many diseases, including mental illness, developmental disorders and cancer. CNVs form at a faster rate than other types of mutation, and seem to do so by similar mechanisms in bacteria, yeast and humans. Here we review current models of the mechanisms that cause copy number variation. Non-homologous end-joining mechanisms are well known, but recent models focus on perturbation of DNA replication and replication of non-contiguous DNA segments. For example, cellular stress might induce repair of broken replication forks to switch from high-fidelity homologous recombination to non-homologous repair, thus promoting copy number change.

Genomic disorders ten years on

It is now becoming generally accepted that a significant amount of human genetic variation is due to structural changes of the genome rather than to base-pair changes in the DNA. As for base-pair changes, knowledge of gene and genome function has been informed by structural alterations that convey clinical phenotypes. Genomic disorders are a class of human conditions that result from structural changes of the human genome that convey traits or susceptibility to traits. The path to the delineation of genomic disorders is intertwined with the evolving technologies that have enabled the resolution of human genome analyses to continue increasing. Similarly, the ability to perform high-resolution human genome analysis has fueled the current and future clinical implementation of such discoveries in the evolving field of genome medicine.

A burst of segmental duplications in the genome of the African great ape ancestor

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

A microhomology-mediated break-induced replication model for the origin of human copy number variation

Chromosome structural changes with nonrecurrent endpoints associated with genomic disorders offer windows into the mechanism of origin of copy number variation (CNV). A recent report of nonrecurrent duplications associated with Pelizaeus-Merzbacher disease identified three distinctive characteristics. First, the majority of events can be seen to be complex, showing discontinuous duplications mixed with deletions, inverted duplications, and triplications. Second, junctions at endpoints show microhomology of 2–5 base pairs (bp). Third, endpoints occur near pre-existing low copy repeats (LCRs). Using these observations and evidence from DNA repair in other organisms, we derive a model of microhomology-mediated break-induced replication (MMBIR) for the origin of CNV and, ultimately, of LCRs. We propose that breakage of replication forks in stressed cells that are deficient in homologous recombination induces an aberrant repair process with features of break-induced replication (BIR). Under these circumstances, single-strand 3′ tails from broken replication forks will anneal with microhomology on any single-stranded DNA nearby, priming low-processivity polymerization with multiple template switches generating complex rearrangements, and eventual re-establishment of processive replication.

Evolutionary analysis of the highly dynamic CHEK2 duplicon in anthropoids

BACKGROUND

Segmental duplications (SDs) are euchromatic portions of genomic DNA (> or = 1 kb) that occur at more than one site within the genome, and typically share a high level of sequence identity (>90%). Approximately 5% of the human genome is composed of such duplicated sequences. Here we report the detailed investigation of CHEK2 duplications. CHEK2 is a multiorgan cancer susceptibility gene encoding a cell cycle checkpoint kinase acting in the DNA-damage response signalling pathway. The continuous presence of the CHEK2 gene in all eukaryotes and its important role in maintaining genome stability prompted us to investigate the duplicative evolution and phylogeny of CHEK2 and its paralogs during anthropoid evolution.

RESULTS

To study CHEK2 duplicon evolution in anthropoids we applied a combination of comparative FISH and in silico analyses. Our comparative FISH results with a CHEK2 fosmid probe revealed the single-copy status of CHEK2 in New World monkeys, Old World monkeys and gibbons. Whereas a single CHEK2 duplication was detected in orangutan, a multi-site signal pattern indicated a burst of duplication in African great apes and human. Phylogenetic analysis of paralogous and ancestral CHEK2 sequences in human, chimpanzee and rhesus macaque confirmed this burst of duplication, which occurred after the radiation of orangutan and African great apes. In addition, we used inter-species quantitative PCR to determine CHEK2 copy numbers. An amplification of CHEK2 was detected in African great apes and the highest CHEK2 copy number of all analysed species was observed in the human genome. Furthermore, we detected variation in CHEK2 copy numbers within the analysed set of human samples.

CONCLUSION

Our detailed analysis revealed the highly dynamic nature of CHEK2 duplication during anthropoid evolution. We determined a burst of CHEK2 duplication after the radiation of orangutan and African great apes and identified the highest CHEK2 copy number in human. In conclusion, our analysis of CHEK2 duplicon evolution revealed that SDs contribute to inter-species variation. Furthermore, our qPCR analysis led us to presume CHEK2 copy number variation in human, and molecular diagnostics of the cancer susceptibility gene CHEK2 inside the duplicated region might be hampered by the individual-specific set of duplicons.

Chromosomal evolution of the PKD1 gene family in primates

BACKGROUND

The autosomal dominant polycystic kidney disease (ADPKD) is mostly caused by mutations in the PKD1 (polycystic kidney disease 1) gene located in 16p13.3. Moreover, there are six pseudogenes of PKD1 that are located proximal to the master gene in 16p13.1. In contrast, no pseudogene could be detected in the mouse genome, only a single copy gene on chromosome 17. The question arises how the human situation originated phylogenetically. To address this question we applied comparative FISH-mapping of a human PKD1-containing genomic BAC clone and a PKD1-cDNA clone to chromosomes of a variety of primate species and the dog as a non-primate outgroup species.

RESULTS

Comparative FISH with the PKD1-cDNA clone clearly shows that in all primate species studied distinct single signals map in subtelomeric chromosomal positions orthologous to the short arm of human chromosome 16 harbouring the master PKD1 gene. Only in human and African great apes, but not in orangutan, FISH with both BAC and cDNA clones reveals additional signal clusters located proximal of and clearly separated from the PKD1 master genes indicating the chromosomal position of PKD1 pseudogenes in 16p of these species, respectively. Indeed, this is in accordance with sequencing data in human, chimpanzee and orangutan. Apart from the master PKD1 gene, six pseudogenes are identified in both, human and chimpanzee, while only a single-copy gene is present in the whole-genome sequence of orangutan. The phylogenetic reconstruction of the PKD1-tree reveals that all human pseudogenes are closely related to the human PKD1 gene, and all chimpanzee pseudogenes are closely related to the chimpanzee PKD1 gene. However, our statistical analyses provide strong indication that gene conversion events may have occurred within the PKD1 family members of human and chimpanzee, respectively.

CONCLUSION

PKD1 must have undergone amplification very recently in hominid evolution. Duplicative transposition of the PKD1 gene and further amplification and evolution of the PKD1 pseudogenes may have arisen in a common ancestor of Homo, Pan and Gorilla approximately 8 MYA. Reticulate evolutionary processes such as gene conversion and non-allelic homologous recombination (NAHR) may have resulted in concerted evolution of PKD1 family members in human and chimpanzee and, thus, simulate an independent evolution of the PKD1 pseudogenes from their master PKD1 genes in human and chimpanzee.

The genomic distribution of intraspecific and interspecific sequence divergence of human segmental duplications relative to human/chimpanzee chromosomal rearrangements

Background

It has been suggested that chromosomal rearrangements harbor the molecular footprint of the biological phenomena which they induce, in the form, for instance, of changes in the sequence divergence rates of linked genes. So far, all the studies of these potential associations have focused on the relationship between structural changes and the rates of evolution of single-copy DNA and have tried to exclude segmental duplications (SDs). This is paradoxical, since SDs are one of the primary forces driving the evolution of structure and function in our genomes and have been linked not only with novel genes acquiring new functions, but also with overall higher DNA sequence divergence and major chromosomal rearrangements.

Results

Here we take the opposite view and focus on SDs. We analyze several of the features of SDs, including the rates of intraspecific divergence between paralogous copies of human SDs and of interspecific divergence between human SDs and chimpanzee DNA. We study how divergence measures relate to chromosomal rearrangements, while considering other factors that affect evolutionary rates in single copy DNA.

Conclusion

We find that interspecific SD divergence behaves similarly to divergence of single-copy DNA. In contrast, old and recent paralogous copies of SDs do present different patterns of intraspecific divergence. Also, we show that some relatively recent SDs accumulate in regions that carry inversions in sister lineages.

Structural differences of orthologous genes: insights from human-primate comparisons

The genomic basis of phenotypic distinctions between humans and nonhuman primates remains insufficiently explained. We hypothesized that interspecies structural differences of orthologous genes can cause such distinctions and searched protein-coding genes conserved between humans and nonhuman primates for species-specific initial and terminal exons. We inferred gene structure differences from genomic locations where portions of primate transcripts aligned with the human genome outside of any human exons. Of 22,466 high-confidence FANTOM3 human transcriptional units, 7424 (33%) had nonhuman primate full-length cDNA support. One hundred eighty-three of the loci contained 68,424 bp of sequence exonic in nonhuman primates but not humans. Fifty-four of 183 included species-specific portions of protein-coding regions. Six genes had evidence of intergenic splicing in a nonhuman primate but not in human. It is imperative that primate transcriptome projects be accelerated on par with genome projects to understand better interspecies gene structure distinctions.

Molecular mechanisms of chromosomal rearrangement during primate evolution

Breakpoint analysis of the large chromosomal rearrangements which have occurred during primate evolution promises to yield new insights into the underlying mechanisms of mutagenesis. Comparison of these evolutionary breakpoints with those that are disease-associated in humans, and which occur during either meiotic or mitotic cell division, should help to identify basic mechanistic similarities as well as differences. It has recently become clear that segmental duplications (SDs) have had a very significant impact on genome plasticity during primate evolution. In comparisons of the human and chimpanzee genomes, SDs have been found in flanking regions of 70-80% of inversions and approximately 40% of deletions/duplications. A strong spatial association between primate-specific breakpoints and SDs has also become evident from comparisons of human with other mammalian genomes. The lineage-specific hyperexpansion of certain SDs observed in the genomes of human, chimpanzee, gorilla and gibbon is indicative of the intrinsic instability of some SDs in primates. However, since many primate-specific breakpoints map to regions lacking SDs, but containing interspersed high-copy repetitive sequence elements such as SINEs, LINEs, LTRs, alpha-satellites and (AT)( n ) repeats, we may infer that a range of different molecular mechanisms have probably been involved in promoting chromosomal breakage during the evolution of primate genomes.

An evolution revolution provides further revelation

The extent of copy-number variation (CNV) in the human genome has been appreciated only recently. Nevertheless, for almost four decades, gene duplication has been a prevailing hypothesis for evolutionary change. Recently, gene CNV spanning 60 million years of human and primate evolution has been determined enabling lineage-specific gene CNV to be identified. Primate lineage-specific gene CNV studies reveal that almost one third of all human genes exhibit a copy-number change in one or more primate species. Intriguingly, human lineage-specific gene amplification can be correlated to the emergence of human-specific traits such as cognition and endurance running.

Segmental duplications and evolutionary plasticity at tumor chromosome break-prone regions

We have previously found that the borders of evolutionarily conserved chromosomal regions often coincide with tumor-associated deletion breakpoints within human 3p12-p22. Moreover, a detailed analysis of a frequently deleted region at 3p21.3 (CER1) showed associations between tumor breaks and gene duplications. We now report on the analysis of 54 chromosome 3 breaks by multipoint FISH (mpFISH) in 10 carcinoma-derived cell lines. The centromeric region was broken in five lines. In lines with highly complex karyotypes, breaks were clustered near known fragile sites, FRA3B, FRA3C, and FRA3D (three lines), and in two other regions: 3p12.3-p13 ( approximately 75 Mb position) and 3q21.3-q22.1 ( approximately 130 Mb position) (six lines). All locations are shown based on NCBI Build 36.1 human genome sequence. The last two regions participated in three of four chromosome 3 inversions during primate evolution. Regions at 75, 127, and 131 Mb positions carry a large ( approximately 250 kb) segmental duplication (tumor break-prone segmental duplication [TBSD]). TBSD homologous sequences were found at 15 sites on different chromosomes. They were located within bands frequently involved in carcinoma-associated breaks. Thirteen of them have been involved in inversions during primate evolution; 10 were reused by breaks during mammalian evolution; 14 showed copy number polymorphism in man. TBSD sites showed an increase in satellite repeats, retrotransposed sequences, and other segmental duplications. We propose that the instability of these sites stems from specific organization of the chromosomal region, associated with location at a boundary between different CG-content isochores and with the presence of TBSDs and "instability elements," including satellite repeats and retroviral sequences.

Pharmacogenomics in Alzheimer's disease

Pharmacological treatment in Alzheimer's disease (AD) accounts for 10-20% of direct costs, and fewer than 20% of AD patients are moderate responders to conventional drugs (donepezil, rivastigmine, galantamine, memantine), with doubtful cost-effectiveness. Both AD pathogenesis and drug metabolism are genetically regulated complex traits in which hundreds of genes cooperatively participate. Structural genomics studies demonstrated that more than 200 genes might be involved in AD pathogenesis regulating dysfunctional genetic networks leading to premature neuronal death. The AD population exhibits a higher genetic variation rate than the control population, with absolute and relative genetic variations of 40-60% and 0.85-1.89%, respectively. AD patients also differ in their genomic architecture from patients with other forms of dementia. Functional genomics studies in AD revealed that age of onset, brain atrophy, cerebrovascular hemodynamics, brain bioelectrical activity, cognitive decline, apoptosis, immune function, lipid metabolism dyshomeostasis, and amyloid deposition are associated with AD-related genes. Pioneering pharmacogenomics studies also demonstrated that the therapeutic response in AD is genotype-specific, with apolipoprotein E (APOE) 4/4 carriers the worst responders to conventional treatments. About 10-20% of Caucasians are carriers of defective cytochrome P450 (CYP) 2D6 polymorphic variants that alter the metabolism and effects of AD drugs and many psychotropic agents currently administered to patients with dementia. There is a moderate accumulation of AD-related genetic variants of risk in CYP2D6 poor metabolizers (PMs) and ultrarapid metabolizers (UMs), who are the worst responders to conventional drugs. The association of the APOE-4 allele with specific genetic variants of other genes (e.g., CYP2D6, angiotensin-converting enzyme [ACE]) negatively modulates the therapeutic response to multifactorial treatments affecting cognition, mood, and behavior. Pharmacogenetic and pharmacogenomic factors may account for 60-90% of drug variability in drug disposition and pharmacodynamics. The incorporation of pharmacogenetic/pharmacogenomic protocols to AD research and clinical practice can foster therapeutics optimization by helping to develop cost-effective pharmaceuticals and improving drug efficacy and safety.

Evolution of man in the light of molecular genetics: a review. Part I. Our evolutionary history and genomics

The discovery in the mid 1970s of efficient methods of DNA sequencing and their subsequent development into more and more rapid procedures followed by sequencing the genomes of many species, including man in 2001, revolutionised the whole of biology. Remarkably, new light could be cast on the evolutionary relations of different species, and the tempo and mode of evolution within a given species, notably man, could quantitatively be illuminated including ongoing evolution possibly involving also the size of the brains. This review is a short summary of the results of the molecular genetic investigations of human evolution including the time and place of the formation of our species, our evolutionary relation to the closest living species relatives as well as extinct forms of the genus Homo. The nature and amount of genetic polymorphism in man is also considered with special emphasis on the causes of this variation, and the role of natural selection in human evolution. A consensus about the mosaic nature of our genome and the rather dynamic structure of our ancestral population is gradually emerging. The modern gene pool has most likely been contributed to several different ancestral demes either before or after the emergence of the anatomically modern human phenotype in the extent that even the nature of the evolutionary lineage leading to the anatomically modern man as a distinct biological species is disputable. Regulation of the function of genes, as well as the evolution of brains will be dealt with in the second part of this review.