Sandra Louzada

Punctuated bursts in human male demography inferred from 1,244 worldwide Y-chromosome sequences

We report the sequences of 1,244 human Y chromosomes randomly ascertained from 26 worldwide populations by the 1000 Genomes Project. We discovered more than 65,000 variants, including single-nucleotide variants, multiple-nucleotide variants, insertions and deletions, short tandem repeats, and copy number variants. Of these, copy number variants contribute the greatest predicted functional impact. We constructed a calibrated phylogenetic tree on the basis of binary single-nucleotide variants and projected the more complex variants onto it, estimating the number of mutations for each class. Our phylogeny shows bursts of extreme expansion in male numbers that have occurred independently among each of the five continental superpopulations examined, at times of known migrations and technological innovations.

Multiplexed pancreatic genome engineering and cancer induction by transfection-based CRISPR/Cas9 delivery in mice

Mouse transgenesis has provided fundamental insights into pancreatic cancer, but is limited by the long duration of allele/model generation. Here we show transfection-based multiplexed delivery of CRISPR/Cas9 to the pancreas of adult mice, allowing simultaneous editing of multiple gene sets in individual cells. We use the method to induce pancreatic cancer and exploit CRISPR/Cas9 mutational signatures for phylogenetic tracking of metastatic disease. Our results demonstrate that CRISPR/Cas9-multiplexing enables key applications, such as combinatorial gene-network analysis, in vivo synthetic lethality screening and chromosome engineering. Negative-selection screening in the pancreas using multiplexed-CRISPR/Cas9 confirms the vulnerability of pancreatic cells to Brca2-inactivation in a Kras-mutant context. We also demonstrate modelling of chromosomal deletions and targeted somatic engineering of inter-chromosomal translocations, offering multifaceted opportunities to study complex structural variation, a hallmark of pancreatic cancer. The low-frequency mosaic pattern of transfection-based CRISPR/Cas9 delivery faithfully recapitulates the stochastic nature of human tumorigenesis, supporting wide applicability for biological/preclinical research.

Quantitative Genetics of CTCF Binding Reveal Local Sequence Effects and Different Modes of X-Chromosome Association

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

The pig X and Y Chromosomes: structure, sequence, and evolution.

We have generated an improved assembly and gene annotation of the pig X Chromosome, and a first draft assembly of the pig Y Chromosome, by sequencing BAC and fosmid clones from Duroc animals and incorporating information from optical mapping and fiber-FISH. The X Chromosome carries 1033 annotated genes, 690 of which are protein coding. Gene order closely matches that found in primates (including humans) and carnivores (including cats and dogs), which is inferred to be ancestral. Nevertheless, several protein-coding genes present on the human X Chromosome were absent from the pig, and 38 pig-specific X-chromosomal genes were annotated, 22 of which were olfactory receptors. The pig Y-specific Chromosome sequence generated here comprises 30 megabases (Mb). A 15-Mb subset of this sequence was assembled, revealing two clusters of male-specific low copy number genes, separated by an ampliconic region including the HSFY gene family, which together make up most of the short arm. Both clusters contain palindromes with high sequence identity, presumably maintained by gene conversion. Many of the ancestral X-related genes previously reported in at least one mammalian Y Chromosome are represented either as active genes or partial sequences. This sequencing project has allowed us to identify genes--both single copy and amplified--on the pig Y Chromosome, to compare the pig X and Y Chromosomes for homologous sequences, and thereby to reveal mechanisms underlying pig X and Y Chromosome evolution.

Defining the Sister Rat Mammary Tumor Cell Lines HH-16 cl.2/1 and HH-16.cl.4 as an In Vitro Cell Model for Erbb2

Cancer cell lines have been shown to be reliable tools in genetic studies of breast cancer, and the characterization of these lines indicates that they are good models for studying the biological mechanisms underlying this disease. Here, we describe the molecular cytogenetic/genetic characterization of two sister rat mammary tumor cell lines, HH-16 cl.2/1 and HH-16.cl.4, for the first time. Molecular cytogenetic analysis using rat and mouse chromosome paint probes and BAC/PAC clones allowed the characterization of clonal chromosome rearrangements; moreover, this strategy assisted in revealing detected breakpoint regions and complex chromosome rearrangements. This comprehensive cytogenetic analysis revealed an increase in the number of copies of the Mycn and Erbb2 genes in the investigated cell lines. To analyze its possible correlation with expression changes, relative RNA expression was assessed by real-time reverse transcription quantitative PCR and RNA FISH. Erbb2 was found to be overexpressed in HH-16.cl.4, but not in the sister cell line HH-16 cl.2/1, even though these lines share the same initial genetic environment. Moreover, the relative expression of Erbb2 decreased after global genome demethylation in the HH-16.cl.4 cell line. As these cell lines are commercially available and have been used in previous studies, the present detailed characterization improves their value as an in vitro cell model. We believe that the development of appropriate in vitro cell models for breast cancer is of crucial importance for revealing the genetic and cellular pathways underlying this neoplasy and for employing them as experimental tools to assist in the generation of new biotherapies.

Different evolutionary trails in the related genomes Cricetus cricetus and Peromyscus eremicus (Rodentia, Cricetidae) uncovered by orthologous satellite DNA repositioning

Constitutive heterochromatin comprises a substantial fraction of the eukaryotic genomes and is mainly composed of tandemly arrayed satellite DNAs (satDNA). These repetitive sequences represent a very dynamic and fast evolving component of genomes. In the present work we report the isolation of Cricetus cricetus (CCR, Cricetidae, Rodentia) centromeric repetitive sequences from chromosome 4 (CCR4/10sat), using the laser microdissection and laser pressure catapulting procedure, followed by DOP-PCR amplification and labelling. Physical mapping by fluorescent in situ hybridisation of these sequences onto C. cricetus and another member of Cricetidae, Peromyscus eremicus, displayed quite interesting patterns. Namely, the centromeric sequences showed to be present in another C. cricetus chromosome (CCR10) besides CCR4. Moreover, these almost chromosome-specific sequences revealed to be present in the P. eremicus genome, and most interestingly, displaying a ubiquitous scattered distribution throughout this karyotype. Finally and in both species, a co-localisation of CCR4/10sat with constitutive heterochromatin was found, either by classical C-banding or C-banding sequential to in situ endonuclease restriction. The presence of these orthologous sequences in both genomes is suggestive of a phylogenetic proximity. Furthermore, the existence of common repetitive DNA sequences with a different chromosomal location foresees the occurrence of an extensive process of karyotype restructuring somehow related with intragenomic movements of these repetitive sequences during the evolutionary process of C. cricetus and P. eremicus species.

Alcohol and endogenous aldehydes damage chromosomes and mutate stem cells

Haematopoietic stem cells renew blood. Accumulation of DNA damage in these cells promotes their decline, while misrepair of this damage initiates malignancies. Here we describe the features and mutational landscape of DNA damage caused by acetaldehyde, an endogenous and alcohol-derived metabolite. This damage results in DNA double-stranded breaks that, despite stimulating recombination repair, also cause chromosome rearrangements. We combined transplantation of single haematopoietic stem cells with whole-genome sequencing to show that this damage occurs in stem cells, leading to deletions and rearrangements that are indicative of microhomology-mediated end-joining repair. Moreover, deletion of p53 completely rescues the survival of aldehyde-stressed and mutated haematopoietic stem cells, but does not change the pattern or the intensity of genome instability within individual stem cells. These findings characterize the mutation of the stem-cell genome by an alcohol-derived and endogenous source of DNA damage. Furthermore, we identify how the choice of DNA-repair pathway and a stringent p53 response limit the transmission of aldehyde-induced mutations in stem cells.

Evolution of the rapidly mutating human salivary agglutinin gene (DMBT1) and population subsistence strategy.

Significance Humans have undergone an evolutionary very recent change in environment of their own making. The development of agriculture profoundly altered diet and exposure to pathogens, and yet the evolutionary response to this is still poorly understood. Here, we characterize extensive copy number variation (CNV) of the gene encoding salivary agglutinin (deleted in malignant brain tumors 1, DMBT1). Salivary agglutinin comprises 10% of salivary protein and binds bacteria, including mediating the attachment of the causative agent of dental caries, Streptococcus mutans, to teeth. We show that DMBT1 is a very fast-mutating protein-coding locus, and DMBT1 CNV correlates with a population history of agriculture. Furthermore, we examine the relationship between variation of the S. mutans region that binds salivary agglutinin and CNV of the DMBT1 gene. The dietary change resulting from the domestication of plant and animal species and development of agriculture at different locations across the world was one of the most significant changes in human evolution. An increase in dietary carbohydrates caused an increase in dental caries following the development of agriculture, mediated by the cariogenic oral bacterium Streptococcus mutans. Salivary agglutinin [SAG, encoded by the deleted in malignant brain tumors 1 (DMBT1) gene] is an innate immune receptor glycoprotein that binds a variety of bacteria and viruses, and mediates attachment of S. mutans to hydroxyapatite on the surface of the tooth. In this study we show that multiallelic copy number variation (CNV) within DMBT1 is extensive across all populations and is predicted to result in between 7–20 scavenger–receptor cysteine-rich (SRCR) domains within each SAG molecule. Direct observation of de novo mutation in multigeneration families suggests these CNVs have a very high mutation rate for a protein-coding locus, with a mutation rate of up to 5% per gamete. Given that the SRCR domains bind S. mutans and hydroxyapatite in the tooth, we investigated the association of sequence diversity at the SAG-binding gene of S. mutans, and DMBT1 CNV. Furthermore, we show that DMBT1 CNV is also associated with a history of agriculture across global populations, suggesting that dietary change as a result of agriculture has shaped the pattern of CNV at DMBT1, and that the DMBT1-S. mutans interaction is a promising model of host-pathogen-culture coevolution in humans.

Evolutionary routes and KRAS dosage define pancreatic cancer phenotypes

The poor correlation of mutational landscapes with phenotypes limits our understanding of the pathogenesis and metastasis of pancreatic ductal adenocarcinoma (PDAC). Here we show that oncogenic dosage-variation has a critical role in PDAC biology and phenotypic diversification. We find an increase in gene dosage of mutant KRAS in human PDAC precursors, which drives both early tumorigenesis and metastasis and thus rationalizes early PDAC dissemination. To overcome the limitations posed to gene dosage studies by the stromal richness of PDAC, we have developed large cell culture resources of metastatic mouse PDAC. Integration of cell culture genomes, transcriptomes and tumour phenotypes with functional studies and human data reveals additional widespread effects of oncogenic dosage variation on cell morphology and plasticity, histopathology and clinical outcome, with the highest KrasMUT levels underlying aggressive undifferentiated phenotypes. We also identify alternative oncogenic gains (Myc, Yap1 or Nfkb2), which collaborate with heterozygous KrasMUT in driving tumorigenesis, but have lower metastatic potential. Mechanistically, different oncogenic gains and dosages evolve along distinct evolutionary routes, licensed by defined allelic states and/or combinations of hallmark tumour suppressor alterations (Cdkn2a, Trp53, Tgfβ-pathway). Thus, evolutionary constraints and contingencies direct oncogenic dosage gain and variation along defined routes to drive the early progression of PDAC and shape its downstream biology. Our study uncovers universal principles of Ras-driven oncogenesis that have potential relevance beyond pancreatic cancer.

Praomys tullbergi (Muridae, Rodentia) genome architecture decoded by comparative chromosome painting with Mus and Rattus

The order Rodentia and in particular the Muridae are characterised by extremely high rates of chromosome evolution and remarkable chromosome diversity. The Praomys group (Murinae, Muridae and Rodentia) constitutes a diverse and abundant group divided into two complexes, the jacksoni complex and the tullbergi complex which includes the species Praomys tullbergi. Comparative chromosome painting using the two index genomes, Mus musculus and Rattus norvegicus, was performed resulting in a high resolution chromosome map for P. tullbergi. The combined use of rat and mouse probes and the assistance of the assembly of all the available sequencing data from Ensembl genome browser allowed a great dissection of P. tullbergi genome, the detection of inversion events and ultimately the refinement of P. tullbergi comparative map. A key achievement was the reconstruction of a high precision Muroidea ancestral karyotype (Muridae/Cricetidae and Murine) based in a broad species analysis combining previous reported comparative maps together with the presented data. This permitted the reconstruction of the evolutionary history of chromosome changes since the ancestral Muroidea genome and enlightened the phylogenetic relationships with the related species mouse and rat. The analysis of constitutive heterochromatin and its co-localisation with the identified evolutionary breakpoints regions was performed suggesting the involvement of repetitive sequences in the chromosome rearrangements that originated the present P. tullbergi genome architecture.