Skittle: A 2-Dimensional Genome Visualization Tool

<p>Abstract</p> <p>Background</p> <p>It is increasingly evident that there are multiple and overlapping patterns within the genome, and that these patterns contain different types of information - regarding both genome function and genome history. In order to discover additional genomic patterns which may have biological significance, novel strategies are required. To partially address this need, we introduce a new data visualization tool entitled Skittle.</p> <p>Results</p> <p>This program first creates a 2-dimensional nucleotide display by assigning four colors to the four nucleotides, and then text-wraps to a user adjustable width. This nucleotide display is accompanied by a "repeat map" which comprehensively displays all local repeating units, based upon analysis of all possible local alignments. Skittle includes a smooth-zooming interface which allows the user to analyze genomic patterns at any scale.</p> <p>Skittle is especially useful in identifying and analyzing tandem repeats, including repeats not normally detectable by other methods. However, Skittle is also more generally useful for analysis of any genomic data, allowing users to correlate published annotations and observable visual patterns, and allowing for sequence and construct quality control.</p> <p>Conclusions</p> <p>Preliminary observations using Skittle reveal intriguing genomic patterns not otherwise obvious, including structured variations inside tandem repeats. The striking visual patterns revealed by Skittle appear to be useful for hypothesis development, and have already led the authors to theorize that imperfect tandem repeats could act as information carriers, and may form tertiary structures within the interphase nucleus.</p

Refinement of Bos taurus sequence assembly based on BAC-FISH experiments

<p>Abstract</p> <p>Background</p> <p>The sequencing of the cow genome was recently published (Btau_4.0 assembly). A second, alternate cow genome assembly (UMD2), based on the same raw sequence data, was also published. The two assemblies have been subsequently updated to Btau_4.2 and UMD3.1, respectively.</p> <p>Results</p> <p>We compared the Btau_4.2 and UMD3.1 alternate assemblies. Inconsistencies were grouped into three main categories: (i) DNA segments showing almost coincidental chromosomal mapping but discordant orientation (inversions); (ii) DNA segments showing a discordant map position along the same chromosome; and (iii) sequences present in one chromosomal assembly but absent in the corresponding chromosome of the other assembly. The latter category mainly consisted of large amounts of scaffolds that were unassigned in Btau_4.2 but successfully mapped in UMD3.1. We sampled 70 inconsistencies and identified appropriate cow BACs for each of them. These clones were then utilized in FISH experiments on cow metaphase or interphase nuclei in order to disambiguate the discrepancies. In almost all instances the FISH results agreed with the UMD3.1 assembly. Occasionally, however, the mapping data of both assemblies were discordant with the FISH results.</p> <p>Conclusions</p> <p>Our work demonstrates how FISH, which is assembly independent, can be efficiently used to solve assembly problems frequently encountered using the shotgun approach.</p

Triple-gated motion and blood pool clearance corrections improve reproducibility of coronary 18F-NaF PET

PurposeTo improve the test-retest reproducibility of coronary plaque 18F-sodium fluoride (18F-NaF) positron emission tomography (PET) uptake measurements.MethodsWe recruited 20 patients with coronary artery disease who underwent repeated hybrid PET/CT angiography (CTA) imaging within 3&nbsp;weeks. All patients had 30-min PET acquisition and CTA during a single imaging session. Five PET image-sets with progressive motion correction were reconstructed: (i) a static dataset (no-MC), (ii) end-diastolic PET (standard), (iii) cardiac motion corrected (MC), (iv) combined cardiac and gross patient motion corrected (2 × MC) and, (v) cardiorespiratory and gross patient motion corrected (3 × MC). In addition to motion correction, all datasets were corrected for variations in the background activities which are introduced by variations in the injection-to-scan delays (background blood pool clearance correction, BC). Test-retest reproducibility of PET target-to-background ratio (TBR) was assessed by Bland-Altman analysis and coefficient of reproducibility.ResultsA total of 47 unique coronary lesions were identified on CTA. Motion correction in combination with BC improved the PET TBR test-retest reproducibility for all lesions (coefficient of reproducibility: standard = 0.437, no-MC = 0.345 (27% improvement), standard&nbsp;+&nbsp;BC = 0.365 (20% improvement), no-MC + BC = 0.341 (27% improvement), MC + BC = 0.288 (52% improvement), 2 × MC + BC = 0.278 (57% improvement) and 3 × C + BC = 0.254 (72% improvement), all p &lt; 0.001). Importantly, in a sub-analysis of 18F-NaF-avid lesions with gross patient motion &gt;&nbsp;10&nbsp;mm following corrections, reproducibility was improved by 133% (coefficient of reproducibility: standard = 0.745, 3 × MC = 0.320).ConclusionJoint corrections for cardiac, respiratory, and gross patient motion in combination with background blood pool corrections markedly improve test-retest reproducibility of coronary 18F-NaF PET

The Evolutionary Origin of Man Can Be Traced in the Layers of Defunct Ancestral Alpha Satellites Flanking the Active Centromeres of Human Chromosomes

Alpha satellite domains that currently function as centromeres of human chromosomes are flanked by layers of older alpha satellite, thought to contain dead centromeres of primate progenitors, which lost their function and the ability to homogenize satellite repeats, upon appearance of a new centromere. Using cladistic analysis of alpha satellite monomers, we elucidated complete layer patterns on chromosomes 8, 17, and X and related them to each other and to primate alpha satellites. We show that discrete and chronologically ordered alpha satellite layers are partially symmetrical around an active centromere and their succession is partially shared in non-homologous chromosomes. The layer structure forms a visual representation of the human evolutionary lineage with layers corresponding to ancestors of living primates and to entirely fossil taxa. Surprisingly, phylogenetic comparisons suggest that alpha satellite arrays went through periods of unusual hypermutability after they became “dead” centromeres. The layer structure supports a model of centromere evolution where new variants of a satellite repeat expanded periodically in the genome by rounds of inter-chromosomal transfer/amplification. Each wave of expansion covered all or many chromosomes and corresponded to a new primate taxon. Complete elucidation of the alpha satellite phylogenetic record would give a unique opportunity to number and locate the positions of major extinct taxa in relation to human ancestors shared with extant primates. If applicable to other satellites in non-primate taxa, analysis of centromeric layers could become an invaluable tool for phylogenetic studies

A comparison of genomic copy number calls by Partek Genomics Suite, Genotyping Console and Birdsuite algorithms to quantitative PCR

<p>Abstract</p> <p>Background</p> <p>Copy number variants are >1 kb genomic amplifications or deletions that can be identified using array platforms. However, arrays produce substantial background noise that contributes to high false discovery rates of variants. We hypothesized that quantitative PCR could finitely determine copy number and assess the validity of calling algorithms.</p> <p>Results</p> <p>Using data from 29 Affymetrix SNP 6.0 arrays, we determined copy numbers using three programs: Partek Genomics Suite, Affymetrix Genotyping Console 2.0 and Birdsuite. We compared array calls at 25 chromosomal regions to those determined by qPCR and found nearly identical calls in regions of copy number 2. Conversely, agreement differed in regions called variant by at least one method. The highest overall agreement in calls, 91%, was between Birdsuite and quantitative PCR. Partek Genomics Suite calls agreed with quantitative PCR 76% of the time while the agreement of Affymetrix Genotyping Console 2.0 with quantitative PCR was 79%.</p> <p>Conclusions</p> <p>In 38 independent samples, 96% of Birdsuite calls agreed with quantitative PCR. Analysis of three copy number calling programs and quantitative PCR showed Birdsuite to have the greatest agreement with quantitative PCR.</p

Large Tandem, Higher Order Repeats and Regularly Dispersed Repeat Units Contribute Substantially to Divergence Between Human and Chimpanzee Y Chromosomes

Comparison of human and chimpanzee genomes has received much attention, because of paramount role for understanding evolutionary step distinguishing us from our closest living relative. In order to contribute to insight into Y chromosome evolutionary history, we study and compare tandems, higher order repeats (HORs), and regularly dispersed repeats in human and chimpanzee Y chromosome contigs, using robust Global Repeat Map algorithm. We find a new type of long-range acceleration, human-accelerated HOR regions. In peripheral domains of 35mer human alphoid HORs, we find riddled features with ten additional repeat monomers. In chimpanzee, we identify 30mer alphoid HOR. We construct alphoid HOR schemes showing significant human-chimpanzee difference, revealing rapid evolution after human-chimpanzee separation. We identify and analyze over 20 large repeat units, most of them reported here for the first time as: chimpanzee and human ~1.6 kb 3mer secondary repeat unit (SRU) and ~23.5 kb tertiary repeat unit (~0.55 kb primary repeat unit, PRU); human 10848, 15775, 20309, 60910, and 72140 bp PRUs; human 3mer SRU (~2.4 kb PRU); 715mer and 1123mer SRUs (5mer PRU); chimpanzee 5096, 10762, 10853, 60523 bp PRUs; and chimpanzee 64624 bp SRU (10853 bp PRU). We show that substantial human-chimpanzee differences are concentrated in large repeat structures, at the level of as much as ~70% divergence, sizably exceeding previous numerical estimates for some selected noncoding sequences. Smeared over the whole sequenced assembly (25 Mb) this gives ~14% human--chimpanzee divergence. This is significantly higher estimate of divergence between human and chimpanzee than previous estimates.Comment: 22 pages, 7 figures, 12 tables. Published in Journal of Molecular Evolutio

Genome-Wide Analysis of Copy Number Variation in Type 1 Diabetes

Type 1 diabetes (T1D) tends to cluster in families, suggesting there may be a genetic component predisposing to disease. However, a recent large-scale genome-wide association study concluded that identified genetic factors, single nucleotide polymorphisms, do not account for overall familiality. Another class of genetic variation is the amplification or deletion of >1 kilobase segments of the genome, also termed copy number variations (CNVs). We performed genome-wide CNV analysis on a cohort of 20 unrelated adults with T1D and a control (Ctrl) cohort of 20 subjects using the Affymetrix SNP Array 6.0 in combination with the Birdsuite copy number calling software. We identified 39 CNVs as enriched or depleted in T1D versus Ctrl. Additionally, we performed CNV analysis in a group of 10 monozygotic twin pairs discordant for T1D. Eleven of these 39 CNVs were also respectively enriched or depleted in the Twin cohort, suggesting that these variants may be involved in the development of islet autoimmunity, as the presently unaffected twin is at high risk for developing islet autoimmunity and T1D in his or her lifetime. These CNVs include a deletion on chromosome 6p21, near an HLA-DQ allele. CNVs were found that were both enriched or depleted in patients with or at high risk for developing T1D. These regions may represent genetic variants contributing to development of islet autoimmunity in T1D

Propagation of RML Prions in Mice Expressing PrP Devoid of GPI Anchor Leads to Formation of a Novel, Stable Prion Strain

PrPC, a host protein which in prion-infected animals is converted to PrPSc, is linked to the cell membrane by a GPI anchor. Mice expressing PrPC without GPI anchor (tgGPI- mice), are susceptible to prion infection but accumulate anchorless PrPSc extra-, rather than intracellularly. We investigated whether tgGPI− mice could faithfully propagate prion strains despite the deviant structure and location of anchorless PrPSc. We found that RML and ME7, but not 22L prions propagated in tgGPI− brain developed novel cell tropisms, as determined by the Cell Panel Assay (CPA). Surprisingly, the levels of proteinase K-resistant PrPSc (PrPres) in RML- or ME7-infected tgGPI− brain were 25–50 times higher than in wild-type brain. When returned to wild-type brain, ME7 prions recovered their original properties, however RML prions had given rise to a novel prion strain, designated SFL, which remained unchanged even after three passages in wild-type mice. Because both RML PrPSc and SFL PrPSc are stably propagated in wild-type mice we propose that the two conformations are separated by a high activation energy barrier which is abrogated in tgGPI− mice