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.     

An analysis approach to identify specific functional sites in orthologous proteins using sequence and structural information: Application to neuroserpin reveals regions that differentially regulate inhibitory activity

The analysis of sequence conservation is commonly used to predict functionally important sites in proteins. We have developed an approach that first identifies highly conserved sites in a set of orthologous sequences using a weighted substitution‐matrix‐based conservation score and then filters these conserved sites based on the pattern of conservation present in a wider alignment of sequences from the same family and structural information to identify surface‐exposed sites. This allows us to detect specific functional sites in the target protein and exclude regions that are likely to be generally important for the structure or function of the wider protein family. We applied our method to two members of the serpin family of serine protease inhibitors. We first confirmed that our method successfully detected the known heparin binding site in antithrombin while excluding residues known to be generally important in the serpin family. We next applied our sequence analysis approach to neuroserpin and used our results to guide site‐directed polyalanine mutagenesis experiments. The majority of the mutant neuroserpin proteins were found to fold correctly and could still form inhibitory complexes with tissue plasminogen activator (tPA). Kinetic analysis of tPA inhibition, however, revealed altered inhibitory kinetics in several of the mutant proteins, with some mutants showing decreased association with tPA and others showing more rapid dissociation of the covalent complex. Altogether, these results confirm that our sequence analysis approach is a useful tool that can be used to guide mutagenesis experiments for the detection of specific functional sites in proteins. Proteins 2015; 83:135–152. © 2014 Wiley Periodicals, Inc.     

酵母双杂交系统筛选与RapGAP相互作用的蛋白质

目的筛选与Rap GAP相互作用的蛋白质,为进一步研究人源Rap1GAP介导的信号转导通路、揭示其与肿瘤的关系提供实验依据。方法选用与Rap1GAP同源的来自美丽线虫的Rap GAP作为饵蛋白,以来源于美丽线虫的c DNA文库作为靶蛋白,应用p PC97、p PC86组成的酵母双杂交系统筛选c DNA文库中与Rap GAP相互作用的蛋白质。结果通过营养缺陷平板(-LTH)筛选出63个拟似阳性菌落。经过Lac Z鉴定,19个菌落为阳性,其中7个为强阳性。提取来自19个酵母菌落中的重组DNA,经PCR扩增,12个菌落出现阳性结果。将该19个重组DNA分别电转化入DH5α细菌,涂板培养后,每板挑取4~10个克隆,通过Sal I和Not I双酶切鉴定进行阳性克隆筛选。将阳性克隆的重组DNA进行序列测定。测序结果与Gen Bank比较,其中4个克隆的DNA片段为Y39b6a基因片段、2个为Rap GAP、1个为苯丙氨酸-4-羟化酶、1个为细胞色素C氧化酶,还有1个DNA片段编码美丽线虫特有的小分子蛋白的基因片段,其余11个DNA片段不编码已知蛋白质。结论初步筛选出与Rap GAP相互作用的蛋白质,特别是其中有2个克隆为Rap GAP,提示Rap GAP可能以二聚体的方式存在。     

Genomic and experimental evidence for multiple metabolic functions in the RidA/YjgF/YER057c/UK114 (Rid) protein family

Background

It is now recognized that enzymatic or chemical side-reactions can convert normal metabolites to useless or toxic ones and that a suite of enzymes exists to mitigate such metabolite damage. Examples are the reactive imine/enamine intermediates produced by threonine dehydratase, which damage the pyridoxal 5''-phosphate cofactor of various enzymes causing inactivation. This damage is pre-empted by RidA proteins, which hydrolyze the imines before they do harm. RidA proteins belong to the YjgF/YER057c/UK114 family (here renamed the Rid family). Most other members of this diverse and ubiquitous family lack defined functions.

Results

Phylogenetic analysis divided the Rid family into a widely distributed, apparently archetypal RidA subfamily and seven other subfamilies (Rid1 to Rid7) that are largely confined to bacteria and often co-occur in the same organism with RidA and each other. The Rid1 to Rid3 subfamilies, but not the Rid4 to Rid7 subfamilies, have a conserved arginine residue that, in RidA proteins, is essential for imine-hydrolyzing activity. Analysis of the chromosomal context of bacterial RidA genes revealed clustering with genes for threonine dehydratase and other pyridoxal 5''-phosphate-dependent enzymes, which fits with the known RidA imine hydrolase activity. Clustering was also evident between Rid family genes and genes specifying FAD-dependent amine oxidases or enzymes of carbamoyl phosphate metabolism. Biochemical assays showed that Salmonella enterica RidA and Rid2, but not Rid7, can hydrolyze imines generated by amino acid oxidase. Genetic tests indicated that carbamoyl phosphate overproduction is toxic to S. enterica cells lacking RidA, and metabolomic profiling of Rid knockout strains showed ten-fold accumulation of the carbamoyl phosphate-related metabolite dihydroorotate.

Conclusions

Like the archetypal RidA subfamily, the Rid2, and probably the Rid1 and Rid3 subfamilies, have imine-hydrolyzing activity and can pre-empt damage from imines formed by amine oxidases as well as by pyridoxal 5''-phosphate enzymes. The RidA subfamily has an additional damage pre-emption role in carbamoyl phosphate metabolism that has yet to be biochemically defined. Finally, the Rid4 to Rid7 subfamilies appear not to hydrolyze imines and thus remain mysterious.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1584-3) contains supplementary material, which is available to authorized users.     

Surface exploration of Amphibalanus amphitrite cyprids on microtextured surfaces

Microtopography is one of several strategies used by marine organisms to inhibit colonization by fouling organisms. While replicates of natural microtextures discourage settlement, details of larval interactions with the structured surfaces remain scarce. Close-range microscopy was used to quantify the exploration of cyprids of Amphibalanus amphitrite on cylindrical micropillars with heights of 5 and 30 μm and diameters ranging from 5 to 100 μm. While 5 μm-high structures had little impact, 30 μm-high pillars significantly influenced cyprid exploration. An observed step length decrease and step duration increase on 5 μm diameter pillars is attributed to the small dimensions of the voids excluding the cyprid's attachment disc and consequently reducing the area of adhesive contact. When exploring larger diameter pillars, cyprids preferred using the voids to form temporary attachment points. This may enhance their resistance to flow. No-choice assay settlement patterns mirrored this exploration behaviour, albeit in a pattern counter to what was predicted.     

Optimizing Human Synovial Fluid Preparation for Two-Dimensional Gel Electrophoresis

Background

Prenatal screening for Down Syndrome (DS) would benefit from an increased number of biomarkers to improve sensitivity and specificity. Improving sensitivity and specificity would decrease the need for potentially risky invasive diagnostic procedures.

Results

We have performed an in depth two-dimensional difference gel electrophoresis (2D DIGE) study to identify potential biomarkers. We have used maternal plasma samples obtained from first and second trimesters from mothers carrying DS affected fetuses compared with mothers carrying normal fetuses. Plasma samples were albumin/IgG depleted and expanded pH ranges of pH 4.5 - 5.5, pH 5.3 - 6.5 and pH 6 - 9 were used for two-dimensional gel electrophoresis (2DE). We found no differentially expressed proteins in the first trimester between the two groups. Significant up-regulation of ceruloplasmin, inter-alpha-trypsin inhibitor heavy chain H4, complement proteins C1s subcomponent, C4-A, C5, and C9 and kininogen 1 were detected in the second trimester in maternal plasma samples where a DS affected fetus was being carried. However, ceruloplasmin could not be confirmed as being consistently up-regulated in DS affected pregnancies by Western blotting.

Conclusions

Despite the in depth 2DE approach used in this study the results underline the deficiencies of gel-based proteomics for detection of plasma biomarkers. Gel-free approaches may be more productive to increase the number of plasma biomarkers for DS for non-invasive prenatal screening and diagnosis.     

Transgenic expression of trehalulose synthase results in high concentrations of the sucrose isomer trehalulose in mature stems of field‐grown sugarcane

Sugarcane plants were developed that produce the sucrose isomers trehalulose and isomaltulose through expression of a vacuole‐targeted trehalulose synthase modified from the gene in ‘Pseudomonas mesoacidophila MX‐45’ and controlled by the maize ubiquitin (Ubi‐1) promoter. Trehalulose concentration in juice increased with internode maturity, reaching about 600 mm , with near‐complete conversion of sucrose in the most mature internodes. Plants remained vigorous, and trehalulose production in selected lines was retained over multiple vegetative generations under glasshouse and field conditions.     

Underwater application of quantitative PCR on an ocean mooring

The Environmental Sample Processor (ESP) is a device that allows for the underwater, autonomous application of DNA and protein probe array technologies as a means to remotely identify and quantify, in situ, marine microorganisms and substances they produce. Here, we added functionality to the ESP through the development and incorporation of a module capable of solid-phase nucleic acid extraction and quantitative PCR (qPCR). Samples collected by the instrument were homogenized in a chaotropic buffer compatible with direct detection of ribosomal RNA (rRNA) and nucleic acid purification. From a single sample, both an rRNA community profile and select gene abundances were ascertained. To illustrate this functionality, we focused on bacterioplankton commonly found along the central coast of California and that are known to vary in accordance with different oceanic conditions. DNA probe arrays targeting rRNA revealed the presence of 16S rRNA indicative of marine crenarchaea, SAR11 and marine cyanobacteria; in parallel, qPCR was used to detect 16S rRNA genes from the former two groups and the large subunit RuBisCo gene (rbcL) from Synecchococcus. The PCR-enabled ESP was deployed on a coastal mooring in Monterey Bay for 28 days during the spring-summer upwelling season. The distributions of the targeted bacterioplankon groups were as expected, with the exception of an increase in abundance of marine crenarchaea in anomalous nitrate-rich, low-salinity waters. The unexpected co-occurrence demonstrated the utility of the ESP in detecting novel events relative to previously described distributions of particular bacterioplankton groups. The ESP can easily be configured to detect and enumerate genes and gene products from a wide range of organisms. This study demonstrated for the first time that gene abundances could be assessed autonomously, underwater in near real-time and referenced against prevailing chemical, physical and bulk biological conditions.     

Patterns of diversifying selection in the phytotoxin-like scr74 gene family of Phytophthora infestans

Phytophthora infestans, the organism responsible for the Irish famine, causes late blight, a re-emerging disease of potato and tomato. Little is known about the molecular evolution of P. infestans genes. To identify candidate effector genes (virulence or avirulence genes) that may have co-evolved with the host, we mined expressed sequence tag (EST) data from infection stages of P. infestans for secreted and potentially polymorphic genes. This led to the identification of scr74, a gene that encodes a predicted 74-amino acid secreted cysteine-rich protein with similarity to the Phytophthora cactorum phytotoxin PcF. The expression of scr74 was upregulated approximately 60-fold 2 to 4 days after inoculation of tomato and was also significantly induced during early stages of colonization of potato. The scr74 gene was found to belong to a highly polymorphic gene family within P. infestans with 21 different sequences identified. Using the approximate and maximum likelihood (ML) methods, we found that diversifying selection likely caused the extensive polymorphism observed within the scr74 gene family. Pairwise comparisons of 17 scr74 sequences revealed elevated ratios of nonsynonymous to synonymous nucleotide-substitution rates, particularly in the mature region of the proteins. Using ML, all 21 polymorphic amino acid sites were identified to be under diversifying selection. Of these 21 amino acids, 19 are located in the mature protein region, suggesting that selection may have acted on the functional portions of the proteins. Further investigation of gene copy number and organization revealed that the scr74 gene family comprises at least three copies located in a region of no more than 300 kb of the P. infestans genome. We found evidence that recombination contributed to sequence divergence within at least one gene locus. These results led us to propose an evolutionary model that involves gene duplication and recombination, followed by functional divergence of scr74 genes. This study provides support for using diversifying selection as a criterion for identifying candidate effector genes from sequence databases.