Mapping and identification of a potential candidate gene for a novel maturity locus, E10, in soybean.

KEY MESSAGE: E10 is a new maturity locus in soybean and FT4 is the predicted/potential functional gene underlying the locus. Flowering and maturity time traits play crucial roles in economic soybean production. Early maturity is critical for north and west expansion of soybean in Canada. To date, 11 genes/loci have been identified which control time to flowering and maturity; however, the molecular bases of almost half of them are not yet clear. We have identified a new maturity locus called "E10" located at the end of chromosome Gm08. The gene symbol E10e10 has been approved by the Soybean Genetics Committee. The e10e10 genotype results in 5-10 days earlier maturity than E10E10. A set of presumed E10E10 and e10e10 genotypes was used to identify contrasting SSR and SNP haplotypes. These haplotypes, and their association with maturity, were maintained through five backcross generations. A functional genomics approach using a predicted protein-protein interaction (PPI) approach (Protein-protein Interaction Prediction Engine, PIPE) was used to investigate approximately 75 genes located in the genomic region that SSR and SNP analyses identified as the location of the E10 locus. The PPI analysis identified FT4 as the most likely candidate gene underlying the E10 locus. Sequence analysis of the two FT4 alleles identified three SNPs, in the 5'UTR, 3'UTR and fourth exon in the coding region, which result in differential mRNA structures. Allele-specific markers were developed for this locus and are available for soybean breeders to efficiently develop earlier maturing cultivars using molecular marker assisted breeding.

Efficient prediction of human protein-protein interactions at a global scale.

BACKGROUND: Our knowledge of global protein-protein interaction (PPI) networks in complex organisms such as humans is hindered by technical limitations of current methods. RESULTS: On the basis of short co-occurring polypeptide regions, we developed a tool called MP-PIPE capable of predicting a global human PPI network within 3 months. With a recall of 23% at a precision of 82.1%, we predicted 172,132 putative PPIs. We demonstrate the usefulness of these predictions through a range of experiments. CONCLUSIONS: The speed and accuracy associated with MP-PIPE can make this a potential tool to study individual human PPI networks (from genomic sequences alone) for personalized medicine.

High-throughput design of bacterial anti-sense RNAs using CAREng.

Summary: Short RNA (sRNA) modulation of gene expression is an increasingly popular tool for bacterial functional genomics. Antisense pairing between an sRNA and a target messenger RNA results in post-transcriptional down-regulation of a specific gene and can thus be used both for investigating individual gene function and for large-scale genetic screens. sRNAs have several advantages over knockout libraries in studies of gene function, including inducibility, the capacity to interrogate essential genes and easy portability to multiple genetic backgrounds. High-throughput, systematic design of antisense RNAs will increase the efficiency and repeatability of sRNA screens. To this end, we present CAREng, the Computer-Automated sRNA Engineer. CAREng designs antisense RNAs for all coding sequences in a given genome, while checking for potential off-targets. Availability and implementation: CAREng is available as a Python script and through a web portal (https://caren.carleton.ca). Supplementary information: Supplementary data are available at Bioinformatics Advances online.

Binding site prediction for protein-protein interactions and novel motif discovery using re-occurring polypeptide sequences.

BACKGROUND: While there are many methods for predicting protein-protein interaction, very few can determine the specific site of interaction on each protein. Characterization of the specific sequence regions mediating interaction (binding sites) is crucial for an understanding of cellular pathways. Experimental methods often report false binding sites due to experimental limitations, while computational methods tend to require data which is not available at the proteome-scale. Here we present PIPE-Sites, a novel method of protein specific binding site prediction based on pairs of re-occurring polypeptide sequences, which have been previously shown to accurately predict protein-protein interactions. PIPE-Sites operates at high specificity and requires only the sequences of query proteins and a database of known binary interactions with no binding site data, making it applicable to binding site prediction at the proteome-scale. RESULTS: PIPE-Sites was evaluated using a dataset of 265 yeast and 423 human interacting proteins pairs with experimentally-determined binding sites. We found that PIPE-Sites predictions were closer to the confirmed binding site than those of two existing binding site prediction methods based on domain-domain interactions, when applied to the same dataset. Finally, we applied PIPE-Sites to two datasets of 2347 yeast and 14,438 human novel interacting protein pairs predicted to interact with high confidence. An analysis of the predicted interaction sites revealed a number of protein subsequences which are highly re-occurring in binding sites and which may represent novel binding motifs. CONCLUSIONS: PIPE-Sites is an accurate method for predicting protein binding sites and is applicable to the proteome-scale. Thus, PIPE-Sites could be useful for exhaustive analysis of protein binding patterns in whole proteomes as well as discovery of novel binding motifs. PIPE-Sites is available online at http://pipe-sites.cgmlab.org/.

Whole Genome Sequencing of the Mutamouse Model Reveals Strain- and Colony-Level Variation, and Genomic Features of the Transgene Integration Site.

The MutaMouse transgenic rodent model is widely used for assessing in vivo mutagenicity. Here, we report the characterization of MutaMouse's whole genome sequence and its genetic variants compared to the C57BL/6 reference genome. High coverage (>50X) next-generation sequencing (NGS) of whole genomes from multiple MutaMouse animals from the Health Canada (HC) colony showed ~5 million SNVs per genome, ~20% of which are putatively novel. Sequencing of two animals from a geographically separated colony at Covance indicated that, over the course of 23 years, each colony accumulated 47,847 (HC) and 17,677 (Covance) non-parental homozygous single nucleotide variants. We found no novel nonsense or missense mutations that impair the MutaMouse response to genotoxic agents. Pairing sequencing data with array comparative genomic hybridization (aCGH) improved the accuracy and resolution of copy number variants (CNVs) calls and identified 300 genomic regions with CNVs. We also used long-read sequence technology (PacBio) to show that the transgene integration site involved a large deletion event with multiple inversions and rearrangements near a retrotransposon. The MutaMouse genome gives important genetic context to studies using this model, offers insight on the mechanisms of structural variant formation, and contributes a framework to analyze aCGH results alongside NGS data.

Insights into the suitability of utilizing brown rats (Rattus norvegicus) as a model for healing spinal cord injury with epidermal growth factor and fibroblast growth factor-II by predicting protein-protein interactions.

The stimulation of the proliferation and differentiation of neural stem cells (NSCs) offers the possibility of a renewable source of replacement cells to treat numerous neurological diseases including spinal cord injury, traumatic brain injury and stroke. Epidermal growth factor (EGF) and fibroblast growth factor-2 (FGF-2) have been used to stimulate NSCs to renew, expand, and produce precursors for neural repair within an adult brown rat (Rattus norvegicus). To provide greater insight into the interspecies protein-protein interactions between human FGF-2 and EGF proteins and native R. norvegicus proteins, we have utilized the Massively Parallel Protein-Protein Interaction Prediction Engine (MP-PIPE) in an attempt to computationally shed light on the pathways potentially driving neurosphere proliferation. This study determined similar and differing protein interaction pathways between the two growth factors and the proteins in R. norvegicus compared with the proteins in H. sapiens. The protein-protein interactions predicted that EGF and FGF-2 may behave differently in rats than in humans. The identification and improved understanding of these differences may help to improve the clinical translation of NSC therapies from rats to humans.

In Silico Engineering of Synthetic Binding Proteins from Random Amino Acid Sequences.

Synthetic proteins with high affinity and selectivity for a protein target can be used as research tools, biomarkers, and pharmacological agents, but few methods exist to design such proteins de novo. To this end, the In-Silico Protein Synthesizer (InSiPS) was developed to design synthetic binding proteins (SBPs) that bind pre-determined targets while minimizing off-target interactions. InSiPS is a genetic algorithm that refines a pool of random sequences over hundreds of generations of mutation and selection to produce SBPs with pre-specified binding characteristics. As a proof of concept, we design SBPs against three yeast proteins and demonstrate binding and functional inhibition of two of three targets in vivo. Peptide SPOT arrays confirm binding sites, and a permutation array demonstrates target specificity. Our foundational approach will support the field of de novo design of small binding polypeptide motifs and has robust applicability while offering potential advantages over the limited number of techniques currently available.

Fitness Tradeoffs of Antibiotic Resistance in Extraintestinal Pathogenic Escherichia coli.

Evolutionary trade-offs occur when selection on one trait has detrimental effects on other traits. In pathogenic microbes, it has been hypothesized that antibiotic resistance trades off with fitness in the absence of antibiotic. Although studies of single resistance mutations support this hypothesis, it is unclear whether trade-offs are maintained over time, due to compensatory evolution and broader effects of genetic background. Here, we leverage natural variation in 39 extraintestinal clinical isolates of Escherichia coli to assess trade-offs between growth rates and resistance to fluoroquinolone and cephalosporin antibiotics. Whole-genome sequencing identifies a broad range of clinically relevant resistance determinants in these strains. We find evidence for a negative correlation between growth rate and antibiotic resistance, consistent with a persistent trade-off between resistance and growth. However, this relationship is sometimes weak and depends on the environment in which growth rates are measured. Using in vitro selection experiments, we find that compensatory evolution in one environment does not guarantee compensation in other environments. Thus, even in the face of compensatory evolution and other genetic background effects, resistance may be broadly costly, supporting the use of drug restriction protocols to limit the spread of resistance. Furthermore, our study demonstrates the power of using natural variation to study evolutionary trade-offs in microbes.

Uncharacterized ORF HUR1 influences the efficiency of non-homologous end-joining repair in Saccharomyces cerevisiae.

Non-Homologous End Joining (NHEJ) is a highly conserved pathway that repairs Double-Strand Breaks (DSBs) within DNA. Here we show that the deletion of yeast uncharacterized ORF HUR1, Hydroxyurea Resistance1 affects the efficiency of NHEJ. Our findings are supported by Protein-Protein Interaction (PPI), genetic interaction and drug sensitivity analyses. To assess the activity of HUR1 in DSB repair, we deleted its non-overlapping region with PMR1, referred to as HUR1-A. We observed that similar to deletion of TPK1 and NEJ1, and unlike YKU70 (important for NHEJ of DNA with overhang and not blunt end), deletion of HUR1-A reduced the efficiency of NHEJ in both overhang and blunt end plasmid repair assays. Similarly, a chromosomal repair assay showed a reduction for repair efficiency when HUR1-A was deleted. In agreement with a functional connection for Hur1p with Tpk1p and NEJ1p, double mutant strains Δhur1-A/Δtpk1, and Δhur1-A/Δnej1 showed the same reduction in the efficiency of plasmid repair, compared to both single deletion strains. Also, using a Homologous Recombination (HR) specific plasmid-based DSB repair assay we observed that deletion of HUR1-A influenced the efficiency of HR repair, suggesting that HUR1 might also play additional roles in other DNA repair pathways.

Evolution of protein-protein interaction networks in yeast.

Interest in the evolution of protein-protein and genetic interaction networks has been rising in recent years, but the lack of large-scale high quality comparative datasets has acted as a barrier. Here, we carried out a comparative analysis of computationally predicted protein-protein interaction (PPI) networks from five closely related yeast species. We used the Protein-protein Interaction Prediction Engine (PIPE), which uses a database of known interactions to make sequence-based PPI predictions, to generate high quality predicted interactomes. Simulated proteomes and corresponding PPI networks were used to provide null expectations for the extent and nature of PPI network evolution. We found strong evidence for conservation of PPIs, with lower than expected levels of change in PPIs for about a quarter of the proteome. Furthermore, we found that changes in predicted PPI networks are poorly predicted by sequence divergence. Our analyses identified a number of functional classes experiencing fewer PPI changes than expected, suggestive of purifying selection on PPIs. Our results demonstrate the added benefit of considering predicted PPI networks when studying the evolution of closely related organisms.

Designing anti-Zika virus peptides derived from predicted human-Zika virus protein-protein interactions.

The production of anti-Zika virus (ZIKV) therapeutics has become increasingly important as the propagation of the devastating virus continues largely unchecked. Notably, a causal relationship between ZIKV infection and neurodevelopmental abnormalities has been widely reported, yet a specific mechanism underlying impaired neurological development has not been identified. Here, we report on the design of several synthetic competitive inhibitory peptides against key pathogenic ZIKV proteins through the prediction of protein-protein interactions (PPIs). Often, PPIs between host and viral proteins are crucial for infection and pathogenesis, making them attractive targets for therapeutics. Using two complementary sequence-based PPI prediction tools, we first produced a comprehensive map of predicted human-ZIKV PPIs (involving 209 human protein candidates). We then designed several peptides intended to disrupt the corresponding host-pathogen interactions thereby acting as anti-ZIKV therapeutics. The data generated in this study constitute a foundational resource to aid in the multi-disciplinary effort to combat ZIKV infection, including the design of additional synthetic proteins.

Phosphatase complex Pph3/Psy2 is involved in regulation of efficient non-homologous end-joining pathway in the yeast Saccharomyces cerevisiae.

One of the main mechanisms for double stranded DNA break (DSB) repair is through the non-homologous end-joining (NHEJ) pathway. Using plasmid and chromosomal repair assays, we showed that deletion mutant strains for interacting proteins Pph3p and Psy2p had reduced efficiencies in NHEJ. We further observed that this activity of Pph3p and Psy2p appeared linked to cell cycle Rad53p and Chk1p checkpoint proteins. Pph3/Psy2 is a phosphatase complex, which regulates recovery from the Rad53p DNA damage checkpoint. Overexpression of Chk1p checkpoint protein in a parallel pathway to Rad53p compensated for the deletion of PPH3 or PSY2 in a chromosomal repair assay. Double mutant strains Δpph3/Δchk1 and Δpsy2/Δchk1 showed additional reductions in the efficiency of plasmid repair, compared to both single deletions which is in agreement with the activity of Pph3p and Psy2p in a parallel pathway to Chk1p. Genetic interaction analyses also supported a role for Pph3p and Psy2p in DNA damage repair, the NHEJ pathway, as well as cell cycle progression. Collectively, we report that the activity of Pph3p and Psy2p further connects NHEJ repair to cell cycle progression.

Short Co-occurring Polypeptide Regions Can Predict Global Protein Interaction Maps.

A goal of the post-genomics era has been to elucidate a detailed global map of protein-protein interactions (PPIs) within a cell. Here, we show that the presence of co-occurring short polypeptide sequences between interacting protein partners appears to be conserved across different organisms. We present an algorithm to automatically generate PPI prediction method parameters for various organisms and illustrate that global PPIs can be predicted from previously reported PPIs within the same or a different organism using protein primary sequences. The PPI prediction code is further accelerated through the use of parallel multi-core programming, which improves its usability for large scale or proteome-wide PPI prediction. We predict and analyze hundreds of novel human PPIs, experimentally confirm protein functions and importantly predict the first genome-wide PPI maps for S. pombe (∼9,000 PPIs) and C. elegans (∼37,500 PPIs).