Complex Segmental Duplications Research Articles

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus.

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Chasing perfection: validation and polishing strategies for telomere-to-telomere genome assemblies.

Advances in long-read sequencing technologies and genome assembly methods have enabled the recent completion of the first telomere-to-telomere human genome assembly, which resolves complex segmental duplications and large tandem repeats, including centromeric satellite arrays in a complete hydatidiform mole (CHM13). Although derived from highly accurate sequences, evaluation revealed evidence of small errors and structural misassemblies in the initial draft assembly. To correct these errors, we designed a new repeat-aware polishing strategy that made accurate assembly corrections in large repeats without overcorrection, ultimately fixing 51% of the existing errors and improving the assembly quality value from 70.2 to 73.9 measured from PacBio high-fidelity and Illumina k-mers. By comparing our results to standard automated polishing tools, we outline common polishing errors and offer practical suggestions for genome projects with limited resources. We also show how sequencing biases in both high-fidelity and Oxford Nanopore Technologies reads cause signature assembly errors that can be corrected with a diverse panel of sequencing technologies.

Evidence for opposing selective forces operating on human-specific duplicated TCAF genes in Neanderthals and humans

TRP channel-associated factor 1/2 (TCAF1/TCAF2) proteins antagonistically regulate the cold-sensor protein TRPM8 in multiple human tissues. Understanding their significance has been complicated given the locus spans a gap-ridden region with complex segmental duplications in GRCh38. Using long-read sequencing, we sequence-resolve the locus, annotate full-length TCAF models in primate genomes, and show substantial human-specific TCAF copy number variation. We identify two human super haplogroups, H4 and H5, and establish that TCAF duplications originated ~1.7 million years ago but diversified only in Homo sapiens by recurrent structural mutations. Conversely, in all archaic-hominin samples the fixation for a specific H4 haplotype without duplication is likely due to positive selection. Here, our results of TCAF copy number expansion, selection signals in hominins, and differential TCAF2 expression between haplogroups and high TCAF2 and TRPM8 expression in liver and prostate in modern-day humans imply TCAF diversification among hominins potentially in response to cold or dietary adaptations.

The origin and evolution of fibromelanosis in domesticated chickens: Genomic comparison of Indonesian Cemani and Chinese Silkie breeds.

Like Chinese Silkie, Indonesian Ayam Cemani exhibits fibromelanosis or dermal hyperpigmentation and possesses complex segmental duplications on chromosome 20 that involve the endothelin 3 gene, EDN3. A genomic region, DR1 of 127 kb, together with another region, DR2 of 171 kb, was duplicated by unequal crossing over, accompanied by inversion of one DR2. Quantitative PCR and copy number variation analyses on the Cemani genome sequence confirmed the duplication of EDN3. These genetic arrangements are identical in Cemani and Silkie, indicating a single origin of the genetic cause of Fm. The two DR1s harbor two distinct EDN3 haplotypes in a form of permanent heterozygosity, although they remain allelic in the ancestral Red Jungle Fowl population and some domesticated chicken breeds, with their allelic divergence time being as recent as 0.3 million years ago. In Cemani and Silkie breeds, artificial selection favoring the Fm phenotype has left an unambiguous record for selective sweep that extends in both directions from tandemly duplicated EDN3 loci. This highly homozygous tract is different in length between Cemani and Silkie, reflecting their distinct breeding histories. It is estimated that the Fm phenotype came into existence at least 6600–9100 years ago, prior to domestication of Cemani and Silkie, and that throughout domestication there has been intense artificial selection with strength s > 50% in each breed.

Common structural features characterize interstitial intrachromosomal Xp and 18q triplications

Rare intrachromosomal triplications producing partial tetrasomies have been reported for a number of chromosomes. A detailed molecular characterization, necessary to define the mechanism of their formation, has so far been lacking. We report on the detailed clinical, cytogenetic, and molecular characterization of two triplications, one de novo involving chromosome 18q, the other familial on chromosome Xp. The clinical phenotype of the patient with 18q triplication, very likely due to overexpression of one or more of the genes in the region, consists mainly of facial dysmorphisms and developmental delay. The familial Xp triplication does not cause an increase in the number of copies of any gene and is almost certainly a polymorphism. The rearrangements are actually complex duplications/triplications. In both patients, their proximal breakpoints are located within complex segmental duplications, one containing the VCX gene cluster on chromosome Xp, the other the TCEB3 genes on chromosome 18q. A proximal duplicated region is also present in both patients. All junctions we analyzed were formed by non-homologous end joining (NHEJ). The structural features shared between our patients suggest the involvement of a common mechanism in the genesis of interstitial intrachromosomal triplications.

Complex Segmental Duplications Mediate a Recurrent dup(X)(p11.22-p11.23) Associated with Mental Retardation, Speech Delay, and EEG Anomalies in Males and Females

Submicroscopic copy-number variations make a considerable contribution to the genetic etiology of human disease. We have analyzed subjects with idiopathic mental retardation (MR) by using whole-genome oligonucleotide-based array comparative genomic hybridization (aCGH) and identified familial and de novo recurrent Xp11.22-p11.23 duplications in males and females with MR, speech delay, and a peculiar electroencephalographic (EEG) pattern in childhood. The size of the duplications ranges from 0.8-9.2 Mb. Most affected females show preferential activation of the duplicated X chromosome. Carriers of the smallest duplication show X-linked recessive inheritance. All other affected individuals present dominant expression and comparable clinical phenotypes irrespective of sex, duplication size, and X-inactivation pattern. The majority of the rearrangements are mediated by recombination between flanking complex segmental duplications. The identification of common clinical features, including the typical EEG pattern, predisposing genomic structure, and peculiar X-inactivation pattern, suggests that duplication of Xp11.22-p11.23 constitutes a previously undescribed syndrome.

Cost-effective sequence-based nonhuman primate MHC class I genotyping from RNA

The major histocompatibility complex (MHC) class I region of nonhuman primates (NHPs) is highly polymorphic and has undergone complex segmental duplications such that the number of expressed genes differs between individuals. In addition the relative abundance of transcripts varies more than 100-fold between NHP class I genes. This unparalleled complexity makes rapid, efficient class I genotyping difficult for NHPs. The ‘gold standard’ of cDNA library construction, screening and sequencing is both costly and labor-intensive. Several rapid genotyping methods have been utilized, but all require some degree of prior sequence knowledge. Here, we describe a method for sequence-based MHC class I genotyping which reduces cost by (1) pooling molecularly barcoded class I cDNA-PCR amplicons for cloning and (2) targeting sequencing of a region of concentrated polymorphism spanning the two exons encoding the peptide binding domain. This method can efficiently genotype both known and novel MHC class I alleles. In addition, full-length cDNA amplicons with novel sequences can be resequenced in their entireties to expand the repertoire of characterized MHC class I sequences for NHPs.

Characterization of six human disease-associated inversion polymorphisms

The human genome is a highly dynamic structure that shows a wide range of genetic polymorphic variation. Unlike other types of structural variation, little is known about inversion variants within normal individuals because such events are typically balanced and are difficult to detect and analyze by standard molecular approaches. Using sequence-based, cytogenetic and genotyping approaches, we characterized six large inversion polymorphisms that map to regions associated with genomic disorders with complex segmental duplications mapping at the breakpoints. We developed a metaphase FISH-based assay to genotype inversions and analyzed the chromosomes of 27 individuals from three HapMap populations. In this subset, we find that these inversions are less frequent or absent in Asians when compared with European and Yoruban populations. Analyzing multiple individuals from outgroup species of great apes, we show that most of these large inversion polymorphisms are specific to the human lineage with two exceptions, 17q21.31 and 8p23 inversions, which are found to be similarly polymorphic in other great ape species and where the inverted allele represents the ancestral state. Investigating linkage disequilibrium relationships with genotyped SNPs, we provide evidence that most of these inversions appear to have arisen on at least two different haplotype backgrounds. In these cases, discovery and genotyping methods based on SNPs may be confounded and molecular cytogenetics remains the only method to genotype these inversions.

Discovery of a previously unrecognized microdeletion syndrome of 16p11.2–p12.2

We have identified a recurrent de novo pericentromeric deletion in 16p11.2-p12.2 in four individuals with developmental disabilities by microarray-based comparative genomic hybridization analysis. The identification of common clinical features in these four individuals along with the characterization of complex segmental duplications flanking the deletion regions suggests that nonallelic homologous recombination mediated these rearrangements and that deletions in 16p11.2-p12.2 constitute a previously undescribed syndrome.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5'-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Analysis of copy number variation in the normal human population within a region containing complex segmental duplications on 22q11 using high-resolution array-CGH

A previously detected copy number polymorphism (Ep CNP) in patients affected with neuroectodermal tumors led us to investigate its frequency and length in the normal population. For this purpose, a program called Sequence Allocator was developed and applied for the construction of an array that consisted of unique and duplicated fragments, allowing the assessment of copy number variation within regions of segmental duplications. The average resolution of this array was 11 kb and we determined the size of the Ep CNP to be 290 kb. Analysis of normal controls identified 7.7 and 7.1% gains in peripheral blood and lymphoblastoid cell line (LCL) DNA, respectively, while deletions were found only in the LCL group (7.1%). This array platform allows the detection of DNA copy number variation within regions of pronounced genomic complexity, which constitutes an improvement over available technologies.

Evolutionary mechanisms shaping the genomic structure of the Williams-Beuren syndrome chromosomal region at human 7q11.23.

About 5% of the human genome consists of segmental duplications or low-copy repeats, which are large, highly homologous (>95%) fragments of sequence. It has been estimated that these segmental duplications emerged during the past approximately 35 million years (Myr) of human evolution and that they correlate with chromosomal rearrangements. Williams-Beuren syndrome (WBS) is a segmental aneusomy syndrome that is the result of a frequent de novo deletion at 7q11.23, mediated by large (approximately 400-kb) region-specific complex segmental duplications composed of different blocks. We have precisely defined the structure of the segmental duplications on human 7q11.23 and characterized the copy number and structure of the orthologous regions in other primates (macaque, orangutan, gorilla, and chimpanzee). Our data indicate a recent origin and rapid evolution of the 7q11.23 segmental duplications, starting before the diversification of hominoids (approximately 12-16 million years ago [Mya]), with species-specific duplications and intrachromosomal rearrangements that lead to significant differences among those genomes. Alu sequences are located at most edges of the large hominoid-specific segmental duplications, suggesting that they might have facilitated evolutionary rearrangements. We propose a mechanistic model based on Alu-mediated duplicated transposition along with nonallelic homologous recombination for the generation and local expansion of the segmental duplications. The extraordinary rate of evolutionary turnover of this region, rich in segmental duplications, results in important genomic variation among hominoid species, which could be of functional relevance and predispose to disease.