Investigating loss of heterozygosity in a SCRaMbLEd yeast genome

Science China Life Sciences -     

GPHMM: an integrated hidden Markov model for identification of copy number alteration and loss of heterozygosity in complex tumor samples using whole genome SNP arrays

There is an increasing interest in using single nucleotide polymorphism (SNP) genotyping arrays for profiling chromosomal rearrangements in tumors, as they allow simultaneous detection of copy number and loss of heterozygosity with high resolution. Critical issues such as signal baseline shift due to aneuploidy, normal cell contamination, and the presence of GC content bias have been reported to dramatically alter SNP array signals and complicate accurate identification of aberrations in cancer genomes. To address these issues, we propose a novel Global Parameter Hidden Markov Model (GPHMM) to unravel tangled genotyping data generated from tumor samples. In contrast to other HMM methods, a distinct feature of GPHMM is that the issues mentioned above are quantitatively modeled by global parameters and integrated within the statistical framework. We developed an efficient EM algorithm for parameter estimation. We evaluated performance on three data sets and show that GPHMM can correctly identify chromosomal aberrations in tumor samples containing as few as 10% cancer cells. Furthermore, we demonstrated that the estimation of global parameters in GPHMM provides information about the biological characteristics of tumor samples and the quality of genotyping signal from SNP array experiments, which is helpful for data quality control and outlier detection in cohort studies.     

Telomeres, hidden mosaicism, loss of heterozygosity, and complex genetic traits

Telomeres appear to function as an endogenous timing mechanism in human beings. Telomere attrition not only provides a satisfactory explanation for some aspects of aging, it might also resolve enigmatic features of complex genetic traits that are age-dependent. If, with the passage of time, telomere attrition in human beings leads to genomic instability and particularly the loss of chromosomes, then the age dependency of phenotypic expressions of complex genetic traits might result from the temporal loss of heterozygosity and the consequent expression of disease-causing genes. In this way, telomere attrition might play a role not only in aging, but also in the diverse expression of complex genetic traits, such as essential hypertension, non-insulin-dependent diabetes mellitus, atherosclerosis, and cancer.     

Somatic mosaicism for copy-neutral loss of heterozygosity and DNA copy number variations in the human genome

Background

Somatic mosaicism denotes the presence of genetically distinct populations of somatic cells in one individual who has developed from a single fertilised oocyte. Mosaicism may result from a mutation that occurs during postzygotic development and is propagated to only a subset of the adult cells. Our aim was to investigate both somatic mosaicism for copy-neutral loss of heterozygosity (cn-LOH) events and DNA copy number variations (CNVs) in fully differentiated tissues.

Results

We studied panels of tissue samples (11–12 tissues per individual) from four autopsy subjects using high-resolution Illumina HumanOmniExpress-12 BeadChips to reveal the presence of possible intra-individual tissue-specific cn-LOH and CNV patterns.We detected five mosaic cn-LOH regions >5 Mb in some tissue samples in three out of four individuals. We also detected three CNVs that affected only a portion of the tissues studied in one out of four individuals. These three somatic CNVs range from 123 to 796 kb and are also found in the general population. An attempt was made to explain the succession of genomic events that led to the observed somatic genetic mosaicism under the assumption that the specific mosaic patterns of CNV and cn-LOH changes reflect their formation during the postzygotic embryonic development of germinal layers and organ systems.

Conclusions

Our results give further support to the idea that somatic mosaicism for CNVs, and also cn-LOHs, is a common phenomenon in phenotypically normal humans. Thus, the examination of only a single tissue might not provide enough information to diagnose potentially deleterious CNVs within an individual. During routine CNV and cn-LOH analysis, DNA derived from a buccal swab can be used in addition to blood DNA to get information about the CNV/cn-LOH content in tissues of both mesodermal and ectodermal origin. Currently, the real frequency and possible phenotypic consequences of both CNVs and cn-LOHs that display somatic mosaicism remain largely unknown. To answer these questions, future studies should involve larger cohorts of individuals and a range of tissues.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1916-3) contains supplementary material, which is available to authorized users.     

Extreme diversification driven by parallel events of massive loss of heterozygosity in the hybrid lineage of Candida albicans

Candida albicans is the most commonly reported species causing candidiasis. The taxonomic classification of C. albicans and related lineages is controversial, with Candida africana (syn. C. albicans var. africana) and Candida stellatoidea (syn. C. albicans var. stellatoidea) being considered different species or C. albicans varieties depending on the authors. Moreover, recent genomic analyses have suggested a shared hybrid origin of C. albicans and C. africana, but the potential parental lineages remain unidentified. Although the genomes of C. albicans and C. africana have been extensively studied, the genome of C. stellatoidea has not been sequenced so far. In order to get a better understanding of the evolution of the C. albicans clade, and to assess whether C. stellatoidea could represent one of the unknown C. albicans parental lineages, we sequenced C. stellatoidea type strain (CBS 1905). This genome was compared to that of C. albicans and of the closely related lineage C. africana. Our results show that, similarly to C. africana, C. stellatoidea descends from the same hybrid ancestor as other C. albicans strains and that it has undergone a parallel massive loss of heterozygosity.     

Synaptonemal complex karyotype of zebrafish

Meiotic cells of zebrafish have been prepared to show synaptonemal complexes (SCs) by light and electron microscopy. Completely paired SCs from both spermatocytes and oocytes were measured to produce an SC karyotype. The SC karyotype resembles the somatic karyotype of zebrafish and has no recognisable sex bivalent. Measurements of total SC length indicate that SCs grow longer and develop centromeres during pachytene. Oocytes consistently have longer SCs than spermatocytes, presumably correlated with the reported higher recombination frequency in females than in males.     

The evolution of vertebrate tetraspanins: gene loss,retention, and massive positive selection after whole genome duplications

Background  

The vertebrate tetraspanin family has many features which make it suitable for preserving the imprint of ancient sequence evolution and amenable for phylogenomic analysis. So we believe that an in-depth analysis of the tetraspanin evolution not only provides more complete understanding of tetraspanin biology, but offers new insights into the influence of the two rounds of whole genome duplication (2R-WGD) at the origin of vertebrates.     

Fugu: a compact vertebrate reference genome

At 400 Mb, the Japanese pufferfish, Fugu rubripes, has the smallest vertebrate genome but has a similar gene repertoire to other vertebrates. Its genes are densely packed with short intergenic and intronic sequences devoid of repetitive elements. It likely has a mutational bias towards DNA elimination and is probably close to a 'minimal' vertebrate genome. As such it is a useful reference genome for gene discovery and gene validation in other vertebrates. Its usefulness in the discovery of conserved regulatory elements has already been demonstrated. The Fugu genome sequence is a good complement to genetic studies in other vertebrates.     

Synaptonemal complex karyotype of Eimeria tenella

In most organisms, biological variability rests on the behaviour of the chromosomes in the meiotic context. Despite the importance of meiosis, very little is known about the meiotic behaviour of the Eimeria chromosomes. The aim of the present study is to describe the standard synaptonemal complex karyotype from Eimeria tenella oocyst spreads by electron microscopy. For that purpose, complete sets of pachytene synaptonemal complexes were obtained and the morphological pachytene karyotype was determined. The authors used a previously reported method that overcomes the difficulty of the extreme resistance of protozoan oocysts to disruption and permits the release of intact meiotic chromosomes. The chromosomes were selected under a light microscope and those selected were stained with phosphotungtic acid and studied by transmission electron microscopy. The authors confirmed 14 chromosomes, which were observed as synaptonemal complexes, and the karyotype was constructed by arranging synaptonemal complexes according to their relative lengths and kinetochore position. Components of the synaptonemal complex, lateral elements, central element, recombination nodules and kinetochore were observed. Measures of the kynetochore, width of the synaptonemal complex, diameter of the recombination nodule and length of the telomeres are given. Minimal and no significant differences were found between measures of chromosomes isolated from different Eimeria tenella strains. To the best of our knowledge, the present investigation for the first time identifies and describes the morphological characteristics of the synaptonemal complex of Eimeria tenella during the meiosis that occurs within the oocysts. In addition, the authors provide evidence of the presence of recombination nodules, suggesting that the recombination process may play an important role in the molecular evolution of this parasite.     

Mechanisms of loss of heterozygosity in retinoblastoma.

Retinoblastoma (RB) tumors arise when both alleles of the RB1 gene are inactivated by two mutational events (M1 and M2). M1 can be an initial germline or somatic mutation; M2 is frequently loss of heterozygosity (LOH), which makes the cell homozygous or hemizygous for the original mutation. LOH is the major mechanism by which many cancers are initiated. To further delineate the mechanism of LOH, we screened a total of 37 RB tumors for LOH by Southern blot analysis. The tumors were from 17 bilaterally and 17 unilaterally affected patients. Nineteen of 30 informative tumors (63%) from 27 patients showed LOH. Proximal and distal flanking markers on chromosome 13 were informative in 13 tumors, allowing evaluation of the mechanisms by which LOH occurred. Mitotic recombination was implicated in 6 (46%) of the 13 tumors.     

SeqMap: mapping massive amount of oligonucleotides to the genome

SeqMap is a tool for mapping large amount of short sequences to the genome. It is designed for finding all the places in a reference genome where each sequence may come from. This task is essential to the analysis of data from ultra high-throughput sequencing machines. With a carefully designed index-filtering algorithm and an efficient implementation, SeqMap can map tens of millions of short sequences to a genome of several billions of nucleotides. Multiple substitutions and insertions/deletions of the nucleotide bases in the sequences can be tolerated and therefore detected. SeqMap supports FASTA input format and various output formats, and provides command line options for tuning almost every aspect of the mapping process. A typical mapping can be done in a few hours on a desktop PC. Parallel use of SeqMap on a cluster is also very straightforward.     

Molecular cytogenetic definition of the chicken genome: the first complete avian karyotype

Chicken genome mapping is important for a range of scientific disciplines. The ability to distinguish chromosomes of the chicken and other birds is thus a priority. Here we describe the molecular cytogenetic characterization of each chicken chromosome using chromosome painting and mapping of individual clones by FISH. Where possible, we have assigned the chromosomes to known linkage groups. We propose, on the basis of size, that the NOR chromosome is approximately the size of chromosome 22; however, we suggest that its original assignment of 16 should be retained. We also suggest a definitive chromosome classification system and propose that the probes developed here will find wide utility in the fields of developmental biology, DT40 studies, agriculture, vertebrate genome organization, and comparative mapping of avian species.     

Spontaneous loss of heterozygosity in diploid Saccharomyces cerevisiae cells

To obtain a broad perspective of the events leading to spontaneous loss of heterozygosity (LOH), we have characterized the genetic alterations that functionally inactivated the URA3 marker hemizygously or heterozygously situated either on chromosome III or chromosome V in diploid Saccharomyces cerevisiae cells. Analysis of chromosome structure in a large number of LOH clones by pulsed-field gel electrophoresis and PCR showed that chromosome loss, allelic recombination, and chromosome aberration were the major classes of genetic alterations leading to LOH. The frequencies of chromosome loss and chromosome aberration were significantly affected when the marker was located in different chromosomes, suggesting that chromosome-specific elements may affect the processes that led to these alterations. Aberrant-sized chromosomes were detected readily in approximately 8% of LOH events when the URA3 marker was placed in chromosome III. Molecular mechanisms underlying the chromosome aberrations were further investigated by studying the fate of two other genetic markers on chromosome III. Chromosome aberration caused by intrachromosomal rearrangements was predominantly due to a deletion between the MAT and HMR loci that occurred at a frequency of 3.1 x 10(-6). Another type of chromosome aberration, which occurred at a frequency slightly higher than that of the intrachromosomal deletion, appeared to be caused by interchromosomal rearrangement, including unequal crossing over between homologous chromatids and translocation with another chromosome.     

Circumventing heterozygosity: sequencing the amplified genome of a single haploid Drosophila melanogaster embryo

Heterozygosity is a major challenge to efficient, high-quality genomic assembly and to the full genomic survey of polymorphism and divergence. In Drosophila melanogaster lines derived from equatorial populations are particularly resistant to inbreeding, thus imposing a major barrier to the determination and analyses of genomic variation in natural populations of this model organism. Here we present a simple genome sequencing protocol based on the whole-genome amplification of the gynogenetically derived haploid genome of a progeny of females mated to males homozygous for the recessive male sterile mutation, ms(3)K81. A single "lane" of paired-end sequences (2 × 76 bp) provides a good syntenic assembly with >95% high-quality coverage (more than five reads). The amplification of the genomic DNA moderately inflates the variation in coverage across the euchromatic portion of the genome. It also increases the frequency of chimeric clones. But the low frequency and random genomic distribution of the chimeric clones limits their impact on the final assemblies. This method provides a solid path forward for population genomic sequencing and offers applications to many other systems in which small amounts of genomic DNA have unique experimental relevance.     

Pooled analysis of loss of heterozygosity in breast cancer: a genome scan provides comparative evidence for multiple tumor suppressors and identifies novel candidate regions

Somatic loss of heterozygosity (LOH) has been widely reported in breast cancer as a means of identifying putative tumor-suppressor genes. However, individual studies have rarely spanned more than a single chromosome, and the varying criteria used to declare LOH complicate efforts to formally differentiate regions of consistent versus sporadic (random) loss. We report here the compilation of an extensive database from 151 published LOH studies of breast cancer, with summary data from >15,000 tumors and primary allelotypes from >4,300 tumors. Allelic loss was evaluated at 1,168 marker loci, with large variation in the density of informative observations across the genome. Using studies in which primary allelotype information was available, we employed a likelihood-based approach with a formal chromosomal instability and selection model. The approach seeks direct evidence for preferential loss at each locus compared with nearby loci, accounts for heterogeneity across studies, and enables the direct comparison of candidate regions across the genome. Striking preferential loss was observed (in descending order of significance) in specific regions of chromosomes 7q, 16q, 13q, 17p, 8p, 21q, 3p, 18q, 2q, and 19p, as well as other regions, in many cases coinciding with previously identified candidate genes or known fragile sites. Many of these observations were not possible from any single LOH study, and our results suggest that many previously reported LOH results are not systematic or reproducible. Our approach provides a comparative framework for further investigation of regions exhibiting LOH and identifies broad genomic regions for which there exist few data.     

Soft selection reduces loss of heterozygosity in asexual reproduction

The adaptive value of sexual reproduction is still debated in evolutionary theory. It has been proposed that the advantage of sexual reproduction over asexual reproduction is to promote genetic diversity, to prevent the accumulation of harmful mutations or to preserve heterozygosity. Since these hypothetical advantages depend on the type of asexual reproduction, understanding how selection affects the taxonomic distribution of each type could help us discriminate between existing hypotheses. Here, I argue that soft selection, competition among embryos or offspring in selection arenas prior to the hard selection of the adult phase, reduces loss of heterozygosity in certain types of asexual reproduction. Since loss of heterozygosity leads to the unmasking of recessive deleterious mutations in the progeny of asexual individuals, soft selection facilitates the evolution of these types of asexual reproduction. Using a population genetics model, I calculate how loss of heterozygosity affects fitness for different types of apomixis and automixis, and I show that soft selection significantly reduces loss of heterozygosity, hence increases fitness, in apomixis with suppression of the first meiotic division and in automixis with central fusion, the most common types of asexual reproduction. Therefore, if sexual reproduction evolved to preserve heterozygosity, soft selection should be associated with these types of asexual reproduction. I discuss the evidence for this prediction and how this and other observations on the distribution of different types of asexual reproduction in nature is consistent with the heterozygosity hypothesis.     

An estimate of unique DNA sequence heterozygosity in the human genome

Summary Patients with recessive X-linked ichthyosis Patients with recessive X-linked ichthyosis (RXLI), one hereditary form of scaly skin, lack activity of the enzyme steroid sulfatase in all tissues studied. To investigate the molecular defect underlying the lack of enzyme activity, we prepared antisera against normal enzyme by injecting normal placental microsomal suspensions or partially purified steroid sulfatase into rabbits. Antibody activity was assessed by immunoprecipitation of detergent solubilized steroid sulfatase. In addition, we prepared rabbit antisera against RXLI placental microsomal suspensions. To detect immunologically cross-reactive material in patients' placentas, extracts were studied by immunoblot techniques and by competition with normal enzyme for antibody binding. Patients' extracts did not contain immunoreactive material co-migrating on electrophoresis with purified enzyme nor did they inhibit immunoprecipitation of normal enzyme. Sera from rabbits immunized with RXLI placental microsomes contain no antibodies to normal steroid sulfatase, as judged by their failure to immunoprecipitate normal enzyme or to react with normal steroid sulfatase on immunoblot. Thus the mutation in RXLI appears to reduce steroid sulfatase enzyme protein as well as enzyme activity. Portions of this material have appeared in abstract form in Clinical Research 31:564A, 1983 and 32:138A, 1984     

High-resolution detection of loss of heterozygosity of dinucleotide microsatellite markers

Dinucleotide microsatellite markers are frequently investigated to study inheritance, genetic stability, and allele frequency distribution in a wide variety of genetic disorders. Previous studies have encountered significant problems regarding resolution and detection of dinucleotide, microsatellites. In this study, a useful method to investigate loss of heterozygosity (LOH) of dinucleotide microsatellite markers is described that involves the use of nondenaturing (Spreadex) submerged gel electrophoresis and SYBR Green I nucleic acid staining. This method omits the gel casting step and the use of hazardous radioactive materials frequently used in many microsatellite studies that employ polyacrylamide gel nucleic acid denaturation analysis. Using this method, 62 patients' paired tumor and normal samples were investigated to detect allele deletions in a region of chromosome 7q31.1, which is believed to harbor a tumor suppressor gene. Interpretable results were obtained in all cases. These results were compared to those attained using ABI Prism Genetic Analyzer 310 and Gene-Scan. There were no discrepancies in results obtained between the two assays. The Spreadex system is cheap, does not require larger equipment costs, and may prove to be a useful system for high-throughput investigation of microsatellites. It may have diagnostic significance and also prove useful if applied to population-based genomic screening and linkage analysis.