Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Predicting Heterogeneity in Patient Response to Morphine Treatment for Neonatal Opioid Withdrawal Syndrome

Infants with neonatal opioid withdrawal syndrome commonly receive morphine treatment to manage their withdrawal signs. However, the effectiveness of this pharmacotherapy in managing the infants' withdrawal signs vary widely. We sought to understand how information available early in infant monitoring can anticipate this treatment response, focusing on early modified Finnegan Neonatal Abstinence Scoring System (FNASS) scores, polygenic risk for opioid dependence (polygenic risk score (PRS)), and drug exposure. Using k-means clustering, we divided the 213 infants in our cohort into 3 groups based on their FNASS scores in the 12 hours before and after the initiation of pharmacotherapy. We found that these groups were pairwise significantly different for risk factors, including methadone exposure, and for in-hospital outcomes, including total morphine received, length of stay, and highest FNASS score. Whereas PRS was not predictive of receipt of treatment, PRS was pairwise significantly different between a subset of the groups. Using tree-based machine learning methods, we then constructed network graphs of the relationships among these groups, FNASS scores, PRS, drug exposures, and in-hospital outcomes. The resulting networks also showed meaningful connection between early FNASS scores and PRS, as well as between both of those and later in-hospital outcomes. These analyses present clinicians with the opportunity to better anticipate infant withdrawal progression and prepare accordingly, whether with expedited morphine treatment or non-pharmacotherapeutic alternative treatments.

Aetiology of diarrhoeal disease and evaluation of viral-bacterial coinfection in children under 5 years old in China: a matched case-control study

Globally, diarrhoeal diseases are the second leading cause of death among children under 5 years old. Few case-control studies on the aetiology of diarrhoea have been conducted in China. A case-control study on 922 children under 5 years old who presented with diarrhoea and individually matched controls was conducted in China between May 2011 and January 2013. Quantitative PCR was used to analyze stool samples for 10 diarrhoeal pathogens. Potential enteric pathogens were detected in 377 (81.8%) of 461 children with diarrhoea and 215 controls (46.6%, p <0.001). Rotavirus, norovirus GII, Shigella and adenovirus were qualitatively associated with diarrhoea. Using receiver operating characteristic curves, the optimal cutoff threshold for defining a symptomatic individual was 72, 5840, and 10(4) copies per reaction for rotavirus (odds ratio 259), norovirus GII (odds ratio 10.6) and Shigella (odds ratio 5.1). The attributable fractions were 0.18 for rotavirus, 0.08 for norovirus GII, 0.01 for Shigella and 0.04 for adenovirus. Coinfections between pathogens were common. Two pairs, rotavirus and adenovirus, and norovirus GII and Salmonella were positively associated. The co-occurrence of rotavirus and sapovirus, astrovirus, enterotoxigenic Escherichia coli or Campylobacter jejuni only occurred in children with disease. Coinfection was not correlated with clinical symptoms. Quantitative data are critical. Our results indicate that increased pathogen loads increase the OR between diarrhoea and rotavirus, norovirus GII and Shigella. Coinfections with rotavirus and norovirus GII are common and occur in a nonrandom distribution. Despite testing for ten diarrhoeal pathogens, over two-thirds of cases do not have a recognized attributable cause.

Discord reigns among nuclear, mitochondrial and phenotypic estimates of divergence in nine lineages of trans-Beringian birds

Proposals for genetic thresholds for species delimitation assume that simple genetic data sets (e.g. mitochondrial sequence data) are correlated with speciation; i.e. such data sets accurately reflect organismal lineage divergence. We used taxonomically stratified phenotypic levels of differentiation (populations, subspecies and species) among nine avian lineages using paired, trans-Beringian samples from three lineages each in three orders (Anseriformes, Charadriiformes, and Passeriformes) to test this assumption. Using mitochondrial DNA sequence data and nuclear genomic data (amplified fragment length polymorphisms), we found a lack of concordance between these two genomes in their respective estimates of divergence and little or no relationship between phenotype (taxonomic relatedness) and genetic differentiation between taxon pairs. There are several possible reasons for the discord observed (e.g. selection on one of the genomes or perhaps lineage sorting), but the implications are that genetic estimates of lineage divergence may not be correlated with estimates from other parts of the genome, are not well correlated with the speciation process and are thus not reliable indicators of species limits.

Working through polytomies: auklets revisited

Polytomies, or phylogenetic "bushes", are the result of a series of internodes occurring in a short period of evolutionary time (which can result in data that do not contain enough information), or data that have too much homoplasy to resolve a bifurcating branching pattern. In this study we used the Aethia auklet polytomy to explore the effectiveness of different methods for resolving polytomies: mitochondrial DNA gene choice, number of individuals per species sampled, model of molecular evolution, and AFLP loci. We recovered a fully-resolved phylogeny using NADH dehydrogenase subunit 2 (ND2) sequence data under two different Bayesian models. We were able to corroborate this tree under one model with an expanded mtDNA dataset. Effectiveness of additional intraspecific sampling varied with node, and fully 20% of the subsampled datasets failed to return a congruent phylogeny when we sampled only one or two individuals per species. We did not recover a resolved phylogeny using AFLP data. Conflict in the AFLP dataset showed that nearly all possible relationships were supported at low levels of confidence, suggesting that either AFLPs are not useful at the genetic depth of the Aethia auklet radiation (7-9% divergent in the mtDNA ND2 gene), perhaps resulting in too much homoplasy, or that the Aethia auklets have experienced incomplete lineage sorting at many nuclear loci.

Multilocus phylogeography of a holarctic duck: colonization of north america from eurasia by gadwall (Anas strepera)

More than 100 species of birds have Holarctic distributions extending across Eurasia and North America, and many of them likely achieved these distributions by recently colonizing one continent from the other. Mitochondrial DNA (mtDNA) and five nuclear introns were sequenced to test the direction and timing of colonization for a Holarctic duck, the gadwall (Anas strepera). Three lines of evidence suggest gadwalls colonized North America from Eurasia. First, New World (NW) gadwalls had fewer alleles at every locus and 61% of the allelic richness found in Old World (OW) gadwalls. Second, NW gadwalls had lower mtDNA allelic richness than other NW ducks. Third, coalescent analysis suggested that less than 5% of the ancestral population contributed to NW gadwalls at the time of divergence. Gadwalls likely colonized North America during the Late Pleistocene (approximately 81,000 years ago), but the confidence interval on that estimate was large (8500-450,000 years ago). Intercontinental gene flow and selection also likely contributed to genetic diversity in gadwalls. This study illustrates the use of multiple loci and coalescent analyses for critically testing a priori hypotheses regarding dispersal and colonization and provides an independent datapoint supporting an OW to NW bias in the direction of colonization.