Search

Insects have sensitive olfactory systems to interact with environment and respond to the change in host plant conditions. Key genes in the system can be potential targets for developing new and efficient pest behaviour control methods. Riptortus pedestris is an important soybean pest in East Asia and has caused serious damage to the soybean plants in Huang-Huai-Hai region of China. However, the current treatment of pests is dominated by chemical insecticides and lacks efficient sustainable prevention and control technologies. In this study, we identified 49 putative odorant-binding proteins (OBPs) (43 were new genes) and 25 chemosensory proteins (CSPs) (17 were new genes) in R. pedestris genome. These OBP and CSP genes are clustered in highly conserved groups from other hemipteran species in phylogenetic trees. Most RpedOBPs displayed antennal-biased expression. Among the 49 RpedOBPs, 33 were significantly highly expressed in the antennae, including three male-biased and nine female-biased. While many RpedCSPs were detected both in the antennae and in non-antennal tissues, only 11 RpedCSPs displayed antennal-biased expression, in which four RpedCSPs were male-biased and five RpedCSPs were female-biased. Some OBP and CSP genes showed sex-biased expression profiles. Our results not only provide a foundation for future exploration of the functions of RpedOBPs and RpedCSPs but also aid in developing environmentally friendly insecticides in the future.

Reversal reactions (RRs) in leprosy are characterized by a reduction in the number of bacilli in lesions associated with an increase in cell-mediated immunity against the intracellular bacterium Mycobacterium leprae, the causative pathogen of leprosy. To identify the mechanisms that contribute to cell-mediated immunity in leprosy, we measured changes in the whole blood-derived transcriptome of patients with leprosy before, during and after RR. We identified an 'RR signature' of 1017 genes that were upregulated at the time of the clinical diagnosis of RR. Using weighted gene correlated network analysis (WGCNA), we detected a module of 794 genes, bisque4, that was significantly correlated with RR, of which 434 genes were part of the RR signature. An enrichment for both IFN-γ and IFN-β downstream gene pathways was present in the RR signature as well as the RR upregulated genes in the bisque4 module, including those encoding proteins of the guanylate binding protein (GBP) family that contributes to antimicrobial responses against mycobacteria. Specifically, GBP1, GBP2, GBP3 and GBP5 mRNAs were upregulated in the RR peripheral blood transcriptome, with GBP1, GBP2 and GBP5 mRNAs also upregulated in the RR disease lesion transcriptome. These data indicate that RRs involve a systemic upregulation of IFN-γ downstream genes including GBP family members as part of the host antimicrobial response against mycobacteria. [ABSTRACT FROM AUTHOR]

Peripheral Blood gene expression is widely used in the discovery of biomarkers and development of therapeutics. Recently, a spate of commercial blood collection and preservation systems have been introduced with proprietary variations that may differentially impact the transcriptomic profiles. Comparative analysis of these collection platforms will help optimize protocols to detect, identify, and reproducibly validate true biological variance among subjects. In the current study, we tested two recently introduced whole blood collection methods, RNAgard® and PAXgene® RNA, in addition to the traditional method of peripheral blood mononuclear cells (PBMCs) separated from whole blood and preserved in Trizol reagent. Study results revealed striking differences in the transcriptomic profiles from the three different methods that imply ex vivo changes in gene expression occurred during the blood collection, preservation, and mRNA extraction processes. When comparing the ability of the three preservation methods to accurately capture individuals' expression differences, RNAgard® outperformed PAXgene® RNA, and both showed better individual separation of transcriptomic profiles than PBMCs. Hence, our study recommends using a single blood collection platform, and strongly cautions against combining methods during the course of a defined study. [ABSTRACT FROM AUTHOR]

Introduction: The prevalence of bacteremia caused by Gram negative non-fermentative (GNNF) bacteria has been increasing globally over the past decade. Many studies have investigated their epidemiology but focus on the common GNNF including Pseudomonas aeruginosa and Acinetobacter baumannii. Knowledge of the uncommon GNNF bacteremias is very limited. This study explores invasive bloodstream infection GNNF isolates that were initially unidentified after testing with standard microbiological techniques. All isolations were made during laboratory-based surveillance activities in two rural provinces of Thailand between 2006 and 2014. Methods: A subset of GNNF clinical isolates (204/947), not identified by standard manual biochemical methodologies were run on the BD Phoenix automated identification and susceptibility testing system. If an organism was not identified (12/204) DNA was extracted for whole genome sequencing (WGS) on a MiSeq platform and data analysis performed using 3 web-based platforms: Taxonomer, CGE KmerFinder and One Codex. Results: The BD Phoenix automated identification system recognized 92% (187/204) of the GNNF isolates, and because of their taxonomic complexity and high phenotypic similarity 37% (69/187) were only identified to the genus level. Five isolates grew too slowly for identification. Antimicrobial sensitivity (AST) data was not obtained for 93/187 (50%) identified isolates either because of their slow growth or their taxa were not in the AST database associated with the instrument. WGS identified the 12 remaining unknowns, four to genus level only. Conclusion: The GNNF bacteria are of increasing concern in the clinical setting, and our inability to identify these organisms and determine their AST profiles will impede treatment. Databases for automated identification systems and sequencing annotation need to be improved so that opportunistic organisms are better covered. [ABSTRACT FROM AUTHOR]

Background: The full scope of the genes expressed by schistosomes during intramolluscan development has yet to be characterized. Understanding the gene products deployed by larval schistosomes in their snail hosts will provide insights into their establishment, maintenance, asexual reproduction, ability to castrate their hosts, and their prolific production of human-infective cercariae. Using the Illumina platform, the intramolluscan transcriptome of Schistosoma mansoni was investigated in field-derived specimens of the prominent vector species Biomphalaria pfeifferi at 1 and 3 days post infection (d) and from snails shedding cercariae. These S. mansoni samples were derived from the same snails used in our complementary B. pfeifferi transcriptomic study. We supplemented this view with microarray analyses of S. mansoni from B. glabrata at 2d, 4d, 8d, 16d, and 32d to highlight robust features of S. mansoni transcription, even when a different technique and vector species was used. Principal findings: Transcripts representing at least 7,740 (66%) of known S. mansoni genes were expressed during intramolluscan development, with the greatest number expressed in snails shedding cercariae. Many transcripts were constitutively expressed throughout development featuring membrane transporters, and metabolic enzymes involved in protein and nucleic acid synthesis and cell division. Several proteases and protease inhibitors were expressed at all stages, including some proteases usually associated with cercariae. Transcripts associated with G-protein coupled receptors, germ cell perpetuation, and stress responses and defense were well represented. We noted transcripts homologous to planarian anti-bacterial factors, several neural development or neuropeptide transcripts including neuropeptide Y, and receptors that may be associated with schistosome germinal cell maintenance that could also impact host reproduction. In at least one snail the presence of larvae of another digenean species (an amphistome) was associated with repressed S. mansoni transcriptional activity. Conclusions/Significance: This in vivo study, emphasizing field-derived snails and schistosomes, but supplemented with observations from a lab model, provides a distinct view from previous studies of development of cultured intramolluscan stages from lab-maintained organisms. We found many highly represented transcripts with suspected or unknown functions, with connection to intramolluscan development yet to be elucidated. [ABSTRACT FROM AUTHOR]

Background: Excess weight gain throughout adulthood can lead to adverse clinical outcomes and are influenced by complex factors that are difficult to measure in free-living individuals. Metabolite profiling offers an opportunity to systematically discover new predictors for weight gain that are relatively easy to measure compared to traditional approaches. Methods and results: Using baseline metabolite profiling data of middle-aged individuals from the Framingham Heart Study (FHS; n = 1,508), we identified 42 metabolites associated (p < 0.05) with longitudinal change in body mass index (BMI). We performed stepwise linear regression to select 8 of these metabolites to build a metabolite risk score (MRS) for predicting future weight gain. We replicated the MRS using data from the Mexico City Diabetes Study (MCDS; n = 768), in which one standard deviation increase in the MRS corresponded to ~0.03 increase in BMI (kg/m2) per year (i.e. ~0.09 kg/year for a 1.7 m adult). We observed that none of the available anthropometric, lifestyle, and glycemic variables fully account for the MRS prediction of weight gain. Surprisingly, we found the MRS to be strongly correlated with baseline insulin sensitivity in both cohorts and to be negatively predictive of T2D in MCDS. Genome-wide association study of the MRS identified 2 genome-wide (p < 5 × 10−8) and 5 suggestively (p < 1 × 10−6) significant loci, several of which have been previously linked to obesity-related phenotypes. Conclusions: We have constructed and validated a generalizable MRS for future weight gain that is an independent predictor distinct from several other known risk factors. The MRS captures a composite biological picture of weight gain, perhaps hinting at the anabolic effects of preserved insulin sensitivity. Future investigation is required to assess the relationships between MRS-predicted weight gain and other obesity-related diseases. [ABSTRACT FROM AUTHOR]

Genomic resources for the domestic dog have improved with the widespread adoption of a 173k SNP array platform and updated reference genome. SNP arrays of this density are sufficient for detecting genetic associations within breeds but are underpowered for finding associations across multiple breeds or in mixed-breed dogs, where linkage disequilibrium rapidly decays between markers, even though such studies would hold particular promise for mapping complex diseases and traits. Here we introduce an imputation reference panel, consisting of 365 diverse, whole-genome sequenced dogs and wolves, which increases the number of markers that can be queried in genome-wide association studies approximately 130-fold. Using previously genotyped dogs, we show the utility of this reference panel in identifying potentially novel associations, including a locus on CFA20 significantly associated with cranial cruciate ligament disease, and fine-mapping for canine body size and blood phenotypes, even when causal loci are not in strong linkage disequilibrium with any single array marker. This reference panel resource will improve future genome-wide association studies for canine complex diseases and other phenotypes. [ABSTRACT FROM AUTHOR]

Objective: Moyamoya disease (MMD) is a chronic occlusive cerebrovascular disease with unknown etiology, sharing many similar clinical symptoms with other vascular disorders. This study aimed to investigate gene dysregulation in peripheral blood of MMD and compare it with other vascular disorders. Methods: Transcriptomic profiles of 12 MMD patients and 8 healthy controls were obtained using RNA sequencing. Differentially expressed genes (DEGs) were identified and several were validated by quantitative real-time PCR in independent samples. Biological pathway enrichment analysis of DEGs and deconvolution of leukocyte subsets in peripheral blood were performed. Expression profiles for other vascular diseases were downloaded from public database and consistent DEGs were calculated. Gene set enrichment analysis (GSEA) was conducted to compare gene dysregulation pattern between MMD and other vascular diseases. Results: A total of 533 DEGs were identified for MMD. Up-regulated genes were mainly involved in extracellular matrix (ECM) organization, whereas down-regulated genes were primarily associated with inflammatory and immune responses. As for cell populations, significantly increased naïve B cells and naïve CD4 cells as well as obviously decreased resting natural killer cells were observed in peripheral blood of MMD patients. GSEA analysis indicated that only up-regulated genes of ischemic stroke and down-regulated genes of coronary artery disease and myocardial infarction were enriched in up-regulated and down-regulated genes of MMD, respectively. Conclusion: Dysregulated genes in peripheral blood of MMD mainly played key roles in ECM organization, inflammatory and immune responses. This gene dysregulation pattern was specific compared with other vascular diseases. Besides, naïve B cells, naïve CD4 cells and resting natural killer cells were aberrantly disrupted in peripheral blood of MMD patients. These results will help elucidate the complicated pathogenic mechanism of MMD. [ABSTRACT FROM AUTHOR]

In haplodiplontic lineages, sexual reproduction occurs in haploid parents without meiosis. Although widespread in multicellular lineages such as brown algae (Phaeophyceae), haplodiplontic gametogenesis has been little studied at the molecular level. We addressed this by generating an annotated reference transcriptome for the gametophytic phase of the sugar kelp, Saccharina latissima. Transcriptional profiles of microscopic male and female gametophytes were analysed at four time points during the transition from vegetative growth to gametogenesis. Gametogenic signals resulting from a switch in culture irradiance from red to white light activated a core set of genes in a sex-independent manner, involving rapid activation of ribosome biogenesis, transcription and translation related pathways, with several acting at the post-transcriptional or post-translational level. Additional genes regulating nutrient acquisition and key carbohydrate-energy pathways were also identified. Candidate sex-biased genes under gametogenic conditions had potentially key roles in controlling female- and male-specific gametogenesis. Among these were several sex-biased or -specific E3 ubiquitin-protein ligases that may have important regulatory roles. Females specifically expressed several genes that coordinate gene expression and/or protein degradation, and the synthesis of inositol-containing compounds. Other female-biased genes supported parallels with oogenesis in divergent multicellular lineages, in particular reactive oxygen signalling via an NADPH-oxidase. Males specifically expressed the hypothesised brown algal sex-determining factor. Male-biased expression mainly involved upregulation of genes that control mitotic cell proliferation and spermatogenesis in other systems, as well as multiple flagella-related genes. Our data and results enhance genome-level understanding of gametogenesis in this ecologically and economically important multicellular lineage. [ABSTRACT FROM AUTHOR]

Background: Tyrosine kinase inhibitors targeted to anaplastic lymphoma kinase (ALK) have been demonstrated to be effective for lung cancer patients with an ALK fusion gene. Application of liquid biopsy, i.e., detection and quantitation of the fusion product in plasma cell-free DNA (cfDNA), could improve clinical practice. To detect ALK fusions, because fusion breakpoints occur somewhere in intron 19 of the ALK gene, sequencing of the entire intron is required to locate breakpoints. Results: We constructed a target sequencing system using an adapter and a set of primers that cover the entire ALK intron 19. This system can amplify fragments, including breakpoints, regardless of fusion partners. The data analysis pipeline firstly detected fusions by alignment to selected target sequences, and then quantitated the fusion alleles aligning to the identified breakpoint sequences. Performance was validated using 20 cfDNA samples from ALK-positive non-small cell lung cancer patients and samples from 10 healthy volunteers. Sensitivity and specificity were 50 and 100%, respectively. Conclusions: We demonstrated that PCR-based target sequencing using a tiling primer set and two-step mapping/alignment quantitatively detected ALK fusions in cfDNA from lung cancer patients. The system offers an alternative to existing approaches based on hybridization capture. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWAS) have been successful in identifying genetic variants associated with complex diseases. However, association analyses between genotypes and phenotypes are not straightforward due to the complex relationships between genetic and environmental factors. Moreover, multiple correlated phenotypes further complicate such analyses. To resolve this complexity, we present an analysis using structural equation modeling (SEM). Unlike current methods that focus only on identifying direct associations between diseases and genetic variants such as single-nucleotide polymorphisms (SNPs), our method introduces the effects of intermediate phenotypes, which are related phenotypes distinct from the target, into the systematic genetic study of diseases. Moreover, we consider multiple diseases simultaneously in a single model. The procedure can be summarized in four steps: 1) selection of informative SNPs, 2) extraction of latent variables from the selected SNPs, 3) investigation of the relationships among intermediate phenotypes and diseases, and 4) construction of an SEM. As a result, a quantitative map can be drawn that simultaneously shows the relationship among multiple SNPs, phenotypes, and diseases. In this study, we considered two correlated diseases, hypertension and type 2 diabetes (T2D), which are known to have a substantial overlap in their disease mechanism and have significant public health implications. As intermediate phenotypes for these diseases, we considered three obesity-related phenotypes—subscapular skin fold thickness, body mass index, and waist circumference—as traits representing subcutaneous adiposity, overall adiposity, and abdominal adiposity, respectively. Using GWAS data collected from the Korea Association Resource (KARE) project, we applied the proposed SEM process. Among 327,872 SNPs, 24 informative SNPs were selected in the first step (p<1.0E-05). Ten latent variables were generated in step 2. After an exploratory analysis, we established a path diagram among phenotypes and diseases in step 3. Finally, in step 4, we produced a quantitative map with paths moving from specific SNPs to hypertension through intermediate phenotypes and T2D. The resulting model had high goodness-of-fit measures (χ2 = 536.52, NFI = 0.997, CFI = 0.998, GFI = 0.995, AGFI = 0.993, RMSEA = 0.012). [ABSTRACT FROM AUTHOR]

Background: Lower birth weight is associated with diabetes although the underlying mechanisms are unclear. Muscle mass could be a modifiable link and hence a target of intervention. We assessed the associations of birth weight with muscle and fat mass observationally in a population with little socio-economic patterning of birth weight and using Mendelian randomization (MR) for validation. Methods: In the population-representative “Children of 1997” birth cohort (n = 8,327), we used multivariable linear regression to assess the adjusted associations of birth weight (kg) with muscle mass (kg) and body fat (%) at ~17.5 years. Genetically predicted birth weight (effect size) was applied to summary genetic associations with fat-free mass and fat mass (kg) from the UK Biobank (n = ~331,000) to obtain unconfounded estimates using inverse-variance weighting. Results: Observationally, birth weight was positively associated with muscle mass (3.29 kg per kg birth weight, 95% confidence interval (CI) 2.83 to 3.75) and body fat (1.09% per kg birth weight, 95% CI 0.54 to 1.65). Stronger associations with muscle mass were observed in boys than in girls (p for interaction 0.004). Using MR, birth weight was positively associated with fat-free mass (0.77 kg per birth weight z-score, 95% CI 0.22 to 1.33) and fat mass (0.58, 95% CI 0.01 to 1.15). No difference by sex was evident. Conclusion: Higher birth weight increasing muscle mass may be relevant to lower birth weight increasing the risk of diabetes and suggests post-natal muscle mass as a potential target of intervention. [ABSTRACT FROM AUTHOR]

Parkinson’s disease (PD) is an irreversible and incurable multigenic neurodegenerative disorder. It involves progressive loss of mid brain dopaminergic neurons in the substantia nigra pars compacta (SN). We compared brain gene expression profiles with those from the peripheral blood cells of a separate sample of PD patients to identify disease-associated genes. Here, we demonstrate the use of gene expression profiling of brain and blood for detecting valid targets and identifying early PD biomarkers. Implementing this systematic approach, we discovered putative PD risk genes in brain, delineated biological processes and molecular functions that may be particularly disrupted in PD and also identified several putative PD biomarkers in blood. 20 of the differentially expressed genes in SN were also found to be differentially expressed in the blood. Further application of this methodology to other brain regions and neurological disorders should facilitate the discovery of highly reliable and reproducible candidate risk genes and biomarkers for PD. The identification of valid peripheral biomarkers for PD may ultimately facilitate early identification, intervention, and prevention efforts as well. [ABSTRACT FROM AUTHOR]

The confluence of deep sequencing and powerful machine learning is providing an unprecedented peek at the darkest of the dark genomic matter, the non-coding genomic regions lacking any functional annotation. While deep sequencing uncovers rare tumor variants, the heterogeneity of the disease confounds the best of machine learning (ML) algorithms. Here we set out to answer if the dark-matter of the genome encompass signals that can distinguish the fine subtypes of disease that are otherwise gnomically indistinguishable. We introduce a novel stochastic regularization, ReVeaL, that empowers ML to discriminate subtle cancer subtypes even from the same ‘cell of origin’. Analogous to heritability, implicitly defined on whole genome, we use predictability (F1 score) definable on portions of the genome. In an effort to distinguish cancer subtypes using dark-matter DNA, we applied ReVeaL to a new WGS dataset from 727 patient samples with seven forms of hematological cancers and assessed the predictivity over several genomic regions including genic, non-dark, non-coding, non-genic, and dark. ReVeaL enabled improved discrimination of cancer subtypes for all segments of the genome. The non-genic, non-coding and dark-matter had the highest F1 scores, with dark-matter having the highest level of predictability. Based on ReVeaL’s predictability of different genomic regions, dark-matter contains enough signal to significantly discriminate fine subtypes of disease. Hence, the agglomeration of rare variants, even in the hitherto unannotated and ill-understood regions of the genome, may play a substantial role in the disease etiology and deserve much more attention. [ABSTRACT FROM AUTHOR]

The haptophyte alga Emiliania huxleyi is the most abundant coccolithophore in the modern ocean and produces elaborate calcite crystals, called coccolith, in a separate intracellular compartment known as the coccolith vesicle. Despite the importance of biomineralization in coccolithophores, the molecular mechanism underlying it remains unclear. Understanding this precise machinery at the molecular level will provide the knowledge needed to enable further manipulation of biomineralization. In our previous study, altering the calcium concentration modified the calcifying ability of E. huxleyi CCMP371. Therefore in this study, we tested E. huxleyi cells acclimated to three different calcium concentrations (0, 0.1, and 10 mM). To understand the whole transcript profile at different calcium concentrations, RNA-sequencing was performed and used for de novo assembly and annotation. The differentially expressed genes (DEGs) among the three different calcium concentrations were analyzed. The functional classification by gene ontology (GO) revealed that ‘intrinsic component of membrane’ was the most enriched of the GO terms at the ambient calcium concentration (10 mM) compared with the limited calcium concentrations (0 and 0.1 mM). Moreover, the DEGs in those comparisons were enriched mainly in ‘secondary metabolites biosynthesis, transport and catabolism’ and ‘signal transduction mechanisms’ in the KOG clusters and ‘processing in endoplasmic reticulum’, and ‘ABC transporters’ in the KEGG pathways. Furthermore, metabolic pathways involved in protein synthesis were enriched among the differentially expressed proteins. The results of this study provide a molecular profile for understanding the expression of transcripts and proteins in E. huxleyi at different calcium concentrations, which will help to identify the detailed mechanism of its calcification. [ABSTRACT FROM AUTHOR]

Streptococcus agalactiae(GBS) is a seriously threat to the farmed tilapia, and oral vaccination was considered to be the most desirable means which requires deep understanding of virulence mechanism of the fish-derived GBS. Our previous genome study of the fish-derived attenuated strain YM001 showed that there were two large deletions in YM001 compared to its parental virulent strain HN016. In this study, a combined transcriptomic and proteomic analysis was performed on YM001 and HN016 strains, and the important genes were verified by RT-qPCR in bacteria strains and infected-tilapia tissues. Overall, we have shown that a total of 958 genes and 331 proteins were significantly differential expressed between YM001 and HN016. By functional annotation of these DEGs and DEPs, genes that were enriched in pentose phosphate pathway(pgm, ptsG, pgi pfkA, fbaA and FBP3) and pyruvate metabolism pathway(pdhA, pdhB, pdhC and pdhD) were identifed as important candidate genes for leads low growth ability in attenuated strain, which may be an important reasons leading virulence attenuation in the end. The expression levels the candidate genes in pentose phosphate pathway and pyruvate metabolism pathway were significant differential expressed in tilapia’ brain and spleen when infected with YM001 and HN016. Our study indicated that the pentose phosphate pathway and pyruvate metabolism pathway that affecting the growth of the strain may be one of the important reasons for the virulence attenuation in HN016. [ABSTRACT FROM AUTHOR]

Rationale: In virtually all models of heart failure, prognosis is determined by right ventricular (RV) function; thus, understanding the cellular mechanisms contributing to RV dysfunction is critical. Whole organ remodeling is associated with cell-specific changes, including cardiomyocyte dedifferentiation and activation of cardiac fibroblasts (Cfib) which in turn is linked to disorganization of cytoskeletal proteins and loss of sarcomeric structures. However, how these cellular changes contribute to RV function remains unknown. We’ve previously shown significant organ-level RV dysfunction in a large animal model of pulmonary hypertension (PH) which was not mirrored by reduced function of isolated cardiomyocytes. We hypothesized that factors produced by the endogenous Cfib contribute to global RV dysfunction by generating a heterogeneous cellular environment populated by dedifferentiated cells. Objective: To determine the effect of Cfib conditioned media (CM) from the PH calf (PH-CM) on adult rat ventricular myocytes (ARVM) in culture. Methods and results: Brief exposure (<2 days) to PH-CM results in rapid, marked dedifferentiation of ARVM to a neonatal-like phenotype exhibiting spontaneous contractile behavior. Dedifferentiated cells maintain viability for over 30 days with continued expression of cardiomyocyte proteins including TnI and α-actinin yet exhibit myofibroblast characteristics including expression of α-smooth muscle actin. Using a bioinformatics approach to identify factor(s) that contribute to dedifferentiation, we found activation of the PH Cfib results in a unique transcriptome correlating with factors both in the secretome and with activated pathways in the dedifferentiated myocyte. Further, we identified upregulation of periostin in the Cfib and CM, and demonstrate that periostin is sufficient to drive cardiomyocyte dedifferentiation. Conclusions: These data suggest that paracrine factor(s) released by Cfib from the PH calf signal a phenotypic transformation in a population of cardiomyocytes that likely contributes to RV dysfunction. Therapies targeting this process, such as inhibition of periostin, have the potential to prevent RV dysfunction. [ABSTRACT FROM AUTHOR]

CD36 is a platelet membrane glycoprotein whose engagement with oxidized low-density lipoprotein (oxLDL) results in platelet activation. The CD36 gene has been associated with platelet count, platelet volume, as well as lipid levels and CVD risk by genome-wide association studies. Platelet CD36 expression levels have been shown to be associated with both the platelet oxLDL response and an elevated risk of thrombo-embolism. Several genomic variants have been identified as associated with platelet CD36 levels, however none have been conclusively demonstrated to be causative. We screened 81 expression quantitative trait loci (eQTL) single nucleotide polymorphisms (SNPs) associated with platelet CD36 expression by a Massively Parallel Reporter Assay (MPRA) and analyzed the results with a novel Bayesian statistical method. Ten eQTLs located 13kb to 55kb upstream of the CD36 transcriptional start site of transcript ENST00000309881 and 49kb to 92kb upstream of transcript ENST00000447544, demonstrated significant transcription shifts between their minor and major allele in the MPRA assay. Of these, rs2366739 and rs1194196, separated by only 20bp, were confirmed by luciferase assay to alter transcriptional regulation. In addition, electromobility shift assays demonstrated differential DNA:protein complex formation between the two alleles of this locus. Furthermore, deletion of the genomic locus by CRISPR/Cas9 in K562 and Meg-01 cells results in upregulation of CD36 transcription. These data indicate that we have identified a variant that regulates expression of CD36, which in turn affects platelet function. To assess the clinical relevance of our findings we used the PhenoScanner tool, which aggregates large scale GWAS findings; the results reinforce the clinical relevance of our variants and the utility of the MPRA assay. The study demonstrates a generalizable paradigm for functional testing of genetic variants to inform mechanistic studies, support patient management and develop precision therapies. [ABSTRACT FROM AUTHOR]

Background: Although the tprK gene of Treponema pallidum are thought to play a critical role in the pathogenesis of syphilis, the profile of variations in tprK during the development of human syphilis infection have remained unclear. Methods/Principal findings: Through next-generation sequencing, we compared the tprK gene of 14 secondary syphilis patients with that of 14 primary syphilis patients, and the results showed an increased number of variants within the seven V regions of the tprK gene in the secondary syphilis samples. The length of the sequences within each V region also presented a 3-bp changing pattern. Interestingly, the frequencies of predominant sequences within the V regions in the secondary syphilis samples were generally decreased compared with those found in the primary syphilis samples, particularly in the V7 region, where a frequency below 60% was found in up to 57% (8/14) of all secondary samples compared with 7% (1/14) of all primary samples. Moreover, the number of minor variants distributed between frequencies of 10 and 49.9% was increased. The alignment of all amino acid sequences within each V region of the primary and secondary syphilis samples revealed that some amino acid sequences, particularly the amino acid sequences IASDGGAIKH and IASEDGSAGNLKH in V1, were highly stable. Additionally, the amino acid sequences in V6 also exhibited notable intrastrain heterogeneity and were likely to form a strain-specific pattern at the interstrain level. Conclusions: The identification of different profiles of the tprK gene in primary and secondary syphilis patients indicated that the tprK gene of T. pallidum undergoes constant variation to result in the best adaptation to the host. The highly stable peptides found in V1 are likely promising potential vaccine components. The highly heterogenetic regions (e.g., V6) could help to understand the role of tprK in immune evasion. [ABSTRACT FROM AUTHOR]

Major depression is a common and severe psychiatric disorder with a highly polygenic genetic architecture. Genome-wide association studies have successfully identified multiple independent genetic loci that harbour variants associated with major depression, but the exact causal genes and biological mechanisms are largely unknown. Tissue-specific network approaches may identify molecular mechanisms underlying major depression and provide a biological substrate for integrative analyses. We provide a framework for the identification of individual risk genes and gene co-expression networks using genome-wide association summary statistics and gene expression information across multiple human brain tissues and whole blood. We developed a novel gene-based method called eMAGMA that leverages tissue-specific eQTL information to identify 99 biologically plausible risk genes associated with major depression, of which 58 are novel. Among these novel associations is Complement Factor 4A (C4A), recently implicated in schizophrenia through its role in synaptic pruning during postnatal development. Major depression risk genes were enriched in gene co-expression modules in multiple brain tissues and the implicated gene modules contained genes involved in synaptic signalling, neuronal development, and cell transport pathways. Modules enriched with major depression signals were strongly preserved across brain tissues, but were weakly preserved in whole blood, highlighting the importance of using disease-relevant tissues in genetic studies of psychiatric traits. We identified tissue-specific genes and gene co-expression networks associated with major depression. Our novel analytical framework can be used to gain fundamental insights into the functioning of the nervous system in major depression and other brain-related traits. [ABSTRACT FROM AUTHOR]

Neutrophil extracellular traps (NET) formation is part of the neutrophil response to infections, but excessive or inappropriate NETosis may trigger the production of autoantibodies and cause organ damage in autoimmune disorders. Spontaneously netting neutrophils are not frequent and induction of NET in vitro by selected stimuli is necessary to investigate their structure. In the present work, the protein composition and post-translational modifications of NET produced under different stimuli have been studied by means of proteomic analysis. Neutrophils from healthy donors were stimulated by PMA, A23187, Escherichia coli LPS or untreated; after three hours, cells were washed, treated with DNase and supernatants collected for mass spectrometry. Data were analyzed by unsupervised hierarchical clustering analyses. We identified proteins contained in NETs of any source or exclusive of one stimulus: LPS-induced and spontaneous NET diverge in protein composition, while PMA- and A23187-induced NET appear more similar. Among the post-translational modifications we examined, methionine sulfoxidation is frequent especially in PMA- and LPS-induced NETs. Myeloperoxidase is the protein more extensively modified. Thus, proteomic analysis indicates that NETs induced by different stimuli are heterogeneous in terms of both protein composition and post-translational modifications, suggesting that NET induced in different conditions may have different biological effects. [ABSTRACT FROM AUTHOR]

While many disease-associated single nucleotide polymorphisms (SNPs) are associated with gene expression (expression quantitative trait loci, eQTLs), a large proportion of complex disease genome-wide association study (GWAS) variants are of unknown function. Some of these SNPs may contribute to disease by regulating gene splicing. Here, we investigate whether SNPs that are associated with alternative splicing (splice QTL or sQTL) can identify novel functions for existing GWAS variants or suggest new associated variants in chronic obstructive pulmonary disease (COPD). RNA sequencing was performed on whole blood from 376 subjects from the COPDGene Study. Using linear models, we identified 561,060 unique sQTL SNPs associated with 30,333 splice sites corresponding to 6,419 unique genes. Similarly, 708,928 unique eQTL SNPs involving 15,913 genes were detected at 10% FDR. While there is overlap between sQTLs and eQTLs, 55.3% of sQTLs are not eQTLs. Co-localization analysis revealed that 7 out of 21 loci associated with COPD (p<1x10-6) in a published GWAS have at least one shared causal variant between the GWAS and sQTL studies. Among the genes identified to have splice sites associated with top GWAS SNPs was FBXO38, in which a novel exon was discovered to be protective against COPD. Importantly, the sQTL in this locus was validated by qPCR in both blood and lung tissue, demonstrating that splice variants relevant to lung tissue can be identified in blood. Other identified genes included CDK11A and SULT1A2. Overall, these data indicate that analysis of alternative splicing can provide novel insights into disease mechanisms. In particular, we demonstrated that SNPs in a known COPD GWAS locus on chromosome 5q32 influence alternative splicing in the gene FBXO38. [ABSTRACT FROM AUTHOR]

The evolutionarily conserved nature of the few well-known anti-aging interventions that affect lifespan, such as caloric restriction, suggests that aging-related research in model organisms is directly relevant to human aging. Since human lifespan is a complex trait, a systems-level approach will contribute to a more comprehensive understanding of the underlying aging landscape. Here, we integrate evolutionary and functional information of normal aging across human and model organisms at three levels: gene-level, process-level, and network-level. We identify evolutionarily conserved modules of normal aging across diverse taxa, and notably show proteostasis to be conserved in normal aging. Additionally, we find that mechanisms related to protein quality control network are enriched for genes harboring genetic variants associated with 22 age-related human traits and associated to caloric restriction. These results demonstrate that a systems-level approach, combined with evolutionary conservation, allows the detection of candidate aging genes and pathways relevant to human normal aging. [ABSTRACT FROM AUTHOR]

Purpose: The role of insulin resistance (IR) in developing postmenopausal breast cancer has not been thoroughly resolved and may be confounded by lifestyle factors such as obesity. We examined whether genetically determined IR is causally associated with breast cancer risk. Methods: We conducted Mendelian randomization (MR) analyses using individual-level data from our previous meta-analysis of a genome-wide association study (GWAS) (n = 11,109 non-Hispanic white postmenopausal women). Four single-nucleotide polymorphisms were associated with fasting glucose (FG), 2 with fasting insulin (FI), and 6 with homeostatic model assessment–IR (HOMA-IR) but were not associated with obesity. We used this GWAS to employ hazard ratios (HRs) for breast cancer risk by adjusting for potential confounding factors. Results: No direct association was observed between comprising 12 IR genetic instruments and breast cancer risk (HR = 0.93, 95% CI: 0.76–1.14). In phenotype-specific analysis, genetically elevated FG was associated with reduced risk for breast cancer (main contributor of this MR-effect estimate: G6PC2 rs13431652; HR = 0.59, 95% CI: 0.35–0.99). Genetically driven FI and HOMA-IR were not significantly associated. Stratification analyses by body mass index, exercise, and dietary fat intake with combined phenotypes showed that genetically elevated IR was associated with greater breast cancer risk in overall obesity and inactive subgroups (single contributor: MTRR/LOC729506 rs13188458; HR = 2.21, 95% CI: 1.03–4.75). Conclusions: We found complex evidence for causal association between IR and risk of breast cancer, which may support the potential value of intervention trials to lower IR and reduce breast cancer risk. [ABSTRACT FROM AUTHOR]

Background: PPM1D (Protein phosphatase magnesium-dependent 1δ) is known as a damage response regulator, a part of the p53 negative feedback loop. Truncating mutations of PPM1D, resulting in overexpression, are frequently found in the blood of patients with breast or ovarian cancer. To identify whether the PPM1D mutation predisposes patients to such cancers or if it results from the cancer and therapy, somatic PPM1D mutations in association with previous cancer and chemotherapy need to be explored. Methods: We performed next-generation sequencing (NGS) analysis of blood samples from patients suspected to have hereditary cancer. We grouped the patients according to their diagnoses and history of chemotherapy. For the patients with PPM1D mutations in blood, tumor tissue specimens were examined for the PPM1D mutation using conventional sequencing. Results: A total of 1,195 patients, including 719 patients with breast cancer and 240 with ovarian cancer, were tested, and four (~0.3%) had the truncating mutation in PPM1D. All truncating mutations were in exon 6, in mosaic form, with a mean allele fraction of 11.15%. While 395 out of the 1,195 patients had undergone chemotherapy, the four with the truncating mutation had a history of cisplatin-based chemotherapy. No corresponding mutations were identified in the tumor tissues. Conclusions: We investigated the frequency of the somatic mosaic PPM1D mutation, in patients with breast or ovarian cancer, which is suggested to be low and related to a history of cisplatin-based chemotherapy. It may be a marker of previous exposure to selective pressure for cells with an impaired DNA damage response. [ABSTRACT FROM AUTHOR]

In order to further our understanding of how gene expression contributes to key functional properties of neurons, we combined publicly accessible gene expression, electrophysiology, and morphology measurements to identify cross-cell type correlations between these data modalities. Building on our previous work using a similar approach, we distinguished between correlations which were “class-driven,” meaning those that could be explained by differences between excitatory and inhibitory cell classes, and those that reflected graded phenotypic differences within classes. Taking cell class identity into account increased the degree to which our results replicated in an independent dataset as well as their correspondence with known modes of ion channel function based on the literature. We also found a smaller set of genes whose relationships to electrophysiological or morphological properties appear to be specific to either excitatory or inhibitory cell types. Next, using data from PatchSeq experiments, allowing simultaneous single-cell characterization of gene expression and electrophysiology, we found that some of the gene-property correlations observed across cell types were further predictive of within-cell type heterogeneity. In summary, we have identified a number of relationships between gene expression, electrophysiology, and morphology that provide testable hypotheses for future studies. [ABSTRACT FROM AUTHOR]

Background: Systematic reviews of randomised controlled trials (RCTs) have suggested that maternal vitamin D (25[OH]D) and calcium supplementation increase birth weight. However, limitations of many trials were highlighted in the reviews. Our aim was to combine genetic and RCT data to estimate causal effects of these two maternal traits on offspring birth weight.Methods and Findings: We performed two-sample mendelian randomisation (MR) using genetic instrumental variables associated with 25(OH)D and calcium that had been identified in genome-wide association studies (GWAS; sample 1; N = 122,123 for 25[OH]D and N = 61,275 for calcium). Associations between these maternal genetic variants and offspring birth weight were calculated in the UK Biobank (UKB) (sample 2; N = 190,406). We used data on mother-child pairs from two United Kingdom birth cohorts (combined N = 5,223) in sensitivity analyses to check whether results were influenced by fetal genotype, which is correlated with the maternal genotype (r ≈ 0.5). Further sensitivity analyses to test the reliability of the results included MR-Egger, weighted-median estimator, 'leave-one-out', and multivariable MR analyses. We triangulated MR results with those from RCTs, in which we used randomisation to supplementation with vitamin D (24 RCTs, combined N = 5,276) and calcium (6 RCTs, combined N = 543) as an instrumental variable to determine the effects of 25(OH)D and calcium on birth weight. In the main MR analysis, there was no strong evidence of an effect of maternal 25(OH)D on birth weight (difference in mean birth weight -0.03 g [95% CI -2.48 to 2.42 g, p = 0.981] per 10% higher maternal 25[OH]D). The effect estimate was consistent across our MR sensitivity analyses. Instrumental variable analyses applied to RCTs suggested a weak positive causal effect (5.94 g [95% CI 2.15-9.73, p = 0.002] per 10% higher maternal 25[OH]D), but this result may be exaggerated because of risk of bias in the included RCTs. The main MR analysis for maternal calcium also suggested no strong evidence of an effect on birth weight (-20 g [95% CI -44 to 5 g, p = 0.116] per 1 SD higher maternal calcium level). Some sensitivity analyses suggested that the genetic instrument for calcium was associated with birth weight via exposures that are independent of calcium levels (horizontal pleiotropy). Application of instrumental variable analyses to RCTs suggested that calcium has a substantial effect on birth weight (178 g [95% CI 121-236 g, p = 1.43 × 10-9] per 1 SD higher maternal calcium level) that was not consistent with any of the MR results. However, the RCT instrumental variable estimate may have been exaggerated because of risk of bias in the included RCTs. Other study limitations include the low response rate of UK Biobank, which may bias MR estimates, and the lack of suitable data to test whether the effects of genetic instruments on maternal calcium levels during pregnancy were the same as those outside of pregnancy.Conclusions: Our results suggest that maternal circulating 25(OH)D does not influence birth weight in otherwise healthy newborns. However, the effect of maternal circulating calcium on birth weight is unclear and requires further exploration with more research including RCT and/or MR analyses with more valid instruments. [ABSTRACT FROM AUTHOR]

Objective: An early onset of menarche and, later, menopause are well-established risk factors for the development of breast cancer and endometrial cancer. Although the largest GWASs have identified 389 independent signals for age at menarche (AAM) and 44 regions for age at menopause (ANM), GWAS can only identify the associations between variants and traits. The aim of this study was to identify genes whose expression levels were associated with AAM or ANM due to pleiotropy or causality by integrating GWAS data with genome-wide expression quantitative trait loci (eQTLs) data. We also aimed to identify the pleiotropic genes that influenced both phenotypes. Method: We employed GWAS data of AAM and ANM and genome-wide eQTL data from whole blood. The summary data-based Mendelian randomization method was used to prioritize the associated genes for further study. The colocalization analysis was used to identify the pleiotropic genes associated with both phenotypes. Results: We identified 31 genes whose expression was associated with AAM and 24 genes whose expression was associated with ANM due to pleiotropy or causality. Two pleiotropic genes were identified to be associated with both phenotypes. Conclusion: The results point out the most possible genes which were responsible for the association. Our study prioritizes the associated genes for further functional mechanistic study of AAM and ANM and illustrates the benefit of integrating different omics data into the study of complex traits. [ABSTRACT FROM AUTHOR]

The largest ever recorded epidemic of the Chikungunya virus (CHIKV) broke out in 2004 and affected four continents. Acute symptomatic infections are typically associated with the onset of fever and often debilitating polyarthralgia/polyarthritis. In this study, a systems biology approach was adopted to analyze the blood transcriptomes of adults acutely infected with the CHIKV. Gene signatures that were associated with viral RNA levels and the onset of symptoms were identified. Among these genes, the putative role of the Eukaryotic Initiation Factor (eIF) family genes and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC3A) in the CHIKV replication process were displayed. We further compared these signatures with signatures induced by the Dengue virus infection and rheumatoid arthritis. Finally, we demonstrated that the CHIKV in vitro infection of murine bone marrow-derived macrophages induced IL-1 beta production in a mechanism that is significantly dependent on the inflammasome NLRP3 activation. The observations provided valuable insights into virus-host interactions during the acute phase and can be instrumental in the investigation of new and effective therapeutic interventions. [ABSTRACT FROM AUTHOR]

Chronic pain is highly prevalent worldwide and represents a significant socioeconomic and public health burden. Several aspects of chronic pain, for example back pain and a severity-related phenotype ‘chronic pain grade’, have been shown previously to be complex heritable traits with a polygenic component. Additional pain-related phenotypes capturing aspects of an individual’s overall sensitivity to experiencing and reporting chronic pain have also been suggested as a focus for investigation. We made use of a measure of the number of sites of chronic pain in individuals within the UK general population. This measure, termed Multisite Chronic Pain (MCP), is a complex trait and its genetic architecture has not previously been investigated. To address this, we carried out a large-scale genome-wide association study (GWAS) of MCP in ~380,000 UK Biobank participants. Our findings were consistent with MCP having a significant polygenic component, with a Single Nucleotide Polymorphism (SNP) heritability of 10.2%. In total 76 independent lead SNPs at 39 risk loci were associated with MCP. Additional gene-level association analyses identified neurogenesis, synaptic plasticity, nervous system development, cell-cycle progression and apoptosis genes as enriched for genetic association with MCP. Genetic correlations were observed between MCP and a range of psychiatric, autoimmune and anthropometric traits, including major depressive disorder (MDD), asthma and Body Mass Index (BMI). Furthermore, in Mendelian randomisation (MR) analyses a causal effect of MCP on MDD was observed. Additionally, a polygenic risk score (PRS) for MCP was found to significantly predict chronic widespread pain (pain all over the body), indicating the existence of genetic variants contributing to both of these pain phenotypes. Overall, our findings support the proposition that chronic pain involves a strong nervous system component with implications for our understanding of the physiology of chronic pain. These discoveries may also inform the future development of novel treatment approaches. [ABSTRACT FROM AUTHOR]

The nosocomial pathogen Acinetobacter baumannii acquired clinical significance due to the rapid development of its multi-drug resistant (MDR) phenotype. A. baumannii strains have the ability to colonize several ecological niches including soil, water, and animals, including humans. They also survive under extremely harsh environmental conditions thriving on rare and recalcitrant carbon compounds. However, the molecular basis behind such extreme adaptability of A. baumannii is unknown. We have therefore determined the complete genome sequence of A. baumannii DS002, which was isolated from agricultural soils, and compared it with 78 complete genome sequences of A. baumannii strains having complete information on the source of their isolation. Interestingly, the genome of A. baumannii DS002 showed high similarity to the genome of A. baumannii SDF isolated from the body louse. The environmental and clinical strains, which do not share a monophyletic origin, showed the existence of a strain-specific unique gene pool that supports niche-specific survival. The strains isolated from infected samples contained a genetic repertoire with a unique gene pool coding for iron acquisition machinery, particularly those required for the biosynthesis of acinetobactin. Interestingly, these strains also contained genes required for biofilm formation. However, such gene sets were either partially or completely missing in the environmental isolates, which instead harbored genes required for alternate carbon catabolism and a TonB-dependent transport system involved in the acquisition of iron via siderophores or xenosiderophores. [ABSTRACT FROM AUTHOR]

In this study, we perform a full genome-wide association study (GWAS) to identify statistically significantly associated single nucleotide polymorphisms (SNPs) with three red blood cell (RBC) components and follow it with two independent PheWASs to examine associations between phenotypic data (case-control status of diagnoses or disease), significant SNPs, and RBC component levels. We first identified associations between the three RBC components: mean platelet volume (MPV), mean corpuscular volume (MCV), and platelet counts (PC), and the genotypes of approximately 500,000 SNPs on the Illumina Infimum DNA Human OmniExpress-24 BeadChip using a single cohort of 4,673 Northern Nevadans. Twenty-one SNPs in five major genomic regions were found to be statistically significantly associated with MPV, two regions with MCV, and one region with PC, with p<5x10-8. Twenty-nine SNPs and nine chromosomal regions were identified in 30 previous GWASs, with effect sizes of similar magnitude and direction as found in our cohort. The two strongest associations were SNP rs1354034 with MPV (p = 2.4x10-13) and rs855791 with MCV (p = 5.2x10-12). We then examined possible associations between these significant SNPs and incidence of 1,488 phenotype groups mapped from International Classification of Disease version 9 and 10 (ICD9 and ICD10) codes collected in the extensive electronic health record (EHR) database associated with Healthy Nevada Project consented participants. Further leveraging data collected in the EHR, we performed an additional PheWAS to identify associations between continuous red blood cell (RBC) component measures and incidence of specific diagnoses. The first PheWAS illuminated whether SNPs associated with RBC components in our cohort were linked with other hematologic phenotypic diagnoses or diagnoses of other nature. Although no SNPs from our GWAS were identified as strongly associated to other phenotypic components, a number of associations were identified with p-values ranging between 1x10-3 and 1x10-4 with traits such as respiratory failure, sleep disorders, hypoglycemia, hyperglyceridemia, GERD and IBS. The second PheWAS examined possible phenotypic predictors of abnormal RBC component measures: a number of hematologic phenotypes such as thrombocytopenia, anemias, hemoglobinopathies and pancytopenia were found to be strongly associated to RBC component measures; additional phenotypes such as (morbid) obesity, malaise and fatigue, alcoholism, and cirrhosis were also identified to be possible predictors of RBC component measures. [ABSTRACT FROM AUTHOR]

Background: Ischemic retinopathies (IRs) are leading causes of visual impairment. They are characterized by an initial phase of microvascular degeneration and a second phase of aberrant pre-retinal neovascularization (NV). microRNAs (miRNAs) regulate gene expression, and a number play a role in normal and pathological NV. But, post-transcriptional modulation of miRNAs in the eye during the development of IRs has not been systematically evaluated. Aims & methods: Using Next Generation Sequencing (NGS) we profiled miRNA expression in the retina and choroid during vasodegenerative and NV phases of oxygen-induced retinopathy (OIR). Results: Approximately 20% of total miRNAs exhibited altered expression (up- or down-regulation); 6% of miRNA were found highly expressed in retina and choroid of rats subjected to OIR. During OIR-induced vessel degeneration phase, miR-199a-3p, -199a-5p, -1b, -126a-3p displayed a robust decreased expression (> 85%) in the retina. While in the choroid, miR-152-3p, -142-3p, -148a-3p, -532-3p were upregulated (>200%) and miR-96-5p, -124-3p, -9a-3p, -190b-5p, -181a-1-3p, -9a-5p, -183-5p were downregulated (>70%) compared to controls. During peak pathological NV, miR-30a-5p, -30e-5p and 190b-5p were markedly reduced (>70%), and miR-30e-3p, miR-335, -30b-5p strongly augmented (by up to 300%) in the retina. Whereas in choroid, miR-let-7f-5p, miR-126a-5p and miR-101a-3p were downregulated by (>81%), and miR-125a-5p, let-7e-5p and let-7g-5p were upregulated by (>570%) during NV. Changes in miRNA observed using NGS were validated using qRT-PCR for the 24 most modulated miRNAs. In silico approach to predict miRNA target genes (using algorithms of miRSystem database) identified potential new target genes with pro-inflammatory, apoptotic and angiogenic properties. Conclusion: The present study is the first comprehensive description of retinal/choroidal miRNAs profiling in OIR (using NGS technology). Our results provide a valuable framework for the characterization and possible therapeutic potential of specific miRNAs involved in ocular IR-triggered inflammation, angiogenesis and degeneration. [ABSTRACT FROM AUTHOR]

Loin muscle area (LMA) and loin muscle depth (LMD) are important traits influencing the production performance of breeding pigs. However, the genetic architecture of these two traits is still poorly understood. To discern the genetic architecture of LMA and LMD, a material consisting of 6043 Duroc pigs belonging to two populations with different genetic backgrounds was collected and applied in genome-wide association studies (GWAS) with a genome-wide distributed panel of 50K single nucleotide polymorphisms (SNPs). To improve the power of detection for common SNPs, we conducted a meta-analysis in these two pig populations and uncovered additional significant SNPs. As a result, we identified 75 significant SNPs for LMA and LMD on SSC6, 7, 12, 16, and 18. Among them, 25 common SNPs were associated with LMA and LMD. One pleiotropic quantitative trait locus (QTL), which was located on SSC7 with a 283 kb interval, was identified to affect LMA and LMD. Marker ALGA0040260 is a key SNP for this QTL, explained 1.77% and 2.48% of the phenotypic variance for LMA and LMD, respectively. Another genetic region on SSC16 (709 kb) was detected and displayed prominent association with LMA and the peak SNP, WU_10.2_16_35829257, contributed 1.83% of the phenotypic variance for LMA. Further bioinformatics analysis determined eight promising candidate genes (GCLC, GPX8, DAXX, FGF21, TAF11, SPDEF, NUDT3, and PACSIN1) with functions in glutathione metabolism, adipose and muscle tissues development and lipid metabolism. This study provides the first GWAS for the LMA and LMD of Duroc breed to analyze the underlying genetic variants through a large sample size. The findings further advance our understanding and help elucidate the genetic architecture of LMA, LMD and growth-related traits in pigs. [ABSTRACT FROM AUTHOR]

Objective: The special status accorded to elite athletes stems from their uncommon genetic potential to excel in world-class power sports (PS). Genetic polymorphisms have been reported to influence elite PS status. Reports of associations between the α-actinin-3 gene (ACTN3) R577X polymorphism and PS have been inconsistent. In light of new published studies, we perform a meta-analysis to further explore the roles of this polymorphism in PS performance among elite athletes. Methods: Multi-database literature search yielded 44 studies from 38 articles. Pooled odds ratios (ORs) and 95% confidence intervals (CIs) were used in estimating associations (significance threshold was set at Pa ≤ 0.05) using the allele-genotype model (R and X alleles, RX genotype). Outlier analysis was used to examine its impact on association and heterogeneity outcomes. Subgroup analysis was race (Western and Asian) and gender (male/female)-based. Interaction tests were applied to differential outcomes between the subgroups, P-values of which were Bonferroni corrected (Pinteraction BC). Tests for sensitivity and publication bias were performed. Results: Significant overall R allele effects (OR 1.21, 95% CI 1.07–1.37, Pa = 0.002) were confirmed in the Western subgroup (OR 1.11, 95% CI 1.01–1.22, Pa = 0.02) and with outlier treatment (ORs 1.12–1.20, 95% CIs 1.02–1.30, Pa < 10−5–0.01). This treatment resulted in acquired significance of the RX effect in Asian athletes (OR 1.91, 95% CI 1.25–2.92, Pa = 0.003). Gender analysis dichotomized the RX genotype and R allele effects as significantly higher in male (OR 1.14, 95% CI 1.02–1.28, Pa = 0.02) and female (OR 1.58, 95% CI 1.21–2.06, Pa = 0.0009) athletes, respectively, when compared with controls. Significant R female association was improved with a test of interaction (Pinteraction BC = 0.03). The overall, Asian and female outcomes were robust. The R allele results were more robust than the RX genotype outcomes. No evidence of publication bias was found. Conclusions: In this meta-analysis, we present clear associations between the R allele/RX genotype in the ACTN3 polymorphism and elite power athlete status. Significant effects of the R allele (overall analysis, Western and female subgroups) and RX genotype (Asians and males) were for the most part, results of outlier treatment. Interaction analysis improved the female outcome. These robust findings were free of publication bias. [ABSTRACT FROM AUTHOR]

Backgrounds: A recent GWAS Study found a new locus (rs9810888 in CACNA1D) was associated with blood pressure (BP) in Chinese adults. But whether the association exists in children is unknown. Whether lifestyle behaviors could interact with rs9810888 on BP is not clear. This study aimed to identify the association between rs9810888 and BP in children, and also explore the gene-lifestyle interaction. Methods: A case-control study was conducted among 2030 Chinese children aged 7 to 18 years. Genotyping was conducted by using the matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Lifestyle behaviors were investigated with questionnaire. Results: With adjustment for age, age square, sex, study group and body mass index (BMI), rs9810888 was significantly associated with diastolic BP (DBP) (b = 1.69, p = 0.021) and mean arterial BP (MAP) (b = 1.56, p = 0.010). Stratified analysis showed that the rs9810888 GG genotype carriers had higher DBP than GT/TT carriers (b = 3.78, p = 0.023) in the subgroup having protein intake (meat/fish/soybeans/egg)

Several horse breeds have been specifically selected for the ability to exhibit alternative patterns of locomotion, or gaits. A premature stop codon in the gene DMRT3 is permissive for “gaitedness” across breeds. However, this mutation is nearly fixed in both American Standardbred trotters and pacers, which perform a diagonal and lateral gait, respectively, during harness racing. This suggests that modifying alleles must influence the preferred gait at racing speeds in these populations. A genome-wide association analysis for the ability to pace was performed in 542 Standardbred horses (n = 176 pacers, n = 366 trotters) with genotype data imputed to ~74,000 single nucleotide polymorphisms (SNPs). Nineteen SNPs on nine chromosomes (ECA1, 2, 6, 9, 17, 19, 23, 25, 31) reached genome-wide significance (p < 1.44 x 10−6). Variant discovery in regions of interest was carried out via whole-genome sequencing. A set of 303 variants from 22 chromosomes with putative modifying effects on gait was genotyped in 659 Standardbreds (n = 231 pacers, n = 428 trotters) using a high-throughput assay. Random forest classification analysis resulted in an out-of-box error rate of 0.61%. A conditional inference tree algorithm containing seven SNPs predicted status as a pacer or trotter with 99.1% accuracy and subsequently performed with 99.4% accuracy in an independently sampled population of 166 Standardbreds (n = 83 pacers, n = 83 trotters). This highly accurate algorithm could be used by owners/trainers to identify Standardbred horses with the potential to race as pacers or as trotters, according to the genotype identified, prior to initiating training and would enable fine-tuning of breeding programs with designed matings. Additional work is needed to determine both the algorithm’s utility in other gaited breeds and whether any of the predictive SNPs play a physiologically functional role in the tendency to pace or tag true functional alleles. [ABSTRACT FROM AUTHOR]

The African spiny mouse, Acomys spp., is capable of scar-free dermal wound healing. Here, we have performed a comprehensive analysis of gene expression throughout wound healing following full-thickness excisional dermal wounds in both Acomys cahirinus and Mus musculus. Additionally, we provide an annotated, de novo transcriptome assembly of A. cahirinus skin and skin wounds. Using a novel computational comparative RNA-Seq approach along with pathway and co-expression analyses, we identify enrichment of regeneration associated genes as well as upregulation of genes directly related to muscle development or function. Our RT-qPCR data reveals induction of the myogenic regulatory factors, as well as upregulation of embryonic myosin, starting between days 14 and 18 post-wounding in A. cahirinus. In contrast, the myogenic regulatory factors remain downregulated, embryonic myosin is only modestly upregulated, and no new muscle fibers of the panniculus carnosus are generated in M. musculus wounds. Additionally, we show that Col6a1, a key component of the satellite cell niche, is upregulated in A. cahirinus compared to M. musculus. Our data also demonstrate that the macrophage profile and inflammatory response is different between species, with A. cahirinus expressing significantly higher levels of Il10. We also demonstrate differential expression of the upstream regulators Wnt7a, Wnt2 and Wnt6 during wound healing. Our analyses demonstrate that A. cahirinus is capable of de novo skeletal muscle regeneration of the panniculus carnosus following removal of the extracellular matrix. We believe this study represents the first detailed analysis of de novo skeletal muscle regeneration observed in an adult mammal. [ABSTRACT FROM AUTHOR]

CD4+ effector/memory T cells (Tem) represent a leading edge of the adaptive immune system responsible for protecting the body from infection, cancer, and other damaging processes. However, a subset of Tem cells with low expression of CD45Rb (RbLoTem) has been shown to suppress inflammation despite their effector surface phenotype and the lack of FoxP3 expression, the canonical transcription factor found in most regulatory T cells. In this report, we show that RbLoTem cells can suppress inflammation by influencing Treg behavior. Co-culturing activated RbLoTem and Tregs induced high expression of IL-10 in vitro, and conditioned media from RbLoTem cells induced IL-10 expression in FoxP3+ Tregs in vitro and in vivo, indicating that RbLoTem cells communicate with Tregs in a cell-contact independent fashion. Transcriptomic and multi-analyte Luminex data identified both IL-2 and IL-4 as potential mediators of RbLoTem-Treg communication, and antibody-mediated neutralization of either IL-4 or CD124 (IL-4Rα) prevented IL-10 induction in Tregs. Moreover, isolated Tregs cultured with recombinant IL-2 and IL-4 strongly induced IL-10 production. Using house dust mite (HDM)-induced airway inflammation as a model, we confirmed that the in vivo suppressive activity of RbLoTem cells was lost in IL-4-ablated RbLoTem cells. These data support a model in which RbLoTem cells communicate with Tregs using a combination of IL-2 and IL-4 to induce robust expression of IL-10 and suppression of inflammation. [ABSTRACT FROM AUTHOR]

The mechanisms underlying the transition from acute to chronic pain remain unclear. Here, we sought to characterize the transcriptome associated with chronic low back pain as well as the transcriptome of the transition from acute to chronic low back pain. For the analysis, we compared the whole blood transcriptome of: (a) patients at the onset of low back pain who no longer had pain within 6 weeks after onset (acute) with patients who developed chronic low back pain at 6 months (chronic T5); and, (b) patients at the onset of low back pain (chronic T1) who developed chronic pain at 6 months with healthy pain-free (normal) controls. The majority of differentially expressed genes were protein coding. We illustrate a unique chronic low back pain transcriptome characterized by significant enrichment for known pain genes, extracellular matrix genes, and genes from the extended major histocompatibility complex (MHC) genomic locus. The transcriptome of the transition from acute to chronic low back pain was characterized by significant upregulation of antigen presentation pathway (MHC class I and II) genes and downregulation of mitochondrial genes associated with oxidative phosphorylation, suggesting a unique genomic signature of vulnerability to low back pain chronicity. [ABSTRACT FROM AUTHOR]

The genetic programs responsible for pulmonary lymphatic maturation prior to birth are not known. To address this gap in knowledge, we developed a novel cell sorting strategy to collect fetal pulmonary lymphatic endothelial cells (PLECs) for global transcriptional profiling. We identified PLECs based on their unique cell surface immunophenotype (CD31+/Vegfr3+/Lyve1+/Pdpn+) and isolated them from murine lungs during late gestation (E16.5, E17.5, E18.5). Gene expression profiling was performed using whole-genome microarrays, and 1,281 genes were significantly differentially expressed with respect to time (FDR q < 0.05) and grouped into six clusters. Two clusters containing a total of 493 genes strongly upregulated at E18.5 were significantly enriched in genes with functional annotations corresponding to innate immune response, positive regulation of angiogenesis, complement & coagulation cascade, ECM/cell-adhesion, and lipid metabolism. Gene Set Enrichment Analysis identified several pathways coordinately upregulated during late gestation, the strongest of which was the type-I IFN-α/β signaling pathway. Upregulation of canonical interferon target genes was confirmed by qRT-PCR and in situ hybridization in E18.5 PLECs. We also identified transcriptional events consistent with a prenatal PLEC maturation program. This PLEC-specific program included individual genes (Ch25h, Itpkc, Pcdhac2 and S1pr3) as well as a set of chemokines and genes containing an NF-κB binding site in their promoter. Overall, this work reveals transcriptional insights into the genes, signaling pathways and biological processes associated with pulmonary lymphatic maturation in the fetal lung. [ABSTRACT FROM AUTHOR]

Broodiness, one of the maternal behaviors and instincts for natural breeding in birds, is an interesting topic in reproductive biology. Broodiness in poultry is characterized by persistent nesting, usually associated with cessation of egg laying. The study of avian broodiness is essential for bird conservation breeding and commercial poultry industry. In this study, we examined the hypothalamus transcriptome of Muscovy duck in three reproductive stages, including egg-laying anaphase (LA), brooding prophase (BP) and brooding metaphase (BM). Differences in gene expression during the transition from egg-laying to broodiness were examined, and 155, 379, 292 differently expressed genes (DEGs) were obtained by pairwise comparisons of LA-vs-BP, LA-vs-BM and BP-vs-BM, respectively (fold change≥1.5, P < 0.05). Gene Ontology Term (GO) enrichment analysis suggested a possible role of oxidative stress in the hypothalamus might invoke reproductive costs that potentially change genes expression. KEGG analysis revealed glutamatergic synapse, dopaminergic synapse, serotonergic synapse and GABAergic synapse pathway were significantly enriched, and regulator genes were identified. Eight gene expression patterns were illustrated by trend analysis and further clustered into three clusters. Additional six hub genes were identified through combining trend analysis and protein-protein interaction (PPI) analysis. Our results suggested that the cyclical mechanisms of reproductive function conversion include effects of oxidative stress, biosynthesis of neurotransmitters or their receptors, and interactions between glucocorticoids and thyroid hormones and regulatory genes. These candidate genes and biological pathways may be used as targets for artificial manipulation and marker-assisted breeding in the reproductive behavior. [ABSTRACT FROM AUTHOR]

The Philadelphia chromosome-negative myeloproliferative neoplasms (MPN) share similar molecular characteristics in that they frequently harbor hotspot mutations in JAK2, CALR or MPL, leading to activated JAK/STAT signaling. However, these MPN have distinct symptoms, morphology, and natural histories, including different tendencies to progress to fibrosis. Although the role of inflammation in tissue fibrosis is well recognized, inflammatory gene expression in bone marrows involved by MPN has been understudied. We analyzed the expression of inflammatory genes by directly measuring RNA transcript abundance in bone marrow biopsies of 108 MPN patients. Unsupervised analyses identified gene expression patterns that distinguish prefibrotic (grade 1–2) MPN from overtly fibrotic (grade 2–3) MPN with high sensitivity and specificity in two independent cohorts. Furthermore, prefibrotic and overtly fibrotic MPN are separable into subsets with different activities in biological pathways linked to inflammation, including cytokines, chemokines, interferon response, and toll-like receptor signaling. In conclusion, this study demonstrates that MPN with overt fibrosis is associated with significant inflammatory gene upregulation in the bone marrow, revealing potential roles for multiple pro-inflammatory signaling networks in the development of myelofibrosis and suggesting potential therapeutic mechanisms to alleviate this process in the bone marrow microenvironment. [ABSTRACT FROM AUTHOR]

Sertoli cell-only (SCO) syndrome is a severe form of human male infertility seemingly characterized by the lack all spermatogenic cells. However, tubules of some SCO testes contain small patches of active spermatogenesis and thus spermatogonial stem cells. We hypothesized that these stem cells cannot replicate and seed spermatogenesis in barren areas of tubule because as-of-yet unrecognized deficits in Sertoli cell gene expression disable most stem cell niches. Performing the first thorough comparison of the transcriptomes of human testes exhibiting complete spermatogenesis with the transcriptomes of testes with SCO syndrome, we defined transcripts that are both predominantly expressed by Sertoli cells and expressed at aberrant levels in SCO testes. Some of these transcripts encode proteins required for the proper assembly of adherent and gap junctions at sites of contact with other cells, including spermatogonial stem cells (SSCs). Other transcripts encode GDNF, FGF8 and BMP4, known regulators of mouse SSCs. Thus, most SCO Sertoli cells can neither organize junctions at normal sites of cell-cell contact nor stimulate SSCs with adequate levels of growth factors. We propose that the critical deficits in Sertoli cell gene expression we have identified contribute to the inability of spermatogonial stem cells within small patches of spermatogenesis in some SCO testes to seed spermatogenesis to adjacent areas of tubule that are barren of spermatogenesis. Furthermore, we predict that one or more of these deficits in gene expression are primary causes of human SCO syndrome. [ABSTRACT FROM AUTHOR]

Barnacle settlement involves sensing of a variety of exogenous cues. A pair of antennules is the main sensory organ that the cyprid larva uses to explore the surface. Antennules are equipped with a number of setae that have both chemo- and mechanosensing function. The current study explores the repertoire of sensory receptors in Balanus improvisus cyprid antennules with the goal to better understand sensory systems involved in the settling behavior of this species. We carried out transcriptome sequencing of dissected B. improvisus cyprid antennules. The generated transcriptome assembly was used to search for sensory receptors using HMM models. Among potential chemosensory genes, we identified the ionotropic receptors IR25a, IR8a and IR93a, and several divergent IR candidates to be expressed in the cyprid antennules. We found one gustatory-like receptor but no odorant receptors, chemosensory or odorant-binding proteins. Apart from chemosensory receptors, we also identified 13 potential mechanosensory genes represented by several transient receptor potential channels (TRP) subfamilies. Furthermore, we analyzed changes in expression profiles of IRs and TRPs during the B. improvisus settling process. Several of the sensory genes were differentially expressed during the course of larval settlement. This study gives expanded knowledge about the sensory systems present in barnacles, a taxonomic group for which only limited information about receptors is currently available. It furthermore serves as a starting point for more in depth studies of how sensory signaling affects settling behavior in barnacles with implications for preventing biofouling. [ABSTRACT FROM AUTHOR]

Glycine receptors (GlyRs) are ligand-gated chloride channels mediating inhibitory neurotransmission in the brain stem and spinal cord. They function as pentamers composed of alpha and beta subunits for which 5 genes have been identified in human (GLRA1, GLRA2, GLRA3, GLRA4, GLRB). Several in vitro studies showed that the pentameric subtype composition as well as its stoichiometry influence the distribution and the molecular function of the receptor. Moreover, mutations in some of these genes are involved in different human conditions ranging from tinnitus to epilepsy and hyperekplexia, suggesting distinct functions of the different subunits. Although the beta subunit is essential for synaptic clustering of the receptor, the specific role of each alpha subtype is still puzzling in vivo. The zebrafish genome encodes for five glycine receptor alpha subunits (glra1, glra2, glra3, glra4a, glra4b) thus offering a model of choice to investigate the respective role of each subtype on general motor behaviour. After establishing a phylogeny of GlyR subunit evolution between human and zebrafish, we checked the temporal expression pattern of these transcripts during embryo development. Interestingly, we found that glra1 is the only maternally transmitted alpha subunit. We also showed that the expression of the different GlyR subunits starts at different time points during development. Lastly, in order to decipher the role of each alpha subunit on the general motor behaviour of the fish, we knocked out individually each alpha subunit by CRISPR/Cas9-targeted mutagenesis. Surprisingly, we found that knocking out any of the alpha2, 3, a4a or a4b subunit did not lead to any obvious developmental or motor phenotype. However, glra1-/- (hitch) embryos depicted a strong motor dysfunction from 3 days, making them incapable to swim and thus leading to their premature death. Our results infer a strong functional redundancy between alpha subunits and confirm the central role played by glra1 for proper inhibitory neurotransmission controlling locomotion. The genetic tools we developed here will be of general interest for further studies aiming at dissecting the role of GlyRs in glycinergic transmission in vivo and the hitch mutant (hic) is of specific relevance as a new model of hyperekplexia. [ABSTRACT FROM AUTHOR]

For most complex traits, the majority of SNPs identified through genome-wide association studies (GWAS) reside within noncoding regions that have no known function. However, these regions are enriched for the regulatory enhancers specific to the cells relevant to the specific trait. Indeed, many of the GWAS loci that have been functionally characterized lie within enhancers that regulate expression levels of key genes. In order to identify polymorphisms with potential allele-specific regulatory effects, we developed a bioinformatics pipeline that harnesses epigenetic signatures as well as transcription factor (TF) binding motifs to identify putative enhancers containing a SNP with potential allele-specific TF binding in linkage disequilibrium (LD) with a GWAS-identified SNP. We applied the approach to GWAS findings for blood lipids, revealing 7 putative enhancers harboring associated SNPs, 3 of which lie within the introns of LCAT and ABCA1, genes that play crucial roles in cholesterol biogenesis and lipoprotein metabolism. All 3 enhancers demonstrated allele-specific in vitro regulatory activity in liver-derived cell lines. We demonstrated that these putative enhancers are in close physical proximity to the promoters of their respective genes, in situ, likely through chromatin looping. In addition, the associated alleles altered the likelihood of transcription activator STAT3 binding. Our results demonstrate that through our approach, the LD blocks that contain GWAS signals, often hundreds of kilobases in size with multiple SNPs serving as statistical proxies to the true functional site, can provide an experimentally testable hypothesis for the underlying regulatory mechanism linking genetic variants to complex traits. [ABSTRACT FROM AUTHOR]

Mosquitoes act as vectors of numerous pathogens that cause human diseases. Dengue virus (DENV) transmitted by mosquito, Aedes aegypti, is responsible for dengue fever epidemics worldwide with a serious impact on human health. Currently, disease control mainly relies on vector targeted intervention strategies. Therefore, it is imperative to understand the molecular mechanisms underlying the innate immune response of mosquitoes against pathogens. In the present study, the expression profiles of immunity-related genes in the midgut responding to DENV infection by feeding were analyzed by transcriptome and quantitative real-time PCR. The level of Antimicrobial peptides (AMPs) increased seven days post-infection (d.p.i.), which could be induced by the Toll immune pathway. The expression of reactive oxygen species (ROS) genes, including antioxidant genes, such as HPX7, HPX8A, HPX8B, HPX8C were induced at one d.p.i. and peaked again at ten d.p.i. in the midgut. Interestingly, down-regulation of the antioxidant gene HPX8C by RNA interference led to reduction in the virus titer in the mosquito, probably due to the elevated levels of ROS. Application of a ROS inhibitor and scavenger molecules further established the role of oxygen free radicals in the modulation of the immune response to DENV infection. Overall, our comparative transcriptome analyses provide valuable information about the regulation of immunity related genes in the transmission vector in response to DENV infection. It further allows us to identify novel molecular mechanisms underlying the host-virus interaction, which might aid in the development of novel strategies to control mosquito-borne diseases. [ABSTRACT FROM AUTHOR]

Functional studies suggest that promoter polymorphisms of the Prostaglandin D Receptor (PTGDR) gene can be involved in asthma. All-trans Retinoic acid (ATRA) has also been linked to allergic diseases. We have previously described the PTGDR promoter activation mediated by ATRA through response elements (RARE) at position -549T> C. In this study we aimed to analyze the effect of retinoic acid (RA) on the expression of PTGDR, the production of cytokines as well as to evaluate the binding of RA receptors to RA-Response Elements (RARE) sequences. A549 cells were transfected with vectors carrying different PTGDR haplotypes and treated with all-Trans Retinoic Acid (ATRA). PTGDR expression was measured by qPCR. Chromatin Immunoprecipitation assays (ChIP) were performed in ATRA stimulated KU812 cells and in PBMCs of patients carrying CTCT, CCCC or CCCT haplotypes. In addition, a broad panel of cytokines was analyzed by cytometric bead assay in A549 cells. The expression of PTGDR increased in A549 cells transfected with PTGDR-variants. The CCCC haplotype showed a significantly higher expression compared with CTCT. However, we found that RA up-regulated PTGDR expression through RARα mainly in the CTCT variant. Experiments on PBMCs from allergic patients carrying the -549T and -549C variant of the PTGDR promoter after ATRA and RAR antagonist administration confirmed the modulation of PTGDR by ATRA. The cytokine analysis showed that IL4 and IL6 levels were significantly increased in A549 cells transfected with PTGDR. In addition, ATRA treatment decreased the levels of IL4, IL6 and TNFα in A549 cells, whereas it increased IL4 and TNFα levels in PTGDR-transfected cells. We observed genetic differences in the regulation of PTGDR by ATRA that could contribute to the phenotypic differences observed in allergic patients. Our findings showed that RAR modulation by PTGDR might have an impact on Th2 responses, suggesting that RAR could be a potential therapeutic target in allergic inflammation. [ABSTRACT FROM AUTHOR]