Comparative DNA sequence and methylation analyses of orthologous genes in humans and non-human primates: Genetic and epigenetic footprints of evolution

Welche genetische Unterschiede machen uns verschieden von unseren nächsten Verwandten, den Schimpansen, und andererseits so ähnlich zu den Schimpansen? Was wir untersuchen und auch verstehen wollen, ist die komplexe Beziehung zwischen den multiplen genetischen und epigenetischen Unterschieden, deren Interaktion mit diversen Umwelt- und Kulturfaktoren in den beobachteten phänotypischen Unterschieden resultieren. Um aufzuklären, ob chromosomale Rearrangements zur Divergenz zwischen Mensch und Schimpanse beigetragen haben und welche selektiven Kräfte ihre Evolution geprägt haben, habe ich die kodierenden Sequenzen von 2 Mb umfassenden, die perizentrischen Inversionsbruchpunkte flankierenden Regionen auf den Chromosomen 1, 4, 5, 9, 12, 17 und 18 untersucht. Als Kontrolle dienten dabei 4 Mb umfassende kollineare Regionen auf den rearrangierten Chromosomen, welche mindestens 10 Mb von den Bruchpunktregionen entfernt lagen. Dabei konnte ich in den Bruchpunkten flankierenden Regionen im Vergleich zu den Kontrollregionen keine höhere Proteinevolutionsrate feststellen. Meine Ergebnisse unterstützen nicht die chromosomale Speziationshypothese für Mensch und Schimpanse, da der Anteil der positiv selektierten Gene (5,1% in den Bruchpunkten flankierenden Regionen und 7% in den Kontrollregionen) in beiden Regionen ähnlich war. Durch den Vergleich der Anzahl der positiv und negativ selektierten Gene per Chromosom konnte ich feststellen, dass Chromosom 9 die meisten und Chromosom 5 die wenigsten positiv selektierten Gene in den Bruchpunkt flankierenden Regionen und Kontrollregionen enthalten. Die Anzahl der negativ selektierten Gene (68) war dabei viel höher als die Anzahl der positiv selektierten Gene (17). Eine bioinformatische Analyse von publizierten Microarray-Expressionsdaten (Affymetrix Chip U95 und U133v2) ergab 31 Gene, die zwischen Mensch und Schimpanse differentiell exprimiert sind. Durch Untersuchung des dN/dS-Verhältnisses dieser 31 Gene konnte ich 7 Gene als negativ selektiert und nur 1 Gen als positiv selektiert identifizieren. Dieser Befund steht im Einklang mit dem Konzept, dass Genexpressionslevel unter stabilisierender Selektion evolvieren. Die meisten positiv selektierten Gene spielen überdies eine Rolle bei der Fortpflanzung. Viele dieser Speziesunterschiede resultieren eher aus Änderungen in der Genregulation als aus strukturellen Änderungen der Genprodukte. Man nimmt an, dass die meisten Unterschiede in der Genregulation sich auf transkriptioneller Ebene manifestieren. Im Rahmen dieser Arbeit wurden die Unterschiede in der DNA-Methylierung zwischen Mensch und Schimpanse untersucht. Dazu wurden die Methylierungsmuster der Promotor-CpG-Inseln von 12 Genen im Cortex von Menschen und Schimpansen mittels klassischer Bisulfit-Sequenzierung und Bisulfit-Pyrosequenzierung analysiert. Die Kandidatengene wurden wegen ihrer differentiellen Expressionsmuster zwischen Mensch und Schimpanse sowie wegen Ihrer Assoziation mit menschlichen Krankheiten oder dem genomischen Imprinting ausgewählt. Mit Ausnahme einiger individueller Positionen zeigte die Mehrzahl der analysierten Gene keine hohe intra- oder interspezifische Variation der DNA-Methylierung zwischen den beiden Spezies. Nur bei einem Gen, CCRK, waren deutliche intraspezifische und interspezifische Unterschiede im Grad der DNA-Methylierung festzustellen. Die differentiell methylierten CpG-Positionen lagen innerhalb eines repetitiven Alu-Sg1-Elements. Die Untersuchung des CCRK-Gens liefert eine umfassende Analyse der intra- und interspezifischen Variabilität der DNA-Methylierung einer Alu-Insertion in eine regulatorische Region. Die beobachteten Speziesunterschiede deuten darauf hin, dass die Methylierungsmuster des CCRK-Gens wahrscheinlich in Adaption an spezifische Anforderungen zur Feinabstimmung der CCRK-Regulation unter positiver Selektion evolvieren. Der Promotor des CCRK-Gens ist anfällig für epigenetische Modifikationen durch DNA-Methylierung, welche zu komplexen Transkriptionsmustern führen können. Durch ihre genomische Mobilität, ihren hohen CpG-Anteil und ihren Einfluss auf die Genexpression sind Alu-Insertionen exzellente Kandidaten für die Förderung von Veränderungen während der Entwicklungsregulation von Primatengenen. Der Vergleich der intra- und interspezifischen Methylierung von spezifischen Alu-Insertionen in anderen Genen und Geweben stellt eine erfolgversprechende Strategie dar.

A linkage and cross-sectional study of hypertension and obesity using a poly(A) Alu-repeat polymorphism at the glucagon receptor gene locus (17q25)

1. Previous glucagon receptor gene (GCGR) studies have shown a Gly40Ser mutation to be more prevalent in essential hypertension and to affect glucagon binding affinity to its receptor. An Alu-repeat poly(A) polymorphism colocalized to GCGR was used in the present study to test for association and linkage in hypertension as well as association in obesity development. 2. Using a cross-sectional approach, 85 hypertensives and 95 normotensives were genotyped using polymerase chain reaction primers flanking the Alu-repeat. Both hypertensive and normotensive populations were subdivided into lean and obese categories based on body mass index (BMI) to determine involvement of this variant in obesity. For the linkage study, 89 Australian Caucasian hypertension affected sibships (174 sibpairs) were genotyped and the results were analysed using GENE-HUNTER, Mapmaker Sibs, ERPA and SPLINK (all freely available from http://linlkage.rockefeller. edu/soft/list.html). 3. Cross-sectional results for both hypertension and obesity were analysed using Chi-squared and Monte Carlo analyses. Results did not show an association of this variant with either hypertension (χ2 = 6.9, P = 0.14; Monte Carlo χ2 = 7.0, P = 0.11; n = 5000) or obesity (χ2 = 3.3, P = 0.35; Monte Carlo χ2 = 3.26, P = 0.34; n = 5000). In addition, results from the linkage study using hypertensive sib-pairs did not indicate linkage of the poly(A) repent with hypertension. Hence, results did not indicate a role far the Alu-repeat in either hypertension or obesity. However, as the heterozygosity of this poly(A) repeat is low (35%), a larger number of hypertensive sib-pairs may be required to draw definitive conclusions.

Mutations in ADAR1 cause Aicardi-Goutieres syndrome associated with a type I interferon signature

Adenosine deaminases acting on RNA (ADARs) catalyze the hydrolytic deamination of adenosine to inosine in double-stranded RNA (dsRNA) and thereby potentially alter the information content and structure of cellular RNAs. Notably, although the overwhelming majority of such editing events occur in transcripts derived from Alu repeat elements, the biological function of non-coding RNA editing remains uncertain. Here, we show that mutations in ADAR1 (also known as ADAR) cause the autoimmune disorder Aicardi-Goutieres syndrome (AGS). As in Adar1-null mice, the human disease state is associated with upregulation of interferon-stimulated genes, indicating a possible role for ADAR1 as a suppressor of type I interferon signaling. Considering recent insights derived from the study of other AGS-related proteins, we speculate that ADAR1 may limit the cytoplasmic accumulation of the dsRNA generated from genomic repetitive elements.

Transcriptional regulation and functional analysis of the Wilms' tumor gene WT1

The Wilms' tumor gene, WT1, encodes a zinc finger transcription factor which functions as a tumor suppressor. Defects in the WT1 gene can result in the development of nephroblastoma. WT1 is expressed during development, primarily in the metanephric kidney, the mesothelial lining of the abdomen and thorax, and the developing gonads. WT1 expression is tightly regulated and is essential for renal development. The WT1 gene encodes a protein with a proline-rich N-terminus which functions as a transcriptional repressor and C-terminus contains 4 zinc fingers that mediate DNA binding. WT1 represses transcription from a number of growth factors and growth factor receptors. WT1 mRNA undergoes alternative splicing at two sites, resulting in 4 mRNA species and polypeptide products. Exon 5, encoding 17 amino acids is alternatively spliced, and is located between the transcriptional repression domain and the DNA binding domain. The second alternative splice is the terminal 9 nucleotides of zinc finger 3, encoding the tripeptide Lys-Thr-Ser (KTS). The presence or absence of KTS within the zinc fingers of WT1 alters DNA binding.^ I have investigated transcriptional regulation of WT1, characterizing two means of repressing WT1 transcription. I have cloned a transcriptional silencer of the WT1 promoter which is located in the third intron of the WT1 gene. The silencer is 460 bp in length and contains an Alu repeat. The silencer functions in cells of non-renal origin.^ I have found that WT1 protein can autoregulate the WT1 promoter. Using the autoregulation of the WT1 promoter as a functional assay, I have defined differential consensus DNA binding motifs of WT1 isoforms lacking and containing the KTS tripeptide insertion. With these refined consensus DNA binding motifs, I have identified two additional targets of WT1 transcriptional repression, the proto-oncogenes bcl-2 and c-myc.^ I have investigated the ability of the alternatively spliced exon 5 to influence cell growth. In cell proliferation assays, isoforms of WT1 lacking exon 5 repress cell growth. WT1 isoforms containing exon 5 fail to repress cell growth to the same extent, but alter the morphology of the cells. These experiments demonstrate that the alternative splice isoforms of WT1 have differential effects on the function of WT1. These findings suggest a role for the alternative splicing of WT1 in metanephric development. ^

Functional copies of a human gene can be directly isolated by transformation-associated recombination cloning with a small 3′ end target sequence

Unique, small sequences (sequence tag sites) have been identified at the 3′ ends of most human genes that serve as landmarks in genome mapping. We investigated whether a single copy gene could be isolated directly from total human DNA by transformation-associated recombination (TAR) cloning in yeast using a short, 3′ unique target. A TAR cloning vector was constructed that, when linearized, contained a small amount (381 bp) of 3′ hypoxanthine phosphoribosyltransferase (HPRT) sequence at one end and an 189-bp Alu repeat at the other end. Transformation with this vector along with human DNA led to selective isolations of the entire HPRT gene as yeast artificial chromosomes (YACs) that extended from the 3′ end sequence to various Alu positions as much as 600 kb upstream. These YACs were retrofitted with a NeoR and a bacterial artificial chromosome (BAC) sequence to transfer the YACs to bacteria and subsequently the BACs to mouse cells by using a Neo selection. Most of the HPRT isolates were functional, demonstrating that TAR cloning retains the functional integrity of the isolated material. Thus, this modified version of TAR cloning, which we refer to as radial TAR cloning, can be used to isolate large segments of the human genome accurately and directly with only a small amount of sequence information.

"Pseudo-Beijing": evidence for convergent evolution in the direct repeat region of Mycobacterium tuberculosis

Background Mycobacterium tuberculosis has a global population structure consisting of six main phylogenetic lineages associated with specific geographic regions and human populations. One particular M. tuberculosis genotype known as “Beijing” has repeatedly been associated with drug resistance and has been emerging in some parts of the world. “Beijing” strains are traditionally defined based on a characteristic spoligotyping pattern. We used three alternative genotyping techniques to revisit the phylogenetic classification of M. tuberculosis complex (MTBC) strains exhibiting the typical “Beijing” spoligotyping pattern. Methods and Findings MTBC strains were obtained from an ongoing molecular epidemiological study in Switzerland and Nepal. MTBC genotyping was performed based on SNPs, genomic deletions, and 24-loci MIRU-VNTR. We identified three MTBC strains from patients originating from Tibet, Portugal and Nepal which exhibited a spoligotyping patterns identical to the classical Beijing signature. However, based on three alternative molecular markers, these strains were assigned to Lineage 3 (also known as Delhi/CAS) rather than to Lineage 2 (also known as East-Asian lineage). Sequencing of the RD207 in one of these strains showed that the deletion responsible for this “Pseudo-Beijing” spoligotype was about 1,000 base pairs smaller than the usual deletion of RD207 in classical “Beijing” strains, which is consistent with an evolutionarily independent deletion event in the direct repeat (DR) region of MTBC. Conclusions We provide an example of convergent evolution in the DR locus of MTBC, and highlight the limitation of using spoligotypes for strain classification. Our results indicate that a proportion of “Beijing” strains may have been misclassified in the past. Markers that are more phylogenetically robust should be used when exploring strain-specific differences in experimental or clinical phenotypes.

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.

Multiple evolutionary rate classes in animal genome evolution

The proportion of functional sequence in the human genome is currently a subject of debate. The most widely accepted figure is that approximately 5% is under purifying selection. In Drosophila, estimates are an order of magnitude higher, though this corresponds to a similar quantity of sequence. These estimates depend on the difference between the distribution of genomewide evolutionary rates and that observed in a subset of sequences presumed to be neutrally evolving. Motivated by the widening gap between these estimates and experimental evidence of genome function, especially in mammals, we developed a sensitive technique for evaluating such distributions and found that they are much more complex than previously apparent. We found strong evidence for at least nine well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least seven classes in an alignment of four mammals, including human. We also identified at least three rate classes in human ancestral repeats. By positing that the largest of these ancestral repeat classes is neutrally evolving, we estimate that the proportion of nonneutrally evolving sequence is 30% of human ancestral repeats and 45% of the aligned portion of the genome. However, we also question whether any of the classes represent neutrally evolving sequences and argue that a plausible alternative is that they reflect variable structure-function constraints operating throughout the genomes of complex organisms.

High-efficiency silencing of a β-glucuronidase gene in rice is correlated with repetitive transgene structure but is independent of DNA methylation

Two transgenic callus lines of rice, stably expressing a β-glucuronidase (GUS) gene, were supertransformed with a set of constructs designed to silence the resident GUS gene. An inverted-repeat (i/r) GUS construct, designed to produce mRNA with self-complementarity, was much more effective than simple sense and antisense constructs at inducing silencing. Supertransforming rice calluses with a direct-repeat (d/r) construct, although not as effective as those with the i/r construct, was also substantially more effective in silencing the resident GUS gene than the simple sense and antisense constructs. DNA hybridisation analyses revealed that every callus line supertransformed with either simple sense or antisense constructs, and subsequently showing GUS silencing, had the silence-inducing transgenes integrated into the plant genome in inverted-repeat configurations. The silenced lines containing i/r and d/r constructs did not necessarily have inverted-repeat T-DNA insertions. There was significant methylation of the GUS sequences in most of the silenced lines but not in the unsilenced lines. However, demethylation treatment of silenced lines with 5-azacytidine did not reverse the post-transcriptional gene silencing (PTGS) of GUS. Whereas the levels of RNA specific to the resident GUS gene were uniformly low in the silenced lines, RNA specific to the inducer transgenes accumulated to a substantial level, and the majority of the i/r RNA was unpolyadenylated. Altogether, these results suggest that both sense- and antisense-mediated gene suppression share a similar molecular basis, that unpolyadenylated RNA plays an important role in PTGS, and that methylation is not essential for PTGS.

Targeting DNA repeat sequences with Py-Im polyamides

Hairpin pyrrole-imdazole polyamides are cell-permeable, sequence-programmable oligomers that bind in the minor groove of DNA. This thesis describes studies of Py-Im polyamides targeted to biologically important DNA repeat sequences for the purpose of modulating disease states. Design of a hairpin polyamide that binds the CG dyad, a site of DNA methylation that can become dysregulated in cancer, is described. We report the synthesis of a DNA methylation antagonist, its sequence specificity and affinity informed by Bind-n-Seq and iteratively designed, which improves inhibitory activity in a cell-free assay by 1000-fold to low nanomolar IC50. Additionally, a hairpin polyamide targeted to the telomeric sequence is found to trigger a slow necrotic-type cell death with the release of inflammatory molecules in a model of B cell lymphoma. The effects of the polyamide are unique in this class of oligomers; its effects are characterized and a functional assay of phagocytosis by macrophages is described. Additionally, hairpin polyamides targeted to pathologically expanded CTG•CAG triplet repeat DNA sequences, the molecular cause of myotonic dystrophy type 1, are synthesized and assessed for toxicity. Lastly, ChIP-seq of Hypoxia-Inducible Factor is performed under hypoxia-induced conditions. The study results show that ChIP-seq can be employed to understand the genome-wide perturbation of Hypoxia-Inducible Factor occupancy by a Py-Im polyamide.

Molecular evolution of a primate-specific microRNA family

Lineage-specific microRNA (miRNA) families may contribute to developmental novelties during evolution. However, little is known about the origin and evolution of new miRNA families. We report evidence of an Alu-mediated rapid expansion of miRNA genes in a

High-resolution organellar genome analysis of Triticum and Aegilops sheds new light on cytoplasm evolution in wheat

We have utilised polymorphic chloroplast microsatellites to analyse cytoplasmic relationships between accessions in the genera Triticum and Aegilops. Sequencing of PCR products revealed point mutations and insertions/deletions in addition to the standard repeat length expansion/contraction which most likely represent ancient synapomorphies. Phylogenetic analyses revealed three distinct groups of accessions. One of these contained all the non-Aegilops speltoides S-type cytoplasm species, another comprised almost exclusively A, C, D, M, N, T and U cytoplasm-type accessions and the third contained the polyploid Triticum species and all the Ae. speltoides accessions, further confirming that Ae. speltoides or a closely related but now extinct species was the original B-genome donor of cultivated polyploid wheat. Successive decreases in levels of genetic diversity due to domestication were also observed. Finally, we highlight the importance of elucidating longer-term evolutionary processes operating at microsatellite repeat loci.

The Fasciola hepatica genome: gene duplication and polymorphism reveals adaptation to the host environment and the capacity for rapid evolution

BACKGROUND: The liver fluke Fasciola hepatica is a major pathogen of livestock worldwide, causing huge economic losses to agriculture, as well as 2.4 million human infections annually.

RESULTS: Here we provide a draft genome for F. hepatica, which we find to be among the largest known pathogen genomes at 1.3 Gb. This size cannot be explained by genome duplication or expansion of a single repeat element, and remains a paradox given the burden it may impose on egg production necessary to transmit infection. Despite the potential for inbreeding by facultative self-fertilisation, substantial levels of polymorphism were found, which highlights the evolutionary potential for rapid adaptation to changes in host availability, climate change or to drug or vaccine interventions. Non-synonymous polymorphisms were elevated in genes shared with parasitic taxa, which may be particularly relevant for the ability of the parasite to adapt to a broad range of definitive mammalian and intermediate molluscan hosts. Large-scale transcriptional changes, particularly within expanded protease and tubulin families, were found as the parasite migrated from the gut, across the peritoneum and through the liver to mature in the bile ducts. We identify novel members of anti-oxidant and detoxification pathways and defined their differential expression through infection, which may explain the stage-specific efficacy of different anthelmintic drugs.

CONCLUSIONS: The genome analysis described here provides new insights into the evolution of this important pathogen, its adaptation to the host environment and external selection pressures. This analysis also provides a platform for research into novel drugs and vaccines.