nitt/bio-med-arxiv huggingface dataset:
Longitudinal studies at birth and age 7 reveal strong effects of genetic variation on ancestry-associated DNA methylation patterns in blood cells from ethnically admixed children
Authors: McKennan, C.; Naughton, K.; Stanhope, C.; Kattan, M.; O'Connor, G.; Sandel, M.; Visness, C.; Wood, R.; Bacharier, L.; Beigelman, A.; Lovisky-Desir, S.; Togias, A.; Gern, J.; Nicolae, D.; Ober, C.
Category: genomics
Date: 2018-08-24

Epigenetic architecture is influenced by genetic and environmental factors, but little is known about their relative contributions or longitudinal dynamics. Here, we studied DNA methylation (DNAm) at over 750,000 CpG sites in mononuclear blood cells collected at birth and age 7 from 196 children of primarily self-reported Black and Hispanic ethnicities to study race-associated DNAm patterns. We developed a novel Bayesian method for high dimensional longitudinal data and showed that race-associated DNAm patterns at birth and age 7 are nearly identical. Additionally, we estimated that up to 51% of all self-reported race-associated CpGs had race-dependent DNAm levels that were mediated through local genotype and, quite surprisingly, found that genetic factors explained an overwhelming majority of the variation in DNAm levels at other, previously identified, environmentally-associated CpGs. These results not only indicate that race-associated DNAm patterns in blood are present at birth and are primarily genetically, and not environmentally, determined, but also that DNAm in blood cells overall is robust to many environmental exposures during the first 7 years of life.
Cerebrovascular Injuries Induce Lymphatic Invasion into Brain Parenchyma to Guide Vascular Regeneration in Zebrafish
Authors: Chen, J.; He, J.; Yang, Q.; Zhang, Y.; Luo, L.
Category: developmental biology
Date: 2018-08-23

Damage to regional cerebrovascular network and neuronal tissues occurs during acute cerebrovascular diseases, such as ischemic stroke. The promotion of vascular regeneration is the most promising therapeutic approach. To understand cellular and molecular mechanisms underlying brain vascular regeneration, we developed two zebrafish cerebrovascular injury models using genetic ablation and photochemical thrombosis. Although brain parenchyma is physiologically devoid of lymphatic vasculature, we found that cerebrovascular injuries induce rapid ingrowth of meningeal lymphatics into the injured parenchyma. The ingrown lymphatics on one hand become lumenized drain interstitial fluid to resolve brain edema, on the other hand act as \"growing tracks\" for nascent blood vessels. The ingrown lymphatic vessels undergo apoptosis and clearance after cerebrovascular regeneration. This study reveals a pathological function of meningeal lymphatics, through previously unexpected ingrowth into brain parenchyma and a newly identified lymphatic function as vascular \"growing tracks\".\n\nHIGHLIGHTSCerebrovascular injuries induce lymphatic ingrowth into the injured brain parenchyma The ingrown lymphatics drain interstitial fluid to resolve brain edema Nascent blood vessels use the ingrown lymphatic vessels as \"growing tracks\" The ingrown lymphatic vessels undergo apoptosis after vascular regeneration completes
Gut bacterial tyrosine decarboxylases restrict the bioavailability of levodopa, the primary treatment in Parkinson’s disease
Authors: van Kessel, S. P.; Frye, A. K.; El-Gendy, A. O.; Castejon, M.; Keshavarzian, A.; van Dijk, G.; El Aidy, S.
Category: molecular biology
Date: 2018-08-21

Human gut bacteria play a critical role in the regulation of immune and metabolic systems, as well as in the function of the nervous system. The microbiota senses its environment and responds by releasing metabolites, some of which are key regulators of human health and disease. In this study, we identify and characterize gut-associated bacteria in their ability to decarboxylate L-DOPA (also known as Levodopa or L-3,4-dihydroxyphenylalanine) to dopamine via the tyrosine decarboxylases, which are mainly present in the class Bacilli. Although the bacterial tyrosine decarboxylases have a higher affinity for tyrosine compared to L-DOPA, this does not affect their ability to decarboxylate L-DOPA, nor does any inhibitor of the human decarboxylase. This study indicates that in situ bioavailability of L-DOPA is compromised by the gut bacterial tyrosine decarboxylase abundance in Parkinsons patients. Finally, we show that the tyrosine decarboxylase abundance in the microbiota at the site of L-DOPA absorption, the proximal small intestine, significantly influences L-DOPA bioavailability in the plasma of rats. Our results highlight the role of microbial metabolism in drug bioavailability, and specifically, that small intestinal abundance of bacterial tyrosine decarboxylase can explain the highly variable L-DOPA dosage regimens required in the treatment of individual Parkinsons patients.\n\nHighlightsO_LISmall intestinal bacteria is able to convert L-DOPA to dopamine\nC_LIO_LIL-DOPA metabolism by gut bacteria reduce the bioavailability of L-DOPA in the body, thus is a significant explanatory factor of the highly variable L-DOPA dosage regimens required in the treatment of individual Parkinsons patients.\nC_LIO_LIInhibitors of the human DOPA decarboxylase are not potent inhibitors for bacterial tyrosine decarboxylases\nC_LI
DataRemix: a universal data transformation for optimal inference from gene expression datasets
Authors: Mao, W.; Hausler, R.; Chikina, M.
Category: genomics
Date: 2018-08-27

MotivationRNAseq technology provides unprecedented power in the assesment of the transcription abundance and can be used to perform a variety of downstream tasks such as inference of gene-correlation network and eQTL discovery. However, raw gene expression values have to be normalized for nuisance biological variation and technical covariates, and different normalization strategies can lead to dramatically different results in the downstream study.

ResultsWe describe a generalization of SVD-based reconstruction for which the common techniques of whitening, rank-k approximation, and removing the top k principle components are special cases. Our simple three-parameter transformation, DataRemix, can be tuned to reweight the contribution of hidden factors and reveal otherwise hidden biological signals. In particular, we demonstrate that the method can effectively prioritize biological signals over noise without leveraging external dataset-specific knowledge, and can outperform normalization methods that make explicit use of known technical factors. We also show that DataRemix can be efficiently optimized via Thompson Sampling approach, which makes it feasible for computationally expensive objectives such as eQTL analysis. Finally we reanalyze the Depression Gene Networks (DGN) dataset, and we highlight new trans-eQTL networks which were not reported in the initial study.

AvailabilityDataRemix is an R package which is freely available at GitHub (https://github.com/wgmao/DataRemix).

Contactmchikina@pitt.edu
Phylogeographic analysis and identification of factors impacting the diffusion of Foot-and-Mouth disease virus in Africa
Authors: Duchatel, F.; Bronsvoort, M.; Lycett, S.
Category: evolutionary biology
Date: 2018-08-23

Foot and mouth disease (FMD) is endemic in sub-Saharan Africa and can lead to important and continuous economic losses for affected countries. Due to the complexity of the disease epidemiology and the lack of data there is a need to use inferential computational approaches to fill the gaps in our understanding of the circulation of FMD virus on this continent. Using a phylogeographic approach we reconstructed the circulation of FMD virus serotypes A, O and SAT2 in Africa and evaluated the influence of potential environmental and anthropological predictors of virus diffusion. Our results show that over the last hundred year the continental circulation of the tree serotypes was mainly driven by livestock trade. Whilst our analyses show that the serotypes A and O were introduced in Africa trough livestock trades, the SAT2 serotype probably originates from African wildlife population. The circulation of serotype O in eastern Africa is impacted by both indirect transmission through persistence in the environment and anthropological activities such as cattle movements.
Genetic variation in a complex polyploid: unveiling the dynamic allelic features of sugarcane
Authors: Sforca, D. A.; Vautrin, S.; Cardoso-Silva, C. B.; Mancini, M. C.; Cruz, M. V. R.; Pereira, G. d. S.; Conte, M.; Bellec, A.; Dahmer, N.; Fourment, J.; Rodde, N.; Van-Sluys, M.-A.; Vicentini, R.; Garcia, A. A. F.; Forni-Martins, E. R.; Carneiro, M. S.; Hoffmann, H. P.; Pinto, L. R.; Landell, M. G. d. A.; Vincentz, M.; Berges, H.; Souza, A. P.
Category: plant biology
Date: 2018-08-21

BackgroundSugarcane (Saccharum spp.) is highly polyploid and aneuploid. Modern cultivars are derived from hybridization between S. officinarum and S. spontaneum. This combination results in a genome exhibiting variable ploidy among different loci, a huge genome size (approximately 10 Gb) and a high content of repetitive regions. Gene expression mechanisms are poorly understood in these cultivars. An approach using genomic, transcriptomic and genetic mapping can improve our knowledge of the behavior of genetics in sugarcane.\n\nResultsThe hypothetical HP600 and centromere protein C (CENP-C) genes from sugarcane were used to elucidate the allelic expression and genomic and genetic behavior of this complex polyploid. The genomically side-by-side genes HP600 and CENP-C were found in two different homeologous chromosome groups with ploidies of eight and ten. The first region (Region01) was a Sorghum bicolor ortholog with all haplotypes of HP600 and CENP- C expressed, but HP600 exhibited an unbalanced haplotype expression. The second region (Region02) was a scrambled sugarcane sequence formed from different noncollinear genes containing duplications of HP600 and CENP-C (paralogs). This duplication occurred before the Saccharum genus formation and after the separation of sorghum and sugarcane, resulting in a nonexpressed HP600 pseudogene and a recombined fusion version of CENP-C and orthologous gene Sobic.003G299500 with at least two chimerical gene haplotypes expressed. The genetic map construction supported the difficulty of mapping markers located in duplicated regions of complex polyploid genomes.\n\nConclusionAll these findings describe a low synteny region in sugarcane, formed by events occurring in all members of the Saccharum genus. Additionally, evidence of duplicated and truncate gene expression and the behavior of genetic markers in a duplicated region was found. Thus, we describe the complexity involved in sugarcane genetics and genomics and allelic dynamics, which can be useful for understanding the complex polyploid genome.
Continual exit of human cutaneous resident memory CD4 T cells that seed distant skin sites
Authors: Klicznik, M. M.; Morawski, P. A.; Höllbacher, B.; Varkhande, S.; Motley, S.; Rosenblum, M. D.; Long, S. A.; Brachtl, G.; Duhen, T.; Campbell, D.; Gratz, I. K.
Category: immunology
Date: 2018-08-23

Tissue-resident memory T cells (TRM) persist locally in non-lymphoid tissues where they provide front-line defense against recurring insults. TRM at barrier surfaces express the markers CD103 and/or CD69 which function to retain them in epithelial tissues. In humans, neither the long-term migratory behavior of TRM nor their ability to re-enter the circulation and potentially migrate to distant tissue sites have been investigated. Using tissue explant cultures, we found that CD4+CD69+CD103+ TRM in human skin can downregulate CD69 and exit the tissue.

Additionally, we identified a skin-tropic CD4+CD69-CD103+ population in human lymph and blood that is transcriptionally, functionally and clonally related to the CD4+CD69+CD103+ TRM population in the skin. Using a skin xenograft model, we confirmed that a fraction of the human cutaneous CD4+CD103+ TRM population can re-enter circulation, and migrate to secondary human skin sites where they re-assume a TRM phenotype. Thus, our data challenge current concepts regarding the strict tissue compartmentalization of CD4+ T cell memory in humans.

One Sentence SummaryHuman CD4+CD103+ cutaneous resident memory T cells are found in the circulation of healthy subjects, and these cells can seed distant skin sites.
Using the Weibull accelerated failure time regression model to predict time to health events
Authors: Liu, E.; Lim, K.
Category: epidemiology
Date: 2018-08-27

We describe a statistical method protocol to use a Weibull accelerated failure time (AFT) model to predict time to a health-related event. This prediction method is quite common in engineering reliability research but rarely used for medical predictions such as survival time. A worked example for how to perform the prediction using a published dataset is included.
Parallel adaptations to different host plants despite clonal reproduction in the world most devastating nematode pest
Authors: Koutsovoulos, G.; Marques, E.; Arguel, M.-J.; Machado, A. C.; Carneiro, R. M.; Castagnone-Sereno, P.; Danchin, E. G.; Albuquerque, E. V.
Category: evolutionary biology
Date: 2018-08-23

The most devastating nematodes to worldwide agriculture are the root-knot nematodes with Meloidogyne incognita being the most widely distributed and damaging species. This parasitic and ecological success seem surprising given its supposed obligatory clonal reproduction. Clonal reproduction has been suspected based on cytological observations but, so far, never confirmed by population genomics data. At the species level, M. incognita is highly polyphagous with thousands of host plants. However, the host range varies among different M. incognita isolates that may present distinct and more restricted host compatibilities. Historically, four  host races had been defined as a function of ranges of compatible and incompatible plants. We sequenced the genomes of 11 isolates across Brazil, covering these four distinct races to assess (i) how clonal reproduction is and (ii) how the level of genome variability associates with biological traits such as the host races, affected agronomic culture, and geographical distribution. By aligning the genomic reads of the isolates to the M. incognita reference genome assembly, we identified SNV and small-scale insertions/deletions. Analysis of linkage disequilibrium and 4-gametes test, showed no sign of recombination, confirming the clonal mode of reproduction of M. incognita. We showed that there are relatively few point variations between the different isolates, and these variations show no significant association with either the host races, the geographical origin of the samples or the host plant on which they have been collected. Due to the lack of phylogenetic signal underlying their existence, we recommend the discontinuation of the terminology  race. Overall, these results suggest that multiple gains and losses of parasitic abilities and adaptations to different environmental conditions account for the broad host spectrum and wide geographic distribution of M. incognita. Hence, this nematode constitutes a model species to study adaptability without sexual recombination and overall low genomic variations in animals.
Detection of instability in 18S and 5.8S rDNA copy numbers and 18S:5.8S ratios in longitudinally-collected human DNA samples detected by monochrome multiplex qPCR
Authors: Ligerman, C. B.; Baird, L.; Leppert, M. F.; Cawthon, R. M.
Category: genomics
Date: 2018-08-27

Ribosomal DNA encodes the structural RNAs of the ribosomes. Ribosomal DNA instability is a major contributor to aging in yeast, but its role in human aging and longevity is largely unknown. Human 45S rDNA encodes the 18S, 5.8S, and 28S ribosomal RNAs; ranges in copy number from 60 to > 800 per cell, distributed as tandem repeats along the short arms of the five acrocentric chromosomes (p13, p14, p15, p21, and p22); and is prone to frequent homologous and non-homologous recombination. Here we present two multiplex quantitative PCR assays, one for 18S rDNA normalized to the single copy gene beta-globin (HBB), and the other for 5.8S rDNA normalized to the single copy gene albumin (ALB). Longitudinally-collected pairs of DNAs from bloods drawn approximately 16 years apart from 40 females and 39 males, aged < 1 to 77 years at first blood draw, from the Utah CEPH families were assayed. Ribosomal DNA copy number varied over a four-fold range between subjects and increased approximately 14% across the lifespan. Longitudinal within-individual gains in copy number up to +68% and losses down to -25% were observed, while repeated assays of single DNA samples varied approximately +/- 10%. While 18S and 5.8S rDNA tended to be gained and lost together, the 18S:5.8S ratio was also unstable longitudinally, with increases up to +19% and decreases down to -15% observed. The 18S:5.8S ratio at second draw, relative to that at first draw, increased significantly across the lifespan in males, but not in females. To our knowledge this is the first report of within-individual longitudinal changes in the human 18S and 5.8S rDNA copy numbers. These assays will facilitate investigations of the biology of ribosomal RNA genes and their roles in health and disease across the human lifespan.
Instability in 18S and 5.8S rDNA copy numbers and 18S:5.8S ratios in longitudinally-collected human DNA samples detected by monochrome multiplex qPCR
Authors: Ligerman, C. B.; Baird, L.; Leppert, M. F.; Cawthon, R. M.
Category: genomics
Date: 2018-08-28

Ribosomal DNA encodes the structural RNAs of the ribosomes. Ribosomal DNA instability is a major contributor to aging in yeast, but its role in human aging and longevity is largely unknown. Human 45S rDNA encodes the 18S, 5.8S, and 28S ribosomal RNAs; ranges in copy number from 60 to > 800 per cell, distributed as tandem repeats along the short arms of the five acrocentric chromosomes (p13, p14, p15, p21, and p22); and is prone to frequent homologous and non-homologous recombination. Here we present two multiplex quantitative PCR assays, one for 18S rDNA normalized to the single copy gene beta-globin (HBB), and the other for 5.8S rDNA normalized to the single copy gene albumin (ALB). Longitudinally-collected pairs of DNAs from bloods drawn approximately 16 years apart from 40 females and 39 males, aged < 1 to 77 years at first blood draw, from the Utah CEPH families were assayed. Ribosomal DNA copy number varied over a four-fold range between subjects and increased approximately 14% across the lifespan. Longitudinal within-individual gains in copy number up to +68% and losses down to -25% were observed, while repeated assays of single DNA samples varied approximately +/- 10%. While 18S and 5.8S rDNA tended to be gained and lost together, the 18S:5.8S ratio was also unstable longitudinally, with increases up to +19% and decreases down to -15% observed. The 18S:5.8S ratio at second draw, relative to that at first draw, increased significantly across the lifespan in males, but not in females. To our knowledge this is the first report of within-individual longitudinal changes in the human 18S and 5.8S rDNA copy numbers. These assays will facilitate investigations of the biology of ribosomal RNA genes and their roles in health and disease across the human lifespan.
New tools for diet analyses: nanopore sequencing of metagenomic DNA from rat stomach contents to quantify diet
Authors: Pearman, W.; Smith, A. N. H.; Breckell, G.; Dale, J.; Freed, N. E.; Silander, O. K.
Category: ecology
Date: 2018-08-28

BackgroundUsing metagenomics to determine animal diet offers a new and promising alternative to current methods. Here we show that rapid and inexpensive diet quantification is possible through metagenomic sequencing with the portable Oxford Nanopore Technologies (ONT) MinION. Using an amplification-free approach, we profiled the stomach contents from wild-caught rats.

ResultsWe conservatively identified diet items from over 50 taxonomic orders, ranging across nine phyla that include plants, vertebrates, invertebrates, and fungi. This highlights the wide range of taxa that can be identified using this simple approach. We calibrate the accuracy of this method by comparing the characteristics of reads matching the ground-truth host genome (rat) to those matching diet items, and show that at the family-level, false positive taxon assignments are approximately 97.5% accurate. We also suggest a way to mitigate for database biases in metagenomic approaches. Finally, we implement a constrained ordination analysis and show that we can identify the sampling location of an individual rat within tens of kilometres based on diet content alone.

ConclusionsThis work establishes proof-of-principle for long-read metagenomic methods in quantitative diet analysis. We show that diet content can be quantified even with limited expertise, using a simple, amplification free workflow and a relatively inexpensive and accessible next generation sequencing method. Continued increases in the accuracy and throughput of ONT sequencing, along with improved genomic databases, suggests that a metagenomic approach to quantification of animal diets will become an important method in the future.
Systems Biology Understanding of the Effects of Lithium on Cancer.
Authors: Ge, W.; Jakobsson, E.
Category: cancer biology
Date: 2018-08-24

Lithium has many widely varying biochemical and phenomenological effects, suggesting that a systems biology approach is required to understand its action. Multiple lines of evidence point to lithium as a significant factor in development of cancer, showing that understanding lithium action is of high importance. In this paper we undertake first steps towards a systems approach by analyzing mutual enrichment between the interactomes of lithium-sensitive enzymes and the pathways associated with cancer. This work integrates information from two important databases, STRING and KEGG pathways. We find that for the majority of cancer pathways the mutual enrichment is many times greater than chance, reinforcing previous lines of evidence that lithium is an important influence on cancer.
Resource limitation modulates the fate of dissimilated nitrogen in a dual-pathway Actinobacterium
Authors: Vuono, D. C.; Read, R. W.; Hemp, J.; Sullivan, B. W.; Arnone, J. A.; Neveux, I.; Blank, B.; Staub, C.; Loney, E.; Miceli, D.; Winkler, M.; Chakraborty, R.; Stahl, D. A.; Grzymski, J. J.
Category: microbiology
Date: 2018-08-23

Respiratory ammonification and denitrification are two evolutionarily unrelated dissimilatory nitrogen (N) processes central to the global N cycle, the activity of which is thought to be controlled by carbon (C) to nitrate (NO3-) ratio. Here we find that Intrasporangium calvum C5, a novel menaquinone-based dual-pathway denitrifier/respiratory ammonifier, disproportionately utilizes ammonification rather than denitrification when grown under carbon or nitrate limitation, not C:NO3- ratio. Higher growth rates are promoted by ammonification and metabolite and transcriptional profiles during growth show that the bacterium produces its own formate from a fermentable carbon source (lactate) to further generate a proton motive force for the ammonification pathway. Transcript abundances encoding for nitrite reducing enzymes, NrfAH and NirK, also significantly increase in response to nitrite production. Mechanistically, our results suggest that pathway selection is driven by intracellular redox potential (redox poise), which may be lowered during resource limitation, thereby decreasing catalytic activity of upstream electron transport steps needed for denitrification enzymes. Our work advances our understanding of the biogeochemical flexibility of N-cycling organisms and pathway evolution.