Genetic transformation of the frog-killing chytrid fungus Batrachochytrium dendrobatidis

Batrachochytrium dendrobatidis (Bd), a causative agent of chytridiomycosis, is decimating amphibian populations around the world. Bd belongs to the chytrid lineage, a group of early-diverging fungi that are widely used to study fungal evolution. Like all chytrids, Bd develops from a motile form into a sessile, growth form, a transition that involves drastic changes in its cytoskeletal architecture. Efforts to study Bd cell biology, development, and pathogenicity have been limited by the lack of genetic tools with which to test hypotheses about underlying molecular mechanisms. Here, we report the development of a transient genetic transformation system for Bd. We used electroporation to deliver exogenous DNA into Bd cells and detected transgene expression for up to three generations under both heterologous and native promoters. We also adapted the transformation protocol for selection using an antibiotic resistance marker. Finally, we used this system to express fluorescent protein fusions and, as a proof of concept, expressed a genetically encoded probe for the actin cytoskeleton. Using live-cell imaging, we visualized the distribution and dynamics of polymerized actin at each stage of the Bd life cycle, as well as during key developmental transitions. This transformation system enables direct testing of key hypotheses regarding mechanisms of Bd pathogenesis. This technology also paves the way for answering fundamental questions of chytrid cell, developmental, and evolutionary biology

Prevalence and Risk Factors for Early, Undesired Weaning Attributed to Lactation Dysfunction

Background: Breastfeeding durations in the United States fall short of public health objectives. We sought to quantify the prevalence and identify risk factors for early, undesired weaning that mothers attribute to physiologic difficulties with breastfeeding

Variant-specific inflation factors for assessing population stratification at the phenotypic variance level

In modern Whole Genome Sequencing (WGS) epidemiological studies, participant-level data from multiple studies are often pooled and results are obtained from a single analysis. We consider the impact of differential phenotype variances by study, which we term \u27variance stratification\u27. Unaccounted for, variance stratification can lead to both decreased statistical power, and increased false positives rates, depending on how allele frequencies, sample sizes, and phenotypic variances vary across the studies that are pooled. We develop a procedure to compute variant-specific inflation factors, and show how it can be used for diagnosis of genetic association analyses on pooled individual level data from multiple studies. We describe a WGS-appropriate analysis approach, implemented in freely-available software, which allows study-specific variances and thereby improves performance in practice. We illustrate the variance stratification problem, its solutions, and the proposed diagnostic procedure, in simulations and in data from the Trans-Omics for Precision Medicine Whole Genome Sequencing Program (TOPMed), used in association tests for hemoglobin concentrations and BMI

Perinatal eating disorder symptoms, parenting stress, and infant feeding

Introduction: We sought to quantify associations between postpartum eating disorder symptoms, maternal mood, parenting, and infant feeding.Study Design: We analyzed participants in an ongoing longitudinal study of women intending to breastfeed. At baseline all participants underwent a Structured Clinical Interview (SCID) to determine history of anorexia, bulimia, or binge eating disorder (EDHx). At 2 months postpartum, we quantified current symptoms using the Eating Disorder Examination Questionnaire (EDE-Q4). We measured outcomes using the Beck Depression Inventory (BDI), the Edinburgh Postnatal Depression Scale (EPDS), the Parenting Stress Index – Short Form, the Breastfeeding Self-Efficacy Scale, and the modified Differential Emotions Scale (mDES). We used multivariable linear regression to model associations between outcomes and EDE-Q4, EDHx and interactions between EDHx and EDE-Q4.Results: Of the 62 women, 1 had a SCID-verified current eating disorder and 8 had an EDHx. At 2 months postpartum, women with EDHx were more likely to meet clinical thresholds for depression indexed by BDI ≥ 17 (33 vs. 6%, Fisher’s exact p=.03). Adjusting for EDHx, higher EDE-Q4 score was associated with higher BDI (R2 = .45, main effect p<.0001) and reduced breastfeeding self-efficacy (R2=.10, main effect p=.02). EDHx moderated other associations: among women without EDHx, higher EDE-Q4 score was associated with higher EPDS score (R2= .27, interaction p.<05), greater parenting stress (R2=.26, interaction p =.06) and more negative emotions during feeding (R2=.32, interaction p=.03).Conclusion: Maternal eating disorder symptoms were associated with depression symptoms, parenting stress, and negative breastfeeding outcomes. Parenting support may improve outcomes for mothers with disordered eating symptoms

Forecasting the spread of raccoon rabies using a purpose-specific group decisionmaking process

The Centers for Disease Control (CDC) and USDA Wildlife Services (WS) have been involved in an oral rabies vaccination (ORV) program for raccoons (Procyon lotor) that has slowed the westward spread of raccoon rabies. The objective of this study was to forecast the spread of the disease if an ORV zone was not maintained. A group decision-making process was designed to address the forecasting problem and was implemented using a group of 15 experts and 4 support personnel at a meeting at the USDA National Wildlife Research Center. Ten expansion regions were constructed that described the spread of disease at 2-year intervals. This forecast may provide for more accurate cost-benefit analysis of the ORV barrier

ADAM19 and HTR4 Variants and Pulmonary Function: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study

The pulmonary function measures of forced expiratory volume in one second (FEV1) and its ratio to forced vital capacity (FVC) are used in the diagnosis and monitoring of lung diseases and predict cardiovascular mortality in the general population. Genome wide association studies (GWAS) have identified numerous loci associated with FEV1 and FEV1/FVC but the causal variants remain uncertain. We hypothesized that novel or rare variants poorly tagged by GWAS may explain the significant associations between FEV1/FVC and two genes: ADAM19 and HTR4

Identification of Circulating Proteins associated With General Cognitive Function among Middle-Aged and Older adults

Identifying circulating proteins associated with cognitive function may point to biomarkers and molecular process of cognitive impairment. Few studies have investigated the association between circulating proteins and cognitive function. We identify 246 protein measures quantified by the SomaScan assay as associated with cognitive function (p \u3c 4.9E-5, n up to 7289). Of these, 45 were replicated using SomaScan data, and three were replicated using Olink data at Bonferroni-corrected significance. Enrichment analysis linked the proteins associated with general cognitive function to cell signaling pathways and synapse architecture. Mendelian randomization analysis implicated higher levels of NECTIN2, a protein mediating viral entry into neuronal cells, with higher Alzheimer\u27s disease (AD) risk (p = 2.5E-26). Levels of 14 other protein measures were implicated as consequences of AD susceptibility (p \u3c 2.0E-4). Proteins implicated as causes or consequences of AD susceptibility may provide new insight into the potential relationship between immunity and AD susceptibility as well as potential therapeutic targets

Epigenetic and integrative cross-omics analyses of cerebral white matter hyperintensities on MRI

Cerebral white matter hyperintensities on MRI are markers of cerebral small vessel disease, a major risk factor for dementia and stroke. Despite the successful identification of multiple genetic variants associated with this highly heritable condition, its genetic architecture remains incompletely understood. More specifically, the role of DNA methylation has received little attention. We investigated the association between white matter hyperintensity burden and DNA methylation in blood at approximately 450,000 CpG sites in 9,732 middle-aged to older adults from 14 community-based studies. Single-CpG and region-based association analyses were carried out. Functional annotation and integrative cross-omics analyses were performed to identify novel genes underlying the relationship between DNA methylation and white matter hyperintensities. We identified 12 single-CpG and 46 region-based DNA methylation associations with white matter hyperintensity burden. Our top discovery single CpG, cg24202936 (P = 7.6 × 10-8), was associated with F2 expression in blood (P = 6.4 × 10-5), and colocalized with FOLH1 expression in brain (posterior probability =0.75). Our top differentially methylated regions were in PRMT1 and in CCDC144NL-AS1, which were also represented in single-CpG associations (cg17417856 and cg06809326, respectively). Through Mendelian randomization analyses cg06809326 was putatively associated with white matter hyperintensity burden (P = 0.03) and expression of CCDC144NL-AS1 possibly mediated this association. Differentially methylated region analysis, joint epigenetic association analysis, and multi-omics colocalization analysis consistently identified a role of DNA methylation near SH3PXD2A, a locus previously identified in genome-wide association studies of white matter hyperintensities. Gene set enrichment analyses revealed functions of the identified DNA methylation loci in the blood-brain barrier and in the immune response. Integrative cross-omics analysis identified 19 key regulatory genes in two networks related to extracellular matrix organization, and lipid and lipoprotein metabolism. A drug repositioning analysis indicated antihyperlipidemic agents, more specifically peroxisome proliferator-activated receptor alpha, as possible target drugs for white matter hyperintensities. Our epigenome-wide association study and integrative cross-omics analyses implicate novel genes influencing white matter hyperintensity burden, which converged on pathways related to the immune response and to a compromised blood brain barrier possibly due to disrupted cell-cell and cell-extracellular matrix interactions. The results also suggest that antihyperlipidemic therapy may contribute to lowering risk for white matter hyperintensities possibly through protection against blood brain barrier disruption

Integrative pathway genomics of lung function and airflow obstruction

Chronic respiratory disorders are important contributors to the global burden of disease. Genome-wide association studies (GWASs) of lung function measures have identified several trait-associated loci, but explain only a modest portion of the phenotypic variability. We postulated that integrating pathway-based methods with GWASs of pulmonary function and airflow obstruction would identify a broader repertoire of genes and processes influencing these traits. We performed two independent GWASs of lung function and applied gene set enrichment analysis to one of the studies and validated the results using the second GWAS. We identified 131 significantly enriched gene sets associated with lung function and clustered them into larger biological modules involved in diverse processes including development, immunity, cell signaling, proliferation and arachidonic acid. We found that enrichment of gene sets was not driven by GWAS-significant variants or loci, but instead by those with less stringent association P-values. Next, we applied pathway enrichment analysis to a meta-analyzed GWAS of airflow obstruction. We identified several biologic modules that functionally overlapped with those associated with pulmonary function. However, differences were also noted, including enrichment of extracellular matrix (ECM) processes specifically in the airflow obstruction study. Network analysis of the ECM module implicated a candidate gene, matrix metalloproteinase 10 (MMP10), as a putative disease target. We used a knockout mouse model to functionally validate MMP10's role in influencing lung's susceptibility to cigarette smoke-induced emphysema. By integrating pathway analysis with population-based genomics, we unraveled biologic processes underlying pulmonary function traits and identified a candidate gene for obstructive lung diseas

Multiethnic Meta-Analysis Identifies Ancestry-Specific and Cross-Ancestry Loci for Pulmonary Function

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N = 60,552), African (N = 8429), Asian (N = 9959), and Hispanic/Latino (N = 11,775) ethnicities. We identify over 50 additional loci at genome-wide significance in ancestry-specific or multiethnic meta-analyses. Using recent fine-mapping methods incorporating functional annotation, gene expression, and differences in linkage disequilibrium between ethnicities, we further shed light on potential causal variants and genes at known and newly identified loci. Several of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12. Our study highlights the utility of multiethnic and integrative genomics approaches to extend existing knowledge of the genetics of lung function and clinical relevance of implicated loci