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.

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.

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 .

Bioinformatics approaches for genomics and post genomics applications of next-generation sequencing.

Technical advances such as the development of molecular cloning, Sanger sequencing, PCR and oligonucleotide microarrays are key to our current capacity to sequence, annotate and study complete organismal genomes. Recent years have seen the development of a variety of so-called 'next-generation' sequencing platforms, with several others anticipated to become available shortly. The previously unimaginable scale and economy of these methods, coupled with their enthusiastic uptake by the scientific community and the potential for further improvements in accuracy and read length, suggest that these technologies are destined to make a huge and ongoing impact upon genomic and post-genomic biology. However, like the analysis of microarray data and the assembly and annotation of complete genome sequences from conventional sequencing data, the management and analysis of next-generation sequencing data requires (and indeed has already driven) the development of informatics tools able to assemble, map, and interpret huge quantities of relatively or extremely short nucleotide sequence data. Here we provide a broad overview of bioinformatics approaches that have been introduced for several genomics and functional genomics applications of next-generation sequencing.

The evolutionary history of the polyQ tract in huntingtin sheds light on its functional pro-neural activities.

Huntington's disease is caused by a pathologically long (>35) CAG repeat located in the first exon of the Huntingtin gene (HTT). While pathologically expanded CAG repeats are the focus of extensive investigations, non-pathogenic CAG tracts in protein-coding genes are less well characterized. Here, we investigated the function and evolution of the physiological CAG tract in the HTT gene. We show that the poly-glutamine (polyQ) tract encoded by CAGs in the huntingtin protein (HTT) is under purifying selection and subjected to stronger selective pressures than CAG-encoded polyQ tracts in other proteins. For natural selection to operate, the polyQ must perform a function. By combining genome-edited mouse embryonic stem cells and cell assays, we show that small variations in HTT polyQ lengths significantly correlate with cells' neurogenic potential and with changes in the gene transcription network governing neuronal function. We conclude that during evolution natural selection promotes the conservation and purity of the CAG-encoded polyQ tract and that small increases in its physiological length influence neural functions of HTT. We propose that these changes in HTT polyQ length contribute to evolutionary fitness including potentially to the development of a more complex nervous system.

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.

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.