Genome-wide association study of early-onset and late-onset postpartum depression: the IGEDEPP prospective study

Postpartum depression (PPD) appears at two peak periods: early-onset prior to 2 months after delivery and late-onset (2 months after delivery and beyond). The aim of our study is to evaluate the different genetic factors associated with early- and late-onset PPD. With the French multicenter interaction of gene and environment of depression during postpartum (IGEDEPP) cohort, we conducted a genome-wide association study (GWAS) on 234 women with early-onset PPD and 223 women with late-onset PPD, as well as 1,204 controls with no history of lifetime depression. We performed post-GWAS analyses: functional mapping and annotation of GWAS results using MAGMA thanks to Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA), expression quantitative trait loci (QTL) analyses, mapping using data from the PsychENCODE and GTEx, and polygenic risk score (PRS) analysis based on published GWAS. We found two new significant candidate loci for early-onset PPD, rs6436132 in PTPRN gene on chromosome 2 and rs184644645 in RAD18 on chromosome 14, respectively, and one region of interest with five significant associated SNPs in chromosome 20 for late-onset PPD. Variant rs6436132 is the most significant associated with early-onset PPD, and it is a QTL that significantly modifies the expression and splicing of the PTPRN gene in different brain tissues. We also found an enrichment of uterus tissue in the early expression of PPD genes. PRS analysis showed a genetic overlap between both early and late-onset PPD and major depressive disorder, but only early-onset PPD overlaps with bipolar disorder. Our study presents two GWAS separately, highlighting two candidate loci for early-onset PPD and one different region of interest for late-onset PPD. These results have important consequences in our understanding of these disorders, especially since our data reinforce the hormonal pathophysiological hypotheses for early-onset PPD.

Meta-Analyses of Genome-Wide Association Studies for Postpartum Depression

OBJECTIVE

Postpartum depression (PPD) is a common subtype of major depressive disorder (MDD) that is more heritable, yet is understudied in psychiatric genetics. The authors conducted meta-analyses of genome-wide association studies (GWASs) to investigate the genetic architecture of PPD.

METHOD

Meta-analyses were conducted on 18 cohorts of European ancestry (17,339 PPD cases and 53,426 controls), one cohort of East Asian ancestry (975 cases and 3,780 controls), and one cohort of African ancestry (456 cases and 1,255 controls), totaling 18,770 PPD cases and 58,461 controls. Post-GWAS analyses included 1) single-nucleotide polymorphism (SNP)-based heritability ([Formula: see text]), 2) genetic correlations between PPD and other phenotypes, and 3) enrichment of the PPD GWAS findings in 27 human tissues and 265 cell types from the mouse central and peripheral nervous system.

RESULTS

No SNP achieved genome-wide significance in the European or the trans-ancestry meta-analyses. The [Formula: see text] of PPD was 0.14 (SE=0.02). Significant genetic correlations were estimated for PPD with MDD, bipolar disorder, anxiety disorders, posttraumatic stress disorder, insomnia, age at menarche, and polycystic ovary syndrome. Cell-type enrichment analyses implicate inhibitory neurons in the thalamus and cholinergic neurons within septal nuclei of the hypothalamus, a pattern that differs from MDD.

CONCLUSIONS

While more samples are needed to reach genome-wide levels of significance, the results presented confirm PPD as a polygenic and heritable phenotype. There is also evidence that despite a high correlation with MDD, PPD may have unique genetic components. Cell enrichment results suggest GABAergic neurons, which converge on a common mechanism with the only medication approved by the U.S. Food and Drug Administration for PPD (brexanolone).

Early control of viral load by favipiravir promotes survival to Ebola virus challenge and prevents cytokine storm in non-human primates

Ebola virus has been responsible for two major epidemics over the last several years and there has been a strong effort to find potential treatments that can improve the disease outcome. Antiviral favipiravir was thus tested on non-human primates infected with Ebola virus. Half of the treated animals survived the Ebola virus challenge, whereas the infection was fully lethal for the untreated ones. Moreover, the treated animals that did not survive died later than the controls. We evaluated the hematological, virological, biochemical, and immunological parameters of the animals and performed proteomic analysis at various timepoints of the disease. The viral load strongly correlated with dysregulation of the biological functions involved in pathogenesis, notably the inflammatory response, hemostatic functions, and response to stress. Thus, the management of viral replication in Ebola virus disease is of crucial importance in preventing the immunopathogenic disorders and septic-like shock syndrome generally observed in Ebola virus-infected patients.

Ebola virus was responsible for several epidemics in the recent years and is now considered as a major public health concern in Central and West African countries. We and others demonstrated that pathogenic events observed during Ebola virus disease are linked to a deleterious immune response. However, the mechanisms implicated are not fully understood. Here, we studied immune responses depending on the viral loads observed in infected cynomolgus monkeys. An antiviral treatment allowed the reduction of viral load in some animals and we observed that these animals did not experience deleterious immune response and the loss of hemostasis. The release of pathogen-associated molecular patterns may thus be limited by the inhibition of viral replication, avoiding the overstimulation of the immune system and consequently the pathogenic events observed in Ebola virus disease.

General regression model: A "model-free" association test for quantitative traits allowing to test for the underlying genetic model

Most genome-wide association studies used genetic-model-based tests assuming an additive mode of inheritance, leading to underpowered association tests in case of departure from additivity. The general regression model (GRM) association test proposed by Fisher and Wilson in 1980 makes no assumption on the genetic model. Interestingly, it also allows formal testing of the underlying genetic model. We conducted a simulation study of quantitative traits to compare the power of the GRM test to the classical linear regression tests, the maximum of the three statistics (MAX), and the allele-based (allelic) tests. Simulations were performed on two samples sizes, using a large panel of genetic models, varying genetic models, minor allele frequencies, and the percentage of explained variance. In case of departure from additivity, the GRM was more powerful than the additive regression tests (power gain reaching 80%) and had similar power when the true model is additive. GRM was also as or more powerful than the MAX or allelic tests. The true simulated model was mostly retained by the GRM test. Application of GRM to HbA1c illustrates its gain in power. To conclude, GRM increases power to detect association for quantitative traits, allows determining the genetic model and is easily applicable.

Immune parameters and outcomes during Ebola virus disease

BACKGROUND

The West African Ebola virus epidemic from 2014-2016 highlighted the lack of knowledge about the pathogenicity of the virus and the factors responsible for outcome. A performant and rapid diagnosis is of crucial importance, as is overcoming the difficulty of providing high-quality patient management during such an extensive outbreak. Here, we propose to study the role of the immune mediators during Ebola virus disease and to define some molecules of importance in the outcome.

METHODS

Plasma from Guinean patients sampled during the outbreak were analyzed using RT-qPCR, magnetic bead assay, ELISA, and high-quality statistical analyses. We also performed a transcriptomic analysis in leukocytes samples. Therefore, we deeply characterized the immune responses involved in Ebola virus disease.

RESULTS

We evaluated the immune patterns depending on the outcome of the disease. Survivors presented an efficient and well-balanced immune response, whereas fatalities were characterized by an intense inflammatory response, overexpression of multiple cytokines, and a "chemokine storm." The plasma concentration of most of the parameters tested increased until death. Statistical analyses also allowed us to define a panel of markers highly predictive of outcome.

CONCLUSION

The immune response observed in fatalities was highly similar to that characterizing septic shock syndrome. Our results suggest that immune responses can play a major pathogenic role during severe Ebola virus infection and argue in favor of therapeutic approaches that act on both viral replication and the induction of shock syndrome.

FUNDING

French Ministry of Foreign Affairs, the Agence Française de Développement, and the Institut Pasteur.

A systematic review of non-genetic predictors and genetic factors of glycated haemoglobin in type 1 diabetes one year after diagnosis

Type 1 diabetes (T1D) results from autoimmune destruction of the pancreatic βcells. Although all T1D patients require daily administration of exogenous insulin, their insulin requirement to achieve good glycaemic control may vary significantly. Glycated haemoglobin (HbA1c) level represents a stable indicator of glycaemic control and is a reliable predictor of long-term complications of T1D. The purpose of this article is to systematically review the role of non-genetic predictors and genetic factors of HbA1c level in T1D patients after the first year of T1D, to exclude the honeymoon period. A total of 1974 articles published since January 2011 were identified and 78 were finally included in the analysis of non-genetic predictors. For genetic factors, a total of 277 articles were identified and 14 were included. The most significantly associated factors with HbA1c level are demographic (age, ethnicity, and socioeconomic status), personal (family characteristics, parental care, psychological traits...) and features related to T1D (duration of T1D, adherence to treatment …). Only a few studies have searched for genetic factors influencing HbA1c level, most of which focused on candidate genes using classical genetic statistical methods, with generally limited power and incomplete adjustment for confounding factors and multiple testing. Our review shows the complexity of explaining HbA1c level variations, which involves numerous correlated predictors. Overall, our review underlines the lack of studies investigating jointly genetic and non-genetic factors and their interactions to better understand factors influencing glycaemic control for T1D patients.

On the turn-inducing properties of asparagine: the structuring role of the amide side chain, from isolated model peptides to crystallized proteins

The anchoring properties of an asparagine (Asn) residue to its local backbone environment in turn model peptides is characterized using gas phase laser spectroscopy and compared to crystallized protein structures.

Reactive and Inelastic Channels in the Ca*···FCH3 Transition State: A Simple Branching Mechanism

To study the excited state dynamics between a calcium atom and the CH3F molecule, a Ca···CH3F 1:1 complex has been prepared by a supersonic expansion with laser ablation of calcium metal in the gas phase. Tunable laser excitation of these complexes in molecular states correlating to Ca (1)P1(4s4p) + CH3F allows observing two competitive channels: the direct dissociation and the reactive channel into CaF* + CH3. The translational recoil, as well as the alignment of the fragments Ca* and CaF* have been analyzed by velocity map imaging and time-of-flight mass spectrometry. This revealed that both the dissociation and reaction processes are quasi direct and are of comparable intensity. We provide a simple interpretation for this process: the electronically excited potential surface of the Ca*···FCH3 complex initiates a fast predissociation from a suspended well to two repulsive surfaces that lead either to Ca (1)P1(4s4p) (Ω = 1) + CH3F or to CaF((2)Δ) + CH3.

Unraveling non-covalent interactions within flexible biomolecules: from electron density topology to gas phase spectroscopy

The NCI (Non-Covalent Interactions) method, a recently-developed theoretical strategy to visualize weak non-covalent interactions from the topological analysis of the electron density and of its reduced gradient, is applied in the present paper to document intra- and inter-molecular interactions in flexible molecules and systems of biological interest in combination with IR spectroscopy. We first describe the conditions of application of the NCI method to the specific case of intramolecular interactions. Then we apply it to a series of stable conformations of isolated molecules as an interpretative technique to decipher the different physical interactions at play in these systems. Examples are chosen among neutral molecular systems exhibiting a large diversity of interactions, for which an extensive spectroscopic characterization under gas-phase isolation conditions has been obtained using state-of-the-art conformer-specific experimental techniques. The interactions presently documented range from weak intra-molecular H-bonds in simple amino-alcohols, to more complex patterns, with interactions of various strengths in model peptides, as well as in chiral bimolecular systems, where invaluable hints for the understanding of chiral recognition are revealed. We also provide a detailed technical appendix, which discusses the choices of cut-offs as well as the applicability of the NCI analysis to specific constrained systems, where local effects require attention. Finally, the NCI technique provides IR spectroscopists with an elegant visualization of the interactions that potentially impact their vibrational probes, namely the OH and NH stretching motions. This contribution illustrates the power and the conditions of use of the NCI technique, with the aim of providing an easy tool for all chemists, experimentalists and theoreticians, for the visualization and characterization of the interactions shaping complex molecular systems.

Gas-phase folding of a two-residue model peptide chain: on the importance of an interplay between experiment and theory

In order to assess the ability of theory to describe properly the dispersive interactions that are ubiquitous in peptide and protein systems, an isolated short peptide chain has been studied using both gas-phase laser spectroscopy and quantum chemistry. The experimentally observed coexistence of an extended form and a folded form in the supersonic expansion was found to result from comparable Gibbs free energies for the two species under the high-temperature conditions (< or = 320 K) resulting from the laser desorption technique used to vaporize the molecules. These data have been compared to results obtained using a series of quantum chemistry methods, including DFT, DFT-D, and post-Hartree-Fock methods, which give rise to a wide range of relative stabilities predicted for these two forms. The experimental observation was best reproduced by an empirically dispersion-corrected functional (B97-D) and a hybrid functional with a significant Hartree-Fock exchange term (M06-2X). In contrast, the popular post-Hartree-Fock method MP2, which is often used for benchmarking these systems, had to be discarded because of a very large basis-set superposition error. The applicability of the atomic counterpoise correction (ACP) is also discussed. This work also introduces the mandatory theoretical examination of experimental abundances. DeltaH(0 K) predictions are clearly not sufficient for discussion of folding, as the conformation inversion temperature is crucial to the conformation determination and requires taking into account thermodynamical corrections (DeltaG) in order to computationally isolate the most stable conformation.

Transition-state spectroscopy of the photoinduced Ca + CH3F reaction. 3. Reaction following the local excitation to Ca(4s3d 1D)

The Ca* + CH3F --> CaF* + CH3 reaction was studied both experimentally and theoretically. The reaction was photoinduced in Ca...CH3F complexes, which were illuminated by a tunable laser in the range 18 000-24 000 cm-1. The absorption band that leads to the reaction extends between 19 000 and 23 000 cm-1. It is formed of three broad overlapping structures corresponding to the excitation of different electronic states of the complex. The two structures of lowest energy were considered in detail. They are associated with two series of respectively 2 and 3 molecular states correlating to Ca(4s3d 1D) + CH3F at infinite separation between Ca and CH3F. The assignment of these structures to specific electronic transitions of the complex stemmed from theoretical calculations where the Ca...CH3F complex is described by a linear Ca-F-C backbone. 2D potential energy surfaces were calculated by associating a pseudopotential description of the [Ca2+] and [F7+] cores, a core polarization operator on calcium, an extensive Gaussian basis, and a treatment of the electronic problem at the CI-MRCI level. All the excited levels correlating to the 4s2 1S, 4s3d 1D, and 4s4p 1P levels of Ca in the Ca + CH3F channel were documented in a calculation that explored the rearrangement channels where either Ca + CH3F or CaF + CH3 are formed. Then, wavepacket calculations on the 2D-PES's allowed one to simulate the absorption spectrum of the complex, in an approximation where the various electronic states of the complex are not coupled together. The assignment above stemmed from this. The second outcome of the calculation was that whatever the excited level of the complex that is considered, the reaction has to proceed through energy barriers. The electronic excitation of the complex on the red side of the absorption band does not seem to deposit enough energy in the system to overcome these barriers (even the lowest one) or to stimulate tunneling reactions. An alternative reaction mechanism involving a transfer to triplet PES's is proposed.

Intramolecular recognition in a jet-cooled short peptide chain: γ-turn helicity probed by a neighbouring residue

gamma-Turn, the shortest secondary structure of peptides, exists as two helical forms gamma(l) and gamma(d) of opposite handedness. The present gas phase study of capped l-Phe-Xxx peptides (Xxx = l-Ala, d-Ala or Aib: aminoisobutyric acid) provides a unique example of intramolecular chiral recognition of the gamma-turn helicity on Ala or Aib by the neighbouring residue Phe within the chain. With the chiral l- or d-Ala residues, the presence of a side-chain operates a discrimination between the two helical forms: one of them is widely favoured over the other (gamma(l) or gamma(d), respectively). This enables us to validate and calibrate the recognition capabilities of the nearby l-Phe residue. The discriminating interactions have been precisely characterized from their spectroscopic UV and IR signatures and identified by comparison with quantum chemistry calculations. Then, in the case of the non-chiral residue Aib, the two helical forms of the gamma-turn, which are simultaneously observed in the jet, have been discriminated and assigned by comparison with the chiral residues. The relative abundances of the diastereomeric forms l-Phe-Aib(gamma(l)) and l-Phe-Aib(gamma(d)) enable us to determine the most efficient recognition configuration.

Transition State Spectroscopy of the Photoinduced Ca + CH3F Reaction. 2. Experimental and Ab Initio Studies of the Free Ca···FCH3 Complex

The Ca* + CH3F --> CaF + CH3 reaction was photoinduced in 1:1 Ca...CH3F complexes formed in a supersonic expansion. The transition state of the reaction was explored by monitoring the electronically excited product, CaF, while scanning the laser that turns on the reaction. Moreover, the electronic structure of the Ca...FCH3 system was studied using ab initio methods by associating a pseudopotential description of the [Ca2+] and [F7+] cores, a core polarization operator on calcium, an extensive Gaussian basis and a treatment of the electronic problem at the CCSD(T) (ground state) and RSPT2 (excited states) level. In this contribution we present experimental results for the free complex and a comparison with the results of a previous experiment where the Ca...CH3F complexes are deposited at the surface of large argon clusters. The ab initio calculations allowed an interpretation of the experimental data in terms of two reaction mechanisms, one involving a partial charge transfer state, the other involving the excitation of the C-F stretch in the CH3F moiety prior to charge transfer.

Transition State Spectroscopy of the Photoinduced Ca + CH3F Reaction. 1. A Cluster Isolated Chemical Reaction Study

The "cluster isolated chemical reactions" technique is used to examine the dynamics of the photoinduced reaction producing electronically excited CaF when 1:1 Ca-CH(3)F complexes are deposited at the surface of large argon clusters. This technique ensures quantitatively that 1:1 complexes are actually at the origin of the observed signals. The reaction is monitored by observing the CaF chemiluminescence while scanning the photoexcitation laser. The resulting action spectrum contains information about the absorption bands of the complex, filtered by the dynamics of the reaction. The observations suggest a profound alteration of the calcium electronic structure and a control of the reaction by the CF stretch in CH(3)F.