Nanopore adaptive sampling of a metagenomic sample derived from a human monkeypox case.

In 2022, a series of human monkeypox cases in multiple countries led to the largest and most widespread outbreak outside the known endemic areas. Setup of proper genomic surveillance is of utmost importance to control such outbreaks. To this end, we performed Nanopore (PromethION P24) and Illumina (NextSeq. 2000) Whole Genome Sequencing (WGS) of a monkeypox sample. Adaptive sampling was applied for in silico depletion of the human host genome, allowing for the enrichment of low abundance viral DNA without a priori knowledge of sample composition. Nanopore sequencing allowed for high viral genome coverage, tracking of sample composition during sequencing, strain determination, and preliminary assessment of mutational pattern. In addition to that, only Nanopore data allowed us to resolve the entire monkeypox virus genome, with respect to two structural variants belonging to the genes OPG015 and OPG208. These SVs in important host range genes seem stable throughout the outbreak and are frequently misassembled and/or misannotated due to the prevalence of short read sequencing or short read first assembly. Ideally, standalone standard Illumina sequencing should not be used for Monkeypox WGS and de novo assembly, since it will obfuscate the structure of the genome, which has an impact on the quality and completeness of the genomes deposited in public databases and thus possibly on the ability to evaluate the complete genetic reason for the host range change of monkeypox in the current pandemic.

Knockout of Factor-Inhibiting HIF (Hif1an) in Colon Epithelium Attenuates Chronic Colitis but Does Not Reduce Colorectal Cancer in Mice.

Inflammatory bowel disease such as chronic colitis promotes colorectal cancer, which is a common cause of cancer mortality worldwide. Hypoxia is a characteristic of inflammation as well as of solid tumors and enforces a gene expression response controlled by hypoxia-inducible factors (HIFs). Once established, solid tumors are immunosuppressive to escape their abatement through immune cells. Although HIF activity is known to 1) promote cancer development and 2) drive tumor immune suppression through the secretion of adenosine, both prolyl hydroxylases and an asparaginyl hydroxylase termed factor-inhibiting HIF (FIH) negatively regulate HIF. Thus, FIH may act as a tumor suppressor in colorectal cancer development. In this study, we examined the role of colon epithelial FIH in a mouse model of colitis-induced colorectal cancer. We recapitulated colitis-associated colorectal cancer development in mice using the azoxymethane/dextran sodium sulfate model in Vil1-Cre/FIH+f/+f and wild-type siblings. Colon samples were analyzed regarding RNA and protein expression and histology. Vil1-Cre/FIH+f/+f mice showed a less severe colitis progress compared with FIH+f/+f animals and a lower number of infiltrating macrophages in the inflamed tissue. RNA sequencing analyses of colon tissue revealed a lower expression of genes associated with the immune response in Vil1-Cre/FIH+f/+f mice. However, tumor occurrence did not significantly differ between Vil1-Cre/FIH+f/+f and wild-type mice. Thus, FIH knockout in colon epithelial cells did not modulate colorectal cancer development but reduced the inflammatory response in chronic colitis.

Gene expression variability in long-term survivors of childhood cancer and cancer-free controls in response to ionizing irradiation.

BACKGROUND: Differential expression analysis is usually adjusted for variation. However, most studies that examined the expression variability (EV) have used computations affected by low expression levels and did not examine healthy tissue. This study aims to calculate and characterize an unbiased EV in primary fibroblasts of childhood cancer survivors and cancer-free controls (N0) in response to ionizing radiation. METHODS: Human skin fibroblasts of 52 donors with a first primary neoplasm in childhood (N1), 52 donors with at least one second primary neoplasm (N2 +), as well as 52 N0 were obtained from the KiKme case-control study and exposed to a high (2 Gray) and a low dose (0.05 Gray) of X-rays and sham- irradiation (0 Gray). Genes were then classified as hypo-, non-, or hyper-variable per donor group and radiation treatment, and then examined for over-represented functional signatures. RESULTS: We found 22 genes with considerable EV differences between donor groups, of which 11 genes were associated with response to ionizing radiation, stress, and DNA repair. The largest number of genes exclusive to one donor group and variability classification combination were all detected in N0: hypo-variable genes after 0 Gray (n = 49), 0.05 Gray (n = 41), and 2 Gray (n = 38), as well as hyper-variable genes after any dose (n = 43). While after 2 Gray positive regulation of cell cycle was hypo-variable in N0, (regulation of) fibroblast proliferation was over-represented in hyper-variable genes of N1 and N2+. In N2+, 30 genes were uniquely classified as hyper-variable after the low dose and were associated with the ERK1/ERK2 cascade. For N1, no exclusive gene sets with functions related to the radiation response were detected in our data. CONCLUSION: N2+ showed high degrees of variability in pathways for the cell fate decision after genotoxic insults that may lead to the transfer and multiplication of DNA-damage via proliferation, where apoptosis and removal of the damaged genome would have been appropriate. Such a deficiency could potentially lead to a higher vulnerability towards side effects of exposure to high doses of ionizing radiation, but following low-dose applications employed in diagnostics, as well.

Androglobin gene expression patterns and FOXJ1-dependent regulation indicate its functional association with ciliogenesis.

Androglobin (ADGB) represents the latest addition to the globin superfamily in metazoans. The chimeric protein comprises a calpain domain and a unique circularly permutated globin domain. ADGB expression levels are most abundant in mammalian testis, but its cell-type-specific expression, regulation, and function have remained unexplored. Analyzing bulk and single-cell mRNA-Seq data from mammalian tissues, we found that-in addition to the testes-ADGB is prominently expressed in the female reproductive tract, lungs, and brain, specifically being associated with cell types forming motile cilia. Correlation analysis suggested coregulation of ADGB with FOXJ1, a crucial transcription factor of ciliogenesis. Investigating the transcriptional regulation of the ADGB gene, we characterized its promoter using epigenomic datasets, exogenous promoter-dependent luciferase assays, and CRISPR/dCas9-VPR-mediated activation approaches. Reporter gene assays revealed that FOXJ1 indeed substantially enhanced luciferase activity driven by the ADGB promoter. ChIP assays confirmed binding of FOXJ1 to the endogenous ADGB promoter region. We dissected the minimal sequence required for FOXJ1-dependent regulation and fine mapped the FOXJ1 binding site to two evolutionarily conserved regions within the ADGB promoter. FOXJ1 overexpression significantly increased endogenous ADGB mRNA levels in HEK293 and MCF-7 cells. Similar results were observed upon RFX2 overexpression, another key transcription factor in ciliogenesis. The complex transcriptional regulation of the ADGB locus was illustrated by identifying a distal enhancer, responsible for synergistic regulation by RFX2 and FOXJ1. Finally, cell culture studies indicated an ADGB-dependent increase in the number of ciliated cells upon overexpression of the full-length protein, confirming a ciliogenesis-associated role of ADGB in mammals.

Identification of antiparasitic drug targets using a multi-omics workflow in the acanthocephalan model.

BACKGROUND: With the expansion of animal production, parasitic helminths are gaining increasing economic importance. However, application of several established deworming agents can harm treated hosts and environment due to their low specificity. Furthermore, the number of parasite strains showing resistance is growing, while hardly any new anthelminthics are being developed. Here, we present a bioinformatics workflow designed to reduce the time and cost in the development of new strategies against parasites. The workflow includes quantitative transcriptomics and proteomics, 3D structure modeling, binding site prediction, and virtual ligand screening. Its use is demonstrated for Acanthocephala (thorny-headed worms) which are an emerging pest in fish aquaculture. We included three acanthocephalans (Pomphorhynchus laevis, Neoechinorhynchus agilis, Neoechinorhynchus buttnerae) from four fish species (common barbel, European eel, thinlip mullet, tambaqui). RESULTS: The workflow led to eleven highly specific candidate targets in acanthocephalans. The candidate targets showed constant and elevated transcript abundances across definitive and accidental hosts, suggestive of constitutive expression and functional importance. Hence, the impairment of the corresponding proteins should enable specific and effective killing of acanthocephalans. Candidate targets were also highly abundant in the acanthocephalan body wall, through which these gutless parasites take up nutrients. Thus, the candidate targets are likely to be accessible to compounds that are orally administered to fish. Virtual ligand screening led to ten compounds, of which five appeared to be especially promising according to ADMET, GHS, and RO5 criteria: tadalafil, pranazepide, piketoprofen, heliomycin, and the nematicide derquantel. CONCLUSIONS: The combination of genomics, transcriptomics, and proteomics led to a broadly applicable procedure for the cost- and time-saving identification of candidate target proteins in parasites. The ligands predicted to bind can now be further evaluated for their suitability in the control of acanthocephalans. The workflow has been deposited at the Galaxy workflow server under the URL tinyurl.com/yx72rda7 .

Radiation-response in primary fibroblasts of long-term survivors of childhood cancer with and without second primary neoplasms: the KiKme study.

BACKGROUND: The etiology and most risk factors for a sporadic first primary neoplasm in childhood or subsequent second primary neoplasms are still unknown. One established causal factor for therapy-associated second primary neoplasms is the exposure to ionizing radiation during radiation therapy as a mainstay of cancer treatment. Second primary neoplasms occur in 8% of all cancer survivors within 30 years after the first diagnosis in Germany, but the underlying factors for intrinsic susceptibilities have not yet been clarified. Thus, the purpose of this nested case-control study was the investigation and comparison of gene expression and affected pathways in primary fibroblasts of childhood cancer survivors with a first primary neoplasm only or with at least one subsequent second primary neoplasm, and controls without neoplasms after exposure to a low and a high dose of ionizing radiation. METHODS: Primary fibroblasts were obtained from skin biopsies from 52 adult donors with a first primary neoplasm in childhood (N1), 52 with at least one additional primary neoplasm (N2+), as well as 52 without cancer (N0) from the KiKme study. Cultured fibroblasts were exposed to a high [2 Gray (Gy)] and a low dose (0.05 Gy) of X-rays. Messenger ribonucleic acid was extracted 4 h after exposure and Illumina-sequenced. Differentially expressed genes (DEGs) were computed using limma for R, selected at a false discovery rate level of 0.05, and further analyzed for pathway enrichment (right-tailed Fisher's Exact Test) and (in-) activation (z ≥|2|) using Ingenuity Pathway Analysis. RESULTS: After 0.05 Gy, least DEGs were found in N0 (n = 236), compared to N1 (n = 653) and N2+ (n = 694). The top DEGs with regard to the adjusted p-value were upregulated in fibroblasts across all donor groups (SESN1, MDM2, CDKN1A, TIGAR, BTG2, BLOC1S2, PPM1D, PHLDB3, FBXO22, AEN, TRIAP1, and POLH). Here, we observed activation of p53 Signaling in N0 and to a lesser extent in N1, but not in N2+. Only in N0, DNA (excision-) repair (involved genes: CDKN1A, PPM1D, and DDB2) was predicted to be a downstream function, while molecular networks in N2+ were associated with cancer, as well as injury and abnormalities (among others, downregulation of MSH6, CCNE2, and CHUK). After 2 Gy, the number of DEGs was similar in fibroblasts of all donor groups and genes with the highest absolute log2 fold-change were upregulated throughout (CDKN1A, TIGAR, HSPA4L, MDM2, BLOC1SD2, PPM1D, SESN1, BTG2, FBXO22, PCNA, and TRIAP1). Here, the p53 Signaling-Pathway was activated in fibroblasts of all donor groups. The Mitotic Roles of Polo Like Kinase-Pathway was inactivated in N1 and N2+. Molecular Mechanisms of Cancer were affected in fibroblasts of all donor groups. P53 was predicted to be an upstream regulator in fibroblasts of all donor groups and E2F1 in N1 and N2+. Results of the downstream analysis were senescence in N0 and N2+, transformation of cells in N0, and no significant effects in N1. Seven genes were differentially expressed in reaction to 2 Gy dependent on the donor group (LINC00601, COBLL1, SESN2, BIN3, TNFRSF10A, EEF1AKNMT, and BTG2). CONCLUSION: Our results show dose-dependent differences in the radiation response between N1/N2+ and N0. While mechanisms against genotoxic stress were activated to the same extent after a high dose in all groups, the radiation response was impaired after a low dose in N1/N2+, suggesting an increased risk for adverse effects including carcinogenesis, particularly in N2+.

Transcriptomes of Clusterin- and S100B-transfected neuronal cells elucidate protective mechanisms against hypoxia and oxidative stress in the hooded seal (Cystophora cristata) brain.

BACKGROUND: The hooded seal (Cystophora cristata) exhibits impressive diving skills and can tolerate extended durations of asphyxia, hypoxia and oxidative stress, without suffering from irreversible neuronal damage. Thus, when exposed to hypoxia in vitro, neurons of fresh cortical and hippocampal tissue from hooded seals maintained their membrane potential 4-5 times longer than neurons of mice. We aimed to identify the molecular mechanisms underlying the intrinsic neuronal hypoxia tolerance. Previous comparative transcriptomics of the visual cortex have revealed that S100B and clusterin (apolipoprotein J), two stress proteins that are involved in neurological disorders characterized by hypoxic conditions, have a remarkably high expression in hooded seals compared to ferrets. When overexpressed in murine neuronal cells (HN33), S100B and clusterin had neuroprotective effects when cells were exposed to hypoxia. However, their specific roles in hypoxia have remained largely unknown. METHODS: In order to shed light on potential molecular pathways or interaction partners, we exposed HN33 cells transfected with either S100B, soluble clusterin (sCLU) or nuclear clusterin (nCLU) to normoxia, hypoxia and oxidative stress for 24 h. We then determined cell viability and compared the transcriptomes of transfected cells to control cells. Potential pathways and upstream regulators were identified via Gene Ontology (GO) and Ingenuity Pathway Analysis (IPA). RESULTS: HN33 cells transfected with sCLU and S100B demonstrated improved glycolytic capacity and reduced aerobic respiration at normoxic conditions. Additionally, sCLU appeared to enhance pathways for cellular homeostasis to counteract stress-induced aggregation of proteins. S100B-transfected cells sustained lowered energy-intensive synaptic signaling. In response to hypoxia, hypoxia-inducible factor (HIF) pathways were considerably elevated in nCLU- and sCLU-transfected cells. In a previous study, S100B and sCLU decreased the amount of reactive oxygen species and lipid peroxidation in HN33 cells in response to oxidative stress, but in the present study, these functional effects were not mirrored in gene expression changes. CONCLUSIONS: sCLU and S100B overexpression increased neuronal survival by decreasing aerobic metabolism and synaptic signaling in advance to hypoxia and oxidative stress conditions, possibly to reduce energy expenditure and the build-up of deleterious reactive oxygen species (ROS). Thus, a high expression of CLU isoforms and S100B is likely beneficial during hypoxic conditions.

Pro-Apoptotic and Anti-Invasive Properties Underscore the Tumor-Suppressing Impact of Myoglobin on a Subset of Human Breast Cancer Cells.

The expression of myoglobin (MB), well known as the oxygen storage and transport protein of myocytes, is a novel hallmark of the luminal subtype in breast cancer patients and correlates with better prognosis. The mechanisms by which MB impacts mammary tumorigenesis are hitherto unclear. We aimed to unravel this role by using CRISPR/Cas9 technology to generate MB-deficient clones of MCF7 and SKBR3 breast cancer cell lines and subsequently characterize them by transcriptomics plus molecular and functional analyses. As main findings, loss of MB at normoxia upregulated the expression of cell cyclins and increased cell survival, while it prevented apoptosis in MCF7 cells. Additionally, MB-deficient cells were less sensitive to doxorubicin but not ionizing radiation. Under hypoxia, the loss of MB enhanced the partial epithelial to mesenchymal transition, thus, augmenting the migratory and invasive behavior of cells. Notably, in human invasive mammary ductal carcinoma tissues, MB and apoptotic marker levels were positively correlated. In addition, MB protein expression in invasive ductal carcinomas was associated with a positive prognostic value, independent of the known tumor suppressor p53. In conclusion, we provide multiple lines of evidence that endogenous MB in cancer cells by itself exerts novel tumor-suppressive roles through which it can reduce cancer malignancy.

Genomics and transcriptomics of epizoic Seisonidea (Rotifera, syn. Syndermata) reveal strain formation and gradual gene loss with growing ties to the host.

BACKGROUND: Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensu lato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea. RESULTS: Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria). CONCLUSION: Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits.

RainDrop: Rapid activation matrix computation for droplet-based single-cell RNA-seq reads.

BACKGROUND: Obtaining data from single-cell transcriptomic sequencing allows for the investigation of cell-specific gene expression patterns, which could not be addressed a few years ago. With the advancement of droplet-based protocols the number of studied cells continues to increase rapidly. This establishes the need for software tools for efficient processing of the produced large-scale datasets. We address this need by presenting RainDrop for fast gene-cell count matrix computation from single-cell RNA-seq data produced by 10x Genomics Chromium technology. RESULTS: RainDrop can process single-cell transcriptomic datasets consisting of 784 million reads sequenced from around 8.000 cells in less than 40 minutes on a standard workstation. It significantly outperforms the established Cell Ranger pipeline and the recently introduced Alevin tool in terms of runtime by a maximal (average) speedup of 30.4 (22.6) and 3.5 (2.4), respectively, while keeping high agreements of the generated results. CONCLUSIONS: RainDrop is a software tool for highly efficient processing of large-scale droplet-based single-cell RNA-seq datasets on standard workstations written in C++. It is available at https://gitlab.rlp.net/stnieble/raindrop .

A big data approach to metagenomics for all-food-sequencing.

BACKGROUND: All-Food-Sequencing (AFS) is an untargeted metagenomic sequencing method that allows for the detection and quantification of food ingredients including animals, plants, and microbiota. While this approach avoids some of the shortcomings of targeted PCR-based methods, it requires the comparison of sequence reads to large collections of reference genomes. The steadily increasing amount of available reference genomes establishes the need for efficient big data approaches. RESULTS: We introduce an alignment-free k-mer based method for detection and quantification of species composition in food and other complex biological matters. It is orders-of-magnitude faster than our previous alignment-based AFS pipeline. In comparison to the established tools CLARK, Kraken2, and Kraken2+Bracken it is superior in terms of false-positive rate and quantification accuracy. Furthermore, the usage of an efficient database partitioning scheme allows for the processing of massive collections of reference genomes with reduced memory requirements on a workstation (AFS-MetaCache) or on a Spark-based compute cluster (MetaCacheSpark). CONCLUSIONS: We present a fast yet accurate screening method for whole genome shotgun sequencing-based biosurveillance applications such as food testing. By relying on a big data approach it can scale efficiently towards large-scale collections of complex eukaryotic and bacterial reference genomes. AFS-MetaCache and MetaCacheSpark are suitable tools for broad-scale metagenomic screening applications. They are available at https://muellan.github.io/metacache/afs.html (C++ version for a workstation) and https://github.com/jmabuin/MetaCacheSpark (Spark version for big data clusters).

Comparison of time and dose dependent gene expression and affected pathways in primary human fibroblasts after exposure to ionizing radiation.

BACKGROUND: Exposure to ionizing radiation induces complex stress responses in cells, which can lead to adverse health effects such as cancer. Although a variety of studies investigated gene expression and affected pathways in human fibroblasts after exposure to ionizing radiation, the understanding of underlying mechanisms and biological effects is still incomplete due to different experimental settings and small sample sizes. Therefore, this study aims to identify the time point with the highest number of differentially expressed genes and corresponding pathways in primary human fibroblasts after irradiation at two preselected time points. METHODS: Fibroblasts from skin biopsies of 15 cell donors were exposed to a high (2Gy) and a low (0.05Gy) dose of X-rays. RNA was extracted and sequenced 2 h and 4 h after exposure. Differentially expressed genes with an adjusted p-value < 0.05 were flagged and used for pathway analyses including prediction of upstream and downstream effects. Principal component analyses were used to examine the effect of two different sequencing runs on quality metrics and variation in expression and alignment and for explorative analysis of the radiation dose and time point of analysis. RESULTS: More genes were differentially expressed 4 h after exposure to low and high doses of radiation than after 2 h. In experiments with high dose irradiation and RNA sequencing after 4 h, inactivation of the FAT10 cancer signaling pathway and activation of gluconeogenesis I, glycolysis I, and prostanoid biosynthesis was observed taking p-value (< 0.05) and (in) activating z-score (≥2.00 or ≤ - 2.00) into account. Two hours after high dose irradiation, inactivation of small cell lung cancer signaling was observed. For low dose irradiation experiments, we did not detect any significant (p < 0.05 and z-score ≥ 2.00 or ≤ - 2.00) activated or inactivated pathways for both time points. CONCLUSIONS: Compared to 2 h after irradiation, a higher number of differentially expressed genes were found 4 h after exposure to low and high dose ionizing radiation. Differences in gene expression were related to signal transduction pathways of the DNA damage response after 2 h and to metabolic pathways, that might implicate cellular senescence, after 4 h. The time point 4 h will be used to conduct further irradiation experiments in a larger sample.

MetaCache: context-aware classification of metagenomic reads using minhashing.

MOTIVATION: Metagenomic shotgun sequencing studies are becoming increasingly popular with prominent examples including the sequencing of human microbiomes and diverse environments. A fundamental computational problem in this context is read classification, i.e. the assignment of each read to a taxonomic label. Due to the large number of reads produced by modern high-throughput sequencing technologies and the rapidly increasing number of available reference genomes corresponding software tools suffer from either long runtimes, large memory requirements or low accuracy. RESULTS: We introduce MetaCache-a novel software for read classification using the big data technique minhashing. Our approach performs context-aware classification of reads by computing representative subsamples of k-mers within both, probed reads and locally constrained regions of the reference genomes. As a result, MetaCache consumes significantly less memory compared to the state-of-the-art read classifiers Kraken and CLARK while achieving highly competitive sensitivity and precision at comparable speed. For example, using NCBI RefSeq draft and completed genomes with a total length of around 140 billion bases as reference, MetaCache's database consumes only 62 GB of memory while both Kraken and CLARK fail to construct their respective databases on a workstation with 512 GB RAM. Our experimental results further show that classification accuracy continuously improves when increasing the amount of utilized reference genome data. AVAILABILITY AND IMPLEMENTATION: MetaCache is open source software written in C ++ and can be downloaded at http://github.com/muellan/metacache. CONTACT: bertil.schmidt@uni-mainz.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

AFS: identification and quantification of species composition by metagenomic sequencing.

Summary: DNA-based methods to detect and quantify taxon composition in biological materials are often based on species-specific polymerase chain reaction, limited to detecting species targeted by the assay. Next-generation sequencing overcomes this drawback by untargeted shotgun sequencing of whole metagenomes at affordable cost. Here we present AFS, a software pipeline for quantification of species composition in food. AFS uses metagenomic shotgun sequencing and sequence read counting to infer species proportions. Using Illumina data from a reference sausage comprising four species, we reveal that AFS is independent of the sequencing assay and library preparation protocol. Cost-saving short (50-bp) single-end reads and Nextera ® library preparation yield reliable results. Availability and Implementation: Datasets, binaries and usage instructions are available under http://all-food-seq.sourceforge.net. Raw data is available at NCBI's SRA with accession number PRJNA271645. Contact: hankeln@uni-mainz.de. Supplementary information: Supplementary data are available at Bioinformatics online.

The genomic footprint of climate adaptation in Chironomus riparius.

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool-Seq data and population genetic modelling. Common-garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an FST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.

Chironomus riparius (Diptera) genome sequencing reveals the impact of minisatellite transposable elements on population divergence.

Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky-Muller incompatibilities (DMI), thus contributing to the evolution of endogenous genetic barriers and eventually causing population divergence. We investigated differential TE abundance among conspecific populations of the nonbiting midge Chironomus riparius and evaluated their potential role in causing endogenous genetic incompatibilities between these populations. We focussed on a Chironomus-specific TE, the minisatellite-like Cla-element, whose activity is associated with speciation in the genus. Using a newly generated and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population-specific TE insertion patterns with many private insertions. A significant correlation of the pairwise FST estimated from genomewide single-nucleotide polymorphisms (SNPs) and the FST estimated from TEs is consistent with drift as the major force driving TE population differentiation. However, the significantly higher Cla-element FST level due to a high proportion of differentially fixed Cla-element insertions also indicates selection against segregating (i.e. heterozygous) insertions. With reciprocal crossing experiments and fluorescent in situ hybridization of Cla-elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings. We propose that the inferred negative selection on heterozygous Cla-element insertions may cause endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs may thus have a direct impact on population divergence that is operationally different from drift and local adaptation.

Extensive transcriptional complexity during hypoxia-regulated expression of the myoglobin gene in cancer.

Recently, the ectopic expression of myoglobin (MB) was reported in human epithelial cancer cell lines and breast tumor tissues, where MB expression increased with hypoxia. The better prognosis of MB-positive breast cancer patients suggested that the globin exerts a tumor-suppressive role, possibly by impairing mitochondrial activity in hypoxic breast carcinoma cells. To better understand MB gene regulation in cancer, we systematically investigated the architecture of the human MB gene, its transcripts and promoters. In silico analysis of transcriptome data from normal human tissues and cancer cell lines, followed by RACE-PCR verification, revealed seven novel exons in the MB gene region, most of which are untranslated exons located 5'-upstream of the coding DNA sequence (CDS). Sixteen novel alternatively spliced MB transcripts were detected, most of which predominantly occur in tumor tissue or cell lines. Quantitative RT-PCR analyses of MB expression in surgical breast cancer specimen confirmed the preferential usage of a hitherto unknown, tumor-associated MB promoter, which was functionally validated by luciferase reporter gene assays. In line with clinical observations of MB up-regulation in avascular breast tumors, the novel cancer-associated MB splice variants exhibited increased expression in tumor cells subjected to experimental hypoxia. The novel gene regulatory mechanisms unveiled in this study support the idea of a non-canonical role of MB during carcinogenesis.

Phylogeny of Syndermata (syn. Rotifera): Mitochondrial gene order verifies epizoic Seisonidea as sister to endoparasitic Acanthocephala within monophyletic Hemirotifera.

A monophyletic origin of endoparasitic thorny-headed worms (Acanthocephala) and wheel-animals (Rotifera) is widely accepted. However, the phylogeny inside the clade, be it called Syndermata or Rotifera, has lacked validation by mitochondrial (mt) data. Herein, we present the first mt genome of the key taxon Seison and report conflicting results of phylogenetic analyses: while mt sequence-based topologies showed monophyletic Lemniscea (Bdelloidea+Acanthocephala), gene order analyses supported monophyly of Pararotatoria (Seisonidea+Acanthocephala) and Hemirotifera (Bdelloidea+Pararotatoria). Sequence-based analyses obviously suffered from substitution saturation, compositional bias, and branch length heterogeneity; however, we observed no compromising effects in gene order analyses. Moreover, gene order-based topologies were robust to changes in coding (genes vs. gene pairs, two-state vs. multistate, aligned vs. non-aligned), tree reconstruction methods, and the treatment of the two monogonont mt genomes. Thus, mt gene order verifies seisonids as sister to acanthocephalans within monophyletic Hemirotifera, while deviating results of sequence-based analyses reflect artificial signal. This conclusion implies that the complex life cycle of extant acanthocephalans evolved from a free-living state, as retained by most monogononts and bdelloids, via an epizoic state with a simple life cycle, as shown by seisonids. Hence, Acanthocephala represent a rare example where ancestral transitional stages have counterparts amongst the closest relatives.

Isolation of Microsatellite Loci by Next-Generation Sequencing of the Critically Endangered Northern Bald Ibis, Geronticus eremita.

The Northern Bald Ibis is one of the rarest bird species, extinct in Europe for 400 years and critically endangered worldwide. The European Union-co-financed LIFE+ project "Reason for Hope - Reintroduction of the Northern Bald Ibis in Europe" aims to reintroduce the species in Europe (Germany, Austria, Italy). In order to obtain information on the genetic diversity within zoo colonies and the reintroduced population, 15 polymorphic microsatellite markers, specific for the Northern Bald Ibis, Geronticus eremita (Linnaeus, 1785), have been isolated from next-generation sequencing (Illumina MiSeq) and are described here. The microsatellite primers were tested in 30 individuals and measures of genetic variability were calculated. Values for the observed heterozygosity ranged from 0.393 to 0.867, while expected heterozygosity ranged from 0.573 to 0.718. Ten out of 15 loci were in Hardy-Weinberg equilibrium and only one showed indication for the presence of null alleles. The newly developed PCR primers can be used to examine population genetic parameters, e.g. for future conservation genetic studies of this critically endangered bird species.

The globin gene repertoire of lampreys: convergent evolution of hemoglobin and myoglobin in jawed and jawless vertebrates.

Agnathans (jawless vertebrates) occupy a key phylogenetic position for illuminating the evolution of vertebrate anatomy and physiology. Evaluation of the agnathan globin gene repertoire can thus aid efforts to reconstruct the origin and evolution of the globin genes of vertebrates, a superfamily that includes the well-known model proteins hemoglobin and myoglobin. Here, we report a comprehensive analysis of the genome of the sea lamprey (Petromyzon marinus) which revealed 23 intact globin genes and two hemoglobin pseudogenes. Analyses of the genome of the Arctic lamprey (Lethenteron camtschaticum) identified 18 full length and five partial globin gene sequences. The majority of the globin genes in both lamprey species correspond to the known agnathan hemoglobins. Both genomes harbor two copies of globin X, an ancient globin gene that has a broad phylogenetic distribution in the animal kingdom. Surprisingly, we found no evidence for an ortholog of neuroglobin in the lamprey genomes. Expression and phylogenetic analyses identified an ortholog of cytoglobin in the lampreys; in fact, our results indicate that cytoglobin is the only orthologous vertebrate-specific globin that has been retained in both gnathostomes and agnathans. Notably, we also found two globins that are highly expressed in the heart of P. marinus, thus representing functional myoglobins. Both genes have orthologs in L. camtschaticum. Phylogenetic analyses indicate that these heart-expressed globins are not orthologous to the myoglobins of jawed vertebrates (Gnathostomata), but originated independently within the agnathans. The agnathan myoglobin and hemoglobin proteins form a monophyletic group to the exclusion of functionally analogous myoglobins and hemoglobins of gnathostomes, indicating that specialized respiratory proteins for O2 transport in the blood and O2 storage in the striated muscles evolved independently in both lineages. This dual convergence of O2-transport and O2-storage proteins in agnathans and gnathostomes involved the convergent co-option of different precursor proteins in the ancestral globin repertoire of vertebrates.