Automated scoring of AFLPs using RawGeno v 2.0, a free R CRAN library.

Amplified Fragment Length Polymorphisms (AFLPs) are a cheap and efficient protocol for generating large sets of genetic markers. This technique has become increasingly used during the last decade in various fields of biology, including population genomics, phylogeography, and genome mapping. Here, we present RawGeno, an R library dedicated to the automated scoring of AFLPs (i.e., the coding of electropherogram signals into ready-to-use datasets). Our program includes a complete suite of tools for binning, editing, visualizing, and exporting results obtained from AFLP experiments. RawGeno can either be used with command lines and program analysis routines or through a user-friendly graphical user interface. We describe the whole RawGeno pipeline along with recommendations for (a) setting the analysis of electropherograms in combination with PeakScanner, a program freely distributed by Applied Biosystems; (b) performing quality checks; (c) defining bins and proceeding to scoring; (d) filtering nonoptimal bins; and (e) exporting results in different formats.

Blind docking of 260 protein-ligand complexes with EADock 2.0.

Molecular docking softwares are one of the important tools of modern drug development pipelines. The promising achievements of the last 10 years emphasize the need for further improvement, as reflected by several recent publications (Leach et al., J Med Chem 2006, 49, 5851; Warren et al., J Med Chem 2006, 49, 5912). Our initial approach, EADock, showed a good performance in reproducing the experimental binding modes for a set of 37 different ligand-protein complexes (Grosdidier et al., Proteins 2007, 67, 1010). This article presents recent improvements regarding the scoring and sampling aspects over the initial implementation, as well as a new seeding procedure based on the detection of cavities, opening the door to blind docking with EADock. These enhancements were validated on 260 complexes taken from the high quality Ligand Protein Database [LPDB, (Roche et al., J Med Chem 2001, 44, 3592)]. Two issues were identified: first, the quality of the initial structures cannot be assumed and a manual inspection and/or a search in the literature are likely to be required to achieve the best performance. Second the description of interactions involving metal ions still has to be improved. Nonetheless, a remarkable success rate of 65% was achieved for a large scale blind docking assay, when considering only the top ranked binding mode and a success threshold of 2 A RMSD to the crystal structure. When looking at the five-top ranked binding modes, the success rate increases up to 76%. In a standard local docking assay, success rates of 75 and 83% were obtained, considering only the top ranked binding mode, or the five top binding modes, respectively.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Affymetrix Whole-Transcript Human Gene 1.0 ST array is highly concordant with standard 3' expression arrays.

The recently released Affymetrix Human Gene 1.0 ST array has two major differences compared with standard 3' based arrays: (i) it interrogates the entire mRNA transcript, and (ii) it uses DNA targets. To assess the impact of these differences on array performance, we performed a series of comparative hybridizations between the Human Gene 1.0 ST and the Affymetrix HG-U133 Plus 2.0 and the Illumina HumanRef-8 BeadChip arrays. Additionally, both RNA and DNA targets were hybridized on HG-U133 Plus 2.0 arrays. The results show that the overall reproducibility of the Gene 1.0 ST array is best. When looking only at the high intensity probes, the reproducibility of the Gene 1.0 ST array and the Illumina BeadChip array is equally good. Concordance of array results was assessed using different inter-platform mappings. Agreements are best between the two labeling protocols using HG-U133 Plus 2.0 array. The Gene 1.0 ST array is most concordant with the HG-U133 array hybridized with cDNA targets. This may reflect the impact of the target type. Overall, the high degree of correspondence provides strong evidence for the reliability of the Gene 1.0 ST array.

Robust volume-targeted balanced steady-state free-precession coronary magnetic resonance angiography in a breathhold at 3.0 Tesla: a reproducibility study.

BACKGROUND: Transient balanced steady-state free-precession (bSSFP) has shown substantial promise for noninvasive assessment of coronary arteries but its utilization at 3.0 T and above has been hampered by susceptibility to field inhomogeneities that degrade image quality. The purpose of this work was to refine, implement, and test a robust, practical single-breathhold bSSFP coronary MRA sequence at 3.0 T and to test the reproducibility of the technique. METHODS: A 3D, volume-targeted, high-resolution bSSFP sequence was implemented. Localized image-based shimming was performed to minimize inhomogeneities of both the static magnetic field and the radio frequency excitation field. Fifteen healthy volunteers and three patients with coronary artery disease underwent examination with the bSSFP sequence (scan time = 20.5 ± 2.0 seconds), and acquisitions were repeated in nine subjects. The images were quantitatively analyzed using a semi-automated software tool, and the repeatability and reproducibility of measurements were determined using regression analysis and intra-class correlation coefficient (ICC), in a blinded manner. RESULTS: The 3D bSSFP sequence provided uniform, high-quality depiction of coronary arteries (n = 20). The average visible vessel length of 100.5 ± 6.3 mm and sharpness of 55 ± 2% compared favorably with earlier reported navigator-gated bSSFP and gradient echo sequences at 3.0 T. Length measurements demonstrated a highly statistically significant degree of inter-observer (r = 0.994, ICC = 0.993), intra-observer (r = 0.894, ICC = 0.896), and inter-scan concordance (r = 0.980, ICC = 0.974). Furthermore, ICC values demonstrated excellent intra-observer, inter-observer, and inter-scan agreement for vessel diameter measurements (ICC = 0.987, 0.976, and 0.961, respectively), and vessel sharpness values (ICC = 0.989, 0.938, and 0.904, respectively). CONCLUSIONS: The 3D bSSFP acquisition, using a state-of-the-art MR scanner equipped with recently available technologies such as multi-transmit, 32-channel cardiac coil, and localized B0 and B1+ shimming, allows accelerated and reproducible multi-segment assessment of the major coronary arteries at 3.0 T in a single breathhold. This rapid sequence may be especially useful for functional imaging of the coronaries where the acquisition time is limited by the stress duration and in cases where low navigator-gating efficiency prohibits acquisition of a free breathing scan in a reasonable time period.

A comparative phenotypic study of kallmann syndrome patients carrying monoallelic and biallelic mutations in the prokineticin 2 or prokineticin receptor 2 genes.

Context: Both biallelic and monoallelic mutations in PROK2 or PROKR2 have been found in Kallmann syndrome (KS). Objective: The objective of the study was to compare the phenotypes of KS patients harboring monoallelic and biallelic mutations in these genes. Design and Patients: We studied clinical and endocrine features that reflect the functioning of the pituitary-gonadal axis, and the nonreproductive phenotype, in 55 adult KS patients (42 men and 13 women), of whom 41 had monoallelic mutations and 14 biallelic mutations in PROK2 or PROKR2. Results: Biallelic mutations were associated with more frequent cryptorchidism (70% vs. 34%, P < 0.05) and microphallus (90% vs. 28%, P < 0.001) and lower mean testicular volume (1.2 +/- 0.4 vs. 4.5 +/- 6.0 ml; P < 0.01) in male patients. Likewise, the testosterone level as well as the basal FSH level and peak LH level under GnRH-stimulation were lower in males with biallelic mutations (0.2 +/- 0.1 vs. 0.7 +/- 0.8 ng/ml; P = 0.05, 0.3 +/- 0.1 vs. 1.8 +/- 3.0 IU/liter; P < 0.05, and 0.8 +/- 0.8 vs. 5.2 +/- 5.5 IU/liter; P < 0.05, respectively). Nonreproductive, nonolfactory anomalies were rare in both sexes and were never found in patients with biallelic mutations. The mean body mass index of the patients (23.9 +/- 4.2 kg/m(2) in males and 26.3 +/- 6.6 kg/m(2) in females) did not differ significantly from that of gender-, age-, and treatment-matched KS individuals who did not carry a mutation in PROK2 or PROKR2. Finally, circadian cortisol levels evaluated in five patients, including one with biallelic PROKR2 mutations, were normal in all cases. Conclusion: Male patients carrying biallelic mutations in PROK2 or PROKR2 have a less variable and on average a more severe reproductive phenotype than patients carrying monoallelic mutations in these genes. Nonreproductive, nonolfactory clinical anomalies associated with KS seem to be restricted to patients with monoallelic mutations.

Defining the Effect of the 16p11.2 Duplication on Cognition, Behavior, and Medical Comorbidities.

IMPORTANCE: The 16p11.2 BP4-BP5 duplication is the copy number variant most frequently associated with autism spectrum disorder (ASD), schizophrenia, and comorbidities such as decreased body mass index (BMI). OBJECTIVES: To characterize the effects of the 16p11.2 duplication on cognitive, behavioral, medical, and anthropometric traits and to understand the specificity of these effects by systematically comparing results in duplication carriers and reciprocal deletion carriers, who are also at risk for ASD. DESIGN, SETTING, AND PARTICIPANTS: This international cohort study of 1006 study participants compared 270 duplication carriers with their 102 intrafamilial control individuals, 390 reciprocal deletion carriers, and 244 deletion controls from European and North American cohorts. Data were collected from August 1, 2010, to May 31, 2015 and analyzed from January 1 to August 14, 2015. Linear mixed models were used to estimate the effect of the duplication and deletion on clinical traits by comparison with noncarrier relatives. MAIN OUTCOMES AND MEASURES: Findings on the Full-Scale IQ (FSIQ), Nonverbal IQ, and Verbal IQ; the presence of ASD or other DSM-IV diagnoses; BMI; head circumference; and medical data. RESULTS: Among the 1006 study participants, the duplication was associated with a mean FSIQ score that was lower by 26.3 points between proband carriers and noncarrier relatives and a lower mean FSIQ score (16.2-11.4 points) in nonproband carriers. The mean overall effect of the deletion was similar (-22.1 points; P < .001). However, broad variation in FSIQ was found, with a 19.4- and 2.0-fold increase in the proportion of FSIQ scores that were very low (≤40) and higher than the mean (>100) compared with the deletion group (P < .001). Parental FSIQ predicted part of this variation (approximately 36.0% in hereditary probands). Although the frequency of ASD was similar in deletion and duplication proband carriers (16.0% and 20.0%, respectively), the FSIQ was significantly lower (by 26.3 points) in the duplication probands with ASD. There also were lower head circumference and BMI measurements among duplication carriers, which is consistent with the findings of previous studies. CONCLUSIONS AND RELEVANCE: The mean effect of the duplication on cognition is similar to that of the reciprocal deletion, but the variance in the duplication is significantly higher, with severe and mild subgroups not observed with the deletion. These results suggest that additional genetic and familial factors contribute to this variability. Additional studies will be necessary to characterize the predictors of cognitive deficits.