Three linked enzyme loci in fishes: implications in the evolution of vertebrate chromosomes.

A three-point linkage group comprised of loci coding for adenosine deaminase (ADA), glucose-6-phosphate dehydrogenase (G6PDH), and 6-phospho-gluconate dehydrogenase (6PGD) is described in fish of the genus Xiphophorus (Poeciliidae). The alleles at loci in this group were shown to assort independently from the alleles at three other loci--isocitrate dehydrogenase 1 and 2, and glyceraldehyde-3-phosphate dehydrogenase 1. Alleles at the latter three loci also assort independently from each other. Data were obtained by observing the segregation of electrophoretically variant alleles in reciprocal backcross hybrids derived from crosses between either X. helleri guentheri or X. h. strigatus and X. maculatus. The linkage component of chi2 was significant (less than 0.01) in all crosses, indicating that the linkage group is conserved in all populations of both species of Xiphophorus examined. While data from X. h. guentheri backcrosses indicate the linkage relationship ADA--6%--G6PDH--24%--6PGD, and ADA--29%--6PGD (30% when corrected for double crossovers), data from backcrosses involving strigatus, while supporting the same gene order, yielded significantly different recombination frequencies. The likelihood of the difference being due to an inversion could not be separated from the possibility of a sex effect on recombination in the present data. The linkage of 6PGD and G6PDH has been shown to exist in species of at least three classes of vertebrates, indicating the possibility of evolutionary conservation of this linkage.

Clonal Evolution and Blast Crisis Correlate with Enhanced Proteolytic Activity of Separase in BCR-ABL b3a2 Fusion Type CML under Imatinib Therapy.

Unbalanced (major route) additional cytogenetic aberrations (ACA) at diagnosis of chronic myeloid leukemia (CML) indicate an increased risk of progression and shorter survival. Moreover, newly arising ACA under imatinib treatment and clonal evolution are considered features of acceleration and define failure of therapy according to the European LeukemiaNet (ELN) recommendations. On the basis of 1151 Philadelphia chromosome positive chronic phase patients of the randomized CML-study IV, we examined the incidence of newly arising ACA under imatinib treatment with regard to the p210BCR-ABL breakpoint variants b2a2 and b3a2. We found a preferential acquisition of unbalanced ACA in patients with b3a2 vs. b2a2 fusion type (ratio: 6.3 vs. 1.6, p = 0.0246) concurring with a faster progress to blast crisis for b3a2 patients (p = 0.0124). ESPL1/Separase, a cysteine endopeptidase, is a key player in chromosomal segregation during mitosis. Separase overexpression and/or hyperactivity has been reported from a wide range of cancers and cause defective mitotic spindles, chromosome missegregation and aneuploidy. We investigated the influence of p210BCR-ABL breakpoint variants and imatinib treatment on expression and proteolytic activity of Separase as measured with a specific fluorogenic assay on CML cell lines (b2a2: KCL-22, BV-173; b3a2: K562, LAMA-84). Despite a drop in Separase protein levels an up to 5.4-fold increase of Separase activity under imatinib treatment was observed exclusively in b3a2 but not in b2a2 cell lines. Mimicking the influence of imatinib on BV-173 and LAMA-84 cells by ESPL1 silencing stimulated Separase proteolytic activity in both b3a2 and b2a2 cell lines. Our data suggest the existence of a fusion type-related feedback mechanism that posttranslationally stimulates Separase proteolytic activity after therapy-induced decreases in Separase protein levels. This could render b3a2 CML cells more prone to aneuploidy and clonal evolution than b2a2 progenitors and may therefore explain the cytogenetic results of CML patients.

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Fluorescence in situ hybridization analysis of keratinocyte growth factor gene amplification and dispersion in evolution of great apes and humans

Keratinocyte growth factor (KGF) is a member of the fibroblast growth factor family. Portions of the gene encoding KGF were amplified during primate evolution and are present in multiple nonprocessed copies in the human genome. Nucleotide analysis of a representative sampling of these KGF-like sequences indicated that they were at least 95% identical to corresponding regions of the KGF gene. To localize these sequences to specific chromosomal sites in human and higher primates, we used fluorescence in situ hybridization. In human, using a cosmid probe encoding KGF exon 1, we assigned the location of the KGF gene to chromosome 15q15–21.1. In addition, copies of KGF-like sequences hybridizing only with a cosmid probe encoding exons 2 and 3 were localized to dispersed sites on chromosome 2q21, 9p11, 9q12–13, 18p11, 18q11, 21q11, and 21q21.1. The distribution of KGF-like sequences suggests a role for alphoid DNA in their amplification and dispersion. In chimpanzee, KGF-like sequences were observed at five chromosomal sites, which were each homologous to sites in human, while in gorilla, a subset of four of these homologous sites was identified; in orangutan two sites were identified, while gibbon exhibited only a single site. The chromosomal localization of KGF sequences in human and great ape genomes indicates that amplification and dispersion occurred in multiple discrete steps, with initial KGF gene duplication and dispersion taking place in gibbon and involving loci corresponding to human chromosomes 15 and 21. These findings support the concept of a closer evolutionary relationship of human and chimpanzee and a possible selective pressure for such dispersion during the evolution of higher primates.

Complete sequencing of the Fugu WAGR region from WT1 to PAX6: Dramatic compaction and conservation of synteny with human chromosome 11p13

The pufferfish Fugu rubripes has a genome ≈7.5 times smaller than that of mammals but with a similar number of genes. Although conserved synteny has been demonstrated between pufferfish and mammals across some regions of the genome, there is some controversy as to what extent Fugu will be a useful model for the human genome, e.g., [Gilley, J., Armes, N. & Fried, M. (1997) Nature (London) 385, 305–306]. We report extensive conservation of synteny between a 1.5-Mb region of human chromosome 11 and <100 kb of the Fugu genome in three overlapping cosmids. Our findings support the idea that the majority of DNA in the region of human chromosome 11p13 is intergenic. Comparative analysis of three unrelated genes with quite different roles, WT1, RCN1, and PAX6, has revealed differences in their structural evolution. Whereas the human WT1 gene can generate 16 protein isoforms via a combination of alternative splicing, RNA editing, and alternative start site usage, our data predict that Fugu WT1 is capable of generating only two isoforms. This raises the question of the extent to which the evolution of WT1 isoforms is related to the evolution of the mammalian genitourinary system. In addition, this region of the Fugu genome shows a much greater overall compaction than usual but with significant noncoding homology observed at the PAX6 locus, implying that comparative genomics has identified regulatory elements associated with this gene.

Evolution of the Friedreich’s ataxia trinucleotide repeat expansion: Founder effect and premutations

Friedreich’s ataxia, the most frequent inherited ataxia, is caused, in the vast majority of cases, by large GAA repeat expansions in the first intron of the frataxin gene. The normal sequence corresponds to a moderately polymorphic trinucleotide repeat with bimodal size distribution. Small normal alleles have approximately eight to nine repeats whereas a more heterogeneous mode of large normal alleles ranges from 16 to 34 GAA. The latter class accounts for ≈17% of normal alleles. To identify the origin of the expansion mutation, we analyzed linkage disequilibrium between expansion mutations or normal alleles and a haplotype of five polymorphic markers within or close to the frataxin gene; 51% of the expansions were associated with a single haplotype, and the other expansions were associated with haplotypes that could be related to the major one by mutation at a polymorphic marker or by ancient recombination. Of interest, the major haplotype associated with expansion is also the major haplotype associated with the larger alleles in the normal size range and was almost never found associated with the smaller normal alleles. The results indicate that most if not all large normal alleles derive from a single founder chromosome and that they represent a reservoir for larger expansion events, possibly through “premutation” intermediates. Indeed, we found two such alleles (42 and 60 GAA) that underwent cataclysmic expansion to pathological range in a single generation. This stepwise evolution to large trinucleotide expansions already was suggested for myotonic dystrophy and fragile X syndrome and may relate to a common mutational mechanism, despite sequence motif differences.

Novel gene conversion between X-Y homologues located in the nonrecombining region of the Y chromosome in Felidae (Mammalia)

Genes located on the mammalian Y chromosome outside of the pseudoautosomal region do not recombine with those on the X and are predicted to either undergo selection for male function or gradually degenerate because of an accumulation of deleterious mutations. Here, phylogenetic analyses of X-Y homologues, Zfx and Zfy, among 26 felid species indicate two ancestral episodes of directed genetic exchange (ectopic gene conversion) from X to Y: once during the evolution of pallas cat and once in a common predecessor of ocelot lineage species. Replacement of the more rapidly evolving Y homologue with the evolutionarily constrained X copy may represent a mechanism for adaptive editing of functional genes on the nonrecombining region of the mammalian Y chromosome.

Molecular cytogenetic dissection of human chromosomes 3 and 21 evolution

Chromosome painting in placental mammalians illustrates that genome evolution is marked by chromosomal synteny conservation and that the association of chromosomes 3 and 21 may be the largest widely conserved syntenic block known for mammals. We studied intrachromosomal rearrangements of the syntenic block 3/21 by using probes derived from chromosomal subregions with a resolution of up to 10–15 Mbp. We demonstrate that the rearrangements visualized by chromosome painting, mostly translocations, are only a fraction of the actual chromosomal changes that have occurred during evolution. The ancestral segment order for both primates and carnivores is still found in some species in both orders. From the ancestral primate/carnivore condition an inversion is needed to derive the pig homolog, and a fission of chromosome 21 and a pericentric inversion is needed to derive the Bornean orangutan condition. Two overlapping inversions in the chromosome 3 homolog then would lead to the chromosome form found in humans and African apes. This reconstruction of the origin of human chromosome 3 contrasts with the generally accepted scenario derived from chromosome banding in which it was proposed that only one pericentric inversion was needed. From the ancestral form for Old World primates (now found in the Bornean orangutan) a pericentric inversion and centromere shift leads to the chromosome ancestral for all Old World monkeys. Intrachromosomal rearrangements, as shown here, make up a set of potentially plentiful and informative markers that can be used for phylogenetic reconstruction and a more refined comparative mapping of the genome.

Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics

The evolutionary relationships of 46 Shigella strains representing each of the serotypes belonging to the four traditional Shigella species (subgroups), Dysenteriae, Flexneri, Boydii, and Sonnei, were determined by sequencing of eight housekeeping genes in four regions of the chromosome. Analysis revealed a very similar evolutionary pattern for each region. Three clusters of strains were identified, each including strains from different subgroups. Cluster 1 contains the majority of Boydii and Dysenteriae strains (B1–4, B6, B8, B10, B14, and B18; and D3–7, D9, and D11–13) plus Flexneri 6 and 6A. Cluster 2 contains seven Boydii strains (B5, B7, B9, B11, B15, B16, and B17) and Dysenteriae 2. Cluster 3 contains one Boydii strain (B12) and the Flexneri serotypes 1–5 strains. Sonnei and three Dysenteriae strains (D1, D8, and D10) are outside of the three main clusters but, nonetheless, are clearly within Escherichia coli. Boydii 13 was found to be distantly related to E. coli. Shigella strains, like the other pathogenic forms of E. coli, do not have a single evolutionary origin, indicating convergent evolution of Shigella phenotypic properties. We estimate the three main Shigella clusters to have evolved within the last 35,000 to 270,000 years, suggesting that shigellosis was one of the early infectious diseases of humans.

Spectral karyotyping suggests additional subsets of colorectal cancers characterized by pattern of chromosome rearrangement

The abundant chromosome abnormalities in most carcinomas are probably a reflection of genomic instability present in the tumor, so the pattern and variability of chromosome abnormalities will reflect the mechanism of instability combined with the effects of selection. Chromosome rearrangement was investigated in 17 colorectal carcinoma-derived cell lines. Comparative genomic hybridization showed that the chromosome changes were representative of those found in primary tumors. Spectral karyotyping (SKY) showed that translocations were very varied and mostly unbalanced, with no translocation occurring in more than three lines. At least three karyotype patterns could be distinguished. Some lines had few chromosome abnormalities: they all showed microsatellite instability, the replication error (RER)+ phenotype. Most lines had many chromosome abnormalities: at least seven showed a surprisingly consistent pattern, characterized by multiple unbalanced translocations and intermetaphase variation, with chromosome numbers around triploid, 6–16 structural aberrations, and similarities in gains and losses. Almost all of these were RER−, but one, LS411, was RER+. The line HCA7 showed a novel pattern, suggesting a third kind of genomic instability: multiple reciprocal translocations, with little numerical change or variability. This line was also RER+. The coexistence in one tumor of two kinds of genomic instability is to be expected if the underlying defects are selected for in tumor evolution.

Incipient speciation by sexual isolation in Drosophila: Concurrent evolution at multiple loci

Drosophila melanogaster from Zimbabwe and nearby regions shows strong but asymmetric sexual isolation from its cosmopolitan counterparts. By creating stable chromosome-substitution lines, earlier studies were able to show that the two major autosomes have very large effects on both male mating success and female mating preference. In this study, we genetically dissect this sexual isolation by recombination analysis between a whole-chromosome substitution line (which carries a Zimbabwe-derived third chromosome) and a strain with seven visible markers on that chromosome. Four loci are responsible for male mating success and three others are found to control female mating preference. Because male and female traits are not closely linked, their strong association among isofemale lines is most likely a reflection of sexual selection in nature. The results suggest that a large number of behavioral loci may evolve concurrently in the incipient stage of speciation before other aspects of reproductive isolation (such as hybrid sterility) have become evident. The results shed light on the population genetic processes underlying the formation of nascent species, as well as modes of speciation.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Centromere identification and inactivation on a neo-Y chromosome fusion in threespine stickleback fish (Gasterosteus aculeatus)

Thesis (Ph.D.)--University of Washington, 2016-05

Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA Microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design microarray probes with high specificities and sensitivities. Highly specific probes correspond to probes having unique DNA sequences; whereas highly sensitive probes correspond to those with melting temperature within a desired range and having no secondary structure. The selection of these probes from a set of functional DNA sequences (exons) constitutes a computationally expensive discrete non-linear search problem. We delegate the search task to a simple yet effective Evolution Strategy algorithm. The computational efficiency is also greatly improved by making use of an available bioinformatics tool.

Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function

The mammalian transcriptome contains many nonprotein-coding RNAs (ncRNAs), but most of these are of unclear significance and lack strong sequence conservation, prompting suggestions that they might be non-functional. However, certain long functional ncRNAs such as Air and Xist are also poorly conserved. In this article, we systematically analyzed the conservation of several groups of functional ncRNAs, including miRNAs, snoRNAs and longer ncRNAs whose function has been either documented or confidently predicted. As expected, miRNAs and snoRNAs were highly conserved. By contrast, the longer functional non-micro, non-sno ncRNAs were much less conserved with many displaying rapid sequence evolution. Our findings suggest that longer ncRNAs are under the influence of different evolutionary constraints and that the lack of conservation displayed by the thousands of candidate ncRNAs does not necessarily signify an absence of function.