A MULTILEVEL MODEL TO ADDRESS BATCH EFFECTS IN COPY NUMBER ESTIMATION USING SNP ARRAYS

Submicroscopic changes in chromosomal DNA copy number dosage are common and have been implicated in many heritable diseases and cancers. Recent high-throughput technologies have a resolution that permits the detection of segmental changes in DNA copy number that span thousands of basepairs across the genome. Genome-wide association studies (GWAS) may simultaneously screen for copy number-phenotype and SNP-phenotype associations as part of the analytic strategy. However, genome-wide array analyses are particularly susceptible to batch effects as the logistics of preparing DNA and processing thousands of arrays often involves multiple laboratories and technicians, or changes over calendar time to the reagents and laboratory equipment. Failure to adjust for batch effects can lead to incorrect inference and requires inefficient post-hoc quality control procedures that exclude regions that are associated with batch. Our work extends previous model-based approaches for copy number estimation by explicitly modeling batch effects and using shrinkage to improve locus-specific estimates of copy number uncertainty. Key features of this approach include the use of diallelic genotype calls from experimental data to estimate batch- and locus-specific parameters of background and signal without the requirement of training data. We illustrate these ideas using a study of bipolar disease and a study of chromosome 21 trisomy. The former has batch effects that dominate much of the observed variation in quantile-normalized intensities, while the latter illustrates the robustness of our approach to datasets where as many as 25% of the samples have altered copy number. Locus-specific estimates of copy number can be plotted on the copy-number scale to investigate mosaicism and guide the choice of appropriate downstream approaches for smoothing the copy number as a function of physical position. The software is open source and implemented in the R package CRLMM available at Bioconductor (http:www.bioconductor.org).

A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

Losses of chromosome arms 4q, 8p, 13q and gain of 8q are correlated with increasing chromosomal instability in hepatocellular carcinoma

OBJECTIVE: Chromosomal instability is a key feature in hepatocellular carcinoma (HCC). Array comparative genomic hybridization (aCGH) revealed recurring structural aberrations, whereas fluorescence in situ hybridization (FISH) indicated an increasing number of numerical aberrations in dedifferentiating HCC. Therefore, we examined whether there was a correlation between structural and numerical aberrations of chromosomal instability in HCC. METHODS AND RESULTS: 27 HCC (5 well, 10 moderately, 12 lower differentiated) already cytogenetically characterized by aCGH were analyzed. FISH analysis using probes for chromosomes 1, 3, 7, 8 and 17 revealed 1.46-4.24 signals/nucleus, which correlated with the histological grade (well vs. moderately,p < 0.0003; moderately vs. lower, p < 0.004). The number of chromosomes to each other was stable with exceptions only seen for chromosome 8. Loss of 4q and 13q, respectively, were correlated with the number of aberrations detected by aCGH (p < 0.001, p < 0.005; Mann-Whitney test). Loss of 4q and gain of 8q were correlated with an increasing number of numerical aberrations detected by FISH (p < 0.020, p < 0.031). Loss of 8p was correlated with the number of structural imbalances seen in aCGH (p < 0.048), but not with the number of numerical changes seen in FISH. CONCLUSION: We found that losses of 4q, 8p and 13q were closely correlated with an increasing number of aberrations detected by aCGH, whereas a loss of 4q and a gain of 8q were also observed in the context of polyploidization, the cytogenetic correlate of morphological dedifferentiation.

Molecular phylogenetic analyses identify Alpine differentiation and dysploid chromosome number changes as major forces for the evolution of the European endemic Phyteuma (Campanulaceae)

Phyteuma is a chromosomally and ecologically diverse vascular plant genus and constitutes an excellent system for studying both the role of chromosomal change for species diversification and the evolution of high-mountain biota. This kind of research is, however, hampered by the lack of a sound phylogenetic framework exacerbated by the notoriously low predictive power of traditional taxonomy with respect to phylogenetic relationships in Campanulaceae. Based on a comprehensive taxon sampling and analyses of nuclear and plastid sequence and AFLP fingerprint data, Phyteuma is confirmed as a monophyletic group sister to the monotypic Physoplexis, which is in line with their peculiar flower morphologies. Within Phyteuma two clades, largely corresponding to previously recognized sections, are consistently found. The traditional circumscription of taxonomic series is largely rejected. Whereas distinctness of the currently recognized species is mostly corroborated, some interspecific relationships remain ambiguous due to incongruences between nuclear and plastid data. Major forces for diversification and evolution of Phyteuma are descending dysploidy (i.e., a decrease in chromosome base number) as well as allopatric and ecological differentiation within the Alps, the genus' center of species diversity. (C) 2013 Elsevier Inc. All rights reserved.

Chromosomal organization and evolutionary conservation of the human chromosome 19q linkage group containing three DNA repair genes

A series of human-rodent somatic cell hybrids were investigated by Southern blot analysis for the presence or absence of twenty-six molecular markers and three isozyme loci from human chromosome 19. Based on the co-retention of these markers in the various independent hybrid clones containing portions of human chromosome 19 and on pulsed field mapping, chromosome 19 is divided into twenty ordered regions. The most likely marker order for the chromosome is: (LDLR, C3)-(cen-MANNB)-D19S7-PEPD-D19S9-GPI-TGF$ \beta$-(CYP2A, NCA, CGM2, BCKAD)-PSG1a-(D19S8, XRCC1)-(D19S19, ATP1A3)-(D19S37, APOC2)-CKMM-ERCC2-ERCC1-(D19S62, D19S51)-D19S6-D19S50-D19S22-(CGB, FTL)-qter.^ The region of 19q between the proximal marker D19S7 and the distal gene coding for the beta subunit of chorionic gonadotropin (CGB) is about 37 Mb in size and covers about 37 cM genetic distance. The ration of genetic to physical distance on 19q is therefore very close to the genomic average OF 1 cM/Mb. Estimates of physical distances for intervals between chromosome 19 markers were calculated using a mapping function which estimates distances based on the number of breaks in hybrid clone panels. The consensus genetic distances between individual markers (established at HBM10) were compared to these estimates of physical distances. The close agreement between the two estimates suggested that spontaneously broken hybrids are as appropriate for this type of study as radiation hybrids.^ All three DNA repair genes located on chromosome 19 were found to have homologues on Chinese hamster chromosome 9, which is hemizygous in CHO cells, providing an explanation for the apparent ease with which mutations at these loci were identified in CHO cells. Homologues of CKMM and TGF$\beta$ (from human chromosome 19q) and a mini-satellite DNA specific to the distal region of human chromosome 19q were also mapped to Chinese hamster 9. Markers from 19p did not map to this hamster chromosome. Thus the q-arm of chromosome 19, at least between the genes PEPD and ERCC1, appears to be a linkage group which is conserved intact between humans and Chinese hamsters. ^

Melanotic tumors of the nervous system are characterized by distinct mutational, chromosomal and epigenomic profiles.

Melanotic tumors of the nervous system show overlapping histological characteristics but differ substantially in their biological behavior. In order to achieve a better delineation of such tumors, we performed an in-depth molecular characterization. Eighteen melanocytomas, 12 melanomas, and 14 melanotic and 14 conventional schwannomas (control group) were investigated for methylome patterns (450k array), gene mutations associated with melanotic tumors and copy number variants (CNVs). The methylome fingerprints assigned tumors to entity-specific groups. Methylation groups also showed a substantial overlap with histology-based diagnosis suggesting that they represent true biological entities. On the molecular level, melanotic schwannomas were characterized by a complex karyotype with recurrent monosomy of chromosome 22q and variable whole chromosomal gains and recurrent losses commonly involving chromosomes 1, 17p and 21. Melanocytomas carried GNAQ/11 mutations and presented with CNV involving chromosomes 3 and 6. Melanomas were frequently mutated in the TERT promoter, harbored additional oncogene mutations and showed recurrent chromosomal losses involving chromosomes 9, 10 and 6q, as well as gains of 22q. Together, melanotic nervous system tumors have several distinct mutational and chromosomal alterations and can reliably be distinguished by methylome profiling.

Shorter telomeres correlate with an increase in the number of uniparental disomies in patients with chronic lymphocytic leukemia.

This study investigated the correlation of the extent of chromosomal aberrations including uniparental disomies (UPDs) by SNP-chip analysis and FISH to telomere length in 46 patients with CLL. CLL harboring high risk aberrations, i.e. deletions of 11q22-23 or 17p13, had significantly shorter telomeres (higher ΔTL) compared to patients with CLL without such abnormalities. Patients with high chromosomal aberration rates had a worse overall survival compared to cases with lower aberration rates. Interestingly, however, an increase was found in the number of UPDs with shorter telomeres. These findings support the idea that telomeres in CLL cells play a role in the overall chromosome stability and could be involved in the occurrence of UPDs.

Chromosomal breakage-fusion-bridge events cause genetic intratumor heterogeneity

It has long been known that rearrangements of chromosomes through breakage-fusion-bridge (BFB) cycles may cause variability of phenotypic and genetic traits within a cell population. Because intercellular heterogeneity is often found in neoplastic tissues, we investigated the occurrence of BFB events in human solid tumors. Evidence of frequent BFB events was found in malignancies that showed unspecific chromosome aberrations, including ring chromosomes, dicentric chromosomes, and telomeric associations, as well as extensive intratumor heterogeneity in the pattern of structural changes but not in tumors with tumor-specific aberrations and low variability. Fluorescence in situ hybridization analysis demonstrated that chromosomes participating in anaphase bridge formation were involved in a significantly higher number of structural aberrations than other chromosomes. Tumors with BFB events showed a decreased elimination rate of unstable chromosome aberrations after irradiation compared with normal cells and other tumor cells. This result suggests that a combination of mitotically unstable chromosomes and an elevated tolerance to chromosomal damage leads to constant genomic reorganization in many malignancies, thereby providing a flexible genetic system for clonal evolution and progression.

Heat effects on DNA repair after ionising radiation: hyperthermia commonly increases the number of non-repaired double-strand breaks and structural rearrangements

After ionising radiation double-strand breaks (dsb) are lethal if not repaired or misrepaired. Cell killing is greatly enhanced by hyperthermia and it is questioned here whether heat not only affects dsb repair capacity but also fidelity in a chromosomal context. dsb repair experiments were designed so as to mainly score non-homologous end joining, while homologous recombination was largely precluded. Human male G0 fibroblasts were either preheated (45°C, 20 min) or not before X-irradiation. dsb induction and repair were measured by conventional gel electrophoresis and an assay combining restriction digestion using a rare cutting enzyme (NotI) and Southern hybridisation, which detects large chromosomal rearrangements (>100 kb). dsb induction rate in an X-chromosomal NotI fragment was 4.8 × 10–3 dsb/Gy/Mb. Similar values were found for the genome overall and also when cells were preheated. After 50 Gy, fibroblasts were competent to largely restore the original restriction fragment size. Five per cent of dsb remained non-rejoined and 14% were misrejoined. Correct restitution of restriction fragments occurred preferably during the first hour but continued at a slow rate for 12–16 h. In addition, dsb appeared to misrejoin throughout the entire repair period. After hyperthermia the fractions of non-rejoined and misrejoined dsb were similarly increased to 13 and 51%, respectively. It is suggested that heat increases the probability of dsb being incorrectly rejoined but it is not likely to interfere with one dsb repair pathway in particular.

An Escherichia coli chromosomal "addiction module" regulated by guanosine [corrected] 3',5'-bispyrophosphate: a model for programmed bacterial cell death.

"Addiction modules" consist of two genes. In most of them the product of one is long lived and toxic while the product of the second is short lived and antagonizes the toxic effect; so far, they have been described mainly in a number of prokaryotic extrachromosomal elements responsible for the postsegregational killing effect. Here we show that the chromosomal genes mazE and mazF, located in the Escherichia coli rel operon, have all of the properties required for an addiction module. Furthermore, the expression of mazEF is regulated by the cellular level of guanosine [corrected] 3',5'-bispyrophosphate, the product of the RelA protein under amino acid starvation. These properties suggest that the mazEF system may be responsible for programmed cell death in E. coli and thus may have a role in the physiology of starvation.

Chromosomal localization of the mammalian peptide-methionine sulfoxide reductase gene and its differential expression in various tissues.

Peptide methionine sulfoxide reductase (MsrA; EC 1.8.4.6) is a ubiquitous protein that can reduce methionine sulfoxide residues in proteins as well as in a large number of methyl sulfoxide compounds. The expression of MsrA in various rat tissues was determined by using immunocytochemical staining. Although the protein was found in all tissues examined, it was specifically localized to renal medulla and retinal pigmented epithelial cells, and it was prominent in neurons and throughout the nervous system. In addition, blood and alveolar macrophages showed high expression of the enzyme. The msrA gene was mapped to the central region of mouse chromosome 14, in a region of homology with human chromosomes 13 and 8p21.

Influence of alkyltransferase activity and chromosomal locus on mutational hotspots in Chinese hamster ovary cells.

High-density mutational spectra have been established for exon 3 of the gene encoding adenine phosphoribosyltransferase (APRT) of the Chinese hamster ovary (CHO) cell line derivative D422 and closely related and/or modified lines by using the mutagen ethyl methanesulfonate (EMS). The total number of selectable sites (GC-->AT transitions yielding a selectable APRT- phenotype) was estimated at 31 based on our own accumulated data base of 136 sequenced exon 3 mutations and on literature reports. D422 and two other APRT hemizygous lines each yielded very similar spectra and showed two populations of mutable sites: (i) 24 "baseline" sites that followed the Poisson distribution and therefore were equally susceptible to mutation and (ii) two hotspots, one comprising a cluster at nucleotides 1293-1309 and the other at nucleotide 1365. Collectively, the latter sites were about 10-fold more frequently mutated than the others. CHO cells are mer- as they lack the repair enzyme O6-methylguanidine methyltransferase (EC 2.1.1.63). In modified repair-proficient CHO cells, the distribution of mutations among all of the 31 sites was random, with only 3 of the 19 GC-->AT transitions in the above hotspots. To determine whether the distribution was locus-dependent, two independent lines carrying single copies of transfected APRT genes were generated from a derivative of D422 carrying a deletion in the endogenous APRT gene. Nucleotides 1293-1309 were again no longer preferentially mutated, but the site at nucleotide 1365 was still a hotspot. We conclude that mutational spectra in mer- cells are at least in part locus dependent and that some sequences are particularly susceptible to EMS mutagenesis and perhaps also to methyltransferase repair.

Alignment of a 1.2-Mb chromosomal region from three strains of Rhodobacter capsulatus reveals a significantly mosaic structure.

High-resolution physical maps of the genomes of three Rhodobacter capsulatus strains, derived from ordered cosmid libraries, were aligned. The 1.2-Mb segment of the SB1003 genome studied here is adjacent to a 1-Mb region analyzed previously [Fonstein, M., Nikolskaya, T. & Haselkorn, H. (1995) J. Bacteriol. 177, 2368-2372]. Probes derived from the ordered cosmid set of R. capsulatus SB1003 were used to link cosmids from the St. Louis and 2.3.1 strain libraries. Cosmids selected this way did not merge into a single contig but formed several unlinked groups. EcoRV restriction maps of the ordered cosmids were then constructed using lambda terminase and fused to derive fragments of the chromosomal map. In order to link these fragments, their ends were transcribed to produce secondary probes for hybridization to gridded cosmid libraries of the same strains. This linking reduced the number of subcontigs to three for the St. Louis strain and one for the 2.3.1 strain. Hybridization of the same probes back to the ordered cosmid set of SB1003 positioned the subcontigs on the high-resolution physical map of SB1003. The final alignment of the restriction maps shows numerous large and small translocations in this 1.2-Mb chromosomal region of the three Rhodobacter strains. In addition, the chromosomes of the three strains, whose fine-structure maps can now be compared over 2.2 Mb, are seen to contain regions of 15-80 kb in which restriction sites are highly polymorphic, interspersed among regions in which the positions of restriction sites are highly conserved.

Petroleomics By Electrospray Ionization Ft-icr Mass Spectrometry Coupled To Partial Least Squares With Variable Selection Methods: Prediction Of The Total Acid Number Of Crude Oils.

Negative-ion mode electrospray ionization, ESI(-), with Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) was coupled to a Partial Least Squares (PLS) regression and variable selection methods to estimate the total acid number (TAN) of Brazilian crude oil samples. Generally, ESI(-)-FT-ICR mass spectra present a power of resolution of ca. 500,000 and a mass accuracy less than 1 ppm, producing a data matrix containing over 5700 variables per sample. These variables correspond to heteroatom-containing species detected as deprotonated molecules, [M - H](-) ions, which are identified primarily as naphthenic acids, phenols and carbazole analog species. The TAN values for all samples ranged from 0.06 to 3.61 mg of KOH g(-1). To facilitate the spectral interpretation, three methods of variable selection were studied: variable importance in the projection (VIP), interval partial least squares (iPLS) and elimination of uninformative variables (UVE). The UVE method seems to be more appropriate for selecting important variables, reducing the dimension of the variables to 183 and producing a root mean square error of prediction of 0.32 mg of KOH g(-1). By reducing the size of the data, it was possible to relate the selected variables with their corresponding molecular formulas, thus identifying the main chemical species responsible for the TAN values.

The Basic Reproduction Number Obtained From Jacobian And Next Generation Matrices - A Case Study Of Dengue Transmission Modelling.

The basic reproduction number is a key parameter in mathematical modelling of transmissible diseases. From the stability analysis of the disease free equilibrium, by applying Routh-Hurwitz criteria, a threshold is obtained, which is called the basic reproduction number. However, the application of spectral radius theory on the next generation matrix provides a different expression for the basic reproduction number, that is, the square root of the previously found formula. If the spectral radius of the next generation matrix is defined as the geometric mean of partial reproduction numbers, however the product of these partial numbers is the basic reproduction number, then both methods provide the same expression. In order to show this statement, dengue transmission modelling incorporating or not the transovarian transmission is considered as a case study. Also tuberculosis transmission and sexually transmitted infection modellings are taken as further examples.