Therapeutic benefits of central LH receptor agonism in the APP/PS1 AD model involve trophic and immune regulation and are reproductive status dependent

The mechanisms that underly reproductive hormone effects on cognition, neuronal plasticity, and AD risk, particularly in relation to gonadotropin LH receptor (LHCGR) signaling, remain poorly understood. To address this gap in knowledge and clarify the impact of circulating steroid hormones on the therapeutic effects of CNS LHCGR activation, we delivered the LHCGR agonist human chorionic gonadotropin (hCG) intracerebroventricularly (ICV) and evaluated functional, structural, plasticity-related signaling cascades, Aβ pathology, and transcriptome differences in reproductively intact and ovariectomized (OVX) APP/PS1 AD female mice. Here we demonstrate that CNS hCG delivery restored function to wild-type levels only in OVX APP/PS1 mice. Spine density was increased in all hCG treated groups independently of reproductive status. Notably, increases in BDNF signaling and cognition, were selectively upregulated only in the OVX hCG-treated group. RNA sequencing analyses identified a significant increase in peripheral myeloid and pro-inflammatory genes within the hippocampi of the OVX group that were completely reversed by hCG treatment, identifying a potential mechanism underlying the selective therapeutic effect of LHCGR activation. Interestingly, in intact mice, hCG administration mimicked the effects of gonadectomy. Together, our findings indicate that CNS LHCGR agonism in the post-menopausal context is beneficial through trophic and immune mechanisms. Our findings also underscore the presence of a steroid-LHCGR mechanistic interaction that is unexplored yet potentially meaningful to fully understand "post-menopausal" brain function and CNS hormone treatment response.

The LH receptor regulates hippocampal spatial memory and restores dendritic spine density in ovariectomized APP/PS1 AD mice

Activation of the luteinizing hormone receptor (LHCGR) rescues spatial memory function and spine density losses associated with gonadectomy and high circulating gonadotropin levels in females. However, whether this extends to the AD brain or the mechanisms that underlie these benefits remain unknown. To address this question, we delivered the LHCGR agonist human chorionic gonadotropin (hCG) intracerebroventricularly (ICV), under reproductively intact and ovariectomized conditions to mimic the post-menopausal state in the APP/PS1mouse brain. Cognitive function was tested using the Morris water maze task, and hippocampal dendritic spine density, Aβ pathology, and signaling changes associated with these endpoints were determined to address mechanisms. Here we show that central LHCGR activation restored function in ovariectomized APP/PS1 female mice to wild-type levels without altering Aβ pathology. LHCGR activation increased hippocampal dendritic spine density regardless of reproductive status, and this was mediated by BDNF-dependent and independent signaling. We also show that ovariectomy in the APP/PS1 brain elicits an increase in peripherally derived pro-inflammatory genes which are inhibited by LHCGR activation. This may mediate reproductive status specific effects of LHCGR agonism on cognitive function and BDNF expression. Together, this work highlights the relevance of the LHCGR on cognition and its therapeutic potential in the "menopausal" AD brain.

Alterations in RNA editing in skeletal muscle following exercise training in individuals with Parkinson's disease

Parkinson's Disease (PD) is the second most common neurodegenerative disease behind Alzheimer's Disease, currently affecting more than 10 million people worldwide and 1.5 times more males than females. The progression of PD results in the loss of function due to neurodegeneration and neuroinflammation. The etiology of PD is multifactorial, including both genetic and environmental origins. Here we explored changes in RNA editing, specifically editing through the actions of the Adenosine Deaminases Acting on RNA (ADARs), in the progression of PD. Analysis of ADAR editing of skeletal muscle transcriptomes from PD patients and controls, including those that engaged in a rehabilitative exercise training program revealed significant differences in ADAR editing patterns based on age, disease status, and following rehabilitative exercise. Further, deleterious editing events in protein coding regions were identified in multiple genes with known associations to PD pathogenesis. Our findings of differential ADAR editing complement findings of changes in transcriptional networks identified by a recent study and offer insights into dynamic ADAR editing changes associated with PD pathogenesis.

Changes in ADAR RNA editing patterns in CMV and ZIKV congenital infections

BACKGROUND

RNA editing is a process that increases transcriptome diversity, often through Adenosine Deaminases Acting on RNA (ADARs) that catalyze the deamination of adenosine to inosine. ADAR editing plays an important role in regulating brain function and immune activation, and is dynamically regulated during brain development. Additionally, the ADAR1 p150 isoform is induced by interferons in viral infection and plays a role in antiviral immune response. However, the question of how virus-induced ADAR expression affects host transcriptome editing remains largely unanswered. This question is particularly relevant in the context of congenital infections, given the dynamic regulation of ADAR editing during brain development, the importance of this editing for brain function, and subsequent neurological symptoms of such infections, including microcephaly, sensory issues, and other neurodevelopmental abnormalities. Here, we begin to address this question, examining ADAR expression in publicly available datasets of congenital infections of human cytomegalovirus (HCMV) microarray expression data, as well as mouse cytomegalovirus (MCMV) and mouse/ human induced pluripotent neuroprogenitor stem cell (hiNPC) Zika virus (ZIKV) RNA-seq data.

RESULTS

We found that in all three datasets, ADAR1 was overexpressed in infected samples compared to uninfected samples. In the RNA-seq datasets, editing rates were also analyzed. In all mouse infections cases, the number of editing sites was significantly increased in infected samples, albeit this was not the case for hiNPC ZIKV samples. Mouse ZIKV samples showed altered editing of well-established protein-recoding sites such as Gria3, Grik5, and Nova1, as well as editing sites that may impact miRNA binding.

CONCLUSIONS

Our findings provide evidence for changes in ADAR expression and subsequent dysregulation of ADAR editing of host transcriptomes in congenital infections. These changes in editing patterns of key neural genes have potential significance in the development of neurological symptoms, thus contributing to neurodevelopmental abnormalities. Further experiments should be performed to explore the full range of editing changes that occur in different congenital infections, and to confirm the specific functional consequences of these editing changes.

Comparative transcriptome analysis of Peromyscus leucopus and C3H mice infected with the Lyme disease pathogen

Lyme disease (LD), the most prevalent tick-borne disease of humans in the Northern Hemisphere, is caused by the spirochetal bacterium of Borreliella burgdorferi (Bb) sensu lato complex. In nature, Bb spirochetes are continuously transmitted between Ixodes ticks and mammalian or avian reservoir hosts. Peromyscus leucopus mice are considered the primary mammalian reservoir of Bb in the United States. Earlier studies demonstrated that experimentally infected P. leucopus mice do not develop disease. In contrast, C3H mice, a widely used laboratory strain of Mus musculus in the LD field, develop severe Lyme arthritis. To date, the exact tolerance mechanism of P. leucopus mice to Bb-induced infection remains unknown. To address this knowledge gap, the present study has compared spleen transcriptomes of P. leucopus and C3H/HeJ mice infected with Bb strain 297 with those of their respective uninfected controls. Overall, the data showed that the spleen transcriptome of Bb-infected P. leucopus mice was much more quiescent compared to that of the infected C3H mice. To date, the current investigation is one of the few that have examined the transcriptome response of natural reservoir hosts to Borreliella infection. Although the experimental design of this study significantly differed from those of two previous investigations, the collective results of the current and published studies have consistently demonstrated very limited transcriptomic responses of different reservoir hosts to the persistent infection of LD pathogens.

Importance

The bacterium Borreliella burgdorferi (Bb) causes Lyme disease, which is one of the emerging and highly debilitating human diseases in countries of the Northern Hemisphere. In nature, Bb spirochetes are maintained between hard ticks of Ixodes spp. and mammals or birds. In the United States, the white-footed mouse, Peromyscus leucopus, is one of the main Bb reservoirs. In contrast to humans and laboratory mice (e.g., C3H mice), white-footed mice rarely develop clinical signs (disease) despite being (persistently) infected with Bb. How the white-footed mouse tolerates Bb infection is the question that the present study has attempted to address. Comparisons of genetic responses between Bb-infected and uninfected mice demonstrated that, during a long-term Bb infection, C3H mice reacted much stronger, whereas P. leucopus mice were relatively unresponsive.

Genetic characterization of Staphylococcus aureus isolated from Norway rats in Boston, Massachusetts

BACKGROUND

Despite the importance of domesticated animals in the generation and transmission of antibiotic-resistant Staphylococcus aureus, the role of wild animals, specifically rodents, in the ecology of S. aureus remains unclear. We recovered and genotyped S. aureus isolates from wild Norway rats (Rattus norvegicus) in Boston, Massachusetts to examine genetic relationships between common human and animal S. aureus isolates in a large US metropolitan area.

METHODS

We collected and necropsied 63 rats from June 2016 to June 2017. Nasal, foot pad, fur, and fecal swabs were collected. Staphylococcus aureus was isolated using culture-based methods and polymerase chain reaction confirmation. S. aureus isolates were spa typed, tested for antibiotic susceptibility, and whole genome sequenced. Assembled sequences were uploaded to the Comprehensive Antibiotic Resistance Database to identify antibiotic resistance elements. A phylogenetic tree was constructed using the neighbor-joining method with the maximum composite likelihood distance in MEGA7.

RESULTS

We recovered 164 Gram-positive bacterial isolates from Norway rats. Nineteen isolates from eight individual rats were confirmed as S. aureus (prevalence: 12.9% (8/63)). All S. aureus isolates were methicillin-susceptible S. aureus (MSSA), pvl-negative, and resistant to penicillin. Two isolates displayed resistance to erythromycin. Four different S. aureus spa types were detected (t933, t10751, t18202, and t189). Thirteen unique antibiotic resistance elements were identified, and all isolates shared genes mepR, mgrA, arlR, and S. aureus norA. Phylogenetic analysis if the 19 S. aureus isolates revealed they were genetically similar to four clades of S. aureus with similar resistance gene profiles isolated from both human- and animal-derived S. aureus, as well as formed a distinct phylogenetic cluster composed only of rat isolates.

CONCLUSIONS

Wild rodents may serve as a reservoir or vector of antibiotic resistance genes in the urban environment with relevance for human and animal health.

Neuroprotective Mechanisms of Amylin Receptor Activation, Not Antagonism, in the APP/PS1 Mouse Model of Alzheimer's Disease

BACKGROUND

Amylin, a pancreatic amyloid peptide involved in energy homeostasis, is increasingly studied in the context of Alzheimer's disease (AD) etiology. To date, conflicting pathogenic and neuroprotective roles for this peptide and its analogs for AD pathogenesis have been described.

OBJECTIVE

Whether the benefits of amylin are associated with peripheral improvement of metabolic tone/function or directly through the activation of central amylin receptors is also unknown and downstream signaling mechanisms of amylin receptors are major objectives of this study.

METHODS

To address these questions more directly we delivered the amylin analog pramlintide systemically (IP), at previously identified therapeutic doses, while centrally (ICV) inhibiting the receptor using an amylin receptor antagonist (AC187), at doses known to impact CNS function.

RESULTS

Here we show that pramlintide improved cognitive function independently of CNS receptor activation and provide transcriptomic data that highlights potential mechanisms. Furthermore, we show than inhibition of the amylin receptor increased amyloid-beta pathology in female APP/PS1 mice, an effect than was mitigated by peripheral delivery of pramlintide. Through transcriptomic analysis of pramlintide therapy in AD-modeled mice we found sexual dimorphic modulation of neuroprotective mechanisms: oxidative stress protection in females and membrane stability and reduced neuronal excitability markers in males.

CONCLUSION

These data suggest an uncoupling of functional and pathology-related events and highlighting a more complex receptor system and pharmacological relationship that must be carefully studied to clarify the role of amylin in CNS function and AD.

Metabolic labeling unveils alterations in the turnover of HDL-associated proteins during diabetes progression in mice

Several aspects of diabetes pathophysiology and complications result from hyperglycemia-induced alterations in the structure and function of plasma proteins. Furthermore, insulin has a significant influence on protein metabolism by affecting both the synthesis and degradation of proteins in various tissues. To understand the role of progressive hyperglycemia on plasma proteins, in this study, we measured the turnover rates of high-density lipoprotein (HDL)-associated proteins in control (chow diet), prediabetic [a high-fat diet (HFD) for 8 wk] or diabetic [HFD for 8 wk with low-dose streptozotocin (HFD + STZ) in weeks 5-8 of HFD] C57BL/6J mice using heavy water (2H2O)-based metabolic labeling approach. Compared with control mice, HFD and HFD + STZ mice showed elevations of fasting plasma glucose levels in the prediabetic and diabetic range, respectively. Furthermore, the HFD and HFD + STZ mice showed increased hepatic triglyceride (TG) levels, total plasma cholesterol, and plasma TGs. The kinetics of 40 proteins were quantified using the proteome dynamics method, which revealed an increase in the fractional synthesis rate (FSR) of HDL-associated proteins in the prediabetic mice compared with control mice, and a decrease in FSR in the diabetic mice. The pathway analysis revealed that proteins with altered turnover rates were involved in acute-phase response, lipid metabolism, and coagulation. In conclusion, prediabetes and diabetes have distinct effects on the turnover rates of HDL proteins. These findings suggest that an early dysregulation of the HDL proteome dynamics can provide mechanistic insights into the changes in protein levels in these conditions.NEW & NOTEWORTHY This study is the first to examine the role of gradual hyperglycemia during diabetes disease progression on HDL-associated protein dynamics in the prediabetes and diabetic mice. Our results show that the fractional synthesis rate of HDL-associated proteins increased in the prediabetic mice whereas it decreased in the diabetic mice compared with control mice. These kinetic changes can help to elucidate the mechanism of altered protein levels and HDL dysfunction during diabetes disease progression.

Amylin Pharmacology in Alzheimer's Disease Pathogenesis and Treatment

The metabolic peptide hormone amylin, in concert with other metabolic peptides like insulin and leptin, has an important role in metabolic homeostasis and has been intimately linked to Alzheimer's disease (AD). Interestingly, this pancreatic amyloid peptide is known to self-aggregate much like amyloid-beta and has been reported to be a source of pathogenesis in both Type II diabetes mellitus (T2DM) and Alzheimer's disease. The traditional "gain of toxic function" properties assigned to amyloid proteins are, however, contrasted by several reports highlighting neuroprotective effects of amylin and a recombinant analog, pramlintide, in the context of these two diseases. This suggests that pharmacological therapies aimed at modulating the amylin receptor may be therapeutically beneficial for AD development, as they already are for T2DMM. However, the nature of amylin receptor signaling is highly complex and not well studied in the context of CNS function. Therefore, to begin to address this pharmacological paradox in amylin research, the goal of this review is to summarize the current research on amylin signaling and CNS functions and critically address the paradoxical nature of this hormone's signaling in the context of AD pathogenesis.

Proteome dynasmics and bioinformatics reveal major alterations in the turnover rate of functionally related cardiac and plasma proteins in a dog model of congestive heart failure

Protein pool turnover is a critically important cellular homeostatic component, yet it has been little explored in the context of heart failure (HF) pathophysiology. We employed in vivo ²H labeling/ proteome dynamics for non-biased discovery of turnover alterations involving functionally linked cardiac and plasma proteins in canine tachypacing-induced HF, an established preclinical model of dilated cardiomyopathy. Compared to control, dogs with congestive HF displayed bidirectional turnover changes of 28 cardiac proteins, i.e. reduced half-life of several key enzymes involved in glycolysis, homocysteine metabolism and glycogenesis, and increased half-life of proteins involved in proteolysis. Changes in plasma proteins were more modest: only 5 proteins, involved in various functions including proteolysis inhibition, hemoglobin, calcium and ferric-iron binding, displayed increased or decreased turnover rates. In other dogs undergoing cardiac tachypacing, we infused for 2 weeks the myokine Follistatin-like protein 1 (FSTL1), known for its ameliorative effects on HF-induced alterations. Proteome dynamics proved very sensitive in detecting the partial or complete prevention, by FSTL1, of cardiac and plasma protein turnover alterations. In conclusion, our study unveiled, for the first time in a large mammal, numerous HF-related alterations that may serve as the basis for future mechanistic research and/or as conceptually new molecular markers.

ADAR Editing in Viruses: An Evolutionary Force to Reckon with

Adenosine Deaminases that Act on RNA (ADARs) are RNA editing enzymes that play a dynamic and nuanced role in regulating transcriptome and proteome diversity. This editing can be highly selective, affecting a specific site within a transcript, or nonselective, resulting in hyperediting. ADAR editing is important for regulating neural functions and autoimmunity, and has a key role in the innate immune response to viral infections, where editing can have a range of pro- or antiviral effects and can contribute to viral evolution. Here we examine the role of ADAR editing across a broad range of viral groups. We propose that the effect of ADAR editing on viral replication, whether pro- or antiviral, is better viewed as an axis rather than a binary, and that the specific position of a given virus on this axis is highly dependent on virus- and host-specific factors, and can change over the course of infection. However, more research needs to be devoted to understanding these dynamic factors and how they affect virus-ADAR interactions and viral evolution. Another area that warrants significant attention is the effect of virus-ADAR interactions on host-ADAR interactions, particularly in light of the crucial role of ADAR in regulating neural functions. Answering these questions will be essential to developing our understanding of the relationship between ADAR editing and viral infection. In turn, this will further our understanding of the effects of viruses such as SARS-CoV-2, as well as many others, and thereby influence our approach to treating these deadly diseases.

Transcriptional response of key metabolic and stress response genes of a nuculanid bivalve, Lembulus bicuspidatus from an oxygen minimum zone exposed to hypoxia-reoxygenation

Benthic animals inhabiting the edges of marine oxygen minimum zones (OMZ) are exposed to unpredictable large fluctuations of oxygen levels. Sessile organisms including bivalves must depend on physiological adaptations to withstand these conditions. However, as habitats are rather inaccessible, physiological adaptations of the OMZ margin inhabitants to oxygen fluctuations are not well understood. We therefore investigated the transcriptional responses of selected key genes involved in energy metabolism and stress protection in a dominant benthic species of the northern edge of the Namibian OMZ, the nuculanid clam Lembulus bicuspidatus,. We exposed clams to normoxia (~5.8 ml O₂ l^-1), severe hypoxia (36 h at ~0.01 ml O₂ l^-1) and post-hypoxic recovery (24 h of normoxia following 36 h of severe hypoxia). Using newly identified gene sequences, we determined the transcriptional responses to hypoxia and reoxygenation of the mitochondrial aerobic energy metabolism (pyruvate dehydrogenase E1 complex, cytochrome c oxidase, citrate synthase, and adenine nucleotide translocator), anaerobic glycolysis (hexokinase (HK), phosphoenolpyruvate carboxykinase (PEPCK), phosphofructokinase, and aldolase), mitochondrial antioxidants (glutaredoxin, peroxiredoxin, and uncoupling protein UCP2) and stress protection mechanisms (a molecular chaperone HSP70 and a mitochondrial quality control protein MIEAP) in the gills and the labial palps of L. bicuspidatus. Exposure to severe hypoxia transcriptionally stimulated anaerobic glycolysis (including HK and PEPCK), antioxidant protection (UCP2), and quality control mechanisms (HSP70 and MIEAP) in the gills of L. bicuspidatus. Unlike UCP2, mRNA levels of the thiol-dependent mitochondrial antioxidants were not affected by hypoxia-reoxygenation stress. Transcript levels of marker genes for aerobic energy metabolism were not responsive to oxygen fluctuations in L. bicuspidatus. Our findings highlight the probable importance of anaerobic succinate production (via PEPCK) and mitochondrial and proteome quality control mechanisms in responses to oxygen fluctuations of the OMZ bivalve L.bicuspidatus. The reaction of L.bicuspidatus to oxygen fluctuations implies parallels to that of other hypoxia-tolerant bivalves, such as intertidal species.

Unintended consequences of innate immune activation by coronaviruses

Background

Methods

Results

Conclusions

Remote Conversations To Enhance Class Experience in the Time of COVID: Co-Teaching with a Wingman Model

Remote delivery poses numerous challenges, particularly for discussion-based classes. Here, we describe a class management model that uses a conversational approach between instructional coleaders to introduce topics and facilitate student discussions. The conversational approach enabled us to model various aspects of the scientific process, such as asking questions and generating hypotheses, aimed at an emergent subject, such as COVID-19. Students' feedback was positive, noting how the class enabled them to connect their existing knowledge with the daily influx of new insights about SARS-CoV-2. While our class focuses on COVID-19-related topics, this model of a conversational style class can be easily implemented for other course topics.

Automated Isoform Diversity Detector (AIDD): a pipeline for investigating transcriptome diversity of RNA-seq data

Background

As the number of RNA-seq datasets that become available to explore transcriptome diversity increases, so does the need for easy-to-use comprehensive computational workflows. Many available tools facilitate analyses of one of the two major mechanisms of transcriptome diversity, namely, differential expression of isoforms due to alternative splicing, while the second major mechanism—RNA editing due to post-transcriptional changes of individual nucleotides—remains under-appreciated. Both these mechanisms play an essential role in physiological and diseases processes, including cancer and neurological disorders. However, elucidation of RNA editing events at transcriptome-wide level requires increasingly complex computational tools, in turn resulting in a steep entrance barrier for labs who are interested in high-throughput variant calling applications on a large scale but lack the manpower and/or computational expertise.

Results

Here we present an easy-to-use, fully automated, computational pipeline (Automated Isoform Diversity Detector, AIDD) that contains open source tools for various tasks needed to map transcriptome diversity, including RNA editing events. To facilitate reproducibility and avoid system dependencies, the pipeline is contained within a pre-configured VirtualBox environment. The analytical tasks and format conversions are accomplished via a set of automated scripts that enable the user to go from a set of raw data, such as fastq files, to publication-ready results and figures in one step. A publicly available dataset of Zika virus-infected neural progenitor cells is used to illustrate AIDD’s capabilities.

Conclusions

AIDD pipeline offers a user-friendly interface for comprehensive and reproducible RNA-seq analyses. Among unique features of AIDD are its ability to infer RNA editing patterns, including ADAR editing, and inclusion of Guttman scale patterns for time series analysis of such editing landscapes. AIDD-based results show importance of diversity of ADAR isoforms, key RNA editing enzymes linked with the innate immune system and viral infections. These findings offer insights into the potential role of ADAR editing dysregulation in the disease mechanisms, including those of congenital Zika syndrome. Because of its automated all-inclusive features, AIDD pipeline enables even a novice user to easily explore common mechanisms of transcriptome diversity, including RNA editing landscapes.

Diversifying selection detected in only a minority of xenobiotic-metabolizing CYP1-3 genes among primate species

1. Primates exhibit a high degree of among-species dietary diversity, which likely exposes them to varying levels of xenobiotic compounds. Here, we examined the evolution of primate CYP1-3 gene families, and we classified the 15 CYP1-3 gene subfamilies as either xenobiotic-metabolizing (XM) or endogenous-metabolizing (EM) based on sources in the P450 literature. 2. We predicted that XM P450s would show (1) greater variability in gene-copy number and (2) more evidence of diversifying selection and, especially on codons that encode the substrate-recognition sites (SRSs) for the final enzymes. 3. Counter to our first prediction, EM and XM P450s showed similar levels of variation in gene-copy number. We did find, however, that four XM P450 subfamilies (CYP2C, CYP2D, CYP2E, and CYP3A) showed evidence of diversifying selection while no EM subfamilies demonstrated any consistent signal of diversifying selection. Of these four, CYP2C, CYP2D, and CYP3A showed significant links between SRSs and diversifying selection. 4. These results reveal an amount of evolutionary dynamism that would not be expected when viewing P450 subfamilies along a simple binary EM/XM spectrum. We recommend that comparative studies of cytochrome P450 evolution should focus on the CYP2C, CYP2D, CYP2E, and CYP3A subfamilies.

Machine Learning Maps Research Needs in COVID-19 Literature

As of August 2020, thousands of COVID-19 (coronavirus disease 2019) publications have been produced. Manual assessment of their scope is an overwhelming task, and shortcuts through metadata analysis (e.g., keywords) assume that studies are properly tagged. However, machine learning approaches can rapidly survey the actual text of publication abstracts to identify research overlap between COVID-19 and other coronaviruses, research hotspots, and areas warranting exploration. We propose a fast, scalable, and reusable framework to parse novel disease literature. When applied to the COVID-19 Open Research Dataset, dimensionality reduction suggests that COVID-19 studies to date are primarily clinical, modeling, or field based, in contrast to the vast quantity of laboratory-driven research for other (non-COVID-19) coronavirus diseases. Furthermore, topic modeling indicates that COVID-19 publications have focused on public health, outbreak reporting, clinical care, and testing for coronaviruses, as opposed to the more limited number focused on basic microbiology, including pathogenesis and transmission.

•

AI/machine learning techniques can analyze coronavirus research at massive scale

•

COVID-19 research has so far focused on non-lab-based (e.g., observational) research

•

COVID-19 lab-based/basic microbiological research is less prevalent than expected

The impact of the COVID-19 pandemic has led scientists to produce a vast quantity of research aimed at understanding, monitoring, and containing the disease; however, it remains unclear whether the research that has been produced to date sufficiently addresses existing knowledge gaps. We use artificial intelligence (AI)/machine learning techniques to analyze this massive amount of information at scale. We find key discrepancies between literature about COVID-19 and what we would expect based on research on other coronaviruses. These discrepancies—namely, the lack of basic microbiological research, which is often expensive and time-consuming—may negatively impact efforts to mitigate the pandemic and raise questions regarding the research community's ability to quickly respond to future crises. Continually measuring what is being produced, both now and in the future, is key to making better resource allocation and goal prioritization decisions as a society moving forward.

Lipid-protein interactions for ECA1 an N-ANTH domain protein involved in stress signaling in plants

Epsin-like Clathrin Adaptor 1 (ECA1/ PICALM1A) is an A/ENTH domain protein that acts as an adaptor protein in clathrin-mediated endocytosis. ECA1 is recruited to the membrane during salt stress signaling in plants in a phosphatidic acid (PA)-dependent manner. PA is a lipid second messenger that rapidly and transiently increases in concentration under stress stimuli. Upon an increase in PA concentration another lipid, diacylglycerol pyrophosphate (DGPP), starts to accumulate. The accumulation of DGPP is suggested to be a cue for attenuating PA signaling during stress in plants. We showed in vitro that ECA1-PA binding is modulated as a function of membrane curvature stress and charge. In this work, we investigate ECA1 binding to DGPP in comparison with PA. We show that ECA1 has more affinity for the less charged PA, and this binding is pH dependent. Additionally, plant PA binding proteins SnRK2.10, TGD2C, and PDK1-PH2 were investigated for their interaction with DGPP, since no known DGPP binding proteins are available in the literature to date. Our results shed further light on DGPP and its interactions with membrane proteins which brings us closer toward understanding the complexity of protein interactions with anionic lipids, especially the enigmatic anionic lipid DGPP.

Interferons and viruses induce a novel primate-specific isoform dACE2 and not the SARS-CoV-2 receptor ACE2

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which causes COVID-19, utilizes angiotensin-converting enzyme 2 (ACE2) for entry into target cells. ACE2 has been proposed as an interferon-stimulated gene (ISG). Thus, interferon-induced variability in ACE2 expression levels could be important for susceptibility to COVID-19 or its outcomes. Here, we report the discovery of a novel, primate-specific isoform of ACE2, which we designate as deltaACE2 (dACE2). We demonstrate that dACE2, but not ACE2, is an ISG. In vitro, dACE2, which lacks 356 N-terminal amino acids, was non-functional in binding the SARS-CoV-2 spike protein and as a carboxypeptidase. Our results reconcile current knowledge on ACE2 expression and suggest that the ISG-type induction of dACE2 in IFN-high conditions created by treatments, inflammatory tumor microenvironment, or viral co-infections is unlikely to affect the cellular entry of SARS-CoV-2 and promote infection.

COVID-19 and emerging viral infections: The case for interferon lambda

With the first reports on coronavirus disease 2019 (COVID-19), which is caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the scientific community working in the field of type III IFNs (IFN-λ) realized that this class of IFNs could play an important role in this and other emerging viral infections. In this Viewpoint, we present our opinion on the benefits and potential limitations of using IFN-λ to prevent, limit, and treat these dangerous viral infections.

Interactive effects of salinity variation and exposure to ZnO nanoparticles on the innate immune system of a sentinel marine bivalve, Mytilus edulis

ZnO nanoparticles (nZnO) are released into the coastal environment from multiple sources, yet their toxicity to marine organisms is not well understood. We investigated the interactive effects of salinity (normal 15, low 5, and fluctuating 5-15) and nZnO (100 μg l-1) on innate immunity of the blue mussels Mytilus edulis from a brackish area of the Baltic Sea. Exposure to ionic Zn (100 μg l-1) was used to test whether the toxic effects of nZnO can be attributed to the potential release of Zn2+. Functional parameters and the expression of key immune-related genes were investigated in the mussels exposed to nZnO or ionic Zn under different salinity regimes for 21 days. nZnO exposures elevated hemocyte mortality, suppressed adhesion, stimulated phagocytosis, and led to an apparent increase in lysosomal volume. At salinity 15, nZnO suppressed the mRNA expression of the Toll-like receptors TLRb and c, C-lectin, and the complement system component C3q indicating impaired ability for pathogen recognition. In contrast, the mRNA levels of an antimicrobial peptide defensin increased during nZnO exposure at salinity 15. At fluctuating salinity (5-15), nZnO exposure increased expression of multiple immune-related genes in hemocytes including the complement system components C1 and C3q, and the Toll-like receptors TLRa, b and c. Low salinity (5) had strong immunosuppressive effects on the functional and molecular immune traits of M. edulis that overshadowed the effects of nZnO. The salinity-dependent modulation of immune response to nZnO cannot be attributed to the differences in the aggregation or solubility of nZnO, and likely reflects the interaction of the toxic effects of nanoparticles and physiological effects of the osmotic stress. These findings have implications for the environmental risk assessment of nanomaterials and the development of the context-specific biomarker baselines for coastal pollution monitoring.

Effects of intermittent hypoxia on cell survival and inflammatory responses in the intertidal marine bivalves Mytilus edulis and Crassostrea gigas

Hypoxia is a major stressor in estuarine and coastal habitats, leading to adverse effects in aquatic organisms. Estuarine bivalves such as blue mussels (Mytilus edulis) and Pacific oysters (Crassostrea gigas) can survive periodic oxygen deficiency but the molecular mechanisms that underlie cellular injury during hypoxia-reoxygenation are not well understood. We examined the molecular markers of autophagy, apoptosis and inflammation during short-term (1 day) and long-term (6 days) hypoxia and post-hypoxic recovery (1 h) in mussels and oysters by measuring the lysosomal membrane stability, activity of a key autophagic enzyme (cathepsin D) and mRNA expression of the genes involved in the cellular survival and inflammation, including caspase 2, 3 and 8, Bcl-2, BAX, TGF-β-activated kinase 1 (TAK1), nuclear factor kappa B1 (NF-κB) and NF-κB activating kinases IKKα and TBK1. Crassostrea gigas exhibited higher hypoxia tolerance, as well as blunted or delayed inflammatory and apoptotic response to hypoxia and reoxygenation as shown by the later onset and/or the lack of transcriptional activation of caspases, BAX and the inflammatory effector NF-κB, compared with M. edulis Long-term hypoxia resulted in upregulation of Bcl-2 in the oysters and mussels, implying activation of anti-apoptotic mechanisms. Our findings indicate the potential importance of the cell survival pathways in hypoxia tolerance of marine bivalves, and demonstrate the utility of the molecular markers of apoptosis and autophagy for the assessment of sublethal hypoxic stress in bivalve populations.

Explaining Pathogenicity of Congenital Zika and Guillain-Barré Syndromes: Does Dysregulation of RNA Editing Play a Role?

Previous studies of Zika virus (ZIKV) pathogenesis have focused primarily on virus-driven pathology and neurotoxicity, as well as host-related changes in cell proliferation, autophagy, immunity, and uterine function. It is now hypothesized that ZIKV pathogenesis arises instead as an (unintended) consequence of host innate immunity, specifically, as the side effect of an otherwise well-functioning machine. The hypothesis presented here suggests a new way of thinking about the role of host immune mechanisms in disease pathogenesis, focusing on dysregulation of post-transcriptional RNA editing as a candidate driver of a broad range of observed neurodevelopmental defects and neurodegenerative clinical symptoms in both infants and adults linked with ZIKV infections. The authors collect and synthesize existing evidence of ZIKV-mediated changes in the expression of adenosine deaminases acting on RNA (ADARs), known links between abnormal RNA editing and pathogenesis, as well as ideas for future research directions, including potential treatment strategies.

Stable G-quadruplex enabling sequences are selected against by the context-dependent codon bias

The distributions of secondary structural elements appear to differ between coding regions (CDS) of mRNAs compared to the untranslated regions (UTRs), presumably as a mechanism to fine-tune gene expression, including efficiency of translation. However, a systematic and comprehensive analysis of secondary structure avoidance because of potential bias in codon usage is difficult as some of the common secondary structures, such as, hairpins can be formed by numerous sequence combinations. Using G-quadruplex (GQ) as the model secondary structure we studied the impact of codon bias on GQs within the CDS. Because GQs can be predicted using specific consensus sequence motifs, they provide an excellent platform for investigation of the selectivity of such putative structures at the codon level. Using a bioinformatics approach, we calculated the frequencies of putative GQs within the CDS of a variety of species. Our results suggest that the most stable GQs appear to be significantly underrepresented within the CDS, through the use of specific synonymous codon combinations. Furthermore, we identified many peptide sequence motifs in which silent mutations can potentially alter translation via stable GQ formation. This work not only provides a comprehensive analysis on how stable secondary structures appear to be avoided within the CDS of mRNA, but also broadens the current understanding of synonymous codon usage as they relate to the structure-function relationship of RNA.

Potential trade-offs between biomineralization and immunity revealed by shell properties and gene expression profiles of two closely related Crassostrea species

Species of the Ostreidae family are key ecosystem engineers and many of them - including Crassostrea gigas and Crassostreavirginica - are commercially important aquaculture species. Despite similarities in their morphology and ecology, these two species differ in their ability to defend against pathogens, potentially reflecting species-specific differential specialization of hemocytes on immune defense versus biomineralization. To test this hypothesis, we investigated the expression levels of immune- and biomineralization-related genes as well as mineralogical and mechanical properties of the shells and the calcium sequestration ability of the hemocytes of C. gigas and C. virginica The expression of biomineralization-related genes was higher in C. virginica than in C. gigas in multiple tissues including the mantle edge and hemocytes, while the expression of immune genes was higher in the hemocytes of C. gigas Hemocytes of C. virginica contained more calcium (stored intracellularly as calcium carbonate mineral) compared with those of C. gigas Analysis of the adult shells showed that the crystallinity of calcite was higher and the laths of the foliated layer of the shell were thicker in C. virginica than in C. gigas Mechanically, the shells of C. virginica were stiffer, harder and stronger than those of C. gigas Taken together, our results show that the species-specific differences in physiology (such as disease resistance and exoskeleton properties) are reflected at the cellular and molecular levels in the differential specialization of hemocytes on potentially competing functions (immunity and biomineralization) as well as different expression profiles of other tissues involved in biomineralization (such as the mantle edge).

Matchmaking in Bioinformatics

Ever return from a meeting feeling elated by all those exciting talks, yet unsure how all those presented glamorous and/or exciting tools can be useful in your research?  Or do you have a great piece of software you want to share, yet only a handful of people visited your poster? We have all been there, and that is why we organized the Matchmaking for Computational and Experimental Biologists Session at the latest ISCB/GLBIO’2017 meeting in Chicago (May 15-17, 2017). The session exemplifies a novel approach, mimicking “matchmaking”, to encouraging communication, making connections and fostering collaborations between computational and non-computational biologists. More specifically, the session facilitates face-to-face communication between researchers with similar or differing research interests, which we feel are critical for promoting productive discussions and collaborations.  To accomplish this, three short scheduled talks were delivered, focusing on RNA-seq, integration of clinical and genomic data, and chromatin accessibility analyses.  Next, small-table developer-led discussions, modeled after speed-dating, enabled each developer (including the speakers) to introduce a specific tool and to engage potential users or other developers around the table.  Notably, we asked the audience whether any other tool developers would want to showcase their tool and we thus added four developers as moderators of these small-table discussions.  Given the positive feedback from the tool developers, we feel that this type of session is an effective approach for promoting valuable scientific discussion, and is particularly helpful in the context of conferences where the number of participants and activities could hamper such interactions.

Retained duplications and deletions of CYP2C genes among primates

The human genome encodes about 60 functional enzymes of the cytochrome P450 superfamily, including four functional enzymes of the cytochrome P450 2C (CYP2C) subfamily. These enzymes have been shown to metabolize drugs and xenobiotic toxins, such as those in the diet, and are therefore of great importance for biomedical research and applications. While the pharmacology of P450 enzymes has been studied extensively, our understanding of molecular evolution of this gene family is incomplete, in part because a great variation exists in the number of CYP2C homologs across genomes. In humans, the enzymes encoded by these genes are responsible for the metabolism of more than 20% of clinical drugs, but this is not the naturalistic function of these enzymes, which is the metabolism of xenobiotics such as plant secondary metabolites. In this paper, we sought to correlate evolutionary relationships among primate CYP2C genes with known dietary profiles from these species, testing the hypothesis that these genes have evolved under the pressure of dietary toxins. Aside from a small number of deeply divergent genes, primate CYP2C paralogs form three separate clades: CYP2C18, CYP2C9/CYP2C19, and CYP2C8/CYP2C20. Our results showed that the CYP2C18 gene has been separately lost in Nomascus leucogenys and the Panini genomes, and there is no evidence that this gene has been under any positive selection among primates. While CYP2C20 has been retained in cercopithecoids, orthologous loci were separately lost in platyrrhines and hominoids. Notably, nine codons exhibited signature of positive selection. Finally, the CYP2C19 locus was duplicated in basal catarrhines, resulting in the birth of CYP2C9; but the ancestral locus was only retained in hominoid taxa. Overall, our findings support the hypothesis that primate CYP2C genes have evolved in response to selective pressures provided by dietary toxins, although not all gene clusters have evolved in the same manner. Our results may indicate an evolutionarily deep difference in ecology or physiology among higher-order primate taxa.

Seasonality of birth defects in West Africa: could congenital Zika syndrome be to blame?

The link between Zika virus infection during pregnancy and microcephaly and other neurodevelopmental defects in infants, referred to as congenital Zika syndrome (CZS), was recently discovered. One key question that remains is whether such neurodevelopmental abnormalities are limited to the recently evolved Asiatic ZIKV strains or if they can also be induced by endemic African strains. Thus, we examined birth registries from one particular hospital from a country in West Africa, where ZIKV is endemic. Results showed a seasonal pattern of birth defects that is consistent with potential CZS, which corresponds to a range of presumed maternal infection that encompasses both the peak of the warm, rainy season as well as the months immediately following it, when mosquito activity is likely high. While we refrain from definitively linking ZIKV infection and birth defects in West Africa at this time, in part due to scant data available from the region, we hope that this report will initiate broader surveillance efforts that may help shed light onto mechanisms underlying CZS.

Biomineralization-related specialization of hemocytes and mantle tissues of the Pacific oyster Crassostrea gigas

The molluscan exoskeleton (shell) plays multiple important roles including structural support, protection from predators and stressors, and physiological homeostasis. Shell formation is a tightly regulated biological process that allows molluscs to build their shells even in environments unfavorable for mineral precipitation. Outer mantle edge epithelial cells (OME) and hemocytes were implicated in this process; however, the exact functions of these cell types in biomineralization are not clear. Pacific oysters (Crassostrea gigas) were used to study differences in the expression profiles of selected biomineralization-related genes in hemocytes and mantle cells, and the functional characteristics of hemocytes such as adhesion, motility and phagocytosis. The specialized role of OME in shell formation was supported by high expression levels of the extracellular matrix (ECM) related and cell-cell interaction genes. Density gradient separation of hemocytes revealed distinct phenotypes based on the cell morphology, gene expression patterns, motility and adhesion characteristics. These hemocyte fractions can be categorized into two functional groups, i.e. biomineralization and immune response cells. Gene expression profiles of the putative biomineralizing hemocytes indicate that in addition to their proposed role in mineral transport, hemocytes also contribute to the formation of the ECM, thus challenging the current paradigm of the mantle as the sole source of the ECM for shell formation. Our findings corroborate the specialized roles of hemocytes and the OME in biomineralization and emphasize complexity of the biological controls over shell formation in bivalves.

RNA editing by the host ADAR system affects the molecular evolution of the Zika virus

Zika virus (ZIKV) is a mosquito‐transmitted flavivirus, linked to microcephaly and fetal death in humans. Here, we investigate whether host‐mediated RNA editing of adenosines (ADAR) plays a role in the molecular evolution of ZIKV. Using complete coding sequences for the ZIKV polyprotein, we show that potential ADAR substitutions are underrepresented at the ADAR‐resistant GA dinucleotides of both the positive and negative strands, that these changes are spatially and temporally clustered (as expected of ADAR editing) for three branches of the viral phylogeny, and that ADAR mutagenesis can be linked to its codon usage. Furthermore, resistant GA dinucleotides are enriched on the positive (but not negative) strand, indicating that the former is under stronger purifying selection than the latter. ADAR editing also affects the evolution of the rhabdovirus sigma. Our study now documents that host ADAR editing is a mutation and evolutionary force of positive‐ as well as negative‐strand RNA viruses.

Role of Host-Driven Mutagenesis in Determining Genome Evolution of Sigma Virus (DMelSV; Rhabdoviridae) in Drosophila melanogaster

Sigma virus (DMelSV) is ubiquitous in natural populations of Drosophila melanogaster. Host-mediated, selective RNA editing of adenosines to inosines (ADAR) may contribute to control of viral infection by preventing transcripts from being transported into the cytoplasm or being translated accurately; or by increasing the viral genomic mutation rate. Previous PCR-based studies showed that ADAR mutations occur in DMelSV at low frequency. Here we use SOLiD^TM deep sequencing of flies from a single host population from Athens, GA, USA to comprehensively evaluate patterns of sequence variation in DMelSV with respect to ADAR. GA dinucleotides, which are weak targets of ADAR, are strongly overrepresented in the positive strand of the virus, consistent with selection to generate ADAR resistance on this complement of the transient, double-stranded RNA intermediate in replication and transcription. Potential ADAR sites in a worldwide sample of viruses are more likely to be “resistant” if the sites do not vary among samples. Either variable sites are less constrained and hence are subject to weaker selection than conserved sites, or the variation is driven by ADAR. We also find evidence of mutations segregating within hosts, hereafter referred to as hypervariable sites. Some of these sites were variable only in one or two flies (i.e., rare); others were shared by four or even all five of the flies (i.e., common). Rare and common hypervariable sites were indistinguishable with respect to susceptibility to ADAR; however, polymorphism in rare sites were more likely to be consistent with the action of ADAR than in common ones, again suggesting that ADAR is deleterious to the virus. Thus, in DMelSV, host mutagenesis is constraining viral evolution both within and between hosts.

Coevolutionary Analysis Identifies Protein-Protein Interaction Sites between HIV-1 Reverse Transcriptase and Integrase

The replication of human immunodeficiency virus-1 (HIV-1) requires reverse transcription of the viral RNA genome and integration of newly synthesized pro-viral DNA into the host genome. This is mediated by the viral proteins reverse transcriptase (RT) and integrase (IN). The formation and stabilization of the pre-integration complex (PIC), which is an essential step for reverse transcription, nuclear import, chromatin targeting, and subsequent integration, involves direct and indirect modes of interaction between RT and IN proteins. While epitope-based treatments targeting IN-viral DNA and IN-RT complexes appear to be a promising combination for an anti-HIV treatment, the mechanisms of IN-RT interactions within the PIC are not well understood due to the transient nature of the protein complex and the intrinsic flexibility of its components. Here, we identify potentially interacting regions between the IN and RT proteins within the PIC through the coevolutionary analysis of amino acid sequences of the two proteins. Our results show that specific regions in the two proteins have strong coevolutionary signatures, suggesting that these regions either experience direct and prolonged interactions between them that require high affinity and/or specificity or that the regions are involved in interactions mediated by dynamic conformational changes and, hence, may involve both direct and indirect interactions. Other regions were found to exhibit weak, but positive correlations, implying interactions that are likely transient and/or have low affinity. We identified a series of specific regions of potential interactions between the IN and RT proteins (e.g., specific peptide regions within the C-terminal domain of IN were identified as potentially interacting with the Connection domain of RT). Coevolutionary analysis can serve as an important step in predicting potential interactions, thus informing experimental studies. These studies can be integrated with structural data to gain a better understanding of the mechanisms of HIV protein interactions.

Molecular trajectories leading to the alternative fates of duplicate genes

Gene duplication generates extra gene copies in which mutations can accumulate without risking the function of pre-existing genes. Such mutations modify duplicates and contribute to evolutionary novelties. However, the vast majority of duplicates appear to be short-lived and experience duplicate silencing within a few million years. Little is known about the molecular mechanisms leading to these alternative fates. Here we delineate differing molecular trajectories of a relatively recent duplication event between humans and chimpanzees by investigating molecular properties of a single duplicate: DNA sequences, gene expression and promoter activities. The inverted duplication of the Glutathione S-transferase Theta 2 (GSTT2) gene had occurred at least 7 million years ago in the common ancestor of African great apes and is preserved in chimpanzees (Pan troglodytes), whereas a deletion polymorphism is prevalent in humans. The alternative fates are associated with expression divergence between these species, and reduced expression in humans is regulated by silencing mutations that have been propagated between duplicates by gene conversion. In contrast, selective constraint preserved duplicate divergence in chimpanzees. The difference in evolutionary processes left a unique DNA footprint in which dying duplicates are significantly more similar to each other (99.4%) than preserved ones. Such molecular trajectories could provide insights for the mechanisms underlying duplicate life and death in extant genomes.

Evolutionary diversification of the avian fatty acid-binding proteins

Phylogenetic analysis of avian and other vertebrate fatty acid binding proteins (FABPs) supported the hypothesis that several gene duplications within this family occurred prior to the most recent common ancestor (MRCA) of tetrapods and bony fishes. The chicken genome encodes two liver-expressed FABPs: (1) L-FABP or FABP1; and (2) Lb-FABP. We propose that the latter be designated FABP10, because in our phylogenetic analysis it clustered with zebrafish FABP10. Bioinformatic analysis of across-tissue gene expression patterns in the chicken showed some congruence with phylogenetic relationships. On the basis of expression, chicken FABP genes seemed to form two major groups: (1) a cluster of genes many of which showed predominant expression in the digestive system (FABP1, FABP2, FABP6, FABP10, RBP1, and CRABP1); and (2) a cluster of genes most of which had predominant expression in tissues other than those of the digestive system, including muscle and the central nervous system (FABP3, FABP4, FABP5, FABP7, and PMP2). Since these clusters corresponded to major clusters in the phylogenetic tree as well, it seems a plausible hypothesis that the earliest duplication in the vertebrate FABP family led to the divergence of a gut-specialized gene from a gene expressed mainly in nervous and muscular systems. Data on gene expression in livers of two lines of chickens selected for high growth and low growth showed differences between FABP1 and FABP10 expressions in the liver, supporting the hypothesis of functional divergence between the two chicken liver-expressed FABPs related to food intake.

Novel protein genes in animal mtDNA: a new sex determination system in freshwater mussels (Bivalvia: Unionoida)?

Mitochondrial (mt) function depends critically on optimal interactions between components encoded by mt and nuclear DNAs. mitochondrial DNA (mtDNA) inheritance (SMI) is thought to have evolved in animal species to maintain mito-nuclear complementarity by preventing the spread of selfish mt elements thus typically rendering mtDNA heteroplasmy evolutionarily ephemeral. Here, we show that mtDNA intraorganismal heteroplasmy can have deterministic underpinnings and persist for hundreds of millions of years. We demonstrate that the only exception to SMI in the animal kingdom, that is, the doubly uniparental mtDNA inheritance system in bivalves, with its three-way interactions among egg mt-, sperm mt- and nucleus-encoded gene products, is tightly associated with the maintenance of separate male and female sexes (dioecy) in freshwater mussels. Specifically, this mother-through-daughter and father-through-son mtDNA inheritance system, containing highly differentiated mt genomes, is found in all dioecious freshwater mussel species. Conversely, all hermaphroditic species lack the paternally transmitted mtDNA (=possess SMI) and have heterogeneous macromutations in the recently discovered, novel protein-coding gene (F-orf) in their maternally transmitted mt genomes. Using immunoelectron microscopy, we have localized the F-open reading frame (ORF) protein, likely involved in specifying separate sexes, in mitochondria and in the nucleus. Our results support the hypothesis that proteins coded by the highly divergent maternally and paternally transmitted mt genomes could be directly involved in sex determination in freshwater mussels. Concomitantly, our study demonstrates novel features for animal mt genomes: the existence of additional, lineage-specific, mtDNA-encoded proteins with functional significance and the involvement of mtDNA-encoded proteins in extra-mt functions. Our results open new avenues for the identification, characterization, and functional analyses of ORFs in the intergenic regions, previously defined as "noncoding," found in a large proportion of animal mt genomes.

Molecular characterization and mRNA expression of two key enzymes of hypoxia-sensing pathways in eastern oysters Crassostrea virginica (Gmelin): hypoxia-inducible factor α (HIF-α) and HIF-prolyl hydroxylase (PHD)

Oxygen homeostasis is crucial for development, survival and normal function of all metazoans. A family of transcription factors called hypoxia-inducible factors (HIF) is critical in mediating the adaptive responses to reduced oxygen availability. The HIF transcription factor consists of a constitutively expressed β subunit and an oxygen-dependent α subunit; the abundance of the latter determines the activity of HIF and is regulated by a family of O(2)- and Fe(2+)-dependent enzymes prolyl hydroxylases (PHDs). Currently very little is known about the function of this important pathway and the molecular structure of its key players in hypoxia-tolerant intertidal mollusks including oysters, which are among the animal champions of anoxic and hypoxic tolerance and thus can serve as excellent models to study the role of HIF cascade in adaptations to oxygen deficiency. We have isolated transcripts of two key components of the oxygen sensing pathway - the oxygen-regulated HIF-α subunit and PHD - from an intertidal mollusk, the eastern oyster Crassostrea virginica, and determined the transcriptional responses of these two genes to anoxia, hypoxia and cadmium (Cd) stress. HIF-α and PHD homologs from eastern oysters C. virginica show significant sequence similarity and share key functional domains with the earlier described isoforms from vertebrates and invertebrates. Phylogenetic analysis shows that genetic diversification of HIF and PHD isoforms occurred within the vertebrate lineage indicating functional diversification and specialization of the oxygen-sensing pathways in this group, which parallels situation observed for many other important genes. HIF-α and PHD homologs are broadly expressed at the mRNA level in different oyster tissues and show transcriptional responses to prolonged hypoxia in the gills consistent with their putative role in oxygen sensing and the adaptive response to hypoxia. Similarity in amino acid sequence, domain structure and transcriptional responses between HIF-α and PHD homologs from oysters and other invertebrate and vertebrate species implies the highly conserved functions of these genes throughout the evolutionary history of animals, in accordance with their critical role in oxygen sensing and homeostasis.

Frequent associations between CTL and T-Helper epitopes in HIV-1 genomes and implications for multi-epitope vaccine designs

BACKGROUND

Epitope vaccines have been suggested as a strategy to counteract viral escape and development of drug resistance. Multiple studies have shown that Cytotoxic T-Lymphocyte (CTL) and T-Helper (Th) epitopes can generate strong immune responses in Human Immunodeficiency Virus (HIV-1). However, not much is known about the relationship among different types of HIV epitopes, particularly those epitopes that can be considered potential candidates for inclusion in the multi-epitope vaccines.

RESULTS

In this study we used association rule mining to examine relationship between different types of epitopes (CTL, Th and antibody epitopes) from nine protein-coding HIV-1 genes to identify strong associations as potent multi-epitope vaccine candidates. Our results revealed 137 association rules that were consistently present in the majority of reference and non-reference HIV-1 genomes and included epitopes of two different types (CTL and Th) from three different genes (Gag, Pol and Nef). These rules involved 14 non-overlapping epitope regions that frequently co-occurred despite high mutation and recombination rates, including in genomes of circulating recombinant forms. These epitope regions were also highly conserved at both the amino acid and nucleotide levels indicating strong purifying selection driven by functional and/or structural constraints and hence, the diminished likelihood of successful escape mutations.

CONCLUSIONS

Our results provide a comprehensive systematic survey of CTL, Th and Ab epitopes that are both highly conserved and co-occur together among all subtypes of HIV-1, including circulating recombinant forms. Several co-occurring epitope combinations were identified as potent candidates for inclusion in multi-epitope vaccines, including epitopes that are immuno-responsive to different arms of the host immune machinery and can enable stronger and more efficient immune responses, similar to responses achieved with adjuvant therapies. Signature of strong purifying selection acting at the nucleotide level of the associated epitopes indicates that these regions are functionally critical, although the exact reasons behind such sequence conservation remain to be elucidated.

Analysis of gene expression in prostate cancer epithelial and interstitial stromal cells using laser capture microdissection

BACKGROUND

The prostate gland represents a multifaceted system in which prostate epithelia and stroma have distinct physiological roles. To understand the interaction between stroma and glandular epithelia, it is essential to delineate the gene expression profiles of these two tissue types in prostate cancer. Most studies have compared tumor and normal samples by performing global expression analysis using a mixture of cell populations. This report presents the first study of prostate tumor tissue that examines patterns of differential expression between specific cell types using laser capture microdissection (LCM).

METHODS

LCM was used to isolate distinct cell-type populations and identify their gene expression differences using oligonucleotide microarrays. Ten differentially expressed genes were then analyzed in paired tumor and non-neoplastic prostate tissues by quantitative real-time PCR. Expression patterns of the transcription factors, WT1 and EGR1, were further compared in established prostate cell lines. WT1 protein expression was also examined in prostate tissue microarrays using immunohistochemistry.

RESULTS

The two-step method of laser capture and microarray analysis identified nearly 500 genes whose expression levels were significantly different in prostate epithelial versus stromal tissues. Several genes expressed in epithelial cells (WT1, GATA2, and FGFR-3) were more highly expressed in neoplastic than in non-neoplastic tissues; conversely several genes expressed in stromal cells (CCL5, CXCL13, IGF-1, FGF-2, and IGFBP3) were more highly expressed in non-neoplastic than in neoplastic tissues. Notably, EGR1 was also differentially expressed between epithelial and stromal tissues. Expression of WT1 and EGR1 in cell lines was consistent with these patterns of differential expression. Importantly, WT1 protein expression was demonstrated in tumor tissues and was absent in normal and benign tissues.

CONCLUSIONS

The prostate represents a complex mix of cell types and there is a need to analyze distinct cell populations to better understand their potential interactions. In the present study, LCM and microarray analysis were used to identify novel gene expression patterns in prostate cell populations, including identification of WT1 expression in epithelial cells. The relevance of WT1 expression in prostate cancer was confirmed by analysis of tumor tissue and cell lines, suggesting a potential role for WT1 in prostate tumorigenesis.

Effects of cadmium exposure and intermittent anoxia on nitric oxide metabolism in eastern oysters, Crassostrea virginica

Nitric oxide (NO) is an intracellular signaling molecule synthesized by a group of enzymes called nitric oxide synthases (NOS) and involved in regulation of many cellular functions including mitochondrial metabolism and bioenergetics. In invertebrates, the involvement of NO in bioenergetics and metabolic responses to environmental stress is poorly understood. We determined sensitivity of mitochondrial and cellular respiration to NO and the effects of cadmium (Cd) and intermittent anoxia on NO metabolism in eastern oysters, Crassostrea virginica. NOS activity was strongly suppressed by exposure to 50 μg l–1 Cd for 30 days (4.76 vs 1.19 pmol NO min–1 mg–1 protein in control and Cd-exposed oysters, respectively) and further decreased during anoxic exposure in Cd-exposed oysters but not in their control counterparts. Nitrate/nitrite content (indicative of NO levels) decreased during anoxic exposure to less than 10% of the normoxic values and recovered within 1 h of re-oxygenation in control oysters. In Cd-exposed oysters, the recovery of the normoxic NO levels lagged behind, reflecting their lower NOS activity. Oyster mitochondrial respiration was inhibited by exogenous NO, with sensitivity on a par with that of mammalian mitochondria, and ADP-stimulated mitochondrial respiration was significantly more sensitive to NO than resting respiration. In isolated gill cells, manipulations of endogenous NOS activity either with a specific NOS inhibitor (aminoguanidine) or a NOS substrate (l-arginine) had no effect on respiration, likely due to the fact that mitochondria in the resting state are relatively NO insensitive. Likewise, Cd-induced stimulation of cellular respiration did not correlate with decreased NOS activity in isolated gill cells. High sensitivity of phosphorylating (ADP-stimulated) oyster mitochondria to NO suggests that regulation of bioenergetics is an evolutionarily conserved function of NO and that NO-dependent regulation of metabolism may be most prominent under the conditions of high metabolic flux when the ADP-to-ATP ratio is high.

Discovery of novel targets for multi-epitope vaccines: screening of HIV-1 genomes using association rule mining

Background

Studies have shown that in the genome of human immunodeficiency virus (HIV-1) regions responsible for interactions with the host's immune system, namely, cytotoxic T-lymphocyte (CTL) epitopes tend to cluster together in relatively conserved regions. On the other hand, "epitope-less" regions or regions with relatively low density of epitopes tend to be more variable. However, very little is known about relationships among epitopes from different genes, in other words, whether particular epitopes from different genes would occur together in the same viral genome. To identify CTL epitopes in different genes that co-occur in HIV genomes, association rule mining was used.

Results

Using a set of 189 best-defined HIV-1 CTL/CD8+ epitopes from 9 different protein-coding genes, as described by Frahm, Linde & Brander (2007), we examined the complete genomic sequences of 62 reference HIV sequences (including 13 subtypes and sub-subtypes with approximately 4 representative sequences for each subtype or sub-subtype, and 18 circulating recombinant forms). The results showed that despite inclusion of recombinant sequences that would be expected to break-up associations of epitopes in different genes when two different genomes are recombined, there exist particular combinations of epitopes (epitope associations) that occur repeatedly across the world-wide population of HIV-1. For example, Pol epitope LFLDGIDKA is found to be significantly associated with epitopes GHQAAMQML and FLKEKGGL from Gag and Nef, respectively, and this association rule is observed even among circulating recombinant forms.

Conclusion

We have identified CTL epitope combinations co-occurring in HIV-1 genomes including different subtypes and recombinant forms. Such co-occurrence has important implications for design of complex vaccines (multi-epitope vaccines) and/or drugs that would target multiple HIV-1 regions at once and, thus, may be expected to overcome challenges associated with viral escape.

Cross-species mapping of bidirectional promoters enables prediction of unannotated 5' UTRs and identification of species-specific transcripts

Background

Bidirectional promoters are shared regulatory regions that influence the expression of two oppositely oriented genes. This type of regulatory architecture is found more frequently than expected by chance in the human genome, yet many specifics underlying the regulatory design are unknown. Given that the function of most orthologous genes is similar across species, we hypothesized that the architecture and regulation of bidirectional promoters might also be similar across species, representing a core regulatory structure and enabling annotation of these regions in additional mammalian genomes.

Results

By mapping the intergenic distances of genes in human, chimpanzee, bovine, murine, and rat, we show an enrichment for pairs of genes equal to or less than 1,000 bp between their adjacent 5' ends ("head-to-head") compared to pairs of genes that fall in the same orientation ("head-to-tail") or whose 3' ends are side-by-side ("tail-to-tail"). A representative set of 1,369 human bidirectional promoters was mapped to orthologous sequences in other mammals. We confirmed predictions for 5' UTRs in nine of ten manual picks in bovine based on comparison to the orthologous human promoter set and in six of seven predictions in human based on comparison to the bovine dataset. The two predictions that did not have orthology as bidirectional promoters in the other species resulted from unique events that initiated transcription in the opposite direction in only those species. We found evidence supporting the independent emergence of bidirectional promoters from the family of five RecQ helicase genes, which gained their bidirectional promoters and partner genes independently rather than through a duplication process. Furthermore, by expanding our comparisons from pairwise to multispecies analyses we developed a map representing a core set of bidirectional promoters in mammals.

Conclusion

We show that the orthologous positions of bidirectional promoters provide a reliable guide to directly annotate over one thousand regulatory regions in sequences of mammalian genomes, while also serving as a useful tool to predict 5' UTR positions and identify genes that are novel to a single species.

Molecular characterization and expression of a novel homolog of uncoupling protein 5 (UCP5) from the eastern oyster Crassostrea virginica (Bivalvia: Ostreidae)

Uncoupling proteins (UCPs) belong to the mitochondrial anion carrier gene family which has been implicated in diverse physiological functions ranging from thermoregulation to antioxidant defense. In mammals, the UCP family is well characterized and contains five members (UCP1-5). In contrast, invertebrate homologues of uncoupling proteins are much less studied both from the viewpoints of structure and function. In this study we report nucleotide and predicted protein structure of an important member of UCP family, UCP5 from eastern oysters Crassostrea virginica. UCP5 from oysters appears to be a close homolog of the mammalian brain mitochondrial carrier protein (BMCP1, or UCP5) and is the first full-length UCP described from a Lophotrochozoan invertebrate. Evolutionary analysis of UCP sequences indicates at least three monophyletic UCP branches (UCP1-3, UCP4 and UCP5) that have diverged early in the evolution, prior to the divergence of vertebrates and invertebrates. In oysters, two forms of UCP5 transcript are found (UCP5S and UCP5L) that differ by 152 bp in length due to the presence of an intron in UCP5L. UCP5 was expressed in all studied oyster tissues, unlike mammals, where UCP5 is predominantly expressed in brains and male gonads. Hypoxia-reoxygenation stress, sublethal Cd exposure (50 ?g L(?1) Cd for 56 days) and acclimation to different temperatures (12 and 20 °C) had no significant effect on UCP5 mRNA expression in oysters indicative of its relative unimportance in antioxidant defense and temperature adaptation of oyster mitochondria. These data suggest that despite the relatively high degree of evolutionary conservation of the UCP5 amino acid sequence, its functional significance in mitochondria changed in the course of evolution of mollusks and vertebrates.

Evolutionary conservation of zinc finger transcription factor binding sites in promoters of genes co-expressed with WT1 in prostate cancer

BACKGROUND

Gene expression analyses have led to a better understanding of growth control of prostate cancer cells. We and others have identified the presence of several zinc finger transcription factors in the neoplastic prostate, suggesting a potential role for these genes in the regulation of the prostate cancer transcriptome. One of the transcription factors (TFs) identified in the prostate cancer epithelial cells was the Wilms tumor gene (WT1). To rapidly identify coordinately expressed prostate cancer growth control genes that may be regulated by WT1, we used an in silico approach.

RESULTS

Evolutionary conserved transcription factor binding sites (TFBS) recognized by WT1, EGR1, SP1, SP2, AP2 and GATA1 were identified in the promoters of 24 differentially expressed prostate cancer genes from eight mammalian species. To test the relationship between sequence conservation and function, chromatin of LNCaP prostate cancer and kidney 293 cells were tested for TF binding using chromatin immunoprecipitation (ChIP). Multiple putative TFBS in gene promoters of placental mammals were found to be shared with those in human gene promoters and some were conserved between genomes that diverged about 170 million years ago (i.e., primates and marsupials), therefore implicating these sites as candidate binding sites. Among those genes coordinately expressed with WT1 was the kallikrein-related peptidase 3 (KLK3) gene commonly known as the prostate specific antigen (PSA) gene. This analysis located several potential WT1 TFBS in the PSA gene promoter and led to the rapid identification of a novel putative binding site confirmed in vivo by ChIP. Conversely for two prostate growth control genes, androgen receptor (AR) and vascular endothelial growth factor (VEGF), known to be transcriptionally regulated by WT1, regulatory sequence conservation was observed and TF binding in vivo was confirmed by ChIP.

CONCLUSION

Overall, this targeted approach rapidly identified important candidate WT1-binding elements in genes coordinately expressed with WT1 in prostate cancer cells, thus enabling a more focused functional analysis of the most likely target genes in prostate cancer progression. Identifying these genes will help to better understand how gene regulation is altered in these tumor cells.