Analyzing postprandial metabolomics data using multiway models: a simulation study.

BACKGROUND: Analysis of time-resolved postprandial metabolomics data can improve the understanding of metabolic mechanisms, potentially revealing biomarkers for early diagnosis of metabolic diseases and advancing precision nutrition and medicine. Postprandial metabolomics measurements at several time points from multiple subjects can be arranged as a subjects by metabolites by time points array. Traditional analysis methods are limited in terms of revealing subject groups, related metabolites, and temporal patterns simultaneously from such three-way data. RESULTS: We introduce an unsupervised multiway analysis approach based on the CANDECOMP/PARAFAC (CP) model for improved analysis of postprandial metabolomics data guided by a simulation study. Because of the lack of ground truth in real data, we generate simulated data using a comprehensive human metabolic model. This allows us to assess the performance of CP models in terms of revealing subject groups and underlying metabolic processes. We study three analysis approaches: analysis of fasting-state data using principal component analysis, T0-corrected data (i.e., data corrected by subtracting fasting-state data) using a CP model and full-dynamic (i.e., full postprandial) data using CP. Through extensive simulations, we demonstrate that CP models capture meaningful and stable patterns from simulated meal challenge data, revealing underlying mechanisms and differences between diseased versus healthy groups. CONCLUSIONS: Our experiments show that it is crucial to analyze both fasting-state and T0-corrected data for understanding metabolic differences among subject groups. Depending on the nature of the subject group structure, the best group separation may be achieved by CP models of T0-corrected or full-dynamic data. This study introduces an improved analysis approach for postprandial metabolomics data while also shedding light on the debate about correcting baseline values in longitudinal data analysis.

Rice florigens control a common set of genes at the shoot apical meristem including the F-BOX BROADER TILLER ANGLE 1 that regulates tiller angle and spikelet development.

Rice flowering is triggered by transcriptional reprogramming at the shoot apical meristem (SAM) mediated by florigenic proteins produced in leaves in response to changes in photoperiod. Florigens are more rapidly expressed under short days (SDs) compared to long days (LDs) and include the HEADING DATE 3a (Hd3a) and RICE FLOWERING LOCUS T1 (RFT1) phosphatidylethanolamine binding proteins. Hd3a and RFT1 are largely redundant at converting the SAM into an inflorescence, but whether they activate the same target genes and convey all photoperiodic information that modifies gene expression at the SAM is currently unclear. We uncoupled the contribution of Hd3a and RFT1 to transcriptome reprogramming at the SAM by RNA sequencing of dexamethasone-inducible over-expressors of single florigens and wild-type plants exposed to photoperiodic induction. Fifteen highly differentially expressed genes common to Hd3a, RFT1, and SDs were retrieved, 10 of which still uncharacterized. Detailed functional studies on some candidates revealed a role for LOC_Os04g13150 in determining tiller angle and spikelet development and the gene was renamed BROADER TILLER ANGLE 1 (BRT1). We identified a core set of genes controlled by florigen-mediated photoperiodic induction and defined the function of a novel florigen target controlling tiller angle and spikelet development.

Characterizing human postprandial metabolic response using multiway data analysis.

INTRODUCTION: Analysis of time-resolved postprandial metabolomics data can improve our understanding of the human metabolism by revealing similarities and differences in postprandial responses of individuals. Traditional data analysis methods often rely on data summaries or univariate approaches focusing on one metabolite at a time. OBJECTIVES: Our goal is to provide a comprehensive picture in terms of the changes in the human metabolism in response to a meal challenge test, by revealing static and dynamic markers of phenotypes, i.e., subject stratifications, related clusters of metabolites, and their temporal profiles. METHODS: We analyze Nuclear Magnetic Resonance (NMR) spectroscopy measurements of plasma samples collected during a meal challenge test from 299 individuals from the COPSAC2000 cohort using a Nightingale NMR panel at the fasting and postprandial states (15, 30, 60, 90, 120, 150, 240 min). We investigate the postprandial dynamics of the metabolism as reflected in the dynamic behaviour of the measured metabolites. The data is arranged as a three-way array: subjects by metabolites by time. We analyze the fasting state data to reveal static patterns of subject group differences using principal component analysis (PCA), and fasting state-corrected postprandial data using the CANDECOMP/PARAFAC (CP) tensor factorization to reveal dynamic markers of group differences. RESULTS: Our analysis reveals dynamic markers consisting of certain metabolite groups and their temporal profiles showing differences among males according to their body mass index (BMI) in response to the meal challenge. We also show that certain lipoproteins relate to the group difference differently in the fasting vs. dynamic state. Furthermore, while similar dynamic patterns are observed in males and females, the BMI-related group difference is observed only in males in the dynamic state. CONCLUSION: The CP model is an effective approach to analyze time-resolved postprandial metabolomics data, and provides a compact but a comprehensive summary of the postprandial data revealing replicable and interpretable dynamic markers crucial to advance our understanding of changes in the metabolism in response to a meal challenge.

Maternal inflammation during pregnancy is associated with risk of ADHD in children at age 10.

INTRODUCTION: Maternal inflammation during pregnancy may affect early neurodevelopment in offspring as suggested by preclinical and register data. However, clinical evidence for risk of aberrant neurodevelopment later in childhood is scarce. In the population-based COPSAC2010 mother-child cohort, we investigated associations between maternal inflammation levels during pregnancy and the risk of a diagnosis of ADHD as well as the load of ADHD symptoms in the children at age 10. METHODS: The COPSAC2010 cohort consists of 700 mother-child pairs followed prospectively since pregnancy week 24.Maternal high-sensitivity C-Reactive Protein (hs-CRP) level at week 24 of gestation was investigated in relation to child neurodevelopment by age 10 using logistic and linear regression models with extensive confounder adjustment, including socioeconomic status and maternal polygenic risk of ADHD. The children completed a comprehensive examination of neurodevelopment including categorical (i.e., diagnostic) and dimensional (i.e., symptom load) psychopathology using the Kiddie Schedule for Affective Disorders and Schizophrenia Present and Lifetime Version (K-SADS-PL) and parental rated ADHD-Rating Scale (ADHD-RS). RESULTS: A total of 604 (86 %) of the 700 children in the COPSAC2010 cohort participated in the COPSYCH visit at age 10. Sixty-five (10.8 %) fulfilled a research diagnosis of ADHD (16 girls and 49 boys). Higher maternal hs-CRP level in pregnancy at week 24 (median 5.4 mg/L) was significantly associated with increased risk for a diagnosis of ADHD, adjusted OR 1.40, 95 %CI (1.16-1.70), p = 0.001. Additionally, higher maternal hs-CRP was associated with increased ADHD symptom load in the entire cohort, reflected by ADHD-RS raw scores. DISCUSSION: These clinical data demonstrated a robust association of prenatal maternal inflammation assessed by hs-CRP with a diagnosis of ADHD by age 10. Moreover, maternal inflammation was associated with ADHD symptom load in the complete cohort. Identifying inflammation as an important marker will provide a potential target for future increased awareness and prevention during pregnancy thereby ultimately improving neurodevelopmental outcomes in children.

Supplementation With Fish Oil in Pregnancy Reduces Gastroenteritis in Early Childhood.

BACKGROUND: We hypothesized that insufficient intake of fish oil-derived omega-3 long-chain polyunsaturated fatty acids (n-3 LCPUFAs) during pregnancy is a contributing factor to gastroenteritis in early childhood. We examined the effect of n-3 LCPUFA supplementation on gastroenteritis symptoms in the offspring's first 3 years of life. METHODS: This was a double-blinded, randomized controlled trial whereby 736 mothers were administered n-3 LCPUFA or control from pregnancy week 24 until 1 week after birth. We measured the number of days with gastroenteritis, number of episodes with gastroenteritis, and the risk of having a gastroenteritis episode in the first 3 years of life. RESULTS: A median reduction of 2.5 days with gastroenteritis (P = .018) was shown, corresponding to a 14% reduction in the n-3 LCPUFA group compared with controls in the first 3 years of life (P = .037). A reduction in the number of gastroenteritis episodes (P = .027) and a reduced risk of having an episode (hazard ratio, 0.80 [95% confidence interval, .66-.97]; P = .023) were also shown. CONCLUSIONS: Fish oil supplementation from the 24th week of pregnancy led to a reduction in the number of days and episodes with gastroenteritis symptoms in the first 3 years of life. The findings suggest n-3 LCPUFA supplementation as a preventive measure against gastrointestinal infections in early childhood. CLINICAL TRIALS REGISTRATION: NCT00798226.

Next generation sequencing of SARS-CoV-2 genomes: challenges, applications and opportunities.

Various next generation sequencing (NGS) based strategies have been successfully used in the recent past for tracing origins and understanding the evolution of infectious agents, investigating the spread and transmission chains of outbreaks, as well as facilitating the development of effective and rapid molecular diagnostic tests and contributing to the hunt for treatments and vaccines. The ongoing COVID-19 pandemic poses one of the greatest global threats in modern history and has already caused severe social and economic costs. The development of efficient and rapid sequencing methods to reconstruct the genomic sequence of SARS-CoV-2, the etiological agent of COVID-19, has been fundamental for the design of diagnostic molecular tests and to devise effective measures and strategies to mitigate the diffusion of the pandemic. Diverse approaches and sequencing methods can, as testified by the number of available sequences, be applied to SARS-CoV-2 genomes. However, each technology and sequencing approach has its own advantages and limitations. In the current review, we will provide a brief, but hopefully comprehensive, account of currently available platforms and methodological approaches for the sequencing of SARS-CoV-2 genomes. We also present an outline of current repositories and databases that provide access to SARS-CoV-2 genomic data and associated metadata. Finally, we offer general advice and guidelines for the appropriate sharing and deposition of SARS-CoV-2 data and metadata, and suggest that more efficient and standardized integration of current and future SARS-CoV-2-related data would greatly facilitate the struggle against this new pathogen. We hope that our 'vademecum' for the production and handling of SARS-CoV-2-related sequencing data, will contribute to this objective.

Lack of ApoA-I in ApoEKO Mice Causes Skin Xanthomas, Worsening of Inflammation, and Increased Coronary Atherosclerosis in the Absence of Hyperlipidemia.

BACKGROUND: HDL (high-density lipoprotein) and its major protein component, apoA-I (apolipoprotein A-I), play a unique role in cholesterol homeostasis and immunity. ApoA-I deficiency in hyperlipidemic, atheroprone mice was shown to drive cholesterol accumulation and inflammatory cell activation/proliferation. The present study was aimed at investigating the impact of apoA-I deficiency on lipid deposition and local/systemic inflammation in normolipidemic conditions. METHODS: ApoE deficient mice, apoE/apoA-I double deficient (DKO) mice, DKO mice overexpressing human apoA-I, and C57Bl/6J control mice were fed normal laboratory diet until 30 weeks of age. Plasma lipids were quantified, atherosclerosis development at the aortic sinus and coronary arteries was measured, skin ultrastructure was evaluated by electron microscopy. Blood and lymphoid organs were characterized through histological, immunocytofluorimetric, and whole transcriptome analyses. RESULTS: DKO were characterized by almost complete HDL deficiency and by plasma total cholesterol levels comparable to control mice. Only DKO showed xanthoma formation and severe inflammation in the skin-draining lymph nodes, whose transcriptome analysis revealed a dramatic impairment in energy metabolism and fatty acid oxidation pathways. An increased presence of CD4+ T effector memory cells was detected in blood, spleen, and skin-draining lymph nodes of DKO. A worsening of atherosclerosis at the aortic sinus and coronary arteries was also observed in DKO versus apoE deficient. Human apoA-I overexpression in the DKO background was able to rescue the skin phenotype and halt atherosclerosis development. CONCLUSIONS: HDL deficiency, in the absence of hyperlipidemia, is associated with severe alterations of skin morphology, aortic and coronary atherosclerosis, local and systemic inflammation.

VID22 counteracts G-quadruplex-induced genome instability.

Genome instability is a condition characterized by the accumulation of genetic alterations and is a hallmark of cancer cells. To uncover new genes and cellular pathways affecting endogenous DNA damage and genome integrity, we exploited a Synthetic Genetic Array (SGA)-based screen in yeast. Among the positive genes, we identified VID22, reported to be involved in DNA double-strand break repair. vid22Δ cells exhibit increased levels of endogenous DNA damage, chronic DNA damage response activation and accumulate DNA aberrations in sequences displaying high probabilities of forming G-quadruplexes (G4-DNA). If not resolved, these DNA secondary structures can block the progression of both DNA and RNA polymerases and correlate with chromosome fragile sites. Vid22 binds to and protects DNA at G4-containing regions both in vitro and in vivo. Loss of VID22 causes an increase in gross chromosomal rearrangement (GCR) events dependent on G-quadruplex forming sequences. Moreover, the absence of Vid22 causes defects in the correct maintenance of G4-DNA rich elements, such as telomeres and mtDNA, and hypersensitivity to the G4-stabilizing ligand TMPyP4. We thus propose that Vid22 is directly involved in genome integrity maintenance as a novel regulator of G4 metabolism.

Comparative Genomics Reveals Early Emergence and Biased Spatiotemporal Distribution of SARS-CoV-2.

Effective systems for the analysis of molecular data are fundamental for monitoring the spread of infectious diseases and studying pathogen evolution. The rapid identification of emerging viral strains, and/or genetic variants potentially associated with novel phenotypic features is one of the most important objectives of genomic surveillance of human pathogens and represents one of the first lines of defense for the control of their spread. During the COVID 19 pandemic, several taxonomic frameworks have been proposed for the classification of SARS-Cov-2 isolates. These systems, which are typically based on phylogenetic approaches, represent essential tools for epidemiological studies as well as contributing to the study of the origin of the outbreak. Here, we propose an alternative, reproducible, and transparent phenetic method to study changes in SARS-CoV-2 genomic diversity over time. We suggest that our approach can complement other systems and facilitate the identification of biologically relevant variants in the viral genome. To demonstrate the validity of our approach, we present comparative genomic analyses of more than 175,000 genomes. Our method delineates 22 distinct SARS-CoV-2 haplogroups, which, based on the distribution of high-frequency genetic variants, fall into four major macrohaplogroups. We highlight biased spatiotemporal distributions of SARS-CoV-2 genetic profiles and show that seven of the 22 haplogroups (and of all of the four haplogroup clusters) showed a broad geographic distribution within China by the time the outbreak was widely recognized-suggesting early emergence and widespread cryptic circulation of the virus well before its isolation in January 2020. General patterns of genomic variability are remarkably similar within all major SARS-CoV-2 haplogroups, with UTRs consistently exhibiting the greatest variability, with s2m, a conserved secondary structure element of unknown function in the 3'-UTR of the viral genome showing evidence of a functional shift. Although several polymorphic sites that are specific to one or more haplogroups were predicted to be under positive or negative selection, overall our analyses suggest that the emergence of novel types is unlikely to be driven by convergent evolution and independent fixation of advantageous substitutions, or by selection of recombined strains. In the absence of extensive clinical metadata for most available genome sequences, and in the context of extensive geographic and temporal biases in the sampling, many questions regarding the evolution and clinical characteristics of SARS-CoV-2 isolates remain open. However, our data indicate that the approach outlined here can be usefully employed in the identification of candidate SARS-CoV-2 genetic variants of clinical and epidemiological importance.

Critical assessment of bioinformatics methods for the characterization of pathological repeat expansions with single-molecule sequencing data.

A number of studies have reported the successful application of single-molecule sequencing technologies to the determination of the size and sequence of pathological expanded microsatellite repeats over the last 5 years. However, different custom bioinformatics pipelines were employed in each study, preventing meaningful comparisons and somewhat limiting the reproducibility of the results. In this review, we provide a brief summary of state-of-the-art methods for the characterization of expanded repeats alleles, along with a detailed comparison of bioinformatics tools for the determination of repeat length and sequence, using both real and simulated data. Our reanalysis of publicly available human genome sequencing data suggests a modest, but statistically significant, increase of the error rate of single-molecule sequencing technologies at genomic regions containing short tandem repeats. However, we observe that all the methods herein tested, irrespective of the strategy used for the analysis of the data (either based on the alignment or assembly of the reads), show high levels of sensitivity in both the detection of expanded tandem repeats and the estimation of the expansion size, suggesting that approaches based on single-molecule sequencing technologies are highly effective for the detection and quantification of tandem repeat expansions and contractions.

CorGAT: a tool for the functional annotation of SARS-CoV-2 genomes.

SUMMARY: While over 200 000 genomic sequences are currently available through dedicated repositories, ad hoc methods for the functional annotation of SARS-CoV-2 genomes do not harness all currently available resources for the annotation of functionally relevant genomic sites. Here, we present CorGAT, a novel tool for the functional annotation of SARS-CoV-2 genomic variants. By comparisons with other state of the art methods we demonstrate that, by providing a more comprehensive and rich annotation, our method can facilitate the identification of evolutionary patterns in the genome of SARS-CoV-2. AVAILABILITYAND IMPLEMENTATION: Galaxy.http://corgat.cloud.ba.infn.it/galaxy; software: https://github.com/matteo14c/CorGAT/tree/Revision_V1; docker: https://hub.docker.com/r/laniakeacloud/galaxy_corgat. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

VINYL: Variant prIoritizatioN bY survivaL analysis.

MOTIVATION: Clinical applications of genome re-sequencing technologies typically generate large amounts of data that need to be carefully annotated and interpreted to identify genetic variants potentially associated with pathological conditions. In this context, accurate and reproducible methods for the functional annotation and prioritization of genetic variants are of fundamental importance. RESULTS: In this article, we present VINYL, a flexible and fully automated system for the functional annotation and prioritization of genetic variants. Extensive analyses of both real and simulated datasets suggest that VINYL can identify clinically relevant genetic variants in a more accurate manner compared to equivalent state of the art methods, allowing a more rapid and effective prioritization of genetic variants in different experimental settings. As such we believe that VINYL can establish itself as a valuable tool to assist healthcare operators and researchers in clinical genomics investigations. AVAILABILITY AND IMPLEMENTATION: VINYL is available at http://beaconlab.it/VINYL and https://github.com/matteo14c/VINYL. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Aortic Gene Expression Profiles Show How ApoA-I Levels Modulate Inflammation, Lysosomal Activity, and Sphingolipid Metabolism in Murine Atherosclerosis.

OBJECTIVE: HDL (high-density lipoprotein) particles are known to possess several antiatherogenic properties that include the removal of excess cholesterol from peripheral tissues, the maintenance of endothelial integrity, antioxidant, and anti-inflammatory activities. ApoA-I overexpression in apoE-deficient (EKO) mice has been shown to increase HDL levels and to strongly reduce atherosclerosis development. The aim of the study was to investigate gene expression patterns associated with atherosclerosis development in the aorta of EKO mice and how HDL plasma levels relate to gene expression patterns at different stages of atherosclerosis development and with different dietary treatments. Approach and Results: Eight-week-old EKO mice, EKO mice overexpressing human apoA-I, and wild-type mice as controls were fed either normal laboratory or Western diet for 6 or 22 weeks. Cholesterol distribution among lipoproteins was evaluated, and atherosclerosis of the aorta was quantified. High-throughput sequencing technologies were used to analyze the transcriptome of the aorta of the 3 genotypes in each experimental condition. In addition to the well-known activation of inflammation and immune response, the impairment of sphingolipid metabolism, phagosome-lysosome system, and osteoclast differentiation emerged as relevant players in atherosclerosis development. The reduced atherosclerotic burden in the aorta of EKO mice expressing high levels of apoA-I was accompanied by a reduced activation of immune system markers, as well as reduced perturbation of lysosomal activity and a better regulation of the sphingolipid synthesis pathway. CONCLUSIONS: ApoA-I modulates atherosclerosis development in the aorta of EKO mice affecting the expression of pathways additional to those associated with inflammation and immune response.

Evolution of an insect immune barrier through horizontal gene transfer mediated by a parasitic wasp.

Genome sequencing data have recently demonstrated that eukaryote evolution has been remarkably influenced by the acquisition of a large number of genes by horizontal gene transfer (HGT) across different kingdoms. However, in depth-studies on the physiological traits conferred by these accidental DNA acquisitions are largely lacking. Here we elucidate the functional role of Sl gasmin, a gene of a symbiotic virus of a parasitic wasp that has been transferred to an ancestor of the moth species Spodoptera littoralis and domesticated. This gene is highly expressed in circulating immune cells (haemocytes) of larval stages, where its transcription is rapidly boosted by injection of microorganisms into the body cavity. RNAi silencing of Sl gasmin generates a phenotype characterized by a precocious suppression of phagocytic activity by haemocytes, which is rescued when these immune cells are incubated in plasma samples of control larvae, containing high levels of the encoded protein. Proteomic analysis demonstrates that the protein Sl gasmin is released by haemocytes into the haemolymph, where it opsonizes the invading bacteria to promote their phagocytosis, both in vitro and in vivo. Our results show that important physiological traits do not necessarily originate from evolution of pre-existing genes, but can be acquired by HGT events, through unique pathways of symbiotic evolution. These findings indicate that insects can paradoxically acquire selective advantages with the help of their natural enemies.

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.

Transcriptome analysis reveals rice MADS13 as an important repressor of the carpel development pathway in ovules.

In angiosperms, floral homeotic genes encoding MADS-domain transcription factors regulate the development of floral organs. Specifically, members of the SEPALLATA (SEP) and AGAMOUS (AG) subfamilies form higher-order protein complexes to control floral meristem determinacy and to specify the identity of female reproductive organs. In rice, the AG subfamily gene OsMADS13 is intimately involved in the determination of ovule identity, since knock-out mutant plants develop carpel-like structures in place of ovules, resulting in female sterility. Little is known about the regulatory pathways at the base of rice gynoecium development. To investigate molecular mechanisms acting downstream of OsMADS13, we obtained transcriptomes of immature inflorescences from wild-type and Osmads13 mutant plants. Among a total of 476 differentially expressed genes (DEGs), a substantial overlap with DEGs from the SEP-family Osmads1 mutant was found, suggesting that OsMADS1 and OsMADS13 may act on a common set of target genes. Expression studies and preliminary analyses of two up-regulated genes encoding Zinc-finger transcription factors indicated that our dataset represents a valuable resource for the identification of both OsMADS13 target genes and novel players in rice ovule development. Taken together, our study suggests that OsMADS13 is an important repressor of the carpel pathway during ovule development.

ESCRT genes and regulation of developmental signaling.

ESCRT (Endosomal Sorting Complex Required for Transport) proteins have been shown to control an increasing number of membrane-associated processes. Some of these, and prominently regulation of receptor trafficking, profoundly shape signal transduction. Evidence in fungi, plants and multiple animal models support the emerging concept that ESCRTs are main actors in coordination of signaling with the changes in cells and tissues occurring during development and homeostasis. Consistent with their pleiotropic function, ESCRTs are regulated in multiple ways to tailor signaling to developmental and homeostatic needs. ESCRT activity is crucial to correct execution of developmental programs, especially at key transitions, allowing eukaryotes to thrive and preventing appearance of congenital defects.

CONSTANS Imparts DNA Sequence Specificity to the Histone Fold NF-YB/NF-YC Dimer.

Nuclear Factor Y (NF-Y) is a heterotrimeric transcription factor that binds CCAAT elements. The NF-Y trimer is composed of a Histone Fold Domain (HFD) dimer (NF-YB/NF-YC) and NF-YA, which confers DNA sequence specificity. NF-YA shares a conserved domain with the CONSTANS, CONSTANS-LIKE, TOC1 (CCT) proteins. We show that CONSTANS (CO/B-BOX PROTEIN1 BBX1), a master flowering regulator, forms a trimer with Arabidopsis thaliana NF-YB2/NF-YC3 to efficiently bind the CORE element of the FLOWERING LOCUS T promoter. We term this complex NF-CO. Using saturation mutagenesis, electrophoretic mobility shift assays, and RNA-sequencing profiling of co, nf-yb, and nf-yc mutants, we identify CCACA elements as the core NF-CO binding site. CO physically interacts with the same HFD surface required for NF-YA association, as determined by mutations in NF-YB2 and NF-YC9, and tested in vitro and in vivo. The co-7 mutation in the CCT domain, corresponding to an NF-YA arginine directly involved in CCAAT recognition, abolishes NF-CO binding to DNA. In summary, a unifying molecular mechanism of CO function relates it to the NF-YA paradigm, as part of a trimeric complex imparting sequence specificity to HFD/DNA interactions. It is likely that members of the large CCT family participate in similar complexes with At-NF-YB and At-NF-YC, broadening HFD combinatorial possibilities in terms of trimerization, DNA binding specificities, and transcriptional regulation.

Single-cell transcriptomics reveals specific RNA editing signatures in the human brain.

While RNA editing by A-to-I deamination is a requisite for neuronal function in humans, it is under-investigated in single cells. Here we fill this gap by analyzing RNA editing profiles of single cells from the brain cortex of living human subjects. We show that RNA editing levels per cell are bimodally distributed and distinguish between major brain cell types, thus providing new insights into neuronal dynamics.

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.