Identification of the transforming EML4-ALK fusion gene in non-small-cell lung cancer

Improvement in the clinical outcome of lung cancer is likely to be achieved by identification of the molecular events that underlie its pathogenesis. Here we show that a small inversion within chromosome 2p results in the formation of a fusion gene comprising portions of the echinoderm microtubule-associated protein-like 4 (EML4) gene and the anaplastic lymphoma kinase (ALK) gene in non-small-cell lung cancer (NSCLC) cells. Mouse 3T3 fibroblasts forced to express this human fusion tyrosine kinase generated transformed foci in culture and subcutaneous tumours in nude mice. The EML4-ALK fusion transcript was detected in 6.7% (5 out of 75) of NSCLC patients examined; these individuals were distinct from those harbouring mutations in the epidermal growth factor receptor gene. Our data demonstrate that a subset of NSCLC patients may express a transforming fusion kinase that is a promising candidate for a therapeutic target as well as for a diagnostic molecular marker in NSCLC.

Mapping and sequencing of structural variation from eight human genomes

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

Chromosome inversions, local adaptation and speciation

We study the evolution of inversions that capture locally adapted alleles when two populations are exchanging migrants or hybridizing. By suppressing recombination between the loci, a new inversion can spread. Neither drift nor coadaptation between the alleles (epistasis) is needed, so this local adaptation mechanism may apply to a broader range of genetic and demographic situations than alternative hypotheses that have been widely discussed. The mechanism can explain many features observed in inversion systems. It will drive an inversion to high frequency if there is no countervailing force, which could explain fixed differences observed between populations and species. An inversion can be stabilized at an intermediate frequency if it also happens to capture one or more deleterious recessive mutations, which could explain polymorphisms that are common in some species. This polymorphism can cycle in frequency with the changing selective advantage of the locally favored alleles. The mechanism can establish underdominant inversions that decrease heterokaryotype fitness by several percent if the cause of fitness loss is structural, while if the cause is genic there is no limit to the strength of underdominance that can result. The mechanism is expected to cause loci responsible for adaptive species-specific differences to map to inversions, as seen in recent QTL studies. We discuss data that support the hypothesis, review other mechanisms for inversion evolution, and suggest possible tests.

Comparative genome sequencing of Drosophila pseudoobscura: chromosomal, gene, and cis-element evolution

We have sequenced the genome of a second Drosophila species, Drosophila pseudoobscura, and compared this to the genome sequence of Drosophila melanogaster, a primary model organism. Throughout evolution the vast majority of Drosophila genes have remained on the same chromosome arm, but within each arm gene order has been extensively reshuffled, leading to a minimum of 921 syntenic blocks shared between the species. A repetitive sequence is found in the D. pseudoobscura genome at many junctions between adjacent syntenic blocks. Analysis of this novel repetitive element family suggests that recombination between offset elements may have given rise to many paracentric inversions, thereby contributing to the shuffling of gene order in the D. pseudoobscura lineage. Based on sequence similarity and synteny, 10,516 putative orthologs have been identified as a core gene set conserved over 25-55 million years (Myr) since the pseudoobscura/melanogaster divergence. Genes expressed in the testes had higher amino acid sequence divergence than the genome-wide average, consistent with the rapid evolution of sex-specific proteins. Cis-regulatory sequences are more conserved than random and nearby sequences between the species--but the difference is slight, suggesting that the evolution of cis-regulatory elements is flexible. Overall, a pattern of repeat-mediated chromosomal rearrangement, and high coadaptation of both male genes and cis-regulatory sequences emerges as important themes of genome divergence between these species of Drosophila.

Hyperphagia, severe obesity, impaired cognitive function, and hyperactivity associated with functional loss of one copy of the brain-derived neurotrophic factor (BDNF) gene

The neurotrophin brain-derived neurotrophic factor (BDNF) inhibits food intake, and rodent models of BDNF disruption all exhibit increased food intake and obesity, as well as hyperactivity. We report an 8-year-old girl with hyperphagia and severe obesity, impaired cognitive function, and hyperactivity who harbored a de novo chromosomal inversion, 46,XX,inv(11)(p13p15.3), a region encompassing the BDNF gene. We have identified the proximal inversion breakpoint that lies 850 kb telomeric of the 5' end of the BDNF gene. The patient's genomic DNA was heterozygous for a common coding polymorphism in BDNF, but monoallelic expression was seen in peripheral lymphocytes. Serum concentration of BDNF protein was reduced compared with age- and BMI-matched subjects. Haploinsufficiency for BDNF was associated with increased ad libitum food intake, severe early-onset obesity, hyperactivity, and cognitive impairment. These findings provide direct evidence for the role of the neurotrophin BDNF in human energy homeostasis, as well as in cognitive function, memory, and behavior.

The barley pan-genome reveals the hidden legacy of mutation breeding

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Highly multiplexed molecular inversion probe genotyping: over 10,000 targeted SNPs genotyped in a single tube assay

Large-scale genetic studies are highly dependent on efficient and scalable multiplex SNP assays. In this study, we report the development of Molecular Inversion Probe technology with four-color, single array detection, applied to large-scale genotyping of up to 12,000 SNPs per reaction. While generating 38,429 SNP assays using this technology in a population of 30 trios from the Centre d'Etude Polymorphisme Humain family panel as part of the International HapMap project, we established SNP conversion rates of approximately 90% with concordance rates >99.6% and completeness levels >98% for assays multiplexed up to 12,000plex levels. Furthermore, these individual metrics can be "traded off" and, by sacrificing a small fraction of the conversion rate, the accuracy can be increased to very high levels. No loss of performance is seen when scaling from 6,000plex to 12,000plex assays, strongly validating the ability of the technology to suppress cross-reactivity at high multiplex levels. The results of this study demonstrate the suitability of this technology for comprehensive association studies that use targeted SNPs in indirect linkage disequilibrium studies or that directly screen for causative mutations.

Functional Correction of Large Factor VIII Gene Chromosomal Inversions in Hemophilia A Patient-Derived iPSCs Using CRISPR-Cas9

Hemophilia A is an X-linked genetic disorder caused by mutations in the F8 gene, which encodes the blood coagulation factor VIII. Almost half of all severe hemophilia A cases result from two gross (140-kbp or 600-kbp) chromosomal inversions that involve introns 1 and 22 of the F8 gene, respectively. We derived induced pluripotent stem cells (iPSCs) from patients with these inversion genotypes and used CRISPR-Cas9 nucleases to revert these chromosomal segments back to the WT situation. We isolated inversion-corrected iPSCs with frequencies of up to 6.7% without detectable off-target mutations based on whole-genome sequencing or targeted deep sequencing. Endothelial cells differentiated from corrected iPSCs expressed the F8 gene and functionally rescued factor VIII deficiency in an otherwise lethal mouse model of hemophilia. Our results therefore provide a proof of principle for functional correction of large chromosomal rearrangements in patient-derived iPSCs and suggest potential therapeutic applications.

Efficient sorting of genomic permutations by translocation, inversion and block interchange

MOTIVATION

Finding genomic distance based on gene order is a classic problem in genome rearrangements. Efficient exact algorithms for genomic distances based on inversions and/or translocations have been found but are complicated by special cases, rare in simulations and empirical data. We seek a universal operation underlying a more inclusive set of evolutionary operations and yielding a tractable genomic distance with simple mathematical form.

RESULTS

We study a universal double-cut-and-join operation that accounts for inversions, translocations, fissions and fusions, but also produces circular intermediates which can be reabsorbed. The genomic distance, computable in linear time, is given by the number of breakpoints minus the number of cycles (b-c) in the comparison graph of the two genomes; the number of hurdles does not enter into it. Without changing the formula, we can replace generation and re-absorption of a circular intermediate by a generalized transposition, equivalent to a block interchange, with weight two. Our simple algorithm converts one multi-linear chromosome genome to another in the minimum distance.

Genome sequencing and comparative analysis of Saccharomyces cerevisiae strain YJM789

We sequenced the genome of Saccharomyces cerevisiae strain YJM789, which was derived from a yeast isolated from the lung of an AIDS patient with pneumonia. The strain is used for studies of fungal infections and quantitative genetics because of its extensive phenotypic differences to the laboratory reference strain, including growth at high temperature and deadly virulence in mouse models. Here we show that the approximately 12-Mb genome of YJM789 contains approximately 60,000 SNPs and approximately 6,000 indels with respect to the reference S288c genome, leading to protein polymorphisms with a few known cases of phenotypic changes. Several ORFs are found to be unique to YJM789, some of which might have been acquired through horizontal transfer. Localized regions of high polymorphism density are scattered over the genome, in some cases spanning multiple ORFs and in others concentrated within single genes. The sequence of YJM789 contains clues to pathogenicity and spurs the development of more powerful approaches to dissecting the genetic basis of complex hereditary traits.

Two chloroplast DNA inversions originated simultaneously during the early evolution of the sunflower family (Asteraceae)

The chloroplast DNA (cpDNA) inversion in the Asteraceae has been cited as a classic example of using genomic rearrangements for defining major lineages of plants. We further characterize cpDNA inversions in the Asteraceae using extensive sequence comparisons among 56 species, including representatives of all major clades of the family and related families. We determine the boundaries of the 22-kb (now known as 22.8 kb) inversion that defines a major split within the Asteraceae, and in the process, we characterize the second and a new, smaller 3.3-kb inversion that occurs at one end of the larger inversion. One end point of the smaller inversion is upstream of the trnE-UUC gene, and the other end point is located between the trnC-GCA and rpoB genes. Although a diverse sampling of Asteraceae experienced substantial length variation and base substitution during the long evolutionary history subsequent to the inversion events, the precise locations of the inversion end points are identified using comparative sequence alignments in the inversion regions. The phylogenetic distribution of two inversions is identical among the members of Asteraceae, suggesting that the inversion events likely occurred simultaneously or within a short time period shortly after the origin of the family. Estimates of divergence times based on ndhF and rbcL sequences suggest that two inversions originated during the late Eocene (38-42 MYA). The divergence time estimates also suggest that the Asteraceae originated in the mid Eocene (42-47 MYA).

Population genomics of inversion polymorphisms in Drosophila melanogaster

Chromosomal inversions have been an enduring interest of population geneticists since their discovery in Drosophila melanogaster. Numerous lines of evidence suggest powerful selective pressures govern the distributions of polymorphic inversions, and these observations have spurred the development of many explanatory models. However, due to a paucity of nucleotide data, little progress has been made towards investigating selective hypotheses or towards inferring the genealogical histories of inversions, which can inform models of inversion evolution and suggest selective mechanisms. Here, we utilize population genomic data to address persisting gaps in our knowledge of D. melanogaster's inversions. We develop a method, termed Reference-Assisted Reassembly, to assemble unbiased, highly accurate sequences near inversion breakpoints, which we use to estimate the age and the geographic origins of polymorphic inversions. We find that inversions are young, and most are African in origin, which is consistent with the demography of the species. The data suggest that inversions interact with polymorphism not only in breakpoint regions but also chromosome-wide. Inversions remain differentiated at low levels from standard haplotypes even in regions that are distant from breakpoints. Although genetic exchange appears fairly extensive, we identify numerous regions that are qualitatively consistent with selective hypotheses. Finally, we show that In(1)Be, which we estimate to be ∼60 years old (95% CI 5.9 to 372.8 years), has likely achieved high frequency via sex-ratio segregation distortion in males. With deeper sampling, it will be possible to build on our inferences of inversion histories to rigorously test selective models-particularly those that postulate that inversions achieve a selective advantage through the maintenance of co-adapted allele complexes.

Recombination-restarted replication makes inverted chromosome fusions at inverted repeats

Impediments to DNA replication are known to induce gross chromosomal rearrangements (GCRs) and copy-number variations (CNVs). GCRs and CNVs underlie human genomic disorders and are a feature of cancer. During cancer development, environmental factors and oncogene-driven proliferation promote replication stress. Resulting GCRs and CNVs are proposed to contribute to cancer development and therapy resistance. When stress arrests replication, the replisome remains associated with the fork DNA (stalled fork) and is protected by the inter-S-phase checkpoint. Stalled forks efficiently resume when the stress is relieved. However, if the polymerases dissociate from the fork (fork collapse) or the fork structure breaks (broken fork), replication restart can proceed either by homologous recombination or microhomology-primed re-initiation. Here we ascertain the consequences of replication with a fork restarted by homologous recombination in fission yeast. We identify a new mechanism of chromosomal rearrangement through the observation that recombination-restarted forks have a considerably high propensity to execute a U-turn at small inverted repeats (up to 1 in 40 replication events). We propose that the error-prone nature of restarted forks contributes to the generation of GCRs and gene amplification in cancer, and to non-recurrent CNVs in genomic disorders.

The Rose-comb mutation in chickens constitutes a structural rearrangement causing both altered comb morphology and defective sperm motility

Rose-comb, a classical monogenic trait of chickens, is characterized by a drastically altered comb morphology compared to the single-combed wild-type. Here we show that Rose-comb is caused by a 7.4 Mb inversion on chromosome 7 and that a second Rose-comb allele arose by unequal crossing over between a Rose-comb and wild-type chromosome. The comb phenotype is caused by the relocalization of the MNR2 homeodomain protein gene leading to transient ectopic expression of MNR2 during comb development. We also provide a molecular explanation for the first example of epistatic interaction reported by Bateson and Punnett 104 years ago, namely that walnut-comb is caused by the combined effects of the Rose-comb and Pea-comb alleles. Transient ectopic expression of MNR2 and SOX5 (causing the Pea-comb phenotype) occurs in the same population of mesenchymal cells and with at least partially overlapping expression in individual cells in the comb primordium. Rose-comb has pleiotropic effects, as homozygosity in males has been associated with poor sperm motility. We postulate that this is caused by the disruption of the CCDC108 gene located at one of the inversion breakpoints. CCDC108 is a poorly characterized protein, but it contains a MSP (major sperm protein) domain and is expressed in testis. The study illustrates several characteristic features of the genetic diversity present in domestic animals, including the evolution of alleles by two or more consecutive mutations and the fact that structural changes have contributed to fast phenotypic evolution.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

The latitudinal cline in the In(3R)Payne inversion polymorphism has shifted in the last 20 years in Australian Drosophila melanogaster populations

Clinal variation has been described in a number of inversions in Drosophila but these clines are often characterized by cytological techniques using small sample sizes, and associations with specific genes are rarely considered. Here we have developed a molecular assay for In(3R)Payne in Drosophila melanogaster from eastern Australia populations. It shows in repeated samples that the inversion cline is very tightly associated with latitude and is almost fixed in tropical populations while relatively rare in temperate populations. This steep cline has shifted in position in the last 20 years. The heat shock gene, hsr-omega, located centrally inside the inversion sequence, shows a different clinal pattern to In(3R)Payne. These results suggest strong ongoing selection on In(3R)Payne over the last 100 years since the colonization of Australia that is partly independent of hsr-omega.

Evaluation of the genomic extent of effects of fixed inversion differences on intraspecific variation and interspecific gene flow in Drosophila pseudoobscura and D. persimilis

There is increasing evidence that chromosomal inversions may facilitate the formation or persistence of new species by allowing genetic factors conferring species-specific adaptations or reproductive isolation to be inherited together and by reducing or eliminating introgression. However, the genomic domain of influence of the inverted regions on introgression has not been carefully studied. Here, we present a detailed study on the consequences that distance from inversion breakpoints has had on the inferred level of gene flow and divergence between Drosophila pseudoobscura and D. persimilis. We identified the locations of the inversion breakpoints distinguishing D. pseudoobscura and D. persimilis in chromosomes 2, XR, and XL. Population genetic data were collected at specific distances from the inversion breakpoints of the second chromosome and at two loci inside the XR and XL inverted regions. For loci outside the inverted regions, we found that distance from the nearest inversion breakpoint had a significant effect on several measures of divergence and gene flow between D. pseudoobscura and D. persimilis. The data fitted a logarithmic relationship, showing that the suppression of crossovers in inversion heterozygotes also extends to loci located outside the inversion but close to it (within 1-2 Mb). Further, we detected a significant reduction in nucleotide variation inside the inverted second chromosome region of D. persimilis and near one breakpoint, consistent with a scenario in which this inversion arose and was fixed in this species by natural selection.

DECKO: Single-oligo, dual-CRISPR deletion of genomic elements including long non-coding RNAs

BACKGROUND

CRISPR genome-editing technology makes it possible to quickly and cheaply delete non-protein-coding regulatory elements. We present a vector system adapted for this purpose called DECKO (Double Excision CRISPR Knockout), which applies a simple two-step cloning to generate lentiviral vectors expressing two guide RNAs (gRNAs) simultaneously. The key feature of DECKO is its use of a single 165 bp starting oligonucleotide carrying the variable sequences of both gRNAs, making it fully scalable from single-locus studies to complex library cloning.

RESULTS

We apply DECKO to deleting the promoters of one protein-coding gene and two oncogenic lncRNAs, UCA1 and the highly-expressed MALAT1, focus of many previous studies employing RNA interference approaches. DECKO successfully deleted genomic fragments ranging in size from 100 to 3000 bp in four human cell lines. Using a clone-derivation workflow lasting approximately 20 days, we obtained 9 homozygous and 17 heterozygous promoter knockouts in three human cell lines. Frequent target region inversions were observed. These clones have reductions in steady-state MALAT1 RNA levels of up to 98 % and display reduced proliferation rates.

CONCLUSIONS

We present a dual CRISPR tool, DECKO, which is cloned using a single starting oligonucleotide, thereby affording simplicity and scalability to CRISPR knockout studies of non-coding genomic elements, including long non-coding RNAs.

L1 integration in a transgenic mouse model

To study integration of the human LINE-1 retrotransposon (L1) in vivo, we developed a transgenic mouse model of L1 retrotransposition that displays de novo somatic L1 insertions at a high frequency, occasionally several insertions per mouse. We mapped 3' integration sites of 51 insertions by Thermal Asymmetric Interlaced PCR (TAIL-PCR). Analysis of integration locations revealed a broad genomic distribution with a modest preference for intergenic regions. We characterized the complete structures of 33 de novo retrotransposition events. Our results highlight the large number of highly truncated L1s, as over 52% (27/51) of total integrants were <1/3 the length of a full-length element. New integrants carry all structural characteristics typical of genomic L1s, including a number with inversions, deletions, and 5'-end microhomologies to the target DNA sequence. Notably, at least 13% (7/51) of all insertions contain a short stretch of extra nucleotides at their 5' end, which we postulate result from template-jumping by the L1-encoded reverse transcriptase. We propose a unified model of L1 integration that explains all of the characteristic features of L1 retrotransposition, such as 5' truncations, inversions, extra nucleotide additions, and 5' boundary and inversion point microhomologies.

Human inversions and their functional consequences

Polymorphic inversions are a type of structural variants that are difficult to analyze owing to their balanced nature and the location of breakpoints within complex repeated regions. So far, only a handful of inversions have been studied in detail in humans and current knowledge about their possible functional effects is still limited. However, inversions have been related to phenotypic changes and adaptation in multiple species. In this review, we summarize the evidences of the functional impact of inversions in the human genome. First, given that inversions have been shown to inhibit recombination in heterokaryotes, chromosomes displaying different orientation are expected to evolve independently and this may lead to distinct gene-expression patterns. Second, inversions have a role as disease-causing mutations both by directly affecting gene structure or regulation in different ways, and by predisposing to other secondary arrangements in the offspring of inversion carriers. Finally, several inversions show signals of being selected during human evolution. These findings illustrate the potential of inversions to have phenotypic consequences also in humans and emphasize the importance of their inclusion in genome-wide association studies.

Comparative geno-plasticity analysis of Mycoplasma bovis HB0801 (Chinese isolate)

Mycoplasma bovis pneumonia in cattle has been epidemic in China since 2008. To investigate M. bovis pathogenesis, we completed genome sequencing of strain HB0801 isolated from a lesioned bovine lung from Hubei, China. The genomic plasticity was determined by comparing HB0801 with M. bovis strain ATCC® 25523™/PG45 from cow mastitis milk, Chinese strain Hubei-1 from lesioned lung tissue, and 16 other Mycoplasmas species. Compared to PG45, the genome size of HB0801 was reduced by 11.7 kb. Furthermore, a large chromosome inversion (580 kb) was confirmed in all Chinese isolates including HB0801, HB1007, a strain from cow mastitis milk, and Hubei-1. In addition, the variable surface lipoproteins (vsp) gene cluster existed in HB0801, but contained less than half of the genes, and had poor identity to that in PG45, but they had conserved structures. Further inter-strain comparisons revealed other mechanisms of gene acquisition and loss in HB0801 that primarily involved insertion sequence (IS) elements, integrative conjugative element, restriction and modification systems, and some lipoproteins and transmembrane proteins. Subsequently, PG45 and HB0801 virulence in cattle was compared. Results indicated that both strains were pathogenic to cattle. The scores of gross pathological assessment for the control group, and the PG45- and HB0801-infected groups were 3, 13 and 9, respectively. Meanwhile the scores of lung lesion for these three groups were 36, 70, and 69, respectively. In addition, immunohistochemistry detection demonstrated that both strains were similarly distributed in lungs and lymph nodes. Although PG45 showed slightly higher virulence in calves than HB0801, there was no statistical difference between the strains (P>0.05). Compared to Hubei-1, a total of 122 SNP loci were disclosed in HB0801. In conclusion, although genomic plasticity was thought to be an evolutionary advantage, it did not apparently affect virulence of M. bovis strains in cattle.

Recurrent DNA inversion rearrangements in the human genome

Several lines of evidence suggest that reiterated sequences in the human genome are targets for nonallelic homologous recombination (NAHR), which facilitates genomic rearrangements. We have used a PCR-based approach to identify breakpoint regions of rearranged structures in the human genome. In particular, we have identified intrachromosomal identical repeats that are located in reverse orientation, which may lead to chromosomal inversions. A bioinformatic workflow pathway to select appropriate regions for analysis was developed. Three such regions overlapping with known human genes, located on chromosomes 3, 15, and 19, were analyzed. The relative proportion of wild-type to rearranged structures was determined in DNA samples from blood obtained from different, unrelated individuals. The results obtained indicate that recurrent genomic rearrangements occur at relatively high frequency in somatic cells. Interestingly, the rearrangements studied were significantly more abundant in adults than in newborn individuals, suggesting that such DNA rearrangements might start to appear during embryogenesis or fetal life and continue to accumulate after birth. The relevance of our results in regard to human genomic variation is discussed.

Selection for Chromosome Architecture in Bacteria

Bacterial chromosomes are immense polymers whose faithful replication and segregation are crucial to cell survival. The ability of proteins such as FtsK to move unidirectionally toward the replication terminus, and direct DNA translocation into the appropriate daughter cell during cell division, requires that bacterial genomes maintain an architecture for the orderly replication and segregation of chromosomes. We suggest that proteins that locate the replication terminus exploit strand-biased sequences that are overrepresented on one DNA strand, and that selection increases with decreased distance to the replication terminus. We report a generalized method for detecting these architecture imparting sequences (AIMS) and have identified AIMS in nearly all bacterial genomes. Their increased abundance on leading strands and decreased abundance on lagging strands toward replication termini are not the result of changes in mutational bias; rather, they reflect a gradient of long-term positive selection for AIMS. The maintenance of the pattern of AIMS across the genomes of related bacteria independent of their positions within individual genes suggests a well-conserved role in genome biology. The stable gradient of AIMS abundance from replication origin to terminus suggests that the replicore acts as a target of selection, where selection for chromosome architecture results in the maintenance of gene order and in the lack of high-frequency DNA inversion within replicores.

Patterns of diversity and linkage disequilibrium within the cosmopolitan inversion In(3R)Payne in Drosophila melanogaster are indicative of coadaptation

The cosmopolitan inversion In(3R)Payne in Drosophila melanogaster decreases in frequency with increasing distance from the equator on three continents, indicating it is subject to strong natural selection. We investigated patterns of genetic variation and linkage disequilibrium (LD) in 24 molecular markers located within and near In(3R)Payne to determine if different parts of the inversion responded to selection the same way. We found reduced variation in the markers we used compared to others distributed throughout the genome, consistent with the inversion having a relatively recent origin (<N(e) generations). LD between markers and In(3R)Payne varied significantly among markers within the inversion, with regions of high association interspersed by regions of low association. Several factors indicate that these patterns were not due to demographic factors such as admixture and bottlenecks associated with colonization, but instead reflected strong epistatic selection. Furthermore, we found that nonadjacent regions with high association to the inversion contained markers with the strongest clinal patterns in allele frequency; in most cases, the level of clinal variation was beyond what could be explained by hitchhiking with In(3R)Payne, indicating that genes within these regions are targets of selection. Our results provide some support for the hypothesis that inversions persist in natural populations because they hold together favorable combinations of alleles that act together to facilitate adaptive shifts.

Genomic architecture of human 17q21 linked to frontotemporal dementia uncovers a highly homologous family of low-copy repeats in the tau region

Familial frontotemporal dementia (FTD), characterized by tau-negative, ubiquitin-positive inclusions at autopsy, is linked to a chromosomal region at 17q21 (FTDU-17), encompassing the gene encoding the microtubule associated protein tau, MAPT. Mutations in MAPT were previously identified in familial FTD with parkinsonism (FTDP-17); however, in FTDU-17 patients, no pathogenic mutations were found in exonic regions consistent with the lack of tauopathy in FTDU-17 brains. Here, we excluded mutations in MAPT by genomic sequencing of 138.5 kb in FTDU-17 patients. Next, to facilitate the identification of the actual underlying genetic defect, we assembled the 6.5 Mb FTDU-17 sequence. Annotation demonstrated that MAPT is surrounded by three highly homologous low-copy repeats (LCRs) in a region of 1.7 Mb. Using evolutionary studies, short tandem repeat-based linkage disequilibrium (LD) and macro-restriction mapping, we demonstrated that these LCRs are at the basis of a series of rearrangements in the MAPT genomic region. One is an inversion that occurred 3 million years ago and resulted in a common polymorphism in humans to date. This inversion plus flanking LCRs spanned approximately 1.3 Mb and was shown to underlie the extended LD and haplotypes H1 and H2 across MAPT. However, in the FTDU-17 families, we ascertained segregation analysis precluding a relationship between the FTDU-17 and the H1/H2 inversion. The presence of multiple homologous LCRs in the region predicts that other potentially more complex genomic rearrangements might be underlying FTDU-17.