A conserved domain of Drosophila RNA-binding protein Pumilio interacts with multiple CCR4-NOT deadenylase complex subunits to repress target mRNAs.

Pumilio is a sequence-specific RNA-binding protein that controls development, stem cell fate, and neurological functions in Drosophila. Pumilio represses protein expression by destabilizing target mRNAs in a manner dependent on the CCR4-NOT deadenylase complex. Three unique repression domains in the N-terminal region of Pumilio were previously shown to recruit CCR4-NOT, but how they do so was not well understood. In this study, we identified the motifs that are necessary and sufficient for the activity of the third repression domain of Pumilio, designated RD3, which is present in all isoforms and has conserved regulatory function. We identified multiple conserved regions of RD3 that are important for repression activity in cell-based reporter gene assays. Using yeast two-hybrid assays, we show that RD3 contacts specific regions of the Not1, Not2, and Not3 subunits of the CCR4-NOT complex. Our results indicate that RD3 makes multivalent interactions with CCR4-NOT mediated by conserved short linear interaction motifs. Specifically, two phenylalanine residues in RD3 make crucial contacts with Not1 that are essential for its repression activity. Using reporter gene assays, we also identify three new target mRNAs that are repressed by Pumilio and show that RD3 contributes to their regulation. Together, these results provide important insights into the mechanism by which Pumilio recruits CCR4-NOT to regulate the expression of target mRNAs.

Evolution of respiratory syncytial virus genotype BA in Kilifi, Kenya, 15 years on.

Respiratory syncytial virus (RSV) is recognised as a leading cause of severe acute respiratory disease and deaths among infants and vulnerable adults. Clinical RSV isolates can be divided into several known genotypes. RSV genotype BA, characterised by a 60-nucleotide duplication in the G glycoprotein gene, emerged in 1999 and quickly disseminated globally replacing other RSV group B genotypes. Continual molecular epidemiology is critical to understand the evolutionary processes maintaining the success of the BA viruses. We analysed 735 G gene sequences from samples collected from paediatric patients in Kilifi, Kenya, between 2003 and 2017. The virus population comprised of several genetically distinct variants (n = 56) co-circulating within and between epidemics. In addition, there was consistent seasonal fluctuations in relative genetic diversity. Amino acid changes increasingly accumulated over the surveillance period including two residues (N178S and Q180R) that mapped to monoclonal antibody 2D10 epitopes, as well as addition of putative N-glycosylation sequons. Further, switching and toggling of amino acids within and between epidemics was observed. On a global phylogeny, the BA viruses from different countries form geographically isolated clusters suggesting substantial localized variants. This study offers insights into longitudinal population dynamics of a globally endemic RSV genotype within a discrete location.

CRISPR-based assessment of genomic structure in the conserved SQUAMOSA promoter-binding-like gene clusters in rice.

Although SQUAMOSA promoter-binding-like (SPL) transcription factors are important regulators of development in rice (Oryza sativa), prior assessments of the SPL family have been limited to single genes. A functional comparison across the full gene family in standardized genetic backgrounds has not been reported previously. Here, we demonstrate that the SPL gene family in rice is enriched due to the most recent whole genome duplication (WGD). Notably, 10 of 19 rice SPL genes (52%) cluster in four units that have persisted for at least 50 million years. We show that SPL gene grouping and retention following WGD is widespread in angiosperms, suggesting the conservatism and importance of this gene arrangement. We used Cas9 editing to generate transformation lines for all 19 SPL genes in a common set of backgrounds, and found that knockouts of 14 SPL genes exhibited defects in plant height, 10 exhibited defects in panicle size, and nine had altered grain lengths. We observed subfunctionalization of genes in the paleoduplicated pairs, but little evidence of neofunctionalization. Expression of OsSPL3 was negatively correlated with that of its closest neighbor in its synteny group, OsSPL4, and its sister paired gene, OsSPL12, in the opposing group. Nucleotide diversity was lower in eight of the nine singleton genes in domesticated rice, relative to wild rice, whereas the reverse was true for the paired genes. Together, these results provide functional information on eight previously unexamined OsSPL family members and suggest that paleoduplicate pair redundancy benefits plant survival and innovation.

Comparative analysis of the genetic variability within the Q-type C2H2 zinc-finger transcription factors in the economically important cabbage, canola and Chinese cabbage genomes.

BACKGROUND: Brassica oleracea, B. rapa and B. napus encompass many economically important vegetable and oil crops; such as cabbage, broccoli, canola and Chinese cabbage. The genome sequencing of these species allows for gene discovery with an eye towards discerning the natural variability available for future breeding. The Q-type C2H2 zinc-finger protein (ZFP) transcription factors contain zinc finger motifs with a conserved QALGGH as part of the motif and they may play a critical role in the plants response to stress. While they may contain from one to five ZF domains (ZFD) this work focuses on the ZFPs that contain two zinc-fingers, which bind to the promoter of genes, and negatively regulate transcription via the EAR motif. B. oleracea and rapa are diploid and evolved into distinct species about 3.7 million years ago. B. napus is polyploid and formed by fusion of the diploids about 7500 years ago. RESULTS: This work identifies a total of 146 Q-type C2H2-ZFPs with 37 in B. oleracea, 35 in B. rapa and 74 in B. napus. The level of sequence similarity and arrangement of these genes on their chromosomes have mostly remained intact in B. napus, when compared to the chromosomes inherited from either B. rapa or oleracea. In contrast, the difference between the protein sequences of the orthologs of B. rapa and oleracea is greater and their organization on the chromosomes is much more divergent. In general, the 146 proteins are highly conserved especially within the known motifs. Differences within subgroups of ZFPs were identified. Considering that B. napus has twice the number of these proteins in its genome, RNA-Seq data was mined and the expression of 68 of the 74 genes was confirmed. CONCLUSION: Alignment of these proteins gives a snapshot of the variability that may be available naturally in Brassica species. The aim is to study how different ZFPs bind different genes or how dissimilar EAR motifs alter the negative regulation of the genes bound to the ZFP. Results from such studies could be used to enhance tolerance in future Brassica breeding programs.

Risks and Recommendations in Prenatally Detected De Novo Balanced Chromosomal Rearrangements from Assessment of Long-Term Outcomes.

The 6%-9% risk of an untoward outcome previously established by Warburton for prenatally detected de novo balanced chromosomal rearrangements (BCRs) does not account for long-term morbidity. We performed long-term follow-up (mean 17 years) of a registry-based nationwide cohort of 41 individuals carrying a prenatally detected de novo BCR with normal first trimester screening/ultrasound scan. We observed a significantly higher frequency of neurodevelopmental and/or neuropsychiatric disorders than in a matched control group (19.5% versus 8.3%, p = 0.04), which was increased to 26.8% upon clinical follow-up. Chromosomal microarray of 32 carriers revealed no pathogenic imbalances, illustrating a low prognostic value when fetal ultrasound scan is normal. In contrast, mate-pair sequencing revealed disrupted genes (ARID1B, NPAS3, CELF4), regulatory domains of known developmental genes (ZEB2, HOXC), and complex BCRs associated with adverse outcomes. Seven unmappable autosomal-autosomal BCRs with breakpoints involving pericentromeric/heterochromatic regions may represent a low-risk group. We performed independent phenotype-aware and blinded interpretation, which accurately predicted benign outcomes (specificity = 100%) but demonstrated relatively low sensitivity for prediction of the clinical outcome in affected carriers (sensitivity = 45%-55%). This sensitivity emphasizes the challenges associated with prenatal risk prediction for long-term morbidity in the absence of phenotypic data given the still immature annotation of the morbidity genome and poorly understood long-range regulatory mechanisms. In conclusion, we upwardly revise the previous estimates of Warburton to a morbidity risk of 27% and recommend sequencing of the chromosomal breakpoints as the first-tier diagnostic test in pregnancies with a de novo BCR.

Expression and function of microRNA-9 in the mid-hindbrain area of embryonic chick.

BACKGROUND: MiR-9 is a small non-coding RNA that is highly conserved between species and primarily expressed in the central nervous system (CNS). It is known to influence proliferation and neuronal differentiation in the brain and spinal cord of different vertebrates. Different studies have pointed to regional and species-specific differences in the response of neural progenitors to miR-9. METHODS: In ovo and ex ovo electroporation was used to overexpress or reduce miR-9 followed by mRNA in situ hybridisation and immunofluorescent stainings to evaluate miR- expression and the effect of changed miR-9 expression. RESULTS: We have investigated the expression and function of miR-9 during early development of the mid-hindbrain region (MH) in chick. Our analysis reveals a closer relationship of chick miR-9 to mammalian miR-9 than to fish and a dynamic expression pattern in the chick neural tube. Early in development, miR-9 is diffusely expressed in the entire brain, bar the forebrain, and it becomes more restricted to specific areas of the CNS at later stages. MiR-9 overexpression at HH9-10 results in a reduction of FGF8 expression and premature neuronal differentiation in the mid-hindbrain boundary (MHB). Within the midbrain miR-9 does not cause premature neuronal differentiation it rather reduces proliferation in the midbrain. CONCLUSION: Our findings indicate that miR-9 has regional specific effects in the developing mid-hindbrain region with a divergence of response of regional progenitors.

A conserved motif within cox 2 allows broad detection of economically important fruit flies (Diptera: Tephritidae).

The genera Anastrepha, Bactrocera, Ceratitis, Dacus and Rhagoletis in the family Tephritidae order Diptera are economically important, worldwide distributed and cause damage to a large number of commercially produced fruits and vegetables. China had regulated these five genera as quarantine pests, including the species Carpomya vesuviana. An accurate molecular method not depending on morphology able to detect all the quarantine fruit flies simultaneously is required for quarantine monitoring. This study contributes a comparative analysis of 146 mitochondrial genomes of Diptera species and found variable sites at the mt DNA cox2 gene only conserved in economically important fruit flies species. Degenerate primers (TephFdeg/TephR) were designed specific for the economically important fruit flies. A 603 bp fragment was amplified after testing each of the 40 selected representative species belonging to each economically important Tephritid genera, no diagnostic fragments were detected/amplified in any of the other Tephritidae and Diptera species examined. PCR sensitivity assays demonstrated the limit of detection of targeted DNA was 0.1 ng/µl. This work contributes an innovative approach for detecting all reported economically important fruit flies in a single-step PCR specific for reported fruit fly species of quarantine concern in China.

Network analysis of pseudogene-gene relationships: from pseudogene evolution to their functional potentials.

Pseudogenes are fossil relatives of genes. Pseudogenes have long been thought of as "junk DNAs", since they do not code proteins in normal tissues. Although most of the human pseudogenes do not have noticeable functions, ∼20% of them exhibit transcriptional activity. There has been evidence showing that some pseudogenes adopted functions as lncRNAs and work as regulators of gene expression. Furthermore, pseudogenes can even be "reactivated" in some conditions, such as cancer initiation. Some pseudogenes are transcribed in specific cancer types, and some are even translated into proteins as observed in several cancer cell lines. All the above have shown that pseudogenes could have functional roles or potentials in the genome. Evaluating the relationships between pseudogenes and their gene counterparts could help us reveal the evolutionary path of pseudogenes and associate pseudogenes with functional potentials. It also provides an insight into the regulatory networks involving pseudogenes with transcriptional and even translational activities.In this study, we develop a novel approach integrating graph analysis, sequence alignment and functional analysis to evaluate pseudogene-gene relationships, and apply it to human gene homologs and pseudogenes. We generated a comprehensive set of 445 pseudogene-gene (PGG) families from the original 3,281 gene families (13.56%). Of these 438 (98.4% PGG, 13.3% total) were non-trivial (containing more than one pseudogene). Each PGG family contains multiple genes and pseudogenes with high sequence similarity. For each family, we generate a sequence alignment network and phylogenetic trees recapitulating the evolutionary paths. We find evidence supporting the evolution history of olfactory family (both genes and pseudogenes) in human, which also supports the validity of our analysis method. Next, we evaluate these networks in respect to the gene ontology from which we identify functions enriched in these pseudogene-gene families and infer functional impact of pseudogenes involved in the networks. This demonstrates the application of our PGG network database in the study of pseudogene function in disease context.

Technological properties assessment and two component systems distribution of Streptococcus thermophilus strains isolated from fermented milk.

Streptococcus thermophilus is one of the economically most representatives of lactic acid bacteria, which is widely used as a starter to produce fermented milk products. In this study, 22 S. thermophilus strains were isolated from 26 fermented milk samples. Most isolates showed the ability to ferment a broad range of carbohydrates. Interestingly, eight strains are galactose positive, which is a desirable property in various industrial dairy fermentations. Four different nucleotide sequences were found in the galR-galK intergenic regions. The 16S-23S intergenic spacer region sequences of most isolates were determined as ITS-St-II type, which are related with protease positive and fast acidification. CS18 presented excellent technological performances, and showed potential as a promising starter candidate. To gain a comprehensive view of stress response mechanisms of strains, the distribution of all the two-component systems (TCSs) in strains were investigated. TCS analysis indicated that the nucleotide sequence of TCSs have obvious differences in different strains. And the strains with the special nucleotide sequences of TCS have distinctive traits. Therefore, it was speculated that there is a certain connection between the traits' difference and the TCS difference of strains.

Concerted millisecond timescale dynamics in the intrinsically disordered carboxyl terminus of γ-tubulin induced by mutation of a conserved tyrosine residue.

Tubulins are an ancient family of eukaryotic proteins characterized by an amino-terminal globular domain and disordered carboxyl terminus. These carboxyl termini play important roles in modulating the behavior of microtubules in living cells. However, the atomic-level basis of their function is not well understood. These regions contain multiple acidic residues and their overall charges are modulated in vivo by post-translational modifications, for example, phosphorylation. In this study, we describe an application of NMR and computer Monte Carlo simulations to investigate how the modification of local charge alters the conformational sampling of the γ-tubulin carboxyl terminus. We compared the dynamics of two 39-residue polypeptides corresponding to the carboxyl-terminus of yeast γ-tubulin. One polypeptide comprised the wild-type amino acid sequence while the second contained a Y > D mutation at Y11 in the polypeptide (Y445 in the full protein). This mutation introduces additional negative charge at a site that is phosphorylated in vivo and produces a phenotype with perturbed microtubule function. NMR relaxation measurements show that the Y11D mutation produces dramatic changes in the millisecond-timescale motions of the entire polypeptide. This observation is supported by Monte Carlo simulations that-similar to NMR-predict the WT γ-CT is largely unstructured and that the substitution of Tyr 11 with Asp causes the sampling of extended conformations that are unique to the Y11D polypeptide.

Characteristics of 29 novel atypical solute carriers of major facilitator superfamily type: evolutionary conservation, predicted structure and neuronal co-expression.

Solute carriers (SLCs) are vital as they are responsible for a major part of the molecular transport over lipid bilayers. At present, there are 430 identified SLCs, of which 28 are called atypical SLCs of major facilitator superfamily (MFS) type. These are MFSD1, 2A, 2B, 3, 4A, 4B, 5, 6, 6 L, 7, 8, 9, 10, 11, 12, 13A, 14A and 14B; SV2A, SV2B and SV2C; SVOP and SVOPL; SPNS1, SPNS2 and SPNS3; and UNC93A and UNC93B1. We studied their fundamental properties, and we also included CLN3, an atypical SLC not yet belonging to any protein family (Pfam) clan, because its involvement in the same neuronal degenerative disorders as MFSD8. With phylogenetic analyses and bioinformatic sequence comparisons, the proteins were divided into 15 families, denoted atypical MFS transporter families (AMTF1-15). Hidden Markov models were used to identify orthologues from human to Drosophila melanogaster and Caenorhabditis elegans Topology predictions revealed 12 transmembrane segments (for all except CLN3), corresponding to the common MFS structure. With single-cell RNA sequencing and in situ proximity ligation assay on brain cells, co-expressions of several atypical SLCs were identified. Finally, the transcription levels of all genes were analysed in the hypothalamic N25/2 cell line after complete amino acid starvation, showing altered expression levels for several atypical SLCs.

Genome-wide analysis of basic/helix-loop-helix gene family in peanut and assessment of its roles in pod development.

The basic/helix-loop-helix (bHLH) proteins constitute a superfamily of transcription factors that are known to play a range of regulatory roles in eukaryotes. Over the past few decades, many bHLH family genes have been well-characterized in model plants, such as Arabidopsis, rice and tomato. However, the bHLH protein family in peanuts has not yet been systematically identified and characterized. Here, 132 and 129 bHLH proteins were identified from two wild ancestral diploid subgenomes of cultivated tetraploid peanuts, Arachis duranensis (AA) and Arachis ipaensis (BB), respectively. Phylogenetic analysis indicated that these bHLHs could be classified into 19 subfamilies. Distribution mapping results showed that peanut bHLH genes were randomly and unevenly distributed within the 10 AA chromosomes and 10 BB chromosomes. In addition, 120 bHLH gene pairs between the AA-subgenome and BB-subgenome were found to be orthologous and 101 of these pairs were highly syntenic in AA and BB chromosomes. Furthermore, we confirmed that 184 bHLH genes expressed in different tissues, 22 of which exhibited tissue-specific expression. Meanwhile, we identified 61 bHLH genes that may be potentially involved in peanut-specific subterranean. Our comprehensive genomic analysis provides a foundation for future functional dissection and understanding of the regulatory mechanisms of bHLH transcription factors in peanuts.

A computational assessment of the predicted structures of Human Macrophage Migration Inhibitory Factor 1 orthologs in parasites and its affinity to human CD74 receptor.

The human macrophage migration inhibitory factor 1 (Hu-MIF-1) is a protein involved in the inflammatory and immunology response to parasite infection. In the present study, the existence of Hu-MIF-1 from parasites have been explored by mining WormBase. A total of 35 helminths were found to have Hu-MIF-1 homologs, including some parasites of importance for public health. Physicochemical, structural, and biological properties of Hu-MIF-1 were compared with its orthologs in parasites showing that most of these are secretory proteins, with positive net charge and presence of the Cys-Xaa-Xaa-Cys motif that is critical for its oxidoreductase activity. The inhibitor-binding site present in Hu-MIF-1 is well conserved among parasite MIFs suggesting that Hu-MIF inhibitors may target orthologs in pathogens. The binding of Hu-MIF-1 to its cognate receptor CD74 was predicted by computer-assisted docking, and it resulted to be very similar to the predicted complexes formed by parasite MIFs and human CD74. More than 1 plausible conformation of MIFs in the extracellular loops of CD74 may be possible as demonstrated by the different predicted conformations of MIF orthologs in complex with CD74. Parasite MIFs in complex with CD74 resulted with some charged residues oriented to CD74, which was not observed in the Hu-MIF-1/CD74 complex. Our findings predict the binding mode of Hu-MIF-1 and orthologs with CD74, which can assist in the design of novel MIF inhibitors. Whether the parasite MIFs function specifically subvert host immune responses to suit the parasite is an open question that needs to be further investigated. Future research should lead to a better understanding of parasite MIF action in the parasite biology.

The MiiA motif is a common marker present in polytopic surface proteins of oral and urinary tract invasive bacteria.

Many surface virulence factors of bacterial pathogens show mosaicism and confounding phylogenetic origin. The Streptococcus gordonii platelet-binding GspB protein, the Streptococcus sanguinis SrpA adhesin and the Streptococcus pneumoniae DiiA protein, share an imperfect 27-residue motif. Given the disparate domain architectures of these proteins and its association to invasive disease, this motif was named MiiA from Multiarchitecture invasion-involved motif A. MiiA is predicted to adopt a beta-sheet folding, probably related to the Ig-like fold, with a symmetrical positioning of two conserved aspartic residues. A specific hidden Markov model profiling MiiA was built, which specifically detected the motif in proteins from 58 species, mainly in cell-wall proteins from Gram-positive bacteria. These proteins contained one to ten MiiA motifs, which were embedded within larger repeat units of 70-82 residues. MiiA motifs combined to other domains and elements such as coiled-coils and low-complexity regions. The species carrying MiiA-proteins included commensals from the urogenital tract and the oral cavity, which can cause opportunistic endocarditis and sepsis. Intra-protein MiiA repeats showed a complex mixture of orthologal, paralogal and inter-species relationships, suggestive of a multistep origin. Presence of these repeats in proteins involved in oligosaccharide recognition and lifestyle of species suggest a putative function for MiiA repeats in sugars binding, probably those present in receptors of epithelial and blood cells. MiiA modules appear to have been transferred horizontally between species co-habiting in the same niche to create their own MiiA-containing determinants. The present work provides a global study and a catalog of potential MiiA virulence factors that should be analyzed experimentally.

The complex evolutionary history of aminoacyl-tRNA synthetases.

Aminoacyl-tRNA synthetases (AARSs) are a superfamily of enzymes responsible for the faithful translation of the genetic code and have lately become a prominent target for synthetic biologists. Our large-scale analysis of >2500 prokaryotic genomes reveals the complex evolutionary history of these enzymes and their paralogs, in which horizontal gene transfer played an important role. These results show that a widespread belief in the evolutionary stability of this superfamily is misconceived. Although AlaRS, GlyRS, LeuRS, IleRS, ValRS are the most stable members of the family, GluRS, LysRS and CysRS often have paralogs, whereas AsnRS, GlnRS, PylRS and SepRS are often absent from many genomes. In the course of this analysis, highly conserved protein motifs and domains within each of the AARS loci were identified and used to build a web-based computational tool for the genome-wide detection of AARS coding sequences. This is based on hidden Markov models (HMMs) and is available together with a cognate database that may be used for specific analyses. The bioinformatics tools that we have developed may also help to identify new antibiotic agents and targets using these essential enzymes. These tools also may help to identify organisms with alternative pathways that are involved in maintaining the fidelity of the genetic code.

Low-cost functional plasticity of TRPV1 supports heat tolerance in squirrels and camels.

The ability to sense heat is crucial for survival. Increased heat tolerance may prove beneficial by conferring the ability to inhabit otherwise prohibitive ecological niches. This phenomenon is widespread and is found in both large and small animals. For example, ground squirrels and camels can tolerate temperatures more than 40 °C better than many other mammalian species, yet a molecular mechanism subserving this ability is unclear. Transient receptor potential vanilloid 1 (TRPV1) is a polymodal ion channel involved in the detection of noxious thermal and chemical stimuli by primary afferents of the somatosensory system. Here, we show that thirteen-lined ground squirrels (Ictidomys tridecemlineatus) and Bactrian camels (Camelus ferus) express TRPV1 orthologs with dramatically reduced temperature sensitivity. The loss of sensitivity is restricted to temperature and does not affect capsaicin or acid responses, thereby maintaining a role for TRPV1 as a detector of noxious chemical cues. We show that heat sensitivity can be reengineered in both TRPV1 orthologs by a single amino acid substitution in the N-terminal ankyrin-repeat domain. Conversely, reciprocal mutations suppress heat sensitivity of rat TRPV1, supporting functional conservation of the residues. Our studies suggest that squirrels and camels co-opt a common molecular strategy to adapt to hot environments by suppressing the efficiency of TRPV1-mediated heat detection at the level of somatosensory neurons. Such adaptation is possible because of the remarkable functional flexibility of the TRPV1 molecule, which can undergo profound tuning at the minimal cost of a single amino acid change.

Conservation in Mammals of Genes Associated with Aggression-Related Behavioral Phenotypes in Honey Bees.

The emerging field of sociogenomics explores the relations between social behavior and genome structure and function. An important question is the extent to which associations between social behavior and gene expression are conserved among the Metazoa. Prior experimental work in an invertebrate model of social behavior, the honey bee, revealed distinct brain gene expression patterns in African and European honey bees, and within European honey bees with different behavioral phenotypes. The present work is a computational study of these previous findings in which we analyze, by orthology determination, the extent to which genes that are socially regulated in honey bees are conserved across the Metazoa. We found that the differentially expressed gene sets associated with alarm pheromone response, the difference between old and young bees, and the colony influence on soldier bees, are enriched in widely conserved genes, indicating that these differences have genomic bases shared with many other metazoans. By contrast, the sets of differentially expressed genes associated with the differences between African and European forager and guard bees are depleted in widely conserved genes, indicating that the genomic basis for this social behavior is relatively specific to honey bees. For the alarm pheromone response gene set, we found a particularly high degree of conservation with mammals, even though the alarm pheromone itself is bee-specific. Gene Ontology identification of human orthologs to the strongly conserved honey bee genes associated with the alarm pheromone response shows overrepresentation of protein metabolism, regulation of protein complex formation, and protein folding, perhaps associated with remodeling of critical neural circuits in response to alarm pheromone. We hypothesize that such remodeling may be an adaptation of social animals to process and respond appropriately to the complex patterns of conspecific communication essential for social organization.

In silico assessment of phosphorylation and O-ß-GlcNAcylation sites in human NPC1 protein critical for Ebola virus entry.

Ebola is a highly pathogenic enveloped virus responsible for deadly outbreaks of severe hemorrhagic fever. It enters human cells by binding a multifunctional cholesterol transporter Niemann-Pick C1 (NPC1) protein. Post translational modification (PTM) information for NPC1 is crucial to understand Ebola virus (EBOV) entry and action due to changes in phosphorylation or glycosylation at the binding site. It is difficult and costly to experimentally assess this type of interaction, so in silico strategy was employed. Identification of phosphorylation sites, including conserved residues that could be possible targets for 21 predicted kinases was followed by interplay study between phosphorylation and O-ß-GlcNAc modification of NPC1. Results revealed that only 4 out of 48 predicted phosphosites exhibited O-ß-GlcNAc activity. Predicted outcomes were integrated with residue conservation and 3D structural information. Three Yin Yang sites were located in the α-helix regions and were conserved in studied vertebrate and mammalian species. Only one modification site S425 was found in ß-turn region located near the N-terminus of NPC1 and was found to differ in pig, mouse, cobra and humans. The predictions suggest that Yin Yang sites may not be important for virus attachment to NPC1, whereas phosphosite 473 may be important for binding and hence entry of Ebola virus. This information could be useful in addressing further experimental studies and therapeutic strategies targeting PTM events in EBOV entry.

In silico identification of AMPylating enzymes and study of their divergent evolution.

AMPylation is a novel post-translational modification (PTM) involving covalent attachment of an AMP moiety to threonine/tyrosine side chains of a protein. AMPylating enzymes belonging to three different families, namely Fic/Doc, GS-ATase and DrrA have been experimentally characterized. Involvement of these novel enzymes in a myriad of biological processes makes them interesting candidates for genome-wide search. We have used SVM and HMM to develop a computational protocol for identification of AMPylation domains and their classification into various functional subfamilies catalyzing AMPylation, deAMPylation, phosphorylation and phosphocholine transfer. Our analysis has not only identified novel PTM catalyzing enzymes among unannotated proteins, but has also revealed how this novel enzyme family has evolved to generate functional diversity by subtle changes in sequence/structures of the proteins. Phylogenetic analysis of Fic/Doc has revealed three new isofunctional subfamilies, thus adding to their functional divergence. Also, frequent occurrence of Fic/Doc proteins on highly mobile and unstable genomic islands indicated their evolution via extensive horizontal gene transfers. On the other hand phylogenetic analyses indicate lateral evolution of GS-ATase family and an early duplication event responsible for AMPylation and deAMPylation activity of GS-ATase. Our analysis also reveals molecular basis of substrate specificity of DrrA proteins.

Genome-wide assessment of recurrent genomic imbalances in canine leukemia identifies evolutionarily conserved regions for subtype differentiation.

Leukemia in dogs is a heterogeneous disease with survival ranging from days to years, depending on the subtype. Strides have been made in both human and canine leukemia to improve classification and understanding of pathogenesis through immunophenotyping, yet classification and choosing appropriate therapy remains challenging. In this study, we assessed 123 cases of canine leukemia (28 ALLs, 24 AMLs, 25 B-CLLs, and 46 T-CLLs) using high-resolution oligonucleotide array comparative genomic hybridization (oaCGH) to detect DNA copy number alterations (CNAs). For the first time, such data were used to identify recurrent CNAs and inclusive genes that may be potential drivers of subtype-specific pathogenesis. We performed predictive modeling to identify CNAs that could reliably differentiate acute subtypes (ALL vs. AML) and chronic subtypes (B-CLL vs. T-CLL) and used this model to differentiate cases with up to 83.3 and 95.8 % precision, respectively, based on CNAs at only one to three genomic regions. In addition, CGH datasets for canine and human leukemia were compared to reveal evolutionarily conserved copy number changes between species, including the shared gain of HSA 21q in ALL and ∼25 Mb of shared gain of HSA 12 and loss of HSA 13q14 in CLL. These findings support the use of canine leukemia as a relevant in vivo model for human leukemia and justify the need to further explore the conserved genomic regions of interest for their clinical impact.