I analysed data from Native Elongating Transcript Sequencing (NET-seq), a method for sequencing RNAs that are still in the process of transcription (“nascent” RNAs). NET-seq captures nascent RNAs by targeting the enzyme that transcribes the genes into RNA, known as RNA Polymerase II (Pol II). From the nascent RNAs, one can verify whether the introns contained within have been spliced out. The splicing efficiency of each intron is then estimated using a proxy metric called the “splicing ratio” (SR). The higher the SR, the faster the intron is presumed to be spliced. The data can also be used to map the locations of Pol IIs and thus to infer the relative speed with which different gene regions were transcribed (although this may be confounded by Pol II phosphorylation patterns). This is important, as there is evidence for links between splicing and transcriptional pausing. However, these inferences can be marred by technical biases related to the nucleotide content of the sequences. I developed a simulation method to test our biological conclusions against such biases. I also designed a "peak calling" algorithm to determine which putative instances of Pol II pausing were the most reliable. This methodological work was disseminated through a blog post (
https://imm.medicina.ulisboa.pt/news/the-peaks-and-valleys-of-the-nascent-transcriptome-in-drosophila-embryos(si apre in una nuova finestra))
a seminar to students at the Faculty of Science of the University of Lisbon, and two practicals on an introductory bioinformatics course organized by Egypt Scholars, directed at students in Egypt.
I proceeded with a detailed characterisation of co-transcriptional splicing in the fruit fly, published as a co-first author paper in the RNA Journal (rna.078933.121v1) and presented at two national and three international conferences. SR varied drastically between introns, and correlated with properties of the intron in unexpected ways. Moreover, Pol II tended to pause at different locations depending on SR. I used Bayesian modelling to explore different hypotheses for mechanisms underlying these patterns. I concluded that the data could only be accounted for by a model where the same intron can stochastically switch between different modes of splicing kinetics.
Next, I checked whether the frequency or evolutionary conservation of splicing-related sequence signals depended on SR. I failed to uncover any significant patterns. This could be because the data was insufficient for such data-hungry analyses, or because variation in SR is either not functionally relevant or not controlled through sequence signals.
In addition, I used three types of interventions to improve statistical thinking skills among biomedical researchers. Firstly, I designed an eight-week introductory statistics course. The course emphasized conceptual understanding, hands-on practice on real data and group work. I taught this course at the iMM in 2021, training a total of ca. 50 early career researchers. The course was repeated in the spring of 2022, as an online Arabic-language version, delivered in collaboration with Egypt Scholars. Both iterations of the course received overwhelmingly positive feedback.
Secondly, in the June of 2022, I organized an international summer school on applying modelling techniques to biological data. The summer school, funded by a Horizon 2020 grant, was attended by researchers from the iMM in Portugal, the Max-Delbrück Zentrum in Germany, the Weizmann Institute of Science in Israel, and the University of Oxford in the UK. The participants took part in five days of hands-on workshops, delivered by an international group of instructors.
Thirdly, I worked individually with researchers to help them better understand their data. This included the supervision of 2 PhD students, 1 Master’s student and one intern, as well as aiding several other researchers. This work has led to one co-first author publication (10.3390/biomedicines10020199) with at least three other manuscripts in preparation.