GWAS reveals genomic associations with swine inflammation and necrosis syndrome.

The recently identified swine inflammation and necrosis syndrome (SINS) occurs in high prevalence from newborn piglets to fattening pigs and resembles an important concern for animal welfare. The primary endogenous syndrome affects the tail, ears, teats, coronary bands, claws and heels. The basis of clinical inflammation and necrosis has been substantiated by histopathology, metabolomic and liver transcriptomic. Considerable variation in SINS scores is evident in offspring of different boars under the same husbandry conditions. The high complexity of metabolic alterations and the influence of the boar led to the hypothesis of a polygenic architecture of SINS. This should be investigated by a genome-wide association study. For this purpose, 27 sows were simultaneously inseminated with mixed semen from two extreme boars. The mixed semen always contained ejaculate from a Pietrain boar classified as extremely SINS susceptible and additionally either the ejaculate from a Pietrain boar classified as SINS stable or from a Duroc boar classified as SINS stable. The 234 piglets were phenotyped on day 3 of life, sampled and genetically assigned to the respective boar. The piglets showed the expected genetic differentiation with respect to SINS susceptibility. The suspected genetic complexity was confirmed both in the number and genome-wide distribution of 221 significantly associated SNPs, and led to 49 candidate genes. As the SNPs were almost exclusively located in noncoding regions, functional nucleotides have not yet been identified. The results suggest that the susceptibility of piglets to SINS depends not only on environmental conditions but also on genomic variation.

OVarFlow: a resource optimized GATK 4 based Open source Variant calling workFlow.

BACKGROUND: The advent of next generation sequencing has opened new avenues for basic and applied research. One application is the discovery of sequence variants causative of a phenotypic trait or a disease pathology. The computational task of detecting and annotating sequence differences of a target dataset between a reference genome is known as "variant calling". Typically, this task is computationally involved, often combining a complex chain of linked software tools. A major player in this field is the Genome Analysis Toolkit (GATK). The "GATK Best Practices" is a commonly referred recipe for variant calling. However, current computational recommendations on variant calling predominantly focus on human sequencing data and ignore ever-changing demands of high-throughput sequencing developments. Furthermore, frequent updates to such recommendations are counterintuitive to the goal of offering a standard workflow and hamper reproducibility over time. RESULTS: A workflow for automated detection of single nucleotide polymorphisms and insertion-deletions offers a wide range of applications in sequence annotation of model and non-model organisms. The introduced workflow builds on the GATK Best Practices, while enabling reproducibility over time and offering an open, generalized computational architecture. The workflow achieves parallelized data evaluation and maximizes performance of individual computational tasks. Optimized Java garbage collection and heap size settings for the GATK applications SortSam, MarkDuplicates, HaplotypeCaller, and GatherVcfs effectively cut the overall analysis time in half. CONCLUSIONS: The demand for variant calling, efficient computational processing, and standardized workflows is growing. The Open source Variant calling workFlow (OVarFlow) offers automation and reproducibility for a computationally optimized variant calling task. By reducing usage of computational resources, the workflow removes prior existing entry barriers to the variant calling field and enables standardized variant calling.

iCLIP analysis of RNA substrates of the archaeal exosome.

BACKGROUND: The archaeal exosome is an exoribonucleolytic multiprotein complex, which degrades single-stranded RNA in 3' to 5' direction phosphorolytically. In a reverse reaction, it can add A-rich tails to the 3'-end of RNA. The catalytic center of the exosome is in the aRrp41 subunit of its hexameric core. Its RNA-binding subunits aRrp4 and aDnaG confer poly(A) preference to the complex. The archaeal exosome was intensely characterized in vitro, but still little is known about its interaction with natural substrates in the cell, particularly because analysis of the transcriptome-wide interaction of an exoribonuclease with RNA is challenging. RESULTS: To determine binding sites of the exosome to RNA on a global scale, we performed individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) analysis with antibodies directed against aRrp4 and aRrp41 of the chrenarchaeon Sulfolobus solfataricus. A relatively high proportion (17-19%) of the obtained cDNA reads could not be mapped to the genome. Instead, they corresponded to adenine-rich RNA tails, which are post-transcriptionally synthesized by the exosome, and to circular RNAs (circRNAs). We identified novel circRNAs corresponding to 5' parts of two homologous, transposase-related mRNAs. To detect preferred substrates of the exosome, the iCLIP reads were compared to the transcript abundance using RNA-Seq data. Among the strongly enriched exosome substrates were RNAs antisense to tRNAs, overlapping 3'-UTRs and RNAs containing poly(A) stretches. The majority of the read counts and crosslink sites mapped in mRNAs. Furthermore, unexpected crosslink sites clustering at 5'-ends of RNAs was detected. CONCLUSIONS: In this study, RNA targets of an exoribonuclease were analyzed by iCLIP. The data documents the role of the archaeal exosome as an exoribonuclease and RNA-tailing enzyme interacting with all RNA classes, and underlines its role in mRNA turnover, which is important for adaptation of prokaryotic cells to changing environmental conditions. The clustering of crosslink sites near 5'-ends of genes suggests simultaneous binding of both RNA ends by the S. solfataricus exosome. This may serve to prevent translation of mRNAs dedicated to degradation in 3'-5' direction.

Comparative analyses of the variation of the transcriptome and proteome of Rhodobacter sphaeroides throughout growth.

BACKGROUND: In natural environments, bacteria must frequently cope with extremely scarce nutrients. Most studies focus on bacterial growth in nutrient replete conditions, while less is known about the stationary phase. Here, we are interested in global gene expression throughout all growth phases, including the adjustment to deep stationary phase. RESULTS: We monitored both the transcriptome and the proteome in cultures of the alphaproteobacterium Rhodobacter sphaeroides, beginning with the transition to stationary phase and at different points of the stationary phase and finally during exit from stationary phase (outgrowth) following dilution with fresh medium. Correlation between the transcriptomic and proteomic changes was very low throughout the growth phases. Surprisingly, even in deep stationary phase, the abundance of many proteins continued to adjust, while the transcriptome analysis revealed fewer adjustments. This pattern was reversed during the first 90 min of outgrowth, although this depended upon the duration of the stationary phase. We provide a detailed analysis of proteomic changes based on the clustering of orthologous groups (COGs), and compare these with the transcriptome. CONCLUSIONS: The low correlation between transcriptome and proteome supports the view that post-transcriptional processes play a major role in the adaptation to growth conditions. Our data revealed that many proteins with functions in transcription, energy production and conversion and the metabolism and transport of amino acids, carbohydrates, lipids, and secondary metabolites continually increased in deep stationary phase. Based on these findings, we conclude that the bacterium responds to sudden changes in environmental conditions by a radical and rapid reprogramming of the transcriptome in the first 90 min, while the proteome changes were modest. In response to gradually deteriorating conditions, however, the transcriptome remains mostly at a steady state while the bacterium continues to adjust its proteome. Even long after the population has entered stationary phase, cells are still actively adjusting their proteomes.

Cellular Gene Expression during Hepatitis C Virus Replication as Revealed by Ribosome Profiling.

BACKGROUND: Hepatitis C virus (HCV) infects human liver hepatocytes, often leading to liver cirrhosis and hepatocellular carcinoma (HCC). It is believed that chronic infection alters host gene expression and favors HCC development. In particular, HCV replication in Endoplasmic Reticulum (ER) derived membranes induces chronic ER stress. How HCV replication affects host mRNA translation and transcription at a genome wide level is not yet known. METHODS: We used Riboseq (Ribosome Profiling) to analyze transcriptome and translatome changes in the Huh-7.5 hepatocarcinoma cell line replicating HCV for 6 days. RESULTS: Established viral replication does not cause global changes in host gene expression-only around 30 genes are significantly differentially expressed. Upregulated genes are related to ER stress and HCV replication, and several regulated genes are known to be involved in HCC development. Some mRNAs (PPP1R15A/GADD34, DDIT3/CHOP, and TRIB3) may be subject to upstream open reading frame (uORF) mediated translation control. Transcriptional downregulation mainly affects mitochondrial respiratory chain complex core subunit genes. CONCLUSION: After establishing HCV replication, the lack of global changes in cellular gene expression indicates an adaptation to chronic infection, while the downregulation of mitochondrial respiratory chain genes indicates how a virus may further contribute to cancer cell-like metabolic reprogramming ("Warburg effect") even in the hepatocellular carcinoma cells used here.

Crystal structure of plasmoredoxin, a redox-active protein unique for malaria parasites.

Plasmoredoxin is a 22 â€‹kDa thiol-disulfide oxidoreductase involved in cellular redox regulatory processes and antioxidant defense. The 1.6 â€‹Å structure of the protein, solved via X-ray crystallography, adopts a modified thioredoxin fold. The structure reveals that plasmoredoxin, unique for malarial parasites, forms a new subgroup of thioredoxin-like proteins together with tryparedoxin, unique for kinetoplastids. Unlike most members of this superfamily, Plrx does not have a proline residue within the CxxC redox motif. In addition, the Plrx structure has a distinct C-terminal domain. Similar to human thioredoxin, plasmoredoxin forms monomers and dimers, which are also structurally similar to the human thioredoxin dimer, and, as in humans, plasmoredoxin is inactive as a dimer. Monomer-dimer equilibrium depends on the surrounding redox conditions, which could support the parasite in reacting to oxidative challenges. Based on structural considerations, the residues of the dimer interface are likely to interact with target proteins. In contrast to human and Plasmodium falciparum thioredoxin, however, there is a cluster of positively charged residues at the dimer interface of plasmoredoxin. These intersubunit (lysine) residues might allow binding of the protein to cellular membranes or to plasminogen. Malaria parasites lack catalase and glutathione peroxidase and therefore depend on their other glutathione and thioredoxin-dependent redox relays. Plasmoredoxin could be part of a so far unknown electron transfer system that only occurs in these parasites. Since the surface charge of plasmoredoxin differs significantly from other members of the thioredoxin superfamily, its three-dimensional structure can provide a model for designing selective redox-modulatory inhibitors.

Structural and Functional Characterization of Plasmodium falciparum Nicotinic Acid Mononucleotide Adenylyltransferase.

Nicotinic acid mononucleotide adenylyltransferase (NaMNAT) is an indispensable enzyme for the synthesis of NAD and NAD phosphate. It catalyzes the adenylylation of nicotinic acid mononucleotide (NaMN) to yield nicotinic acid adenine dinucleotide (NaAD). Since NAD(H) and NAD phosphate(H) are essentially involved in metabolic and redox regulatory reactions, NaMNAT is an attractive drug target in the fight against bacterial and parasitic infections. Notably, NaMNAT of the malaria parasite Plasmodium falciparum possesses only 20% sequence identity with the homologous human enzyme. Here, we present for the first time the two X-ray structures of P. falciparum NaMNAT (PfNaMNAT)-in the product-bound state with NaAD and complexed with an α,ß-non-hydrolizable ATP analog-the structures were determined to a resolution of 2.2Å and 2.5Å, respectively. The overall architecture of PfNaMNAT was found to be more similar to its bacterial homologs than its human counterparts although the PPHK motif conserved in bacteria is missing. Furthermore, PfNaMNAT possesses two cysteine residues within the active site that have not been described for any other NaMNATase so far and are likely to be involved in redox regulation of PfNaMNAT activity. Enzymatic studies and surface plasmon resonance data reveal that PfNaMNAT is capable of utilizing NaMN and nicotinamide mononucleotide with a slight preference for NaMN. Surprisingly, a comparison with the active site of Escherichia coli NaMNAT showed very similar architectures, despite different substrate preferences.

Targeting the intersubunit cavity of Plasmodium falciparum glutathione reductase by a novel natural inhibitor: computational and experimental evidence.

Glutathione reductase (GR), a homodimeric FAD-dependent disulfide reductase, is essential for redox homeostasis of the malaria parasite Plasmodium falciparum and has been proposed as an antimalarial drug target. In this study we performed a virtual screening against PfGR, using the structures of about 170,000 natural compounds. Analysis of the two top-scoring molecules, TTB and EPB, indicated that these ligands are likely to interact with the homodimer intersubunit cavity of PfGR with high binding energy scores of -9.67 and -9.60kcal/mol, respectively. Both compounds had a lower affinity for human GR due to differences in structure and electrostatic properties. In order to assess the putative interactions in motion, molecular dynamics simulations were carried out for 30ns, resulting in TTB being more dynamically and structurally favored than EPB. A closely related compound MDPI 21618 was tested on recombinant PfGR and hGR, resulting in IC50 values of 11.3±2.5µM and 10.2±1.7µM, respectively. Kinetic characterization of MDPI 21618 on PfGR revealed a mixed-type inhibition with respect to glutathione disulfide (Ki=9.7±2.3µM) and an uncompetitive inhibition with respect to NADPH. Furthermore, MDPI 21618 was found to inhibit the growth of the chloroquine-sensitive P. falciparum strain 3D7 with an IC50 of 3.2±1.9µM and the chloroquine-resistant Dd2 strain with an IC50 of 3.2+1.6µM. In drug combination assays with chloroquine, artemisinin, or mefloquine MDPI 21618 showed an antagonistic action, which might suggest partially overlapping routes of action. This study further substantiates research on PfGR as a potential antimalarial drug target.