Variation revealed by SNP genotyping and morphology provides insight into the origin of the tomato

Tomato, Solanum lycopersicum, is divided into two widely distributed varieties: the cultivated S. lycopersicum var. lycopersicum, and the weedy S. lycopersicum var. cerasiforme. Solanum pimpinellifolium is the most closely related wild species of tomato. The roles of S. pimpinellifolium and S. l. cerasiforme during the domestication of tomato are still under debate. Some authors consider S. l. cerasiforme to be the ancestor, whereas others think that S. l. cerasiforme is an admixture of S. pimpinellifolium and the cultivated S. l. lycopersicum. It is also not clear whether the domestication occurred in the Andean region or in Mesoamerica. We characterized 272 accessions (63 S. pimpinellifolium, 106 S. l. cerasiforme, 95 S. l. lycopersicum and 8 derived from hybridization processes) were morphologically and genetically using the SolCap platform (7,414 SNPs). The two species were distinguished in a PCA analysis and displayed a rich geographic structure. Solanum lycopersicum var. cerasiforme and S. l. lycopersicum were also differentiated in the PCA and Structure analyses, which supports maintaining them as different varieties. Solanum pimpinellifolium and the Andean S. l. cerasiforme were more diverse than the non-Andean S. lycopersicum. Solanum lycopersicum var. cerasiforme was morphologically and molecularly intermediate between S. pimpinellifolium and tomato. Solanum lycopersicum var. cerasiforme, with the exception of several Ecuadorian and Mexican accessions, is composed of the products of admixture processes according to the Structure analysis. The non-admixtured S. l. cerasiforme might be similar to the ancestral cultivars from which the cultivated tomato originated, and presents remarkable morphological diversity, including fruits of up to 6 cm in diameter. The data obtained would fit a model in which a pre-domestication took place in the Andean region, with the domestication being completed in Mesoamerica. Subsequently, the Spaniards took plants from Mesoamerica to Spain and from there they were exported to the rest of the world.Blanca Postigo, JM.; Cañizares Sales, J.; Cordero Romay, L.; Pascual Bañuls, L.; Díez Niclós, MJTDJ.; Nuez Viñals, F. (2012). Variation revealed by SNP genotyping and morphology provides insight into the origin of the tomato. PLoS ONE. 7(10):1-17. doi:10.1371/journal.pone.0048198S11771

ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence

Background: The possibilities offered by next generation sequencing (NGS) platforms are revolutionizing biotechnological laboratories. Moreover, the combination of NGS sequencing and affordable high-throughput genotyping technologies is facilitating the rapid discovery and use of SNPs in non-model species. However, this abundance of sequences and polymorphisms creates new software needs. To fulfill these needs, we have developed a powerful, yet easy-to-use application. Results: The ngs_backbone software is a parallel pipeline capable of analyzing Sanger, 454, Illumina and SOLiD (Sequencing by Oligonucleotide Ligation and Detection) sequence reads. Its main supported analyses are: read cleaning, transcriptome assembly and annotation, read mapping and single nucleotide polymorphism (SNP) calling and selection. In order to build a truly useful tool, the software development was paired with a laboratory experiment. All public tomato Sanger EST reads plus 14.2 million Illumina reads were employed to test the tool and predict polymorphism in tomato. The cleaned reads were mapped to the SGN tomato transcriptome obtaining a coverage of 4.2 for Sanger and 8.5 for Illumina. 23,360 single nucleotide variations (SNVs) were predicted. A total of 76 SNVs were experimentally validated, and 85% were found to be real. Conclusions: ngs_backbone is a new software package capable of analyzing sequences produced by NGS technologies and predicting SNVs with great accuracy. In our tomato example, we created a highly polymorphic collection of SNVs that will be a useful resource for tomato researchers and breeders. The software developed along with its documentation is freely available under the AGPL license and can be downloaded from http://bioinf. comav.upv.es/ngs_backbone/ or http://github.com/JoseBlanca/franklin.Blanca Postigo, JM.; Pascual Bañuls, L.; Ziarsolo Areitioaurtena, P.; Nuez Viñals, F.; Cañizares Sales, J. (2011). Ngs_backbone: a pipeline for read cleaning, mapping and SNP calling using Next Generation Sequence. BMC Genomics. 12:1-8. doi:10.1186/1471-2164-12-285S1812Metzker ML: Sequencing technologies - the next generation. Nature Reviews Genetics. 2010, 11 (1): 31-46. 10.1038/nrg2626.454 sequencing. [ http://www.454.com/ ]Illumina Inc. [ http://www.illumina.com/ ]Flicek P, Birney E: Sense from sequence reads: methods for alignment and assembly (vol 6, pg S6, 2009). Nature Methods. 2010, 7 (6): 479-479.Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WEG, Wetter T, Suhai S: Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Research. 2004, 14 (6): 1147-1159. 10.1101/gr.1917404.Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology. 2009, 10 (3):Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data P: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.1000 Genomes. A deep Catalog of Human Genetic Variation. [ http://1000genomes.org/wiki/doku.php?id=1000_genomes:analysis:vcf4.0 ]The seqanswers internet forum. [ http://seqanswers.com/ ]Blankenberg D, Taylor J, Schenck I, He JB, Zhang Y, Ghent M, Veeraraghavan N, Albert I, Miller W, Makova KD, Ross CH, Nekrutenko A: A framework for collaborative analysis of ENCODE data: Making large-scale analyses biologist-friendly. Genome Research. 2007, 17 (6): 960-964. 10.1101/gr.5578007.CloVR Automated Sequence Analysis from Your Desktop. [ http://clovr.org/ ]Papanicolaou A, Stierli R, Ffrench-Constant RH, Heckel DG: Next generation transcriptomes for next generation genomes using est2assembly. Bmc Bioinformatics. 2009, 10:Applied Biosystems by life technologies. [ http://www.appliedbiosystems.com/absite/us/en/home/applications-technologies/solid-next-generation-sequencing.html ]Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang HY, Landherr L, Tomsho LP, Hu Y, Carlson JE, Ma H, Schuster SC, Soltis DE, Soltis PS, Altman N, dePamphilis CW: Comparison of next generation sequencing technologies for transcriptome characterization. Bmc Genomics. 2009, 10:Murchison EP, Tovar C, Hsu A, Bender HS, Kheradpour P, Rebbeck CA, Obendorf D, Conlan C, Bahlo M, Blizzard CA, Pyecroft S, Kreiss A, Kellis M, Stark A, Harkins TT, Marshall Graves JA, Woods GM, Hanon GJ, Papenfuss AT: The Tasmanian Devil Transcriptome Reveals Schwann Cell Origins of a Clonally Transmissible Cancer. Science. 2010, 327 (5961): 84-87. 10.1126/science.1180616.Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA: Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. Bmc Genomics. 2010, 11:Babik W, Stuglik M, Qi W, Kuenzli M, Kuduk K, Koteja P, Radwan J: Heart transcriptome of the bank vole (Myodes glareolus): towards understanding the evolutionary variation in metabolic rate. BMC Genomics. 2010, 11: 390-10.1186/1471-2164-11-390.Miller JC, Tanksley SD: RFLP analysis of phylogenetic-relationships and genetic-variation in the genus Lycopersicon. Theoretical and Applied Genetics. 1990, 80 (4): 437-448.Williams CE, Stclair DA: Phenetic relationships and levels of variability detected by restriction-fragment-length-polymorphism and random amplified polymorphic DNA analysis of cultivated and wild accessions of Lycopersicon-esculentum. Genome. 1993, 36 (3): 619-630. 10.1139/g93-083.Rick CM: Tomato, Lycopersicon esculentum (Solanaceae). Evolution of crop plants. Edited by: Simmonds NW. 1976, London: Longman Group, 268-273.Labate JA, Baldo AM: Tomato SNP discovery by EST mining and resequencing. Molecular Breeding. 2005, 16 (4): 343-349. 10.1007/s11032-005-1911-5.Yano K, Watanabe M, Yamamoto N, Maeda F, Tsugane T, Shibata D: Expressed sequence tags (EST) database of a miniature tomato cultivar, Micro-Tom. Plant and Cell Physiology. 2005, 46: S139-S139.Jimenez-Gomez JM, Maloof JN: Sequence diversity in three tomato species: SNPs, markers, and molecular evolution. Bmc Plant Biology. 2009, 9:Yang WC, Bai XD, Kabelka E, Eaton C, Kamoun S, van der Knaap E, Francis D: Discovery of single nucleotide polymorphisms in Lycopersicon esculentum by computer aided analysis of expressed sequence tags. Molecular Breeding. 2004, 14 (1): 21-34.Van Deynze A, Stoffel K, Buell CR, Kozik A, Liu J, van der Knaap E, Francis D: Diversity in conserved genes in tomato. Bmc Genomics. 2007, 8:Sim SC, Robbins MD, Chilcott C, Zhu T, Francis DM: Oligonucleotide array discovery of polymorphisms in cultivated tomato (Solanum lycopersicum L.) reveals patterns of SNP variation associated with breeding. Bmc Genomics. 2009, 10:Bioinformatics at COMAV. [ http://bioinf.comav.upv.es/ngs_backbone/index.html ]Broad institute. [ http://www.broadinstitute.org/igv ]Bioinformatics at COMAV. [ http://bioinf.comav.upv.es/ngs_backbone/install.html ]Github social coding. [ http://github.com/JoseBlanca/franklin ]Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17 (12): 1093-1104. 10.1093/bioinformatics/17.12.1093.Picard. [ http://picard.sourceforge.net/index.shtml ]McKenna A, Hanna M, Banks E, Sivachenko A, Citulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research. 2010, 20: 1297-1303. 10.1101/gr.107524.110.Sol Genomics Network. [ ftp://ftp.solgenomics.net/ ]NCBI Genbank. [ http://www.ncbi.nlm.nih.gov/genbank/ ]Gundry CN, Vandersteen JG, Reed GH, Pryor RJ, Chen J, Wittwer CT: Amplicon melting analysis with labeled primers: A closed-tube method for differentiating homozygotes and heterozygotes. Clinical Chemistry. 2003, 49 (3): 396-406. 10.1373/49.3.396

Genome Wide Association Study to predict severe asthma exacerbations in children using random forests classifiers

<p>Abstract</p> <p>Background</p> <p>Personalized health-care promises tailored health-care solutions to individual patients based on their genetic background and/or environmental exposure history. To date, disease prediction has been based on a few environmental factors and/or single nucleotide polymorphisms (SNPs), while complex diseases are usually affected by many genetic and environmental factors with each factor contributing a small portion to the outcome. We hypothesized that the use of random forests classifiers to select SNPs would result in an improved predictive model of asthma exacerbations. We tested this hypothesis in a population of childhood asthmatics.</p> <p>Methods</p> <p>In this study, using emergency room visits or hospitalizations as the definition of a severe asthma exacerbation, we first identified a list of top Genome Wide Association Study (GWAS) SNPs ranked by Random Forests (RF) importance score for the CAMP (Childhood Asthma Management Program) population of 127 exacerbation cases and 290 non-exacerbation controls. We predict severe asthma exacerbations using the top 10 to 320 SNPs together with age, sex, pre-bronchodilator FEV1 percentage predicted, and treatment group.</p> <p>Results</p> <p>Testing in an independent set of the CAMP population shows that severe asthma exacerbations can be predicted with an Area Under the Curve (AUC) = 0.66 with 160-320 SNPs in comparison to an AUC score of 0.57 with 10 SNPs. Using the clinical traits alone yielded AUC score of 0.54, suggesting the phenotype is affected by genetic as well as environmental factors.</p> <p>Conclusions</p> <p>Our study shows that a random forests algorithm can effectively extract and use the information contained in a small number of samples. Random forests, and other machine learning tools, can be used with GWAS studies to integrate large numbers of predictors simultaneously.</p

The Evolution of a Female Genital Trait Widely Distributed in the Lepidoptera: Comparative Evidence for an Effect of Sexual Coevolution

Sexual coevolution is considered responsible for the evolution of many male genital traits, but its effect on female genital morphology is poorly understood. In many lepidopterans, females become temporarily unreceptive after mating and the length of this refractory period is inversely related to the amount of spermatophore remaining in their genital tracts. Sperm competition can select for males that delay female remating by transferring spermatophores with thick spermatophore envelopes that take more time to be broken. These envelopes could select for signa, sclerotized sharp structures located within the female genital tract, that are used for breaking spermatophores. Thus, this hypothesis predicts that thick spermatophore envelopes and signa evolve in polyandrous species, and that these adaptations are lost when monandry evolves subsequently. Here we test the expected associations between female mating pattern and presence/absence of signa, and review the scant information available on the thickness of spermatophore envelopes.We made a literature review and found information on female mating pattern (monandry/polyandry), presence/absence of signa and phylogenetic position for 37 taxa. We built a phylogenetic supertree for these taxa, mapped both traits on it, and tested for the predicted association by using Pagel's test for correlated evolution. We found that, as predicted by our hypothesis, monandry evolved eight times and in five of them signa were lost; preliminary evidence suggests that at least in two of the three exceptions males imposed monandry on females by means of specially thick spermatophore envelopes. Previously published data on six genera of Papilionidae is in agreement with the predicted associations between mating pattern and the characteristics of spermatophore envelopes and signa.Our results support the hypothesis that signa are a product of sexually antagonistic coevolution with spermatophore envelopes

Calpain and PARP Activation during Photoreceptor Cell Death in P23H and S334ter Rhodopsin Mutant Rats

Retinitis pigmentosa (RP) is a heterogeneous group of inherited neurodegenerative diseases affecting photoreceptors and causing blindness. Many human cases are caused by mutations in the rhodopsin gene. An important question regarding RP pathology is whether different genetic defects trigger the same or different cell death mechanisms. To answer this question, we analysed photoreceptor degeneration in P23H and S334ter transgenic rats carrying rhodopsin mutations that affect protein folding and sorting respectively. We found strong activation of calpain and poly(ADP-ribose) polymerase (PARP) in both mutants, concomitant with calpastatin down-regulation, increased oxidative DNA damage and accumulation of PAR polymers. These parameters were strictly correlated with the temporal progression of photoreceptor degeneration, mirroring earlier findings in the phosphodiesterase-6 mutant rd1 mouse, and suggesting execution of non-apoptotic cell death mechanisms. Interestingly, activation of caspases-3 and -9 and cytochrome c leakage—key events in apoptotic cell death—were observed only in the S334ter mutant, which also showed increased expression of PARP-1. The identification of the same metabolic markers triggered by different mutations in two different species suggests the existence of common cell death mechanisms, which is a major consideration for any mutation independent treatment

Hypermethylated 14-3-3-σ and ESR1 gene promoters in serum as candidate biomarkers for the diagnosis and treatment efficacy of breast cancer metastasis

Background: Numerous hypermethylated genes have been reported in breast cancer, and the silencing of these genes plays an important role in carcinogenesis, tumor progression and diagnosis. These hypermethylated promoters are very rarely found in normal breast. It has been suggested that aberrant hypermethylation may be useful as a biomarker, with implications for breast cancer etiology, diagnosis, and management. The relationship between primary neoplasm and metastasis remains largely unknown. There has been no comprehensive comparative study on the clinical usefulness of tumor-associated methylated DNA biomarkers in primary breast carcinoma and metastatic breast carcinoma. The objective of the present study was to investigate the association between clinical extension of breast cancer and methylation status of Estrogen Receptor1 (ESR1) and Stratifin (14-3-3-σ) gene promoters in disease-free and metastatic breast cancer patients. Methods: We studied two cohorts of patients: 77 patients treated for breast cancer with no signs of disease, and 34 patients with metastatic breast cancer. DNA was obtained from serum samples, and promoter methylation status was determined by using DNA bisulfite modification and quantitative methylation-specific PCR. Results: Serum levels of methylated gene promoter 14-3-3-σ significantly differed between Control and Metastatic Breast Cancer groups (P < 0.001), and between Disease-Free and Metastatic Breast Cancer groups (P < 0.001). The ratio of the 14-3-3-σ level before the first chemotherapy cycle to the level just before administration of the second chemotherapy cycle was defined as the Biomarker Response Ratio [BRR]. We calculated BRR values for the "continuous decline" and "rise-and-fall" groups. Subsequent ROC analysis showed a sensitivity of 75% (95% CI: 47.6 - 86.7) and a specificity of 66.7% (95% CI: 41.0 - 86.7) to discriminate between the groups for a cut-off level of BRR = 2.39. The area under the ROC curve (Z = 0.804 ± 0.074) indicates that this test is a good approach to post-treatment prognosis. Conclusions: The relationship of 14-3-3-σ with breast cancer metastasis and progression found in this study suggests a possible application of 14-3-3-σ as a biomarker to screen for metastasis and to follow up patients treated for metastatic breast cancer, monitoring their disease status and treatment response.This study was supported by a grant from the Ministerio de Ciencia e Innovación: SAF 2004-00889; JL Linares is supported by the Junta de Andalucía (P06-CTS-1385)

Facilitative parenting and children's social, emotional and behavioural adjustment

Facilitative parenting (FP) supports the development of children’s social and emotional competence and effective peer relationships. Previous research has shown that FP discriminates between children bullied by peers from children who are not bullied, according to reports of teachers. This study investigates the association between FP and children’s social, emotional and behavioral problems, over and above the association with dysfunctional parenting (DP). 215 parents of children aged 5–11 years completed questionnaires about parenting and child behavior, and children and teachers completed measures of child bullying victimization. As predicted, FP accounted for variance in teacher reports of children’s bullying victimization as well as parent reports of children’s social and emotional problems and prosocial behavior better than that accounted for by DP. However for children’s reports of peer victimization the whole-scale DP was a better predictor than FP. Contrary to predictions, FP accounted for variance in conduct problems and hyperactivity better than DP. When analyses were replicated substituting subscales of dysfunctional and FP, a sub-set of FP subscales including Warmth, Supports Friendships, Not Conflicting, Child Communicates and Coaches were correlated with low levels of problems on a broad range of children’s adjustment problems. Parent–child conflict accounted for unique variance in children’s peer victimization (teacher report), peer problems, depression, emotional problems, conduct problems and hyperactivity. The potential relevance of FP as a protective factor for children against a wide range of adjustment problems is discussed

Identification of a biomarker panel for colorectal cancer diagnosis

<p>Abstract</p> <p>Background</p> <p>Malignancies arising in the large bowel cause the second largest number of deaths from cancer in the Western World. Despite progresses made during the last decades, colorectal cancer remains one of the most frequent and deadly neoplasias in the western countries.</p> <p>Methods</p> <p>A genomic study of human colorectal cancer has been carried out on a total of 31 tumoral samples, corresponding to different stages of the disease, and 33 non-tumoral samples. The study was carried out by hybridisation of the tumour samples against a reference pool of non-tumoral samples using Agilent Human 1A 60-mer oligo microarrays. The results obtained were validated by qRT-PCR. In the subsequent bioinformatics analysis, gene networks by means of Bayesian classifiers, variable selection and bootstrap resampling were built. The consensus among all the induced models produced a hierarchy of dependences and, thus, of variables.</p> <p>Results</p> <p>After an exhaustive process of pre-processing to ensure data quality--lost values imputation, probes quality, data smoothing and intraclass variability filtering--the final dataset comprised a total of 8, 104 probes. Next, a supervised classification approach and data analysis was carried out to obtain the most relevant genes. Two of them are directly involved in cancer progression and in particular in colorectal cancer. Finally, a supervised classifier was induced to classify new unseen samples.</p> <p>Conclusions</p> <p>We have developed a tentative model for the diagnosis of colorectal cancer based on a biomarker panel. Our results indicate that the gene profile described herein can discriminate between non-cancerous and cancerous samples with 94.45% accuracy using different supervised classifiers (AUC values in the range of 0.997 and 0.955).</p

Validation of the Spanish Version of the ICECAP-O for Nursing Home Residents with Dementia

Background Measurement of health-related quality of life (HRQoL) is important for a chronic disease, such as dementia, which impairs the quality of life of affected patients in addition to their length of life. This is important in the context of economic evaluations when interventions do not (only) affect HRQoL and these other factors also affect overall quality of life. Objective To validate the Spanish translation of the ICECAP-O's capability to measure Health-related quality of life in elderly with dementia who live in nursing homes. Method Cross-sectional study. For 217 residents living in 8 Spanish nursing homes, questionnaires were completed by nursing professionals serving as proxy respondents. We analyzed the internal consistency and other psychometric properties. We investigated the convergent validity of the ICECAP-O with other HRQoL instruments, the EQ-5D extended with a cognitive dimension (EQ-5D+C), the Alzheimer's Disease Related Quality of Life (ADRQL) measures, and the Barthel Index measure of activities of daily living (ADL). Results The ICECAP-O presents satisfactory internal consistency (alpha 0.820). The factorial analysis indicated a structure of five principal dimensions that explain 66.57% of the total variance. Convergent validity between the ICECAP-O, EQ-5D+C, ADRQL, and Barthel Index scores was moderate to good (with correlations of 0.62, 0.61, and 0.68, respectively), but differed between dimensions of the instruments. Discriminant validity was confirmed by finding differences in ICECAP-O scores between subgroups based on ADL scores (0.70 low, 0.59 medium, and 0.39 high level care), dementia severity (0.72 mild, 0.63 medium, and 0.50 severe), and ages (0.59 below 75 years and 0.84 above 75 years). Conclusions This study presented the first use of a Spanish version of the ICECAP-O. The results indicate that the ICECAP-O appears to be a reliable Health-related quality of life measurement instrument showing good convergent and discriminant validity for people with dementia