2LD, GENECOUNTING and HAP: Computer programs for linkage disequilibrium analysis

Computer programs are introduced which calculate pair-wise linkage disequilibrium statistics and conduct haplotype frequency estimation, including X chromosome data, and using a heuristic algorithm to handle multiple genetic markers and missing data. AVAILABILITY: Programs 2LD, GENECOUNTING and HAP are available on Internet from http://www.hgmp.mrc.ac.uk/~jzhao and http://www.iop.kcl.ac.uk/IoP/Departments/PsychMed/GEpiBSt/software.shtml     

Genetic profiling of the critically endangered palm species Medemia argun using newly developed chloroplast DNA markers

Background: Medemia argun is a rare wild palm tree species. Its global existence is assumed to include the main population of about 1000 trees in the Nubian Desert of Sudan and some scattered individuals in southern Egypt. The species had previously been assumed to be extinct, but then reported to be extant about 20 years ago.

Aims: To assess genetic variation and explore population genetic structure of M. argun, through development and analysis of microsatellite markers.

Methods: The genome sequence mining approach was applied in order to identify microsatellites in the chloroplast genome of Bismarckia nobilis, a species closely related of M. argun. A set of 49 markers were designed, and their characteristics are now provided. Seven chloroplast DNA markers were developed for use in the genetic characterisation of this threatened species.

Results: Five markers were found polymorphic in M. argun, which enabled the assessment of the genetic diversity of the species. Significant genetic differentiation was observed among generations and collection sites, accompanied by low genetic variation. The seven markers developed were polymorphic among the wild relatives Hyphaene thebaica and Borassus aethiopum.

Conclusions: This is the first study to report molecular markers for M. argun. Our results suggest that the genetic consequences of population fragmentation in M. argun are beginning to be evident.     

<Emphasis Type="Italic">snpReady</Emphasis>: a tool to assist breeders in genomic analysis

The snpReady R package is a new instrument developed to help breeders in genomic projects such as genomic prediction and association studies. This package offers three different methods to build the genomic relationship matrix, a new imputation method for missing markers based on Wright’s theory, and a population genetic overview. Therefore, we implemented three functions (raw.data, G.matrix, and popgen). Hence, this tool allows the raw data to be transformed from different genotyping platforms to numeric matrices and performs quality control (missing data and allele frequency). Moreover, the package generates and exports four different relationship matrices (proposed by Yang et al. (N 569:565–569, 2010), VanRaden (JDS 91:4414–23, 2008), and the Gaussian kernel) depending on the purpose and software to be used in further analysis. Finally, based on the genotypic matrix, the package estimates the genetic variability, effective population size, and endogamy, among other population genetic parameters. Empirical comparisons between the method of imputation proposed and other well-known approaches have shown a lower accuracy of imputation, however, with no significant impact on the genome prediction accuracies when a lower amount of missing data is allowed. The functions and arguments were designed to carry out the preparation of genomic datasets in a straightforward, fast, and more computationally efficient way.The package and its details are available at CRAN or http://www.github.com/italo-granato/snpReady.     

Single-nucleotide polymorphism,linkage disequilibrium and geographic structure in the malaria parasite <Emphasis Type="Italic">Plasmodium vivax</Emphasis>: prospects for genome-wide association studies

Background

The ideal malaria parasite populations for initial mapping of genomic regions contributing to phenotypes such as drug resistance and virulence, through genome-wide association studies, are those with high genetic diversity, allowing for numerous informative markers, and rare meiotic recombination, allowing for strong linkage disequilibrium (LD) between markers and phenotype-determining loci. However, levels of genetic diversity and LD in field populations of the major human malaria parasite P. vivax remain little characterized.

Results

We examined single-nucleotide polymorphisms (SNPs) and LD patterns across a 100-kb chromosome segment of P. vivax in 238 field isolates from areas of low to moderate malaria endemicity in South America and Asia, where LD tends to be more extensive than in holoendemic populations, and in two monkey-adapted strains (Salvador-I, from El Salvador, and Belem, from Brazil). We found varying levels of SNP diversity and LD across populations, with the highest diversity and strongest LD in the area of lowest malaria transmission. We found several clusters of contiguous markers with rare meiotic recombination and characterized a relatively conserved haplotype structure among populations, suggesting the existence of recombination hotspots in the genome region analyzed. Both silent and nonsynonymous SNPs revealed substantial between-population differentiation, which accounted for ~40% of the overall genetic diversity observed. Although parasites clustered according to their continental origin, we found evidence for substructure within the Brazilian population of P. vivax. We also explored between-population differentiation patterns revealed by loci putatively affected by natural selection and found marked geographic variation in frequencies of nucleotide substitutions at the pvmdr-1 locus, putatively associated with drug resistance.

Conclusion

These findings support the feasibility of genome-wide association studies in carefully selected populations of P. vivax, using relatively low densities of markers, but underscore the risk of false positives caused by population structure at both local and regional levels.See commentary: http://www.biomedcentral.com/1741-7007/8/90

     

ReCount: A multi-experiment resource of analysis-ready RNA-seq gene count datasets

Background

With next-generation sequencing technologies, experiments that were considered prohibitive only a few years ago are now possible. However, while these technologies have the ability to produce enormous volumes of data, the sequence reads are prone to error. This poses fundamental hurdles when genetic diversity is investigated.

Results

We developed ShoRAH, a computational method for quantifying genetic diversity in a mixed sample and for identifying the individual clones in the population, while accounting for sequencing errors. The software was run on simulated data and on real data obtained in wet lab experiments to assess its reliability.

Conclusions

ShoRAH is implemented in C++, Python, and Perl and has been tested under Linux and Mac OS X. Source code is available under the GNU General Public License at http://www.cbg.ethz.ch/software/shorah.     

PlantGM: a database for genetic markers in rice (Oryza sativa) and Chinese cabbage (Brassica rapa)

The Plant Genetic Map Database (PlantGM) has been developed as a web-based system which provides information about genetic markers in rice (Oryza sativa) and Chinese cabbage (Brassica rapa). The database has three major parts and functions; (1) Map Search, (2) Marker Search, and (3) QTL Search. At present, the database provides characterization information for about 3258 genetic markers. It has 2800 RFLP and 112 QTL markers related to rice in addition to 321 RFLP and 25 PCR-based markers for Chinese cabbage. In addition, a genetic linkage map was also constructed by using 1,054 markers from 2,912 markers in rice.

Availability

The database is available for free at http://www.niab.go.kr/nabic/PlantGM     

A physical map of the Chinese chestnut (Castanea mollissima) genome and its integration with the genetic map

Three Chinese chestnut bacterial artificial chromosome (BAC) libraries were developed and used for physical map construction. Specifically, high information content fingerprinting was used to assemble 126,445 BAC clones into 1,377 contigs and 12,919 singletons. Integration of the dense Chinese chestnut genetic map with the physical map was achieved via high-throughput hybridization using overgo probes derived from sequence-based genetic markers. A total of 1,026 probes were anchored to the physical map including 831 probes corresponding to 878 expressed sequence tag-based markers. Within the physical map, three BAC contigs were anchored to the three major fungal blight-resistant quantitative trait loci on chestnut linkage groups B, F, and G. A subset of probes corresponding to orthologous genes in poplar showed only a limited amount of conserved gene order between the poplar and chestnut genomes. The integrated genetic and physical map of Chinese chestnut is available at www.fagaceae.org/physical_maps.     

Estimating Individual Admixture Proportions from Next Generation Sequencing Data

Inference of population structure and individual ancestry is important both for population genetics and for association studies. With next generation sequencing technologies it is possible to obtain genetic data for all accessible genetic variations in the genome. Existing methods for admixture analysis rely on known genotypes. However, individual genotypes cannot be inferred from low-depth sequencing data without introducing errors. This article presents a new method for inferring an individual’s ancestry that takes the uncertainty introduced in next generation sequencing data into account. This is achieved by working directly with genotype likelihoods that contain all relevant information of the unobserved genotypes. Using simulations as well as publicly available sequencing data, we demonstrate that the presented method has great accuracy even for very low-depth data. At the same time, we demonstrate that applying existing methods to genotypes called from the same data can introduce severe biases. The presented method is implemented in the NGSadmix software available at http://www.popgen.dk/software.     

BioTorrents: A File Sharing Service for Scientific Data

The transfer of scientific data has emerged as a significant challenge, as datasets continue to grow in size and demand for open access sharing increases. Current methods for file transfer do not scale well for large files and can cause long transfer times. In this study we present BioTorrents, a website that allows open access sharing of scientific data and uses the popular BitTorrent peer-to-peer file sharing technology. BioTorrents allows files to be transferred rapidly due to the sharing of bandwidth across multiple institutions and provides more reliable file transfers due to the built-in error checking of the file sharing technology. BioTorrents contains multiple features, including keyword searching, category browsing, RSS feeds, torrent comments, and a discussion forum. BioTorrents is available at http://www.biotorrents.net.     

Contextual Cross-Referencing of Species Names for Fiddler Crabs (Genus Uca): An Experiment in Cyber-Taxonomy

Cyber-taxonomy of name usage has focused primarily on producing authoritative lists of names or cross-linking names and data across disparate databases. A feature missing from much of this work is the recording and analysis of the context in which a name was used—context which can be critical for understanding not only what name an author used, but to which currently recognized species they actually refer. An experiment on recording contextual information associated with name usage was conducted for the fiddler crabs (genus Uca). Data from approximately one quarter of all publications that mention fiddler crabs, including 95% of those published prior to 1924 and 67% of those published prior to 1976, have currently been recorded in a database. Approaches and difficulties in recording and analyzing the context of name use are discussed. These results are not meant to be a full solution, rather to highlight problems which have not been previously investigated and may act as a springboard for broader approaches and discussion. Some data on the accessibility of the literature, including in particular electronic forms of publication, are also presented. The resulting data has been integrated for general browsing into the website http://www.fiddlercrab.info; the raw data and code used to construct the website is available at https://github.com/msrosenberg/fiddlercrab.info.     

lociNGS: A Lightweight Alternative for Assessing Suitability of Next-Generation Loci for Evolutionary Analysis

Genomic enrichment methods and next-generation sequencing produce uneven coverage for the portions of the genome (the loci) they target; this information is essential for ascertaining the suitability of each locus for further analysis. lociNGS is a user-friendly accessory program that takes multi-FASTA formatted loci, next-generation sequence alignments and demographic data as input and collates, displays and outputs information about the data. Summary information includes the parameters coverage per locus, coverage per individual and number of polymorphic sites, among others. The program can output the raw sequences used to call loci from next-generation sequencing data. lociNGS also reformats subsets of loci in three commonly used formats for multi-locus phylogeographic and population genetics analyses – NEXUS, IMa2 and Migrate. lociNGS is available at https://github.com/SHird/lociNGS and is dependent on installation of MongoDB (freely available at http://www.mongodb.org/downloads). lociNGS is written in Python and is supported on MacOSX and Unix; it is distributed under a GNU General Public License.     

Genomic Target Database (GTD): A database of potential targets in human pathogenic bacteria

A Genomic Target Database (GTD) has been developed having putative genomic drug targets for human bacterial pathogens. The selected pathogens are either drug resistant or vaccines are yet to be developed against them. The drug targets have been identified using subtractive genomics approaches and these are subsequently classified into

1. Drug targets in pathogen specific unique metabolic pathways,
2. Drug targets in host-pathogen common metabolic pathways, and
3. Membrane localized drug targets.

HTML code is used to link each target to its various properties and other available public resources. Essential resources and tools for subtractive genomic analysis, sub-cellular localization, vaccine and drug designing are also mentioned. To the best of authors knowledge, no such database (DB) is presently available that has listed metabolic pathways and membrane specific genomic drug targets based on subtractive genomics. Listed targets in GTD are readily available resource in developing drug and vaccine against the respective pathogen, its subtypes, and other family members. Currently GTD contains 58 drug targets for four pathogens. Shortly, drug targets for six more pathogens will be listed.

Availability

GTD is available at IIOAB website http://www.iioab.webs.com/GTD.htm. It can also be accessed at http://www.iioabdgd.webs.com.GTD is free for academic research and non-commercial use only. Commercial use is strictly prohibited without prior permission from IIOAB.     

A database for human Y chromosome protein data

The human Y chromosome is the sex determining chromosome. The number of proteins associated with this chromosome is 196 and 107 of the 196 proteins have yet not been characterised. Here, we describe the analysis of these 107 proteins by computing various physicochemical properties using sequence and predicted structural data to elucidate molecular function. We present the derived data in the form a form a database made freely available for download, review, refinement and update.

Availability

http://puratham.googlepages.com/ or http://puratham.googlepages.com/ftpconnection     

Pseudosynapsis and Decreased Stringency of Meiotic Repair Pathway Choice on the Hemizygous Sex Chromosome of Caenorhabditis elegans Males

During meiosis, accurate chromosome segregation relies on homology to mediate chromosome pairing, synapsis, and crossover recombination. Crossovers are dependent upon formation and repair of double-strand breaks (DSBs) by homologous recombination (HR). In males of many species, sex chromosomes are largely hemizygous, yet DSBs are induced along nonhomologous regions. Here we analyzed the genetic requirements for meiotic DSB repair on the completely hemizygous X chromosome of Caenorhabditis elegans males. Our data reveal that the kinetics of DSB formation, chromosome pairing, and synapsis are tightly linked in the male germ line. Moreover, DSB induction on the X is concomitant with a brief period of pseudosynapsis that may allow X sister chromatids to masquerade as homologs. Consistent with this, neither meiotic kleisins nor the SMC-5/6 complex are essential for DSB repair on the X. Furthermore, early processing of X DSBs is dependent on the CtIP/Sae2 homolog COM-1, suggesting that as with paired chromosomes, HR is the preferred pathway. In contrast, the X chromosome is refractory to feedback mechanisms that ensure crossover formation on autosomes. Surprisingly, neither RAD-54 nor BRC-2 are essential for DSB repair on the X, suggesting that unlike autosomes, the X is competent for repair in the absence of HR. When both RAD-54 and the structure-specific nuclease XPF-1 are abrogated, X DSBs persist, suggesting that single-strand annealing is engaged in the absence of HR. Our findings indicate that alteration in sister chromatid interactions and flexibility in DSB repair pathway choice accommodate hemizygosity on sex chromosomes.     

Protein kinase C-α negatively regulates EGF-induced PLC- activity through direct phosphorylation

Phospholipase C- is a PLC isozyme that contains a CDC25 homology domain and a pair of RA domains in addition to a conserved PLC catalytic domain. PLC- is activated by both growth factors and GPCR ligands in a distinct manner. Growth factors such as EGF stimulate PLC- in an RA2 domain-dependent manner through Ras and Rap. On the other hand, several GPCR ligands that are linked with Ga₁₂ or Ga₁₃ can activate PLC- by associating with GTP-RhoA. GTP-RhoA binds with the region in the PLC- Y domain. Gs-linked ligands such as PGE1 and adrenaline stimulate PLC- by cAMP-dependent activation of Epac and Rap2B. PLC- is important for cardiac development and function. In addition, several lines of evidence indicate that PLC- promotes cell growth in an activity-dependent or -independent manner. In particular, PLC--dependent suppression of EGF receptor downregulation contributes to its growth promoting activity. Proper regulation of PLC- activity is essential for preventing tumor formation. Our previous report indicated that EGF-dependent ubiquitination of PLC- is required for the control of PLC--dependent cell growth. Recently, we found that PLC- is phosphorylated by growth factor stimulation, and this is another mechanism of the negative regulation. PLC- is phosphorylated by PKC-α upon stimulation with growth factors such as EGF and PDGF. The EGF-induced phosphorylation of PLC- was abolished by PKC inhibitors and by the expression of the dominant negative mutant of PKC-α. Furthermore, PKC-α was found to phosphorylate PLC- directly in vitro, suggesting that PLC- is a substrate of PKC-α in cells. In addition, PLC- was co-immunoprecipitated with PKC-α in an EGF-dependent manner. Immunocytochemical studies showed that PLC- co-localized with PKC-α in the plasma membrane after EGF stimulation. In addition, inhibition of PKC activity enhanced PLC--mediated PIP₂ hydrolysis, suggesting that PKC-α negatively regulates PLC- activity. Taken together, these results suggest for the first time that PLC- is regulated by PKC-α-dependent phosphorylation.     

Transformer‐4 version 2.0.1, a free multi‐platform software to quickly reformat genotype matrices of any marker type,and archive them in the Demiurge information system

Transformer‐4 version 2.0.1 (T4) is a multi‐platform freeware programmed in java that can transform a genotype matrix in Excel or XML format into the input formats of one or several of the most commonly used population genetic software, for any possible combination of the populations that the matrix contains. T4 also allows the users to (i) draw allozyme gel interpretations for any number of diploid individuals, and then generate a genotype matrix ready to be used by T4; and (ii) produce basic reports about the data in the matrices. Furthermore, T4 is the only way to optionally submit ‘genetic diversity digests’ for publication in the Demiurge online information system ( http://www.demiurge-project.org ). Each such digest undergoes peer‐review, and it consists of a geo‐referenced data matrix in the tfm4 format plus any ancillary document or hyperlink that the digest authors see fit to include. The complementarity between T4 and Demiurge facilitates a free, safe, permanent, and standardized data archival and analysis system for researchers, and may also be a convenient resource for scientific journals, public administrations, or higher educators. T4 and its converters are freely available (at, respectively, http://www.demiurge-project.org/download_t4 and http://www.demiurge-project.org/converterstore ) upon registration in the Demiurge information system ( http://demiurge-project.org/register ). Users have to click on the link provided on an account validation email, and accept Demiurge's terms of use (see http://www.demiurge-project.org/termsofuse ). A thorough user's guide is available within T4. A 3‐min promotional video about T4 and Demiurge can be seen at http://vimeo.com/29828406 .     

Proteogenomics in the context of the Human Proteome Project (HPP)

Introduction: The technological and scientific progress performed in the Human Proteome Project (HPP) has provided to the scientific community a new set of experimental and bioinformatic methods in the challenging field of shotgun and SRM/MRM-based Proteomics. The requirements for a protein to be considered experimentally validated are now well-established, and the information about the human proteome is available in the neXtProt database, while targeted proteomic assays are stored in SRMAtlas. However, the study of the missing proteins continues being an outstanding issue.

Areas covered: This review is focused on the implementation of proteogenomic methods designed to improve the detection and validation of the missing proteins. The evolution of the methodological strategies based on the combination of different omic technologies and the use of huge publicly available datasets is shown taking the Chromosome 16 Consortium as reference.

Expert commentary: Proteogenomics and other strategies of data analysis implemented within the C-HPP initiative could be used as guidance to complete in a near future the catalog of the human proteins. Besides, in the next years, we will probably witness their use in the B/D-HPP initiative to go a step forward on the implications of the proteins in the human biology and disease.     

FCDD: A Database for Fruit Crops Diseases

Fruit Crops Diseases Database (FCDD) requires a number of biotechnology and bioinformatics tools. The FCDD is a unique bioinformatics resource that compiles information about 162 details on fruit crops diseases, diseases type, its causal organism, images, symptoms and their control. The FCDD contains 171 phytochemicals from 25 fruits, their 2D images and their 20 possible sequences. This information has been manually extracted and manually verified from numerous sources, including other electronic databases, textbooks and scientific journals. FCDD is fully searchable and supports extensive text search. The main focus of the FCDD is on providing possible information of fruit crops diseases, which will help in discovery of potential drugs from one of the common bioresource-fruits. The database was developed using MySQL. The database interface is developed in PHP, HTML and JAVA. FCDD is freely available.

Availability

http://www.fruitcropsdd.com/     

Micronuclei as a marker for medical screening of subjects continuously occupationally exposed to low doses of ionizing radiation

Context: Genotoxicity assays are widely employed in human biomonitoring studies to assess genetic damage inflicted by genotoxic agents.

Objective: Evaluation of micronuclei (MN) as a screening marker of occupational ionizing radiation (IR) exposure.

Materials and methods: Using micronucleus test, peripheral blood lymphocytes (PBL) of 402 control and exposed subjects were screened for genetic damage.

Results: The mean frequencies of micronucleus test parameters were significantly higher in exposed persons. Increase of micronucleus yield with duration of exposure (DOE) by 0.303MN/year was revealed.

Discussion and conclusion: The obtained data encourage us to consider MN as valuable markers for preventive medical screening of occupationally exposed groups.     

ModuleRole: A Tool for Modulization,Role Determination and Visualization in Protein-Protein Interaction Networks

Rapidly increasing amounts of (physical and genetic) protein-protein interaction (PPI) data are produced by various high-throughput techniques, and interpretation of these data remains a major challenge. In order to gain insight into the organization and structure of the resultant large complex networks formed by interacting molecules, using simulated annealing, a method based on the node connectivity, we developed ModuleRole, a user-friendly web server tool which finds modules in PPI network and defines the roles for every node, and produces files for visualization in Cytoscape and Pajek. For given proteins, it analyzes the PPI network from BioGRID database, finds and visualizes the modules these proteins form, and then defines the role every node plays in this network, based on two topological parameters Participation Coefficient and Z-score. This is the first program which provides interactive and very friendly interface for biologists to find and visualize modules and roles of proteins in PPI network. It can be tested online at the website http://www.bioinfo.org/modulerole/index.php, which is free and open to all users and there is no login requirement, with demo data provided by “User Guide” in the menu Help. Non-server application of this program is considered for high-throughput data with more than 200 nodes or user’s own interaction datasets. Users are able to bookmark the web link to the result page and access at a later time. As an interactive and highly customizable application, ModuleRole requires no expert knowledge in graph theory on the user side and can be used in both Linux and Windows system, thus a very useful tool for biologist to analyze and visualize PPI networks from databases such as BioGRID.

Availability

ModuleRole is implemented in Java and C, and is freely available at http://www.bioinfo.org/modulerole/index.php. Supplementary information (user guide, demo data) is also available at this website. API for ModuleRole used for this program can be obtained upon request.