Characterizing Genetic Drivers of Lymphoma through High-Throughput Sequencing

The advent of next-generation sequencing, now nearing a decade in age, has enabled, among other capabilities, measurement of genome-wide sequence features at unprecedented scale and resolution.

In this dissertation, I describe work to understand the genetic underpinnings of non-Hodgkin’s lymphoma through exploration of the epigenetics of its cell of origin, initial characterization and interpretation of driver mutations, and finally, a larger-scale, population-level study that incorporates mutation interpretation with clinical outcome.

In the first research chapter, I describe genomic characteristics of lymphomas through the lens of their cells of origin. Just as many other cancers, such as breast cancer or lung cancer, are categorized based on their cell of origin, lymphoma subtypes can be examined through the context of their normal B Cells of origin, Naïve, Germinal Center, and post-Germinal Center. By applying integrative analysis of the epigenetics of normal B Cells of origin through chromatin-immunoprecipitation sequencing, we find that differences in normal B Cell subtypes are reflected in the mutational landscapes of the cancers that arise from them, namely Mantle Cell, Burkitt, and Diffuse Large B-Cell Lymphoma.

In the next research chapter, I describe our first endeavor into understanding the genetic heterogeneity of Diffuse Large B Cell Lymphoma, the most common form of non-Hodgkin’s lymphoma, which affects 100,000 patients in the world. Through whole-genome sequencing of 1 case as well as whole-exome sequencing of 94 cases, we characterize the most recurrent genetic features of DLBCL and lay the groundwork for a larger study.

In the last research chapter, I describe work to characterize and interpret the whole exomes of 1001 cases of DLBCL in the largest single-cancer study to date. This highly-powered study enabled sub-gene, gene-level, and gene-network level understanding of driver mutations within DLBCL. Moreover, matched genomic and clinical data enabled the connection of these driver mutations to clinical features such as treatment response or overall survival. As sequencing costs continue to drop, whole-exome sequencing will become a routine clinical assay, and another diagnostic dimension in addition to existing methods such as histology. However, to unlock the full utility of sequencing data, we must be able to interpret it. This study undertakes a first step in developing the understanding necessary to uncover the genomic signals of DLBCL hidden within its exomes. However, beyond the scope of this one disease, the experimental and analytical methods can be readily applied to other cancer sequencing studies.

Thus, this dissertation leverages next-generation sequencing analysis to understand the genetic underpinnings of lymphoma, both by examining its normal cells of origin as well as through a large-scale study to sensitively identify recurrently mutated genes and their relationship to clinical outcome.

Synergy between histone deacetylase inhibitors and DNA- damaging agents is mediated by histone deacetylase 2 in colorectal cancer

Previous studies have associated the overexpression of histone deacetylase 2 (HDAC2) and the presence of TP53 mutations with the progression to advanced stage drug resistant colorectal cancer (CRC). However, the mechanistic link between HDAC2 expression and the TP53 mutational status has remained unexplored. Here, we investigated the function of HDAC2 in drug resistance by assessing the synergistic effects of DNA-targeted chemotherapeutic agents and HDAC inhibitors (HDACis) on two TP53-mutated colorectal adenocarcinoma CRC cell lines (SW480 and HT-29) and on the TP53-wild type carcinoma cell line (HCT116 p53+/+) and its TP53 deficient sub-line (HCT116 p53-/-). We showed that in the untreated SW480 and HT-29 cells the steady-state level of HDAC2 was low compared to a TP53-wild type carcinoma cell line (HCT116 p53+/+). Increased expression of HDAC2 correlated with drug resistance, and depletion by shRNA sensitised the multi-drug resistance cell line HT-29 to CRC chemotherapeutic drugs such as 5-fluorouracil (5-FU) and oxaliplatin (Oxa). Combined treatment with the HDACi suberoylanilide hydroxamic acid plus 5-FU or Oxa reduced the level of HDAC2 expression, modified chromatin structure and induced mitotic cell death in HT-29 cells. Non-invasive bioluminescence imaging revealed significant reductions in xenograft tumour growth with HDAC2 expression level reduced to <50% in treated animals. Elevated levels of histone acetylation on residues H3K9, H4K12 and H4K16 were also found to be associated with resistance to VPA/Dox or SAHA/Dox treatment. Our results suggest that HDAC2 expression rather than the p53 mutation status influences the outcome of combined treatment with a HDACi and DNA-damaging agents in CRC.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND: Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk. METHOD: We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation. RESULTS: Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05. CONCLUSION: This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk

Epigenetic Therapy for the Treatment of Hypertension-Induced Cardiac Hypertrophy and Fibrosis

BACKGROUND: The development of heart failure is associated with changes in the size, shape, and structure of the heart that has a negative impact on cardiac function. These pathological changes involve excessive extracellular matrix deposition within the myocardial interstitium and myocyte hypertrophy. Alterations in fibroblast phenotype and myocyte activity are associated with reprogramming of gene transcriptional profiles that likely requires epigenetic alterations in chromatin structure. The aim of our work was to investigate the potential of a currently licensed anticancer epigenetic modifier as a treatment option for cardiac diseases associated with hypertension-induced cardiac hypertrophy and fibrosis.

METHODS AND RESULTS: The effects of DNA methylation inhibition with 5-azacytidine (5-aza) were examined in a human primary fibroblast cell line and in a spontaneously hypertensive rat (SHR) model. The results from this work allude to novel in vivo antifibrotic and antihypertrophic actions of 5-aza. Administration of the DNA methylation inhibitor significantly improved several echocardiographic parameters associated with hypertrophy and diastolic dysfunction. Myocardial collagen levels and myocyte size were reduced in 5-aza-treated SHRs. These findings are supported by beneficial in vitro effects in cardiac fibroblasts. Collagen I, collagen III, and α-smooth muscle actin were reduced in a human ventricular cardiac fibroblast cell line treated with 5-aza.

CONCLUSION: These findings suggest a role for epigenetic modifications in contributing to the profibrotic and hypertrophic changes evident during disease progression. Therapeutic intervention with 5-aza demonstrated favorable effects highlighting the potential use of this epigenetic modifier as a treatment option for cardiac pathologies associated with hypertrophy and fibrosis.

Hypoxia-induced epigenetic modifications are associated with cardiac tissue fibrosis and the development of a myofibroblast-like phenotype

Ischemia caused by coronary artery disease and myocardial infarction leads to aberrant ventricular remodeling and cardiac fibrosis. This occurs partly through accumulation of gene expression changes in resident fibroblasts, resulting in an overactive fibrotic phenotype. Long-term adaptation to a hypoxic insult is likely to require significant modification of chromatin structure in order to maintain the fibrotic phenotype. Epigenetic changes may play an important role in modulating hypoxia-induced fibrosis within the heart. Therefore, the aim of the study was to investigate the potential pro-fibrotic impact of hypoxia on cardiac fibroblasts and determine whether alterations in DNA methylation could play a role in this process. This study found that within human cardiac tissue, the degree of hypoxia was associated with increased expression of collagen 1 and alpha-smooth muscle actin (ASMA). In addition, human cardiac fibroblast cells exposed to prolonged 1% hypoxia resulted in a pro-fibrotic state. These hypoxia-induced pro-fibrotic changes were associated with global DNA hypermethylation and increased expression of the DNA methyltransferase (DNMT) enzymes DNMT1 and DNMT3B. Expression of these methylating enzymes was shown to be regulated by hypoxia-inducible factor (HIF)-1α. Using siRNA to block DNMT3B expression significantly reduced collagen 1 and ASMA expression. In addition, application of the DNMT inhibitor 5-aza-2'-deoxycytidine suppressed the pro-fibrotic effects of TGFβ. Epigenetic modifications and changes in the epigenetic machinery identified in cardiac fibroblasts during prolonged hypoxia may contribute to the pro-fibrotic nature of the ischemic milieu. Targeting up-regulated expression of DNMTs in ischemic heart disease may prove to be a valuable therapeutic approach.

Mapping Protein-DNA Interactions Using ChIP-exo and Illumina-Based Sequencing

Chromatin immunoprecipitation (ChIP) provides a means of enriching DNA associated with transcription factors, histone modifications, and indeed any other proteins for which suitably characterized antibodies are available. Over the years, sequence detection has progressed from quantitative real-time PCR and Southern blotting to microarrays (ChIP-chip) and now high-throughput sequencing (ChIP-seq). This progression has vastly increased the sequence coverage and data volumes generated. This in turn has enabled informaticians to predict the identity of multi-protein complexes on DNA based on the overrepresentation of sequence motifs in DNA enriched by ChIP with a single antibody against a single protein. In the course of the development of high-throughput sequencing, little has changed in the ChIP methodology until recently. In the last three years, a number of modifications have been made to the ChIP protocol with the goal of enhancing the sensitivity of the method and further reducing the levels of nonspecific background sequences in ChIPped samples. In this chapter, we provide a brief commentary on these methodological changes and describe a detailed ChIP-exo method able to generate narrower peaks and greater peak coverage from ChIPped material.

Changes of 5-hydroxymethylcytosine distribution during myeloid and lymphoid differentiation of CD34+ cells

BACKGROUND: Hematopoietic stem cell renewal and differentiation are regulated through epigenetic processes. The conversion of 5-methylcytosine into 5-hydroxymethylcytosine (5hmC) by ten-eleven-translocation enzymes provides new insights into the epigenetic regulation of gene expression during development. Here, we studied the potential gene regulatory role of 5hmC during human hematopoiesis.

RESULTS: We used reduced representation of 5-hydroxymethylcytosine profiling (RRHP) to characterize 5hmC distribution in CD34+ cells, CD4+ T cells, CD19+ B cells, CD14+ monocytes and granulocytes. In all analyzed blood cell types, the presence of 5hmC at gene bodies correlates positively with gene expression, and highest 5hmC levels are found around transcription start sites of highly expressed genes. In CD34+ cells, 5hmC primes for the expression of genes regulating myeloid and lymphoid lineage commitment. Throughout blood cell differentiation, intragenic 5hmC is maintained at genes that are highly expressed and required for acquisition of the mature blood cell phenotype. Moreover, in CD34+ cells, the presence of 5hmC at enhancers associates with increased binding of RUNX1 and FLI1, transcription factors essential for hematopoiesis.

CONCLUSIONS: Our study provides a comprehensive genome-wide overview of 5hmC distribution in human hematopoietic cells and new insights into the epigenetic regulation of gene expression during human hematopoiesis.

Intraclonal heterogeneity and distinct molecular mechanisms characterize the development of t(4;14) and t(11;14) myeloma.

We have used whole exome sequencing to compare a group of presentation t(4;14) with t(11;14) cases of myeloma to define the mutational landscape. Each case was characterized by a median of 24.5 exonic nonsynonymous single-nucleotide variations, and there was a consistently higher number of mutations in the t(4;14) group, but this number did not reach statistical significance. We show that the transition and transversion rates in the 2 subgroups are similar, suggesting that there was no specific mechanism leading to mutation differentiating the 2 groups. Only 3% of mutations were seen in both groups, and recurrently mutated genes include NRAS, KRAS, BRAF, and DIS3 as well as DNAH5, a member of the axonemal dynein family. The pattern of mutation in each group was distinct, with the t(4;14) group being characterized by deregulation of chromatin organization, actin filament, and microfilament movement. Recurrent RAS pathway mutations identified subclonal heterogeneity at a mutational level in both groups, with mutations being present as either dominant or minor subclones. The presence of subclonal diversity was confirmed at a single-cell level using other tumor-acquired mutations. These results are consistent with a distinct molecular pathogenesis underlying each subgroup and have important impacts on targeted treatment strategies. The Medical Research Council Myeloma IX trial is registered under ISRCTN68454111.

A global expression-based analysis of the consequences of the t(4;14) translocation in myeloma

Purpose: Our purpose in this report was to define genes and pathways dysregulated as a consequence of the t(4;14) in myeloma, and to gain insight into the downstream functional effects that may explain the different prognosis of this subgroup.Experimental Design: Fibroblast growth factor receptor 3 (FGFR3) overexpression, the presence of immunoglobulin heavy chain-multiple myeloma SET domain (IgH-MMSET) fusion products and the identification of t(4;14) breakpoints were determined in a series of myeloma cases. Differentially expressed genes were identified between cases with (n = 55) and without (n = 24) a t(4;14) by using global gene expression analysis.Results: Cases with a t(4;14) have a distinct expression pattern compared with other cases of myeloma. A total of 127 genes were identified as being differentially expressed including MMSET and cyclin D2, which have been previously reported as being associated with this translocation. Other important functional classes of genes include cell signaling, apoptosis and related genes, oncogenes, chromatin structure, and DNA repair genes. Interestingly, 25% of myeloma cases lacking evidence of this translocation had up-regulation of the MMSET transcript to the same level as cases with a translocation.Conclusions: t(4;14) cases form a distinct subgroup of myeloma cases with a unique gene signature that may account for their poor prognosis. A number of non-t(4;14) cases also express MMSET consistent with this gene playing a role in myeloma pathogenesis.

Bcl-xL (S49) and (S62) sequential phosphorylation/dephosphorylation during mitosis prevents chromosome instability and aneuploidy

Une caractéristique intéressante de la protéine Bcl-xL est la présence d'un domaine en boucle non-structurée entre les hélices α1 and α2 de la protéine. Ce domaine protéique n'est pas essentiel pour sa fonction anti-apoptotique et absent chez CED-9, la protéine orthologue chez Caenorhabditis elegans. A l'intérieur de ce domaine, Bcl-xL subit une phosphorylation et déphosphorylation dynamique sur les résidus Ser49 et Ser62 en phase G2 du cycle cellulaire et lors de la mitose. Lorsque ces résidus sont mutés et les protéines exprimées dans des cellules cancéreuses, les cellules démontrent plusieurs défauts mitotiques liés à l'instabilité chromosomique. Pour analyser les effets de Bcl-xL Ser49 et Ser62 dans les cellules normales, les présentes études ont été réalisées dans des cellules diploïdes humaines normales, et in vivo chez Caenorhabditis elegans. Dans une première étude, nous avons utilisé la lignée cellulaire de cellules fibroblastiques diploïdes humaines normales BJ, exprimant Bcl-xL (type sauvage), (S49A), (S49D), (S62A), (S62D) et les double (S49/62A) et (S49/62D) mutants. Les cellules exprimant les mutants de phosphorylation ont montré des cinétiques de doublement de la population cellulaire réduites. Ces effets sur la cinétique de doublement de la population cellulaire corrèle avec l'apparition de la sénescence cellulaire, sans impact sur les taux de mort cellulaire. Ces cellules sénescentes affichent des phénotypes typiques de sénescence associés notamment à haut niveau de l'activité β-galactosidase associée à la sénescence, la sécrétion d' interleukine-6, l'activation de p53 et de p21WAF1/ Cip1, un inhibiteur des complexes kinase cycline-dépendant, ainsi que la formation de foyers de chromatine nucléaire associés à γH2A.X. Les analyses de fluorescence par hybridation in situ et des caryotypes par coloration au Giemsa ont révélé que l'expression des mutants de phosphorylation de Bcl-xL provoquent de l'instabilité chromosomique et l'aneuploïdie. Ces résultats suggèrent que les cycles de phosphorylation et déphosphorylation dynamiques de Bcl-xL Ser49 et Ser62 sont importants dans le maintien de l'intégrité des chromosomes lors de la mitose dans les cellules normales. Dans une deuxième étude, nous avons entrepris des expériences chez Caenorhabditis elegans pour comprendre l'importance des résidus Ser49 et Ser62 de Bcl-xL in vivo. Les vers transgéniques portant les mutations de Bcl-xL (S49A, S62A, S49D, S62D et S49/62A) ont été générés et leurs effets ont été analysés sur les cellules germinales des jeunes vers adultes. Les vers portant les mutations de Bcl-xL ont montré une diminution de ponte et d'éclosion des oeufs, des variations de la longueur de leurs régions mitotiques et des zones de transition, des anomalies chromosomiques à leur stade de diplotène, et une augmentation de l'apoptose des cellules germinales. Certaines de ces souches transgéniques, en particulier les variants Ser/Ala, ont également montré des variations de durée de vie par rapport aux vers témoins. Ces observations in vivo ont confirmé l'importance de Ser49 et Ser62 à l'intérieur du domaine à boucle de Bcl-xL pour le maintien de la stabilité chromosomique. Ces études auront une incidence sur les futures stratégies visant à développer et à identifier des composés qui pourraient cibler non seulement le domaine anti-apoptotique de la protéine Bcl-xL, mais aussi son domaine mitotique pour la thérapie du cancer.

Étude fonctionnelle d’un nouveau complexe multi-enzymatique régulant l’épigénome

L’ubiquitination, une modification post-traductionnelle importante pour le contrôle de nombreux processus cellulaires, est une réaction réversible. La réaction inverse, nommée déubiquitination est catalysée par les déubiquitinases (DUB). Nous nous sommes intéressés dans nos travaux à étudier l’ubiquitination de l’histone H2A (H2Aub), au niveau des résidus lysines 118 et 119 (K118/K119), une marque épigénétique impliquée dans la régulation de la prolifération cellulaire et la réparation de l’ADN. Le régulateur transcriptionnel BAP1, une déubiquitinase nucléaire, a été initialement identifié pour sa capacité à promouvoir la fonction suppressive de tumeurs de BRCA1. BAP1 forme un complexe multi-protéique avec plusieurs facteurs transcriptionnels et sa fonction principale est la déubiquitination de H2Aub. Plusieurs études ont démontré que BAP1 est un gène suppresseur de tumeurs majeur et qu’il est largement muté et inactivé dans une multitude de cancers. En effet, BAP1 émerge comme étant la DUB la plus mutée au niveau des cancers. Cependant, le ou les mécanismes d’action et de régulation du complexe BAP1 restent très peu connus. Dans cette étude nous nous sommes intéressés à la caractérisation moléculaire et fonctionnelle des partenaires protéiques de BAP1. De manière significative nous avons caractérisé un mécanisme unique de régulation entre deux composants majeurs du complexe BAP1 à savoir, HCF-1 et OGT. En effet, nous avons démontré que HCF-1 est requis pour maintenir le niveau protéique de OGT et que cette dernière est indispensable pour la maturation protéolytique de HCF-1 en promouvant son clivage par O-GlcNAcylation, une signalisation cellulaire nécessaire au bon fonctionnement de HCF-1. Également, nous avons découvert un nouveau mécanisme de régulation de BAP1 par l’ubiquitine ligase atypique UBE2O. En effet, UBE2O agit comme un régulateur négatif de BAP1 puisque l’ubiquitination de ce dernier induit sa séquestration dans le cytoplasme et l’inhibition de sa fonction suppressive de tumeurs. D’autre part nous nous sommes penchés sur la caractérisation de l’association de BAP1 avec deux facteurs de la famille des protéines Polycombes nommés ASXL1 et ASXL2 (ASXL1/2). Nous avons investigué le rôle de BAP1/ASXL1/2, particulièrement dans les mécanismes de déubiquitination et suppression de tumeurs. Nous avons démontré que BAP1 interagit directement iii via son domaine C-terminale avec le même domaine ASXM de ASXL1/2 formant ainsi deux complexes mutuellement exclusifs indispensables pour induire l’activité déubiquitinase de BAP1. De manière significative, ASXM s’associe avec BAP1 pour créer un nouveau domaine composite de liaison à l’ubiquitine. Ces interactions BAP1/ASXL1/2 régulent la progression harmonieuse du cycle cellulaire. De plus, la surexpression de BAP1 et de ASXL2 au niveau des fibroblastes induit la sénescence de manière dépendante de leurs interactions. D’autre part, nous avons identifié des mutations de cancers au niveau de BAP1 le rendant incapable de lier ASXL1/2, d’exercer sa fonction d’autodéubiquitination et de ce fait d’agir comme suppresseur de tumeurs. Ainsi nous avons révélé un lien étroit entre le gène suppresseur de tumeurs BAP1, son activité déubiquitinase et le contrôle de la prolifération cellulaire.

Flexible strategies for association analysis with genomic phenotypes

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

Repeated evolution of heat responsiveness among Brassicaceae COPIA transposable elements

Eukaryotic genomes contain repetitive DNA sequences. This includes simple repeats and more complex transposable elements (TEs). Many TEs reach high copy numbers in the host genome, owing to their amplification abilities by specific mechanisms. There is growing evidence that TEs contribute to gene transcriptional regulation. However, excess of TE activity may lead to reduced genome stability. Therefore, TEs are suppressed by the transcriptional gene silencing machinery via specific chromatin modifications. In contrary, effectiveness of the epigenetic silencing mechanisms imposes risk for TE survival in the host genome. Therefore, TEs may have evolved specific strategies for bypassing epigenetic control and allowing the emergence of new TE copies. Recent studies suggested that the epigenetic silencing can be, at least transiently, attenuated by heat stress in A. thaliana. Heat stress induced strong transcriptional activation of COPIA78 family LTR-retrotransposons named ONSEN, and even their transposition in mutants deficient in siRNA-biogenesis. ONSEN transcriptional activation was facilitated by the presence of heat responsive elements (HREs) within the long terminal repeats, which serve as a binding platform for the HEAT SHOCK FACTORs (HSFs). This thesis focused on the evolution of ONSEN heat responsiveness in Brassicaceae. By using whole-transcriptome sequencing approach, multiple Arabidopsis lyrata ONSENs with conserved heat response were found and together with ONSENs from other Brassicaceae were used to reconstruct the evolution of ONSEN HREs. This indicated ancestral situation with two, in palindrome organized, HSF binding motifs. In the genera Arabidopsis and Ballantinia, a local duplication of this locus increased number of HSF binding motifs to four, forming a high-efficiency HRE. In addition, whole transcriptome analysis revealed novel heat-responsive TE families COPIA20, COPIA37 and HATE. Notably, HATE represents so far unknown COPIA family which occurs in several Brassicaceae species but is absent in A. thaliana. Putative HREs were identified within the LTRs of COPIA20, COPIA37 and HATE of A. lyrata, and could be preliminarily validated by transcriptional analysis upon heat induction in subsequent survey of Brassicaeae species. Subsequent phylogenetic analysis indicated a repeated evolution of heat responsiveness within Brassicaceae COPIA LTR-retrotransposons. This indicates that acquisition of heat responsiveness may represent a successful strategy for survival of TEs within the host genome.

Reorganisation ofHoxdregulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

FACS-based purification of Arabidopsis microspores, sperm cells and vegetative nuclei

Background: The male germline in flowering plants differentiates by asymmetric division of haploid uninucleated microspores, giving rise to a vegetative cell enclosing a smaller generative cell, which eventually undergoes a second mitosis to originate two sperm cells. The vegetative cell and the sperm cells activate distinct genetic and epigenetic mechanisms to control pollen tube growth and germ cell specification, respectively. Therefore, a comprehensive characterization of these processes relies on efficient methods to isolate each of the different cell types throughout male gametogenesis. Results: We developed stable transgenic Arabidopsis lines and reliable purification tools based on Fluorescence-Activated Cell Sorting (FACS) in order to isolate highly pure and viable fractions of each cell/nuclei type before and after pollen mitosis. In the case of mature pollen, this was accomplished by expressing GFP and RFP in the sperm and vegetative nuclei, respectively, resulting in 99% pure sorted populations. Microspores were also purified by FACS taking advantage of their characteristic small size and autofluorescent properties, and were confirmed to be 98% pure. Conclusions: We provide simple and efficient FACS-based purification protocols for Arabidopsis microspores, vegetative nuclei and sperm cells. This paves the way for subsequent molecular analysis such as transcriptomics, DNA methylation analysis and chromatin immunoprecipitation, in the developmental context of microgametogenesis in Arabidopsis.