Transcriptomic profiling of fetal membranes of mice deficient in biglycan and decorin as a model of preterm birth†.

Approximately, 25% of all preterm births are due to preterm premature rupture of membranes. Mice deficient in proteoglycans biglycan (Bgn) and decorin (Dcn) display abnormal fetal membranes and increased incidence of preterm birth. We conducted RNA-Seq to profile fetal membranes and identify molecular pathways that may lead to preterm birth in double knockout (DKO) mice (Bgn-/-; Dcn-/-) compared to wild-type (WT) at two different gestational stages, E12 and E18 (n = 3 in each group). 3264 transcripts were differentially regulated in E18 DKO vs. WT fetal membranes, and 96 transcripts differentially regulated in E12 DKO vs. WT fetal membranes (FDR < 0.05, log 2 FC ≥ 1). Differentially regulated transcripts in E18 DKO fetal membranes were significantly enriched for genes involved in cell cycle regulation, extracellular matrix-receptor interaction, and the complement cascade. Fifty transcripts involved in the cell cycle were altered in E18 DKO fetal membranes (40↓, 10↑, FDR < 0.05), including p21 and p57 (↑), and Tgfb2, Smad3, CycA, Cdk1, and Cdk2(↓). Thirty-one transcripts involved in the complement cascade were altered (11↓, 20↑, FDR < 0.05) in E18 DKO fetal membranes, including C1q, C2, and C3 (↑). Differentially expressed genes in the top three molecular pathways (1) showed evidence of negative or purifying selection, and (2) were significantly enriched (Z-score > 10) for transcription factor binding sites for Nr2f1 at E18. We propose that in DKO mice, cell cycle arrest results in lack of cell proliferation in fetal membranes, inability to contain the growing fetus, and preterm birth.

Maternal weight affects placental DNA methylation of genes involved in metabolic pathways in the common marmoset monkey (Callithrix jacchus).

Accumulating evidence suggests that dysregulation of placental DNA methylation (DNAm) is a mechanism linking maternal weight during pregnancy to metabolic programming outcomes. The common marmoset, Callithrix jaccus, is a platyrrhine primate species that has provided much insight into studies of the primate placenta, maternal condition, and metabolic programming, yet the relationships between maternal weight and placental DNAm are unknown. Here, we report genome-wide DNAm from term marmoset placentas using reduced representation bisulfite sequencing. We identified 74 genes whose DNAm pattern is associated with maternal weight during gestation. These genes are predominantly involved in energy metabolism and homeostasis, including the regulation of glycolytic and lipid metabolic processes pathways.

Methylomic profiles reveal sex-specific differences in leukocyte composition associated with post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) is a debilitating mental disorder precipitated by trauma exposure. However, only some persons exposed to trauma develop PTSD. There are sex differences in risk; twice as many women as men develop a lifetime diagnosis of PTSD. Methylomic profiles derived from peripheral blood are well-suited for investigating PTSD because DNA methylation (DNAm) encodes individual response to trauma and may play a key role in the immune dysregulation characteristic of PTSD pathophysiology. In the current study, we leveraged recent methodological advances to investigate sex-specific differences in DNAm-based leukocyte composition that are associated with lifetime PTSD. We estimated leukocyte composition on a combined methylation array dataset (483 participants, ∼450 k CpG sites) consisting of two civilian cohorts, the Detroit Neighborhood Health Study and Grady Trauma Project. Sex-stratified Mann-Whitney U test and two-way ANCOVA revealed that lifetime PTSD was associated with significantly higher monocyte proportions in males, but not in females (Holm-adjusted p-val < 0.05). No difference in monocyte proportions was observed between current and remitted PTSD cases in males, suggesting that this sex-specific difference may reflect a long-standing trait of lifetime history of PTSD, rather than current state of PTSD. Associations with lifetime PTSD or PTSD status were not observed in any other leukocyte subtype and our finding in monocytes was confirmed using cell estimates based on a different deconvolution algorithm, suggesting that our sex-specific findings are robust across cell estimation approaches. Overall, our main finding of elevated monocyte proportions in males, but not in females with lifetime history of PTSD provides evidence for a sex-specific difference in peripheral blood leukocyte composition that is detectable in methylomic profiles and that may reflect long-standing changes associated with PTSD diagnosis.

Developing reproducible bioinformatics analysis workflows for heterogeneous computing environments to support African genomics.

BACKGROUND: The Pan-African bioinformatics network, H3ABioNet, comprises 27 research institutions in 17 African countries. H3ABioNet is part of the Human Health and Heredity in Africa program (H3Africa), an African-led research consortium funded by the US National Institutes of Health and the UK Wellcome Trust, aimed at using genomics to study and improve the health of Africans. A key role of H3ABioNet is to support H3Africa projects by building bioinformatics infrastructure such as portable and reproducible bioinformatics workflows for use on heterogeneous African computing environments. Processing and analysis of genomic data is an example of a big data application requiring complex interdependent data analysis workflows. Such bioinformatics workflows take the primary and secondary input data through several computationally-intensive processing steps using different software packages, where some of the outputs form inputs for other steps. Implementing scalable, reproducible, portable and easy-to-use workflows is particularly challenging. RESULTS: H3ABioNet has built four workflows to support (1) the calling of variants from high-throughput sequencing data; (2) the analysis of microbial populations from 16S rDNA sequence data; (3) genotyping and genome-wide association studies; and (4) single nucleotide polymorphism imputation. A week-long hackathon was organized in August 2016 with participants from six African bioinformatics groups, and US and European collaborators. Two of the workflows are built using the Common Workflow Language framework (CWL) and two using Nextflow. All the workflows are containerized for improved portability and reproducibility using Docker, and are publicly available for use by members of the H3Africa consortium and the international research community. CONCLUSION: The H3ABioNet workflows have been implemented in view of offering ease of use for the end user and high levels of reproducibility and portability, all while following modern state of the art bioinformatics data processing protocols. The H3ABioNet workflows will service the H3Africa consortium projects and are currently in use. All four workflows are also publicly available for research scientists worldwide to use and adapt for their respective needs. The H3ABioNet workflows will help develop bioinformatics capacity and assist genomics research within Africa and serve to increase the scientific output of H3Africa and its Pan-African Bioinformatics Network.

Systemic lupus erythematosus-associated neutrophil cytosolic factor 2 mutation affects the structure of NADPH oxidase complex.

In a case-control association study with 3716 North Americans of Hispanic descent and 4867 North Americans of European descent, we show that the associations of rs17849502 (NCF2 His-389 → Gln) and rs13306575 (NCF2 Arg-395 → Trp) with systemic lupus erythematosus are independent. We have shown that His-389 → Gln disrupts the binding of NCF2 to the ZF domain of VAV1, resulting in decreased NADPH oxidase activity. With respect to Arg-395 → Trp, using protein docking and structure analyses, we provide a model for the involvement of this mutation in the structure and function of the NADPH oxidase complex. This model assigns a central role to Arg-395 in the structure and stability of the quaternary NCF2/NCF4/VAV1/RAC1 NADPH oxidase complex. Arg-395 stabilizes the C-terminal tail of NCF4 and the conformation of NCF2 loop 395-402, which in turn stabilize the evolutionarily conserved interactions of NCF2/NCF4 with the DH domain of VAV1 and RAC1 region 120-137. Our findings are consistent with the high levels of conservation of all of the residues involved in these interactions.

Lupus-associated causal mutation in neutrophil cytosolic factor 2 (NCF2) brings unique insights to the structure and function of NADPH oxidase.

Systemic lupus erythematosus (SLE), the prototypic systemic autoimmune disease, is a debilitating multisystem autoimmune disorder characterized by chronic inflammation and extensive immune dysregulation in multiple organ systems, resulting in significant morbidity and mortality. Here, we present a multidisciplinary approach resulting in the identification of neutrophil cytosolic factor 2 (NCF2) as an important risk factor for SLE and the detailed characterization of its causal variant. We show that NCF2 is strongly associated with increased SLE risk in two independent populations: childhood-onset SLE and adult-onset SLE. The association between NCF2 and SLE can be attributed to a single nonsynonymous coding mutation in exon 12, the effect of which is the substitution of histidine-389 with glutamine (H389Q) in the PB1 domain of the NCF2 protein, with glutamine being the risk allele. Computational modeling suggests that the NCF2 H389Q mutation reduces the binding efficiency of NCF2 with the guanine nucleotide exchange factor Vav1. The model predicts that NCF2/H389 residue interacts with Vav1 residues E509, N510, E556, and G559 in the ZF domain of Vav1. Furthermore, replacing H389 with Q results in 1.5 kcal/mol weaker binding. To examine the effect of the NCF2 H389Q mutation on NADPH oxidase function, site-specific mutations at the 389 position in NCF2 were tested. Results show that an H389Q mutation causes a twofold decrease in reactive oxygen species production induced by the activation of the Vav-dependent Fcγ receptor-elicited NADPH oxidase activity. Our study completes the chain of evidence from genetic association to specific molecular function.

Role of the MEOX2 homeobox gene in neurovascular dysfunction in Alzheimer disease.

Neurovascular dysfunction substantially contributes to Alzheimer disease. Here, we show that transcriptional profiling of human brain endothelial cells (BECs) defines a subset of genes whose expression is age-independent but is considerably altered in Alzheimer disease, including the homeobox gene MEOX2 (also known as GAX), a regulator of vascular differentiation, whose expression is low in Alzheimer disease. By using viral-mediated MEOX2 gene silencing and transfer, we show that restoring expression of the protein it encodes, GAX, in BECs from individuals with Alzheimer disease stimulates angiogenesis, transcriptionally suppresses AFX1 forkhead transcription factor-mediated apoptosis and increases the levels of a major amyloid-beta peptide (Abeta) clearance receptor, the low-density lipoprotein receptor-related protein 1 (LRP), at the blood-brain barrier. In mice, deletion of Meox2 (also known as Gax) results in reductions in brain capillary density and resting cerebral blood flow, loss of the angiogenic response to hypoxia in the brain and an impaired Abeta efflux from brain caused by reduced LRP levels. The link of MEOX2 to neurovascular dysfunction in Alzheimer disease provides new mechanistic and therapeutic insights into this illness.

Identification of IRAK1 as a risk gene with critical role in the pathogenesis of systemic lupus erythematosus.

A combined forward and reverse genetic approach was undertaken to test the candidacy of IRAK1 (interleukin-1 receptor associated kinase-1) as an X chromosome-encoded risk factor for systemic lupus erythematosus (SLE). In studying approximately 5,000 subjects and healthy controls, 5 SNPs spanning the IRAK1 gene showed disease association (P values reaching 10(-10), odds ratio >1.5) in both adult- and childhood-onset SLE, in 4 different ethnic groups, with a 4 SNP haplotype (GGGG) being strongly associated with the disease. The functional role of IRAK1 was next examined by using congenic mouse models bearing the disease loci: Sle1 or Sle3. IRAK1 deficiency abrogated all lupus-associated phenotypes, including IgM and IgG autoantibodies, lymphocytic activation, and renal disease in both models. In addition, the absence of IRAK1 reversed the dendritic cell "hyperactivity" associated with Sle3. Collectively, the forward genetic studies in human SLE and the mechanistic studies in mouse models establish IRAK1 as a disease gene in lupus, capable of modulating at least 2 key checkpoints in disease development. This demonstration of an X chromosome gene as a disease susceptibility factor in human SLE raises the possibility that the gender difference in SLE may in part be attributed to sex chromosome genes.

Replication of simulated prebiotic amphiphile vesicles controlled by experimental lipid physicochemical properties.

We present a new embodiment of the graded autocatalysis replication domain (GARD) for the growth, replication and evolution of lipid vesicles based on a semi-empirical foundation using experimentally measured kinetic values of selected extant lipid species. Extensive simulations using this formalism elucidated the details of the dependence of the replication and properties of the vesicles on the physicochemical properties and concentrations of the lipids, both in the environment and in the vesicle. As expected, the overall concentration and number of amphiphilic components strongly affect average replication time. Furthermore, variations in acyl chain length and unsaturation of vesicles also influence replication rate, as do the relative concentrations of individual lipid types. Understanding of the dependence of replication rates on physicochemical parameters opens a new direction in the study of prebiotic vesicles and lays the groundwork for future studies involving the competition between lipid vesicles for available amphiphilic monomers.

RAGE mediates amyloid-beta peptide transport across the blood-brain barrier and accumulation in brain.

Amyloid-beta peptide (Abeta) interacts with the vasculature to influence Abeta levels in the brain and cerebral blood flow, providing a means of amplifying the Abeta-induced cellular stress underlying neuronal dysfunction and dementia. Systemic Abeta infusion and studies in genetically manipulated mice show that Abeta interaction with receptor for advanced glycation end products (RAGE)-bearing cells in the vessel wall results in transport of Abeta across the blood-brain barrier (BBB) and expression of proinflammatory cytokines and endothelin-1 (ET-1), the latter mediating Abeta-induced vasoconstriction. Inhibition of RAGE-ligand interaction suppresses accumulation of Abeta in brain parenchyma in a mouse transgenic model. These findings suggest that vascular RAGE is a target for inhibiting pathogenic consequences of Abeta-vascular interactions, including development of cerebral amyloidosis.

Associations between childhood family emotional health, fronto-limbic grey matter volume, and saliva 5mC in young adulthood.

BACKGROUND: Poor family emotional health (FEH) during childhood is prevalent and impactful, and likely confers similar neurodevelopmental risks as other adverse social environments. Pointed FEH study efforts are underdeveloped, and the mechanisms by which poor FEH are biologically embedded are unclear. The current exploratory study examined whether variability in 5-methyl-cytosine (5mC) and fronto-limbic grey matter volume may represent pathways through which FEH may become biologically embedded. RESULTS: In 98 university students aged 18-22 years, retrospective self-reported childhood FEH was associated with right hemisphere hippocampus (b = 10.4, p = 0.005), left hemisphere amygdala (b = 5.3, p = 0.009), and right hemisphere amygdala (b = 5.8, p = 0.016) volumes. After pre-processing and filtering to 5mC probes correlated between saliva and brain, analyses showed that childhood FEH was associated with 49 5mC principal components (module eigengenes; MEs) (prange = 3 × 10-6 to 0.047). Saliva-derived 5mC MEs partially mediated the association between FEH and right hippocampal volume (Burlywood ME indirect effect b = - 111, p = 0.014), and fully mediated the FEH and right amygdala volume relationship (Pink4 ME indirect effect b = - 48, p = 0.026). Modules were enriched with probes falling in genes with immune, central nervous system (CNS), cellular development/differentiation, and metabolic functions. CONCLUSIONS: Findings extend work highlighting neurodevelopmental variability associated with adverse social environment exposure during childhood by specifically implicating poor FEH, while informing a mechanism of biological embedding. FEH-associated epigenetic signatures could function as proxies of altered fronto-limbic grey matter volume associated with poor childhood FEH and inform further investigation into primarily affected tissues such as endocrine, immune, and CNS cell types.

RANDOMIZE: A Web Server for Data Randomization.

The microarray-based Illumina Infinium MethylationEpic BeadChip (Epic 850k) has become a useful and standard tool for epigenome wide deoxyribonucleic acid (DNA) methylation profiling. Data from this technology may suffer from batch effects due to improper handling of the samples during the plating process. Batch effects are a significant issue and can give rise to spurious and inaccurate results and reduction in power to detect real biological differences. Careful study design, such as randomizing the samples to uniformly distribute the samples across the factors responsible for batch effects, is crucial to address batch effects and other technical artifacts. Randomization helps to reduce the likelihood of bias and impact of difference among groups. This process of randomizing the samples can be a tedious, error-prone, and time-consuming task without a user-friendly and efficient tool. We present RANDOMIZE, a web-based application designed to perform randomization of relevant metadata to evenly distribute samples across the factors typically responsible for batch effects in DNA methylation microarrays, such as rows, chips and plates. We demonstrate that the tool is efficient, fast and easy to use. The tool is freely available online at https://coph-usf.shinyapps.io/RANDOMIZE/ and can be accessed using any web browser. Sample data and tutorial is also available with the tool.

Early Experiences of Threat, but Not Deprivation, Are Associated With Accelerated Biological Aging in Children and Adolescents.

BACKGROUND: Recent conceptual models argue that early life adversity (ELA) accelerates development, which may contribute to poor mental and physical health outcomes. Evidence for accelerated development in youths comes from studies of telomere shortening or advanced pubertal development following circumscribed ELA experiences and neuroimaging studies of circuits involved in emotional processing. It is unclear whether all ELA is associated with accelerated development across global metrics of biological aging or whether this pattern emerges following specific adversity types. METHODS: In 247 children and adolescents 8 to 16 years of age with wide variability in ELA exposure, we evaluated the hypothesis that early environments characterized by threat, but not deprivation, would be associated with accelerated development across two global biological aging metrics: DNA methylation (DNAm) age and pubertal stage relative to chronological age. We also examined whether accelerated development explained associations of ELA with depressive symptoms and externalizing problems. RESULTS: Exposure to threat-related ELA (e.g., violence) was associated with accelerated DNAm age and advanced pubertal stage, but exposure to deprivation (e.g., neglect, food insecurity) was not. In models including both ELA types, threat-related ELA was uniquely associated with accelerated DNAm age (ß = .18) and advanced pubertal stage (ß = .28), whereas deprivation was uniquely associated with delayed pubertal stage (ß = -.21). Older DNAm age was related to greater depressive symptoms, and a significant indirect effect of threat exposure on depressive symptoms was observed through DNAm age. CONCLUSIONS: Early threat-related experiences are particularly associated with accelerated biological aging in youths, which may be a mechanism linking ELA with depressive symptoms.

Amplification of diacylglycerol activation of protein kinase C by cholesterol.

The combined effects of cholesterol, a major cell membrane component, and the lipid second messenger diacylglycerol on the activity of protein kinase C (PK-C) and the structure of phosphatidylcholine/phosphatidylserine bilayers were investigated using specific PK-C assays and (2)H NMR. Whereas the classical activation of PK-C was observed as an effect of diacylglycerol, in the absence of this second messenger, cholesterol did not affect PK-C activity. A novel effect of amplified PK-C activation was observed in the presence of both cholesterol and diacylglycerol concentrations within the physiological range of each of these components. (2)H NMR results suggest that this phenomenon is due to cholesterol- and diacylglycerol-induced increased propensity of the lipids to adopt nonbilayer phases, effectively destabilizing the bilayer structure. The magnitude of the effect was a function of cholesterol concentration, implying that laterally separated cell membrane domains with distinct cholesterol concentrations have the capacity to differ in their sensitivity to extracellular stimuli.

Function2Gene: a gene selection tool to increase the power of genetic association studies by utilizing public databases and expert knowledge.

BACKGROUND: Many common disorders have multiple genetic components which convey increased susceptibility. SNPs have been used to identify genetic components which are associated with a disease. Unfortunately, many studies using these methods suffer from low reproducibility due to lack of power. RESULTS: We present a set of programs which implement a novel method for searching for disease-associated genes using prior information to select and order genes from publicly available databases by their prior likelihood of association with the disease. These programs were used in a published study of childhood-onset SLE which yielded novel associations with modest sample size. CONCLUSION: Using prior information to decrease the size of the problem space to an amount commensurate with available samples and resources while maintaining appropriate power enables researchers to increase their likelihood of discovering reproducible associations.

Extracellular Vesicles and the Promise of Continuous Liquid Biopsies.

The rapid and accurate diagnosis of patients with minimally invasive procedures was once only found in science fiction. However, the discovery of extracellular vesicles (EVs) and their near ubiquity in body fluids, coupled with the advent of inexpensive next generation sequencing techniques and EV purification protocols, promises to make science fiction a reality. Purifying and sequencing the RNA content of EV from routine blood draws and urine samples are likely to enable pathologists and physicians to diagnose and track the progress of diseases in many inaccessible tissues in the near future. Here we present the evolutionary background of EV, summarize the biology of EV formation and cargo selection, and discuss the current barriers to making continuous liquid biopsies through the use of EV a science reality.

Replication of Simulated Prebiotic Amphiphilic Vesicles in a Finite Environment Exhibits Complex Behavior That Includes High Progeny Variability and Competition.

We studied the simulated replication and growth of prebiotic vesicles composed of 140 phospholipids and cholesterol using our R-GARD (Real Graded Autocatalysis Replication Domain) formalism that utilizes currently extant lipids that have known rate constants of lipid-vesicle interactions from published experimental data. R-GARD normally modifies kinetic parameters of lipid-vesicle interactions based on vesicle composition and properties. Our original R-GARD model tracked the growth and division of one vesicle at a time in an environment with unlimited lipids at a constant concentration. We explore here a modified model where vesicles compete for a finite supply of lipids. We observed that vesicles exhibit complex behavior including initial fast unrestricted growth, followed by intervesicle competition for diminishing resources, then a second growth burst driven by better-adapted vesicles, and ending with a final steady state. Furthermore, in simulations without kinetic parameter modifications ("invariant kinetics"), the initial replication was an order of magnitude slower, and vesicles' composition variability at the final steady state was much lower. The complex kinetic behavior was not observed either in the previously published R-GARD simulations or in additional simulations presented here with only one lipid component. This demonstrates that both a finite environment (inducing selection) and multiple components (providing variation for selection to act upon) are crucial for portraying evolution-like behavior. Such properties can improve survival in a changing environment by increasing the ability of early protocellular entities to respond to rapid environmental fluctuations likely present during abiogenesis both on Earth and possibly on other planets. This in silico simulation predicts that a relatively simple in vitro chemical system containing only lipid molecules might exhibit properties that are relevant to prebiotic processes. Key Words: Phospholipid vesicles-Prebiotic compartments-Prebiotic vesicle competition-Prebiotic vesicle variability. Astrobiology 18, 419-430.

Epigenetic profiles associated with major depression in the human brain.

We conducted an epigenome-wide association study of Major Depressive Disorder (MDD) in brain-derived DNA using two analytic approaches. DNA methylation data (GSE41826) was used in differential methylation (DM) analyses controlling for age, sex, suicide status, and post-mortem interval; and in weighted gene co-methylation network analyses (WGCNA) in probes mapping to transcription start sites. No probes in the DM analysis survived FDR correction. Nominally significant DM probes were enriched in synaptic function-related genes. WGCNA revealed one module correlated with MDD, enriched in genes associated with mitochondrial function. DM and WGCNA both showed enrichment of genes involved in transcription and DNA binding.

Epigenetic meta-analysis across three civilian cohorts identifies NRG1 and HGS as blood-based biomarkers for post-traumatic stress disorder.

AIM: Trauma exposure is a necessary, but not deterministic, contributor to post-traumatic stress disorder (PTSD). Epigenetic factors may distinguish between trauma-exposed individuals with versus without PTSD. MATERIALS & METHODS: We conducted a meta-analysis of PTSD epigenome-wide association studies in trauma-exposed cohorts drawn from civilian contexts. Whole blood-derived DNA methylation levels were analyzed in 545 study participants, drawn from the three civilian cohorts participating in the PTSD working group of the Psychiatric Genomics Consortium. RESULTS: Two CpG sites significantly associated with current PTSD in NRG1 (cg23637605) and in HGS (cg19577098). CONCLUSION: PTSD is associated with differential methylation, measured in blood, within HGS and NRG1 across three civilian cohorts.

Organizing and running bioinformatics hackathons within Africa: The H3ABioNet cloud computing experience.

The need for portable and reproducible genomics analysis pipelines is growing globally as well as in Africa, especially with the growth of collaborative projects like the Human Health and Heredity in Africa Consortium (H3Africa). The Pan-African H3Africa Bioinformatics Network (H3ABioNet) recognized the need for portable, reproducible pipelines adapted to heterogeneous compute environments, and for the nurturing of technical expertise in workflow languages and containerization technologies. To address this need, in 2016 H3ABioNet arranged its first Cloud Computing and Reproducible Workflows Hackathon, with the purpose of building key genomics analysis pipelines able to run on heterogeneous computing environments and meeting the needs of H3Africa research projects. This paper describes the preparations for this hackathon and reflects upon the lessons learned about its impact on building the technical and scientific expertise of African researchers. The workflows developed were made publicly available in GitHub repositories and deposited as container images on quay.io.