Members of the olfactory receptor gene family are contained in large blocks of DNA duplicated polymorphically near the ends of human chromosomes

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Characterization of Host-Specific Genes from Pine- and Grass-Associated Species of the Fusarium fujikuroi Species Complex

The Fusarium fujikuroi species complex (FFSC) includes socioeconomically important pathogens that cause disease for numerous crops and synthesize a variety of secondary metabolites that can contaminate feedstocks and food. Here, we used comparative genomics to elucidate processes underlying the ability of pine-associated and grass-associated FFSC species to colonize tissues of their respective plant hosts. We characterized the identity, possible functions, evolutionary origins, and chromosomal positions of the host-range-associated genes encoded by the two groups of fungi. The 72 and 47 genes identified as unique to the respective genome groups were potentially involved in diverse processes, ranging from transcription, regulation, and substrate transport through to virulence/pathogenicity. Most genes arose early during the evolution of Fusarium/FFSC and were only subsequently retained in some lineages, while some had origins outside Fusarium. Although differences in the densities of these genes were especially noticeable on the conditionally dispensable chromosome of F. temperatum (representing the grass-associates) and F. circinatum (representing the pine-associates), the host-range-associated genes tended to be located towards the subtelomeric regions of chromosomes. Taken together, these results demonstrate that multiple mechanisms drive the emergence of genes in the grass- and pine-associated FFSC taxa examined. It also highlighted the diversity of the molecular processes potentially underlying niche-specificity in these and other Fusarium species.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Keeping Balance Between Genetic Stability and Plasticity at the Telomere and Subtelomere of Trypanosoma brucei

Telomeres, the nucleoprotein complexes at chromosome ends, are well-known for their essential roles in genome integrity and chromosome stability. Yet, telomeres and subtelomeres are frequently less stable than chromosome internal regions. Many subtelomeric genes are important for responding to environmental cues, and subtelomeric instability can facilitate organismal adaptation to extracellular changes, which is a common theme in a number of microbial pathogens. In this review, I will focus on the delicate and important balance between stability and plasticity at telomeres and subtelomeres of a kinetoplastid parasite, Trypanosoma brucei, which causes human African trypanosomiasis and undergoes antigenic variation to evade the host immune response. I will summarize the current understanding about T. brucei telomere protein complex, the telomeric transcript, and telomeric R-loops, focusing on their roles in maintaining telomere and subtelomere stability and integrity. The similarities and differences in functions and underlying mechanisms of T. brucei telomere factors will be compared with those in human and yeast cells.

The Power of Stress: The Telo-Hormesis Hypothesis

Adaptative response to stress is a strategy conserved across evolution to promote survival. In this context, the groundbreaking findings of Miroslav Radman on the adaptative value of changing mutation rates opened new avenues in our understanding of stress response. Inspired by this work, we explore here the putative beneficial effects of changing the ends of eukaryotic chromosomes, the telomeres, in response to stress. We first summarize basic principles in telomere biology and then describe how various types of stress can alter telomere structure and functions. Finally, we discuss the hypothesis of stress-induced telomere signaling with hormetic effects.

Cell wall protein variation, break-induced replication, and subtelomere dynamics in Candida glabrata

Candida glabrata is an opportunistic pathogen of humans, responsible for up to 30% of disseminated candidiasis. Adherence of C. glabrata to host cells is mediated by adhesin-like proteins (ALPs), about half of which are encoded in the subtelomeres. We performed a de novo assembly of two C. glabrata strains, BG2 and BG3993, using long single-molecule real-time (SMRT) reads, and constructed high-quality telomere-to-telomere assemblies of all 13 chromosomes to assess differences between C. glabrata strains. We documented variation between strains, and in agreement with earlier studies, found high (~0.5%-1%) frequencies of SNVs across the genome, including within subtelomeric regions. We documented changes in ALP gene structure and complement: there are large length differences in ALP genes in different strains, resulting from copy number variation in tandem repeats. We compared strains to characterize chromosome rearrangement events including within the poorly characterized subtelomeric regions. We show that rearrangements within the subtelomere regions all affect ALP-encoding genes, and 14/16 involve just the most terminal ALP gene. We present evidence that these rearrangements are mediated by break-induced replication. This study highlights the constrained nature of subtelomeric changes impacting ALP gene complement and subtelomere structure.

Hotspots of Human Mutation

Mutation of the human genome results in three classes of genomic variation: single nucleotide variants; short insertions or deletions; and large structural variants (SVs). Some mutations occur during normal processes, such as meiotic recombination or B cell development, and others result from DNA replication or aberrant repair of breaks in sequence-specific contexts. Regardless of mechanism, mutations are subject to selection, and some hotspots can manifest in disease. Here, we discuss genomic regions prone to mutation, mechanisms contributing to mutation susceptibility, and the processes leading to their accumulation in normal and somatic genomes. With further, more accurate human genome sequencing, additional mutation hotspots, mechanistic details of their formation, and the relevance of hotspots to evolution and disease are likely to be discovered.

Comprehensive Analysis of Human Subtelomeres by Whole Genome Mapping

Detailed comprehensive knowledge of the structures of individual long-range telomere-terminal haplotypes are needed to understand their impact on telomere function, and to delineate the population structure and evolution of subtelomere regions. However, the abundance of large evolutionarily recent segmental duplications and high levels of large structural variations have complicated both the mapping and sequence characterization of human subtelomere regions. Here, we use high throughput optical mapping of large single DNA molecules in nanochannel arrays for 154 human genomes from 26 populations to present a comprehensive look at human subtelomere structure and variation. The results catalog many novel long-range subtelomere haplotypes and determine the frequencies and contexts of specific subtelomeric duplicons on each chromosome arm, helping to clarify the currently ambiguous nature of many specific subtelomere structures as represented in the current reference sequence (HG38). The organization and content of some duplicons in subtelomeres appear to show both chromosome arm and population-specific trends. Based upon these trends we estimate a timeline for the spread of these duplication blocks.

Genome maps across 26 human populations reveal population-specific patterns of structural variation

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

Large structural variants (SV) are understudied in human genetics research because of the difficulty to detect them in the routinely generated short-read sequencing data. Here, the authors generate optical genome maps of 154 individuals from 26 populations that allow comprehensive examination of large SVs.

Telomeres and genomic evolution

The terminal regions of eukaryotic chromosomes, composed of telomere repeat sequences and sub-telomeric sequences, represent some of the most variable and rapidly evolving regions of the genome. The sub-telomeric regions are characterized by segmentally duplicated repetitive DNA elements, interstitial telomere repeat sequences and families of variable genes. Sub-telomeric repeat sequence families are shared among multiple chromosome ends, often rendering detailed sequence characterization difficult. These regions are composed of constitutive heterochromatin and are subjected to high levels of meiotic recombination. Dysfunction within telomere repeat arrays, either due to disruption in the chromatin structure or because of telomere shortening, can lead to chromosomal fusion and the generation of large-scale genomic rearrangements across the genome. The dynamic nature of telomeric regions, therefore, provides functionally useful variation to create genetic diversity, but also provides a mechanism for rapid genomic evolution that can lead to reproductive isolation and speciation. This article is part of the theme issue 'Understanding diversity in telomere dynamics'.This article is part of the theme issue 'Understanding diversity in telomere dynamics'.

The DNA damage response at dysfunctional telomeres, and at interstitial and subtelomeric DNA double-strand breaks

In mammals, DNA double-strand breaks (DSBs) are primarily repaired by classical non-homologous end joining (C-NHEJ), although homologous recombination repair and alternative NHEJ (A-NHEJ), which involve DSB processing, can also occur. These pathways are tightly regulated to maintain chromosome integrity. The ends of chromosomes, called telomeres, contain telomeric DNA that forms a cap structure in cooperation with telomeric proteins to prevent the activation of the DNA damage response and chromosome fusion at chromosome termini. Telomeres and subtelomeric regions are poor substrates for DNA replication; therefore, regions near telomeres are prone to replication fork stalling and chromosome breakage. Moreover, DSBs near telomeres are poorly repaired. As a result, when DSBs occur near telomeres in normal cells, the cells stop proliferating, while in cancer cells, subtelomeric DSBs induce rearrangements due to the absence of cell cycle checkpoints. The sensitivity of subtelomeric regions to DSBs is due to the improper regulation of processing, because although C-NHEJ is functional at subtelomeric DSBs, excessive processing results in an increased frequency of large deletions and chromosome rearrangements involving A-NHEJ.

Diversification of spatiotemporal expression and copy number variation of the echinoid hbox12/pmar1/micro1 multigene family

Changes occurring during evolution in the cis-regulatory landscapes of individual members of multigene families might impart diversification in their spatiotemporal expression and function. The archetypal member of the echinoid hbox12/pmar1/micro1 family is hbox12-a, a homeobox-containing gene expressed exclusively by dorsal blastomeres, where it governs the dorsal/ventral gene regulatory network during embryogenesis of the sea urchin Paracentrotus lividus. Here we describe the inventory of the hbox12/pmar1/micro1 genes in P. lividus, highlighting that gene copy number variation occurs across individual sea urchins of the same species. We show that the various hbox12/pmar1/micro1 genes group into three subfamilies according to their spatiotemporal expression, which ranges from broad transcription throughout development to transient expression in either the animal hemisphere or micromeres of the early embryo. Interestingly, the promoter regions of those genes showing comparable expression patterns are highly similar, while differing from those of the other subfamilies. Strikingly, phylogenetic analysis suggests that the hbox12/pmar1/micro1 genes are species-specific, exhibiting extensive divergence in their noncoding, but not in their coding, sequences across three distinct sea urchin species. In spite of this, two micromere-specific genes of P. lividus possess a TCF/LEF-binding motif in a similar position, and their transcription relies on Wnt/β-catenin signaling, similar to the pmar1 and micro1 genes, which in other sea urchin species are involved in micromere specification. Altogether, our findings suggest that the hbox12/pmar1/micro1 gene family evolved rather rapidly, generating paralogs whose cis-regulatory sequences diverged following multiple rounds of duplication from a common ancestor.

Trypanosoma brucei TIF2 suppresses VSG switching by maintaining subtelomere integrity

Subtelomeres consist of sequences adjacent to telomeres and contain genes involved in important cellular functions, as subtelomere instability is associated with several human diseases. Balancing between subtelomere stability and plasticity is particularly important for Trypanosoma brucei, a protozoan parasite that causes human African trypanosomiasis. T. brucei regularly switches its major variant surface antigen, variant surface glycoprotein (VSG), to evade the host immune response, and VSGs are expressed exclusively from subtelomeres in a strictly monoallelic fashion. Telomere proteins are important for protecting chromosome ends from illegitimate DNA processes. However, whether they contribute to subtelomere integrity and stability has not been well studied. We have identified a novel T. brucei telomere protein, T. brucei TRF-Interacting Factor 2 (TbTIF2), as a functional homolog of mammalian TIN2. A transient depletion of TbTIF2 led to an elevated VSG switching frequency and an increased amount of DNA double-strand breaks (DSBs) in both active and silent subtelomeric bloodstream form expression sites (BESs). Therefore, TbTIF2 plays an important role in VSG switching regulation and is important for subtelomere integrity and stability. TbTIF2 depletion increased the association of TbRAD51 with the telomeric and subtelomeric chromatin, and TbRAD51 deletion further increased subtelomeric DSBs in TbTIF2-depleted cells, suggesting that TbRAD51-mediated DSB repair is the underlying mechanism of subsequent VSG switching. Surprisingly, significantly more TbRAD51 associated with the active BES than with the silent BESs upon TbTIF2 depletion, and TbRAD51 deletion induced much more DSBs in the active BES than in the silent BESs in TbTIF2-depleted cells, suggesting that TbRAD51 preferentially repairs DSBs in the active BES.

Sex bias in copy number variation of olfactory receptor gene family depends on ethnicity

Gender plays a pivotal role in the human genetic identity and is also manifested in many genetic disorders particularly mental retardation. In this study its effect on copy number variation (CNV), known to cause genetic disorders was explored. As the olfactory receptor (OR) repertoire comprises the largest human gene family, it was selected for this study, which was carried out within and between three populations, derived from 150 individuals from the 1000 Genome Project. Analysis of 3872 CNVs detected among 791 OR loci, in which 307 loci showed CNV, revealed the following novel findings: Sex bias in CNV was significantly more prevalent in uncommon than common CNV variants of OR pseudogenes, in which the male genome showed more CNVs; and in one-copy number loss compared to complete deletion of OR pseudogenes; both findings implying a more recent evolutionary role for gender. Sex bias in copy number gain was also detected. Another novel finding was that the observed sex bias was largely dependent on ethnicity and was in general absent in East Asians. Using a CNV public database for sick children (International Standard Cytogenomic Array Consortium) the application of these findings for improving clinical molecular diagnostics is discussed by showing an example of sex bias in CNV among kids with autism. Additional clinical relevance is discussed, as the most polymorphic CNV-enriched OR cluster in the human genome, located on chr 15q11.2, is found near the Prader–Willi syndrome/Angelman syndrome bi-directionally imprinted region associated with two well-known mental retardation syndromes. As olfaction represents the primitive cognition in most mammals, arguably in competition with the development of a larger brain, the extensive retention of OR pseudogenes in females of this study, might point to a parent-of-origin indirect regulatory role for OR pseudogenes in the embryonic development of human brain. Thus any perturbation in the temporal regulation of olfactory system could lead to developmental delay disorders including mental retardation.

Structural variation in subtelomeres

Subtelomeres are an incredibly dynamic part of the human genome located at the ends of chromosomes just proximal to telomere repeats. Although subtelomeric variation contributes to normal polymorphism in the human genome and is a by-product of rapid evolution in these regions, rearrangements in subtelomeres can also cause intellectual disabilities and birth defects, making robust methods of detecting copy number variation in chromosome ends a must for cytogenetics labs. In recent years, methods for detecting structural variation in subtelomeres have moved from fluorescence in situ hybridization (FISH) to array technology; however, FISH is still necessary to determine the chromosomal structure of subtelomeric gains and losses identified by arrays.

Introduction to telomeres and telomerase

Telomeres are ends of chromosomes that play an important part in the biology of eukaryotic cells. Through the coordinated action of the telomerase and networks of other proteins and factors, the length and integrity of telomeres are maintained to prevent telomere dysfunction that has been linked to senescence, aging, diseases, and cancer. The tools and assays being used to study telomeres are being broadened, which has allowed us to derive a more detailed, high-resolution picture of the various players and pathways at work at the telomeres.

Accidental amplification and inactivation of a methyltransferase gene eliminates cytosine methylation in Mycosphaerella graminicola

A de novo search for repetitive elements in the genome sequence of the wheat pathogen Mycosphaerella graminicola identified a family of repeats containing a DNA cytosine methyltransferase sequence (MgDNMT). All 23 MgDNMT sequences identified carried signatures of repeat induced point mutation (RIP). All copies were subtelomeric in location except for one on chromosome 6. Synteny with M. fijiensis implied that the nontelomeric copy on chromosome 6 served as a template for subsequent amplifications. Southern analysis revealed that the MgDNMT sequence also was amplified in 15 additional M. graminicola isolates from various geographical regions. However, this amplification event was specific to M. graminicola; a search for MgDNMT homologs identified only a single, unmutated copy in the genomes of 11 other ascomycetes. A genome-wide methylation assay revealed that M. graminicola lacks cytosine methylation, as expected if its MgDNMT gene is inactivated. Methylation was present in several other species tested, including the closest known relatives of M. graminicola, species S1 and S2. Therefore, the observed changes most likely occurred within the past 10,500 years since the divergence between M. graminicola and S1. Our data indicate that the recent amplification of a single-copy MgDNMT gene made it susceptible to RIP, resulting in complete loss of cytosine methylation in M. graminicola.

Replication timing of human telomeres is chromosome arm-specific, influenced by subtelomeric structures and connected to nuclear localization

The mechanisms governing telomere replication in humans are still poorly understood. To fill this gap, we investigated the timing of replication of single telomeres in human cells. Using in situ hybridization techniques, we have found that specific telomeres have preferential time windows for replication during the S-phase and that these intervals do not depend upon telomere length and are largely conserved between homologous chromosomes and between individuals, even in the presence of large subtelomeric segmental polymorphisms. Importantly, we show that one copy of the 3.3 kb macrosatellite repeat D4Z4, present in the subtelomeric region of the late replicating 4q35 telomere, is sufficient to confer both a more peripheral localization and a later-replicating property to a de novo formed telomere. Also, the presence of β-satellite repeats next to a newly created telomere is sufficient to delay its replication timing. Remarkably, several native, non-D4Z4–associated, late-replicating telomeres show a preferential localization toward the nuclear periphery, while several early-replicating telomeres are associated with the inner nuclear volume. We propose that, in humans, chromosome arm–specific subtelomeric sequences may influence both the spatial distribution of telomeres in the nucleus and their replication timing.

Functional telomeres are essential for genome stability. While replication of telomeres has been extensively studied in model organisms such as the baker's yeast, little is known about the mechanisms that govern the replication of human telomeres. In this study, we have determined the timing of replication of telomeres of individual human chromosomes and its association with potential modulating factors such as particular subtelomeric sequences, the presence of heterochromatic regions, and nuclear localization. We have found that native telomeres associated with D4Z4 sequences—a macrosatellite naturally located in the subtelomeric regions of 4q, 10q, and acrocentric chromosome extremities—replicate later than others. We also present descriptive and experimental evidence indicating that nuclear localization influences the timing of telomere replication. These results contribute to our understanding of telomere metabolism in humans.

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.

High-resolution SNP arrays in mental retardation diagnostics: how much do we gain?

We used Affymetrix 6.0 GeneChip SNP arrays to characterize copy number variations (CNVs) in a cohort of 70 patients previously characterized on lower-density oligonucleotide arrays affected by idiopathic mental retardation and dysmorphic features. The SNP array platform includes approximately 900,000 SNP probes and 900,000 non-SNP oligonucleotide probes at an average distance of 0.7 Kb, which facilitates coverage of the whole genome, including coding and noncoding regions. The high density of probes is critical for detecting small CNVs, but it can lead to data interpretation problems. To reduce the number of false positives, parameters were set to consider only imbalances >75 Kb encompassing at least 80 probe sets. The higher resolution of the SNP array platform confirmed the increased ability to detect small CNVs, although more than 80% of these CNVs overlapped to copy number 'neutral' polymorphism regions and 4.4% of them did not contain known genes. In our cohort of 70 patients, of the 51 previously evaluated as 'normal' on the Agilent 44K array, the SNP array platform disclosed six additional CNV changes, including three in three patients, which may be pathogenic. This suggests that about 6% of individuals classified as 'normal' using the lower-density oligonucleotide array could be found to be affected by a genomic disorder when evaluated with the higher-density microarray platforms.

DDX11L: a novel transcript family emerging from human subtelomeric regions

BACKGROUND

The subtelomeric regions of human chromosomes exhibit an extraordinary plasticity. To date, due to the high GC content and to the presence of telomeric repeats, the subtelomeric sequences are underrepresented in the genomic libraries and consequently their sequences are incomplete in the finished human genome sequence, and still much remains to be learned about subtelomere organization, evolution and function. Indeed, only in recent years, several studies have disclosed, within human subtelomeres, novel gene family members.

RESULTS

During a project aimed to analyze genes located in the telomeric region of the long arm of the human X chromosome, we have identified a novel transcript family, DDX11L, members of which map to 1pter, 2q13/14.1, 2qter, 3qter, 6pter, 9pter/9qter, 11pter, 12pter, 15qter, 16pter, 17pter, 19pter, 20pter/20qter, Xpter/Xqter and Yqter. Furthermore, we partially sequenced the underrepresented subtelomeres of human chromosomes showing a common evolutionary origin.

CONCLUSION

Our data indicate that an ancestral gene, originated as a rearranged portion of the primate DDX11 gene, and propagated along many subtelomeric locations, is emerging within subtelomeres of human chromosomes, defining a novel gene family. These findings support the possibility that the high plasticity of these regions, sites of DNA exchange among different chromosomes, could trigger the emergence of new genes.

Human copy number polymorphic genes

Recent large-scale genomic studies within human populations have identified numerous genomic regions as copy number variant (CNV). As these CNV regions often overlap coding regions of the genome, large lists of potentially copy number polymorphic genes have been produced that are candidates for disease association. Most of the current data regarding normal genic variation, however, has been generated using BAC or SNP microarrays, which lack precision especially with respect to exons. To address this, we assessed 2,790 candidate CNV genes defined from available studies in nine well-characterized HapMap individuals by designing a customized oligonucleotide microarray targeted specifically to exons. Using exon array comparative genomic hybridization (aCGH), we detected 255 (9%) of the candidates as true CNVs including 134 with evidence of variation over the entire gene. Individuals differed in copy number from the control by an average of 100 gene loci. Both partial- and whole-gene CNVs were strongly associated with segmental duplications (55 and 71%, respectively) as well as regions of positive selection. We confirmed 37% of the whole-gene CNVs using the fosmid end sequence pair (ESP) structural variation map for these same individuals. If we modify the end sequence pair mapping strategy to include low-sequence identity ESPs (98-99.5%) and ESPs with an everted orientation, we can capture 82% of the missed genes leading to more complete ascertainment of structural variation within duplicated genes. Our results indicate that segmental duplications are the source of the majority of full-length copy number polymorphic genes, most of the variant genes are organized as tandem duplications, and a significant fraction of these genes will represent paralogs with levels of sequence diversity beyond thresholds of allelic variation. In addition, these data provide a targeted set of CNV genes enriched for regions likely to be associated with human phenotypic differences due to copy number changes and present a source of copy number responsive oligonucleotide probes for future association studies.

Human subtelomeric copy number variations

Copy number variation is a defining characteristic of human subtelomeres. Human subtelomeric segmental duplication regions ('Subtelomeric Repeats') comprise about 25% of the most distal 500 kb and 80% of the most distal 100 kb in human DNA. Huge allelic disparities seen in subtelomeric DNA sequence content and organization are postulated to have an impact on the dosage of transcripts embedded within the duplicated sequences, on the transcription of genes in adjacent single copy DNA regions, and on the chromatin structures mediating telomere functions including chromosome stability. In addition to the complex duplicon substructure and huge allelic variations in extended subtelomere regions, both copy number variation and alternative sequence organizations for DNA characterize the sequences immediately adjacent to terminal (TTAGGG)n tracts ('subterminal DNA'). The structural variation in subterminal DNA is likely to have important consequences for expression of subterminal transcripts such as a newly-discovered gene family encoding actin-interacting proteins and a non-coding telomeric repeat containing RNA (TERRA) transcript family critical for telomere integrity. Major immediate challenges include discovering the full extent and nature of subtelomeric structural and copy number variation in humans, and developing methods for tracking individual allelic variants in the context of total genomic DNA.

The evolutionary significance of copy number variation in the human genome

Copy number variation provides the raw material for gene family expansion and diversification, which is an important evolutionary force. Moreover, copy number variants (CNVs) can influence gene transcriptional and translational levels and have been associated with complex disease susceptibility. Therefore, natural selection may have affected at least some of the greater than one thousand CNVs thus far discovered among the genomes of phenotypically normal humans. While identifying and understanding particular instances of natural selection may shed light on important aspects of human evolutionary history, our ability to analyze CNVs in traditional population genetic frameworks has been limited. However, progress has been made by adapting some of these frameworks for use with copy number data. Moving forward, these efforts will be aided by non-human organism studies of the population genetics of copy number variation, and by more direct comparisons of within-species copy number variation and between-species copy number fixation.

Human genes involved in copy number variation: mechanisms of origin, functional effects and implications for disease

Copy number variants (CNVs) overlap over 7000 genes, many of which are pivotal in biological pathways. The implications of this are profound, with consequences for evolutionary studies, population genetics, gene function and human phenotype, including elucidation of genetic susceptibility to major common diseases, the heritability of which has thus far defied full explanation. Even though this research is still in its infancy, CNVs have already been associated with a number of monogenic, syndromic and complex diseases: the development of high throughput and high resolution techniques for CNV screening is likely to bring further new insights into the contribution of copy number variation to common diseases. Amongst genes overlapped by CNVs, significant enrichments for certain gene ontology categories have been identified, including those related to immune responses and interactions with the environment. Genes in both of these categories are thought to be important in evolutionary adaptation and to be particular targets of natural selection. Thus, a full appreciation of copy number variation may be important for our understanding of human evolution.

Extraordinary diversity of chemosensory receptor gene repertoires among vertebrates

Chemosensation (smell and taste) is important to the survival and reproduction of vertebrates and is mediated by specific bindings of odorants, pheromones, and tastants by chemoreceptors that are encoded by several large gene families. This review summarizes recent comparative genomic and evolutionary studies of vertebrate chemoreceptor genes. It focuses on the remarkable diversity of chemoreceptor gene repertoires in terms of gene number and gene sequence across vertebrates and the evolutionary mechanisms that are responsible for generating this diversity. We argue that the great among-species variation of chemoreceptor gene repertoires is a result of adaptations of individual species to their environments and diets.

Characterization of chromosome ends in the filamentous fungus Neurospora crassa

Telomeres and subtelomere regions have vital roles in cellular homeostasis and can facilitate niche adaptation. However, information on telomere/subtelomere structure is still limited to a small number of organisms. Prior to initiation of this project, the Neurospora crassa genome assembly contained only 3 of the 14 telomeres. The missing telomeres were identified through bioinformatic mining of raw sequence data from the genome project and from clones in new cosmid and plasmid libraries. Their chromosomal locations were assigned on the basis of paired-end read information and/or by RFLP mapping. One telomere is attached to the ribosomal repeat array. The remaining chromosome ends have atypical structures in that they lack distinct subtelomere domains or other sequence features that are associated with telomeres in other organisms. Many of the chromosome ends terminate in highly AT-rich sequences that appear to be products of repeat-induced point mutation, although most are not currently repeated sequences. Several chromosome termini in the standard Oak Ridge wild-type strain were compared to their counterparts in an exotic wild type, Mauriceville. This revealed that the sequences immediately adjacent to the telomeres are usually genome specific. Finally, despite the absence of many features typically found in the telomere regions of other organisms, the Neurospora chromosome termini still retain the dynamic nature that is characteristic of chromosome ends.

Human telomere structure and biology

Human telomeric DNA is complex and highly variable. Subterminal sequences are associated with cis-acting determinants of allele-specific (TTAGGG)n tract length regulation and may modulate susceptibility of (TTAGGG)n tracts to rapid deletion events. More extensive subtelomeric DNA tracts are filled with segmental duplications and segments that vary in copy number, leading to highly variable subtelomeric allele structures in the human population. RNA transcripts encoded in telomere regions include multicopy protein-encoding gene families and a variety of noncoding RNAs. One recently described family of (UUAGGG)n-containing subterminal RNAs appears to be critical for telomere integrity; these RNAs associate with telomeric chromatin and are regulated by RNA surveillance factors including human homologs of the yeast Est1p protein. An increasingly detailed and complete picture of telomeric DNA sequence organization and structural variation is essential for understanding and tracking allele-specific subterminal and subtelomeric features critical for human biology.

Comparative sequence analysis of primate subtelomeres originating from a chromosome fission event

Subtelomeres are concentrations of interchromosomal segmental duplications capped by telomeric repeats at the ends of chromosomes. The nature of the segments shared by different sets of human subtelomeres reflects their high rate of recent interchromosomal exchange. Here, we characterize the rearrangements incurred by the 15q subtelomere after it arose from a chromosome fission event in the common ancestor of great apes. We used FISH, sequencing of genomic clones, and PCR to map the breakpoint of this fission and track the fate of flanking sequence in human, chimpanzee, gorilla, orangutan, and macaque genomes. The ancestral locus, a cluster of olfactory receptor (OR) genes, lies internally on macaque chromosome 7. Sequence originating from this fission site is split between the terminus of 15q and the pericentromere of 14q in the great apes. Numerous structural rearrangements, including interstitial deletions and transfers of material to or from other subtelomeres, occurred subsequent to the fission, such that each species has a unique 15q structure and unique collection of ORs derived from the fission locus. The most striking rearrangement involved transfer of at least 200 kb from the fission-site region to the end of chromosome 4q, where much still resides in chimpanzee and gorilla, but not in human. This gross structural difference places the subtelomeric defect underlying facioscapulohumeral muscular dystrophy (FSHD) much closer to the telomere in human 4q than in the hybrid 4q-15q subtelomere of chimpanzee.

Extensive copy-number variation of the human olfactory receptor gene family

As much as a quarter of the human genome has been reported to vary in copy number between individuals, including regions containing about half of the members of the olfactory receptor (OR) gene family. We have undertaken a detailed study of copy-number variation of ORs to elucidate the selective and mechanistic forces acting on this gene family and the true impact of copy-number variation on human OR repertoires. We argue that the properties of copy-number variants (CNVs) and other sets of large genomic regions violate the assumptions of statistical methods that are commonly used in the assessment of gene enrichment. Using more appropriate methods, we provide evidence that OR enrichment in CNVs is not due to positive selection but is because of OR preponderance in segmentally duplicated regions, which are known to be frequently copy-number variable, and because purifying selection against CNVs is lower in OR-containing regions than in regions containing essential genes. We also combine multiplex ligation-dependent probe amplification (MLPA) and PCR to assay the copy numbers of 37 candidate CNV ORs in a panel of approximately 50 human individuals. We confirm copy-number variation of 18 ORs but find no variation in this human-diversity panel for 16 other ORs, highlighting the caveat that reported intervals often overrepresent true CNVs. The copy-number variation we describe is likely to underpin significant variation in olfactory abilities among human individuals. Finally, we show that both homology-based and homology-independent processes have played a recent role in remodeling the OR family.

Structural variation of the human genome

There is growing appreciation that the human genome contains significant numbers of structural rearrangements, such as insertions, deletions, inversions, and large tandem repeats. Recent studies have defined approximately 5% of the human genome as structurally variant in the normal population, involving more than 800 independent genes. We present a detailed review of the various structural rearrangements identified to date in humans, with particular reference to their influence on human phenotypic variation. Our current knowledge of the extent of human structural variation shows that the human genome is a highly dynamic structure that shows significant large-scale variation from the currently published genome reference sequence.

Molecular population genetics of Drosophila subtelomeric DNA

DNA sequence surveys in yeast and humans suggest that the forces shaping telomeric polymorphism and divergence are distinctly more dynamic than those in the euchromatic, gene-rich regions of the chromosomes. However, the generality of this pattern across outbreeding, multicellular eukaryotes has not been determined. To characterize the structure and evolution of Drosophila telomeres, we collected and analyzed molecular population genetics data from the X chromosome subtelomere in 58 lines of North American Drosophila melanogaster and 29 lines of African D. melanogaster. We found that Drosophila subtelomeres exhibit high levels of both structural and substitutional polymorphism relative to linked euchromatic regions. We also observed strikingly different patterns of variation in the North American and African samples. Moreover, our analyses of the polymorphism data identify a localized hotspot of recombination in the most-distal portion of the X subtelomere. While the levels of polymorphism decline sharply and in parallel with rates of crossing over per physical length over the distal first euchromatic megabase pairs of the X chromosome, our data suggest that they rise again sharply in the subtelomeric region (approximately 80 kbp). These patterns of historical recombination and geographic differentiation indicate that, similar to yeast and humans, Drosophila subtelomeric DNA is evolving very differently from euchromatic DNA.

The fine-scale and complex architecture of human copy-number variation

Despite considerable excitement over the potential functional significance of copy-number variants (CNVs), we still lack knowledge of the fine-scale architecture of the large majority of CNV regions in the human genome. In this study, we used a high-resolution array-based comparative genomic hybridization (aCGH) platform that targeted known CNV regions of the human genome at approximately 1 kb resolution to interrogate the genomic DNAs of 30 individuals from four HapMap populations. Our results revealed that 1020 of 1153 CNV loci (88%) were actually smaller in size than what is recorded in the Database of Genomic Variants based on previously published studies. A reduction in size of more than 50% was observed for 876 CNV regions (76%). We conclude that the total genomic content of currently known common human CNVs is likely smaller than previously thought. In addition, approximately 8% of the CNV regions observed in multiple individuals exhibited genomic architectural complexity in the form of smaller CNVs within larger ones and CNVs with interindividual variation in breakpoints. Future association studies that aim to capture the potential influences of CNVs on disease phenotypes will need to consider how to best ascertain this previously uncharacterized complexity.

Human subtelomeric duplicon structure and organization

The sequence divergence within subtelomeric duplicon families varies considerably, as does the organization of duplicon blocks at subtelomere alleles; a class of duplicon blocks was identified that are subtelomere-specific.

Background

Human subtelomeric segmental duplications ('subtelomeric repeats') comprise about 25% of the most distal 500 kb and 80% of the most distal 100 kb in human DNA. A systematic analysis of the duplication substructure of human subtelomeric regions was done in order to develop a detailed understanding of subtelomeric sequence organization and a nucleotide sequence-level characterization of subtelomeric duplicon families.

Results

The extent of nucleotide sequence divergence within subtelomeric duplicon families varies considerably, as does the organization of duplicon blocks at subtelomere alleles. Subtelomeric internal (TTAGGG)n-like tracts occur at duplicon boundaries, suggesting their involvement in the generation of the complex sequence organization. Most duplicons have copies at both subtelomere and non-subtelomere locations, but a class of duplicon blocks is identified that are subtelomere-specific. In addition, a group of six subterminal duplicon families are identified that, together with six single-copy telomere-adjacent segments, include all of the (TTAGGG)n-adjacent sequence identified so far in the human genome.

Conclusion

Identification of a class of duplicon blocks that is subtelomere-specific will facilitate high-resolution analysis of subtelomere repeat copy number variation as well as studies involving somatic subtelomere rearrangements. The significant levels of nucleotide sequence divergence within many duplicon families as well as the differential organization of duplicon blocks on subtelomere alleles may provide opportunities for allele-specific subtelomere marker development; this is especially true for subterminal regions, where divergence and organizational differences are the greatest. These subterminal sequence families comprise the immediate cis-elements for (TTAGGG)n tracts, and are prime candidates for subtelomeric sequences regulating telomere-specific (TTAGGG)n tract length in humans.

Genetic elucidation of human hyperosmia to isovaleric acid

The genetic basis of odorant-specific variations in human olfactory thresholds, and in particular of enhanced odorant sensitivity (hyperosmia), remains largely unknown. Olfactory receptor (OR) segregating pseudogenes, displaying both functional and nonfunctional alleles in humans, are excellent candidates to underlie these differences in olfactory sensitivity. To explore this hypothesis, we examined the association between olfactory detection threshold phenotypes of four odorants and segregating pseudogene genotypes of 43 ORs genome-wide. A strong association signal was observed between the single nucleotide polymorphism variants in OR11H7P and sensitivity to the odorant isovaleric acid. This association was largely due to the low frequency of homozygous pseudogenized genotype in individuals with specific hyperosmia to this odorant, implying a possible functional role of OR11H7P in isovaleric acid detection. This predicted receptor–ligand functional relationship was further verified using the Xenopus oocyte expression system, whereby the intact allele of OR11H7P exhibited a response to isovaleric acid. Notably, we also uncovered another mechanism affecting general olfactory acuity that manifested as a significant inter-odorant threshold concordance, resulting in an overrepresentation of individuals who were hyperosmic to several odorants. An involvement of polymorphisms in other downstream transduction genes is one possible explanation for this observation. Thus, human hyperosmia to isovaleric acid is a complex trait, contributed to by both receptor and other mechanisms in the olfactory signaling pathway.

Humans can accurately discern thousands of odors, yet there is considerable inter-individual variation in the ability to detect different odors, with individuals exhibiting low sensitivity (hyposmia), high sensitivity (hyperosmia), or even “blindness” (anosmia) to particular odors. Such differences are thought to stem from genetic differences in olfactory receptor (OR) genes, which encode proteins that initiate olfactory signaling. OR segregating pseudogenes, which have both functional and inactive alleles in the population, are excellent candidates for producing this olfactory phenotype diversity. Here, we provide evidence that a particular segregating OR gene is related to sensitivity to a sweaty odorant, isovaleric acid. We show that hypersensitivity towards this odorant is seen predominantly in individuals who carry at least one copy of the intact allele. Furthermore, we demonstrate that this hyperosmia is a complex trait, being driven by additional factors affecting general olfactory acuity. Our results highlight a functional role of segregating pseudogenes in human olfactory variability, and constitute a step towards deciphering the genetic basis of human olfactory variability.

Genetic epidemiology analysis reveals a multifaceted mechanism underlying enhanced olfactory sensitivity to the sweaty odor of isovaleric acid in humans.

Mutational and selective effects on copy-number variants in the human genome

Comprehensive descriptions of large insertion/deletion or segmental duplication polymorphisms (SDs) in the human genome have recently been generated. These annotations, known collectively as structural or copy-number variants (CNVs), include thousands of discrete genomic regions and span hundreds of millions of nucleotides. Here we review the genomic distribution of CNVs, which is strongly correlated with gene, repeat and segmental duplication content. We explore the evolutionary mechanisms giving rise to this nonrandom distribution, considering the available data on both human polymorphisms and the fixed changes that differentiate humans from other species. It is likely that mutational biases, selective effects and interactions between these forces all contribute substantially to the spectrum of human copy-number variation. Although defining these variants with nucleotide-level precision remains a largely unmet but critical challenge, our understanding of their potential medical impact and evolutionary importance is rapidly emerging.

Co-regulation of a large and rapidly evolving repertoire of odorant receptor genes

The olfactory system meets niche- and species-specific demands by an accelerated evolution of its odorant receptor repertoires. In this review, we describe evolutionary processes that have shaped olfactory and vomeronasal receptor gene families in vertebrate genomes. We emphasize three important periods in the evolution of the olfactory system evident by comparative genomics: the adaptation to land in amphibian ancestors, the decline of olfaction in primates, and the delineation of putative pheromone receptors concurrent with rodent speciation. The rapid evolution of odorant receptor genes, the sheer size of the repertoire, as well as their wide distribution in the genome, presents a developmental challenge: how are these ever-changing odorant receptor repertoires coordinated within the olfactory system? A central organizing principle in olfaction is the specialization of sensory neurons resulting from each sensory neuron expressing only ~one odorant receptor allele. In this review, we also discuss this mutually exclusive expression of odorant receptor genes. We have considered several models to account for co-regulation of odorant receptor repertoires, as well as discussed a new hypothesis that invokes important epigenetic properties of the system.

Hominoid lineage specific amplification of low-copy repeats on 22q11.2 (LCR22s) associated with velo-cardio-facial/digeorge syndrome

Segmental duplications or low-copy repeats (LCRs) constitute approximately 5% of the sequenced portion of the human genome and are associated with many human congenital anomaly disorders. The low-copy repeats on chromosome 22q11.2 (LCR22s) mediate chromosomal rearrangements resulting in deletions, duplications and translocations. The evolutionary mechanisms leading to LCR22 formation is unknown. Four genes, USP18, BCR, GGTLA and GGT, map adjacent to the LCR22s and pseudogene copies are located within them. It has been hypothesized that gene duplication occurred during primate evolution, followed by recombination events, forming pseudogene copies. We investigated whether gene duplication could be detected in non-human hominoid species. FISH mapping was performed using probes to the four functional gene loci. There was evidence for a single copy in humans but additional copies in hominoid species. We then compared LCR22 copy number using LCR22 FISH probes. Lineage specific LCR22 variation was detected in the hominoid species supporting the hypothesis. To independently validate initial findings, real time PCR, and screening of gorilla BAC library filters were performed. This was compared to array comparative genome hybridization data available. The most striking finding was a dramatic amplification of LCR22s in the gorilla. The LCR22s localized to the telomeric or subtelomeric bands of gorilla chromosomes. The most parsimonious explanation is that the LCR22s became amplified by inter-chromosomal recombination between telomeric bands. In summary, our results are consistent with a lineage specific coupling between gene and LCR22 duplication events. The LCR22s thus serve as an important model for evolution of genome variation.

Toward accurate high-throughput SNP genotyping in the presence of inherited copy number variation

Background

The recent discovery of widespread copy number variation in humans has forced a shift away from the assumption of two copies per locus per cell throughout the autosomal genome. In particular, a SNP site can no longer always be accurately assigned one of three genotypes in an individual. In the presence of copy number variability, the individual may theoretically harbor any number of copies of each of the two SNP alleles.

Results

To address this issue, we have developed a method to infer a "generalized genotype" from raw SNP microarray data. Here we apply our approach to data from 48 individuals and uncover thousands of aberrant SNPs, most in regions that were previously unreported as copy number variants. We show that our allele-specific copy numbers follow Mendelian inheritance patterns that would be obscured in the absence of SNP allele information. The interplay between duplication and point mutation in our data shed light on the relative frequencies of these events in human history, showing that at least some of the duplication events were recurrent.

Conclusion

This new multi-allelic view of SNPs has a complicated role in disease association studies, and further work will be necessary in order to accurately assess its importance. Software to perform generalized genotyping from SNP array data is freely available online [1].

What a difference copy number variation makes

DNA copy number variation (CNV) represents a considerable source of human genetic diversity. Recently,1 a global map of copy number variation in the human genome has been drawn up which reveals not only the ubiquity but also the complexity of this type of variation. Thus, two human genomes may differ by more than 20 Mb and it is likely that the full extent of CNV still remains to be discovered. Nearly 3000 genes are associated with CNV. This high degree of variability with regard to gene copy number between two individuals challenges definitions of normality. Many CNVs are located in regions of complex genomic structure and this currently limits the extent to which these variants can be genotyped by using tagging SNPs. However, some CNVs are already amenable to genome-wide association studies so that their influence on human phenotypic diversity and disease susceptibility may soon be determined.

Understanding the recent evolution of the human genome: insights from human-chimpanzee genome comparisons

The sequencing of the chimpanzee genome and the comparison with its human counterpart have begun to reveal the spectrum of genetic changes that has accompanied human evolution. In addition to gross karyotypic rearrangements such as the fusion that formed human chromosome 2 and the human-specific pericentric inversions of chromosomes 1 and 18, there is considerable submicroscopic structural variation involving deletions, duplications, and inversions. Lineage-specific segmental duplications, detected by array comparative genomic hybridization and direct sequence comparison, have made a very significant contribution to this structural divergence, which is at least three-fold greater than that due to nucleotide substitutions. Since structural genomic changes may have given rise to irreversible functional differences between the diverging species, their detailed analysis could help to identify the biological processes that have accompanied speciation. To this end, interspecies comparisons have revealed numerous human-specific gains and losses of genes as well as changes in gene expression. The very considerable structural diversity (polymorphism) evident within both lineages has, however, hampered the analysis of the structural divergence between the human and chimpanzee genomes. The concomitant evaluation of genetic divergence and diversity at the nucleotide level has nevertheless served to identify many genes that have evolved under positive selection and may thus have been involved in the development of human lineage-specific traits. Genes that display signs of weak negative selection have also been identified and could represent candidate loci for complex genomic disorders. Here, we review recent progress in comparing the human and chimpanzee genomes and discuss how the differences detected have improved our understanding of the evolution of the human genome.

Structural divergence between the human and chimpanzee genomes

The structural microheterogeneity evident between the human and chimpanzee genomes is quite considerable and includes inversions and duplications as well as deletions, ranging in size from a few base-pairs up to several megabases (Mb). Insertions and deletions have together given rise to at least 150 Mb of genomic DNA sequence that is either present or absent in humans as compared to chimpanzees. Such regions often contain paralogous sequences and members of multigene families thereby ensuring that the human and chimpanzee genomes differ by a significant fraction of their gene content. There is as yet no evidence to suggest that the large chromosomal rearrangements which serve to distinguish the human and chimpanzee karyotypes have influenced either speciation or the evolution of lineage-specific traits. However, the myriad submicroscopic rearrangements in both genomes, particularly those involving copy number variation, are unlikely to represent exclusively neutral changes and hence promise to facilitate the identification of genes that have been important for human-specific evolution.

Organization of chromosome ends in the rice blast fungus, Magnaporthe oryzae

Eukaryotic pathogens of humans often evade the immune system by switching the expression of surface proteins encoded by subtelomeric gene families. To determine if plant pathogenic fungi use a similar mechanism to avoid host defenses, we sequenced the 14 chromosome ends of the rice blast pathogen, Magnaporthe oryzae. One telomere is directly joined to ribosomal RNA-encoding genes, at the end of the ∼2 Mb rDNA array. Two are attached to chromosome-unique sequences, and the remainder adjoin a distinct subtelomere region, consisting of a telomere-linked RecQ-helicase (TLH) gene flanked by several blocks of tandem repeats. Unlike other microbes, M.oryzae exhibits very little gene amplification in the subtelomere regions—out of 261 predicted genes found within 100 kb of the telomeres, only four were present at more than one chromosome end. Therefore, it seems unlikely that M.oryzae uses switching mechanisms to evade host defenses. Instead, the M.oryzae telomeres have undergone frequent terminal truncation, and there is evidence of extensive ectopic recombination among transposons in these regions. We propose that the M.oryzae chromosome termini play more subtle roles in host adaptation by promoting the loss of terminally-positioned genes that tend to trigger host defenses.

A 360-kb interchromosomal duplication of the human HYDIN locus

The HYDIN gene located in human chromosome band 16q22.2 is a large gene encompassing 423 kb of genomic DNA that has been suggested as a candidate for an autosomal recessive form of congenital hydrocephalus. We have found that the human HYDIN locus has been very recently duplicated, with a nearly identical 360-kb paralogous segment inserted on chromosome 1q21.1. The duplication, among the largest interchromosomal segmental duplications described in humans, is not accounted for in the current human genome assembly and appears to be part of a greater than 550-kb contig that must lie within 1 of the 11 sequence gaps currently remaining in 1q21.1. Both copies of the HYDIN gene are expressed in alternatively spliced transcripts. Elucidation of the role of HYDIN in human disease susceptibility will require careful discrimination among the paralogous copies.