Comparative genomics

A phylogenetic approach to comparative genomics

Comparative genomics, whereby the genomes of different species are compared, has the potential to address broad and fundamental questions at the intersection of genetics and evolution. However, species, genomes and genes cannot be considered as independent data points within statistical tests. Closely related species tend to be similar because they share genes by common descent, which must be accounted for in analyses. This problem of non-independence may be exacerbated when examining genomes or genes but can be addressed by applying phylogeny-based methods to comparative genomic analyses. Here, we review how controlling for phylogeny can change the conclusions of comparative genomics studies. We address common questions on how to apply these methods and illustrate how they can be used to test causal hypotheses. The combination of rapidly expanding genomic datasets and phylogenetic comparative methods is set to revolutionize the biological insights possible from comparative genomic studies.

Data mining and molecular dynamics analysis to detect HIV-1 reverse transcriptase RNase H activity inhibitor

HIV-1 is a deadly virus that affects millions of people worldwide. In this study, we aimed to inhibit viral replication by targeting one of the HIV-1 proteins and identifying a new drug candidate. We used data mining and molecular dynamics methods on HIV-1 genomes. Based on MAUVE analysis, we selected the RNase H activity of the reverse transcriptase (R.T) enzyme as a potential target due to its low mutation rate and high conservation level. We screened about 94,000 small molecule inhibitors by virtual screening. We validated the hit compounds' stability and binding free energy through molecular dynamics simulations and MM/PBSA. Phomoarcherin B, known for its anticancer properties, emerged as the best candidate and showed potential as an HIV-1 reverse transcriptase RNase H activity inhibitor. This study presents a new target and drug candidate for HIV-1 treatment. However, in vitro and in vivo tests are required. Also, the effect of RNase H activity on viral replication and the interaction of Phomoarcherin B with other HIV-1 proteins should be investigated.

Self-cascade deoxynivalenol detoxification by an artificial enzyme with bifunctions of dehydrogenase and aldo/keto reductase from genome mining

Due to the severe health risks for human and animal caused by the intake of toxic deoxynivalenol (DON) derived from Fusarium species, elimination DON in food and feed has been initiated as a critical issue. Enzymatic cascade catalysis by dehydrogenase and aldo-keto reductase represents a fascinating strategy for DON detoxification. Here, one quinone-dpendent alcohol dehydrogenase DADH oxidized DON into less-toxic 3-keto-DON and NADPH-dependent aldo-keto reductase AKR13B3 reduced 3-keto-DON into relatively non-toxic 3-epi-DON were identified from Devosia strain A6-243, indicating that degradation of DON on C3 are two-step sequential cascade processes. To establish the bifunctions, fusion enzyme linking DADH and AKR13B3 was successfully assembled to promote one-step DON degradations with accelerated specific activity and efficiency, resulting 93.29 % of DON removal rate in wheat sample. Three-dimensional simulation analysis revealed that the bifunctional enzyme forms an artificial intramolecular channel to minimize the distance of intermediate from DADH to AKR13B3 for two-step enzymatic reactions, and thereby accelerates this enzymatic process. As the first report of directing single step DON detoxification by an interesting bifunctional artificial enzyme, this work revealed a facile and eco-friendly approach to detoxify DON with application potential and gave valuable insights into execute other mycotoxin detoxification for ensuring food safety.

A vast evolutionarily transient translatome contributes to phenotype and fitness

Translation is the process by which ribosomes synthesize proteins. Ribosome profiling recently revealed that many short sequences previously thought to be noncoding are pervasively translated. To identify protein-coding genes in this noncanonical translatome, we combine an integrative framework for extremely sensitive ribosome profiling analysis, iRibo, with high-powered selection inferences tailored for short sequences. We construct a reference translatome for Saccharomyces cerevisiae comprising 5,400 canonical and almost 19,000 noncanonical translated elements. Only 14 noncanonical elements were evolving under detectable purifying selection. A representative subset of translated elements lacking signatures of selection demonstrated involvement in processes including DNA repair, stress response, and post-transcriptional regulation. Our results suggest that most translated elements are not conserved protein-coding genes and contribute to genotype-phenotype relationships through fast-evolving molecular mechanisms.

Comparative Genomics Analysis of Habitat Adaptation by Lactobacillus kefiranofaciens

Lactobacillus kefiranofaciens is often found in fermented dairy products. Many strains of this species have probiotic properties, contributing to the regulation of immune metabolism and intestinal flora. This species was added to the list of lactic acid bacteria that can be added to food in China, in 2020. However, research on the genomics of this species is scarce. In this study we undertook whole genome sequencing analysis of 82 strains of L. kefiranofaciens from different habitats, of which 9 strains were downloaded from the NCBI RefSeq (National Center for Biotechnology Information RefSeq). The mean genome size of the 82 strains was 2.05 ± 0.25 Mbp, and the mean DNA G + C content was 37.47 ± 0.42%. The phylogenetic evolutionary tree for the core genes showed that all strains belonged to five clades with clear aggregation in relation to the isolation habitat; this indicated that the genetic evolution of L. kefiranofaciens was correlated to the isolation habitat. Analysis of the annotation results identified differences in the functional genes, carbohydrate active enzymes (CAZy) and bacteriocins amongst different isolated strains, which were related to the environment. Isolates from kefir grains had more enzymes for cellulose metabolism and a better ability to use vegetative substrates for fermentation, which could be used in feed production. Isolates from kefir grains also had fewer kinds of bacteriocin than isolates from sour milk and koumiss; helveticin J and lanthipeptide class I were not found in the isolates from kefir grains. The genomic characteristics and evolutionary process of L. kefiranofaciens were analyzed by comparative genomics and this paper explored the differences in the functional genes amongst the strains, aiming to provide a theoretical basis for the research and development of L. kefiranofaciens.

Integrative comparative analysis of avian chromosome evolution by in-silico mapping of the gene ontology of homologous synteny blocks and evolutionary breakpoint regions

Avian chromosomes undergo more intra- than interchromosomal rearrangements, which either induce or are associated with genome variations among birds. Evolving from a common ancestor with a karyotype not dissimilar from modern chicken, two evolutionary elements characterize evolutionary change: homologous synteny blocks (HSBs) constitute common conserved parts at the sequence level, while evolutionary breakpoint regions (EBRs) occur between HSBs, defining the points where rearrangement occurred. Understanding the link between the structural organization and functionality of HSBs and EBRs provides insight into the mechanistic basis of chromosomal change. Previously, we identified gene ontology (GO) terms associated with both; however, here we revisit our analyses in light of newly developed bioinformatic algorithms and the chicken genome assembly galGal6. We aligned genomes available for six birds and one lizard species, identifying 630 HSBs and 19 EBRs. We demonstrate that HSBs hold vast functionality expressed by GO terms that have been largely conserved through evolution. Particularly, we found that genes within microchromosomal HSBs had specific functionalities relevant to neurons, RNA, cellular transport and embryonic development, and other associations. Our findings suggest that microchromosomes may have conserved throughout evolution due to the specificity of GO terms within their HSBs. The detected EBRs included those found in the genome of the anole lizard, meaning they were shared by all saurian descendants, with others being unique to avian lineages. Our estimate of gene richness in HSBs supported the fact that microchromosomes contain twice as many genes as macrochromosomes.

Contributions of Adaptive Laboratory Evolution towards the Enhancement of the Biotechnological Potential of Non-Conventional Yeast Species

Changes in biological properties over several generations, induced by controlling short-term evolutionary processes in the laboratory through selective pressure, and whole-genome re-sequencing, help determine the genetic basis of microorganism's adaptive laboratory evolution (ALE). Due to the versatility of this technique and the imminent urgency for alternatives to petroleum-based strategies, ALE has been actively conducted for several yeasts, primarily using the conventional species Saccharomyces cerevisiae, but also non-conventional yeasts. As a hot topic at the moment since genetically modified organisms are a debatable subject and a global consensus on their employment has not yet been attained, a panoply of new studies employing ALE approaches have emerged and many different applications have been exploited in this context. In the present review, we gathered, for the first time, relevant studies showing the ALE of non-conventional yeast species towards their biotechnological improvement, cataloging them according to the aim of the study, and comparing them considering the species used, the outcome of the experiment, and the employed methodology. This review sheds light on the applicability of ALE as a powerful tool to enhance species features and improve their performance in biotechnology, with emphasis on the non-conventional yeast species, as an alternative or in combination with genome editing approaches.

RHOA Therapeutic Targeting in Hematological Cancers

Primarily identified as an important regulator of cytoskeletal dynamics, the small GTPase Ras homolog gene family member A (RHOA) has been implicated in the transduction of signals regulating a broad range of cellular functions such as cell survival, migration, adhesion and proliferation. Deregulated activity of RHOA has been linked to the growth, progression and metastasis of various cancer types. Recent cancer genome-wide sequencing studies have unveiled both RHOA gain and loss-of-function mutations in primary leukemia/lymphoma, suggesting that this GTPase may exert tumor-promoting or tumor-suppressive functions depending on the cellular context. Based on these observations, RHOA signaling represents an attractive therapeutic target for the development of selective anticancer strategies. In this review, we will summarize the molecular mechanisms underlying RHOA GTPase functions in immune regulation and in the development of hematological neoplasms and will discuss the current strategies aimed at modulating RHOA functions in these diseases.

Phylogenomics provides insights into the evolution of cactophily and host plant shifts in Drosophila

Cactophilic species of the Drosophila buzzatii cluster (repleta group) comprise an excellent model group to investigate genomic changes underlying adaptation to extreme climate conditions and host plants. In particular, these species form a tractable system to study the transition from chemically simpler breeding sites (like prickly pears of the genus Opuntia) to chemically more complex hosts (columnar cacti). Here, we report four highly contiguous genome assemblies of three species of the buzzatii cluster. Based on this genomic data and inferred phylogenetic relationships, we identified candidate taxonomically restricted genes (TRGs) likely involved in the evolution of cactophily and cactus host specialization. Functional enrichment analyses of TRGs within the buzzatii cluster identified genes involved in detoxification, water preservation, immune system response, anatomical structure development, and morphogenesis. In contrast, processes that regulate responses to stress, as well as the metabolism of nitrogen compounds, transport, and secretion were found in the set of species that are columnar cacti dwellers. These findings are in line with the hypothesis that those genomic changes brought about key mechanisms underlying the adaptation of the buzzatii cluster species to arid regions in South America.

Molecular Breeding and Drought Tolerance in Chickpea

Cicer arietinum L. is the third greatest widely planted imperative pulse crop worldwide, and it belongs to the Leguminosae family. Drought is the utmost common abiotic factor on plants, distressing their water status and limiting their growth and development. Chickpea genotypes have the natural ability to fight drought stress using certain strategies viz., escape, avoidance and tolerance. Assorted breeding methods, including hybridization, mutation, and marker-aided breeding, genome sequencing along with omics approaches, could be used to improve the chickpea germplasm lines(s) against drought stress. Root features, for instance depth and root biomass, have been recognized as the greatest beneficial morphological factors for managing terminal drought tolerance in the chickpea. Marker-aided selection, for example, is a genomics-assisted breeding (GAB) strategy that can considerably increase crop breeding accuracy and competence. These breeding technologies, notably marker-assisted breeding, omics, and plant physiology knowledge, underlined the importance of chickpea breeding and can be used in future crop improvement programmes to generate drought-tolerant cultivars(s).

Insufficiency in functional genomics studies, data, and applications: A case study of bio-prospecting research in ruminant microbiome

Over the last two decades, biotechnology has advanced at a rapid pace, propelled by the incorporation of bio-products into various aspects of pharmaceuticals, industry, and the environment. These developments have sparked interest in the bioprospecting of microorganisms and their products in a variety of niche environments. Furthermore, the use of omics technologies has greatly aided our analyses of environmental samples by elucidating the microbial ecological framework, biochemical pathways, and bio-products. However, the more often overemphasis on taxonomic identification in most research publications, as well as the data associated with such studies, is detrimental to immediate industrial and commercial applications. This review identifies several factors that contribute to the complexity of sequence data analysis as potential barriers to the pragmatic application of functional genomics, utilizing recent research on ruminants to demonstrate these limitations in the hopes of broadening our horizons and drawing attention to this gap in bioprospecting studies for other niche environments as well. The review also aims to emphasize the importance of routinely incorporating functional genomics into environmental metagenomics analyses in order to improve solutions that drive rapid industrial biocatalysis developments from derived outputs with the aim of achieving potential benefits in energy-use reduction and environmental considerations for current and future applications.

Characterization of Host-Specific Genes from Pine- and Grass-Associated Species of the Fusarium fujikuroi Species Complex

The Fusarium fujikuroi species complex (FFSC) includes socioeconomically important pathogens that cause disease for numerous crops and synthesize a variety of secondary metabolites that can contaminate feedstocks and food. Here, we used comparative genomics to elucidate processes underlying the ability of pine-associated and grass-associated FFSC species to colonize tissues of their respective plant hosts. We characterized the identity, possible functions, evolutionary origins, and chromosomal positions of the host-range-associated genes encoded by the two groups of fungi. The 72 and 47 genes identified as unique to the respective genome groups were potentially involved in diverse processes, ranging from transcription, regulation, and substrate transport through to virulence/pathogenicity. Most genes arose early during the evolution of Fusarium/FFSC and were only subsequently retained in some lineages, while some had origins outside Fusarium. Although differences in the densities of these genes were especially noticeable on the conditionally dispensable chromosome of F. temperatum (representing the grass-associates) and F. circinatum (representing the pine-associates), the host-range-associated genes tended to be located towards the subtelomeric regions of chromosomes. Taken together, these results demonstrate that multiple mechanisms drive the emergence of genes in the grass- and pine-associated FFSC taxa examined. It also highlighted the diversity of the molecular processes potentially underlying niche-specificity in these and other Fusarium species.

Selenium Metabolism and Selenoproteins in Prokaryotes: A Bioinformatics Perspective

Selenium (Se) is an important trace element that mainly occurs in the form of selenocysteine in selected proteins. In prokaryotes, Se is also required for the synthesis of selenouridine and Se-containing cofactor. A large number of selenoprotein families have been identified in diverse prokaryotic organisms, most of which are thought to be involved in various redox reactions. In the last decade or two, computational prediction of selenoprotein genes and comparative genomics of Se metabolic pathways and selenoproteomes have arisen, providing new insights into the metabolism and function of Se and their evolutionary trends in bacteria and archaea. This review aims to offer an overview of recent advances in bioinformatics analysis of Se utilization in prokaryotes. We describe current computational strategies for the identification of selenoprotein genes and generate the most comprehensive list of prokaryotic selenoproteins reported to date. Furthermore, we highlight the latest research progress in comparative genomics and metagenomics of Se utilization in prokaryotes, which demonstrates the divergent and dynamic evolutionary patterns of different Se metabolic pathways, selenoprotein families, and selenoproteomes in sequenced organisms and environmental samples. Overall, bioinformatics analyses of Se utilization, function, and evolution may contribute to a systematic understanding of how this micronutrient is used in nature.

Genomic Fishing and Data Processing for Molecular Evolution Research

Molecular evolution analyses, such as detection of adaptive/purifying selection or ancestral protein reconstruction, typically require three inputs for a target gene (or gene family) in a particular group of organisms: sequence alignment, model of evolution, and phylogenetic tree. While modern advances in high-throughput sequencing techniques have led to rapid accumulation of genomic-scale data in public repositories and databases, mining such vast amount of information often remains a challenging enterprise. Here, we describe a comprehensive, versatile workflow aimed at the preparation of genome-extracted datasets readily available for molecular evolution research. The workflow involves: (1) fishing (searching and capturing) specific gene sequences of interest from taxonomically diverse genomic data available in databases at variable levels of annotation, (2) processing and depuration of retrieved sequences, (3) production of a multiple sequence alignment, (4) selection of best-fit model of evolution, and (5) solid reconstruction of a phylogenetic tree.

From Sequences to Enzymes: Comparative Genomics to Study Evolutionarily Conserved Protein Functions in Marine Microbes

Comparative genomics is a research field that allows comparison between genomes of different life forms providing information on the organization of the compared genomes, both in terms of structure and encoded functions. Moreover, this approach provides a powerful tool to study and understand the evolutionary changes and adaptation among organisms. Comparative genomics can be used to compare phylogenetically close marine organisms showing different vital strategies and lifestyles and obtain information regarding specific adaptations and/or their evolutionary history. Here we report a basic comparative genomics protocol to extrapolate evolutionary information about a protein of interest conserved across diverse marine microbes. The outlined approach can be used in a number of different settings and might help to gain new insights into the evolution and adaptation of marine microorganisms.

Comparative Genomics and Evolution of Avian Specialized Traits

Genomic data are important for understanding the origin and evolution of traits. Under the context of rapidly developing of sequencing technologies and more widely available genome sequences, researchers are able to study evolutionary mechanisms of traits via comparative genomic methods. Compared with other vertebrates, bird genomes are relatively small and exhibit conserved synteny with few repetitive elements, which makes them suitable for evolutionary studies. Increasing genomic progress has been reported on the evolution of powered flight, body size variation, beak morphology, plumage colouration, high-elevation colonization, migration, and vocalization. By summarizing previous studies, we demonstrate the genetic bases of trait evolution, highlighting the roles of small-scale sequence variation, genomic structural variation, and changes in gene interaction networks. We suggest that future studies should focus on improving the quality of reference genomes, exploring the evolution of regulatory elements and networks, and combining genomic data with morphological, ecological, behavioural, and developmental biology data.

Phenotyping in the era of genomics: MaTrics—a digital character matrix to document mammalian phenotypic traits

A new and uniquely structured matrix of mammalian phenotypes, MaTrics (Mammalian Traits for Comparative Genomics) in a digital form is presented. By focussing on mammalian species for which genome assemblies are available, MaTrics provides an interface between mammalogy and comparative genomics.

MaTrics was developed within a project aimed to find genetic causes of phenotypic traits of mammals using Forward Genomics. This approach requires genomes and comprehensive and recorded information on homologous phenotypes that are coded as discrete categories in a matrix. MaTrics is an evolving online resource providing information on phenotypic traits in numeric code; traits are coded either as absent/present or with several states as multistate. The state record for each species is linked to at least one reference (e.g., literature, photographs, histological sections, CT scans, or museum specimens) and so MaTrics contributes to digitalization of museum collections. Currently, MaTrics covers 147 mammalian species and includes 231 characters related to structure, morphology, physiology, ecology, and ethology and available in a machine actionable NEXUS-format*. Filling MaTrics revealed substantial knowledge gaps, highlighting the need for phenotyping efforts. Studies based on selected data from MaTrics and using Forward Genomics identified associations between genes and certain phenotypes ranging from lifestyles (e.g., aquatic) to dietary specializations (e.g., herbivory, carnivory). These findings motivate the expansion of phenotyping in MaTrics by filling research gaps and by adding taxa and traits. Only databases like MaTrics will provide machine actionable information on phenotypic traits, an important limitation to genomics. MaTrics is available within the data repository Morph·D·Base (www.morphdbase.de).

Oxytocin variation and brain region-specific gene expression in a domesticated avian species

The Bengalese finch was domesticated more than 250 years ago from the wild white-rumped munia (WRM). Similar to other domesticated species, Bengalese finches show a reduced fear response and have lower corticosterone levels, compared to WRMs. Bengalese finches and munias also have different song types. Since oxytocin (OT) has been found to be involved in stress coping and auditory processing, we tested whether the OT sequence and brain expression pattern and content differ in wild munias and domesticated Bengalese finches. We sequenced the OT from 10 wild munias and 11 Bengalese finches and identified intra-strain variability in both the untranslated and protein-coding regions of the sequence, with all the latter giving rise to synonymous mutations. Several of these changes fall in specific transcription factor-binding sites, and show either a conserved or a relaxed evolutionary trend in the avian lineage, and in vertebrates in general. Although in situ hybridization in several hypothalamic nuclei did not reveal significant differences in the number of cells expressing OT between the two strains, real-time quantitative PCR showed a significantly higher OT mRNA expression in the cerebrum of the Bengalese finches relative to munias, but a significantly lower expression in their diencephalon. Our study thus points to a brain region-specific pattern of neurochemical expression in domesticated and wild avian strains, which could be linked to domestication and the behavioral changes associated with it.

Genetic Diversity for Accelerating Microbial Adaptive Laboratory Evolution

Adaptive laboratory evolution (ALE) is a widely used and highly effective tool for improving microbial phenotypes and investigating the evolutionary roots of biological phenomena. Serving as the raw materials of evolution, mutations have been extensively utilized to increase the chances of engineering molecules or microbes with tailor-made functions. The generation of genetic diversity is therefore a core technology for accelerating ALE, and a high-quality mutant library is crucial to its success. Because of its importance, technologies for generating genetic diversity have undergone rapid development in recent years. Here, we review the existing techniques for the construction of mutant libraries, briefly introduce their mechanisms and applications, discuss ongoing and emerging efforts to apply engineering technologies in the construction of mutant libraries, and suggest future perspectives for library construction.

Selenium in Human Health and Gut Microflora: Bioavailability of Selenocompounds and Relationship With Diseases

This review covers current knowledge of selenium in the dietary intake, its bioavailability, metabolism, functions, biomarkers, supplementation and toxicity, as well as its relationship with diseases and gut microbiota specifically on the symbiotic relationship between gut microflora and selenium status. Selenium is essential for the maintenance of the immune system, conversion of thyroid hormones, protection against the harmful action of heavy metals and xenobiotics as well as for the reduction of the risk of chronic diseases. Selenium is able to balance the microbial flora avoiding health damage associated with dysbiosis. Experimental studies have shown that inorganic and organic selenocompounds are metabolized to selenomethionine and incorporated by bacteria from the gut microflora, therefore highlighting their role in improving the bioavailability of selenocompounds. Dietary selenium can affect the gut microbial colonization, which in turn influences the host's selenium status and expression of selenoproteoma. Selenium deficiency may result in a phenotype of gut microbiota that is more susceptible to cancer, thyroid dysfunctions, inflammatory bowel disease, and cardiovascular disorders. Although the host and gut microbiota benefit each other from their symbiotic relationship, they may become competitors if the supply of micronutrients is limited. Intestinal bacteria can remove selenium from the host resulting in two to three times lower levels of host's selenoproteins under selenium-limiting conditions. There are still gaps in whether these consequences are unfavorable to humans and animals or whether the daily intake of selenium is also adapted to meet the needs of the bacteria.

AniProtDB: A Collection of Consistently Generated Metazoan Proteomes for Comparative Genomics Studies

To address the void in the availability of high-quality proteomic data traversing the animal tree, we have implemented a pipeline for generating de novo assemblies based on publicly available data from the NCBI Sequence Read Archive, yielding a comprehensive collection of proteomes from 100 species spanning 21 animal phyla. We have also created the Animal Proteome Database (AniProtDB), a resource providing open access to this collection of high-quality metazoan proteomes, along with information on predicted proteins and protein domains for each taxonomic classification and the ability to perform sequence similarity searches against all proteomes generated using this pipeline. This solution vastly increases the utility of these data by removing the barrier to access for research groups who do not have the expertise or resources to generate these data themselves and enables the use of data from non-traditional research organisms that have the potential to address key questions in biomedicine.

Application of the omics sciences to the study of Naegleria fowleri, Acanthamoeba spp., and Balamuthia mandrillaris: current status and future projections

In this review, we focus on the sequenced genomes of the pathogens Naegleria fowleri, Acanthamoeba spp. and Balamuthia mandrillaris, and the remarkable discoveries regarding the pathogenicity and genetic information of these organisms, using techniques related to the various omics branches like genomics, transcriptomics, and proteomics. Currently, novel data produced through comparative genomics analyses and both differential gene and protein expression in these free-living amoebas have allowed for breakthroughs to identify genes unique to N. fowleri, genes with active transcriptional activity, and their differential expression in conditions of modified virulence. Furthermore, orthologous genes of the various nuclear genomes within the Naegleria and Acanthamoeba genera have been clustered. The proteome of B. mandrillaris has been reconstructed through transcriptome data, and its mitochondrial genome structure has been thoroughly described with a unique characteristic that has come to light: a type I intron with the capacity of interrupting genes through its self-splicing ribozymes activity. With the integration of data derived from the diverse omic sciences, there is a potential approximation that reflects the molecular complexity required for the identification of virulence factors, as well as crucial information regarding the comprehension of the molecular mechanisms with which these interact. Altogether, these breakthroughs could contribute to radical advances in both the fields of therapy design and medical diagnosis in the foreseeable future.

Genomic diversity of Mycobacterium avium subsp. paratuberculosis: pangenomic approach for highlighting unique genomic features with newly constructed complete genomes

Mycobacterium avium subsp. paratuberculosis (MAP) is a causative agent of Johne’s disease, which is a chronic granulomatous enteropathy in ruminants. Determining the genetic diversity of MAP is necessary to understand the epidemiology and biology of MAP, as well as establishing disease control strategies. In the present study, whole genome-based alignment and comparative analysis were performed using 40 publicly available MAP genomes, including newly sequenced Korean isolates. First, whole genome-based alignment was employed to identify new genomic structures in MAP genomes. Second, the genomic diversity of the MAP population was described by pangenome analysis. A phylogenetic tree based on the core genome and pangenome showed that the MAP was differentiated into two major types (C- and S-type), which was in keeping with the findings of previous studies. However, B-type strains were discriminated from C-type strains. Finally, functional analysis of the pangenome was performed using three virulence factor databases (i.e., PATRIC, VFDB, and Victors) to predict the phenotypic diversity of MAP in terms of pathogenicity. Based on the results of the pangenome analysis, we developed a real-time PCR technique to distinguish among S-, B- and C-type strains. In conclusion, the results of our study suggest that the phenotypic differences between MAP strains can be explained by their genetic polymorphisms. These results may help to elucidate the diversity of MAP, extending from genomic features to phenotypic traits.

Elevated diversity in loci linked to facial morphology is consistent with the hypothesis that individual facial recognition is important across hominoids

OBJECTIVES

The ability to use visual signals to identify individuals is an important feature of primate social groups, including humans. Sheehan and Nachman (2014) showed that loci linked to facial morphology had elevated levels of diversity and interpreted this as evidence that the human face is under frequency-dependent selection to enhance individual recognition (Nature Communications 5). In our study, we tested whether this pattern is found in non-human ape species, to help understand whether individual recognition might also play a role in species other than humans.

MATERIALS AND METHODS

We examined levels of genetic diversity in an available population genomic dataset of humans, chimpanzees, bonobos, gorillas, and orangutans for three sets of loci, (1) loci linked to facial morphology, (2) loci linked to height, and (3) neutrally evolving regions. We tested whether loci linked to facial morphology were more variable than loci linked to height or neutrally evolving loci in each of these species.

RESULTS

We found significantly elevated diversity in loci linked to facial morphology in chimpanzees, gorillas, and Sumatran and Bornean orangutans.

DISCUSSION

Our findings closely parallel those of Sheehan and Nachman and are consistent with the idea that selection for facial diversity and individual recognition has not only shaped the evolution of the human face, but it has similarly shaped the evolution of most of our closest primate relatives. We also discuss alternative hypotheses for this pattern.

Bioinformatics as a Tool for the Structural and Evolutionary Analysis of Proteins

This chapter deals with the topic of bioinformatics, computational, mathematics, and statistics tools applied to biology, essential for the analysis and characterization of biological molecules, in particular proteins, which play an important role in all cellular and evolutionary processes of the organisms. In recent decades, with the next generation sequencing technologies and bioinformatics, it has facilitated the collection and analysis of a large amount of genomic, transcriptomic, proteomic, and metabolomic data from different organisms that have allowed predictions on the regulation of expression, transcription, translation, structure, and mechanisms of action of proteins as well as homology, mutations, and evolutionary processes that generate structural and functional changes over time. Although the information in the databases is greater every day, all bioinformatics tools continue to be constantly modified to improve performance that leads to more accurate predictions regarding protein functionality, which is why bioinformatics research remains a great challenge.

MapGL: inferring evolutionary gain and loss of short genomic sequence features by phylogenetic maximum parsimony

Background

Comparative genomics studies are growing in number partly because of their unique ability to provide insight into shared and divergent biology between species. Of particular interest is the use of phylogenetic methods to infer the evolutionary history of cis-regulatory sequence features, which contribute strongly to phenotypic divergence and are frequently gained and lost in eutherian genomes. Understanding the mechanisms by which cis-regulatory element turnover generate emergent phenotypes is crucial to our understanding of adaptive evolution. Ancestral reconstruction methods can place species-specific cis-regulatory features in their evolutionary context, thus increasing our understanding of the process of regulatory sequence turnover. However, applying these methods to gain and loss of cis-regulatory features historically required complex workflows, preventing widespread adoption by the broad scientific community.

Results

MapGL simplifies phylogenetic inference of the evolutionary history of short genomic sequence features by combining the necessary steps into a single piece of software with a simple set of inputs and outputs. We show that MapGL can reliably disambiguate the mechanisms underlying differential regulatory sequence content across a broad range of phylogenetic topologies and evolutionary distances. Thus, MapGL provides the necessary context to evaluate how genomic sequence gain and loss contribute to species-specific divergence.

Conclusions

MapGL makes phylogenetic inference of species-specific sequence gain and loss easy for both expert and non-expert users, making it a powerful tool for gaining novel insights into genome evolution.

Characterisation of sperm piRNAs and their correlation with semen quality traits in swine

Piwi-interacting RNAs (piRNAs) are a class of non-coding RNAs that are essential in the transcriptional silencing of transposable elements and warrant genome stability in the mammalian germline. In this study, we have identified piRNAs in porcine sperm using male germline and zygote datasets from human, mice, cow and pig, and evaluated the relation between their abundances and sperm quality traits. In our analysis, we identified 283 382 piRNAs, 1355 of which correlated with P ≤ 0.01 to at least one semen quality trait. Fifty-seven percent of the correlated piRNAs mapped less than 50 kb apart from any other piRNA in the pig genome. Furthermore, piRNA location was significantly enriched near long interspersed nuclear elements. Moreover, some of the significant piRNAs mapped within or close to genes relevant for fertility or spermatogenesis such as CSNK1G2 and PSMF1.

Chromatin Remodelers in the 3D Nuclear Compartment

Chromatin remodeling complexes (CRCs) use ATP hydrolysis to maintain correct expression profiles, chromatin stability, and inherited epigenetic states. More than 20 CRCs have been described to date, which encompass four large families defined by their ATPase subunits. These complexes and their subunits are conserved from yeast to humans through evolution. Their activities depend on their catalytic subunits which through ATP hydrolysis provide the energy necessary to fulfill cellular functions such as gene transcription, DNA repair, and transposon silencing. These activities take place at the first levels of chromatin compaction, and CRCs have been recognized as essential elements of chromatin dynamics. Recent studies have demonstrated an important role for these complexes in the maintenance of higher order chromatin structure. In this review, we present an overview of the organization of the genome within the cell nucleus, the different levels of chromatin compaction, and importance of the architectural proteins, and discuss the role of CRCs and how their functions contribute to the dynamics of the 3D genome organization.

Empowering a Methanol-Dependent Escherichia coli via Adaptive Evolution Using a High-Throughput Microbial Microdroplet Culture System

Recently, a methanol-essential Escherichia coli was constructed; this strain is highly dependent on a supply of gluconate as a co-substrate for growth. Adaptive laboratory evolution is commonly applied to obtain mutants with specific phenotypes under certain selected pressure. However, conventional adaptive evolution approaches are not only laborious and time consuming, but they also come with lower throughput and inefficiency. In order to empower the aforementioned E. coli with reduced gluconate usage and enhanced growth rate, an irrational strategy based on a microbial microdroplet culture (MMC) platform was developed in this study. Given the automatic high-throughput selection of the MMC, a three-stage regime of an adaptive evolution experiment via gradually decreasing the availability of gluconate during the cultivation was performed for 50 days continuously in order to obtain the mutations. Finally, a candidate mutant was obtained with a 3-fold faster growth rate, a 43% shorter lag phase, and 40% less gluconate usage compared with the starting strain. Moreover, the gene mutations of gntU, idnT, edd, and pckA were identified by analyzing the whole-genome sequencing of mutants, which are strongly associated with the efficiency of gluconate uptake and cell growth. In conclusion, we have successfully demonstrated the feasibility of using MMC platform to empower the target strain with specific requirements in a manner of labor, time efficiency, and directed evolution.

Design of a Multiepitope-Based Peptide Vaccine against the E Protein of Human COVID-19: An Immunoinformatics Approach

Background

A new endemic disease has spread across Wuhan City, China, in December 2019. Within few weeks, the World Health Organization (WHO) announced a novel coronavirus designated as coronavirus disease 2019 (COVID-19). In late January 2020, WHO declared the outbreak of a “public-health emergency of international concern” due to the rapid and increasing spread of the disease worldwide. Currently, there is no vaccine or approved treatment for this emerging infection; thus, the objective of this study is to design a multiepitope peptide vaccine against COVID-19 using an immunoinformatics approach.

Method

Several techniques facilitating the combination of the immunoinformatics approach and comparative genomic approach were used in order to determine the potential peptides for designing the T-cell epitope-based peptide vaccine using the envelope protein of 2019-nCoV as a target.

Results

Extensive mutations, insertion, and deletion were discovered with comparative sequencing in the COVID-19 strain. Additionally, ten peptides binding to MHC class I and MHC class II were found to be promising candidates for vaccine design with adequate world population coverage of 88.5% and 99.99%, respectively.

Conclusion

The T-cell epitope-based peptide vaccine was designed for COVID-19 using the envelope protein as an immunogenic target. Nevertheless, the proposed vaccine rapidly needs to be validated clinically in order to ensure its safety and immunogenic profile to help stop this epidemic before it leads to devastating global outbreaks.

Intelligent Design of 14-3-3 Docking Proteins Utilizing Synthetic Evolution Artificial Intelligence (SYN-AI)

The ability to write DNA code from scratch will allow for the discovery of new and interesting chemistries as well as allowing the rewiring of cell signal pathways. Herein, we have utilized synthetic evolution artificial intelligence (SYN-AI) to intelligently design a set of 14-3-3 docking genes. SYN-AI engineers synthetic genes utilizing a parental gene as an evolution template. Wherein, evolution is fast-forwarded by transforming template gene sequences to DNA secondary and tertiary codes based upon gene hierarchical structural levels. The DNA secondary code allows identification of genomic building blocks across an orthologous sequence space comprising multiple genomes. Where, the DNA tertiary code allows engineering of supersecondary structures. SYN-AI constructed a library of 10 million genes that was reduced to three structurally functional 14-3-3 docking genes by applying natural selection protocols. Synthetic protein identity was verified utilizing Clustal Omega sequence alignments and Phylogeny.fr phylogenetic analysis. Wherein, we were able to confirm the three-dimensional structure utilizing I-TASSER and protein-ligand interactions utilizing COACH and Cofactor. The conservation of allosteric communications was confirmed utilizing elastic and anisotropic network models. Whereby, we utilized elNemo and ANM2.1 to confirm conservation of the 14-3-3 ζ amphipathic groove. Notably, to the best of our knowledge, we report the first 14-3-3 docking genes to be written from scratch.

IRAK3 modulates downstream innate immune signalling through its guanylate cyclase activity

Interleukin-1 receptor associated kinase 3 (IRAK3) is a cytoplasmic homeostatic mediator of inflammatory responses and is potentially useful as a prognostic marker in inflammation. IRAK3 inhibits signalling cascades downstream of myddosome complexes associated with toll like receptors. IRAK3 contains a death domain that interacts with other IRAK family members, a pseudokinase domain and a C-terminus domain involved with tumour necrosis factor receptor associated factor 6 (TRAF6). Previous bioinformatic studies revealed that IRAK3 contained a guanylate cyclase centre in its pseudokinase domain but its role in IRAK3 action is unresolved. We demonstrate that wildtype IRAK3 is capable of producing cGMP. Furthermore, we show that a specific point mutation in the guanylate cyclase centre reduced cGMP production. Cells containing toll like receptor 4 and a nuclear factor kappa-light-chain-enhancer of activated B cells (NFĸB) reporter system were transfected with IRAK3 or mutant IRAK3 proteins. Cell-permeable cGMP treatment of untransfected control cells suppresses downstream signalling through modulation of the NFĸB in the presence of lipopolysaccharides. Cells transfected with wildtype IRAK3 also suppress lipopolysaccharide induced NFĸB activity in the absence of exogenous cGMP. Lipopolysaccharide induced NFĸB activity was not suppressed in cells transfected with the IRAK3 mutant with reduced cGMP-generating capacity. Whereas in the presence of exogenously applied cell-permeable cGMP the IRAK3 mutant was able to retain its function by suppressing lipopolysaccharide induced NFĸB activity. Furthermore, increasing the amount of membrane permeable cGMP did not affect IRAK3's ability to reduce NFĸB activity. These results suggest that cGMP generated by IRAK3 may be involved in regulatory function of the protein where the presence of cGMP may selectively affect downstream signalling pathway(s) by modulating binding and/or activity of nearby proteins that interact in the inflammatory signalling cascade.

Understanding Adaptation, Coevolution, Host Specialization, and Mating System in Castrating Anther-Smut Fungi by Combining Population and Comparative Genomics

Anther-smut fungi provide a powerful system to study host-pathogen specialization and coevolution, with hundreds of Microbotryum species specialized on diverse Caryophyllaceae plants, castrating their hosts through manipulation of the hosts' reproductive organs to facilitate disease transmission. Microbotryum fungi have exceptional genomic characteristics, including dimorphic mating-type chromosomes, that make this genus anexcellent model for studying the evolution of mating systems and their influence on population genetics structure and adaptive potential. Important insights into adaptation, coevolution, host specialization, and mating system evolution have been gained using anther-smut fungi, with new insights made possible by the recent advent of genomic approaches. We illustrate with Microbotryum case studies how using a combination of comparative genomics, population genomics, and transcriptomics approaches enables the integration of different evolutionary perspectives across different timescales. We also highlight current challenges and suggest future studies that will contribute to advancing our understanding of the mechanisms underlying adaptive processes in populations of fungal pathogens.

Genomic analysis of the aggressive tree pathogen Ceratocystis albifundus

The overall goal of this study was to determine whether the genome of an important plant pathogen in Africa, Ceratocystis albifundus, is structured into subgenomic compartments, and if so, to establish how these compartments are distributed across the genome. For this purpose, the publicly available genome of C. albifundus was complemented with the genome sequences for four additional isolates using the Illumina HiSeq platform. In addition, a reference genome for one of the individuals was assembled using both PacBio and Illumina HiSeq technologies. Our results showed a high degree of synteny between the five genomes, although several regions lacked detectable long-range synteny. These regions were associated with the presence of accessory genes, lower genetic similarity, variation in read-map depth, as well as transposable elements and genes associated with host-pathogen interactions (e.g. effectors and CAZymes). Such patterns are regarded as hallmarks of accelerated evolution, particularly of accessory subgenomic compartments in fungal pathogens. Our findings thus showed that the genome of C. albifundus is made-up of core and accessory subgenomic compartments, which is an important step towards characterizing its pangenome. This study also highlights the value of comparative genomics for understanding mechanisms that may underly and influence the biology and evolution of pathogens.

Comparative genomics and metagenomics of the metallomes

Recent achievements and advances in comparative genomic and metagenomic analyses of trace metals were reviewed and discussed.

Genomic Structural Variations Within Five Continental Populations of Drosophila melanogaster

Chromosomal structural variations (SV) including insertions, deletions, inversions, and translocations occur within the genome and can have a significant effect on organismal phenotype. Some of these effects are caused by structural variations containing genes. Large structural variations represent a significant amount of the genetic diversity within a population. We used a global sampling of Drosophila melanogaster (Ithaca, Zimbabwe, Beijing, Tasmania, and Netherlands) to represent diverse populations within the species. We used long-read sequencing and optical mapping technologies to identify SVs in these genomes. Among the five lines examined, we found an average of 2,928 structural variants within these genomes. These structural variations varied greatly in size and location, included many exonic regions, and could impact adaptation and genomic evolution.

From zero to hero: the past, present and future of grain amaranth breeding

Grain amaranth is an underutilized crop with high nutritional quality from the Americas. Emerging genomic and biotechnological tools are becoming available that allow the integration of novel breeding techniques for rapid improvement of amaranth and other underutilized crops. Out of thousands of edible plants, only three cereals-maize, wheat and rice-are the major food sources for a majority of people worldwide. While these crops provide high amounts of calories, they are low in protein and other essential nutrients. The dependence on only few crops, with often narrow genetic basis, leads to a high vulnerability of modern cropping systems to the predicted climate change and accompanying weather extremes. Broadening our food sources through the integration of so-called orphan crops can help to mitigate the effects of environmental change and improve qualitative food security. Thousands of traditional crops are known, but have received little attention in the last century and breeding efforts were limited. Amaranth is such an underutilized pseudocereal that is of particular interest because of its balanced amino acid and micronutrient profiles. Additionally, the C₄ photosynthetic pathway and ability to withstand environmental stress make the crop a suitable choice for future agricultural systems. Despite the potential of amaranth, efforts of genetic improvement lag considerably behind those of major crops. The progress in novel breeding methods and molecular techniques developed in model plants and major crops allow a rapid improvement of underutilized crops. Here, we review the history of amaranth and recent advances in genomic tools and give a concrete perspective how novel breeding techniques can be implemented into breeding programs. Our perspectives are transferable to many underutilized crops. The implementation of these could improve the nutritional quality and climate resilience of future cropping systems.

Homology Detection Using Multilayer Maximum Clustering Coefficient

Homologous sequences are widely used to understand the functions of certain genes or proteins. However, there is no consensus to solve the automatic assignment of functions to protein problem and many algorithms have different ways of identifying homologous clusters in a given set of sequences. In this article, we present an algorithm to deal with specific sets, the set of coding sequences obtained from phylogenetically close genomes (of the same species, genus, or family). When modeled as a graph, these sets have their own characteristics: they form more homogeneous and denser clusters. To solve this problem, our algorithm makes use of the clustering coefficient, which maximization can lead to the expected results from the biological point of view. In addition, we also present an algorithm for the identification of sequence domains based on graph topology. We also compared our results with those of the TribeMCL tool, a well-established algorithm of the area.

Pan4Draft: A Computational Tool to Improve the Accuracy of Pan-Genomic Analysis Using Draft Genomes

High-throughput sequencing technologies are a milestone in molecular biology for facilitating great advances in genomics by enabling the deposit of large volumes of biological data to public databases. The availability of such data has made possible the comparative genomic analysis through pipelines, using the entire gene repertoire of genomes. However, a large number of unfinished genomes exist in public databases; their number is approximately 16-fold higher than the number of complete genomes, which creates bias during comparative analyses. Therefore, the present work proposes a new tool called Pan4Drafts, an automated pipeline for pan-genomic analysis of draft prokaryotic genomes to maximize the representation and accuracy of the gene repertoire of unfinished genomes by using reads from sequencing data. Pan4Draft allows to perform comparative analyses using different methodologies such as combining complete and draft genomes, using only draft genomes or only complete genomes. Pan4Draft is available at http://www.computationalbiology.ufpa.br/pan4drafts and the test dataset is available at https://sourceforge.net/projects/pan4drafts .

Incorporating comparative genomics into the design-test-learn cycle of microbial strain engineering

Engineering microbes with new properties is an important goal in industrial engineering, to establish biological factories for production of biofuels, commodity chemicals and pharmaceutics. But engineering microbes to produce new compounds with high yield remains a major challenge toward economically viable production. Incorporating several modern approaches, including synthetic and systems biology, metabolic modeling and regulatory rewiring, has proven to significantly advance industrial strain engineering. This review highlights how comparative genomics can also facilitate strain engineering, by identifying novel genes and pathways, regulatory mechanisms and genetic background effects for engineering. We discuss how incorporating comparative genomics into the design-test-learn cycle of strain engineering can provide novel information that complements other engineering strategies.

A Guide to the PLAZA 3.0 Plant Comparative Genomic Database

PLAZA 3.0 is an online resource for comparative genomics and offers a versatile platform to study gene functions and gene families or to analyze genome organization and evolution in the green plant lineage. Starting from genome sequence information for over 35 plant species, precomputed comparative genomic data sets cover homologous gene families, multiple sequence alignments, phylogenetic trees, and genomic colinearity information within and between species. Complementary functional data sets, a Workbench, and interactive visualization tools are available through a user-friendly web interface, making PLAZA an excellent starting point to translate sequence or omics data sets into biological knowledge. PLAZA is available at http://bioinformatics.psb.ugent.be/plaza/ .

New insights into the factors affecting synonymous codon usage in human infecting Plasmodium species

Codon usage bias is due to the non-random usage of synonymous codons for coding amino acids. The synonymous sites are under weak selection, and codon usage bias is maintained by the equilibrium in mutational bias, genetic drift and selection pressure. The differential codon usage choices are also relevant to human infecting Plasmodium species. Recently, P. knowlesi switches its natural host, long-tailed macaques, and starts infecting humans. This review focuses on the comparative analysis of codon usage choices among human infecting P. falciparum and P. vivax along with P. knowlesi species taking their coding sequence data. The variation in GC content, amino acid frequencies, effective number of codons and other factors plays a crucial role in determining synonymous codon choices. Within species codon choices are more similar for P. vivax and P. knowlesi in comparison with P. falciparum species. This study suggests that synonymous codon choice modulates the gene expression level, mRNA stability, ribosome speed, protein folding, translation efficiency and its accuracy in Plasmodium species, and provides a valuable information regarding the codon usage pattern to facilitate gene cloning as well as expression and transfection studies for malaria causing species.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics

BACKGROUND

With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies.

RESULTS

We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements.

CONCLUSIONS

Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

Molecular Cloning of Porcine SUN5 Gene and Association between a SNP with Litter Size Trait

SUN domain-containing protein 5 (SUN5) is an important reproduction related gene. In this study, we cloned the full-length coding sequence of porcine SUN5 gene through RT-PCR. Sequence analysis of this gene revealed that the pig SUN5 gene encodes a protein of 383 amino acids that has high homology with the SUN5 protein of eight species: wild Bactrian camel (95%), alpaca (95%), Yangtze River dolphin (94%), sperm whale (94%), sheep (93%), black flying fox (93%), goat (92%), and horse (91%). This gene is structured into 13 exons and 12 introns as revealed by computer-assisted analysis. The prediction of transmembrane helices showed that pig SUN5 protein might be a transmembrane protein. PCR-Taq I-RFLP was established to detect the GU475008:c.138 G>A substitution of porcine SUN5 gene coding sequence and eight pig breeds displayed obvious genotype and allele frequency differences at this mutation locus. Association of this SNP with litter size traits was assessed in Large White (n = 200) and Landrace (n = 200) pig populations, and the results demonstrated that this polymorphic locus was significantly associated with the litter size of all parities in Large White and Landrace sows (P < 0.05). Therefore, the SUN5 gene could be a useful candidate gene for increasing the litter size in pigs.

Differential regulation of mouse and human Mu opioid receptor gene depends on the single stranded DNA structure of its promoter and α-complex protein 1

The Mu opioid receptor (MOR) mediates various functions of opioid-induced analgesia, euphoria and respiratory depression, and is a major target of opioid analgesics. Understanding of MOR gene expression among species is important for understanding its analgesic function in humans. In the current study, the polypyrimidine/polypurine (PPy/u) region, a key element of MOR gene expression, was compared in humans and mice. The mouse PPy/u element is highly homologous to its human element (84%), and the mouse MOR (mMOR) reporter drives luciferase activity 35-fold more effectively than the human MOR (hMOR) reporter. The structural study of reporter plasmids using S1 nuclease indicates that the mouse PPy/u element has a particular conformational structure, namely a single-stranded DNA (ssDNA) region that promotes strong promoter activity. DNA electrophoretic mobility shift assays demonstrated that the recombinant α-complex protein 1 (α-CP1) is capable of binding to a single-stranded mouse PPy/u sequence. Furthermore, plasmid-expressing α-CP1 activated the expression of a luciferase reporter when cotransfected with a single-stranded (p336/306) construct. In addition, the α-CP1 gene induced the mMOR gene in mouse neuronal cells and did not induce the human neuronal MOR gene. The current study demonstrates that α-CP1 functions as a transcriptional activator in the mMOR gene, but does not function in the hMOR gene due to species-specific structural differences. The differences in human and mouse MOR gene expression are based on α-CP1 and the ssDNA structure of the MOR promoter. The MOR gene is species-specifically regulated, as the PPy/u element adopts a unique species-specific conformation and α-CP1 recruitment.

Phylogenomic analysis of gene co-expression networks reveals the evolution of functional modules

Molecular evolutionary studies correlate genomic and phylogenetic information with the emergence of new traits of organisms. These traits are, however, the consequence of dynamic gene networks composed of functional modules, which might not be captured by genomic analyses. Here, we established a method that combines large-scale genomic and phylogenetic data with gene co-expression networks to extensively study the evolutionary make-up of modules in the moss Physcomitrella patens, and in the angiosperms Arabidopsis thaliana and Oryza sativa (rice). We first show that younger genes are less annotated than older genes. By mapping genomic data onto the co-expression networks, we found that genes from the same evolutionary period tend to be connected, whereas old and young genes tend to be disconnected. Consequently, the analysis revealed modules that emerged at a specific time in plant evolution. To uncover the evolutionary relationships of the modules that are conserved across the plant kingdom, we added phylogenetic information that revealed duplication and speciation events on the module level. This combined analysis revealed an independent duplication of cell wall modules in bryophytes and angiosperms, suggesting a parallel evolution of cell wall pathways in land plants. We provide an online tool allowing plant researchers to perform these analyses at http://www.gene2function.de.

Orchestrating differential data access for translational research: a pilot implementation

Background

Translational researchers need robust IT solutions to access a range of data types, varying from public data sets to pseudonymised patient information with restricted access, provided on a case by case basis. The reason for this complication is that managing access policies to sensitive human data must consider issues of data confidentiality, identifiability, extent of consent, and data usage agreements. All these ethical, social and legal aspects must be incorporated into a differential management of restricted access to sensitive data.

Methods

In this paper we present a pilot system that uses several common open source software components in a novel combination to coordinate access to heterogeneous biomedical data repositories containing open data (open access) as well as sensitive data (restricted access) in the domain of biobanking and biosample research. Our approach is based on a digital identity federation and software to manage resource access entitlements.

Results

Open source software components were assembled and configured in such a way that they allow for different ways of restricted access according to the protection needs of the data. We have tested the resulting pilot infrastructure and assessed its performance, feasibility and reproducibility.

Conclusions

Common open source software components are sufficient to allow for the creation of a secure system for differential access to sensitive data. The implementation of this system is exemplary for researchers facing similar requirements for restricted access data. Here we report experience and lessons learnt of our pilot implementation, which may be useful for similar use cases. Furthermore, we discuss possible extensions for more complex scenarios.