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.

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.

The core transcriptome of mammalian placentas and the divergence of expression with placental shape.

INTRODUCTION: The placenta is arguably the most anatomically variable organ in mammals even though its primary function is conserved. METHOD: Using RNA-Seq, we measured the expression profiles of 55 term placentas of 14 species of mammals representing all major eutherian superordinal clades and marsupials, and compared the evolution of expression across clades. RESULTS: We identified a set of 115 core genes which is expressed (FPKM ≥10) in all eutherian placentas, including genes with immune-modulating properties (ANXA2, ANXA1, S100A11, S100A10, and LGALS1), cell-cell interactions (LAMC1, LUM, and LGALS1), invasion (GRB2 and RALB) and syncytialization (ANXA5 and ANXA1). We also identified multiple pre-eclampsia associated genes which are differentially expressed in Homo sapiens when compared to the other 13 species. Multiple genes are significantly associated with placenta morphology, including EREG and WNT5A which are both associated with placental shape. DISCUSSION: 115 genes are important for the core functions of the placenta in all eutherian species analyzed. The molecular functions and pathways enriched in the core placenta align with the evolutionarily conserved functionality of the placenta.

Development of Bioinformatics Infrastructure for Genomics Research.

BACKGROUND: Although pockets of bioinformatics excellence have developed in Africa, generally, large-scale genomic data analysis has been limited by the availability of expertise and infrastructure. H3ABioNet, a pan-African bioinformatics network, was established to build capacity specifically to enable H3Africa (Human Heredity and Health in Africa) researchers to analyze their data in Africa. Since the inception of the H3Africa initiative, H3ABioNet's role has evolved in response to changing needs from the consortium and the African bioinformatics community. OBJECTIVES: H3ABioNet set out to develop core bioinformatics infrastructure and capacity for genomics research in various aspects of data collection, transfer, storage, and analysis. METHODS AND RESULTS: Various resources have been developed to address genomic data management and analysis needs of H3Africa researchers and other scientific communities on the continent. NetMap was developed and used to build an accurate picture of network performance within Africa and between Africa and the rest of the world, and Globus Online has been rolled out to facilitate data transfer. A participant recruitment database was developed to monitor participant enrollment, and data is being harmonized through the use of ontologies and controlled vocabularies. The standardized metadata will be integrated to provide a search facility for H3Africa data and biospecimens. Because H3Africa projects are generating large-scale genomic data, facilities for analysis and interpretation are critical. H3ABioNet is implementing several data analysis platforms that provide a large range of bioinformatics tools or workflows, such as Galaxy, the Job Management System, and eBiokits. A set of reproducible, portable, and cloud-scalable pipelines to support the multiple H3Africa data types are also being developed and dockerized to enable execution on multiple computing infrastructures. In addition, new tools have been developed for analysis of the uniquely divergent African data and for downstream interpretation of prioritized variants. To provide support for these and other bioinformatics queries, an online bioinformatics helpdesk backed by broad consortium expertise has been established. Further support is provided by means of various modes of bioinformatics training. CONCLUSIONS: For the past 4 years, the development of infrastructure support and human capacity through H3ABioNet, have significantly contributed to the establishment of African scientific networks, data analysis facilities, and training programs. Here, we describe the infrastructure and how it has affected genomics and bioinformatics research in Africa.

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.

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.

Thrombospondin-1 modulates the angiogenic phenotype of human cerebral arteriovenous malformation endothelial cells.

BACKGROUND: The management of cerebral arteriovenous malformation (AVM) is challenging, and invasive therapies place vital intracranial structures at risk of injury. The development of noninvasive, pharmacologic approaches relies on identifying factors that mediate key angiogenic processes. Previous studies indicate that endothelial cells (ECs) derived from cerebral AVM (AVM-ECs) are distinct from control brain ECs with regard to important angiogenic characteristics. OBJECTIVE: To determine whether thrombospondin-1 (TSP-1), a potent angiostatic factor, regulates critical angiogenic features of AVM-ECs and to identify factors that modulate TSP-1 production in AVM-ECs. METHODS: EC proliferation, migration, and tubule formation were evaluated with bromodeoxyuridine incorporation, Boyden chamber, and Matrigel studies, respectively. TSP-1 and inhibitor of DNA binding/differentiation 1 (Id1) mRNA levels were quantified with microarray and quantitative real-time polymerase chain reaction analyses. TSP-1 protein expression was measured using Western blotting, immunohistochemical, and enzyme-linked immunosorbent assay techniques. The mechanistic link between Id1 and TSP-1 was established through small interfering RNA-mediated knockdown of Id1 in AVM-ECs followed by Western blot and enzyme-linked immunosorbent assay experiments assessing TSP-1 production. RESULTS: AVM-ECs proliferate faster, migrate more quickly, and form disorganized tubules compared with brain ECs. TSP-1 is significantly down-regulated in AVM-ECs. The addition of TSP-1 to AVM-EC cultures normalizes the rate of proliferation and migration and the efficiency of tubule formation, whereas brain ECs are unaffected. Id1 negatively regulates TSP-1 expression in AVM-ECs. CONCLUSION: These data highlight a novel role for TSP-1 in the pathobiology of AVM angiogenesis and provide a context for its use in the clinical management of brain AVMs.

High-density genotyping of STAT4 reveals multiple haplotypic associations with systemic lupus erythematosus in different racial groups.

OBJECTIVE: Systemic lupus erythematosus (SLE) is the prototypical systemic autoimmune disorder, with complex etiology and a strong genetic component. Recently, gene products involved in the interferon pathway have been under intense investigation in terms of the pathogenesis of SLE. STAT-1 and STAT-4 are transcription factors that play key roles in the interferon and Th1 signaling pathways, making them attractive candidates for involvement in SLE susceptibility. METHODS: Fifty-six single-nucleotide polymorphisms (SNPs) across STAT1 and STAT4 on chromosome 2 were genotyped using the Illumina platform, as part of an extensive association study in a large collection of 9,923 lupus patients and control subjects from different racial groups. DNA samples were obtained from the peripheral blood of patients with SLE and control subjects. Principal components analyses and population-based case-control association analyses were performed, and the P values, false discovery rate q values, and odds ratios with 95% confidence intervals were calculated. RESULTS: We observed strong genetic associations with SLE and multiple SNPs located within STAT4 in different ethnic groups (Fisher's combined P = 7.02 x 10(-25)). In addition to strongly confirming the previously reported association in the third intronic region of this gene, we identified additional haplotypic association across STAT4 and, in particular, a common risk haplotype that is found in multiple racial groups. In contrast, only a relatively weak suggestive association was observed with STAT1, probably due to its proximity to STAT4. CONCLUSION: Our findings indicate that STAT4 is likely to be a crucial component in SLE pathogenesis in multiple racial groups. Knowledge of the functional effects of this association, when they are revealed, might improve our understanding of the disease and provide new therapeutic targets.

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.

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.

Identification of novel susceptibility genes in childhood-onset systemic lupus erythematosus using a uniquely designed candidate gene pathway platform.

OBJECTIVE: Childhood-onset systemic lupus erythematosus (SLE) presents a unique subgroup of patients for genetic study. The present study was undertaken to identify susceptibility genes contributing to SLE, using a novel candidate gene pathway microarray platform to investigate gene expression in patients with childhood-onset SLE and both of their parents. METHODS: Utilizing bioinformatic tools, a platform of 9,412 single-nucleotide polymorphisms (SNPs) from 1,204 genes was designed and validated. Molecular inversion probes and high-throughput SNP technologies were used for assay development. Seven hundred fifty three subjects, corresponding to 251 full trios of childhood-onset SLE families, were genotyped and analyzed using transmission disequilibrium testing (TDT) and multitest corrections. RESULTS: Family-based TDT showed a significant association of SLE with a N673S polymorphism in the P-selectin gene (SELP) (P = 5.74 x 10(-6)) and a C203S polymorphism in the interleukin-1 receptor-associated kinase 1 gene (IRAK1) (P = 9.58 x 10(-6)). These 2 SNPs had a false discovery rate for multitest correction of <0.05, and therefore a >95% probability of being considered as proven. Furthermore, 7 additional SNPs showed q values of <0.5, suggesting association with SLE and providing a direction for followup studies. These additional genes notably included TNFRSF6 (Fas) and IRF5, supporting previous findings of their association with SLE pathogenesis. CONCLUSION: SELP and IRAK1 were identified as novel SLE-associated genes with a high degree of significance, suggesting new directions in understanding the pathogenesis of SLE. The overall design and results of this study demonstrate that the candidate gene pathway microarray platform used provides a novel and powerful approach that is generally applicable in identifying genetic foundations of complex diseases.

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.

Transmembrane domains in the functions of Fc receptors.

In the present study, we use a novel method, PHDhtm, to predict the exact locations and extents of the transmembrane (TM) domains of multisubunit immunoglobulin Fc-receptors. Whereas most previous studies have used single residue hydrophobicity plots for characterizing of these domains, PHDhtm utilizes a system of neural networks and the evolutionary information contained in multiple alignments of related sequences to predict the above. Present PHDhtm application predicts TM domains of immunoglobulin Fc-receptors that in many cases differ significantly from those derived by using earlier methods. Comparisons of helical wheel projections of the presently derived TM domains from PHDhtm with those produced earlier reveal different hydrophobic moments as well as hydrophobic and hydrophilic surfaces. These differences probably alter the character of subunit association within the receptor complexes. This new algorithm can also be used for other membrane protein complexes and may advance both understanding the principles underlying such complexes formation and design of peptides that can interfere with such TM domain association so as to modulate specific cellular responses.

Synergistic perturbation of phosphatidylcholine/sphingomyelin bilayers by diacylglycerol and cholesterol.

The perturbations induced by second messenger diacylglycerols (DAGs) into bovine brain phosphatidylcholine (BBPC) bilayers in the presence or absence of bovine brain sphingomyelin (SM) and/or cholesterol were studied by (2)H NMR. Addition of 15 mol% DAG to BBPC bilayers did not induce non-bilayer lipid phases in the temperature range 30-60 degrees C. Similar measurements performed in the presence of cholesterol revealed that cholesterol progressively destabilizes PC bilayers with respect to DAG-induced perturbations. Thus, at 40 mol% cholesterol, addition of 15 mol% DAG induced the formation of non-bilayer (isotropic and inverted hexagonal) phases at 60 degrees C. Whereas some lateral separation of the bilayers into domains of different cholesterol contents was observed in BBPC/cholesterol membranes, such a lateral heterogeneity was greatly facilitated by the addition of SM. Since both a tendency to form non-bilayer lipid phases and lateral heterogeneity of the membranes are associated with increased activation of a number of membrane-associated enzymes, our results suggest that SM- and cholesterol-enriched regions of biological membranes (rafts) provide an environment with increased sensitivity to the generation of lipid second messengers and modified transmembrane signal transduction properties.