Drought-Responsive ZmFDL1/MYB94 Regulates Cuticle Biosynthesis and Cuticle-Dependent Leaf Permeability

In all land plants, the outer surface of aerial parts is covered by the cuticle, a complex lipid layer that constitutes a barrier against damage caused by environmental factors and provides protection against nonstomatal water loss. We show in this study that both cuticle deposition and cuticle-dependent leaf permeability during the juvenile phase of plant development are controlled by the maize (Zea mays) transcription factor ZmFUSED LEAVES 1 (FDL1)/MYB94. Biochemical analysis showed altered cutin and wax biosynthesis and deposition in fdl1-1 mutant seedlings at the coleoptile stage. Among cutin compounds, ω-hydroxy fatty acids and polyhydroxy-fatty acids were specifically affected, while the reduction of epicuticular waxes was mainly observed in primary long chain alcohols and, to a minor extent, in long-chain wax esters. Transcriptome analysis allowed the identification of candidate genes involved in lipid metabolism and the assembly of a proposed pathway for cuticle biosynthesis in maize. Lack of ZmFDL1/MYB94 affects the expression of genes located in different modules of the pathway, and we highlighted the correspondence between gene transcriptional variations and biochemical defects. We observed a decrease in cuticle-dependent leaf permeability in maize seedlings exposed to drought as well as abscisic acid treatment, which implies coordinated changes in the transcript levels of ZmFDL1/MYB94 and associated genes. Overall, our results suggest that the response to water stress implies the activation of wax biosynthesis and the involvement of both ZmFDL1/MYB94 and abscisic acid regulatory pathways.

SMRT long reads and Direct Label and Stain optical maps allow the generation of a high-quality genome assembly for the European barn swallow (Hirundo rustica rustica)

Background

The barn swallow (Hirundo rustica) is a migratory bird that has been the focus of a large number of ecological, behavioral, and genetic studies. To facilitate further population genetics and genomic studies, we present a reference genome assembly for the European subspecies (H. r. rustica).

Findings

As part of the Genome10K effort on generating high-quality vertebrate genomes (Vertebrate Genomes Project), we have assembled a highly contiguous genome assembly using single molecule real-time (SMRT) DNA sequencing and several Bionano optical map technologies. We compared and integrated optical maps derived from both the Nick, Label, Repair, and Stain technology and from the Direct Label and Stain (DLS) technology. As proposed by Bionano, DLS more than doubled the scaffold N50 with respect to the nickase. The dual enzyme hybrid scaffold led to a further marginal increase in scaffold N50 and an overall increase of confidence in the scaffolds. After removal of haplotigs, the final assembly is approximately 1.21 Gbp in size, with a scaffold N50 value of more than 25.95 Mbp.

Conclusions

This high-quality genome assembly represents a valuable resource for future studies of population genetics and genomics in the barn swallow and for studies concerning the evolution of avian genomes. It also represents one of the very first genomes assembled by combining SMRT long-read sequencing with the new Bionano DLS technology for scaffolding. The quality of this assembly demonstrates the potential of this methodology to substantially increase the contiguity of genome assemblies.

Pilot study on circulating miRNA signature in children with obesity born small for gestational age and appropriate for gestational age

BACKGROUND

Children born small for gestational age (SGA) are at increased risk of metabolic dysfunction. Dysregulation of specific microRNAs (miRNAs) contributes to aberrant gene expression patterns underlying metabolic dysfunction.

OBJECTIVE

We aimed to determine and compare circulating miRNA (c-miRNA) profile of SGA and appropriate for gestational age (AGA) children with obesity and with normal weight, in order to identify biomarkers for early detection of increased risk of developing metabolic dysfunction in SGA and AGA children with obesity.

METHODS

Small non-coding RNAs from serum of 15 SGA children with obesity (OB-SGA), 10 SGA children with normal weight (NW-SGA), 17 AGA children with obesity (OB-AGA) and 12 AGA children with normal weight (NW-AGA) (mean age 11.2 ± 2.6) have been extracted and sequenced in order to detect and quantify miRNA expression profiles.

RESULTS

RNA-seq analyses showed 28 miRNAs dysregulated in OB-SGA vs. NW-SGA and 19 miRNAs dysregulated in OB-AGA vs. NW-AGA. Among these, miR-92a-3p, miR-122-5p, miR-423-5p, miR-484, miR-486-3p and miR-532-5p were up regulated, and miR-181b-5p was down regulated in both OB-SGA and OB-AGA compared with normal weight counterparts. Pathway analysis and miRNA target prediction suggested that these miRNAs were particularly involved in insulin signalling, glucose transport, insulin resistance, cholesterol and lipid metabolism.

CONCLUSION

We identified a specific profile of c-miRNAs in SGA and AGA children with obesity compared with SGA and AGA children with normal weight. These c-miRNAs could represent specific biomarkers for early detection of increased risk of developing metabolic dysfunction in SGA and AGA children with obesity.

CoVaCS: a consensus variant calling system

Background

The advent and ongoing development of next generation sequencing technologies (NGS) has led to a rapid increase in the rate of human genome re-sequencing data, paving the way for personalized genomics and precision medicine. The body of genome resequencing data is progressively increasing underlining the need for accurate and time-effective bioinformatics systems for genotyping - a crucial prerequisite for identification of candidate causal mutations in diagnostic screens.

Results

Here we present CoVaCS, a fully automated, highly accurate system with a web based graphical interface for genotyping and variant annotation. Extensive tests on a gold standard benchmark data-set -the NA12878 Illumina platinum genome- confirm that call-sets based on our consensus strategy are completely in line with those attained by similar command line based approaches, and far more accurate than call-sets from any individual tool. Importantly our system exhibits better sensitivity and higher specificity than equivalent commercial software.

Conclusions

CoVaCS offers optimized pipelines integrating state of the art tools for variant calling and annotation for whole genome sequencing (WGS), whole-exome sequencing (WES) and target-gene sequencing (TGS) data. The system is currently hosted at Cineca, and offers the speed of a HPC computing facility, a crucial consideration when large numbers of samples must be analysed. Importantly, all the analyses are performed automatically allowing high reproducibility of the results. As such, we believe that CoVaCS can be a valuable tool for the analysis of human genome resequencing studies. CoVaCS is available at: https://bioinformatics.cineca.it/covacs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4508-1) contains supplementary material, which is available to authorized users.

A-GAME: improving the assembly of pooled functional metagenomics sequence data

Background

Expression screening of environmental DNA (eDNA) libraries is a popular approach for the identification and characterization of novel microbial enzymes with promising biotechnological properties. In such “functional metagenomics” experiments, inserts, selected on the basis of activity assays, are sequenced with high throughput sequencing technologies. Assembly is followed by gene prediction, annotation and identification of candidate genes that are subsequently evaluated for biotechnological applications.

Results

Here we present A-GAME (A GAlaxy suite for functional MEtagenomics), a web service incorporating state of the art tools and workflows for the analysis of eDNA sequence data. We illustrate the potential of A-GAME workflows using real functional metagenomics data, showing that they outperform alternative metagenomics assemblers. Dedicated tools available in A-GAME allow efficient analysis of pooled libraries and rapid identification of candidate genes, reducing sequencing costs and saving the need for laborious manual annotation.

Conclusion

In conclusion, we believe A-GAME will constitute a valuable resource for the functional metagenomics community.

A-GAME is publicly available at http://beaconlab.it/agame

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4369-z) contains supplementary material, which is available to authorized users.

Transcriptomic analysis of nickel exposure in Sphingobium sp. ba1 cells using RNA-seq

Nickel acts as cofactor for a number of enzymes of many bacteria species. Its homeostasis is ensured by proteins working as ion efflux or accumulation systems. These mechanisms are also generally adopted to counteract life-threatening high extra-cellular Ni²⁺ concentrations. Little is known regarding nickel tolerance in the genus Sphingobium. We studied the response of the novel Sphingobium sp. ba1 strain, able to adapt to high Ni²⁺ concentrations. Differential gene expression in cells cultured in 10 mM Ni²⁺, investigated by RNA-seq analysis, identified 118 differentially expressed genes. Among the 90 up-regulated genes, a cluster including genes coding for nickel and other metal ion efflux systems (similar to either cnrCBA, nccCBA or cznABC) and for a NreB-like permease was found. Comparative analyses among thirty genomes of Sphingobium species show that this cluster is conserved only in two cases, while in the other genomes it is partially present or even absent. The differential expression of genes encoding proteins which could also work as Ni²⁺-accumulators (HupE/UreJ-like protein, NreA and components of TonB-associated transport and copper-homeostasis systems) was also detected. The identification of Sphingobium sp. ba1 strain adaptive mechanisms to nickel ions, can foster its possible use for biodegradation of poly-aromatic compounds in metal-rich environments.

A novel group of IncQ1 plasmids conferring multidrug resistance

The IncQ is a group of non-conjugative but mobilisable plasmids that are found and stably maintained in a wide range of bacteria contributing to the spread of antimicrobial resistance genes and to the insurgence of multidrug resistant bacteria. Here we report the identification, in clinical Salmonella Typhimurium strains, of an IncQ1 plasmid (pNUC) which confers resistance to sulfamethoxazole, streptomycin and tetracycline through the presence of sul2, strAB and tetA genes, respectively. pNUC was detected in five multidrug resistant S. Typhimurium strains collected in Southern Italy from various hospitals and years of isolation. Bioinformatics analyses highlighted the presence of pNUC-like plasmids in pathogenic bacteria of various Enterobacteriaceae genera or species. Taken as a whole, these plasmids constitute a novel group of IncQ1 plasmids that might have originated through recombination events between a tetR-tetA gene cluster (possibly derived from a Tn1721) and a recipient IncQ1 plasmid related to RSF1010. Our findings raise concerns regarding the possible contribution of the newly identified group of IncQ1 plasmids to the spread of tetracycline resistance.

RNA sequencing of Populus x canadensis roots identifies key molecular mechanisms underlying physiological adaption to excess zinc

Populus x canadensis clone I-214 exhibits a general indicator phenotype in response to excess Zn, and a higher metal uptake in roots than in shoots with a reduced translocation to aerial parts under hydroponic conditions. This physiological adaptation seems mainly regulated by roots, although the molecular mechanisms that underlie these processes are still poorly understood. Here, differential expression analysis using RNA-sequencing technology was used to identify the molecular mechanisms involved in the response to excess Zn in root. In order to maximize specificity of detection of differentially expressed (DE) genes, we consider the intersection of genes identified by three distinct statistical approaches (61 up- and 19 down-regulated) and validate them by RT-qPCR, yielding an agreement of 93% between the two experimental techniques. Gene Ontology (GO) terms related to oxidation-reduction processes, transport and cellular iron ion homeostasis were enriched among DE genes, highlighting the importance of metal homeostasis in adaptation to excess Zn by P. x canadensis clone I-214. We identified the up-regulation of two Populus metal transporters (ZIP2 and NRAMP1) probably involved in metal uptake, and the down-regulation of a NAS4 gene involved in metal translocation. We identified also four Fe-homeostasis transcription factors (two bHLH38 genes, FIT and BTS) that were differentially expressed, probably for reducing Zn-induced Fe-deficiency. In particular, we suggest that the down-regulation of FIT transcription factor could be a mechanism to cope with Zn-induced Fe-deficiency in Populus. These results provide insight into the molecular mechanisms involved in adaption to excess Zn in Populus spp., but could also constitute a starting point for the identification and characterization of molecular markers or biotechnological targets for possible improvement of phytoremediation performances of poplar trees.

Addressing the role of microRNAs in reprogramming leaf growth during drought stress in Brachypodium distachyon

Plant responses to drought are regulated by complex genetic and epigenetic networks leading to rapid reprogramming of plant growth. miRNAs have been widely indicated as key players in the regulation of growth and development. The role of miRNAs in drought response was investigated in young leaves of Brachypodium distachyon, a drought-tolerant monocot model species. Adopting an in vivo drought assay, shown to cause a dramatic reduction in leaf size, mostly due to reduced cell expansion, small RNA libraries were produced from proliferating and expanding leaf cells. Next-generation sequencing data were analyzed using an in-house bioinformatics pipeline allowing the identification of 66 annotated miRNA genes and 122 new high confidence predictions greatly expanding the number of known Brachypodium miRNAs. In addition, we identified four TAS3 loci and a large number of siRNA-producing loci that show characteristics suggesting that they may represent young miRNA genes. Most miRNAs showed a high expression level, consistent with their involvement in early leaf development and cell identity. Proliferating and expanding leaf cells respond differently to drought treatment and differential expression analyses suggest novel evidence for an miRNA regulatory network controlling cell division in both normal and stressed conditions and demonstrate that drought triggers a genetic reprogramming of leaf growth in which miRNAs are deeply involved.

The estimation of relative site variability among aligned homologous protein sequences

MOTIVATION

Maximum likelihood-based methods to estimate site by site substitution rate variability in aligned homologous protein sequences rely on the formulation of a phylogenetic tree and generally assume that the patterns of relative variability follow a pre-determined distribution. We present a phylogenetic tree-independent method to estimate the relative variability of individual sites within large datasets of homologous protein sequences. It is based upon two simple assumptions. Firstly that substitutions observed between two closely related sequences are likely, in general, to occur at the most variable sites. Secondly that non-conservative amino acid substitutions tend to occur at more variable sites. Our methodology makes no assumptions regarding the underlying pattern of relative variability between sites.

RESULTS

We have compared, using data simulated under a non-gamma distributed model, the performance of this approach to that of a maximum likelihood method that assumes gamma distributed rates. At low mean rates of evolution our method inferred site by site relative substitution rates more accurately than the maximum likelihood approach in the absence of prior assumptions about the relationships between sequences. Our method does not directly account for the effects of mutational saturation, However, we have incorporated an 'ad-hoc' modification that allows the accurate estimation of relative site variability in fast evolving and saturated datasets.

Unique phylogenetic relationships of glucokinase and glucosephosphate isomerase of the amitochondriate eukaryotes Giardia intestinalis, Spironucleus barkhanus and Trichomonas vaginalis

Glucokinase (GK) and glucosephosphate isomerase (GPI), the first two enzymes of the glycolytic pathway of the diplomonads Giardia intestinalis and Spironucleus barkhanus, Type I amitochondriate eukaryotes, were sequenced. GPI of the parabasalid Trichomonas vaginalis was also sequenced. The diplomonad GKs belong to a family of specific GKs present in cyanobacteria, in some proteobacteria and also in T. vaginalis, a Type II amitochondriate protist. These enzymes are not part of the hexokinase family, which is broadly distributed among eukaryotes, including the Type I amitochondriate parasite Entamoeba histolytica. G. intestinalis GK expressed in Escherichia coli was specific for glucose and glucosamine, as are its eubacterial homologs. The sequence of diplomonad and trichomonad GPIs formed a monophyletic group more closely related to cyanobacterial and chloroplast sequences than to cytosolic GPIs of other eukaryotes and prokaryotes. The findings show that certain enzymes of the energy metabolism of these amitochondriate protists originated from sources different than those of other eukaryotes. The observation that the two diplomonads and T. vaginalis share the same unusual GK and GPI is consistent with gene trees that suggest a close relationship between diplomonads and parabasalids. The intriguing relationships of these enzymes to cyanobacterial (and chloroplast) enzymes might reflect horizontal gene transfer between the common ancestor of the diplomonad and parabasalid lineages and the ancestor of cyanobacteria.

Iron hydrogenases and the evolution of anaerobic eukaryotes

Hydrogenases, oxygen-sensitive enzymes that can make hydrogen gas, are key to the function of hydrogen-producing organelles (hydrogenosomes), which occur in anaerobic protozoa scattered throughout the eukaryotic tree. Hydrogenases also play a central role in the hydrogen and syntrophic hypotheses for eukaryogenesis. Here, we show that sequences related to iron-only hydrogenases ([Fe] hydrogenases) are more widely distributed among eukaryotes than reports of hydrogen production have suggested. Genes encoding small proteins which contain conserved structural features unique to [Fe] hydrogenases were identified on all well-surveyed aerobic eukaryote genomes. Longer sequences encoding [Fe] hydrogenases also occur in the anaerobic eukaryotes Entamoeba histolytica and Spironucleus barkhanus, both of which lack hydrogenosomes. We also identified a new [Fe] hydrogenase sequence from Trichomonas vaginalis, bringing the total of [Fe] hydrogenases reported for this organism to three, all of which may function within its hydrogenosomes. Phylogenetic analysis and hypothesis testing using likelihood ratio tests and parametric bootstrapping suggest that the [Fe] hydrogenases in anaerobic eukaryotes are not monophyletic. Iron-only hydrogenases from Entamoeba, Spironucleus, and Trichomonas are plausibly monophyletic, consistent with the hypothesis that a gene for [Fe] hydrogenase was already present on the genome of the common, perhaps also anaerobic, ancestor of these phylogenetically distinct eukaryotes. Trees where the [Fe] hydrogenase from the hydrogenosomal ciliate Nyctotherus was constrained to be monophyletic with the other eukaryote sequences were rejected using a likelihood ratio test of monophyly. In most analyses, the Nyctotherus sequence formed a sister group with a [Fe] hydrogenase on the genome of the eubacterium Desulfovibrio vulgaris. Thus, it is possible that Nyctotherus obtained its hydrogenosomal [Fe] hydrogenase from a different source from Trichomonas for its hydrogenosomes. We find no support for the hypothesis that components of the Nyctotherus [Fe] hydrogenase fusion protein derive from the mitochondrial respiratory chain.

A single eubacterial origin of eukaryotic pyruvate: ferredoxin oxidoreductase genes: implications for the evolution of anaerobic eukaryotes

The iron sulfur protein pyruvate: ferredoxin oxidoreductase (PFO) is central to energy metabolism in amitochondriate eukaryotes, including those with hydrogenosomes. Thus, revealing the evolutionary history of PFO is critical to understanding the origin(s) of eukaryote anaerobic energy metabolism. We determined a complete PFO sequence for Spironucleus barkhanus, a large fragment of a PFO sequence from Clostridium pasteurianum, and a fragment of a new PFO from Giardia lamblia. Phylogenetic analyses of eubacterial and eukaryotic PFO genes suggest a complex history for PFO, including possible gene duplications and horizontal transfers among eubacteria. Our analyses favor a common origin for eukaryotic cytosolic and hydrogenosomal PFOs from a single eubacterial source, rather than from separate horizontal transfers as previously suggested. However, with the present sampling of genes and species, we were unable to infer a specific eubacterial sister group for eukaryotic PFO. Thus, we find no direct support for the published hypothesis that the donor of eukaryote PFO was the common alpha-proteobacterial ancestor of mitochondria and hydrogenosomes. We also report that several fungi and protists encode proteins with PFO domains that are likely monophyletic with PFOs from anaerobic protists. In Saccharomyces cerevisiae, PFO domains combine with fragments of other redox proteins to form fusion proteins which participate in methionine biosynthesis. Our results are consistent with the view that PFO, an enzyme previously considered to be specific to energy metabolism in amitochondriate protists, was present in the common ancestor of contemporary eukaryotes and was retained, wholly or in part, during the evolution of oxygen-dependent and mitochondrion-bearing lineages.

Molecular data suggest an early acquisition of the mitochondrion endosymbiont

The three deepest branching eucaryotic lineages in small subunit ribosomal RNA phylogenies are the amitochondriate Microspora, Metamonada and Parabasala. They are followed by either the Euglenozoa (e.g. Euglena and Trypanosoma) or the Percolozoa as the first mitochondria-containing eucaryotes. To investigate the hypothesis of an even earlier timing of the mitochondrion endosymbiosis we have amplified a partial cpn-60 coding region from the parabasalid Trichomonas vaginalis and the first such sequence from a percolozoan, Naegleria fowleri. Analysis of predicted protein sequences reveals a high degree of sequence similarity (> or = 40%) with a selection of published bacterial and mitochondrial cpn-60s for both taxa. Both sequences were recovered within a strongly supported monophyletic group, otherwise defined by mitochondrial sequences, which systematically clustered with alpha-proteobacteria. These results provide compelling evidence that the ancestor of T. vaginalis once contained the endosymbiont which gave rise to mitochondria, and suggest that this symbiosis probably occurred before the Trichomonas lineage diverged from the main eukaryote trunk. It also makes feasible the published hypothesis that the Trichomonas hydrogenosome might represent a biochemically modified mitochondrion. Analysis of the N. fowleri cpn-60 did not support the hypothesis that the mitochondrion-containing Percolozoa represent an earlier branch in the cpn-60 tree than Trichomonas or Trypanosoma.

Comparison of three-day butoconazole treatment with seven-day miconazole treatment for vulvovaginal candidiasis

In this multicenter, parallel, randomized, investigator-blind trial, we compared the safety and efficacy of a three-day regimen of 2% butoconazole vaginal cream with those of a seven-day regimen of 2% miconazole vaginal cream. Enrolled were 271 nonpregnant women with vulvovaginal candidiasis. Each patient administered her assigned study medication to the posterior vaginal fornix for three or seven consecutive nights. All 271 patients were included in the safety evaluation, and 225 (111 receiving butoconazole and 114 receiving miconazole) were included in the efficacy evaluation. Eight to ten days after treatment completion, 88% of the butoconazole-treated patients and 91% of the miconazole-treated patients were Candida negative; 80% of the butoconazole-treated patients and 82% of the miconazole-treated patients were considered clinically cured. Thirty days after treatment completion, 73% of the butoconazole-treated patients and 69% of the miconazole-treated patients remained Candida negative; 78% of the butoconazole-treated patients and 80% of the miconazole-treated patients remained free of clinical symptoms of vulvovaginitis. None of the differences between the two treatment groups was statistically significant. Six patients (four receiving butoconazole and two receiving miconazole) reported increased symptoms of vulvovaginal irritation, and three of them (two receiving butoconazole and one receiving miconazole) withdrew from the trial. Thus, the efficacy and safety of the three-day butoconazole treatment regimen were equivalent to those of the seven-day miconazole treatment regimen. The advantage of the shorter butoconazole treatment is increased patient compliance with maintenance of high efficacy.