In silico analysis of Simple Sequence Repeats from chloroplast genomes of Solanaceae species

The availaibilty of chloroplast genome (cpDNA) sequences of Atropa belladonna, Nicotiana sylvestris, N tabacum, N tomentosiformis, Solanum bulbocastanum, S lycopersicum and S tuberosum, which are Solanaceae species, allowed us to analyze the organization of cpSSRs in their genic and intergenic regions In general, the number of cpSSRs in cpDNA ranged from 161 in S tuberosum to 226 in N tabacum, and the number of intergenic cpSSRs was higher than genic cpSSRs The mononucleotide repeats were the most frequent in studied species, but we also identified di-, tri-, tetra-, penta- and hexanucleotide repeats Multiple alignments of all cpSSRs sequence from Solanaceae species made the identification of nucleotide variability possible and the phylogeny was estimated by maximum parsimony Our study showed that the plastome database can be exploited for phylogenetic analyses and biotechnological approaches

In silico analysis of candidate genes involved in light sensing and signal transduction pathways in soybean

Several aspects of photoperception and light signal transduction have been elucidated by studies with model plants. However, the information available for economically important crops, such as Fabaceae species, is scarce. In order to incorporate the existing genomic tools into a strategy to advance soybean research, we have investigated publicly available expressed sequence tag ( EST) sequence databases in order to identify Glycine max sequences related to genes involved in light-regulated developmental control in model plants. Approximately 38,000 sequences from open-access databases were investigated, and all bona fide and putative photoreceptor gene families were found in soybean sequence databases. We have identified G. max orthologs for several families of transcriptional regulators and cytoplasmic proteins mediating photoreceptor-induced responses, although some important Arabidopsis phytochrome-signaling components are absent. Moreover, soybean and Arabidopsis gene-family homologs appear to have undergone a distinct expansion process in some cases. We propose a working model of light perception, signal transduction and response-eliciting in G. max, based on the identified key components from Arabidopsis. These results demonstrate the power of comparative genomics between model systems and crop species to elucidate several aspects of plant physiology and metabolism.

In silico analysis of microRNAS targeting the HLA-G 3 ` untranslated region alleles and haplotypes

It has been reported that microRNAs (miRNA) may have allele-specific targeting for the 3` untranslated region (3` UTR) of the HLA-G locus. In a previous study, we reported 11 3`UTR haplotypes encompassing the 14-bp insertion/deletion polymorphism and seven SNPs (+3003 T/C, +3010 C/G, +3027 C/A, +3035 C/T, +3142 C/G, +3187A/G,and +3196 C/G), of which only the +3142 C/G SNP has been reported to influence the binding of miRNAs. Using bioinformatics analyses, we identified putative miRNA-binding sites considering the haplotypes encompassing these eight polymorphic sites, and we ranked the lowest free energies that could potentially lead to an mRNA degradation or translational repression. When a specific haplotype or a particular SNP was associated with a miRNA-binding site, we defined a free energy difference of 4 kcal/mol between alleles to classify them energetically distant. The best results were obtained for the miR-513a-5p, miR-518c*, miR-1262 and miR-92a-1*, miR-92a-2*, miR-661, miR-1224-5p, and miR-433 miRNAs, all influencing one or more of the +3003, +3010, +3027, and +3035 SNPs. The miR-2110, miR-93, miR-508-5p, miR-331-5p, miR-616, miR-513b, and miR-589* miRNAs targeted the 14-bp fragment region, and miR-148a, miR-19a*, miR-152, mir-148b,and miR-218-2 also influenced the +3142C/G polymorphism. These results suggest that these miRNAs might play a relevant role on the HLA-G expression pattern. (C) 2009 Published by Elsevier Inc. on behalf of American Society for Histocompatibility and Immunogenetics.

PK/DB: database for pharmacokinetic properties and predictive in silico ADME models

The study of pharmacokinetic properties (PK) is of great importance in drug discovery and development. In the present work, PK/DB (a new freely available database for PK) was designed with the aim of creating robust databases for pharmacokinetic studies and in silico absorption, distribution, metabolism and excretion (ADME) prediction. Comprehensive, web-based and easy to access, PK/DB manages 1203 compounds which represent 2973 pharmacokinetic measurements, including five models for in silico ADME prediction (human intestinal absorption, human oral bioavailability, plasma protein binding, bloodbrain barrier and water solubility).

In silico screening of HIV-1 non-nucleoside reverse transcriptase and protease inhibitors

Two targets, reverse transcriptase (RT) and protease from HIV-1, were used during the past two decades to the discovery of non-nucleoside reverse transcriptase inhibitors (NNRTI) and protease inhibitors (PI) that belong to the arsenal of the antiretroviral therapy. Herein these enzymes were chosen as templates for conducting a computer-aided ligand design. Ligand and structure-based drug designs were the starting points to select compounds from a database bearing more than five million compounds by means of cheminformatic tools. New promising lead structures are retrieved from the database, which are open to acquisition and test. Classes of molecules already described as NNRTI or PI in the literature also came out and were useful to prove the reliability of the workflow, and thus validating the work carried out so far. (c) 2007 Elsevier Masson SAS. All rights reserved.

High-throughput SNP genotyping in the highly heterozygous genome of Eucalyptus: assay success, polymorphism and transferability across species

Background: High-throughput SNP genotyping has become an essential requirement for molecular breeding and population genomics studies in plant species. Large scale SNP developments have been reported for several mainstream crops. A growing interest now exists to expand the speed and resolution of genetic analysis to outbred species with highly heterozygous genomes. When nucleotide diversity is high, a refined diagnosis of the target SNP sequence context is needed to convert queried SNPs into high-quality genotypes using the Golden Gate Genotyping Technology (GGGT). This issue becomes exacerbated when attempting to transfer SNPs across species, a scarcely explored topic in plants, and likely to become significant for population genomics and inter specific breeding applications in less domesticated and less funded plant genera. Results: We have successfully developed the first set of 768 SNPs assayed by the GGGT for the highly heterozygous genome of Eucalyptus from a mixed Sanger/454 database with 1,164,695 ESTs and the preliminary 4.5X draft genome sequence for E. grandis. A systematic assessment of in silico SNP filtering requirements showed that stringent constraints on the SNP surrounding sequences have a significant impact on SNP genotyping performance and polymorphism. SNP assay success was high for the 288 SNPs selected with more rigorous in silico constraints; 93% of them provided high quality genotype calls and 71% of them were polymorphic in a diverse panel of 96 individuals of five different species. SNP reliability was high across nine Eucalyptus species belonging to three sections within subgenus Symphomyrtus and still satisfactory across species of two additional subgenera, although polymorphism declined as phylogenetic distance increased. Conclusions: This study indicates that the GGGT performs well both within and across species of Eucalyptus notwithstanding its nucleotide diversity >= 2%. The development of a much larger array of informative SNPs across multiple Eucalyptus species is feasible, although strongly dependent on having a representative and sufficiently deep collection of sequences from many individuals of each target species. A higher density SNP platform will be instrumental to undertake genome-wide phylogenetic and population genomics studies and to implement molecular breeding by Genomic Selection in Eucalyptus.

Expression and cellular localization of microRNA-29b and RAX, an activator of the RNA-dependent protein kinase (PKR), in the retina of streptozotocin-induced diabetic rats

Purpose: The apoptosis of retinal neurons plays a critical role in the pathogenesis of diabetic retinopathy (DR), but the molecular mechanisms underlying this phenomenon remain unclear. The purpose of this study was to investigate the cellular localization and the expression of microRNA-29b (miR-29b) and its potential target PKR associated protein X (RAX), an activator of the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway, in the retina of normal and diabetic rats. Methods: Retinas were obtained from normal and diabetic rats within 35 days after streptozotocin (STZ) injection. In silico analysis indicated that RAX is a potential target of miR-29b. The cellular localization of miR-29b and RAX was assessed by in situ hybridization and immunofluorescence, respectively. The expression levels of miR-29b and RAX mRNA were evaluated by quantitative reverse transcription PCR (qRT-PCR), and the expression of RAX protein was evaluated by western blot. A luciferase reporter assay and inhibition of endogenous RAX were performed to confirm whether RAX is a direct target of miR-29b as predicted by the in silico analysis. Results: We found that miR-29b and RAX are localized in the retinal ganglion cells (RGCs) and the cells of the inner nuclear layer (INL) of the retinas from normal and diabetic rats. Thus, the expression of miR-29b and RAX, as assessed in the retina by quantitative RT-PCR, reflects their expression in the RGCs and the cells of the INL. We also revealed that RAX protein is upregulated (more than twofold) at 3, 6, 16, and 22 days and downregulated (70%) at 35 days, whereas miR-29b is upregulated (more than threefold) at 28 and 35 days after STZ injection. We did not confirm the computational prediction that RAX is a direct target of miR-29b. Conclusions: Our results suggest that RAX expression may be indirectly regulated by miR-29b, and the upregulation of this miRNA at the early stage of STZ-induced diabetes may have a protective effect against the apoptosis of RGCs and cells of the INL by the pro-apoptotic RNA-dependent protein kinase (PKR) signaling pathway.

Global gene expression under nitrogen starvation in Xylella fastidiosa: contribution of the sigma(54) regulon

Background: Xylella fastidiosa, a Gram-negative fastidious bacterium, grows in the xylem of several plants causing diseases such as citrus variegated chlorosis. As the xylem sap contains low concentrations of amino acids and other compounds, X. fastidiosa needs to cope with nitrogen limitation in its natural habitat. Results: In this work, we performed a whole-genome microarray analysis of the X. fastidiosa nitrogen starvation response. A time course experiment (2, 8 and 12 hours) of cultures grown in defined medium under nitrogen starvation revealed many differentially expressed genes, such as those related to transport, nitrogen assimilation, amino acid biosynthesis, transcriptional regulation, and many genes encoding hypothetical proteins. In addition, a decrease in the expression levels of many genes involved in carbon metabolism and energy generation pathways was also observed. Comparison of gene expression profiles between the wild type strain and the rpoN null mutant allowed the identification of genes directly or indirectly induced by nitrogen starvation in a sigma(54)-dependent manner. A more complete picture of the sigma(54) regulon was achieved by combining the transcriptome data with an in silico search for potential sigma(54)-dependent promoters, using a position weight matrix approach. One of these sigma(54)-predicted binding sites, located upstream of the glnA gene (encoding glutamine synthetase), was validated by primer extension assays, confirming that this gene has a sigma(54)-dependent promoter. Conclusions: Together, these results show that nitrogen starvation causes intense changes in the X. fastidiosa transcriptome and some of these differentially expressed genes belong to the sigma(54) regulon.

Novel Primate-Specific Genes, RMEL 1, 2 and 3, with Highly Restricted Expression in Melanoma, Assessed by New Data Mining Tool

Melanoma is a highly aggressive and therapy resistant tumor for which the identification of specific markers and therapeutic targets is highly desirable. We describe here the development and use of a bioinformatic pipeline tool, made publicly available under the name of EST2TSE, for the in silico detection of candidate genes with tissue-specific expression. Using this tool we mined the human EST (Expressed Sequence Tag) database for sequences derived exclusively from melanoma. We found 29 UniGene clusters of multiple ESTs with the potential to predict novel genes with melanoma-specific expression. Using a diverse panel of human tissues and cell lines, we validated the expression of a subset of three previously uncharacterized genes (clusters Hs.295012, Hs.518391, and Hs.559350) to be highly restricted to melanoma/melanocytes and named them RMEL1, 2 and 3, respectively. Expression analysis in nevi, primary melanomas, and metastatic melanomas revealed RMEL1 as a novel melanocytic lineage-specific gene up-regulated during melanoma development. RMEL2 expression was restricted to melanoma tissues and glioblastoma. RMEL3 showed strong up-regulation in nevi and was lost in metastatic tumors. Interestingly, we found correlations of RMEL2 and RMEL3 expression with improved patient outcome, suggesting tumor and/or metastasis suppressor functions for these genes. The three genes are composed of multiple exons and map to 2q12.2, 1q25.3, and 5q11.2, respectively. They are well conserved throughout primates, but not other genomes, and were predicted as having no coding potential, although primate-conserved and human-specific short ORFs could be found. Hairpin RNA secondary structures were also predicted. Concluding, this work offers new melanoma-specific genes for future validation as prognostic markers or as targets for the development of therapeutic strategies to treat melanoma.

Characterization of new IS elements and studies of their dispersion in two subspecies of Leifsonia xyli

Background: Leifsonia xyli is a xylem-inhabiting bacterial species comprised of two subspecies: L. xyli subsp. xyli (Lxx) and L. xyli subsp. cynodontis (Lxc). Lxx is the causal agent of ratoon stunting disease in sugarcane commercial fields and Lxc colonizes the xylem of several grasses causing either mild or no symptoms of disease. The completely sequenced genome of Lxx provided insights into its biology and pathogenicity. Since IS elements are largely reported as an important source of bacterial genome diversification and nothing is known about their role in chromosome architecture of L. xyli, a comparative analysis of Lxc and Lxx elements was performed. Results: Sample sequencing of Lxc genome and comparative analysis with Lxx complete DNA sequence revealed a variable number of IS transposable elements acting upon genomic diversity. A detailed characterization of Lxc IS elements and a comparative review with IS elements of Lxx are presented. Each genome showed a unique set of elements although related to same IS families when considering features such as similarity among transposases, inverted and direct repeats, and element size. Most of the Lxc and Lxx IS families assigned were reported to maintain transposition at low levels using translation regulatory mechanisms, consistent with our in silico analysis. Some of the IS elements were found associated with rearrangements and specific regions of each genome. Differences were also found in the effect of IS elements upon insertion, although none of the elements were preferentially associated with gene disruption. A survey of transposases among genomes of Actinobacteria showed no correlation between phylogenetic relatedness and distribution of IS families. By using Southern hybridization, we suggested that diversification of Lxc isolates is also mediated by insertion sequences in probably recent events. Conclusion: Collectively our data indicate that transposable elements are involved in genome diversification of Lxc and Lxx. The IS elements were probably acquired after the divergence of the two subspecies and are associated with genome organization and gene contents. In addition to enhancing understanding of IS element dynamics in general, these data will contribute to our ongoing comparative analyses aimed at understanding the biological differences of the Lxc and Lxx.

HFE gene mutations in patients with primary iron overload: Is there a significant improvement in molecular diagnosis yield with HFE sequencing?

Rare HFE variants have been shown to be associated with hereditary hemochromatosis (HH), an iron overload disease. The low frequency of the HFE p.C282Y mutation in HH-affected Brazilian patients may suggest that other HFE-related mutations may also be implicated in the pathogenesis of HH in this population. The main aim was to screen for new HFE mutations in Brazilian individuals with primary iron overload and to investigate their relationship with HH. Fifty Brazilian patients with primary iron overload (transferrin saturation >50% in females and 60% in males) were selected. Subsequent bidirectional sequencing for each HFE exon was performed. The effect of HFE mutations on protein structure were analyzed by molecular dynamics simulation and free binding energy calculations. p.C282Y in homozygosis or in heterozygosis with p.H63D were the most frequent genotypic combinations associated with HH in our sample population (present in 17 individuals, 34%). Thirty-six (72.0%) out of the 50 individuals presented at least one HFE mutation. The most frequent genotype associated with HH was the homozygous p.C282Y mutation (n = 11, 22.0%). One novel mutation (p.V256I) was indentified in heterozygosis with the p.H63D mutation. In silico modeling analysis of protein behavior indicated that the p.V256I mutation does not reduce the binding affinity between HFE and beta 2-microglobulin ((beta 2M) in the same way the p.C282Y mutation does compared with the native HFE protein. In conclusion, screening of HFE through direct sequencing, as compared to p.C282Y/p.H63D genotyping, was not able to increase the molecular diagnosis yield of HH. The novel p.V256I mutation could not be implicated in the molecular basis of the HH phenotype, although its role cannot be completely excluded in HH-phenotype development. Our molecular modeling analysis can help in the analysis of novel, previously undescribed, HFE mutations. (C) 2010 Elsevier Inc. All rights reserved.

Intrinsically bent DNA sites in the Drosophila melanogaster third chromosome amplified domain

Bent DNA sites promote the curvature of DNA in both eukaryotic and prokaryotic chromosomes. Here, we investigate the localization and structure of intrinsically bent DNA sites in the extensively characterized Drosophila melanogaster third chromosome DAFC-66D segment (Drosophila amplicon in the follicle cells). This region contains the amplification control element ACE3, which is a replication enhancer that acts in cis to activate the major replication origin ori-beta. Through both electrophoretic and in silico analysis, we have identified three major bent DNA sites in DAFC-66D. The bent DNA site (b1) is localized in the ACE3 element, whereas the other two bent DNA sites (b2 and b3) are localized in the ori-beta region. Four additional bent DNA sites were identified in the intron of the S18 gene and near the TATA box of the S15, S19, and S16 genes. The identification of DNA bent sites in genomic regions previously characterized as functionally relevant for DNA amplification further supports a function for DNA bent sites in DNA replication in eukaryotes.

Evaluation of RET polymorphisms in a six-generation family with G533C RET mutation: specific RET variants may modulate age at onset and clinical presentation

P>Context We previously described a six-generation family with G533C RET mutation and medullary thyroid carcinoma, in the largest family reported do date. Of particular interest, phenotype variability regarding the age of onset and clinical presentation of the disease, was observed. Objective We evaluate whether single SNPs within RET oncogene or haplotype comprising the RET variants (defined by Haploview) could predispose to early development of MTC in this family and influence the clinical manifestation. Design Eight SNPs were selected based on their previous association with the clinical course of hereditary or sporadic MTC, in particular promoting an early onset of disease. The variants were initially tested in 77 G533C-carriers and 100 controls using either PCR-direct sequencing or PCR-RFLP. Association between a SNP or haplotype and age at diagnosis or presence of lymph node metastasis was tested in 34 G533C-carries with MTC. Different bioinformatic tools were used to evaluate the potential effects on RNA splicing. Results An association was found between IVS1-126G > T and age at diagnosis. The variant [IVS8 +82A > G; 85-86 insC] was associated with the presence of lymph node metastases at diagnosis. In silico analysis suggested that this variant may induce abnormal splicing. This in silico analysis predicted that the [IVS8 +82A > G; 85-86 insC] could alter the splicing by disrupting and/or creating exonic splicing enhancer motifs. Conclusions We here identified two RET variants that were associated with phenotype variability in G533C-carriers, which highlights the fact that the modifier effect of a variant might depend on the type of mutation.

Messenger RNA expression of Pabpnl1 and Mbd3l2 genes in oocytes and cleavage embryos

Objective: To identify genes specifically expressed in mammalian oocytes using an in silico subtraction, and to characterize the mRNA patterns of selected genes in oocytes, embryos, and adult tissues. Design: Comparison between oocyte groups and between early embryo stages. Setting: Laboratories of embryo manipulation and molecular biology from Departamento de Genetica (FMRP) and Departamento de Ciencias Basicas (FZEA) - University of Sao Paulo. Sample(s): Oocytes were collected from slaughtered cows for measurements, in vitro fertilization, and in vitro embryo culture. Somatic tissue, excluding gonad and uterus tissue, was collected from male and female cattle. Main Outcome Measure(s): Messenger RNA levels of poly(A)-binding protein nuclear-like 1 (Pabpnl1) and methyl-CpG-binding domain protein 3-like 2 (Mbd3l2). Result(s): Pabpnl1 mRNA was found to be expressed in oocytes, and Mbd3l2 transcripts were present in embryos. Quantification of Pabpnl1 transcripts showed no difference in levels between good-and bad-quality oocytes before in vitro maturation (IVM) or between good-quality oocytes before and after IVM. However, Pabpnl1 transcripts were not detected in bad-quality oocytes after IVM. Transcripts of the Mbd3l2 gene were found in 4-cell, 8-cell, and morula-stage embryos, with the highest level observed in 8-cell embryos. Conclusion(s): Pabpnl1 gene expression is restricted to oocytes and Mbd3l2 to embryos. Different Pabpnl1 mRNA levels in oocytes of varying viability suggest an important role in fertility involving the oocyte potential for embryo development. (Fertil Steril (R) 2010; 93: 2507-12. (C) 2010 by American Society for Reproductive Medicine.)

Clinical evaluation and COL2A1 gene analysis in 21 Brazilian families with Stickler syndrome: Identification of novel mutations, further genotype/phenotype correlation, and its implications for the diagnosis

We present clinical and molecular evaluation from a large cohort of patients with Stickler syndrome: 78 individuals from 21 unrelated Brazilian families. The patients were selected in a Hospital with a craniofacial dysmorphology assistance service and clinical diagnosis was based on the presence of cleft palate associated to facial and ocular anomalies of Stickler syndrome. Analysis of COL2A1 gene revealed 9 novel and 4 previously described pathogenic mutations. Except for the mutation c.556G>T (p.Gly186X), all the others were located in the triple helical domain. We did not find genotype/phenotype correlation in relation to type and position of the mutation in the triple helical domain. However, a significantly higher proportion of myopia in patients with mutations located in this domain was observed in relation to those with the mutation in the non-tripe helical domain (c.556G>T; P < 0.04). A trend towards a higher prevalence of glaucoma, although not statistically significant, was observed in the presence of the mutation c.556G>T. It is possible. that this mutation alters the splicing of the mRNA instead of only creating a premature stop codon and therefore it can lead to protein products of different ocular effects. One novel DNA variation (c.1266+7G>C) occurs near a splice site and it was observed to co-segregate with the phenotype in one of the two families with this DNA variation. As in silico analysis predicted that the c.1266+7G>C DNA variation can affect the efficiency of the splicing, we still cannot rule it out as non-pathogenic. Our study also showed that ascertainment through cleft palate associated to other craniofacial signs can be very efficient for identification of Stickler syndrome patients. Still, high frequency of familial cases and high frequency of underdevelopment of distal lateral tibial epiphyses observed in our patients suggested that the inclusion of this information can improve the clinical diagnosis of Stickler syndrome. (C) 2008 Elsevier Masson SAS. All rights reserved.