Search

Fungal infections cause >1 million deaths annually and the emergence of antifungal resistance has prompted the exploration for novel antifungal targets. Quadruplexes are four-stranded nucleic acid secondary structures, which can regulate processes such as transcription, translation, replication and recombination. They are also found in genes linked to virulence in microbes, and ligands that bind to quadruplexes can eliminate drug-resistant pathogens. Using a computational approach, we quantified putative quadruplex-forming sequences (PQS) in 1359 genomes across the fungal kingdom and explored their presence in genes related to virulence, drug resistance and biological processes associated with pathogenicity in Aspergillus fumigatus. Here we present the largest analysis of PQS in fungi and identify significant heterogeneity of these sequences throughout phyla, genera and species. PQS were genetically conserved in Aspergillus spp. and frequently pathogenic species appeared to contain fewer PQS than their lesser/non-pathogenic counterparts. GO-term analysis identified that PQS-containing genes were involved in processes linked with virulence such as zinc ion binding, the biosynthesis of secondary metabolites and regulation of transcription in A. fumigatus. Although the genome frequency of PQS was lower in A. fumigatus, PQS could be found enriched in genes involved in virulence, and genes upregulated during germination and hypoxia. Moreover, PQS were found in genes involved in drug resistance. Quadruplexes could have important roles within fungal biology and virulence, but their roles require further elucidation.

Candida auris is an emerging antifungal resistant human fungal pathogen increasingly reported in healthcare facilities. It persists in hospital environments, and on skin surfaces, and can form biofilms readily. Here, we investigated the cell surface proteins from C. auris biofilms grown in a synthetic sweat medium mimicking human skin conditions. Cell surface proteins from both biofilm and planktonic control cells were extracted with a buffer containing β-mercaptoethanol and resolved by 2-D gel electrophoresis. Some of the differentially expressed proteins were excised and identified by mass spectrometry. C. albicans orthologs Spe3p, Tdh3p, Sod2p, Ywp1p, and Mdh1p were overexpressed in biofilm cells when compared to the planktonic cells of C. auris. Interestingly, several proteins with zinc ion binding activity were detected. Nrg1p is a zinc-binding transcription factor that negatively regulates hyphal growth in C. albicans. C. auris does not produce true hypha under standard in vitro growth conditions, and the role of Nrg1p in C. auris is currently unknown. Western blot analyses of cell surface and cytosolic proteins of C. auris against anti-CalNrg1 antibody revealed the Nrg1p in both locations. Cell surface localization of Nrg1p in C. auris, an unexpected finding, was further confirmed by immunofluorescence microscopy. Nrg1p expression is uniform across all four clades of C. auris and is dependent on growth conditions. Taken together, the data indicate that C. auris produces several unique proteins during its biofilm growth, which may assist in the skin-colonizing lifestyle of the fungus during its pathogenesis.

Dipeptidyl peptidase III (DPP III), a zinc- dependent exopeptidase, is widely distributed in organisms and present in almost all human tissues. In addition to its involvement in protein catabolism, it plays a role in oxidative stress and blood pressure regulation, and there is evidence of its involvement in pain modulation. Excess zinc ions have been found to inhibit its hydrolytic activity, but the binding affinity, binding site geometry, and mechanism of inhibitory activity have been unknown. Using several different computational approaches, we determined the inhibitory zinc ion binding site, its coordination and relative binding affinity. During some simulations the translocation of the zinc ion from the inhibitory to the catalytic binding site was observed, accompanied by movement of the catalytic zinc ion toward the exit of the substrate binding site. The traced behavior suggests an associative type of metal ion exchange, in which the formation of the ternary complex between enzyme and two metal ions precedes the exit of the catalytic metal ion. Differently from our previous findings that binding of a peptide induces partial opening of hDPP III, the globularity of the protein did not change in MD simulations of the hermorphin-like peptide bound to hDPP III with two zinc ions. However, the entrance to the interdomain cleft widens during Zn diffusion into the protein and was found to be the highest energy barrier in the process of metal translocation from the solvent to the active site. Finally, we discuss why excess zinc reduces enzyme activity.

Background Salinity expansion in arable land is a threat to crop plants. Rice is the staple food crop across several countries worldwide; however, its salt sensitive nature severely affects its growth under excessive salinity. FL478 is a salt tolerant indica recombinant inbred line, which can be a good source of salt tolerance at the seedling stage in rice. To learn about the genetic basis of its tolerance to salinity, we compared transcriptome profiles of FL478 and its sensitive parent (IR29) using RNA-seq technique. Results A total of 1714 and 2670 genes were found differentially expressed (DEGs) under salt stress compared to normal conditions in FL478 and IR29, respectively. Gene ontology analysis revealed the enrichment of transcripts involved in salinity response, regulation of gene expression, and transport in both genotypes. Comparative transcriptome analysis revealed that 1063 DEGs were co-expressed, while 338/252 and 572/908 DEGs were exclusively up/down-regulated in FL478 and IR29, respectively. Further, some biological processes (e.g. iron ion transport, response to abiotic stimulus, and oxidative stress) and molecular function terms (e.g. zinc ion binding and cation transmembrane transporter activity) were specifically enriched in FL478 up-regulated transcripts. Based on the metabolic pathways analysis, genes encoding transport and major intrinsic proteins transporter superfamily comprising aquaporin subfamilies and genes involved in MAPK signaling and signaling receptor kinases were specifically enriched in FL478. A total of 1135 and 1894 alternative splicing events were identified in transcripts of FL478 and IR29, respectively. Transcripts encoding two potassium transporters and two major facilitator family transporters were specifically up-regulated in FL478 under salt stress but not in the salt sensitive genotype. Remarkably, 11 DEGs were conversely regulated in the studied genotypes; for example, OsZIFL, OsNAAT, OsGDSL, and OsELIP genes were up-regulated in FL478, while they were down-regulated in IR29. Conclusions The achieved results suggest that FL478 employs more efficient mechanisms (especially in signal transduction of salt stress, influx and transport of k+, ionic and osmotic homeostasis, as well as ROS inhibition) to respond to the salt stress compared to its susceptible parent. Electronic supplementary material The online version of this article (10.1186/s12284-019-0273-2) contains supplementary material, which is available to authorized users.

Although specific interactions between host and pathogen genotypes have been well documented in invertebrates, the identification of host genes involved in discriminating pathogen genotypes remains a challenge. In the mosquito Aedes aegypti, the main dengue virus (DENV) vector worldwide, statistical associations between host genetic markers and DENV types or strains were previously detected, but the host genes underlying this genetic specificity have not been identified. In particular, it is unknown whether DENV type- or strain-specific resistance relies on allelic variants of the same genes or on distinct gene sets. Here, we investigated the genetic architecture of DENV resistance in a population of Ae. aegypti from Bakoumba, Gabon, which displays a stronger resistance phenotype to DENV type 1 (DENV-1) than to DENV type 3 (DENV-3) infection. Following experimental exposure to either DENV-1 or DENV-3, we sequenced the exomes of large phenotypic pools of mosquitoes that are either resistant or susceptible to each DENV type. Using variation in single-nucleotide polymorphism (SNP) frequencies among the pools, we computed empirical p values based on average gene scores adjusted for the differences in SNP counts, to identify genes associated with infection in a DENV type-specific manner. Among the top 5% most significant genes, 263 genes were significantly associated with resistance to both DENV-1 and DENV-3, 287 genes were only associated with DENV-1 resistance and 290 were only associated with DENV-3 resistance. The shared significant genes were enriched in genes with ATP binding activity and sulfur compound transmembrane transporter activity, whereas the genes uniquely associated with DENV-3 resistance were enriched in genes with zinc ion binding activity. Together, these results indicate that specific resistance to different DENV types relies on largely non-overlapping sets of genes in this Ae. aegypti population and pave the way for further mechanistic studies., Author summary Compatibility between hosts and pathogens is often genetically specific in invertebrates but host genes underlying this genetic specificity have not been elucidated. We investigated the genetic architecture of dengue virus type-specific resistance in the mosquito vector Aedes aegypti. We used a natural population of Ae. aegypti from Bakoumba, Gabon, which is differentially resistant to dengue virus type 1 and dengue virus type 3. We surveyed genetic variation in protein-coding regions of the mosquito genome and compared the frequency of genetic polymorphisms between groups of mosquitoes that are either resistant or susceptible to each dengue virus type. We found that the Ae. aegypti genes associated with resistance to dengue virus type 1 or dengue virus type 3 were largely non-overlapping. This finding indicates that different sets of host genes, rather than different variants of the same genes, confer pathogen-specific resistance in this population. This study is an important step towards identification of mechanisms underlying the genetic specificity of invertebrate host-pathogen interactions.

mRNA processing is highly regulated during development through changes in RNA-binding protein (RBP) activities. CUG-BP, Elav-like family member 1 (CELF1, also called CUGBP1) is an RBP, the expression of which decreases in skeletal muscle soon after birth. CELF1 regulates multiple nuclear and cytoplasmic RNA processing events. In the nucleus, CELF1 regulates networks of postnatal alternative splicing (AS) transitions, while in the cytoplasm, CELF1 regulates mRNA stability and translation. Stabilization and misregulation of CELF1 has been implicated in human diseases including myotonic dystrophy type 1, Alzheimer’s disease and multiple cancers. To understand the contribution of nuclear and cytoplasmic CELF1 activity to normal and pathogenic skeletal muscle biology, we generated transgenic mice for doxycycline-inducible and skeletal muscle-specific expression of active CELF1 mutants engineered to be localized predominantly to either the nucleus or the cytoplasm. Adult mice expressing nuclear, but not cytoplasmic, CELF1 are characterized by strong histopathological defects, muscle loss within 10 days and changes in AS. In contrast, mice expressing cytoplasmic CELF1 display changes in protein levels of targets known to be regulated at the level of translation by CELF1, with minimal changes in AS. These changes are in the absence of overt histopathological changes or muscle loss. RNA-sequencing revealed extensive gene expression and AS changes in mice overexpressing nuclear and naturally localized CELF1 protein, with affected genes involved in cytoskeleton dynamics, membrane dynamics, RNA processing and zinc ion binding. These results support a stronger role for nuclear CELF1 functions as compared to cytoplasmic CELF1 functions in skeletal muscle wasting.

BackgroundMicroRNAs (miRNAs) regulate adipose tissue development, which are closely related to subcutaneous and intramuscular fat deposition and adipocyte differentiation. As an important economic and agricultural animal, rabbits have low adipose tissue deposition and are an ideal model to study adipose regulation. However, the miRNAs related to fat deposition during the growth and development of rabbits are poorly defined.MethodsIn this study, miRNA-sequencing and bioinformatics analyses were used to profile the miRNAs in rabbit perirenal adipose tissue at 35, 85 and 120 days post-birth. Differentially expressed (DE) miRNAs between different stages were identified by DEseq in R. Target genes of DE miRNAs were predicted by TargetScan and miRanda. To explore the functions of identified miRNAs, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed.ResultsApproximately 1.6 GB of data was obtained by miRNA-seq. A total of 987 miRNAs (780 known and 207 newly predicted) and 174 DE miRNAs were identified. The miRNAs ranged from 18 nt to 26 nt. GO enrichment and KEGG pathway analyses revealed that the target genes of the DE miRNAs were mainly involved in zinc ion binding, regulation of cell growth, MAPK signaling pathway, and other adipose hypertrophy-related pathways. Six DE miRNAs were randomly selected, and their expression profiles were validated by q-PCR.ConclusionsThis is the first report of the miRNA profiles of adipose tissue during different growth stages of rabbits. Our data provide a theoretical reference for subsequent studies on rabbit genetics, breeding and the regulatory mechanisms of adipose development.

In the biological process of testicular spermatogenesis, the expression and interaction of many genes are regulated by microRNAs (miRNAs). However, comparisons of miRNA expression between descended testes (DTs) and undescended testes (UDTs) are rarely done in horses. In this study, we selected two UDTs (CKY2b and GU4b) from Chakouyi (CKY) and Guanzhong (GU) horses and eight DTs (GU1&ndash, 3, CKY1, CKY3, CKY2a, GU4a, and GU5). Three groups were compared to evaluate expression patterns of testicular miRNA in stallion testes. Group 1 compared normal CKY horses and GU horses (CKY1 and CKY3 vs. GU1&ndash, 3). Group 2 (CKY2a and GU4a (DTs) vs. CKY2b and GU4b (UDTs)) and group 3 (GU1&ndash, 3, CKY1, CKY3 (DTs) vs. CKY2b and GU4b (UDTs)) compared the expression levels in unilateral retained testes to normal testes. The results show that 42 miRNAs (7 upregulated and 35 downregulated) had significantly different expression levels in both comparisons. The expression levels of eca-miR-545, eca-miR-9084, eca-miR-449a, eca-miR-9024, eca-miR-9121, eca-miR-8908e, eca-miR-136, eca-miR-329b, eca-miR-370, and eca-miR-181b were further confirmed by quantitative real-time PCR assay. The target genes of differentially expressed miRNAs in three comparisons were predicted, and the functions were annotated. The putative target genes of the 42 co-differentially expressed miRNAs were annotated to 15 functional terms, including metal ion binding, GTPase activator activity, zinc ion binding, intracellular, cytoplasm, and cancer pathways, and osteoclast differentiation. Our data indicate that the differentially expressed miRNAs in undescended testis suggests a potential role in male fertility and a relationship with cryptorchidism in horses. The discovery of miRNAs in stallion testes might contribute to a new direction in the search for biomarkers of stallion fertility.

Background The transamidase complex is a molecular machine in the endoplasmic reticulum of eukaryotes that attaches a glycosylphosphatidylinositol (GPI) lipid anchor to substrate proteins after cleaving a C-terminal propeptide with a defined sequence signal. Its five subunits are very hydrophobic; thus, solubility, heterologous expression and complex reconstruction are difficult. Therefore, theoretical approaches are currently the main source of insight into details of 3D structure and of the catalytic process. Results In this work, we generated model 3D structures of the lumenal domain of human GPAA1, the M28-type metallo-peptide-synthetase subunit of the transamidase, including zinc ion and model substrate positions. In comparative molecular dynamics (MD) simulations of M28-type structures and our GPAA1 models, we estimated the metal ion binding energies with evolutionary conserved amino acid residues in the catalytic cleft. We find that canonical zinc binding sites 2 and 3 are strongest binders for Zn1 and, where a second zinc is available, sites 2 and 4 for Zn2. Zinc interaction of site 5 with Zn1 enhances upon substrate binding in structures with only one zinc. Whereas a previously studied glutaminyl cyclase structure, the best known homologue to GPAA1, binds only one zinc ion at the catalytic site, GPAA1 can sterically accommodate two. The M28-type metallopeptidases segregate into two independent branches with regard to one/two zinc ion binding modality in a phylogenetic tree where the GPAA1 family is closer to the joint origin of both groups. For GPAA1 models, MD studies revealed two large loops (flaps) surrounding the active site being involved in an anti-correlated, breathing-like dynamics. Conclusions In the light of combined sequence-analytic and phylogenetic arguments as well as 3D structural modelling results, GPAA1 is most likely a single zinc ion metallopeptidase. Two large flaps environ the catalytic site restricting access to large substrates. Reviewers This article was reviewed by Thomas Dandekar (MD) and Michael Gromiha.

The HIV-1 nucleocapsid 7 (NCp7) plays crucial roles in multiple stages of HIV-1 life cycle, and its biological functions rely on the binding of zinc ions. Understanding the molecular mechanism of how the zinc ions modulate the conformational dynamics and functions of the NCp7 is essential for the drug development and HIV-1 treatment. In this work, using a structure-based coarse-grained model, we studied the effects of zinc cofactors on the folding and target RNA(SL3) recognition of the NCp7 by molecular dynamics simulations. After reproducing some key properties of the zinc binding and folding of the NCp7 observed in previous experiments, our simulations revealed several interesting features in the metal ion modulated folding and target recognition. Firstly, we showed that the zinc binding makes the folding transition states of the two zinc fingers less structured, which is in line with the Hammond effect observed typically in mutation, temperature or denaturant induced perturbations to protein structure and stability. Secondly, We showed that there exists mutual interplay between the zinc ion binding and NCp7-target recognition. Binding of zinc ions enhances the affinity between the NCp7 and the target RNA, whereas the formation of the NCp7-RNA complex reshapes the intrinsic energy landscape of the NCp7 and increases the stability and zinc affinity of the two zinc fingers. Thirdly, by characterizing the effects of salt concentrations on the target RNA recognition, we showed that the NCp7 achieves optimal balance between the affinity and binding kinetics near the physiologically relevant salt concentrations. In addition, the effects of zinc binding on the inter-domain conformational flexibility and folding cooperativity of the NCp7 were also discussed.

The 2525 amino acid SMRT corepressor is an intrinsically disordered hub protein responsible for binding and coordinating the activities of multiple transcription factors and chromatin modifying enzymes. Here we have studied its interaction with HDAC7, a class IIa deacetylase that interacts with the corepressor complex together with the highly active class I deacetylase HDAC3. The binding site of class IIa deacetylases was previously mapped to an approximate 500 amino acid region of SMRT, with recent implication of short glycine-serine-isoleucine (GSI) containing motifs. In order to characterize the interaction in detail, we applied a random library screening approach within this region and obtained a range of stable, soluble SMRT fragments. In agreement with an absence of predicted structural domains, these were characterized as intrinsically disordered by NMR spectroscopy. We identified one of them, comprising residues 1255–1452, as interacting with HDAC7 with micromolar affinity. The binding site was mapped in detail by NMR and confirmed by truncation and alanine mutagenesis. Complementing this with mutational analysis of HDAC7, we show that HDAC7, via its surface zinc ion binding site, binds to a 28 residue stretch in SMRT comprising a GSI motif followed by an alpha helix.

Adipose tissue was the major energy deposition site of the mammals and provided the energy for the body and released the external pressure to the internal organs. In animal production, fat deposition in muscle can affect the meat quality, especially the intramuscular fat (IMF) content. Diacylglycerol acyltransferase-1 (DGAT1) was the key enzyme to control the synthesis of the triacylglycerol in adipose tissue. In order to better understand the regulation mechanism of the DGAT1 in the intramuscular fat deposition, the global gene expression profiling was performed in gastrocnemius muscle between DGAT1 transgenic mice and wild-type mice by microarray. 281 differentially expressed transcripts were identified with at least 1.5-fold change and thepvalue < 0.05. 169 transcripts were upregulated and 112 transcripts were downregulated. Ten genes (SREBF1, DUSP1, PLAGL1, FKBP5, ZBTB16, PPP1R3C, CDC14A, GLUL, PDK4, and UCP3) were selected to validate the reliability of the chip’s results by the real-time PCR. The finding of RT-PCR was consistent with the gene chip. Seventeen signal pathways were analyzed using KEGG pathway database and the pathways concentrated mainly on the G-protein coupled receptor protein signaling pathway, signal transduction, oxidation-reduction reaction, olfactory receptor activity, protein binding, and zinc ion binding. This study implied a function role of DGAT1 in the synthesis of TAG, insulin resistance, and IMF deposition.

Cactophilic Drosophila species provide a valuable model to study gene–environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii–D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii–D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches., This work was supported by grants BFU2008-04988 and BFU2011-30476 from Ministerio de Ciencia e Innovación (Spain) to A.R., by an FPI fellowship to Y.G. and a PIF-UAB fellowship to N.R, and by the National Institute of General Medical Sciences of the National Institute of Health under award number R01GM071813 to E.B.

Plants are exposed to antibiotics produced by soil microorganisms, but little is known about their responses at the transcriptional level. Likewise, few endogenous mechanisms of antibiotic resistance have been reported. The Arabidopsis thaliana ATP Binding Cassette (ABC) transporter AtWBC19 (ABCG19) is known to confer kanamycin resistance, but the exact mechanism of resistance is not well understood. Here we examined the transcriptomes of control seedlings and wbc19 mutant seedlings using RNA-seq analysis. Exposure to kanamycin indicated changes in the organization of the photosynthetic apparatus, metabolic fluxes and metal uptake. Elemental analysis showed a 60% and 80% reduction of iron uptake in control and wbc19 mutant seedlings respectively, upon exposure to kanamycin. The drop in iron content was accompanied by the upregulation of the gene encoding for FERRIC REDUCTION OXIDASE 6 (FRO6) in mutant seedlings but not by the differential expression of other transport genes known to be induced by iron deficiency. In addition, wbc19 mutants displayed a distinct expression profile in the absence of kanamycin. Most notably the expression of several zinc ion binding proteins, including ZINC TRANSPORTER 1 PRECURSOR (ZIP1) was increased, suggesting abnormal zinc uptake. Elemental analysis confirmed a 50% decrease of zinc content in wbc19 mutants. Thus, the antibiotic resistance gene WBC19 appears to also have a role in zinc uptake.

The biosynthesis of the active metal-bound form of the nickel-dependent enzyme urease involves the formation of a lysine-carbamate functional group concomitantly with the delivery of two Ni(2+) ions into the precast active site of the apoenzyme and with GTP hydrolysis. In the urease system, this role is performed by UreG, an accessory protein belonging to the group of homologous P-loop GTPases, often required to complete the biosynthesis of nickel-enzymes. This study is focused on UreG from Helicobacter pylori (HpUreG), a bacterium responsible for gastric ulcers and cancer, infecting large part of the human population, and for which urease is a fundamental virulence factor. The soluble HpUreG was expressed in E. coli and purified to homogeneity. On-line size exclusion chromatography and light scattering indicated that apo-HpUreG exists as a monomer in solution. Circular dichroism, which demonstrated the presence of a well-defined secondary structure, and NMR spectroscopy, which revealed a large number of residues that appear structured on the basis of their backbone amide proton chemical shift dispersion, indicated that, at variance with other UreG proteins so far characterized, this protein is significantly folded in solution. The amino acid sequence of HpUreG is 29% identical to that of HypB from Methanocaldococcus jannaschii, a dimeric zinc-binding GTPase involved in the in vivo assembly of [Ni,Fe]-hydrogenase. A homology-based molecular model of HpUreG was calculated, which allowed us to identify structural and functional features of the protein. Isothermal titration microcalorimetry demonstrated that HpUreG specifically binds 0.5 equivalents of Zn(2+) per monomer (K(d) = 0.33 +/- 0.03 microM), whereas it has 20-fold lower affinity for Ni(2+) (K(d) = 10 +/- 1 microM). Zinc ion binding (but not Ni(2+) binding) causes protein dimerization, as confirmed using light scattering measurements. The structural rearrangement occurring upon Zn(2+)-binding and consequent dimerization was evaluated using circular dichroism and fluorescence spectroscopy. Fully conserved histidine and cysteine residues were identified and their role in zinc binding was verified by site-directed mutagenesis and microcalorimetry. The results are analyzed and discussed with respect to analogous examples of GTPases in nickel metabolism.