Tyzzerella nexilis strains enriched in mobile genetic elements are involved in progressive multiple sclerosis

Multiple sclerosis (MS) is an autoimmune-demyelinating disease with an inflammatory pathology formed by self-reactive lymphocytes with activated glial cells. Progressive MS, characterized by resistance to medications, significantly differs from the non-progressive form in gut microbiome profiles. After confirming an increased abundance of "Tyzzerella nexilis" in various cohorts of progressive MS, we identified a distinct cluster of T. nexilis strains enriched in progressive MS based on long-read metagenomics. The distinct T. nexilis cluster is characterized by a large number of mobile genetic elements (MGEs) and a lack of defense systems against MGEs. Microbial genes for sulfate reduction and flagella formation with pathogenic implications are specific to this cluster. Moreover, these flagellar genes are encoded on MGEs. Mono-colonization with MGE-enriched T. nexilis made germ-free mice more susceptible to experimental autoimmune encephalomyelitis. These results indicate that the progression of MS may be promoted by MGE-enriched T. nexilis with potentially pathogenic properties.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome wide

RNA:DNA hybrids compromise replication fork progression and genome integrity in all cells. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate strongly in noncoding RNAs and 5'-UTRs of coding sequences. For ΔrnhB, hybrids accumulate preferentially in untranslated regions and early in coding sequences. We show that hybrid accumulation is particularly sensitive to gene expression in ΔrnhC cells. DNA replication in ΔrnhC cells is disrupted, leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts cause mutagenesis and shape genome organization.

An insight into the functional genomics and species classification of Eudiplozoon nipponicum (Monogenea, Diplozoidae), a haematophagous parasite of the common carp Cyprinus carpio

BACKGROUND

Monogenea (Platyhelminthes, Neodermata) are the most species-rich class within the Neodermata superclass of primarily fish parasites. Despite their economic and ecological importance, monogenean research tends to focus on their morphological, phylogenetic, and population characteristics, while comprehensive omics analyses aimed at describing functionally important molecules are few and far between. We present a molecular characterisation of monogenean representative Eudiplozoon nipponicum, an obligate haematophagous parasite infecting the gills of the common carp. We report its nuclear and mitochondrial genomes, present a functional annotation of protein molecules relevant to the molecular and biochemical aspect of physiological processes involved in interactions with the fish hosts, and re-examinate the taxonomic position of Eudiplozoon species within the Diplozoidae family.

RESULTS

We have generated 50.81 Gbp of raw sequencing data (Illumina and Oxford Nanopore reads), bioinformatically processed, and de novo assembled them into a genome draft 0.94 Gbp long, consisting of 21,044 contigs (N50 = 87 kbp). The final assembly represents 57% of the estimated total genome size (~ 1.64 Gbp), whereby repetitive and low-complexity regions account for ~ 64% of the assembled length. In total, 36,626 predicted genes encode 33,031 proteins and homology-based annotation of protein-coding genes (PCGs) and proteins characterises 14,785 (44.76%) molecules. We have detected significant representation of functional proteins and known molecular functions. The numbers of peptidases and inhibitors (579 proteins), characterised GO terms (16,016 unique assigned GO terms), and identified KEGG Orthology (4,315 proteins) acting in 378 KEGG pathways demonstrate the variety of mechanisms by which the parasite interacts with hosts on a macromolecular level (immunomodulation, feeding, and development). Comparison between the newly assembled E. nipponicum mitochondrial genome (length of 17,038 bp) and other diplozoid monogeneans confirms the existence of two distinct Eudiplozoon species infecting different fish hosts: Cyprinus carpio and Carassius spp.

CONCLUSIONS

Although the amount of sequencing data and characterised molecules of monogenean parasites has recently increased, a better insight into their molecular biology is needed. The E. nipponicum nuclear genome presented here, currently the largest described genome of any monogenean parasite, represents a milestone in the study of monogeneans and their molecules but further omics research is needed to understand these parasites' biological nature.

RNase H genes cause distinct impacts on RNA:DNA hybrid formation and mutagenesis genome-wide

RNA:DNA hybrids such as R-loops affect genome integrity and DNA replication fork progression. The overall impacts of naturally occurring RNA:DNA hybrids on genome integrity, and the relative contributions of ribonucleases H to mitigating the negative effects of hybrids, remain unknown. Here, we investigate the contributions of RNases HII (RnhB) and HIII (RnhC) to hybrid removal, DNA replication, and mutagenesis genome-wide. Deletion of either rnhB or rnhC triggers RNA:DNA hybrid accumulation, but with distinct patterns of mutagenesis and hybrid accumulation. Across all cells, hybrids accumulate most strongly in non-coding RNAs and 5'-UTRs of coding sequences. For Δ rnhB , hybrids accumulate preferentially in untranslated regions and early in coding sequences. Hybrid accumulation is particularly sensitive to gene expression in Δ rnhC ; in cells lacking RnhC, DNA replication is disrupted leading to transversions and structural variation. Our results resolve the outstanding question of how hybrids in native genomic contexts interact with replication to cause mutagenesis and shape genome organization.

Adaptive Radioresistance of Enterohemorrhagic Escherichia coli O157:H7 Results in Genomic Loss of Shiga Toxin-Encoding Prophages

Enterohemorrhagic Escherichia coli (EHEC) O157:H7 is a foodborne pathogen producing Shiga toxins (Stx1 and Stx2), which can cause hemorrhagic diarrhea and life-threatening infections. O157:H7 strain EDL933 carries prophages CP-933V and BP-933W, which encode Shiga toxin genes (stx1 and stx2, respectively). The aim of this work was to investigate the mechanisms of adaptive resistance of EHEC strain EDL933 to a typically lethal dose of gamma irradiation (1.5 kGy). Adaptive selection through six passages of exposure to 1.5 kGy resulted in the loss of CP-933V and BP-933W prophages from the genome and mutations within three genes: wrbA, rpoA, and Wt_02639 (molY). Three selected EHEC clones that became irradiation adapted to the 1.5-kGy dose (C1, C2, and C3) demonstrated increased resistance to oxidative stress, sensitivity to acid pH, and decreased cytotoxicity to Vero cells. To confirm that loss of prophages plays a role in increased radioresistance, clones C1 and C2 were exposed to bacteriophage-containing lysates. Although phage BP-933W could lysogenize C1, C2, and E. coli K-12 strain MG1655, it was not found to have integrated into the bacterial chromosome in C1-Φ and C2-Φ lysogens. Interestingly, for the E. coli K-12 lysogen (K-12-Φ), BP-933W DNA had integrated at the wrbA gene (K-12-Φ). Both C1-Φ and C2-Φ lysogens regained sensitivity to oxidative stress, were more effectively killed by a 1.5-kGy gamma irradiation dose, and had regained cytotoxicity and acid resistance phenotypes. Further, the K-12-Φ lysogen became cytotoxic, more sensitive to gamma irradiation and oxidative stress, and slightly more acid resistant. IMPORTANCE Gamma irradiation of food products can provide an effective means of eliminating bacterial pathogens such as enterohemorrhagic Escherichia coli (EHEC) O157:H7, a significant foodborne pathogen that can cause severe disease due to the production of Stx. To decipher the mechanisms of adaptive resistance of the O157:H7 strain EDL933, we evolved clones of this bacterium resistant to a lethal dose of gamma irradiation by repeatedly exposing bacterial cells to irradiation following a growth restoration over six successive passages. Our findings provide evidence that adaptive selection involved modifications in the bacterial genome, including deletion of the CP-933V and BP-933W prophages. These mutations in EHEC O157:H7 resulted in loss of stx1 and stx2, loss of cytotoxicity to epithelial cells, and decreased resistance to acidity, critical virulence determinants of EHEC, concomitant with increased resistance to lethal irradiation and oxidative stress. These findings demonstrate that the potential adaptation of EHEC to high doses of radiation would involve elimination of the Stx-encoding phages and likely lead to a substantial attenuation of virulence.

Isolation and infection cycle of a polinton-like virus virophage in an abundant marine alga

Virophages are small double stranded DNA (dsDNA) viruses that can only replicate in a host by co-infecting with another virus. Marine algae are commonly associated with virophage-like elements such as Polinton-like viruses (PLVs) that remain largely uncharacterized. Here we isolated a PLV that co-infects the alga Phaeocystis globosa with the Phaeocystis globosa virus-14T (PgV-14T), a close relative of the "Phaeocystis globosa virus-virophage" genomic sequence. We name this PLV 'Gezel-14T. Gezel is phylogenetically distinct from the Lavidaviridae family where all known virophages belong. Gezel-14T co-infection decreases the fitness of its viral host by reducing burst sizes of PgV-14T, yet insufficiently to spare the cellular host population. Genomic screens show Gezel-14T-like PLVs integrated into Phaeocystis genomes, suggesting that these widespread viruses are capable of integration into cellular host genomes. This system presents an opportunity to better understand the evolution of eukaryotic dsDNA viruses as well as the complex dynamics and implications of viral parasitism.

Phylogeny of Saxifraga section Saxifraga subsection Arachnoideae (Saxifragaceae) and the origin of low elevation shade-dwelling species

Saxifraga section Saxifraga subsection Arachnoideae is a lineage of 12 species distributed mainly in the European Alps. It is unusual in terms of ecological diversification by containing both high elevation species from exposed alpine habitats and low elevation species from shady habitats such as overhanging rocks and cave entrances. Our aims are to explore which of these habitat types is ancestral, and to identify the possible drivers of this remarkable ecological diversification. Using a Hybseq DNA-sequencing approach and a complete species sample we reconstructed and dated the phylogeny of subsection Arachnoideae. Using Landolt indicator values, this phylogenetic tree was used for the reconstruction of the evolution of temperature, light and soil pH requirements in this lineage. Diversification of subsection Arachnoideae started in the late Pliocene and continued through the Pleistocene. Both diversification among and within clades was largely allopatric, and species from shady habitats with low light requirements are distributed in well-known refugia. We hypothesize that low light requirements evolved when species persisting in cold-stage refugia were forced into marginal habitats by more competitive warm-stage vegetation. While we do not claim that such competition resulted in speciation, it very likely resulted in adaptive evolution.

Identification of trypsin-degrading commensals in the large intestine

Increased levels of proteases, such as trypsin, in the distal intestine have been implicated in intestinal pathological conditions^1–3. However, the players and mechanisms that underlie protease regulation in the intestinal lumen have remained unclear. Here we show that Paraprevotella strains isolated from the faecal microbiome of healthy human donors are potent trypsin-degrading commensals. Mechanistically, Paraprevotella recruit trypsin to the bacterial surface through type IX secretion system-dependent polysaccharide-anchoring proteins to promote trypsin autolysis. Paraprevotella colonization protects IgA from trypsin degradation and enhances the effectiveness of oral vaccines against Citrobacter rodentium. Moreover, Paraprevotella colonization inhibits lethal infection with murine hepatitis virus-2, a mouse coronavirus that is dependent on trypsin and trypsin-like proteases for entry into host cells^4,5. Consistently, carriage of putative genes involved in trypsin degradation in the gut microbiome was associated with reduced severity of diarrhoea in patients with SARS-CoV-2 infection. Thus, trypsin-degrading commensal colonization may contribute to the maintenance of intestinal homeostasis and protection from pathogen infection.

Colonization of trypsin-degrading commensal bacteria may contribute to the maintenance of intestinal homeostasis and protection against pathogen infection in humans and mice.

An invasive Haemophilus influenzae serotype b infection in an Anglo-Saxon plague victim

BACKGROUND

The human pathogen Haemophilus influenzae was the main cause of bacterial meningitis in children and a major cause of worldwide infant mortality before the introduction of a vaccine in the 1980s. Although the occurrence of serotype b (Hib), the most virulent type of H. influenzae, has since decreased, reports of infections with other serotypes and non-typeable strains are on the rise. While non-typeable strains have been studied in-depth, very little is known of the pathogen's evolutionary history, and no genomes dating prior to 1940 were available.

RESULTS

We describe a Hib genome isolated from a 6-year-old Anglo-Saxon plague victim, from approximately 540 to 550 CE, Edix Hill, England, showing signs of invasive infection on its skeleton. We find that the genome clusters in phylogenetic division II with Hib strain NCTC8468, which also caused invasive disease. While the virulence profile of our genome was distinct, its genomic similarity to NCTC8468 points to mostly clonal evolution of the clade since the 6th century. We also reconstruct a partial Yersinia pestis genome, which is likely identical to a published first plague pandemic genome of Edix Hill.

CONCLUSIONS

Our study presents the earliest genomic evidence for H. influenzae, points to the potential presence of larger genomic diversity in the phylogenetic division II serotype b clade in the past, and allows the first insights into the evolutionary history of this major human pathogen. The identification of both plague and Hib opens questions on the effect of plague in immunocompromised individuals already affected by infectious diseases.

Prediction of the Sex-Associated Genomic Region in Tunas (Thunnus Fishes)

Fish species have a variety of sex determination systems. Tunas (genus Thunnus) have an XY genetic sex determination system. However, the Y chromosome or responsible locus has not yet been identified in males. In a previous study, a female genome of Pacific bluefin tuna (T. orientalis) was sequenced, and candidates for sex-associated DNA polymorphisms were identified by a genome-wide association study using resequencing data. In the present study, we sequenced a male genome of Pacific bluefin tuna by long-read and linked-read sequencing technologies and explored male-specific loci through a comparison with the female genome. As a result, we found a unique region carrying the male-specific haplotype, where a homolog of estrogen sulfotransferase gene was predicted to be encoded. The genome-wide mapping of previously resequenced data indicated that, among the functionally annotated genes, only this gene, named sult1st6y, was paternally inherited in the males of Pacific bluefin tuna. We reviewed the RNA-seq data of southern bluefin tuna (T. maccoyii) in the public database and found that sult1st6y of southern bluefin tuna was expressed in all male testes, but absent or suppressed in the female ovary. Since estrogen sulfotransferase is responsible for the inactivation of estrogens, it is reasonable to assume that the expression of sult1st6y in gonad cells may inhibit female development, thereby inducing the individuals to become males. Thus, our results raise a promising hypothesis that sult1st6y is the sex determination gene in Thunnus fishes or at least functions at a crucial point in the sex-differentiation cascade.

Minirmd: accurate and fast duplicate removal tool for short reads via multiple minimizers

SUMMARY

Removing duplicate and near-duplicate reads, generated by high-throughput sequencing technologies, is able to reduce computational resources in downstream applications. Here we develop minirmd, a de novo tool to remove duplicate reads via multiple rounds of clustering using different length of minimizer. Experiments demonstrate that minirmd removes more near-duplicate reads than existing clustering approaches and is faster than existing multi-core tools. To the best of our knowledge, minirmd is the first tool to remove near-duplicates on reverse-complementary strand.

AVAILABILITY AND IMPLEMENTATION

https://github.com/yuansliu/minirmd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Myoglobin primary structure reveals multiple convergent transitions to semi-aquatic life in the world's smallest mammalian divers

The speciose mammalian order Eulipotyphla (moles, shrews, hedgehogs, solenodons) combines an unusual diversity of semi-aquatic, semi-fossorial, and fossorial forms that arose from terrestrial forbearers. However, our understanding of the ecomorphological pathways leading to these lifestyles has been confounded by a fragmentary fossil record, unresolved phylogenetic relationships, and potential morphological convergence, calling for novel approaches. The net surface charge of the oxygen-storing muscle protein myoglobin (ZMb), which can be readily determined from its primary structure, provides an objective target to address this question due to mechanistic linkages with myoglobin concentration. Here, we generate a comprehensive 71 species molecular phylogeny that resolves previously intractable intra-family relationships and then ancestrally reconstruct ZMb evolution to identify ancient lifestyle transitions based on protein sequence alone. Our phylogenetically informed analyses confidently resolve fossorial habits having evolved twice in talpid moles and reveal five independent secondary aquatic transitions in the order housing the world's smallest endothermic divers.

Accelerating Population Count with a Hardware Co-Processor for MicroBlaze

This paper proposes a Field-Programmable Gate Array (FPGA)-based hardware accelerator for assisting the embedded MicroBlaze soft-core processor in calculating population count. The population count is frequently required to be executed in cyber-physical systems and can be applied to large data sets, such as in the case of molecular similarity search in cheminformatics, or assisting with computations performed by binarized neural networks. The MicroBlaze instruction set architecture (ISA) does not support this operation natively, so the count has to be realized as either a sequence of native instructions (in software) or in parallel in a dedicated hardware accelerator. Different hardware accelerator architectures are analyzed and compared to one another and to implementing the population count operation in MicroBlaze. The achieved experimental results with large vector lengths (up to 217) demonstrate that the best hardware accelerator with DMA (Direct Memory Access) is ~31 times faster than the best software version running on MicroBlaze. The proposed architectures are scalable and can easily be adjusted to both smaller and bigger input vector lengths. The entire system was implemented and tested on a Nexys-4 prototyping board containing a low-cost/low-power Artix-7 FPGA.

Nubeam-dedup: a fast and RAM-efficient tool to de-duplicate sequencing reads without mapping

SUMMARY

We present Nubeam-dedup, a fast and RAM-efficient tool to de-duplicate sequencing reads without reference genome. Nubeam-dedup represents nucleotides by matrices, transforms reads into products of matrices, and based on which assigns a unique number to a read. Thus, duplicate reads can be efficiently removed by using a collisionless hash function. Compared with other state-of-the-art reference-free tools, Nubeam-dedup uses 50-70% of CPU time and 10-15% of RAM.

AVAILABILITY AND IMPLEMENTATION

Source code in C++ and manual are available at https://github.com/daihang16/nubeamdedup and https://haplotype.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

High Dietary Sugar Reshapes Sweet Taste to Promote Feeding Behavior in Drosophila melanogaster

Recent studies find that sugar tastes less intense to humans with obesity, but whether this sensory change is a cause or a consequence of obesity is unclear. To tackle this question, we study the effects of a high sugar diet on sweet taste sensation and feeding behavior in Drosophila melanogaster. On this diet, fruit flies have lower taste responses to sweet stimuli, overconsume food, and develop obesity. Excess dietary sugar, but not obesity or dietary sweetness alone, caused taste deficits and overeating via the cell-autonomous action of the sugar sensor O-linked N-Acetylglucosamine (O-GlcNAc) transferase (OGT) in the sweet-sensing neurons. Correcting taste deficits by manipulating the excitability of the sweet gustatory neurons or the levels of OGT protected animals from diet-induced obesity. Our work demonstrates that the reshaping of sweet taste sensation by excess dietary sugar drives obesity and highlights the role of glucose metabolism in neural activity and behavior.

BioSeqZip: a collapser of NGS redundant reads for the optimization of sequence analysis

MOTIVATION

High-throughput next-generation sequencing can generate huge sequence files, whose analysis requires alignment algorithms that are typically very demanding in terms of memory and computational resources. This is a significant issue, especially for machines with limited hardware capabilities. As the redundancy of the sequences typically increases with coverage, collapsing such files into compact sets of non-redundant reads has the 2-fold advantage of reducing file size and speeding-up the alignment, avoiding to map the same sequence multiple times.

METHOD

BioSeqZip generates compact and sorted lists of alignment-ready non-redundant sequences, keeping track of their occurrences in the raw files as well as of their quality score information. By exploiting a memory-constrained external sorting algorithm, it can be executed on either single- or multi-sample datasets even on computers with medium computational capabilities. On request, it can even re-expand the compacted files to their original state.

RESULTS

Our extensive experiments on RNA-Seq data show that BioSeqZip considerably brings down the computational costs of a standard sequence analysis pipeline, with particular benefits for the alignment procedures that typically have the highest requirements in terms of memory and execution time. In our tests, BioSeqZip was able to compact 2.7 billion of reads into 963 million of unique tags reducing the size of sequence files up to 70% and speeding-up the alignment by 50% at least.

AVAILABILITY AND IMPLEMENTATION

BioSeqZip is available at https://github.com/bioinformatics-polito/BioSeqZip.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Improvement of the Pacific bluefin tuna (Thunnus orientalis) reference genome and development of male-specific DNA markers

The Pacific bluefin tuna, Thunnus orientalis, is a highly migratory species that is widely distributed in the North Pacific Ocean. Like other marine species, T. orientalis has no external sexual dimorphism; thus, identifying sex-specific variants from whole genome sequence data is a useful approach to develop an effective sex identification method. Here, we report an improved draft genome of T. orientalis and male-specific DNA markers. Combining PacBio long reads and Illumina short reads sufficiently improved genome assembly, with a 38-fold increase in scaffold contiguity (to 444 scaffolds) compared to the first published draft genome. Through analysing re-sequence data of 15 males and 16 females, 250 male-specific SNPs were identified from more than 30 million polymorphisms. All male-specific variants were male-heterozygous, suggesting that T. orientalis has a male heterogametic sex-determination system. The largest linkage disequilibrium block (3,174 bp on scaffold_064) contained 51 male-specific variants. PCR primers and a PCR-based sex identification assay were developed using these male-specific variants. The sex of 115 individuals (56 males and 59 females; sex was diagnosed by visual examination of the gonads) was identified with high accuracy using the assay. This easy, accurate, and practical technique facilitates the control of sex ratios in tuna farms. Furthermore, this method could be used to estimate the sex ratio and/or the sex-specific growth rate of natural populations.

Divergent evolutionary trajectories following speciation in two ectoparasitic honey bee mites

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities.

NGSReadsTreatment - A Cuckoo Filter-based Tool for Removing Duplicate Reads in NGS Data

The Next-Generation Sequencing (NGS) platforms provide a major approach to obtaining millions of short reads from samples. NGS has been used in a wide range of analyses, such as for determining genome sequences, analyzing evolutionary processes, identifying gene expression and resolving metagenomic analyses. Usually, the quality of NGS data impacts the final study conclusions. Moreover, quality assessment is generally considered the first step in data analyses to ensure the use of only reliable reads for further studies. In NGS platforms, the presence of duplicated reads (redundancy) that are usually introduced during library sequencing is a major issue. These might have a serious impact on research application, as redundancies in reads can lead to difficulties in subsequent analysis (e.g., de novo genome assembly). Herein, we present NGSReadsTreatment, a computational tool for the removal of duplicated reads in paired-end or single-end datasets. NGSReadsTreatment can handle reads from any platform with the same or different sequence lengths. Using the probabilistic structure Cuckoo Filter, the redundant reads are identified and removed by comparing the reads with themselves. Thus, no prerequisite is required beyond the set of reads. NGSReadsTreatment was compared with other redundancy removal tools in analyzing different sets of reads. The results demonstrated that NGSReadsTreatment was better than the other tools in both the amount of redundancies removed and the use of computational memory for all analyses performed. Available in https://sourceforge.net/projects/ngsreadstreatment/.

Shifts in the Genetic Landscape of the Western Eurasian Steppe Associated with the Beginning and End of the Scythian Dominance

The Early Iron Age nomadic Scythians have been described as a confederation of tribes of different origins, based on ancient DNA evidence [1-3]. It is still unclear how much of the Scythian dominance in the Eurasian Steppe was due to movements of people and how much reflected cultural diffusion and elite dominance. We present new whole-genome sequences of 31 ancient Western and Eastern Steppe individuals, including Scythians as well as samples pre- and postdating them, allowing us to set the Scythians in a temporal context (in the Western, i.e., Ponto-Caspian Steppe). We detect an increase of eastern (Altaian) affinity along with a decrease in eastern hunter-gatherer (EHG) ancestry in the Early Iron Age Ponto-Caspian gene pool at the start of the Scythian dominance. On the other hand, samples of the Chernyakhiv culture postdating the Scythians in Ukraine have a significantly higher proportion of Near Eastern ancestry than other samples of this study. Our results agree with the Gothic source of the Chernyakhiv culture and support the hypothesis that the Scythian dominance did involve a demic component.

MPIGeneNet: Parallel Calculation of Gene Co-Expression Networks on Multicore Clusters

In this work, we present MPIGeneNet, a parallel tool that applies Pearson's correlation and Random Matrix Theory to construct gene co-expression networks. It is based on the state-of-the-art sequential tool RMTGeneNet, which provides networks with high robustness and sensitivity at the expenses of relatively long runtimes for large scale input datasets. MPIGeneNet returns the same results as RMTGeneNet but improves the memory management, reduces the I/O cost, and accelerates the two most computationally demanding steps of co-expression network construction by exploiting the compute capabilities of common multicore CPU clusters. Our performance evaluation on two different systems using three typical input datasets shows that MPIGeneNet is significantly faster than RMTGeneNet. As an example, our tool is up to 175.41 times faster on a cluster with eight nodes, each one containing two 12-core Intel Haswell processors. The source code of MPIGeneNet, as well as a reference manual, are available at https://sourceforge.net/projects/mpigenenet/.

paraGSEA: a scalable approach for large-scale gene expression profiling

More studies have been conducted using gene expression similarity to identify functional connections among genes, diseases and drugs. Gene Set Enrichment Analysis (GSEA) is a powerful analytical method for interpreting gene expression data. However, due to its enormous computational overhead in the estimation of significance level step and multiple hypothesis testing step, the computation scalability and efficiency are poor on large-scale datasets. We proposed paraGSEA for efficient large-scale transcriptome data analysis. By optimization, the overall time complexity of paraGSEA is reduced from O(mn) to O(m+n), where m is the length of the gene sets and n is the length of the gene expression profiles, which contributes more than 100-fold increase in performance compared with other popular GSEA implementations such as GSEA-P, SAM-GS and GSEA2. By further parallelization, a near-linear speed-up is gained on both workstations and clusters in an efficient manner with high scalability and performance on large-scale datasets. The analysis time of whole LINCS phase I dataset (GSE92742) was reduced to nearly half hour on a 1000 node cluster on Tianhe-2, or within 120 hours on a 96-core workstation. The source code of paraGSEA is licensed under the GPLv3 and available at http://github.com/ysycloud/paraGSEA.