BreakTrans: uncovering the genomic architecture of gene fusions

Producing gene fusions through genomic structural rearrangements is a major mechanism for tumor evolution. Therefore, accurately detecting gene fusions and the originating rearrangements is of great importance for personalized cancer diagnosis and targeted therapy. We present a tool, BreakTrans, that systematically maps predicted gene fusions to structural rearrangements. Thus, BreakTrans not only validates both types of predictions, but also provides mechanistic interpretations. BreakTrans effectively validates known fusions and discovers novel events in a breast cancer cell line. Applying BreakTrans to 43 breast cancer samples in The Cancer Genome Atlas identifies 90 genomically validated gene fusions. BreakTrans is available at http://bioinformatics.mdanderson.org/main/BreakTrans     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Interpreting trends in the binding of cyclic ureas to HIV-1 protease

The design of new HIV protease inhibitors requires an improved understanding of the physical basis of inhibitor/protein binding. Here, the binding affinities of seven aliphatic cyclic ureas to HIV-1 protease are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed U-shaped trend of binding free energy as a function of aliphatic chain length. Interestingly, the decrease in affinity for the longest chains is attributable primarily to the energy cost of partly desolvating charged aspartic and arginine groups at the mouths of the active site. Even aliphatic chains too short to contact these charged groups directly are subject to considerable desolvation penalties. We are not aware of other systems where binding affinity trends have been attributed to long-ranged electrostatic desolvation of ionized groups. A generalized Born/surface area solvation model yields a much smaller change in desolvation energy with chain length and, therefore, does not reproduce the experimental binding affinity trends. This result suggests that the generalized Born model should be used with caution for complex, partly desolvated systems like protein binding sites. We also find that changing the assumed protonation state of the active site aspartyl dyad significantly affects the computed binding affinity trends. The protonation state of the aspartyl dyad in the presence of cyclic ureas is discussed in light of the observation that the monoprotonated state reproduces the experimental results best.     

中国的炭疽杆菌DNA分型及其地理分布

炭疽广泛分布于中国各地,特别是西部地区,并经常造成人畜疾病,在一项合作研究中,用多位点VNTR分析（MLVA）对从1952-1998年自中国主要地理流行区域分离的病人,病畜和土壤等来源的炭疽杆菌进行了基因分型,MLVA分析结果揭示了21种新的基因型,其等位基因组合在以前世界范围分离物的研究中未曾发现,此外,分离物的分群显示,A3b组是地理上最广泛分布的基因组,说明该组可能是中国的“地方流行株”。而来自古丝绸之路重要贸易中心新疆的大量分离株其基因型特别分散。     

Viral discovery and sequence recovery using DNA microarrays

Because of the constant threat posed by emerging infectious diseases and the limitations of existing approaches used to identify new pathogens, there is a great demand for new technological methods for viral discovery. We describe herein a DNA microarray-based platform for novel virus identification and characterization. Central to this approach was a DNA microarray designed to detect a wide range of known viruses as well as novel members of existing viral families; this microarray contained the most highly conserved 70mer sequences from every fully sequenced reference viral genome in GenBank. During an outbreak of severe acute respiratory syndrome (SARS) in March 2003, hybridization to this microarray revealed the presence of a previously uncharacterized coronavirus in a viral isolate cultivated from a SARS patient. To further characterize this new virus, approximately 1 kb of the unknown virus genome was cloned by physically recovering viral sequences hybridized to individual array elements. Sequencing of these fragments confirmed that the virus was indeed a new member of the coronavirus family. This combination of array hybridization followed by direct viral sequence recovery should prove to be a general strategy for the rapid identification and characterization of novel viruses and emerging infectious disease.     

Extending assembly of short DNA sequences to handle error

Inexpensive de novo genome sequencing, particularly in organisms with small genomes, is now possible using several new sequencing technologies. Some of these technologies such as that from Illumina's Solexa Sequencing, produce high genomic coverage by generating a very large number of small reads ( approximately 30 bp). While prior work shows that partial assembly can be performed by k-mer extension in error-free reads, this algorithm is unsuccessful with the sequencing error rates found in practice. We present VCAKE (Verified Consensus Assembly by K-mer Extension), a modification of simple k-mer extension that overcomes error by using high depth coverage. Though it is a simple modification of a previous approach, we show significant improvements in assembly results on simulated and experimental datasets that include error. AVAILABILITY: http://152.2.15.114/~labweb/VCAKE     

MicroRNAs in systemic rheumatic diseases

MicroRNAs (miRNAs) are endogenous, non-coding, single-stranded RNAs about 21 nucleotides in length. miRNAs have been shown to regulate gene expression and thus influence a wide range of physiological and pathological processes. Moreover, they are detected in a variety of sources, including tissues, serum, and other body fluids, such as saliva. The role of miRNAs is evident in various malignant and nonmalignant diseases, and there is accumulating evidence also for an important role of miRNAs in systemic rheumatic diseases. Abnormal expression of miRNAs has been reported in autoimmune diseases, mainly in systemic lupus erythematosus and rheumatoid arthritis. miRNAs can be aberrantly expressed even in the different stages of disease progression, allowing miRNAs to be important biomarkers, to help understand the pathogenesis of the disease, and to monitor disease activity and effects of treatment. Different groups have demonstrated a link between miRNA expression and disease activity, as in the case of renal flares in lupus patients. Moreover, miRNAs are emerging as potential targets for new therapeutic strategies of autoimmune disorders. Taken together, recent data demonstrate that miRNAs can influence mechanisms involved in the pathogenesis, relapse, and specific organ involvement of autoimmune diseases. The ultimate goal is the identification of a miRNA target or targets that could be manipulated through specific therapies, aiming at activation or inhibition of specific miRNAs responsible for the development of disease.     

Frequent coexistence of anti-topoisomerase I and anti-U1RNP autoantibodies in African American patients associated with mild skin involvement: a retrospective clinical study

Introduction  

The presence of anti-topoisomerase I (topo I) antibodies is a classic scleroderma (SSc) marker presumably associated with a unique clinical subset. Here the clinical association of anti-topo I was reevaluated in unselected patients seen in a rheumatology clinic setting.