Hidden breakpoints in genome alignments

During the course of evolution, an organism's genome can undergo changes that affect the large-scale structure of the genome. These changes include gene gain, loss, duplication, chromosome fusion, fission, and rearrangement. When gene gain and loss occurs in addition to other types of rearrangement, breakpoints of rearrangement can exist that are only detectable by comparison of three or more genomes. An arbitrarily large number of these "hidden" breakpoints can exist among genomes that exhibit no rearrangements in pairwise comparisons. We present an extension of the multichromosomal breakpoint median problem to genomes that have undergone gene gain and loss. We then demonstrate that the median distance among three genomes can be used to calculate a lower bound on the number of hidden breakpoints present. We provide an implementation of this calculation including the median distance, along with some practical improvements on the time complexity of the underlying algorithm. We apply our approach to measure the abundance of hidden breakpoints in simulated data sets under a wide range of evolutionary scenarios. We demonstrate that in simulations the hidden breakpoint counts depend strongly on relative rates of inversion and gene gain/loss. Finally we apply current multiple genome aligners to the simulated genomes, and show that all aligners introduce a high degree of error in hidden breakpoint counts, and that this error grows with evolutionary distance in the simulation. Our results suggest that hidden breakpoint error may be pervasive in genome alignments.Comment: 13 pages, 4 figure

Explosive Percolation in the Human Protein Homology Network

We study the explosive character of the percolation transition in a real-world network. We show that the emergence of a spanning cluster in the Human Protein Homology Network (H-PHN) exhibits similar features to an Achlioptas-type process and is markedly different from regular random percolation. The underlying mechanism of this transition can be described by slow-growing clusters that remain isolated until the later stages of the process, when the addition of a small number of links leads to the rapid interconnection of these modules into a giant cluster. Our results indicate that the evolutionary-based process that shapes the topology of the H-PHN through duplication-divergence events may occur in sudden steps, similarly to what is seen in first-order phase transitions.Comment: 13 pages, 6 figure

Growth performance, in vitro antioxidant properties and chemical composition of the halophyte Limonium algarvense Erben are strongly influenced by the irrigation salinity

Limonium algarvense Erben (sea lavender) is a halophyte species with potential to provide natural ingredients with in vitro antioxidant, anti-inflammatory, neuroprotective and antidiabetic properties. This study reports for the first time the 1) cultivation of sea lavender in greenhouse conditions under irrigation with freshwater (approx. 0 mM NaCl) and saline aquaculture wastewater (300 and 600 mM NaCl), and 2) the influence of the irrigation salinity on the plant performance (e.g growth, number of produced leaves and flowers), in vitro antioxidant properties [radical scavenging activity (DPPH and ABTS), ferric reducing antioxidant power (FRAP), metal chelating properties on copper (CCA) and iron (ICA)], toxicity (in vitro on three mammalian cell lines) and chemical composition (determined by LC-ESI-HRMS/MS). The freshwater-irrigated plants had better growth performance than those irrigated with saltwater. Extracts from wild plants, had the highest antioxidant activity, but those from cultivated ones kept high in vitro antioxidant properties and interesting chemical profile. The flowers' extracts of plants irrigated with 300 mM NaCl had the highest antioxidant activities against DPPH, whereas those from freshwater-irrigated plants were more active on ABTS, CCA and FRAP. Most of the extracts showed nil toxicity. The flowers' extracts displayed the highest diversity of compounds, mainly quercetin, apigenin, luteolin, naringenin and their glycoside derivatives. Moreover, their abundance varied with the irrigation salinity. These data indicate that sea lavender plants can be successfully cultivated in greenhouse conditions under fresh- and saltwater irrigation, maintaining interesting biological and chemical properties.Funding Agency Portuguese Foundation for Science and Technology Portuguese National Budget CCMAR/Multi/04326/2019 GreenVet project ALG-01-0145-FEDER-028876 XtrerneAquaCrops FA-05-2017-028 Lisboa-01-0145-FEDER-022125-RNEM-IST ID/QUI/00100/201 Portuguese Foundation for Science and Technology SFRH/BD/116604/2016 CEECIND/00425/2017info:eu-repo/semantics/publishedVersio

Species-level functional profiling of metagenomes and metatranscriptomes.

Functional profiles of microbial communities are typically generated using comprehensive metagenomic or metatranscriptomic sequence read searches, which are time-consuming, prone to spurious mapping, and often limited to community-level quantification. We developed HUMAnN2, a tiered search strategy that enables fast, accurate, and species-resolved functional profiling of host-associated and environmental communities. HUMAnN2 identifies a community's known species, aligns reads to their pangenomes, performs translated search on unclassified reads, and finally quantifies gene families and pathways. Relative to pure translated search, HUMAnN2 is faster and produces more accurate gene family profiles. We applied HUMAnN2 to study clinal variation in marine metabolism, ecological contribution patterns among human microbiome pathways, variation in species' genomic versus transcriptional contributions, and strain profiling. Further, we introduce 'contributional diversity' to explain patterns of ecological assembly across different microbial community types

Electrochemical integration of graphene with light absorbing copper-based thin films

We present an electrochemical route for the integration of graphene with light sensitive copper-based alloys used in optoelectronic applications. Graphene grown using chemical vapor deposition (CVD) transferred to glass is found to be a robust substrate on which photoconductive Cu_{x}S films of 1-2 um thickness can be deposited. The effect of growth parameters on the morphology and photoconductivity of Cu_{x}S films is presented. Current-voltage characterization and photoconductivity decay experiments are performed with graphene as one contact and silver epoxy as the other

BPGA- an ultra-fast pan-genome analysis pipeline

Recent advances in ultra-high-throughput sequencing technology and metagenomics have led to a paradigm shift in microbial genomics from few genome comparisons to large-scale pan-genome studies at different scales of phylogenetic resolution. Pan-genome studies provide a framework for estimating the genomic diversity of the dataset, determining core (conserved), accessory (dispensable) and unique (strain-specific) gene pool of a species, tracing horizontal gene-flux across strains and providing insight into species evolution. The existing pan genome software tools suffer from various limitations like limited datasets, difficult installation/requirements, inadequate functional features etc. Here we present an ultra-fast computational pipeline BPGA (Bacterial Pan Genome Analysis tool) with seven functional modules. In addition to the routine pan genome analyses, BPGA introduces a number of novel features for downstream analyses like core/pan/MLST (Multi Locus Sequence Typing) phylogeny, exclusive presence/absence of genes in specific strains, subset analysis, atypical G + C content analysis and KEGG & COG mapping of core, accessory and unique genes. Other notable features include minimum running prerequisites, freedom to select the gene clustering method, ultra-fast execution, user friendly command line interface and high-quality graphics outputs. The performance of BPGA has been evaluated using a dataset of complete genome sequences of 28 Streptococcus pyogenes strains

Large-scale comparative genomic ranking of taxonomically restricted genes (TRGs) in bacterial and archaeal genomes

BACKGROUND: Lineage-specific, or taxonomically restricted genes (TRGs), especially those that are species and strain-specific, are of special interest because they are expected to play a role in defining exclusive ecological adaptations to particular niches. Despite this, they are relatively poorly studied and little understood, in large part because many are still orphans or only have homologues in very closely related isolates. This lack of homology confounds attempts to establish the likelihood that a hypothetical gene is expressed and, if so, to determine the putative function of the protein. METHODOLOGY/PRINCIPAL FINDINGS: We have developed "QIPP" ("Quality Index for Predicted Proteins"), an index that scores the "quality" of a protein based on non-homology-based criteria. QIPP can be used to assign a value between zero and one to any protein based on comparing its features to other proteins in a given genome. We have used QIPP to rank the predicted proteins in the proteomes of Bacteria and Archaea. This ranking reveals that there is a large amount of variation in QIPP scores, and identifies many high-scoring orphans as potentially "authentic" (expressed) orphans. There are significant differences in the distributions of QIPP scores between orphan and non-orphan genes for many genomes and a trend for less well-conserved genes to have lower QIPP scores. CONCLUSIONS: The implication of this work is that QIPP scores can be used to further annotate predicted proteins with information that is independent of homology. Such information can be used to prioritize candidates for further analysis. Data generated for this study can be found in the OrphanMine at http://www.genomics.ceh.ac.uk/orphan_mine

Combined use of gliadins and SSRs to analyse the genetic variability of the Spanish collection of cultivated diploid wheat (Triticum monococcum L. ssp. monococcum)

This work studied the combined use of gliadins and SSRs to analyse inter- and intra-accession variability of the Spanish collection of cultivated einkorn (Triticum monococcum L. ssp. monococcum) maintained at the CRF-INIA. In general, gliadin loci presented higher discrimination power than SSRs, reflecting the high variability of the gliadins. The loci on chromosome 6A were the most polymorphic with similar PIC values for both marker systems, showing that these markers are very useful for genetic variability studies in wheat. The gliadin results indicated that the Spanish einkorn collection possessed high genetic diversity, being the differentiation large between varieties and small within them. Some associations between gliadin alleles and geographical and agro-morphological data were found. Agro-morphological relations were also observed in the clusters of the SSRs dendrogram. A high concordance was found between gliadins and SSRs for genotype identification. In addition, both systems provide complementary information to resolve the different cases of intra-accession variability not detected at the agro-morphological level, and to identify separately all the genotypes analysed. The combined use of both genetic markers is an excellent tool for genetic resource evaluation in addition to agro-morphological evaluation

Characterisation of pathogen-specific regions and novel effector candidates in Fusarium oxysporum f. sp. cepae

A reference-quality assembly of Fusarium oxysporum f. sp. cepae (Foc), the causative agent of onion basal rot has been generated along with genomes of additional pathogenic and non-pathogenic isolates of onion. Phylogenetic analysis confirmed a single origin of the Foc pathogenic lineage. Genome alignments with other F. oxysporum ff. spp. and non pathogens revealed high levels of syntenic conservation of core chromosomes but little synteny between lineage specific (LS) chromosomes. Four LS contigs in Foc totaling 3.9 Mb were designated as pathogen-specific (PS). A two-fold increase in segmental duplication events was observed between LS regions of the genome compared to within core regions or from LS regions to the core. RNA-seq expression studies identified candidate effectors expressed in planta, consisting of both known effector homologs and novel candidates. FTF1 and a subset of other transcription factors implicated in regulation of effector expression were found to be expressed in planta

Bioinformatics for the human microbiome project

Microbes inhabit virtually all sites of the human body, yet we know very little about the role they play in our health. In recent years, there has been increasing interest in studying human-associated microbial communities, particularly since microbial dysbioses have now been implicated in a number of human diseases [1]–[3]. Dysbiosis, the disruption of the normal microbial community structure, however, is impossible to define without first establishing what “normal microbial community structure” means within the healthy human microbiome. Recent advances in sequencing technologies have made it feasible to perform large-scale studies of microbial communities, providing the tools necessary to begin to address this question [4], [5]. This led to the implementation of the Human Microbiome Project (HMP) in 2007, an initiative funded by the National Institutes of Health Roadmap for Biomedical Research and constructed as a large, genome-scale community research project [6]. Any such project must plan for data analysis, computational methods development, and the public availability of tools and data; here, we provide an overview of the corresponding bioinformatics organization, history, and results from the HMP (Figure 1).National Institutes of Health (U.S.) (NIH U54HG004969)National Institutes of Health (U.S.) (grant R01HG004885)National Institutes of Health (U.S.) (grant R01HG005975)National Institutes of Health (U.S.) (grant R01HG005969