Eimeria tenella is a eukaryotic, intracellular parasite that, along with six other Eimeria species, causes coccidiosis in chickens. This disease can result in weight loss or even death and is estimated to cause 2 billion euros of damages to the chicken industry each year. While much is known of the life cycle of E. tenella in the chicken, less is known about molecular mechanisms of infection and the chicken immune response. In this study, we produced a pipeline for dual RNA-sequencing analysis of a mixed chicken and E. tenella dataset.  We then carried out an analysis on an in vitro infection of the chicken macrophage HD-11 cell line.  This was followed by a differential expression analysis across six time points, 2, 4, 12, 24, 48, and 72 hours post-infection, in order to elucidate these mechanisms. The results showed clear patterns of expression for the chicken immune genes, with strong down-regulation of genes across the immune system at 24 hours and a repetition of early patterns at 72 hours, indicating that reinfection by a second generation of parasite cells was occurring. Several genes that may have important roles in the immune reaction of the chicken were identified, such as MRC2, ITGB3 and ITGA9, along with genes with known roles, such as TLR15. The expression of surface antigen genes in E. tenella was also examined, showing a clear upregulation in the late stages of merogony, suggesting important roles for merozoites. Finally, a co-expression analysis was carried out, showing considerable co-expression among the two organisms.  One of the gene co-expression networks identified appeared to be enriched with both infection specific genes from E. tenella and chicken immune genes. These results, along with the pipeline, will be used in further studies on E. tenella infections and bring us closer to the eventual goal of a vaccine for coccidiosis.

Since the discovery of transfer RNAs (tRNAs) as decoders of the genetic code, life science has transformed. Particularly, as soon as the importance of tRNAs in protein synthesis has been established, researchers recognized that the functionality of tRNAs in cellular regulation exceeds beyond this paradigm. A strong impetus for these discoveries came from advances in large-scale RNA sequencing (RNA-seq) and increasingly sophisticated algorithms. Sequencing tRNAs is challenging both experimentally and in terms of the subsequent computational analysis. In RNA-seq data analysis, mapping tRNA reads to a reference genome is an error-prone task. This is in particular true, as chemical modifications introduce systematic reverse transcription errors while at the same time the genomic loci are only approximately identical due to the post-transcriptional maturation of tRNAs. Additionally, their multi-copy nature complicates the precise read assignment to its true genomic origin. In the course of the thesis a computational workflow was established to enable accurate mapping of tRNA reads. The developed method removes most of the mapping artifacts introduced by simpler mapping schemes, as demonstrated by using both simulated and human RNA-seq data. Subsequently, the resulting mapping profiles can be used for reliable identification of specific chemical tRNA modifications with a false discovery rate of only 2%. For that purpose, computational analysis methods were developed that facilitates the sensitive detection and even classification of most tRNA modifications based on their mapping profiles. This comprised both untreated RNA-seq data of various species, as well as treated data of Bacillus subtilis that has been designed to display modifications in a specific read-out in the mapping profile. The discussion focuses on sources of artifacts that complicate the profiling of tRNA modifications and strategies to overcome them. Exemplary studies on the modification pattern of different human tissues and the developmental stages of Dictyostelium discoideum were carried out. These suggested regulatory functions of tRNA modifications in development and during cell differentiation. The main experimental difficulties of tRNA sequencing are caused by extensive, stable secondary structures and the presence of chemical modifications. Current RNA-seq methods do not sample the entire tRNA pool, lose short tRNA fragments, or they lack specificity for tRNAs. Within this thesis, the benchmark and improvement of LOTTE-seq, a method for specific selection of tRNAs for high-throughput sequencing, exhibited that the method solves the experimental challenges and avoids the disadvantages of previous tRNA-seq protocols. Applying the accurate tRNA mapping strategy to LOTTE-seq and other tRNA-specific RNA- seq methods demonstrated that the content of mature tRNAs is highest in LOTTE-seq data, ranging from 90% in Spinacia oleracea to 100% in D. discoideum. Additionally, the thesis addressed the fact that tRNAs are multi-copy genes that undergo concerted evolution which keeps sequences of paralogous genes effectively identical. Therefore, it is impossible to distinguish orthologs from paralogs by sequence similarity alone. Synteny, the maintenance of relative genomic positions, is helpful to disambiguate evolutionary relationships in this situation. During this thesis a workflow was computed for synteny-based orthology identification of tRNA genes. The workflow is based on the use of pre-computed genome-wide multiple sequence alignment blocks as anchors to establish syntenic conservation of sequence intervals. Syntenic clusters of concertedly evolving genes of different tRNA families are then subdivided and processed by cograph editing to recover their duplication histories. A useful outcome of this study is that it highlights the technical problems and difficulties associated with an accurate analysis of the evolution of multi-copy genes. To showcase the method, evolution of tRNAs in primates and fruit flies were reconstructed. In the last decade, a number of reports have described novel aspects of tRNAs in terms of the diversity of their genes. For example, nuclear-encoded mitochondrial-derived tRNAs (nm-tRNAs) have been reported whose presence provokes intriguing questions about their functionality. Within this thesis an annotation strategy was developed that led to the identification of 335 and 43 novel nm-tRNAs in human and mouse, respectively. Interestingly, downstream analyses showed that the localization of several nm-tRNAs in introns and the over-representation of conserved RNA-binding sites of proteins involved in splicing suggest a potential regulatory function of intronic nm-tRNAs in splicing.

RBMS3 belongs to the family of c-myc gene single-strand binding proteins (MSSPs) that play important roles in transcriptional regulation. Here, we show that RBMS3 functions as a tumor promoter in triple-negative breast cancer (TNBC), a highly aggressive BC subtype. Analysis of RBMS3 expression shows that RBMS3 is upregulated at both mRNA and protein levels in TNBC cells. Functionally, overexpression of RBMS3 increases cell migration, invasion and cancer stem cell (CSC) behaviors. Moreover, RBMS3 induces expression of epithelial-mesenchymal transition (EMT) and CSC markers. Conversely, loss of RBMS3 in TNBC BT549 cells inhibits cell proliferation, migration and mesenchymal phenotype. Correlation analysis shows RBMS3 is associated with TGF-β signaling. Mechanistically, RBMS3 interacts with Smad2, Smad3 and Smad4 mRNA and regulates the stability of these transcripts. Importantly, RBMS3 prevents TGF-β-induced cytostasis and apoptosis in premalignant cancer cells. Moreover, RBMS3 inversely correlates with expression of ESRPs, epithelial-specific splicing regulatory proteins that regulate morphogenesis-associated alternative splicing events. ESRPs appear to suppress EMT through distinct mechanisms: ESRP1 restricted cell migration, whereas ESRP2 prevented cell growth. RBMS3 significantly facilitates the EMT process when ESRPs are lost. Collectively, the studies within this dissertation identify RBMS3 as a positive regulator of EMT and breast cancer progression by regulating the TGF-β signaling pathway.

Plant yield is an agronomic trait dependent on the transport of photosynthate from mature source leaves to sink tissues. Manipulating phloem transport may lead to increased yield, however in a previous study, Arabidopsis thaliana overexpressing sucrose transporter AtSUC2 in the phloem resulted in stunted growth and an apparent P-deficiency. In the course of further characterizing the phenotype and identifying the causative mutation, this research included 1) reverse genetics to test genes hypothesized to modulate carbon-phosphate interactions; 2) whole genome sequencing to identify all T-DNA insertions in plants displaying the phenotype; 3) genetic crosses and segregation analysis to isolate the causative mutation; and 4) transcriptomics to capture gene-expression profiles in plants displaying the phenotype. These phenotypes were traced to a T-DNA insertion located on chromosome 4. Transcriptomics by RNA-Seq and data analysis through bioinformatics pipelines suggest disruptions in metabolic and transport pathways that include phosphate, but do not support a direct role of well-established phosphate acquisition mechanisms. Gene At1G78690 is immediately downstream of the T-DNA insertion site and shows modestly increased expression relative to wild type plants. At1G78690 encodes O-acyl transferase, which is involved in processing N-acylphosphotidyl ethanolamine (NAPE) to N-acyl ethanolamine (NAE). Exogenous NAE application causes stunted growth in specific conditions. From the experiments described herein, At1G78690 emerges as the strongest candidate for causing the observed phenotypes.

Transcription factors activated by exogenous or endogenous stimuli alter gene expression with major effects on chromatin accessibility and the epigenome. This thesis investigates that impact of environmental chemical and hormonal exposure on liver chromatin accessibility in a mouse liver model. Exposure to the constitutive androstane receptor (CAR)-specific agonist ligand 1,4-bis-[2-(3,5-dichloropyridyloxy)]benzene (TCPOBOP) induces localized changes in chromatin accessibility at several thousand DNase hypersensitive sites (DHS). Activating histone marks, associated with enhancers and promoters, were induced by TCPOBOP and were highly enriched at opening DHS. Opening DHS were highly enriched for CAR binding sites and nuclear receptor direct repeat-4 motifs. These DHS were also enriched for the CAR heterodimeric partner RXRA, binding by CEBPA and CEBPB, and motifs for other liver-specific factors. Thus, TCPOBOP alters the enhancer landscape through changes in histone marks and by mechanisms linked to induced CAR binding. In other studies, the impact of pituitary growth hormone (GH) secretion patterns on chromatin accessibility changes associated with sex-biased liver gene expression was examined. In adult male liver, the transcription factor STAT5 is directly activated by each successive plasma GH pulse. In female liver, STAT5 is persistently activated by the near-continuous stimulation by plasma GH. A majority of the ~4,000 GH-regulated, sex-biased DHS have chromatin marks characteristic of enhancers and were enriched for proximity to sex-biased gene promoters. Chromatin accessibility is thus a key feature of sex-differential gene expression. Two major classes of male-biased DHS were identified: dynamic male-biased DHS, almost all bound by STAT5, which undergo repeated cycles of chromatin opening and closing induced by each GH pulse; and static male-biased DHS, whose accessibility is unaffected GH/STAT5 pulses and whose sex bias results from these chromatin sites being more closed in female liver. Sites with STAT5 binding showed greater chromatin opening, many of which also contain the STAT5 motif. Finally, the effect of a single GH pulse on hypophysectomized male mouse liver was investigated to identify DHS responsive to the male, pulsatile-GH, secretion pattern. These studies demonstrate that widespread epigenetic changes associated with target gene expression are induced by xenobiotics and hormones regulating liver gene expression. / 2022-01-31T00:00:00Z

Many omics-level studies have been undertaken on Aiptasia, however, our understanding of the genes and processes associated with symbiosis regulation and maintenance is still limited. To gain deeper insights into the molecular processes underlying this association, we investigated this relationship using multipronged approaches combining next generation sequencing with metabolomics and immunohistochemistry. We identified 731 high-confident symbiosis-associated genes using meta-analysis. Coupled with metabolomic profiling, we exposed that symbiont-derived carbon enables host recycling of ammonium into nonessential amino acids, which may serve as a regulatory mechanism to control symbiont growth through a carbon-dependent negative feedback of nitrogen availability to the symbiont. We then characterized two symbiosis-associated ammonium transporters (AMTs). Both of the proteins exhibit gastrodermis-specific localization in symbiotic anemones. Their tissuespecific localization consistent with the higher ammonium assimilation rate in gastrodermis of symbiotic Aiptasia as shown by 15N labeling and nanoscale secondary ion mass spectrometry (NanoSIMS). Inspired by the tissue-specific localization of AMTs, we investigated spatial expression of genes in Aiptasia. Our results suggested that symbiosis with Symbiodiniaceae is the main driver for transcriptional changes in Aiptasia. We focused on the phagosome-associated genes and identified several key factors involved in phagocytosis and the formation of symbiosome. Our study provided the first insights into the tissue specific complexity of gene expression in Aiptasia. To investigate symbiosis-induced response in symbiont and to find further evidence for the hypotheses generated from our host-focused analyses, we explored the growth and gene expression changes of Symbiodiniaceae in response to the limitations of three essential nutrients: nitrogen, phosphate, and iron, respectively. Comparisons of the expression patterns of in hospite Symbiodiniaceae to these nutrient limiting conditions showed a strong and significant correlation of gene expression profiles to the nitrogen-limited culture condition. This confirmed the nitrogen-limited growing condition of Symbiodiniaceae in hospite, and further supported our hypothesis that the host limits the availability of nitrogen, possibly to regulate symbiont cell density. In summary, we investigated different molecular aspects of symbiosis from both the host’s and symbiont’s perspective. This dissertation provides novel insights into the function of nitrogen, and the potential underlying molecular mechanisms, in the metabolic interactions between Aiptasia and Symbiodiniaceae.

The goal of genetic-based cancer research is often to identify which genes behave differently in cancerous and healthy tissue. This difference in behavior, referred to as differential expression, may lead researchers to more targeted preventative care and treatment. One way to measure the expression of genes is though a process called RNA-Seq, that takes physical tissue samples and maps gene products and fragments in the sample back to the gene that created it, resulting in a large read-count matrix with genes in the rows and a column for each sample. The read-counts for tumor and normal samples are then compared in a process called differential expression analysis. However, normalization of these read-counts is a necessary pre-processing step, in order to account for differences in the read-count values due to non-expression related variables. It is common in recent RNA-Seq normalization methods to also account for gene-level covariates, namely gene length in base pairs and GC-content, the proportion of bases in the gene that are Guanine and Cytosine. Here a colorectal cancer RNA-Seq read-count data set comprised of 30,220 genes and 378 samples is examined. Two of the normalization methods that account for gene length and GC-content, CQN and EDASeq, are extended to account for protein coding status as a third gene-level covariate. The binary nature of protein coding status results in unique computation issues. The results of using the normalized read counts from CQN, EDASeq, and four new normalization methods are used for differential expression analysis via the nonparametric Wilcoxon Rank-Sum Test as well as the lme4 pipeline that produces per-gene models based on a negative binomial distribution. The resulting differential expression results are compared for two genes of interest in colorectal cancer, APC and CTNNB1, both of the WNT signaling pathway.

La transcriptómica estudia el nivel de expresión de los genes en distintas condiciones experimentales para tratar de identificar los genes asociados a un fenotipo dado así como las relaciones de regulación entre distintos genes. Los datos ómicos se caracterizan por contener información de miles de variables en una muestra con pocas observaciones. Las tecnologías de alto rendimiento más comunes para medir el nivel de expresión de miles de genes simultáneamente son los microarrays y, más recientemente, la secuenciación de RNA (RNA-seq). Este trabajo de tesis versará sobre la evaluación, adaptación y desarrollo de modelos estadísticos para el análisis de datos de expresión génica, tanto si ha sido estimada mediante microarrays o bien con RNA-seq. El estudio se abordará con herramientas univariantes y multivariantes, así como con métodos tanto univariantes como multivariantes. / Tarazona Campos, S. (2014). Statistical methods for transcriptomics: From microarrays to RNA-seq [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/48485 / TESIS / Premios Extraordinarios de tesis doctorales

No description available.

Soybean crops are vulnerable to a wide range of pathogens that reduce yield and cause extensive losses worldwide. In the United States, the soilborne pathogens Pythium spp., causing soybean damping-off, and Fusarium virguliforme, causing sudden death syndrome (SDS) of soybean, have been among the top diseases that most reduced soybean yields. This study demonstrated that biological control using native fungal antagonists could be a powerful tool to integrate with current management strategies for more efficient control of Pythium damping-off and SDS in soybean. Trichoderma spp. and Clonostachys rosea demonstrated the ability to mycoparasitize and antagonize the pathogens using different mechanisms and exhibited a protective effect on soybean in field conditions. The development of an efficient biological control program for disease management relies on a deep understanding of the BCA-pathogen interaction’s biology. This research also uncovered the molecular mechanisms involved in the F. virguliforme-T. afroharzianum interaction by using a dual RNAseq approach. Significant changes in both fungal organisms’ transcriptomes were discovered at different stages in their interaction. The results provided here can contribute to the future implementation of effective biological control programs for soybean. The benefits may also extend to other crops.