Influence of estradiol on bovine trophectoderm and uterine gene transcripts around maternal recognition of pregnancy†.

Embryo survival and pregnancy success is increased among animals that exhibit estrus prior to fixed time-artificial insemination, but there are no differences in conceptus survival to d16. The objective of this study was to determine effects of preovulatory estradiol on uterine transcriptomes, select trophectoderm (TE) transcripts, and uterine luminal fluid proteins. Beef cows/heifers were synchronized, artificially inseminated (d0), and grouped into either high (highE2) or low (lowE2) preovulatory estradiol. Uteri were flushed (d16); conceptuses and endometrial biopsies (n = 29) were collected. RNA sequencing was performed on endometrium. Real-time polymerase chain reaction (RT-PCR) was performed on TE (n = 21) RNA to measure relative abundance of IFNT, PTGS2, TM4SF1, C3, FGFR2, and GAPDH. Uterine fluid was analyzed using 2D Liquid Chromatography with tandem mass spectrometry-based Isobaric tags for relative and absolute quantitation (iTRAQ) method. RT-PCR data were analyzed using the MIXED procedure in SAS. There were no differences in messenger RNA (mRNA) abundances in TE, but there were 432 differentially expressed genes (253 downregulated, 179 upregulated) in highE2/conceptus versus lowE2/conceptus groups. There were also 48 differentially expressed proteins (19 upregulated, 29 downregulated); 6 of these were differentially expressed (FDR < 0.10) at the mRNA level. Similar pathways for mRNA and proteins included: calcium signaling, protein kinase A signaling, and corticotropin-releasing hormone signaling. These differences in uterine function may be preparing the conceptus for improved likelihood of survival after d16 among highE2 animals.

Hydrotropism in the primary roots of maize.

Recent studies mainly in Arabidopsis have renewed interest and discussion in some of the key issues in hydrotropism of roots, such as the site of water sensing and the involvement of auxin. We examined hydrotropism in maize (Zea mays) primary roots. We determined the site of water sensing along the root using a nonintrusive method. Kinematic analysis was conducted to investigate spatial root elongation during hydrotropic response. Indole-3-acetic acid (IAA) and other hormones were quantified using LC-MS/MS. The transcriptome was analyzed using RNA sequencing. Main results: The very tip of the root is the most sensitive to the hydrostimulant. Hydrotropic bending involves coordinated adjustment of spatial cell elongation and cell flux. IAA redistribution occurred in maize roots, preceding hydrotropic bending. The redistribution is caused by a reduction of IAA content on the side facing a hydrostimulant, resulting in a higher IAA content on the dry side. Transcriptomic analysis of the elongation zone prior to bending identified IAA response and lignin synthesis/wall cross-linking as some of the key processes occurring during the early stages of hydrotropic response. We conclude that maize roots differ from Arabidopsis in the location of hydrostimulant sensing and the involvement of IAA redistribution.

Root isoflavonoids and hairy root transformation influence key bacterial taxa in the soybean rhizosphere.

Rhizodeposits play a key role in shaping rhizosphere microbial communities. In soybean, isoflavonoids are a key rhizodeposit component that aid in plant defense and enable symbiotic associations with rhizobia. However, it is uncertain if and how they influence rhizosphere microbial communities. Isoflavonoid biosynthesis was silenced via RNA interference of isoflavone synthase in soybean hairy root composite plants. Rhizosphere soil fractions tightly associated with roots were isolated, and PCR amplicons from 16S rRNA gene variable regions V1-V3 and V3-V5 from these fractions were sequenced using 454. The resulting data was resolved using MOTHUR and vegan to identify bacterial taxa and evaluate changes in rhizosphere bacterial communities. The soybean rhizosphere was enriched in Proteobacteria and Bacteroidetes, and had relatively lower levels of Actinobacteria and Acidobacteria compared with bulk soil. Isoflavonoids had a small effect on bacterial community structure, and in particular on the abundance of Xanthomonads and Comamonads. The effect of hairy root transformation on rhizosphere bacterial communities was largely similar to untransformed plant roots with approximately 74% of the bacterial families displaying similar colonization underscoring the suitability of this technique to evaluate the influence of plant roots on rhizosphere bacterial communities. However, hairy root transformation had notable influence on Sphingomonads and Acidobacteria.

Comparative analysis of methicillin-sensitive and resistant Staphylococcus aureus exposed to emodin based on proteomic profiling.

Emodin has a strong antibacterial activity, including methicillin-resistant Staphylococcus aureus (MRSA). However, the mechanism by which emodin induces growth inhibition against MRSA remains unclear. In this study, the isobaric tags for relative and absolute quantitation (iTRAQ) proteomics approach was used to investigate the modes of action of emodin on a MRSA isolate and methicillin-sensitive S. aureus ATCC29213(MSSA). Proteomic analysis showed that expression levels of 145 and 122 proteins were changed significantly in MRSA and MSSA, respectively, after emodin treatment. Comparative analysis of the functions of differentially expressed proteins between the two strains was performed via bioinformatics tools blast2go and STRING database. Proteins related to pyruvate pathway imbalance induction, protein synthesis inhibition, and DNA synthesis suppression were found in both methicillin-sensitive and resistant strains. Moreover, Interference proteins related to membrane damage mechanism were also observed in MRSA. Our findings indicate that emodin is a potential antibacterial agent targeting MRSA via multiple mechanisms.

Code Interpreter for Bioinformatics: Are We There Yet?

The Code Interpreter feature in ChatGPT has the potential to democratize data analysis for non-specialists. As bioinformaticians, we are impressed by its performance in data manipulation and visualization. However, bioinformatics tasks often require execution of third-party packages, access to annotation knowledgebase, and handling large datasets. Code Interpreter's exclusive support for Python, no installation option for additional packages, inability to utilize external resources, and limited storage capacity could pose obstacles to its wide adoption in bioinformatics applications. To address these limitations, we advocated for the necessity of locally deployable, API-based systems for chatbot-aided bioinformatics applications.

Influence of conceptus presence and preovulatory estradiol exposure on uterine gene transcripts and proteins around maternal recognition of pregnancy in beef cattle.

The uterine environment must provide sufficient endocrine conditions and nutrients for pregnancy maintenance and conceptus survival. The objective of this study was to determine the effects of preovulatory estradiol and conceptus presence on uterine transcripts and uterine luminal fluid (ULF) proteins. Beef cows/heifers were synchronized and artificially inseminated (d 0). Uteri were flushed (d 16); conceptuses and endometrial biopsies were collected. Total cellular RNA was extracted from endometrium for RNA sequencing and RT-PCR validation. There were two independent ULF pools made for each of the following groups: highE2/conceptus, highE2/noconceptus, lowE2/conceptus, and lowE2/noconceptus that were analyzed using the 2D LC-MS/MS based iTRAQ method. There were 64 differentially expressed genes (DEGs) and 77 differentially expressed proteins (DEPs) in common among the highE2/conceptus vs highE2/noconceptus and lowE2/conceptus vs lowE2/noconceptus groups. In summary, the interaction between preovulatory estradiol and the conceptus induces the expression of genes, proteins, and pathways necessary for pregnancy.

Subfamily-specific differential contribution of individual monomers and the tether sequence to mouse L1 promoter activity.

BACKGROUND: The internal promoter in L1 5'UTR is critical for autonomous L1 transcription and initiating retrotransposition. Unlike the human genome, which features one contemporarily active subfamily, four subfamilies (A_I, Gf_I and Tf_I/II) have been amplifying in the mouse genome in the last one million years. Moreover, mouse L1 5'UTRs are organized into tandem repeats called monomers, which are separated from ORF1 by a tether domain. In this study, we aim to compare promoter activities across young mouse L1 subfamilies and investigate the contribution of individual monomers and the tether sequence. RESULTS: We observed an inverse relationship between subfamily age and the average number of monomers among evolutionarily young mouse L1 subfamilies. The youngest subgroup (A_I and Tf_I/II) on average carry 3-4 monomers in the 5'UTR. Using a single-vector dual-luciferase reporter assay, we compared promoter activities across six L1 subfamilies (A_I/II, Gf_I and Tf_I/II/III) and established their antisense promoter activities in a mouse embryonic fibroblast cell line and a mouse embryonal carcinoma cell line. Using consensus promoter sequences for three subfamilies (A_I, Gf_I and Tf_I), we dissected the differential roles of individual monomers and the tether domain in L1 promoter activity. We validated that, across multiple subfamilies, the second monomer consistently enhances the overall promoter activity. For individual promoter components, monomer 2 is consistently more active than the corresponding monomer 1 and/or the tether for each subfamily. Importantly, we revealed intricate interactions between monomer 2, monomer 1 and tether domains in a subfamily-specific manner. Furthermore, using three-monomer 5'UTRs, we established a complex nonlinear relationship between the length of the outmost monomer and the overall promoter activity. CONCLUSIONS: The laboratory mouse is an important mammalian model system for human diseases as well as L1 biology. Our study extends previous findings and represents an important step toward a better understanding of the molecular mechanism controlling mouse L1 transcription as well as L1's impact on development and disease.

Early gene response of human brain microvascular endothelial cells to Listeria monocytogenes infection.

The gene expression of human brain microvascular endothelial cells (HBMEC) in response to 4 h of infection by Listeria monocytogenes was analyzed. Four hours after infection, the expression of 456 genes of HBMEC had changed (p < 0.05). We noted that many active genes were involved in the formyl-methionyl-leucyl-phenylalanine pathway in infected HBMEC. In the upregulated genes, mRNA levels of interleukin-8 and interleukin-15 in infected cells increased according to microarray and real-time reverse transcription - PCR analyses. Since both cytokines are regarded as potent chemotactic factors, the results suggest that HBMEC are capable of recruiting cells of innate and adaptive immune responses during early L. monocytogenes infection.

iDEP Web Application for RNA-Seq Data Analysis.

RNA sequencing (RNA-seq) has become a routine method for transcriptomic profiling. We developed a user-friendly web app called iDEP (integrated differential expression and pathway analysis) to help biologists interpret read counts or other types of expression matrices derived from read mapping. With iDEP, users can easily conduct exploratory data analysis, identify differentially expressed genes, and perform pathway analysis. Due to its intuitive user interface and massive annotation database, iDEP is being widely adopted for interactive analysis of RNA-seq data. Using a public dataset on the effect of heat shock on mouse with and without functional Hsf1, we demonstrate how users can prepare data files and conduct in-depth analysis. We also discuss the importance of critical interpretion of results (avoid p-hacking and rationalizing) and validation of significant pathways by using different methods and independent annotation databases.

nCov2019: an R package for studying the COVID-19 coronavirus pandemic.

BACKGROUND: The global spreading of the COVID-19 coronavirus is still a serious public health challenge. Although there are a large number of public resources that provide statistics data, tools for retrospective historical data and convenient visualization are still valuable. To provide convenient access to data and visualization on the pandemic we developed an R package, nCov2019 (https://github.com/YuLab-SMU/nCov2019). METHODS: We collect stable and reliable data of COVID-19 cases from multiple authoritative and up-to-date sources, and aggregate the most recent and historical data for each country or even province. Medical progress information, including global vaccine development and therapeutics candidates, were also collected and can be directly accessed in our package. The nCov2019 package provides an R language interfaces and designed functions for data operation and presentation, a set of interfaces to fetch data subset intuitively, visualization methods, and a dashboard with no extra coding requirement for data exploration and interactive analysis. RESULTS: As of January 14, 2021, the global health crisis is still serious. The number of confirmed cases worldwide has reached 91,268,983. Following the USA, India has reached 10 million confirmed cases. Multiple peaks are observed in many countries. Under the efforts of researchers, 51 vaccines and 54 drugs are under development and 14 of these vaccines are already in the pre-clinical phase. DISCUSSION: The nCov2019 package provides detailed statistics data, visualization functions and the Shiny web application, which allows researchers to keep abreast of the latest epidemic spread overview.

Cbf genes of the Fr-A2 allele are differentially regulated between long-term cold acclimated crown tissue of freeze-resistant and - susceptible, winter wheat mutant lines.

BACKGROUND: In order to identify genes that might confer and maintain freeze resistance of winter wheat, a comparative transcriptome analysis was performed between control and 4 wk cold-acclimated crown tissue of two winter wheat lines that differ in field freeze survival. The lines, generated by azide mutagenesis of the winter wheat cultivar 'Winoka' were designated FR (75% survival) and FS (30% survival). Using two winter lines for this comparative analysis removed the influence of differential expression of the vernalization genes and allowed our study to focus on Cbf genes located within the Fr-A2 allele independent of the effect of the closely mapped Vrn allele. RESULTS: Vernalization genes, (Vrn-A1, B1 and D1), and the transcription factor gene, TaVrt-2, were up-regulated to the same extent in FR and FS lines with cold acclimation thus confirming that azide mutagenesis had not modified the winter habitat of the lines. One category of Cbf genes, (Cbf-2, -A22 and B-22) reflected an increase in level of expression with cold acclimation in both FR and FS lines. Another category of Cbf genes (Cbf-3, -5, -6, -12, -14 and -19) were differentially expressed between cold-acclimated FR and FS lines relative to the non-acclimated controls. Comparison of expression patterns of the two categories of Cbf genes with the expression patterns of a set of ABA-dependent and -independent Cor/Lea genes revealed similar patterns of expression for this sample of Cor/Lea genes with that for Cbf-2 and -22. This pattern of expression was also exhibited by the Vrn genes. CONCLUSION: Some Cor/Lea genes may be co-regulated by the Vrn genes during cold acclimation and the Vrn genes may also control the expression of Cbf-2, -A22 and -B22. The increased freeze survival by the FR line and the increase in expression levels of wheat Cbf genes, Cbf-3, -5, -6, -12, -14 and -19 with cold acclimation in the FR line suggests a possible gain of function mutation resulting in higher levels of expression of these Cbf genes and increased freeze survival.

Evaluating the impact of health awareness events on Google search frequency.

Over two hundred health awareness events take place in the United States in order to educate the public about various diseases. It would be informative and instructive for the organizations to know the impact of these events, although such information could be difficult to measure. We investigated whether 46 selected events attract the public attention by increasing the search frequencies of certain keywords. Internet search data from 2004 to 2017 were downloaded from Google Trend (GT). Three statistical methods including Transfer Function Noise modeling, Wilcoxon Rank Sum test, and Binomial inference were conducted. Our study showed that 10 health awareness events resulted in increased search frequencies in the event months, and 28 events did not, with the rest being classified as unclear.

Genome-wide analysis of antisense transcription with Affymetrix exon array.

BACKGROUND: A large number of natural antisense transcripts have been identified in human and mouse genomes. Study of their potential functions clearly requires cost-efficient method for expression analysis. RESULTS: Here we show that Affymetrix Exon arrays, which were designed to detect conventional transcripts in the sense orientation, can be used to monitor antisense expression across all exonic loci in mammalian genomes. Through modification of the cDNA synthesis protocol, we labeled single-strand cDNA in the reverse orientation as in the standard protocol, thus enabling the detection of antisense transcripts using the same array. Applying this technique to human Jurkat cells, we identified antisense transcription at 2,088 exonic loci of 1,516 UniGene clusters. Many of these antisense transcripts were not observed previously and some were validated by orientation-specific RT-PCR. CONCLUSION: Our results suggest that with a modified protocol Affymetrix human, mouse and rat Exon arrays can be used as a routine method for genome-wide analysis of antisense transcription in these genomes.

Identifying nonspecific SAGE tags by context of gene expression.

Many serial analysis of gene expression (SAGE) tags can be matched to multiple genes, leading to difficulty in SAGE data interpretation and analysis. As only a subset of genes in the human genome are transcribed in a certain type of tissue/cell, we used microarray expression data from different tissue types to define contexts of gene expression and to annotate SAGE tags collected from the same or similar tissue sources. To predict the original transcript contributing a nonspecific SAGE tag collected from a particular tissue, we ranked the corresponding genes by their expression levels determined by microarray. We developed a tissue-specific SAGE tag annotation database based on microarray data collected from 73 normal human tissues and 18 cancer tissues and cell lines. The database can be queried online at: http://www.basic.northwestern.edu/SAGE/. The accuracy of this database was confirmed by experimental data.

A large quantity of novel human antisense transcripts detected by LongSAGE.

MOTIVATION: Taking advantage of the high sensitivity and specificity of LongSAGE tag for transcript detection and genome mapping, we analyzed the 632 813 unique human LongSAGE tags deposited in public databases to identify novel human antisense transcripts. RESULTS: Our study identified 45 321 tags that match the antisense strand of 9804 known mRNA sequences, 6606 of which contain antisense ESTs and 3198 are mapped only by SAGE tags. Quantitative analysis showed that the detected antisense transcripts are present at levels lower than their counterpart sense transcripts. Experimental results confirmed the presence of antisense transcripts detected by the antisense tags. We also constructed an antisense tag database that can be used to identify the antisense SAGE tags originated from the antisense strand of known mRNA sequences included in the RefSeq database. CONCLUSIONS: Our study highlights the benefits of exploring SAGE data for comprehensive identification of human antisense transcripts and demonstrates the prevalence of antisense transcripts in the human genome.

ATOH1 Promotes Leptomeningeal Dissemination and Metastasis of Sonic Hedgehog Subgroup Medulloblastomas.

Medulloblastoma arising from the cerebellum is the most common pediatric brain malignancy, with leptomeningeal metastases often present at diagnosis and recurrence associated with poor clinical outcome. In this study, we used mouse medulloblastoma models to explore the relationship of tumor pathophysiology and dysregulated expression of the NOTCH pathway transcription factor ATOH1, which is present in aggressive medulloblastoma subtypes driven by aberrant Sonic Hedgehog/Patched (SHH/PTCH) signaling. In experiments with conditional ATOH1 mouse mutants crossed to Ptch1+/- mice, which develop SHH-driven medulloblastoma, animals with Atoh1 transgene expression developed highly penetrant medulloblastoma at a young age with extensive leptomeningeal disease and metastasis to the spinal cord and brain, resembling xenografts of human SHH medulloblastoma. Metastatic tumors retained abnormal SHH signaling like tumor xenografts. Conversely, ATOH1 expression was detected consistently in recurrent and metastatic SHH medulloblastoma. Chromatin immunoprecipitation sequencing and gene expression profiling identified candidate ATOH1 targets in tumor cells involved in development and tumorigenesis. Among these targets specific to metastatic tumors, there was an enrichment in those implicated in extracellular matrix remodeling activity, cytoskeletal network and interaction with microenvironment, indicating a shift in transcriptomic and epigenomic landscapes during metastasis. Treatment with bone morphogenetic protein or SHH pathway inhibitors decreased tumor cell proliferation and suppressed metastatic tumor growth, respectively. Our work reveals a dynamic ATOH1-driven molecular cascade underlying medulloblastoma metastasis that offers possible therapeutic opportunities. Cancer Res; 77(14); 3766-77. ©2017 AACR.

SAGE detects microRNA precursors.

BACKGROUND: MicroRNAs (miRNAs) have been shown to play important roles in regulating gene expression. Since miRNAs are often evolutionarily conserved and their precursors can be folded into stem-loop hairpins, many miRNAs have been predicted. Yet experimental confirmation is difficult since miRNA expression is often specific to particular tissues and developmental stages. RESULTS: Analysis of 29 human and 230 mouse longSAGE libraries revealed the expression of 22 known and 10 predicted mammalian miRNAs. Most were detected in embryonic tissues. Four SAGE tags detected in human embryonic stem cells specifically match a cluster of four human miRNAs (mir-302a, b, c&d) known to be expressed in embryonic stem cells. LongSAGE data also suggest the existence of a mouse homolog of human and rat mir-493. CONCLUSION: The observation that some orphan longSAGE tags uniquely match miRNA precursors provides information about the expression of some known and predicted miRNAs.

Identification and characterization of lin-28 homolog B (LIN28B) in human hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide. Several studies have identified signature gene sets that may be useful as potential diagnostic tools by global microarray analysis. Here we report the cloning and characterization of a novel gene, lin-28 homolog B (LIN28B), which is overexpressed in hepatocellular carcinoma. The heterochronic gene lin-28 is a key regulator of developmental timing in the nematode Caenorhabditis elegans. Similar with lin-28 proteins, LIN28B conserves a cold shock domain and a pair of CCHC zinc finger domains. Phylogenetic analysis suggests that they might arise as a result of duplication from an ancestral gene. Overexpression of LIN28B was noted in most HCC cell lines and clinical samples. By western blot analysis using a polyclonal antibody against LIN28B, a short LIN28B isoform was also identified in non-tumor liver tissue and fetal liver. Although predominantly localized in the cytoplasm, we found that LIN28B protein shows cell cycle-dependent nuclear translocation in Huh7 cells. Induced expression of exogenous LIN28B in a tet-off cell line promoted cancer cell proliferation. Interestingly, the segment of the unusually long 3'UTR of LIN28B contains complementary sites to let-7 microRNA of mammals. And our studies provided indirect evidence that LIN28B is a possibly natural target for let-7 mediated regulation. These findings strongly implicate a critical role of LIN28B during development and tumorigenesis and suggest a possible novel mechanism.

Two distinct gene expression signatures in pediatric acute lymphoblastic leukemia with MLL rearrangements.

Acute lymphoblastic leukemia (ALL) with 11q23 translocations is usually associated with MLL gene rearrangement, but little is known about its leukemogenesis. We analyzed the gene expression profiles of pediatric ALL samples according to their translocations. Using oligonucleotide microarray analysis, we identified distinct expression profiles for 23 ALL samples with 11q23 translocations, including t(4;11) (n = 15), t(11;19) (n = 6), and t(5;11) (n = 2), compared with 9 ALL samples with other translocations, including t(12;21) (n = 6) and t(1;19) (n = 3). Gene expression scores of FLT3, MeisI, and CD44 for samples with MLL rearrangements were particularly high compared with those for other ALL samples. Statistical analysis of the gene expression profiles for the 21 ALL samples with MLL rearrangements at diagnosis revealed two subgroups that exclusively correlated with prognosis but not with any other clinico-pathological factor. The transcription factors CBF2 and CDP were highly expressed in the poor and good prognosis subgroups, respectively. In addition, their downstream target genes were differentially expressed. These findings provide new insights into the biological mechanisms of leukemogenesis and prognosis for pediatric ALL with MLL rearrangements.

2.45 GHz radiofrequency fields alter gene expression in cultured human cells.

The biological effect of radiofrequency (RF) fields remains controversial. We address this issue by examining whether RF fields can cause changes in gene expression. We used the pulsed RF fields at a frequency of 2.45 GHz that is commonly used in telecommunication to expose cultured human HL-60 cells. We used the serial analysis of gene expression (SAGE) method to measure the RF effect on gene expression at the genome level. We observed that 221 genes altered their expression after a 2-h exposure. The number of affected genes increased to 759 after a 6-h exposure. Functional classification of the affected genes reveals that apoptosis-related genes were among the upregulated ones and the cell cycle genes among the downregulated ones. We observed no significant increase in the expression of heat shock genes. These results indicate that the RF fields at 2.45 GHz can alter gene expression in cultured human cells through non-thermal mechanism.