Maternal undernutrition results in transcript changes in male offspring that may promote resistance to high fat diet induced weight gain

Maternal nutrition during embryonic development and lactation influences multiple aspects of offspring health. Using mice, this study investigates the effects of maternal caloric restriction (CR) during mid-gestation and lactation on offspring neonatal development and on adult metabolic function when challenged by a high fat diet (HFD). The CR maternal model produced male and female offspring that were significantly smaller, in terms of weight and length, and females had delayed puberty. Adult offspring born to CR dams had a sexually dimorphic response to the high fat diet. Compared to offspring of maternal control dams, adult female, but not male, CR offspring gained more weight in response to high fat diet at 10 weeks. In adipose tissue of male HFD offspring, maternal undernutrition resulted in blunted expression of genes associated with weight gain and increased expression of genes that protect against weight gain. Regardless of maternal nutrition status, HFD male offspring showed increased expression of genes associated with progression toward nonalcoholic fatty liver disease (NAFLD). Furthermore, we observed significant, sexually dimorphic differences in serum TSH. These data reveal tissue- and sex-specific changes in gene and hormone regulation following mild maternal undernutrition, which may offer protection against diet induced weight gain in adult male offspring.

Hydrangea arborescens 'Annabelle' Flower Formation and Flowering in the Current Year

The perennial woody plant Hydrangea arborescens 'Annabelle' is of great research value due to its unique mechanism of flower development that occurs in the current year, resulting in decorative flowers that can be enjoyed for a relatively long period of time. However, the mechanisms underlying the regulation of current-year flower development in H. arborescens 'Annabelle' are still not fully understood. In this study, we conducted an associated analysis to explore the core regulating network in H. arborescens 'Annabelle' by combining phenological observations, physiological assays, and transcriptome comparisons across seven flower developmental stages. Through this analysis, we constructed a gene co-expression network (GCN) based on the highest reciprocal rank (HRR), using 509 differentially expressed genes (DEGs) identified from seven flowering-related pathways, as well as the biosynthesis of eight flowering-related phytohormones and signal transduction in the transcriptomic analysis. According to the analysis of the GCN, we identified 14 key genes with the highest functional connectivity that played critical roles in specific development stages. We confirmed that 135 transcription factors (AP2/ERF, bHLH, CO-like, GRAS, MIKC, SBP, WRKY) were highly co-expressed with the 14 key genes, indicating their close associations with the development of current-year flowers. We further proposed a hypothetical model of a gene regulatory network for the development of the whole flower. This model suggested that the photoperiod, aging, and gibberellin pathways, along with the phytohormones abscisic acid (ABA), gibberellin (GA), brassinosteroid (BR), and jasmonic acid (JA), work synergistically to promote the floral transition. Additionally, auxin, GA, JA, ABA, and salicylic acid (SA) regulated the blooming process by involving the circadian clock. Cytokinin (CTK), ethylene (ETH), and SA were key regulators that affected flower senescence. Additionally, several floral integrators (HaLFY, HaSOC1-2, HaAP1, HaFULL, HaAGL24, HaFLC, etc.) were dominant contributors to the development of H. arborescens flowers. Overall, this research provides a comprehensive understanding of the dynamic mechanism underlying the entire process of current-year flower development, thereby offering valuable insights for further studies on the flower development of H. arborescens 'Annabelle'.

Defining solute carrier transporter signatures of murine immune cell subsets

Solute carrier (SLC) transporters are membrane-bound proteins that facilitate nutrient transport, and the movement across cellular membranes of various substrates ranging from ions to amino acids, metabolites and drugs. Recently, SLCs have gained increased attention due to their functional linkage to innate immunological processes such as the clearance of dead cells and anti-microbial defense. Further, the druggable nature of these transporters provides unique opportunities for improving outcomes in different immunological diseases. Although the SLCs represent the largest group of transporters and are often identified as significant hits in omics data sets, their role in immunology has been insufficiently explored. This is partly due to the absence of tools that allow identification of SLC expression in particular immune cell types and enable their comparison before embarking on functional studies. In this study, we used publicly available RNA-Seq data sets to analyze the transcriptome in adaptive and innate immune cells, focusing on differentially and highly expressed SLCs. This revealed several new insights: first, we identify differentially expressed SLC transcripts in phagocytes (macrophages, dendritic cells, and neutrophils) compared to adaptive immune cells; second, we identify new potential immune cell markers based on SLC expression; and third, we provide user-friendly online tools for researchers to explore SLC genes of interest (and the rest of the genes as well), in three-way comparative dot plots among immune cells. We expect this work to facilitate SLC research and comparative transcriptomic studies across different immune cells.

Comparative transcriptome profiling provides insights into the growth promotion activity of Pseudomonas fluorescens strain SLU99 in tomato and potato plants

The use of biocontrol agents with plant growth-promoting activity has emerged as an approach to support sustainable agriculture. During our field evaluation of potato plants treated with biocontrol rhizobacteria, four bacteria were associated with increased plant height. Using two important solanaceous crop plants, tomato and potato, we carried out a comparative analysis of the growth-promoting activity of the four bacterial strains: Pseudomonas fluorescens SLU99, Serratia plymuthica S412, S. rubidaea AV10, and S. rubidaea EV23. Greenhouse and in vitro experiments showed that P. fluorescens SLU99 promoted plant height, biomass accumulation, and yield of potato and tomato plants, while EV23 promoted growth in potato but not in tomato plants. SLU99 induced the expression of plant hormone-related genes in potato and tomato, especially those involved in maintaining homeostasis of auxin, cytokinin, gibberellic acid and ethylene. Our results reveal potential mechanisms underlying the growth promotion and biocontrol effects of these rhizobacteria and suggest which strains may be best deployed for sustainably improving crop yield.

Major transcriptomic differences are induced by warmer temperature conditions experienced during asexual and sexual reproduction in Fragaria vesca ecotypes

A major challenge for plants in a rapidly changing climate is to adapt to rising temperatures. Some plants adapt to temperature conditions by generating an epigenetic memory that can be transmitted both meiotically and mitotically. Such epigenetic memories may increase phenotypic variation to global warming and provide time for adaptation to occur through classical genetic selection. The goal of this study was to understand how warmer temperature conditions experienced during sexual and asexual reproduction affect the transcriptomes of different strawberry (Fragaria vesca) ecotypes. We let four European F. vesca ecotypes reproduce at two contrasting temperatures (18 and 28°C), either asexually through stolon formation for several generations, or sexually by seeds (achenes). We then analyzed the transcriptome of unfolding leaves, with emphasis on differential expression of genes belonging to the epigenetic machinery. For asexually reproduced plants we found a general transcriptomic response to temperature conditions but for sexually reproduced plants we found less significant responses. We predicted several splicing isoforms for important genes (e.g. a SOC1, LHY, and SVP homolog), and found significantly more differentially presented splicing event variants following asexual vs. sexual reproduction. This difference could be due to the stochastic character of recombination during meiosis or to differential creation or erasure of epigenetic marks during embryogenesis and seed development. Strikingly, very few differentially expressed genes were shared between ecotypes, perhaps because ecotypes differ greatly both genetically and epigenetically. Genes related to the epigenetic machinery were predominantly upregulated at 28°C during asexual reproduction but downregulated after sexual reproduction, indicating that temperature-induced change affects the epigenetic machinery differently during the two types of reproduction.

The transcriptome response of astronaut leukocytes to long missions aboard the International Space Station reveals immune modulation

Introduction

Spaceflight leads to the deconditioning of multiple body systems including the immune system. We sought to characterize the molecular response involved by capturing changes in leukocyte transcriptomes from astronauts transitioning to and from long-duration spaceflight.

Methods

Fourteen male and female astronauts with ~6-month- long missions aboard the International Space Station (ISS) had 10 blood samples collected throughout the three phases of the study: one pre-flight (PF), four in-flight (IF) while onboard the ISS, and five upon return to Earth (R). We measured gene expression through RNA sequencing of leukocytes and applied generalized linear modeling to assess differential expression across all 10 time points followed by the analysis of selected time points and functional enrichment of changing genes to identify shifts in biological processes.

Results

Our temporal analysis identified 276 differentially expressed transcripts grouped into two clusters (C) showing opposite profiles of expression with transitions to and from spaceflight: (C1) decrease-then-increase and (C2) increase-then-decrease. Both clusters converged toward average expression between ~2 and ~6 months in space. Further analysis of spaceflight transitions identified the decrease-then-increase pattern with most changes: 112 downregulated genes between PF and early spaceflight and 135 upregulated genes between late IF and R. Interestingly, 100 genes were both downregulated when reaching space and upregulated when landing on Earth. Functional enrichment at the transition to space related to immune suppression increased cell housekeeping functions and reduced cell proliferation. In contrast, egress to Earth is related to immune reactivation.

Conclusion

The leukocytes' transcriptome changes describe rapid adaptations in response to entering space followed by opposite changes upon returning to Earth. These results shed light on immune modulation in space and highlight the major adaptive changes in cellular activity engaged to adapt to extreme environments.

Integrated transcriptomic and metabolomic analyses reveal key metabolic pathways in response to potassium deficiency in coconut (Cocos nucifera L.) seedlings

Potassium ions (K⁺) are important for plant growth and crop yield. However, the effects of K⁺ deficiency on the biomass of coconut seedlings and the mechanism by which K⁺ deficiency regulates plant growth remain largely unknown. Therefore, in this study, we compared the physiological, transcriptome, and metabolite profiles of coconut seedling leaves under K⁺-deficient and K⁺-sufficient conditions using pot hydroponic experiments, RNA-sequencing, and metabolomics technologies. K⁺ deficiency stress significantly reduced the plant height, biomass, and soil and plant analyzer development value, as well as K content, soluble protein, crude fat, and soluble sugar contents of coconut seedlings. Under K⁺ deficiency, the leaf malondialdehyde content of coconut seedlings were significantly increased, whereas the proline (Pro) content was significantly reduced. Superoxide dismutase, peroxidase, and catalase activities were significantly reduced. The contents of endogenous hormones such as auxin, gibberellin, and zeatin were significantly decreased, whereas abscisic acid content was significantly increased. RNA-sequencing revealed that compared to the control, there were 1003 differentially expressed genes (DEGs) in the leaves of coconut seedlings under K⁺ deficiency. Gene Ontology analysis revealed that these DEGs were mainly related to "integral component of membrane," "plasma membrane," "nucleus", "transcription factor activity," "sequence-specific DNA binding," and "protein kinase activity." Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that the DEGs were mainly involved in "MAPK signaling pathway-plant," "plant hormone signal transduction," "starch and sucrose metabolism," "plant-pathogen interaction," "ABC transporters," and "glycerophospholipid metabolism." Metabolomic analysis showed that metabolites related to fatty acids, lipidol, amines, organic acids, amino acids, and flavonoids were generally down-regulated in coconut seedlings under K⁺ deficiency, whereas metabolites related to phenolic acids, nucleic acids, sugars, and alkaloids were mostly up-regulated. Therefore, coconut seedlings respond to K⁺ deficiency stress by regulating signal transduction pathways, primary and secondary metabolism, and plant-pathogen interaction. These results confirm the importance of K⁺ for coconut production, and provide a more in-depth understanding of the response of coconut seedlings to K⁺ deficiency and a basis for improving K⁺ utilization efficiency in coconut trees.

Diagnosis of childhood and adolescent growth hormone deficiency using transcriptomic data

BACKGROUND

Gene expression (GE) data have shown promise as a novel tool to aid in the diagnosis of childhood growth hormone deficiency (GHD) when comparing GHD children to normal children. The aim of this study was to assess the utility of GE data in the diagnosis of GHD in childhood and adolescence using non-GHD short stature children as a control group.

METHODS

GE data was obtained from patients undergoing growth hormone stimulation testing. Data were taken for the 271 genes whose expression was utilized in our previous study. The synthetic minority oversampling technique was used to balance the dataset and a random forest algorithm applied to predict GHD status.

RESULTS

24 patients were recruited to the study and eight subsequently diagnosed with GHD. There were no significant differences in gender, age, auxology (height SDS, weight SDS, BMI SDS) or biochemistry (IGF-I SDS, IGFBP-3 SDS) between the GHD and non-GHD subjects. A random forest algorithm gave an AUC of 0.97 (95% CI 0.93 - 1.0) for the diagnosis of GHD.

CONCLUSION

This study demonstrates highly accurate diagnosis of childhood GHD using a combination of GE data and random forest analysis.

Comprehensive molecular and morphological resolution of blubber stratification in a deep-diving, fasting-adapted seal

Blubber is a modified subcutaneous adipose tissue in marine mammals that provides energy storage, thermoregulation, hydrodynamic locomotion, and buoyancy. Blubber displays vertical stratification by lipid content, fatty acid composition, and vascularization, leading to the assumption that deeper blubber layers are metabolically active, while superficial layers are mainly structural and thermoregulatory. However, few studies have examined functional stratification of marine mammal blubber directly, especially in pinnipeds. We characterized morphological and transcriptional differences across blubber layers in the northern elephant seal, a deep-diving and fasting-adapted phocid. We collected blubber from seals early in their fasting period and divided blubber cores into three similarly sized portions. We hypothesized that the innermost blubber portion would have higher 1) heterogeneity in adipocyte size, 2) microvascular density, and 3) expression of genes associated with metabolism and hormone signaling than outer blubber. We found that adipocyte area and variance increased from outermost (skin-adjacent) to innermost (muscle-adjacent) blubber layers, suggesting that inner blubber has a higher capacity for lipid storage and turnover than outer blubber. Inner blubber had a higher proportion of CD144+ endothelial cells, suggesting higher microvascular density. In contrast, outer blubber had a higher proportion of CD4+ immune cells than inner blubber, suggesting higher capacity for response to tissue injury. Transcriptome analysis identified 61 genes that were differentially expressed between inner and outer blubber layers, many of which have not been studied previously in marine mammals. Based on known functions of these genes in other mammals, we suggest that inner blubber has potentially higher 1) adipogenic capacity, 2) cellular diversity, and 3) metabolic and neuroendocrine signaling activity, while outer blubber may have higher 1) extracellular matrix synthesis activity and 2) responsiveness to pathogens and cell stressors. We further characterized expression of nine genes of interest identified by transcriptomics and two adipokines with higher precision across blubber layers using targeted assays. Our study provides functional insights into stratification of blubber in marine mammals and a molecular key, including CD144, CD4, HMGCS2, GABRG2, HCAR2, and COL1A2, for distinguishing blubber layers for physiological and functional studies in seals.

Dynamic landscape of long noncoding RNAs during leaf aging in Arabidopsis

Leaf senescence, the last stage of leaf development, is essential for whole-plant fitness as it marks the relocation of nutrients from senescing leaves to reproductive or other developing organs. Temporally coordinated physiological and functional changes along leaf aging are fine-tuned by a highly regulated genetic program involving multi-layered regulatory mechanisms. Long noncoding RNAs (lncRNAs) are newly emerging as hidden players in many biological processes; however, their contribution to leaf senescence has been largely unknown. Here, we performed comprehensive analyses of RNA-seq data representing all developmental stages of leaves to determine the genome-wide lncRNA landscape along leaf aging. A total of 771 lncRNAs, including 232 unannotated lncRNAs, were identified. Time-course analysis revealed 446 among 771 developmental age-related lncRNAs (AR-lncRNAs). Intriguingly, the expression of AR-lncRNAs was regulated more dynamically in senescing leaves than in growing leaves, revealing the relevant contribution of these lncRNAs to leaf senescence. Further analyses enabled us to infer the function of lncRNAs, based on their interacting miRNA or mRNA partners. We considered functionally diverse lncRNAs including antisense lncRNAs (which regulate overlapping protein-coding genes), competitive endogenous RNAs (ceRNAs; which regulate paired mRNAs using miRNAs as anchors), and mRNA-interacting lncRNAs (which affect the stability of mRNAs). Furthermore, we experimentally validated the senescence regulatory function of three novel AR-lncRNAs including one antisense lncRNA and two mRNA-interacting lncRNAs through molecular and phenotypic analyses. Our study provides a valuable resource of AR-lncRNAs and potential regulatory networks that link the function of coding mRNA and AR-lncRNAs. Together, our results reveal AR-lncRNAs as important elements in the leaf senescence process.

Exposure of Larvae to Sublethal Thiacloprid Delays Bee Development and Affects Transcriptional Responses of Newly Emerged Honey Bees

Understanding the cause of honey bee (Apis mellifera) population decline has attracted immense attention worldwide in recent years. Exposure to neonicotinoid pesticides is considered one of the most probable factors due to the physiological and behavioral damage they cause to honey bees. However, the influence of thiacloprid, a relatively less toxic cyanogen-substituted form of neonicotinoid, on honey bee (Apis mellifera L.) development is not well studied. The toxicity of sublethal thiacloprid to larvae, pupae, and emerging honey bees was assessed under laboratory conditions. We found that thiacloprid reduced the survival rate of larvae and pupae, and delayed the development of bees which led to lower bodyweight and size. Furthermore, we identified differentially expressed genes involved in metabolism and immunity though RNA-sequencing of newly-emerged adult bees. GO enrichment analysis identified genes involved in metabolism, catalytic activity, and transporter activity. KEGG pathway analysis indicated that thiacloprid induced up-regulation of genes related to glutathione metabolism and Toll-like receptor signaling pathway. Overall, our results suggest that chronic sublethal thiacloprid can affect honey bee colonies by reducing survival and delaying bee development.

Transcriptome and Metabolome Profiling of a Novel Isolate Chlorella sorokiniana G32 (Chlorophyta) Displaying Enhanced Starch Accumulation at High Growth Rate Under Mixotrophic Condition

Chlorella sorokiniana is one of the most productive microalgal species with a high potential for the production of biofuels and other high value-added molecules. Many studies have focused on its capability of mixotrophic growth using reduced organic carbon and growth pattern shift between autotrophic and mixotrophic conditions. In this study, we investigated growth patterns of a novel isolate, C. sorokiniana G32, under mixotrophic growth conditions supplemented with a low level (1.25 g L-1) and a high level (5 g L-1) of glucose. Physiological, transcriptomic (i.e., RNA-seq), and metabolomic (i.e., LC-MS/MS) methods were used. We showed that peak growth based on OD680nm absorbance is ∼4-fold higher with high glucose vs. low glucose supplementation. Photosynthetic efficiency (Fv/Fm) in G32 mixotrophic cultures with high or low glucose supplementation remains identical to that of G32 phototrophic growth. We also found that the conversion rate between absorbance-based cell density and cell dry weight with high glucose supplementation was lower than with low glucose. This suggests that more cell biomass is produced under high glucose treatment than with low glucose. The result was confirmed via sucrose density gradient centrifugation. It is likely that accumulation of high concentration of starch may account for this effect. Transcriptomic analysis of G32 cultures (i.e., via RNA-seq) in response to reciprocal change of glucose levels reveals that expression of a subset of differentially expressed genes (DEGs) is correlated with the amount of glucose supplementation. These DEGs are designated as glucose-specific responsive (GSR) genes. GSR genes are enriched for a number of energy metabolic pathways. Together with metabolomics data (i.e., LC-MS/MS), we show that under high-level supplementation, glucose is preferentially oxidized through an oxidative pentose phosphate pathway. Collectively, our results indicate the mechanism of regulation of glucose assimilation and energy metabolism in G32 under mixotrophic conditions with different levels of glucose supplementation revealed by transcriptomic and metabolomic analyses. We propose that C. sorokiniana G32 has the potential for the production of high value-added molecules.

Comparative transcriptome analysis of Citrus macrophylla tree infected with Citrus tristeza virus stem pitting mutants provides new insight into the role of phloem regeneration in stem pitting disease

Stem pitting is a complex and economically important virus-associated disease of perennial woody plants. Molecular mechanisms and pathways occurring during virus-plant interaction that result in this phenomenon are still obscure. Previous studies indicated that different Citrus tristeza virus (CTV) mutants induce defined stem pitting phenotypes ranging from mild (CTVΔp13) to severe (CTVΔp33) in Citrus macrophylla trees. In this study, we conducted comparative transcriptome analyses of C. macrophylla trees infected with CTV mutants (CTVΔp13 and CTVΔp33) and a full-length virus in comparison to healthy plants as control. The mild CTV stem pitting mutant had very few differentially expressed genes (DEGs) related to plant defense mechanism and plant growth and development. In contrast, substantial gene expression changes were observed in plants infected with the severe mutant and the full-length virus, indicating that both the p13 and p33 proteins of CTV acted as a regulator of symptom production by activating and modulating plant responses, respectively. The analysis of transcriptome data for CTVΔp33 and the full-length virus suggested that xylem specification has been blocked by detecting several genes encoding xylem, cell wall and lignin degradation, and cell wall loosening enzymes. Furthermore, stem pitting was accompanied by downregulation of transcription factors involved in regulation of xylem differentiation and downregulation of some genes involved in lignin biosynthesis, showing that the xylem differentiation and specification program has been shut off. Upregulation of genes encoding transcription factors associated with phloem and cambium development indicated the activation of this program in stem pitting disease. Furthermore, we detected the induction of several DEGs encoding proteins associated with cell cycle re-entry such as chromatin remodeling factors and cyclin, and histone modification. This kind of expression pattern of genes related to xylem differentiation and specification, phloem and cambium development, and cell cycle re-entry is demonstrated during secondary vascular tissue (SVT) regeneration. The microscopy analysis confirmed that the regeneration of new phloem is associated with stem pitting phenotypes. The findings of this study, thus, provide evidence for the association between stem pitting phenotypes and SVT regeneration, suggesting that the expression of these genes might play important roles in development of stem pitting symptoms. Overall, our findings suggest that phloem regeneration contributes to development of stem pitting symptoms.

Transcriptional and Epigenetic Regulation of Monocyte and Macrophage Dysfunction by Chronic Alcohol Consumption

Drinking alcohol, even in moderation, can affect the immune system. Studies have shown disproportionate effects of alcohol on circulating and tissue-resident myeloid cells (granulocytes, monocytes, macrophages, dendritic cells). These cells orchestrate the body's first line of defense against microbial challenges as well as maintain tissue homeostasis and repair. Alcohol's effects on these cells are dependent on exposure pattern, with acute drinking dampening but chronic drinking enhancing production of inflammatory mediators. Although chronic drinking is associated with heightened systemic inflammation, studies on tissue resident macrophage populations in several organs including the spleen, liver, brain, and lung have also shown compromised functional and metabolic capacities of these cells. Many of these effects are thought to be mediated by oxidative stress caused by alcohol and its metabolites which can directly impact the cellular epigenetic landscapes. In addition, since myeloid cells are relatively short-lived in circulation and are under constant repopulation from the bone marrow compartment, alcohol's effects on bone marrow progenitors and hematopoiesis are important for understanding the impact of alcohol systemically on these myeloid populations. Alcohol-induced disruption of progenitor, circulating, and tissue resident myeloid populations contribute to the increased susceptibility of patients with alcohol use disorders to viral and bacterial infections. In this review, we provide an overview of the impact of chronic alcohol consumption on the function of monocytes and macrophages in host defense, tissue repair and inflammation. We then summarize our current understanding of the mechanisms underlying alcohol-induced disruption and examine changes in transcriptome and epigenome of monocytes and mcrophages. Overall, chronic alcohol consumption leads to hyper-inflammation concomitant with decreased microbial and wound healing responses by monocytes/macrophages due to a rewiring of the epigentic and transcriptional landscape. However, in advanced alcoholic liver disease, myeloid cells become immunosuppressed as a response to the surrounding hyper-inflammatory milieu. Therefore, the effect of chronic alcohol on the inflammatory response depends on disease state and the immune cell population.

miRNA transcriptome and myofiber characteristics of lamb skeletal muscle during hypertrophic growth1

Postnatal muscle growth is achieved through hypertrophy of the muscle fibers and is impacted by the activity of satellite cells, the quiescent muscle stem cell. Several miRNAs are preferentially expressed in skeletal muscle and could provide a mechanism for increasing muscle hypertrophy through satellite cell proliferation and/or differentiation. The objectives of this study were to: 1) Characterize the miRNA transcriptome of the longissimus thoracis et lumborum muscle at several developmental timepoints [gestational d 85 (PN1), 110 (PN2), 133 (PN3), postnatal d 42 (PW1), 65 (PW2), 243 (MAT)] during muscle hypertrophy in lambs, and 2) examine miR-29a, identified in sequencing to be differentially regulated across development, loss of function on satellite cell proliferation and differentiation. Muscle fiber characteristics showed drastic increases (p < 0.0001) in fiber size and alterations in muscle fiber type occur during pre and postnatal development. miRNA sequencing comparisons were performed in developmental order (PN1 vs. PN2, PN2 vs. PN3, PN3 vs. PW1, PW1 vs. PW2, PW2 vs. MAT). There were 184 differentially expressed (P _adj < 0.05) miRNA, 142 unique miRNA, from all 5 comparisons made. The transitional stage (PN3 vs. PW1) had the largest number (115) of differentially expressed miRNA. Inhibition of miR-29a in satellite cell culture increased (p < 0.05) cell proliferation and differentiation capacity. Characterization of the miRNA transcriptome provides valuable insights into the miRNA involved in muscle fiber hypertrophy and the potential importance of the transitional period.

Transcriptomic Remodelling of Fetal Endothelial Cells During Establishment of Inflammatory Memory

Inflammatory memory involves the molecular and cellular ‘reprogramming’ of innate immune cells following exogenous stimuli, leading to non-specific protection against subsequent pathogen exposure. This phenomenon has now also been described in non-hematopoietic cells, such as human fetal and adult endothelial cells. In this study we mapped the cell-specific DNA methylation profile and the transcriptomic remodelling during the establishment of inflammatory memory in two distinct fetal endothelial cell types – a progenitor cell (ECFC) and a differentiated cell (HUVEC) population. We show that both cell types have a core transcriptional response to an initial exposure to a viral-like ligand, Poly(I:C), characterised by interferon responsive genes. There was also an ECFC specific response, marked by the transcription factor ELF1, suggesting a non-canonical viral response pathway in progenitor endothelial cells. Next, we show that both ECFCs and HUVECs establish memory in response to an initial viral exposure, resulting in an altered subsequent response to lipopolysaccharide. While the capacity to train or tolerize the induction of specific sets of genes was similar between the two cell types, the progenitor ECFCs show a higher capacity to establish memory. Among tolerized cellular pathways are those involved in endothelial barrier establishment and leukocyte migration, both important for regulating systemic immune-endothelial cell interactions. These findings suggest that the capacity for inflammatory memory may be a common trait across different endothelial cell types but also indicate that the specific downstream targets may vary by developmental stage.

Unraveling the Metabolic Derangements Occurring in Non-infarcted Areas of Pig Hearts With Chronic Heart Failure

Objective: After myocardial infarction (MI), the non-infarcted left ventricle (LV) ensures appropriate contractile function of the heart. Metabolic disturbance in this region greatly exacerbates post-MI heart failure (HF) pathology. This study aimed to provide a comprehensive understanding of the metabolic derangements occurring in the non-infarcted LV that could trigger cardiovascular deterioration. Methods and Results: We used a pig model that progressed into chronic HF over 3 months following MI induction. Integrated gene and metabolite signatures revealed region-specific perturbations in amino acid- and lipid metabolism, insulin signaling and, oxidative stress response. Remote LV, in particular, showed impaired glutamine and arginine metabolism, altered synthesis of lipids, glucose metabolism disorder, and increased insulin resistance. LPIN1, PPP1R3C, PTPN1, CREM, and NR0B2 were identified as the main effectors in metabolism dysregulation in the remote zone and were found differentially expressed also in the myocardium of patients with ischemic and/or dilated cardiomyopathy. In addition, a simultaneous significant decrease in arginine levels and altered PRCP, PTPN1, and ARF6 expression suggest alterations in vascular function in remote area. Conclusions: This study unravels an array of dysregulated genes and metabolites putatively involved in maladaptive metabolic and vascular remodeling in the non-infarcted myocardium and may contribute to the development of more precise therapies to mitigate progression of chronic HF post-MI.

LeuO, a LysR-Type Transcriptional Regulator, Is Involved in Biofilm Formation and Virulence of Acinetobacter baumannii

Acinetobacter baumannii is an important nosocomial pathogen that can survive in different environmental conditions and poses a severe threat to public health due to its multidrug resistance properties. Research on transcriptional regulators, which play an essential role in adjusting to new environments, could provide new insights into A. baumannii pathogenesis. LysR-type transcriptional regulators (LTTRs) are structurally conserved among bacterial species and regulate virulence in many pathogens. We identified a novel LTTR, designated as LeuO encoded in the A. baumannii genome. After construction of LeuO mutant strain, transcriptome analysis showed that LeuO regulates the expression of 194 upregulated genes and 108 downregulated genes responsible for various functions and our qPCR validation of several differentially expressed genes support transcriptome data. Our results demonstrated that disruption of LeuO led to increased biofilm formation and increased pathogenicity in an animal model. However, the adherence and surface motility of the LeuO mutant were reduced compared with those of the wild-type strain. We observed some mutations on amino acids sequence of LeuO in clinical isolates. These mutations in the A. baumannii biofilm regulator LeuO may cause hyper-biofilm in the tested clinical isolates. This study is the first to demonstrate the association between the LTTR member LeuO and virulence traits of A. baumannii.

Transcriptome-Based Network Analysis Reveals Hirudin Potentiates Anti-Renal Fibrosis Efficacy in UUO Rats

Background: Hirudin has been widely used in the treatment of antifibrosis. Previous studies have shown that hirudin can effectively improve the clinical remission rate of chronic kidney disease. However, the mechanism of its renal protection has not been systematically investigated.

Methods: In this study, the reliability of UUO-induced renal interstitial fibrosis was evaluated by histopathological verification. High-throughput transcriptome sequencing was used to elucidate the molecular mechanism of hirudin, differentially expressed mRNAs were identified, and their functions were analyzed by GO analysis and GSEA. In addition, the RNA-seq results were validated by in vitro and vivo experiments.

Results: We found 322 identical differential expressed genes (IDEs) in the UUO hirudin-treated group compared with the sham group. Functional enrichment analysis indicated that cellular amino acid metabolic processes were the most obvious enrichment pathways in biological processes. In terms of molecular functional enrichment analysis, IDEs were mainly enriched in coenzyme binding, pyridoxal phosphate binding and other pathways. In addition, microbody is the most obvious pathway for cellular components. A total of 115 signaling pathways were enriched, and AMPK, JAK-STAT, and PI3K-Akt signaling pathways were the important signaling pathways enriched. We found that PI3K, p-Akt, and mTOR expression were significantly reduced by hirudin treatment. In particular, our results showed that hirudin could induce a decrease in the expression of autophagy-related proteins such as P62, LC3, Beclin-1 in TGF-β1-induced NRK-52E cells.

Conclusion: Our results suggest that hirudin may protect the kidney by ameliorating renal autophagy impairment through modulating the PI3K/Akt pathway.

Transcriptome Analysis Identifies Strategies Targeting Immune Response-Related Pathways to Control Enterotoxigenic Escherichia coli Infection in Porcine Intestinal Epithelial Cells

Enterotoxigenic Escherichia coli (ETEC) is an important cause of post-weaning diarrhea (PWD) worldwide, resulting in huge economic losses to the swine industry worldwide. In this study, to understand the pathogenesis, the transcriptomic analysis was performed to explore the biological processes (BP) in porcine intestinal epithelial J2 cells infected with an emerging ETEC strain isolated from weaned pigs with diarrhea. Under the criteria of |fold change| (FC) ≥ 2 and P < 0.05 with false discovery rate < 0.05, a total of 131 referenced and 19 novel differentially expressed genes (DEGs) were identified after ETEC infection, including 96 upregulated DEGs and 54 downregulated DEGs. The Gene Ontology (GO) analysis of DEGs showed that ETEC evoked BP specifically involved in response to lipopolysaccharide (LPS) and negative regulation of intracellular signal transduction. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that immune response-related pathways were mainly enriched in J2 cells after ETEC infection, in which tumor necrosis factor (TNF), interleukin 17, and mitogen-activated protein kinase (MAPK) signaling pathways possessed the highest rich factor, followed by nucleotide-binding and oligomerization domain-like receptor (NLRs), C-type lectin receptor (CLR), cytokine–cytokine receptor interaction, and Toll-like receptor (TLR), and nuclear factor kappa-B (NF-κB) signaling pathways. Furthermore, 30 of 131 referenced DEGs, especially the nuclear transcription factor AP-1 and NF-κB, participate in the immune response to infection through an integral signal cascade and can be target molecules for prevention and control of enteric ETEC infection by probiotic Lactobacillus reuteri. Our data provide a comprehensive insight into the immune response of porcine intestinal epithelial cells (IECs) to ETEC infection and advance the identification of targets for prevention and control of ETEC-related PWD.

Transcriptional Time Course After Rotator Cuff Tear

Rotator cuff (RC) tears are prevalent in the population above the age of 60. The disease progression leads to muscle atrophy, fibrosis, and fatty infiltration in the chronic state, which is not improved with intervention or surgical repair. This highlights the need to better understand the underlying dysfunction in muscle after RC tendon tear. Contemporary studies aimed at understanding muscle pathobiology after RC tear have considered transcriptional data in mice, rats and sheep models at 2–3 time points (1 to 16 weeks post injury). However, none of these studies observed a transition or resurgence of gene expression after the initial acute time points. In this study, we collected rabbit supraspinatus muscle tissue with high temporal resolution (1, 2, 4, 8, and 16 weeks) post-tenotomy (n = 6/group), to determine if unique, time-dependent transcriptional changes occur. RNA sequencing and analyses were performed to identify a transcriptional timeline of RC muscle changes and related morphological sequelae. At 1-week post-tenotomy, the greatest number of differentially expressed genes was observed (1,069 up/873 down) which decreases through 2 (170/133), 4 (86/41), and 8 weeks (16/18), followed by a resurgence and transition of expression at 16 weeks (1,421/293), a behavior which previously has not been captured or reported. Broadly, 1-week post-tenotomy is an acute time point with expected immune system responses, catabolism, and changes in energy metabolism, which continues into 2 weeks with less intensity and greater contribution from mitochondrial effects. Expression shifts at 4 weeks post-tenotomy to fatty acid oxidation, lipolysis, and general upregulation of adipogenesis related genes. The effects of previous weeks’ transcriptional dysfunction present themselves at 8 weeks post-tenotomy with enriched DNA damage binding, aggresome activity, extracellular matrix-receptor changes, and significant expression of genes known to induce apoptosis. At 16 weeks post-tenotomy, there is a range of enriched pathways including extracellular matrix constituent binding, mitophagy, neuronal activity, immune response, and more, highlighting the chaotic nature of this time point and possibility of a chronic classification. Transcriptional activity correlated significantly with histological changes and were enriched for biologically relevant pathways such as lipid metabolism. These data provide platform for understanding the biological mechanisms of chronic muscle degeneration after RC tears.

Coordinated Expressional Landscape of the Human Placental miRNome and Transcriptome

Placenta is a unique organ that serves its own function, and contributes to maternal gestational adaptation and fetal development. Coordination of its transcriptome to satisfy all the maternal-fetal needs across gestation is not fully understood. MicroRNAs are powerful transcriptome modulators capable to adjust rapidly the expression level and dynamics of large gene sets. This MiR-Seq based study presents a multi-omics investigation of the human placental miRNome and its synergy with the transcriptome. The analysis included 52 placentas representing three trimesters of normal pregnancy, and term cases of late-onset preeclampsia (LO-PE), gestational diabetes and affected fetal growth. Gestational-age dependent differential expression (FDR < 0.05) was detected for 319 of 417 tested miRNAs (76.5%). A shared list of target genes of dynamic miRNAs suggested their coordinated action. The most abundant miR-143-3p revealed as a marker for pregnancy progression. The data suggested critical, but distinct roles of placenta-specific imprinted C19MC and C14MC miRNA clusters. Paternally encoded primate-specific C19MC was highly transcribed during first trimester, potentially fine-tuning the early placental transcriptome in dosage-sensitive manner. Maternally encoded eutherian C14MC showed high expression until term, underlining its key contribution across gestation. A major shift in placental miRNome (16% miRNAs) was observed in LO-PE, but not in other term pregnancy complications. Notably, 13/38 upregulated miRNAs were transcribed from C19MC and only one from C14MC, whereas 11/28 downregulated miRNAs represented C14MC and none C19MC. miR-210-3p, miR-512-5p, miR-32-5p, miR-19a-3p, miR-590-3p, miR-379-5p were differentially expressed in LO-PE and cases of small-for-gestational-age newborns, supporting a shared etiology. Expression correlation analysis with the RNA-Seq data (16,567 genes) of the same samples clustered PE-linked miRNAs into five groups. Large notable clusters of miRNA–gene pairs showing directly and inversely correlated expression dynamics suggested potential functional relationships in both scenarios. The first genome-wide study of placental miR-eQTLs identified 66 placental SNVs associated with the expression of neighboring miRNAs, including PE-linked miRNAs miR-30a-5p, miR-210-3p, miR-490-3p and miR-518-5p. This study provided a rich catalog of miRNAs for further in-depth investigations of their individual and joint effect on placental transcriptome. Several highlighted miRNAs may serve as potential biomarkers for pregnancy monitoring and targets to prevent or treat gestational disorders.

The Transcriptome of Paired Major and Minor Salivary Gland Tissue in Patients With Primary Sjögren's Syndrome

Background

While all salivary glands (SGs) can be involved in primary Sjögren’s syndrome (pSS), their respective role in pathogenesis remains unclear. Our objective was to assess immunopathway activation in paired parotid and labial gland tissue from biopsy-positive and biopsy-negative pSS and non-SS sicca patients.

Methods

Paraffin-embedded, paired parotid and labial salivary gland tissue and peripheral blood mononuclear cells were obtained from 39 pSS and 20 non-SS sicca patients. RNA was extracted, complementary DNA libraries were prepared and sequenced. For analysis of differentially expressed genes (DEGs), patients were subdivided based on fulfillment of ACR-EULAR criteria and histopathology.

Results

With principal component analysis, only biopsy-positive pSS could be separated from non-SS sicca patients based on SG gene expression. When comparing the transcriptome of biopsy-positive pSS and biopsy-negative non-SS sicca patients, 1235 and 624 DEGs (FDR<0.05, log2FC<-1 or >1) were identified for parotid and labial glands, respectively. The number of DEGs between biopsy-negative pSS and non-SS sicca patients was scarce. Overall, transcript expression levels correlated strongly between parotid and labial glands (R² = 0.86, p-value<0.0001). Gene signatures present in both glands of biopsy-positive pSS patients included IFN-α signaling, IL-12/IL-18 signaling, CD3/CD28 T-cell activation, CD40 signaling in B-cells, DN2 B-cells, and FcRL4+ B-cells. Signature scores varied considerably amongst pSS patients.

Conclusion

Transcriptomes of paired major and minor SGs in pSS were overall comparable, although significant inter-individual heterogeneity in immunopathway activation existed. The SG transcriptome of biopsy-negative pSS was indistinguishable from non-SS sicca patients. Different patterns of SG immunopathway activation in pSS argue for personalized treatment approaches.

Scale Drop Disease Virus (SDDV) and Lates calcarifer Herpes Virus (LCHV) Coinfection Downregulate Immune-Relevant Pathways and Cause Splenic and Kidney Necrosis in Barramundi Under Commercial Farming Conditions

Marine farming of barramundi (Lates calcarifer) in Southeast Asia is currently severely affected by viral diseases. To better understand the biological implications and gene expression response of barramundi in commercial farming conditions during a disease outbreak, the presence of pathogens, comparative RNAseq, and histopathology targeting multiple organs of clinically “sick” and “healthy” juveniles were investigated. Coinfection of scale drop disease virus (SDDV) and L. calcarifer herpes virus (LCHV) were detected in all sampled fish, with higher SDDV viral loads in sick than in healthy fish. Histopathology showed that livers in sick fish often had moderate to severe abnormal fat accumulation (hepatic lipidosis), whereas the predominant pathology in the kidneys shows moderate to severe inflammation and glomerular necrosis. The spleen was the most severely affected organ, with sick fish presenting severe multifocal and coalescing necrosis. Principal component analysis (PC1 and PC2) explained 70.3% of the observed variance and strongly associated the above histopathological findings with SDDV loads and with the sick phenotypes, supporting a primary diagnosis of the fish being impacted by scale drop disease (SDD). Extracted RNA from kidney and spleen of the sick fish were also severely degraded likely due to severe inflammation and tissue necrosis, indicating failure of these organs in advanced stages of SDD. RNAseq of sick vs. healthy barramundi identified 2,810 and 556 differentially expressed genes (DEGs) in the liver and muscle, respectively. Eleven significantly enriched pathways (e.g., phagosome, cytokine-cytokine-receptor interaction, ECM-receptor interaction, neuroactive ligand-receptor interaction, calcium signaling, MAPK, CAMs, etc.) and gene families (e.g., tool-like receptor, TNF, lectin, complement, interleukin, chemokine, MHC, B and T cells, CD molecules, etc.) relevant to homeostasis and innate and adaptive immunity were mostly downregulated in sick fish. These DEGs and pathways, also previously identified in L. calcarifer as general immune responses to other pathogens and environmental stressors, suggest a failure of the clinically sick fish to cope and overcome the systemic inflammatory responses and tissue degeneration caused by SDD.

Plasma Testosterone Levels and Atherosclerotic Plaque Gene Expression in Men With Advanced Atherosclerosis

Aims: Low plasma testosterone levels have been shown to predict worse outcome in men with severe atherosclerotic disease. We hypothesized that a low plasma testosterone level affects disease risk through changes in gene expression in atherosclerotic plaques. Therefore, we studied plasma testosterone levels in relation to gene expression levels in atherosclerotic plaque tissue of men with advanced atherosclerotic disease. Methods: Plasma testosterone levels were measured in 203 men undergoing carotid endarterectomy. The corresponding atherosclerotic plaque tissue was used for RNA sequencing. First, we assessed how often the androgen receptor gene was expressed in the plaque. Second, correlations between plasma testosterone levels and pre-selected testosterone-sensitive genes were assessed. Finally, differences within the RNA expression profile of the plaque as a whole, characterized into gene regulatory networks and at individual gene level were assessed in relation to testosterone levels. Results: Testosterone plasma levels were low with a median of 11.6 nmol/L (IQR: 8.6-13.8). RNA-seq of the plaque resulted in reliable expression data for 18,850 genes to be analyzed. Within the RNA seq data, the androgen-receptor gene was expressed in 189 out of 203 (93%) atherosclerotic plaques of men undergoing carotid endarterectomy. The androgen receptor gene expression was not associated with testosterone plasma levels. There were no significant differences in gene expression of atherosclerotic plaques between different endogenous testosterone levels. This remained true for known testosterone-sensitive genes, the complete transcriptomic profile, male-specific gene co-expression modules as well as for individual genes. Conclusion: In men with severe atherosclerotic disease the androgen receptor is highly expressed in plaque tissue. However, plasma testosterone levels were neither associated with pre-selected testosterone sensitive genes, gene expression profiles nor gene regulatory networks in late-stage atherosclerotic plaques. The effect of testosterone on gene expression of the late-stage atherosclerotic plaque appears limited, suggesting that alternate mechanisms explain its effect on clinical outcomes.

Analysis of Gastric Cancer Transcriptome Allows the Identification of Histotype Specific Molecular Signatures With Prognostic Potential

Gastric cancer is the fifth most common malignancy but the third leading cause of cancer-associated mortality worldwide. Therapy for gastric cancer remain largely suboptimal making the identification of novel therapeutic targets an urgent medical need. In the present study we have carried out a high-throughput sequencing of transcriptome expression in patients with gastric cancers. Twenty-four patients, among a series of 53, who underwent an attempt of curative surgery for gastric cancers in a single center, were enrolled. Patients were sub-grouped according to their histopathology into diffuse and intestinal types, and the transcriptome of the two subgroups assessed by RNAseq analysis and compared to the normal gastric mucosa. The results of this investigation demonstrated that the two histopathology phenotypes express two different patterns of gene expression. A total of 2,064 transcripts were differentially expressed between neoplastic and non-neoplastic tissues: 772 were specific for the intestinal type and 407 for the diffuse type. Only 885 transcripts were simultaneously differentially expressed by both tumors. The per pathway analysis demonstrated an enrichment of extracellular matrix and immune dysfunction in the intestinal type including CXCR2, CXCR1, FPR2, CARD14, EFNA2, AQ9, TRIP13, KLK11 and GHRL. At the univariate analysis reduced levels AQP9 was found to be a negative predictor of 4 years survival. In the diffuse type low levels CXCR2 and high levels of CARD14 mRNA were negative predictors of 4 years survival. In summary, we have identified a group of genes differentially regulated in the intestinal and diffuse histotypes of gastric cancers with AQP9, CARD14 and CXCR2 impacting on patients' prognosis, although CXCR2 is the only factor independently impacting overall survival.

Molecular Programming of Drought-Challenged Trichoderma harzianum-Bioprimed Rice (Oryza sativa L.)

Trichoderma biopriming enhances rice growth in drought-stressed soils by triggering various plant metabolic pathways related to antioxidative defense, secondary metabolites, and hormonal upregulation. In the present study, transcriptomic analysis of rice cultivar IR64 bioprimed with Trichoderma harzianum under drought stress was carried out in comparison with drought-stressed samples using next-generation sequencing techniques. Out of the 2,506 significant (p < 0.05) differentially expressed genes (DEGs), 337 (15%) were exclusively expressed in drought-stressed plants, 382 (15%) were expressed in T. harzianum-treated drought-stressed plants, and 1,787 (70%) were commonly expressed. Furthermore, comparative analysis of upregulated and downregulated genes under stressed conditions showed that 1,053 genes (42%) were upregulated and 733 genes (29%) were downregulated in T. harzianum-treated drought-stressed rice plants. The genes exclusively expressed in T. harzianum-treated drought-stressed plants were mostly photosynthetic and antioxidative such as plastocyanin, small chain of Rubisco, PSI subunit Q, PSII subunit PSBY, osmoproteins, proline-rich protein, aquaporins, stress-enhanced proteins, and chaperonins. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis states that the most enriched pathways were metabolic (38%) followed by pathways involved in the synthesis of secondary metabolites (25%), carbon metabolism (6%), phenyl propanoid (7%), and glutathione metabolism (3%). Some of the genes were selected for validation using real-time PCR which showed consistent expression as RNA-Seq data. Furthermore, to establish host-T. harzianum interaction, transcriptome analysis of Trichoderma was also carried out. The Gene Ontology (GO) analysis of T. harzianum transcriptome suggested that the annotated genes are functionally related to carbohydrate binding module, glycoside hydrolase, GMC oxidoreductase, and trehalase and were mainly upregulated, playing an important role in establishing the mycelia colonization of rice roots and its growth. Overall, it can be concluded that T. harzianum biopriming delays drought stress in rice cultivars by a multitude of molecular programming.

Enhanced Glycolysis Is Required for Antileishmanial Functions of Neutrophils Upon Infection With Leishmania donovani

Visceral leishmaniasis (VL) is a fatal parasitic disease if untreated. Treatment options of VL diminish due to emerging drug resistance. Although the principal host cells for the multiplication of Leishmania are macrophages, neutrophils are the first cells infected with the parasites rapidly after parasite inoculation. Leishmania can survive in neutrophils despite the potent antimicrobial effector functions of neutrophils that can eliminate the parasites. Recently, the growing field of immunometabolism provided strong evidence for the therapeutic potential in targeting metabolic processes as a means of controlling immune effector functions. Therefore, the understanding of the immunometabolic profile of neutrophils during Leishmania infection could provide new promising targets for host-directed therapies against VL. To our knowledge, this is the first study addressing the bioenergetics profile of L. donovani-infected primary human neutrophils. Transcriptome analysis of L. donovani-infected neutrophils revealed an early significant upregulation of several glycolytic enzymes. Extracellular flux analysis showed that glycolysis and glycolytic capacity were upregulated in L. donovani-infected neutrophils at 6 h post infection. An increased glucose uptake and accumulation of glycolytic end products were further signs for an elevated glycolytic metabolism in L. donovani-infected neutrophils. At the same time point, oxidative phosphorylation provided NADPH for the oxidative burst but did not contribute to ATP production. Inhibition of glycolysis with 2-DG significantly reduced the survival of L. donovani promastigotes in neutrophils and in culture. However, this reduction was due to a direct antileishmanial effect of 2-DG and not a consequence of enhanced antileishmanial activity of neutrophils. To further address the impact of glucose metabolism during the first days of infection in vivo, we treated C57BL/6 mice with 2-DG prior to infection with L. donovani and assessed the parasite load one day and seven days post infection. Our results show, that seven days post-infection the parasite load of 2-DG treated animals was significantly higher than in mock treated animals. This data indicates that glycolysis serves as major energy source for antimicrobial effector functions against L. donovani. Inhibition of glycolysis abrogates important neutrophil effector functions that are necessary the initial control of Leishmania infection.

Closing the Gap: Mechanisms of Epithelial Fusion During Optic Fissure Closure

A key embryonic process that occurs early in ocular development is optic fissure closure (OFC). This fusion process closes the ventral optic fissure and completes the circumferential continuity of the 3-dimensional eye. It is defined by the coming together and fusion of opposing neuroepithelia along the entire proximal-distal axis of the ventral optic cup, involving future neural retina, retinal pigment epithelium (RPE), optic nerve, ciliary body, and iris. Once these have occurred, cells within the fused seam differentiate into components of the functioning visual system. Correct development and progression of OFC, and the continued integrity of the fused margin along this axis, are important for the overall structure of the eye. Failure of OFC results in ocular coloboma-a significant cause of childhood visual impairment that can be associated with several complex ocular phenotypes including microphthalmia and anterior segment dysgenesis. Despite a large number of genes identified, the exact pathways that definitively mediate fusion have not yet been found, reflecting both the biological complexity and genetic heterogeneity of the process. This review will highlight how recent developmental studies have become focused specifically on the epithelial fusion aspects of OFC, applying a range of model organisms (spanning fish, avian, and mammalian species) and utilizing emerging high-resolution live-imaging technologies, transgenic fluorescent models, and unbiased transcriptomic analyses of segmentally-dissected fissure tissue. Key aspects of the fusion process are discussed, including basement membrane dynamics, unique cell behaviors, and the identities and fates of the cells that mediate fusion. These will be set in the context of what is now known, and how these point the way to new avenues of research.

Attenuated Induction of the Unfolded Protein Response in Adult Human Primary Astrocytes in Response to Recurrent Low Glucose

AIMS/HYPOTHESIS

Recurrent hypoglycaemia (RH) is a major side-effect of intensive insulin therapy for people with diabetes. Changes in hypoglycaemia sensing by the brain contribute to the development of impaired counterregulatory responses to and awareness of hypoglycaemia. Little is known about the intrinsic changes in human astrocytes in response to acute and recurrent low glucose (RLG) exposure.

METHODS

Human primary astrocytes (HPA) were exposed to zero, one, three or four bouts of low glucose (0.1 mmol/l) for three hours per day for four days to mimic RH. On the fourth day, DNA and RNA were collected. Differential gene expression and ontology analyses were performed using DESeq2 and GOseq, respectively. DNA methylation was assessed using the Infinium MethylationEPIC BeadChip platform.

RESULTS

24 differentially expressed genes (DEGs) were detected (after correction for multiple comparisons). One bout of low glucose exposure had the largest effect on gene expression. Pathway analyses revealed that endoplasmic-reticulum (ER) stress-related genes such as HSPA5, XBP1, and MANF, involved in the unfolded protein response (UPR), were all significantly increased following low glucose (LG) exposure, which was diminished following RLG. There was little correlation between differentially methylated positions and changes in gene expression yet the number of bouts of LG exposure produced distinct methylation signatures.

CONCLUSIONS/INTERPRETATION

These data suggest that exposure of human astrocytes to transient LG triggers activation of genes involved in the UPR linked to endoplasmic reticulum (ER) stress. Following RLG, the activation of UPR related genes was diminished, suggesting attenuated ER stress. This may be a consequence of a successful metabolic adaptation, as previously reported, that better preserves intracellular energy levels and a reduced necessity for the UPR.

RNA-Seq Transcriptome Analysis Provides Candidate Genes for Resistance to Tomato Leaf Curl New Delhi Virus in Melon

Tomato leaf curl New Delhi virus (ToLCNDV) emerged in the Mediterranean Basin in 2012 as the first DNA bipartite begomovirus (Geminiviridae family), causing severe yield and economic losses in cucurbit crops. A major resistance locus was identified in the wild melon accession WM-7 (Cucumis melo kachri group), but the mechanisms involved in the resistant response remained unknown. In this work, we used RNA-sequencing to identify disease-associated genes that are differentially expressed in the course of ToLCNDV infection and could contribute to resistance. Transcriptomes of the resistant WM-7 genotype and the susceptible cultivar Piñonet Piel de Sapo (PS) (C. melo ibericus group) in ToLCNDV and mock inoculated plants were compared at four time points during infection (0, 3, 6, and 12 days post inoculation). Different gene expression patterns were observed over time in the resistant and susceptible genotypes in comparison to their respective controls. Differentially expressed genes (DEGs) in ToLCNDV-infected plants were classified using gene ontology (GO) terms, and genes of the categories transcription, DNA replication, and helicase activity were downregulated in WM-7 but upregulated in PS, suggesting that reduced activity of these functions reduces ToLCNDV replication and intercellular spread and thereby contributes to resistance. DEGs involved in the jasmonic acid signaling pathway, photosynthesis, RNA silencing, transmembrane, and sugar transporters entail adverse consequences for systemic infection in the resistant genotype, and lead to susceptibility in PS. The expression levels of selected candidate genes were validated by qRT-PCR to corroborate their differential expression upon ToLCNDV infection in resistant and susceptible melon. Furthermore, single nucleotide polymorphism (SNPs) with an effect on structural functionality of DEGs linked to the main QTLs for ToLCNDV resistance have been identified. The obtained results pinpoint cellular functions and candidate genes that are differentially expressed in a resistant and susceptible melon line in response to ToLCNDV, an information of great relevance for breeding ToLCNDV-resistant melon cultivars.

Whole Transcriptome Analysis Provides Insights Into the Molecular Mechanisms of Chlamydospore-Like Cell Formation in Phanerochaete chrysosporium

Phanerochaete chrysosporium is a white rot fungus naturally isolated from hardwoods and widely used in environmental pollution control because it produces extracellular peroxidases. It forms chlamydospores during nitrogen starvation, which naturally occurs in the habitat of P. chrysosporium. Chlamydospores protect fungi against many stresses; the molecular basis underlying chlamydospore formation in basidiomycetes is poorly explored. Chlamydospores in P. chrysosporium have a different cell wall compared with hyphae, as confirmed by cell wall digestion and microscopy. Furthermore, this study investigated the transcriptome of P. chrysosporium in different life stages, including conidium, hypha, and chlamydospore formation, through RNA sequencing. A total of 2215 differentially expressed genes were identified during these processes. The expression patterns of genes involved in several molecular events critical for chlamydospore formation, including starch and sucrose metabolism, phosphatase and kinase, and transcription factors, were determined. This study serves as a basis for further investigating the function of chlamydospore formation in the biotechnologically relevant fungus P. chrysosporium.

Transcriptome Analysis of Ovarian and Uterine Clear Cell Malignancies

Purpose

Ovarian and uterine clear cell carcinomas (CCCs) are rare but associated with poor prognosis. This study explored RNA transcription patterns characteristic of these tumors.

Experimental Design

RNA sequencing (RNA-seq) of 11 ovarian CCCs and five uterine CCCs was performed and compared to publicly available data from high grade serous ovarian cancers (HGSOCs). Ingenuity Pathway Analyses were performed. CIBERSORT analyses estimated relative fractions of 22 immune cell types in each RNA-seq sample. Sequencing data was correlated with PD-L1 immunohistochemical expression.

Results

RNA-seq revealed 1,613 downregulated and 1,212 upregulated genes (corrected p < 0.05, |FC |≥10) in ovarian CCC versus HGSOC. Two subgroups were identified in the ovarian CCC, characterized by ethnicity and expression differences in ARID1A. There were 3,252 differentially expressed genes between PD-L1+/− ovarian CCCs, revealing immune response, cell death, and DNA repair networks, negatively correlated with PD-L1 expression, whereas cellular proliferation networks positively correlated with expression. In clear cell ovarian versus clear cell uterine cancer, 1,607 genes were significantly upregulated, and 109 genes were significantly downregulated (corrected p < 0.05, |FC|≥10). Comparative pathway analysis of late and early stage ovarian CCCs revealed unique metabolic and PTEN pathways, whereas uterine CCCs had unique Wnt/Ca+, estrogen receptor, and CCR5 signaling. CIBERSORT analysis revealed that activated mast cells and regulatory T cell populations were relatively enriched in uterine CCCs. The PD-L1+ ovarian CCCs had enriched resting NK cells and memory B cell populations, while PD-L1− had enriched CD8 T-cells, monocytes, eosinophils, and activated dendritic cells.

Conclusions

Unique transcriptional expression profiles distinguish clear cell uterine and ovarian cancers from each other and from other more common histologic subtypes. These insights may aid in devising novel therapeutics.

Comparative Analysis of Skeletal Muscle Transcriptional Signatures Associated With Aerobic Exercise Capacity or Response to Training in Humans and Rats

Increasing exercise capacity promotes healthy aging and is strongly associated with lower mortality rates. In this study, we analyzed skeletal muscle transcriptomics coupled to exercise performance in humans and rats to dissect the inherent and response components of aerobic exercise capacity. Using rat models selected for intrinsic and acquired aerobic capacity, we determined that the high aerobic capacity muscle transcriptome is associated with pathways for tissue oxygenation and vascularization. Conversely, the low capacity muscle transcriptome indicated immune response and metabolic dysfunction. Low response to training was associated with an inflammatory signature and revealed a potential link to circadian rhythm. Next, we applied bioinformatics tools to predict potential secreted factors (myokines). The predicted secretome profile for exercise capacity highlighted circulatory factors involved in lipid metabolism and the exercise response secretome was associated with extracellular matrix remodelling. Lastly, we utilized human muscle mitochondrial respiration and transcriptomics data to explore molecular mediators of exercise capacity and response across species. Human transcriptome comparison highlighted epigenetic mechanisms linked to exercise capacity and the damage repair for response. Overall, our findings from this cross-species transcriptome analysis of exercise capacity and response establish a foundation for future studies on the mechanisms that link exercise and health.

Transcriptional and Proteomic Choreography Under Phosphorus Deficiency and Re-supply in the N2 Fixing Cyanobacterium Trichodesmium erythraeum

The N₂ fixing cyanobacterium Trichodesmium is a critically important organism in oligotrophic marine ecosystems, supplying “new” nitrogen (N) to the otherwise N-poor tropical and subtropical regions where it occurs. Low concentrations of phosphorus (P) in these regions can constrain Trichodesmium distribution and N₂ fixation rates. Physiological characterization of a single species in a mixed community can be challenging, and ‘omic approaches are increasingly important tools for tracking nutritional physiology in a taxon-specific manner. As such, studies examining the dynamics of gene and protein markers of physiology (e.g., nutrient stress) are critical for the application and interpretation of such ‘omic data in situ. Here we leveraged combined transcriptomics, proteomics, and enzyme activity assays to track the physiological response of Trichodesmium erythraeum IMS101 to P deficiency and subsequent P re-supply over 72 h of sampling. P deficiency resulted in differential gene expression, protein abundance, and enzyme activity that highlighted a synchronous shift in P physiology with increases in the transcripts and corresponding proteins for hydrolyzing organic phosphorus, taking up phosphate with higher affinity, and modulating intracellular P demand. After P deficiency was alleviated, gene expression of these biomarkers was reduced to replete levels within 4 h of P amendment. A number of these gene biomarkers were adjacent to putative pho boxes and their expression patterns were similar to a sphR response regulator. Protein products of the P deficiency biomarkers were slow to decline, with 84% of the original P deficient protein set still significantly differentially expressed after 72 h. Alkaline phosphatase activity tracked with proteins for this enzyme. With the rapid turnover time of transcripts, they appear to be good biomarkers of a P stress phenotype, whereas proteins, with a slower turnover time, may better reflect cellular activities. These results highlight the importance of validating and pairing transcriptome and proteome data that can be applied to physiological studies of key species in situ.

Global Transcriptomic Analysis of the Canine corpus luteum (CL) During the First Half of Diestrus and Changes Induced by in vivo Inhibition of Prostaglandin Synthase 2 (PTGS2/COX2)

The canine luteal phase exhibits several peculiarities compared with other species. In early diestrus, the corpus luteum (CL) is, at least in part, independent of gonadotropins, and prostaglandins (PGs) appear to be among its main regulators. This was also observed with the inhibition in vivo of COX2, when also transcriptional capacity, vascularization and immune-related factors were affected. Here, we aimed to further investigate the potential effects of PGs withdrawal on the CL transcriptome by performing deep RNA sequencing (RNA-Seq). Samples from a previous in vivo study were used; bitches were treated for 5, 10, 20, or 30 days after ovulation with firocoxib (Previcox®), a PTGS2/COX2 inhibitor, or a placebo. Analysis of results was performed with SUSHI (framework from FGCZ) and with pathways and functional networks analyzers. Time-dependent effects were also investigated and used for quality control. More highly represented differentially expressed genes (DEGs, P < 0.01, FDR < 0.1) in the early CL (days 5 and 10) referred to proliferation and immune system, while in the mature CL (days 20 and 30) they were related with steroidogenesis. The absence of genes concomitantly affected by the treatment at all time-points suggested stage-dependency in the observed effects. Little effect was observed on days 5 and 10. Day 20 had the highest number of DEGs (n = 1,741), related with increased immune response. On day 30, DEGs found (n = 552) referred to decreased steroidogenesis and vascularization. Our results suggest the presence of strong compensatory effects in the early CL and multidirectional effects toward gonadotropin-dependency of the CL after COX2 inhibition.

Enriched Differentiation of Human Otic Sensory Progenitor Cells Derived From Induced Pluripotent Stem Cells

Age-related neurosensory deficit of the inner ear is mostly due to a loss of hair cells (HCs). Development of stem cell-based therapy requires a better understanding of factors and signals that drive stem cells into otic sensory progenitor cells (OSPCs) to replace lost HCs. Human induced pluripotent stem cells (hiPSCs) theoretically represent an unlimited supply for the generation of human OSPCs in vitro. In this study, we developed a monolayer-based differentiation system to generate an enriched population of OSPCs via a stepwise differentiation of hiPSCs. Gene and protein expression analyses revealed the efficient induction of a comprehensive panel of otic/placodal and late otic markers over the course of the differentiation. Furthermore, whole transcriptome analysis confirmed a developmental path of OSPC differentiation from hiPSCs. We found that modulation of WNT and transforming growth factor-β (TGF-β) signaling combined with fibroblast growth factor 3 (FGF3) and FGF10 treatment over a 6-day period drives the expression of early otic/placodal markers followed by late otic sensory markers within 13 days, indicative of a differentiation into embryonic-like HCs. In summary, we report a rapid and efficient strategy to generate an enriched population of OSPCs from hiPSCs, thereby establishing the value of this approach for disease modeling and cell-based therapies of the inner ear.

Meta-analysis and transcriptome profiling reveal hub genes for soybean seed storage composition during seed development

Soybean is an important crop providing edible oil and protein source. Soybean oil and protein contents are quantitatively inherited and significantly affected by environmental factors. In this study, meta-analysis was conducted based on soybean physical maps to integrate quantitative trait loci (QTLs) from multiple experiments in different environments. Meta-QTLs for seed oil, fatty acid composition, and protein were identified. Of them, 11 meta-QTLs were located on hot regions for both seed oil and protein. Next, we selected 4 chromosome segment substitution lines with different seed oil and protein contents to characterize their 3 years of phenotype selection in the field. Using strand-specific RNA-sequencing analysis, we profile the time-course transcriptome patterns of soybean seeds at early maturity, middle maturity, and dry seed stages. Pairwise comparison and K-means clustering analysis revealed 7,482 differentially expressed genes and 45 expression patterns clusters. Weighted gene coexpression network analysis uncovered 46 modules of gene expression patterns. The 2 most significant coexpression networks were visualized, and 7 hub genes were identified that were involved in soybean oil and seed storage protein accumulation processes. Our results provided a transcriptome dataset for soybean seed development, and the candidate hub genes represent a foundation for further research.