Assessment of algorithms for high throughput detection of genomic copy number variation in oligonucleotide microarray data

Background

Genomic deletions and duplications are important in the pathogenesis of diseases, such as cancer and mental retardation, and have recently been shown to occur frequently in unaffected individuals as polymorphisms. Affymetrix GeneChip whole genome sampling analysis (WGSA) combined with 100 K single nucleotide polymorphism (SNP) genotyping arrays is one of several microarray-based approaches that are now being used to detect such structural genomic changes. The popularity of this technology and its associated open source data format have resulted in the development of an increasing number of software packages for the analysis of copy number changes using these SNP arrays.

Results

We evaluated four publicly available software packages for high throughput copy number analysis using synthetic and empirical 100 K SNP array data sets, the latter obtained from 107 mental retardation (MR) patients and their unaffected parents and siblings. We evaluated the software with regards to overall suitability for high-throughput 100 K SNP array data analysis, as well as effectiveness of normalization, scaling with various reference sets and feature extraction, as well as true and false positive rates of genomic copy number variant (CNV) detection.

Conclusion

We observed considerable variation among the numbers and types of candidate CNVs detected by different analysis approaches, and found that multiple programs were needed to find all real aberrations in our test set. The frequency of false positive deletions was substantial, but could be greatly reduced by using the SNP genotype information to confirm loss of heterozygosity.     

Caspase-mediated cleavage of the exosome subunit PM/Scl-75 during apoptosis

Recent studies have implicated the dying cell as a potential reservoir of modified autoantigens that might initiate and drive systemic autoimmunity in susceptible hosts. A number of subunits of the exosome, a complex of 3'→5' exoribonucleases that functions in a variety of cellular processes, are recognized by the so-called anti-PM/Scl autoantibodies, found predominantly in patients suffering from an overlap syndrome of myositis and scleroderma. Here we show that one of these subunits, PM/Scl-75, is cleaved during apoptosis. PM/Scl-75 cleavage is inhibited by several different caspase inhibitors. The analysis of PM/Scl-75 cleavage by recombinant caspase proteins shows that PM/Scl-75 is efficiently cleaved by caspase-1, to a smaller extent by caspase-8, and relatively inefficiently by caspase-3 and caspase-7. Cleavage of the PM/Scl-75 protein occurs in the C-terminal part of the protein at Asp369 (IILD³⁶⁹↓G), and at least a fraction of the resulting N-terminal fragments of PM/Scl-75 remains associated with the exosome. Finally, the implications of PM/Scl-75 cleavage for exosome function and the generation of anti-PM/Scl-75 autoantibodies are discussed.     

Increased incidence of pregnancy complications in women who later develop scleroderma: a case control study

Introduction  

Studies have shown that fetal progenitor cells persist in maternal blood or bone marrow for more than 30 years after delivery. Increased trafficking of fetal cells occurs during pregnancy complications, such as hypertension, preeclampsia, miscarriage and intra-uterine growth restriction (IUGR). Women with these pregnancy complications are significantly more often HLA-class II compatible with their spouses. Women who later develop scleroderma also give birth to an HLA-class II child more often. From these prior studies we hypothesized that preeclampsia and other pregnancy complications could be associated with increased levels of fetal cell trafficking, and later be involved in the development of scleroderma.     

Conceptual approach to renewable barrier film design based on wood hydrolysate

Biomass is converted to oxygen barriers through a conceptually unconventional approach involving the preservation of the biomass native interactions and macromolecular components and enhancing the effect by created interactions with a co-component. A combined calculation/assessment model is elaborated to understand, quantify, and predict which compositions that provide an intermolecular affinity high enough to mediate the molecular packing needed to create a functioning barrier. The biomass used is a wood hydrolysate, a polysaccharide-rich but not highly refined mixture where a fair amount of the native intermolecular and intramolecular hemicelluloses-lignin interactions are purposely preserved, resulting in barriers with very low oxygen permeabilities (OP) both at 50 and 80% relative humidity and considerably lower OPs than coatings based on the corresponding highly purified spruce hemicellulose, O-acetyl galactoglucomannan (AcGGM). The component interactions and mutual affinities effectively mediate an immobilization of the chain segments in a dense disordered structure, modeled through the Hansen's solubility parameter concept and quantified on the nanolength scale by positron annihilation lifetime spectrum (PALS).     

Membrane localization diversity of TPK channels and their physiological role

Potassium (K) is one of the major nutrients that is essential for plant growth and development. The majority of cellular K⁺ resides in the vacuole and tonoplast K⁺ channels of the TPK (Two Pore K) family are main players in cellular K⁺ homeostasis. All TPK channels were previously reported to be expressed in the tonoplast of the large central lytic vacuole (LV) except for one isoform in Arabidopsis that resides in the plasma membrane. However, plant cells often contain more than one type of vacuole that coexist in the same cell. We recently showed that two TPK isoforms (OsTPKa and OsTPKb) from Oryza sativa localize to different vacuoles with OsTPKa predominantly found in the LV tonoplast and OsTPKb primarily in smaller compartments that resemble small vacuoles (SVs). Our study further revealed that it is the C-terminal domain that determines differential targeting of OsTPKa and OsTPKb. Three C-terminal amino acids were particularly relevant for targeting TPKs to their respective endomembranes. In this addendum we further evaluate how the different localization of TPKa and TPKb impact on their physiological role and how TPKs provide a potential tool to study the physiology of different types of vacuole.Key words: TPK channels, small vacuoles, vacuolar targeting, potassiumThe roles of plant vacuolar K⁺ channels are diverse and include potassium homeostasis, turgor regulation and responses to abiotic stress. Vacuolar K⁺-selective channels belong to two-pore K⁺ (TPK) channel families which have been found in genomes of many plant species such as Arabidopsis, poplar, Physcomitrella, Eucalyptus, barley, potato, rice and tobacco (Fig. 1). TPKs have structural similarity to mammalian “tandem P domain” channels with a secondary structure that contains four transmembrane domains and two pore regions (Fig. 2).^1–5 TPK channels have pore regions with a GYGD signature that endows K⁺ selectivity and a variable number of Ca²⁺ binding EF domains in the C terminus.^3–8 One of the best characterized members of the TPK family is AtTPK1 from Arabidopsis thaliana. AtTPK1 activity is voltage independent but sensitive to cytosolic Ca²⁺, cytosolic pH and N-terminal phosphorylation by 14-3-3 proteins.^5,6,8,9 In Arabidopsis, AtTPK1 expresses in the large lytic vacuole (LV) and plays roles in cellular K⁺ homeostasis, K⁺-release during stomatal closure and seed germination.^4,5 Other members of the Arabidopsis TPK family (AtTPK2, AtTPK3, AtTPK5) have been shown to localize to the LV but also showed some expression in smaller, vesicle-like, compartments.⁴ However, none of these isoforms appears to form functional channels in planta although our experiments with heterologous expression of AtTPK3 and AtTPK5 in the K⁺ uptake deficient E. coli LB2003 demonstrates complementation of bacterial growth phenotype (Isayenkov S, et al. unpublished results). Equally intriguing, is the plasma membrane localization of the Arabidopsis TPK4 isoform, in spite of its sequence being very similar to that of other TPKs.¹⁰Open in a separate windowFigure 1Phylogenetic tree of plant TPKs. The three main clusters of TPKs comprise: Cluster 1 with AtTPK1-like channels; Cluster 2 with AtTPK3/TPK5-like channels; Cluster 3 with barley HvTPKb. Bootstrap analysis was performed using ‘Molecular Evolutionary Genetics Analysis, MEGA4’ software available at www.megasoftware.net/mega4/megaOpen in a separate windowFigure 2Two-pore potassium channel secondary structure. TPK channels comprise four transmembrane domains (1–4) and two pore regions (P) per subunit. Functional channels are formed from two subunits. In most TPKs, both P regions contain a K⁺ selectivity signature, GYGD. However, the tobacco NtTPKa isoform has different motifs in the second P domain. In the N terminal region, TPKs have a 14-3-3 binding domain that impact on channel activity, with the binding of 14-3-3 protein leading to channel activation. C-termini of TPKs show a varying number of putative Ca²⁺ binding “EF hands” which may vary from zero to two.     

Genetic dissection of Al tolerance QTLs in the maize genome by high density SNP scan

Background

Aluminum (Al) toxicity is an important limitation to food security in tropical and subtropical regions. High Al saturation on acid soils limits root development, reducing water and nutrient uptake. In addition to naturally occurring acid soils, agricultural practices may decrease soil pH, leading to yield losses due to Al toxicity. Elucidating the genetic and molecular mechanisms underlying maize Al tolerance is expected to accelerate the development of Al-tolerant cultivars.

Results

Five genomic regions were significantly associated with Al tolerance, using 54,455 SNP markers in a recombinant inbred line population derived from Cateto Al237. Candidate genes co-localized with Al tolerance QTLs were further investigated. Near-isogenic lines (NILs) developed for ZmMATE2 were as Al-sensitive as the recurrent line, indicating that this candidate gene was not responsible for the Al tolerance QTL on chromosome 5, qALT5. However, ZmNrat1, a maize homolog to OsNrat1, which encodes an Al³⁺ specific transporter previously implicated in rice Al tolerance, was mapped at ~40 Mbp from qALT5. We demonstrate for the first time that ZmNrat1 is preferentially expressed in maize root tips and is up-regulated by Al, similarly to OsNrat1 in rice, suggesting a role of this gene in maize Al tolerance. The strongest-effect QTL was mapped on chromosome 6 (qALT6), within a 0.5 Mbp region where three copies of the Al tolerance gene, ZmMATE1, were found in tandem configuration. qALT6 was shown to increase Al tolerance in maize; the qALT6-NILs carrying three copies of ZmMATE1 exhibited a two-fold increase in Al tolerance, and higher expression of ZmMATE1 compared to the Al sensitive recurrent parent. Interestingly, a new source of Al tolerance via ZmMATE1 was identified in a Brazilian elite line that showed high expression of ZmMATE1 but carries a single copy of ZmMATE1.

Conclusions

High ZmMATE1 expression, controlled either by three copies of the target gene or by an unknown molecular mechanism, is responsible for Al tolerance mediated by qALT6. As Al tolerant alleles at qALT6 are rare in maize, marker-assisted introgression of this QTL is an important strategy to improve maize adaptation to acid soils worldwide.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-153) contains supplementary material, which is available to authorized users.