Integrating copy number polymorphisms into array CGH analysis using a robust HMM

MOTIVATION: Array comparative genomic hybridization (aCGH) is a pervasive technique used to identify chromosomal aberrations in human diseases, including cancer. Aberrations are defined as regions of increased or decreased DNA copy number, relative to a normal sample. Accurately identifying the locations of these aberrations has many important medical applications. Unfortunately, the observed copy number changes are often corrupted by various sources of noise, making the boundaries hard to detect. One popular current technique uses hidden Markov models (HMMs) to divide the signal into regions of constant copy number called segments; a subsequent classification phase labels each segment as a gain, a loss or neutral. Unfortunately, standard HMMs are sensitive to outliers, causing over-segmentation, where segments erroneously span very short regions. RESULTS: We propose a simple modification that makes the HMM robust to such outliers. More importantly, this modification allows us to exploit prior knowledge about the likely location of "outliers", which are often due to copy number polymorphisms (CNPs). By "explaining away" these outliers with prior knowledge about the locations of CNPs, we can focus attention on the more clinically relevant aberrated regions. We show significant improvements over the current state of the art technique (DNAcopy with MergeLevels) on previously published data from mantle cell lymphoma cell lines, and on published benchmark synthetic data augmented with outliers. AVAILABILITY: Source code written in Matlab is available from http://www.cs.ubc.ca/~sshah/acgh.     

Rapid high resolution western blotting: from gel to image in a single day.

A streamlined protocol is described that allows high sensitivity antigen detection by Western blotting in a single day. The choice of membrane blotting matrix, as well as blocking reagents, has been optimized in order to allow rapid development of the blot with chemiluminescent reagents. The entire process, from gel to blot to a permanent, hard copy image on x-ray film, can be accomplished within six hours.     

Cbf14 copy number variation in the A, B, and D genomes of diploid and polyploid wheat

Freezing tolerance and winter hardiness are complex traits. In the Triticeae, two loci on the group 5 chromosome homoeologs are repeatedly identified as having major effects on these traits. Recently, we found that segments of the genomic region at one of these loci, Frost resistance-2 (Fr-2) is copy number variable in barley. Freezing-tolerant winter-hardy genotypes have greater tandem copy numbers of the genomic region encompassing the C-repeat binding factor genes Cbf2A and Cbf4B at Fr-H2 than the less freezing-tolerant nonwinter-hardy genotypes. Here we report that in wheat the Cbf14 gene at Fr-2 is copy number variable. Using DNA blot hybridizations, we estimated copy numbers of Cbf14 across the different genomes of diploid and polyploid wheat. Copy numbers of Cbf14 are lower in the B genome than in the A and D genomes across all ploidy levels. Among hexaploid red wheats, winter genotypes harbor greater Cbf14 copy numbers than spring genotypes. Cbf14 copy numbers also vary across the red winter wheats such that hard wheats harbor greater copy numbers than soft wheats. Analysis of hexaploid wheat chromosome 5 substitution lines indicates that Cbf14 copy numbers in the introgressions are stable in the different backgrounds. Taken together our data suggest that higher copy number states existed in the diploid wild ancestors prior to the polyploidization events and that the loss of Cbf14 copies occurred in the cultivated germplasm.     

On the origins of tandemly repeated genes: Does histone gene copy number inDrosophila reflect chromosomal location?

Widely regarded beliefs aboutDrosophila histone gene copy numbers and developmental requirements have been generalized from fairly limited data since studies on histone gene arrangements and copy numbers have been largely confined to a single species,D. melanogaster. Histone gene copy numbers and chromosomal locations were examined in three species:D. melanogaster, D. hydei andD. hawaiiensis. Quantitative whole genome blot analysis of DNA from diploid tissues revealed a tenfold variability in histone gene copy numbers for these three species. In situ hybridization to polytene chromosomes showed that the histone DNA (hDNA) chromosomal location is different in all three species. These observations lead us to propose a relationship between histone gene reiteration and chromosomal position.     

Proliferation and copy number variation of BEL-like long terminal repeat retrotransposons within the Diabrotica virgifera virgifera genome

The proliferation of retrotransposons within a genome can contribute to increased size and affect the function of eukaryotic genes. BEL/Pao-like long-terminal repeat (LTR) retrotransposons were annotated from the highly adaptable insect species Diabrotica virgifera virgifera, the Western corn rootworm, using survey sequences from bacterial artificial chromosome (BAC) inserts and contigs derived from a low coverage next-generation genome sequence assembly. Eleven unique D. v. virgifera BEL elements were identified that contained full-length gag–pol coding sequences, whereas 88 different partial coding regions were characterized from partially assembled elements. Estimated genome copy number for full and partial BEL-like elements ranged from ~ 8 to 1582 among individual contigs using a normalized depth of coverage (DOC) among Illumina HiSeq reads (total genome copy number ~ 8821). BEL element copy number was correlated among different D. v. virgifera populations (R² = 0.9846), but individual element numbers varied ≤ 1.68-fold and the total number varied by ~ 527 copies. These data indicate that BEL element proliferation likely contributed to a large genome size, and suggest that differences in copy number are a source of genetic variability among D. v. virgifera.     

Removal of fiber from vines by solid state fermentation/enzymatic degradation: a comparison of flax and kudzu retting

Kinetic data describing the decomposition of the outer sheath of kudzu vines (undergoing a solid fermentation process in a glass beaker of soil) have been analyzed to determine the two constants, K(m)/S(o) and V/S(o), where S(o) is the initial substrate concentration, K(m) the Michaelis constant, and V the maximum product rate. The kinetic data are expressed by a simple time-varying desheathing index, obtained from the number of spatula scrapings required to reach the desired hard cellulosic fibers (vascular bundles) of the plant. A simple relationship between the desheathing index, D.I. and the dimensionless product concentration, P/S(o), is proposed to relate the D.I. data and enzyme kinetic concentration data. Thus, the Michaelis-Menten enzyme kinetic parameters can be estimated from easily obtained physical (non-chemical data; the D.I.(t) measurements). This low energy process for processing vines into valuable fibers is similar to the traditional microbial method for recovering flax fibers for linen cloth, by retting of the flax plant vines, except there is no unbound liquid water present in the soil medium.     

On the use of a general-purpose CAD program for molecular graphics

It is suggested that general-purpose computer-aided design software can be of good use in producing high-quality molecular graphics. Such software is often widely available. Besides being an alternative to existing programs, hard copy creation is facilitated by this approach.     

A fluorescence cytophotometer operated under computer control for multi-parameter cell analysis

Summary A new type of fluorescence cytophotometer has been developed for multi-parameter cell analysis (Olympus BH2-QRFL). For multi-color fluorescence cytophotometry, this instrument is equipped with four sets of interchangeable filters, each consisting of an excitation filter, a dichroic mirror with a barrier filter, and a measuring filter. For permitting automatic operation of the filter sets, the cytophotometer is connected on line with a personal computer (HP 85F). A desired sequence of filter sets can be memorized in the software and multiple cellular constituents can be rapidly and consecutively determined on a single cell basis. All data are stored in the same computer and can be retrieved for further statistical analysis and display either in tabular form, or as histogram, correlation histograms, two-dimensional scatter plot, or two-dimensional frequency distribution histogram, on the CRT (cathode ray tube) with simultaneous hard copy. As an example of multiparameter cell analysis, combined protein and DNA measurements were performed on normal, border-line, and cancerous gynecological cytology specimens by using the ninhydrin-Schiff and Feulgen techniques.     

ArrayCyGHt: a web application for analysis and visualization of array-CGH data

ArrayCyGHt is a web-based application tool for analysis and visualization of microarray-comparative genomic hybridization (array-CGH) data. Full process of array-CGH data analysis, from normalization of raw data to the final visualization of copy number gain or loss, can be straightforwardly achieved on this arrayCyGHt system without the use of any further software. ArrayCyGHt, therefore, provides an easy and fast tool for the analysis of copy number aberrations in any kinds of data format. AVAILABILITY: ArrayCyGHt can be accessed at http://genomics.catholic.ac.kr/arrayCGH/     

Improvement of peanut (Arachis hypogaea L.) transformation efficiency and determination of transgene copy number by relative quantitative real-time PCR

The biolistic method is reliable for delivering genes of interest into various species, but low transformation efficiency can be a limiting factor in its application. To test various conditions that could improve peanut transformation via particle bombardment, embryogenic tissues of the peanut cultivar Georgia Green were co-bombarded with two plasmids: one containing a green fluorescent protein gene and one containing a gene of interest plus a selectable marker. The fluorescence in bombarded embryogenic tissues was measured to evaluate transformation efficiency. A 4.6-fold improvement of transformation efficiency was achieved in stably transformed peanut lines by introducing protamine instead of conventional spermidine in a bombardment mixture with 70 ng/shot plasmid DNA and 50 μg/shot gold. Unexpectedly, the reduction of plasmid DNA from 700 to 70 ng/shot produced transgenic lines with significantly increased numbers of transgene copies. To determine the transgene copy number during plantlet regeneration, relative quantitative real-time polymerase chain reaction (qPCR) was established using fluorescently labeled universal library probes. A correlation of 95% was found for estimation of copy number between Southern blot and qPCR data. Given its speed and high-throughput nature, qPCR can be employed as an effective screening tool to separate high copy number events from low copy number events as early as the shoot formation stage of regeneration.     

Mechanisms regulating the copy numbers of six LTR retrotransposons in the genome of Drosophila melanogaster

There has been debate over the mechanisms that control the copy number of transposable elements in the genome of Drosophila melanogaster. Target sites in D. melanogaster populations are occupied at low frequencies, suggesting that there is some form of selection acting against transposable elements. Three main theories have been proposed to explain how selection acts against transposable elements: insertions of a copy of a transposable element are selected against; chromosomal rearrangements caused by ectopic exchange between element copies are selected against; or the process of transposition itself is selected against. The three theories give different predictions for the pattern of transposable element insertions in the chromosomes of D. melanogaster. We analysed the abundance of six LTR (long terminal repeat) retrotransposons on the X and fourth chromosomes of multiple strains of D. melanogaster, which we compare with the predictions of each theory. The data suggest that no one theory can account for the insertion patterns of all six retrotransposons. Comparing our results with earlier work using these transposable element families, we find a significant correlation between studies in the particular model of copy number regulation supported by the proportion of elements on the X for the different transposable element families. This suggests that different retrotransposon families are regulated by different mechanisms.     

The population genetics of the origin and divergence of the Drosophila simulans complex species

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.     

Control-free calling of copy number alterations in deep-sequencing data using GC-content normalization

SUMMARY: We present a tool for control-free copy number alteration (CNA) detection using deep-sequencing data, particularly useful for cancer studies. The tool deals with two frequent problems in the analysis of cancer deep-sequencing data: absence of control sample and possible polyploidy of cancer cells. FREEC (control-FREE Copy number caller) automatically normalizes and segments copy number profiles (CNPs) and calls CNAs. If ploidy is known, FREEC assigns absolute copy number to each predicted CNA. To normalize raw CNPs, the user can provide a control dataset if available; otherwise GC content is used. We demonstrate that for Illumina single-end, mate-pair or paired-end sequencing, GC-contentr normalization provides smooth profiles that can be further segmented and analyzed in order to predict CNAs. AVAILABILITY: Source code and sample data are available at http://bioinfo-out.curie.fr/projects/freec/.     

Propagation of an amplifiable recombinant plasmid in Saccharomyces cerevisiae: flow cytometry studies and segregated modeling

Efficient expression of a foreign protein product by the yeastSaccharomyces cerevisiaerequires a stable recombinant vector present at a high number of copies per cell. A conditional centromere yeast plasmid was constructed which can be amplified to high copy number by a process of unequal partitioning at cell division, followed by selection for increased copy number. However, in the absence of selection pressure for plasmid amplification, copy number rapidly drops from 25 plasmids/cell to 6 plasmids/cell in less than 10 generations of growth. Copy number subsequently decreases from 6 plasmids/cell to 2 plasmids/cell over a span of 50 generations. A combination of flow cytometric measurement of copy number distributions and segregated mathematical modeling were applied to test the predictions of a conceptual model of conditional centromereplasmid propagation. Measured distributions of plasmid content displayed a significant subpopulation of cells with a copy number of 4-6, evenin a population whose mean copy number was 13.5. This type of copy number distribution was reproduced by a mathematical model which assumes that amaximum of 4-6 centromere plasmids per cell can be stably partitionedat cell division. The model also reproduces the observed biphasic kinetics of plasmid number instability. The agreement between simulation and experimental results provides support for the proposed model and demonstrates the utility of the flow cytometry/segregated modeling approach for the study of multicopy recombinant vector propagation.     

A Novel Simple Method for Determining CYP2D6 Gene Copy Number and Identifying Allele(s) with Duplication/Multiplication

Background

Cytochrome P450 2D6 (CYP2D6) gene duplication and multiplication can result in ultrarapid drug metabolism and therapeutic failure or excessive response in patients. Long range polymerase chain reaction (PCR), restriction fragment length polymorphism (RFLP) and sequencing are usually used for genotyping CYP2D6 duplication/multiplications and identification, but are labor intensive, time consuming, and costly.

Methods

We developed a simple allele quantification-based Pyrosequencing genotyping method that facilitates CYP2D6 copy number variation (CNV) genotyping while also identifying allele-specific CYP2D6 CNV in heterozygous samples. Most routine assays do not identify the allele containing a CNV. A total of 237 clinical and Coriell DNA samples with different known CYP2D6 gene copy numbers were genotyped for CYP2D6 *2, *3, *4, *6, *10, *17, *41 polymorphisms and CNV determination.

Results

The CYP2D6 gene allele quantification/identification were determined simultaneously with CYP2D6*2, *3, *4, *6, *10, *17, *41 genotyping. We determined the exact CYP2D6 gene copy number, identified which allele had the duplication or multiplication, and assigned the correct phenotype and activity score for all samples.

Conclusions

Our method can efficiently identify the duplicated CYP2D6 allele in heterozygous samples, determine its copy number in a fraction of time compared to conventional methods and prevent incorrect ultrarapid phenotype calls. It also greatly reduces the cost, effort and time associated with CYP2D6 CNV genotyping.     

Analysis of extracellular potentials using an IBM PC

1. A system has been developed for using IBM PC-compatible computers in combination with a Grafitek Data Logging Interface to record spike trains on magentic discs for later analysis. 2. The times and amplitudes of spikes detected on two input channels are recorded, together with a third channel containing information on computer-generated stimuli and keyboard-activated event markers. In excess of 50,000 spikes can be recorded with a computer having 640 k of Random Access Memory. 3. The recorded spike trains can be reconstructed on the computer monitor and keyboard-controlled window discriminators can be used to select the spikes for analysis by amplitude. 4. The same recorded data can be analysed to produce displays of spike count against time, amplitude histograms, inter-spike interval histograms, peri-stimulus time histograms(PSTH), raster displays and auto- and cross-correlations between activity on the two channels. Each spike is identified by number, allowing easy location of the start and finish of the section of data to be analysed, and the PSTH, raster and correlation analyses allow pretriggering to investigate event occurring before stimulation. 5. The axes of the displays histograms can be adjusted to produce optimum displays, and hard copy can be produced on dot matrix printers or digital plotters. 6. Quantitative analysis enables comparison between different recordings and treatments.     

Coat colours in the Massese sheep breed are associated with mutations in the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes

Massese is an Italian dairy sheep breed characterized by animals with black skin and horns and black or apparent grey hairs. Owing to the presence of these two coat colour types, this breed can be considered an interesting model to evaluate the effects of coat colour gene polymorphisms on this phenotypic trait. Two main loci have been already shown to affect coat colour in sheep: Agouti and Extension coding for the agouti signalling protein (ASIP) and melanocortin 1 receptor (MC1R) genes, respectively. The Agouti locus is affected by a large duplication including the ASIP gene that may determine the Agouti white and tan allele (A(Wt)). Other disrupting or partially inactivating mutations have been identified in exon 2 (a deletion of 5 bp, D(5); and a deletion of 9 bp, D(9)) and in exon 4 (g.5172T>A, p.C126S) of the ASIP gene. Three missense mutations in the sheep MC1R gene cause the dominant black E(D) allele (p.M73K and p.D121N) and the putative recessive e allele (p.R67C). Here, we analysed these ASIP and MC1R mutations in 161 Massese sheep collected from four flocks. The presence of one duplicated copy allele including the ASIP gene was associated with grey coat colour (P = 9.4E-30). Almost all animals with a duplicated copy allele (37 out of 41) showed uniform apparent grey hair and almost all animals without a duplicated allele (117 out of 120) were completely black. Different forms of duplicated alleles were identified in Massese sheep including, in almost all cases, copies with exon 2 disrupting or partially inactivating mutations making these alleles different from the A(Wt) allele. A few exceptions were observed in the association between ASIP polymorphisms and coat colour: three grey sheep did not carry any duplicated copy allele and four black animals carried a duplicated copy allele. Of the latter four sheep, two carried the E(D) allele of the MC1R gene that may be the cause of their black coat colour. The coat colour of all other black animals may be determined by non-functional ASIP alleles (non-agouti alleles, A(a)) and in a few cases by the E(D) Extension allele. At least three frequent ASIP haplotypes ([D(5):g.5172T], [N:g.5172A] and [D(5):g.5172A]) were detected (organized into six different diplotypes). In conclusion, the results indicated that coat colours in the Massese sheep breed are mainly derived by combining ASIP and MC1R mutations.