Dmitry A. Afonnikov

Comparison of the three-dimensional organization of sperm and fibroblast genomes using the Hi-C approach

BackgroundThe three-dimensional organization of the genome is tightly connected to its biological function. The Hi-C approach was recently introduced as a method that can be used to identify higher-order chromatin interactions genome-wide. The aim of this study was to determine genome-wide chromatin interaction frequencies using the Hi-C approach in mouse sperm cells and embryonic fibroblasts.ResultsThe obtained data demonstrate that the three-dimensional genome organizations of sperm and fibroblast cells show a high degree of similarity both with each other and with the previously described mouse embryonic stem cells. Both A- and B-compartments and topologically associated domains are present in spermatozoa and fibroblasts. Nevertheless, sperm cells and fibroblasts exhibit statistically significant differences between each other in the contact probabilities of defined loci. Tight packaging of the sperm genome results in an enrichment of long-range contacts compared with the fibroblasts. However, only 30% of the differences in the number of contacts are based on differences in the densities of their genome packages; the main source of the differences is the gain or loss of contacts that are specific for defined genome regions. We find that the dependence of the contact probability on genomic distance for sperm is close to the dependence predicted for the fractal globular folding of chromatin.ConclusionsOverall, we can conclude that the three-dimensional structure of the genome is passed through generations without being dramatically changed in sperm cells.

Identification of microRNA Genes in Three Opisthorchiids

Background Opisthorchis felineus, O. viverrini, and Clonorchis sinensis (family Opisthorchiidae) are parasitic flatworms that pose a serious threat to humans in some countries and cause opisthorchiasis/clonorchiasis. Chronic disease may lead to a risk of carcinogenesis in the biliary ducts. MicroRNAs (miRNAs) are small noncoding RNAs that control gene expression at post-transcriptional level and are implicated in the regulation of various cellular processes during the parasite- host interplay. However, to date, the miRNAs of opisthorchiid flukes, in particular those essential for maintaining their complex biology and parasitic mode of existence, have not been satisfactorily described. Methodology/Principal Findings Using a SOLiD deep sequencing-bioinformatic approach, we identified 43 novel and 18 conserved miRNAs for O. felineus (miracidia, metacercariae and adult worms), 20 novel and 16 conserved miRNAs for O. viverrini (adult worms), and 33 novel and 18 conserved miRNAs for C. sinensis (adult worms). The analysis of the data revealed differences in the expression level of conserved miRNAs among the three species and among three the developmental stages of O. felineus. Analysis of miRNA genes revealed two gene clusters, one cluster-like region and one intronic miRNA in the genome. The presence and structure of the two gene clusters were validated using a PCR-based approach in the three flukes. Conclusions This study represents a comprehensive description of miRNAs in three members of the family Opistorchiidae, significantly expands our knowledge of miRNAs in multicellular parasites and provides a basis for understanding the structural and functional evolution of miRNAs in these metazoan parasites. Results of this study also provides novel resources for deeper understanding the complex parasite biology, for further research on the pathogenesis and molecular events of disease induced by the liver flukes. The present data may also facilitate the development of novel approaches for the prevention and treatment of opisthorchiasis/clonorchiasis.

Orthoscape: a cytoscape application for grouping and visualization KEGG based gene networks by taxonomy and homology principles

BackgroundThere are many available software tools for visualization and analysis of biological networks. Among them, Cytoscape (http://cytoscape.org/) is one of the most comprehensive packages, with many plugins and applications which extends its functionality by providing analysis of protein-protein interaction, gene regulatory and gene co-expression networks, metabolic, signaling, neural as well as ecological-type networks including food webs, communities networks etc. Nevertheless, only three plugins tagged ‘network evolution’ found in Cytoscape official app store and in literature. We have developed a new Cytoscape 3.0 application Orthoscape aimed to facilitate evolutionary analysis of gene networks and visualize the results.ResultsOrthoscape aids in analysis of evolutionary information available for gene sets and networks by highlighting: (1) the orthology relationships between genes; (2) the evolutionary origin of gene network components; (3) the evolutionary pressure mode (diversifying or stabilizing, negative or positive selection) of orthologous groups in general and/or branch-oriented mode. The distinctive feature of Orthoscape is the ability to control all data analysis steps via user-friendly interface.ConclusionOrthoscape allows its users to analyze gene networks or separated gene sets in the context of evolution. At each step of data analysis, Orthoscape also provides for convenient visualization and data manipulation.

Extraction of quantitative characteristics describing wheat leaf pubescence with a novel image-processing technique

Leaf pubescence (hairiness) in wheat plays an important biological role in adaptation to the environment. However, this trait has always been methodologically difficult to phenotype. An important step forward has been taken with the use of computer technologies. Computer analysis of a photomicrograph of a transverse fold line of a leaf is proposed for quantitative evaluation of wheat leaf pubescence. The image-processing algorithm is implemented in the LHDetect2 software program accessible as a Web service at http://wheatdb.org/lhdetect2. The results demonstrate that the proposed method is rapid, adequately assesses leaf pubescence density and the length distribution of trichomes and the data obtained using this method are significantly correlated with the density of trichomes on the leaf surface. Thus, the proposed method is efficient for high-throughput analysis of leaf pubescence morphology in cereal genetic collections and mapping populations.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species.

BackgroundThe identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures.ResultsThe current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified.ConclusionsThe results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

Identification of nuclear genes controlling chlorophyll synthesis in barley by RNA-seq

BackgroundAlbinism in plants is characterized by lack of chlorophyll and results in photosynthesis impairment, abnormal plant development and premature death. These abnormalities are frequently encountered in interspecific crosses and tissue culture experiments. Analysis of albino mutant phenotypes with full or partial chlorophyll deficiency can shed light on genetic determinants and molecular mechanisms of albinism. Here we report analysis of RNA-seq transcription profiling of barley (Hordeum vulgare L.) near-isogenic lines, one of which is a carrier of mutant allele of the Alm gene for albino lemma and pericarp phenotype (line i:BwAlm).Results1221 genome fragments have statistically significant changes in expression levels between lines i:BwAlm and Bowman, with 148 fragments having increased expression levels in line i:BwAlm, and 1073 genome fragments, including 42 plastid operons, having decreased levels of expression in line i:BwAlm. We detected functional dissimilarity between genes with higher and lower levels of expression in i:BwAlm line. Genes with lower level of expression in the i:BwAlm line are mostly associated with photosynthesis and chlorophyll synthesis, while genes with higher expression level are functionally associated with vesicle transport. Differentially expressed genes are shown to be involved in several metabolic pathways; the largest fraction of such genes was observed for the Calvin-Benson-Bassham cycle. Finally, de novo assembly of transcriptome contains several transcripts, not annotated in current H. vulgare genome version.ConclusionsOur results provide the new information about genes which could be involved in formation of albino lemma and pericarp phenotype. They demonstrate the interplay between nuclear and chloroplast genomes in this physiological process.

High temperature and pressure influence the interdomain orientation of Nip7 proteins from P. abyssi and P. furiosus: MD simulations.

Interactions between protein domains and their position and movement relative to each other are essential for the stability and normal functioning of a protein molecule. Features of the movement of domains may define the mechanism of enzymatic reactions. Therefore, the description of this motion is an important task in the analysis of the structures and functions of multidomain proteins. In the current work, we investigated the influence of pressure and temperature on changes in the movement of the two domains of the protein Nip7, expressed by deep-water (Pyrococcus abyssi) and shallow-water (Pyrococcus furiosus) archaea. The results of the present study show that the interdomain interfaces of the Nip7 proteins of P. abyssi and P. furiosus are formed by stable hydrophobic interactions. It was shown that high pressure and high temperature significantly changed the orientation of domains in Nip7 proteins which perhaps was connected with functional features of these domains. It was found that increasing the pressure significantly changed the angle of rotation of these domains, to a greater extent in the shallow-water protein, while an increase in temperature slightly reduced the angle of rotation of these domains. Moreover, the results suggest that the type of motion of the domains under study is similar to shear motion.

Residue coevolution reveals functionally important intramolecular interactions in formamidopyrimidine-DNA glycosylase

In protein evolution, functionally important intramolecular interactions, such as polar bridges or hydrophobic interfaces, tend to be conserved. We have analyzed coevolution of physicochemical properties in pairs of amino acid residues in the formamidopyrimidine-DNA glycosylase (Fpg) protein family, identified three conserved polar bridges (Arg54-Glu131, Gln234-Arg244, and Tyr170-Ser208 in the E. coli protein) located in known functional regions of the protein, and analyzed their roles by site-directed mutagenesis. The structure and molecular dynamic modeling showed that the coevolving pairs do not form isolated bridges but rather participate in tight local clusters of hydrogen bonds. The Arg54-Glu131 bridge, connecting the N- and C-terminal domains, was important for DNA binding, as its abolishment or even ion pair reversal inactivated Fpg and greatly decreased the enzymes affinity for DNA. Mutations of the Gln234-Arg244 bridge, located at the base of the single Fpg β-hairpin zinc finger, did not affect the activity but sharply decreased the melting temperature of the protein, with the bridge reversal partially restoring the thermal stability. Finally, Tyr170 mutation to Phe decreased Fpg binding but did not fully inactivate the protein, whereas Ser208 replacement with Ala had no effect; molecular dynamics showed that in both wild-type and S208 A Fpg, Tyr170 quickly re-orients to form an alternative set of hydrogen bonds. Thus, the coevolution analysis approach, combined with biochemical and computational studies, provides a powerful tool for understanding intramolecular interactions important for the function of DNA repair enzymes.

Differential expression of NBS-LRR-encoding genes in the root transcriptomes of two Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis

BackgroundThe characterization of major resistance genes (R genes) in the potato remains an important task for molecular breeding. However, R genes are rapidly evolving and frequently occur in genomes as clusters with complex structures, and their precise mapping and identification are complicated and time consuming.ResultsComparative analysis of root transcriptomes of Solanum phureja genotypes with contrasting resistance to Globodera rostochiensis revealed a number of differentially expressed genes. However, compiling a list of candidate R genes for further segregation analysis was hampered by their scarce annotation. Nevertheless, combination of transcriptomic analysis with data on predicted potato NBS-LRR-encoding genes considerably improved the quality of the results and provided a reasonable number of candidate genes that provide S. phureja with strong resistance to the potato golden cyst nematode.ConclusionCombination of comparative analyses of tissue-specific transcriptomes in resistant and susceptible genotypes may be used as an approach for the rapid identification of candidate potato R genes for co-segregation analysis and may be used in parallel with more sophisticated studies based on genome resequencing.

Metabolic pathways and genes identified by RNA-seq analysis of barley near-isogenic lines differing by allelic state of the Black lemma and pericarp ( Blp ) gene

BackgroundSome plant species have ‘melanin-like’ black seed pigmentation. However, the chemical and genetic nature of this ‘melanin-like’ black pigment have not yet been fully explored due to its complex structure and ability to withstand almost all solvents. Nevertheless, identification of genetic networks participating in trait formation is key to understanding metabolic processes involved in the expression of ‘melanin-like’ black seed pigmentation. The aim of the current study was to identify differentially expressed genes (DEGs) in barley near-isogenic lines (NILs) differing by allelic state of the Blp (black lemma and pericarp) locus.ResultsRNA-seq analysis of six libraries (three replicates for each line) was performed. A total of 957 genome fragments had statistically significant changes in expression levels between lines BLP and BW, with 632 fragments having increased expression levels in line BLP and 325 genome fragments having decreased expression. Among identified DEGs, 191 genes were recognized as participating in known pathways. Among these were metabolic pathways including ‘suberin monomer biosynthesis’, ‘diterpene phytoalexins precursors biosynthesis’, ‘cutin biosynthesis’, ‘cuticular wax biosynthesis’, and ‘phenylpropanoid biosynthesis, initial reactions’. Differential expression was confirmed by real-time PCR analysis of selected genes.ConclusionsMetabolic pathways and genes presumably associated with black lemma and pericarp colour as well as Blp-associated resistance to oxidative stress and pathogens, were revealed. We suggest that the black pigmentation of lemmas and pericarps is related to increased level of phenolic compounds and their oxidation. The effect of functional Blp on the synthesis of ferulic acid and other phenolic compounds can explain the increased antioxidant capacity and biotic and abiotic stress tolerance of black-grained cereals. Their drought tolerance and resistance to diseases affecting the spike may also be related to cuticular wax biosynthesis. In addition, upregulated synthesis of phytoalexins, suberin and universal stress protein (USP) in lemmas and pericarps of the Blp carriers may contribute to their increased disease resistance. Further description of the DEGs haplotypes and study of their association with physiological characteristics may be useful for future application in barley pre-breeding.