Les organes centraux du tractus urogénital sont le testicule et l’ovaire, qui assurent la production de gamètes et d’hormones, et donc la fertilité de l’individu. Ces deux organes, parfaitement distincts et complémentaires, ont pour origine une gonade bipotentielle qui s’engagera vers une trajectoire de différenciation masculine ou féminine au cours de la vie fœtale. Les deux gonades vont par la suite subir plusieurs phases de différenciation et de développement de leurs populations cellulaires, afin d’acquérir leurs fonctions propres qui leur permettront d’assumer leur rôle à l’âge adulte. Depuis plus d’une quinzaine d’années, le concept de syndrome de dysgénésie testiculaire fait état d’un lien entre l’exposition du fœtus à des composés environnementaux et des anomalies du tractus urogénital. Bien que sujette à de vifs débats au sein de la communauté scientifique, cette hypothèse a attiré l’attention de la recherche sur les conséquences de l’exposition des mères aux xénobiotiques sur l’enfant à naître. La différenciation et le développement des gonades fœtales sont gouvernés par des programmes d’expression spécifiques de chaque sexe, dont de nombreuses zones d’ombres subsistent, notamment concernant la fraction non-codante exprimée par le génome humain. La première partie de cette thèse de doctorat présente, pour la première fois, le paysage transcriptionnel contrôlant ces processus complexes entre la 6ième et 17ième semaine de développement chez l’Homme. Grâce à l’avènement des technologies de transcriptomique, il est désormais possible d’identifier et d’observer l’expression des gènes de manière sensible et sans a priori. Le RNA-seq m’a donc permis de décrire de manière exhaustive la dynamique d’expression des gènes, pendant les stades précoces de la différenciation sexuelle, jusqu’aux phénomènes plus tardifs conduisant aux linéages des différentes populations cellulaires spécifiques du testicule et de l’ovaire. Dans une deuxième partie, mon travail de recherche s’est attaché à étudier l’impact de deux perturbateurs endocriniens suspectés, l’ibuprofène et le chlordécone, sur le programme d’expression du testicule fœtal humain. L’utilisation du RNA-seq m’a permis de définir et de comparer la signature toxicogénomique de chaque molécule afin de contribuer à la compréhension de leur mécanisme d’action et d’identifier les populations cellulaires affectées. Enfin, face à l’essor des technologies ultra-haut-débit dans les sciences de la vie, y compris dans les domaines de la reproduction, j’ai activement participé au déploiement d’une nouvelle version du Reprogenomic Viewer dans la dernière partie de ma thèse (http://rgv.genouest.org). Cet outil nternet a pour vocation de centraliser et de rendre accessibles les données de séquençage accumulées au sein de la communauté de la reproduction via des outils de visualisation intuitifs. / Fetal life is a crucial period for sexual reproduction when bipotential gonads differentiate into either a testis or an ovary. Gaining insights into the complex molecular events underlying this process is central to a better understanding of disorders of sexual development. The present work intends to improve the knowledge on molecular pathways at play during gonad development in humans using RNA-sequencing. This project particularly seeks to identify early transcriptional events that may play critical role in the regulatory network driving human sexual differentiation. To address this issue, we defined the transcriptional landscape of fetal human gonads by sequencing total RNA extracted from testes and ovaries between 6 and 17 gestational weeks. The resulting paired-end reads were mapped on the human genome and then assembled into transcripts using the Tuxedo suite. We next defined a high-confidence set of transcripts showing differential expression across samples. Clusters of co-expressed genes were subjected to functional analysis. The analysis of this massive RNA-seq dataset has led to a high-confidence set of 35,194 assembled transcripts; among which 32,391 known and novel isoforms coding genes (mRNAs), 1,209 to long non-coding (lnc) RNAs and 318 to novel unannotated transcripts/genes (NUTs). The dynamic of transcriptional landscape occurring during human fetal gonads development has been described and new genes and interesting candidates, including new genes, have been highlighted as potential key genes governing this biological process. The second interest of this work was the study of the impact of two endocrine disruptors, ibuprofene and chlordecone, on human fetal testis using RNA-seq. The transcriptional alteration induced by these compound in the gonad allowed a deeper understanding of their mechanisms of action of endocrine disruption. The last part of this work was the development of a new version of the ReproGenomics Viewer (http://rgv.genouest.org), a web tool dedicated to the integration and accumulation of sequencing data from studies performed in the field of reproduction.

L’insecte phytophage Riptortus pedestris, appartenant au sous-ordre des Hétéroptères, est un ravageur notoire de cultures agricoles en Asie du sud-est qui se nourrit préférentiellement de plants de soja. Cette punaise est associée à une bactérie symbiotique du genre Burkholderia nommée Burkholderia insecticola, localisée dans une région spécifique de l’intestin de l’insecte appelée la région M4. Cette région M4, organisée en cryptes, constitue l’organe symbiotique dans lequel le symbiote prolifère de manière extracellulaire. Cette interaction favorise la croissance et le développement de la punaise. Récemment, il a été montré que Riptortus produit des peptides antimicrobiens au sein des cryptes, appelés “crypt-specific cysteine-rich peptides” ou peptides CCR pour lesquels le symbiote est particulièrement résistant. Il a été proposé que les peptides antimicrobiens de l’hôte,incluant les peptides CCR, participent à la colonisation spécifique de l’organe symbiotique par B. insecticola. Dans ce travail, une approche Tn-seq a été utilisée pour identifier les gènes bactériens impliqués dans la résistance aux peptides antimicrobiens et dans la symbiose. Dans un premier temps, la robustesse de la méthode Tn-seq a été évaluée en identifiant le génome essentiel de B. insecticola. Puis dans un second temps, les facteurs bactériens impliqués dans la résistance aux peptides antimicrobiens ont été caractérisés via une approche gènes-candidats et l’approche Tn-seq. Dans une dernière partie, une expérience de Tn-seq in vivo a permis d’évaluer l’ampleur du goulot d’étranglement sur la population symbiotique lors de l’infection de l’organe symbiotique et d’identifier les facteurs symbiotiques impliqués dans la colonisation de R. pedestris. / The phytophagous insect Riptortus pedestris, belonging to the Heteroptera suborder, is a notorious crop pest in South-Eastern Asia which feeds preferentially on soybean plants. This bean bug is associated with a bacterial symbiont, a specific Burkholderia species named Burkholderia insecticola, located in the M4 region of the insect’s midgut. This M4 region is organized in crypts and constitutes the symbiotic organ where the symbiont proliferates extracellularly. This interaction promotes the growth and the development of the bean bug. Recently, it was demonstrated that Riptortus produces antimicrobial peptides in the midgut crypts called crypt-specific cysteine-rich peptides (CCR) for which the bacterial symbiont demonstrates a high resistance profile. It was proposed that host antimicrobial peptides, including the CCR peptides, contribute to the specific colonization of the symbiotic organ by B. insecticola. In this work, a Tn-seq approach was used to find bacterial fitness genes involved in antimicrobial peptide resistance and symbiosis. First, the robustness of the Tn-seq method was assessed by identifying the essential genome of B. insecticola. Second, the bacterial factors for antimicrobial peptide resistance were characterized, based on both a candidate-gene and the Tn-seq approach. Finally, a Tn-seq in vivo experiment was performed to reveal the infection bottleneck effect on the symbiotic population and to identify the bacterial symbiosis factors for the colonization of R. pedestris.

Le phytopathogène Agrobacterium tumefaciens est l'agent causal de la maladie appelée galle du collet, et est capable d'infecter plus de 90 familles de plantes dicotylédones. Cette ∝-protéobactérie appartient à la famille Rhizobiaceae. A. tumefaciens est un complexe de différentes espèces regroupées en 10 génomovars (G1 à G8 et G13). A. tumefaciens C58 appartient au groupe du G8. Son génome est constitué de 4 réplicons : 1 chromosome circulaire, 1 chromosome linéaire et des 2 plasmides dispensables : pAt (pour A.tumefaciens) et pTi (pour Tumor inducing, qui est requis pour la virulence). Pour explorer de nouveaux aspects du mode de vie d’A. tumefaciens, et en particulier l'interaction entre la bactérie et sa plante hôte, deux approches différentes ont été utilisées pour identifier, caractériser et analyser les gènes qui pourraient jouer un rôle dans l'adaptation des bactéries à la tumeur. Une expérience de l'évolution par des passages en série de trois souches différentes de l'agent pathogène sur la plante hôte Solanum lycopersicum a été effectuée afin de clarifier la dynamique évolutive du génome au cours de l'infection. Parallèlement, une étude de différents transcriptomes (in planta et in vitro) a été réalisée et étudiée pour élucider des gènes bactériens candidats impliqués dans l'interaction de la bactérie avec la plante et divers composés produits dans la tumeur. Ce travail tente de donner une vue plus générale du processus d'adaptation de la bactérie à la niche écologique qui est la tumeur. / Agrobacterium tumefaciens is the causal agent of the plant disease called crowngall, and it’s able to infect more than 90 families of dicotyledonous plants. It is an α-Proteobacterium and belongs to the Rhizobiaceae family. A. tumefaciens is a complex of different species grouped in 10 genomovars (G1 to G8, and G13). A. tumefaciens C58 belongs to the G8 group. Its genome consists in 4 replicons: 1 chromosome circular, 1 chromosome linear and 2 dispensable plasmids: pAt (for A. tumefaciens) and pTi (for Tumor inducing), which is required for virulence. To explore new aspects of the A. tumefaciens lifestyle, and in particular the interaction between the bacteria and its plant host, two different approaches have been used to identify, characterize and analyze genes that could play a role in the adaptation of the bacteria to tumor lifestyle. An evolution experiment by serial passages of three different strains of thepathogen on the host plant Solanum lycopersicum has been carried out to clarify the evolutionary dynamics of the genome during the course of infection. In parallel, a study of different transcriptomes (in planta and in vitro) was performed and studied to elucidate bacterial candidate genes involved in the interaction of the bacteria with the plant and various compounds produced in the tumor. This work attempts to give a more general view of the process of adaptation of the bacteria to the ecological niche that is the tumor.

Dans l’histoire évolutive, deux gènes paralogues sont issus d’un évènement de duplication de leur ancêtre commun. Les gènes paralogues sont caractérisés par des duplications globales de génome (WGD) ou à petite échelle (SSD) et par leur datation. Les WGDs ont lieu à deux reprises à la base de la lignée des vertébrés. Les évènements de SSD ont lieu à plusieurs moments pouvant être plus récents, plus anciens ou contemporain de la période des évènements de WGD. La rétention des paralogues dans le génome, associée à une divergence de l’expression spatiale est une contribution importante pour l’augmentation de la complexité de l’organisme au cours de l’évolution. Certaines études ont montré que les duplications anciennes seraient plus associées aux maladies. L’objectif de la première partie de la thèse est de créer une ressource sur les paralogues en collectant et en analysant différentes annotations. Nous avons construit une ressource robuste de paralogues humains à partir de listes publiées mais aussi à partir d’annotations externes. L’exploration de différentes annotations nous a permis d’identifier une identité de séquence élevée entre gènes paralogues pouvant biaiser la mesure d’expression des gènes et diminuer leur expression. L’objectif de la seconde partie, est d’explorer l’expression spatiale et la co- expression des paralogues au sein du cerveau humain, à partir des données RNA-seq du consortium GTEx. Les données d’expression GTEx de 13 tissus cérébraux, nous ont permis de montrer que la datation récente mais aussi que le type SSD contribuaient à une expression plus tissu-spécifique. Nous avons utilisé l’analyse de la co-expression (WGCNA) afin de regrouper les paralogues possédant une expression similaire au travers des tissus et nous avons pu suggérer une co-expression des SSD récents. Nos études sur les maladies ont montré que les SSD récents accumulaient des mutations associées à des maladies cérébrales. Finalement, nous avons trouvé que la co-expression des paralogues et leur tissu-spécificité au travers des régions cérébrales pouvaient enrichir nos connaissances sur les gènes associés à des maladies cérébrales. / In evolution history, two paralogous genes originate from the duplication event of a common ancestor gene. Paralogous genes are characterized by whole genome (WGD) or small-scale (SSD) duplications and their duplication date. The WGDs happened twice in the early vertebrate lineage. SSD events take place at any moment in evolutionary history and can be younger, older or dating to the same period than WGD events. Retention of paralogs in the genome associated with divergence of spatial expression is an important contributor to the increase of organism complexity through evolution. Different studies found that old duplications are more associated with diseases. The objective of the first part of the thesis is to create a resource on paralogs by collecting and analyzing annotations. We built a robust resource of human paralogs from published lists of paralogous genes and also from external annotations. Annotation exploration allowed us to identify a high sequence identity between paralogous genes impacting the gene expression measurement from RNA-seq data and decreasing the gene expression. The objective of the second part is to explore spatial expression and co-expression of paralogs in the human brain, from the GTEx consortium RNA-seq data. The GTEx expression data of 13 brain tissues allowed us to show that duplication youth and SSD type contributed to a more tissue-specific expression. We used co-expression analyses (WGCNA) to group paralogs with similar expression across tissues and we suggested the co-expression of younger SSDs. Our disease studies showed the younger SSD accumulation of mutations associated with brain diseases. We finally found that paralog co-expression and their tissue-specificity across brain regions could enrich information of known brain disease-associated genes.

No description available.

Les technologies récentes de séquençage de l’ARN de cellules uniques (scRNA-seq, pour single cell RNA-seq) ont permis de quantifier le niveau d’expression des gènes au niveau de la cellules, alors que les technologies standards de séquençage de l’ARN (RNA-seq, ou bulk RNA-seq) ne permettaient de quantifier que l’expression moyenne des gènes dans un échantillon de cellules. Cette résolution supérieure a permis des avancées majeures dans le domaine de la recherche biomédicale, mais a également posé de nouveaux défis, notamment computationnels. Les données qui découlent des technologies de scRNA-seq sont en effet complexes et plus bruitées que les données de bulk RNA-seq. En outre, les technologies sont nombreuses et leur nombre explose, nécessitant chacune un prétraitement plus ou moins différent. De plus en plus de méthodes sont ainsi développées, mais il n’existe pas encore de norme établie (gold standard) pour le prétraitement et l’analyse de ces données. Le laboratoire du Dr. Mader a récemment fait l’acquisition de la technologie Drop-Seq (une technologie haut débit de scRNA-seq), nécessitant une expertise nouvelle pour le traitement des données qui en découlent. Dans ce mémoire, différentes étapes du prétraitement des données issues de la technologie Drop-Seq sont donc passées en revue et le fonctionnement de certains outils dédiés à cet effet est étudié, permettant d’établir des lignes directrices pour de futures expériences au sein du laboratoire du Dr. Mader. Cette étude est menée sur les premiers jeux de données générés avec la technologie Drop-Seq du laboratoire, issus de lignées cellulaires du cancer du sein. Les méthodes d’analyses, moins spécifiques à la technologie, ne sont pas étudiées dans ce mémoire, mais une analyse exploratoire des jeux de données du laboratoire pose les bases pour une analyse plus poussée. / Recent single cell RNA sequencing technologies (scRNA-seq) have enabled the quantification of gene expression levels at the cellular level, while standard RNA sequencing technologies (RNA-seq, or bulk RNA-seq) have only been able to quantify the average gene expression in a sample of cells. This higher resolution has allowed major advances in biomedical research, but has also raised new challenges, in particular computational ones. The data derived from scRNA-seq technologies are indeed complex and noisier than bulk RNA-seq data. In addition, the number of scRNA-seq technologies is exploding, each of them requiring a rather different pre-processing. More and more methods are thus being developed, but there is still no gold standard for the preprocessing and analysis of these data. Dr. Mader’s laboratory has recently invested in the Drop-Seq technology (a high-throughput scRNAseq technology), requiring new expertise for the processing of the resulting data. In this thesis, different steps for the pre-processing of Drop-Seq data are reviewed and the behavior of some of the dedicated tools are studied, allowing to establish guidelines for future experiments in Dr. Mader’s laboratory. This study is conducted on the first data sets generated with the Drop-Seq technology of the laboratory, derived from breast cancer cell lines. Analytical methods, less specific to the technology, are not investigated in this thesis, but an exploratory analysis of the lab’s datasets lays the foundation for further analysis.

With the rapid progress of next-generation sequencing (NGS) technologies, new tools and methods have emerged to investigate the transcriptomics of various organisms. RNA sequencing (RNA-seq) employs NGS to evaluate the presence and abundance of RNA transcripts in biological samples. This technique offers a comprehensive snapshot of the RNA dynamics within cells. With the ability to profile the entire transcriptome of organisms rapidly and accurately, RNA-seq has become the state-of-the-art method for transcriptome profiling, surpassing the traditional microarray approach. Single-cell RNA sequencing (scRNA-seq) was introduced in 2009 to profile the single-cell gene expression in highly heterogeneous samples such as brain tissue and tumors. The advancement of scRNA-seq technologies enables the in-depth transcriptomic study in each cell subtype. When selecting an scRNA-seq method, researchers must weigh the trade-off between profiling more single cells versus obtaining more comprehensive transcripts per cell, while considering the overall costs. The throughput of full-length scRNA-seq methods is usually lower, as each single cell needs to be processed separately to produce scRNA-seq libraries. However, full-length methods enable the researchers to investigate the splicing variants and allele-specific expression. Non-full-length methods only capture the 3' or 5' ends of transcripts, which limits their application in isoform detection, but as cells are pooled after barcoding for cDNA synthesis, the throughput is 2–3 orders of magnitude higher than full-length methods. We developed a droplet-based platform for full-length single-cell RNA-seq, which enabled the efficient recovery of full-length mRNA from individual cells in a high-throughput manner. The developed platform can process ~8,000 single cells within 2 days and detect ~20% more genes compared to Drop-seq. Besides scRNA-seq technology development, we also applied a low-input RNA-seq method to study the transcriptomics in different biological samples. When handling precious biological samples, a low-input method is necessary to profile the transcriptome of homogeneous cell populations. We first studied the epigenomic and transcriptomic regulations in colorectal cancer (CRC) using MOWChIP-seq, a low-input high-throughput method, in conjunction with our low-input RNA-seq approach. Fusobacterium nucleatum (Fnn) is closely related to the progression of cancers like CRC and pancreatic cancer. However, the molecular mechanisms of how Fnn adjusts the tumor microenvironment (TME) and leads to poor clinical outcomes are still unclear. In this in-vitro study, we characterized how hypoxia, an important TME ignored by previous research, facilitates Fnn infection of CRC and corresponding alterations of global epigenome and transcriptome. We infer that hypoxia has similar effects as Fnn infection alone on the CRC cells. The Fnn infection under hypoxia further boosts the proliferation and progression of CRC. We then applied our low-input RNA-seq method to study brain neuroscience and immunology. Psychedelics like DOI show promising clinical efficacy in patients with psychiatric conditions. Although psychedelics exhibit rapid antidepression action and long-lasting effectiveness compared to conventional treatment, their acute psychotic symptoms and potential for drug abuse discourage their application in clinical practice. In this case, it is important to comprehend the molecular mechanisms responsible for psychedelics' clinical efficacy. This understanding can pave the way for the development of improved treatments that do not rely on psychedelics. After profiling the transcriptome of mouse brain samples exposed to psychedelics with different post-exposure times, we concluded that the psychedelic-induced transcriptomic variations are more transient than epigenomic changes. In the second brain neuroscience project, we first applied 3-color FACS sorting to differentiate four neuron and non-neuron subtypes in human postmortem prefrontal cortex tissues. Then we profiled the gene expression of the four subtypes and validated the FACS sorting by examining the expression of marker genes. Differentially expressed genes between each subtype and the others were extracted and proceeded to gene ontology analysis. We identified unique altered biological pathways related to each subtype. The immunology research focuses on revealing the difference between low-grade inflammation and monocyte exhaustion, as well as the unique biological pathways they regulate. Therefore, we profiled the transcriptome of bone marrow-derived monocytes stimulated by PBS control, a low- or high-dose LPS. In addition to wild-type mice, we also included TRAM-deficient and IRAK-M-deficient mice. We concluded that low-dose LPS specifically regulates the TRAM-dependent pathway of TLR4 signaling, and high-dose LPS exclusively upregulates exhaustion markers by impacting metabolic and proliferative pathways. / Doctor of Philosophy / Transcriptomics is the comprehensive study of RNA transcripts derived from an organism's genome. RNA plays a vital role in maintaining the fundamental functions of cells and organisms. In eukaryotes, the genetic information stored in the DNA of cells is transferred to messenger RNA (mRNA) molecules through a process called transcription. These mRNA molecules serve as a bridge between DNA and proteins, as they carry the instructions encoded in genes to ribosomes for protein synthesis. Studying mRNA transcripts reveals various cellular mechanisms and their impact on overall organism function, gene regulation, and disease pathways. With the aid of next-generation sequencing, various RNA-seq approaches have been developed to study mRNA transcripts quantitatively in the past decades. To better understand the gene expression regulations in biological samples, we first applied bulk RNA-seq to profile the transcriptome of various samples under different conditions. Our in-house bulk RNA-seq protocol has been proven to be both high-performance and cost-effective compared to commercial kits. To better understand cellular diversity and uncover rare cell types in heterogeneous biological samples, we developed a droplet-based scRNA-seq platform that can recover full-length mRNA transcripts in a high throughput manner. It can profile the transcriptome of thousands of single cells within two days. It combines the advantages of the droplet-based scRNA-seq method (high throughput) and the well plate-based scRNA-seq method (full-length mRNA recovery).

No description available.

Cell death shapes multicellular organism development and sustains tissue and organ homeostasis. Over the past decade we have begun to understand the breadth of physiological and biochemical diversity in cell death and clearance pathways, which play vital roles in organismal development and heath. While apoptosis and necrosis have been studied extensively across many model systems and contexts, the discovery of non-apoptotic paradigms of cell death and their roles in disease has greatly expanded the field. Collectively called Regulated Cell Death (RCD), these death pathways are regulated in a tissue and context-dependent manner (e.g. disease state). This dissertation is a culmination of multiple projects investigating cell death and clearance events spanning the ovary and the brain of the model organism, Drosophila melanogaster. We undertook the first multi-modal, high-throughput survey, involving single-cell RNA-seq, TRAP-seq, and proteomics, to compare two different archetypes of germline death in the fly egg chamber - apoptosis and phagoptosis. Our analysis identified several important candidates and pathways that are either unique to or shared between the germline death modalities and affecting oogenesis upon their disruption. We also observed that V-ATPases, proton pumps required for germline phagoptosis, are differentially localized throughout oogenesis, and we identified the specific subunits upregulated in phagoptosis. Furthermore, we identified a novel exon splicing event in the ‘a’ subunit isoform of V-ATPases that may facilitate its sub-cellular localization. Using a novel image analysis method involving image segmentation and spatial statistical inference, we determined that circulating immune cells agglomerate at specific niches within the abdomen, in response to egg chamber degeneration resulting from physiological stress of protein-deprivation. We then turned our focus to phagocytosis in the fly brain, which is essential for pruning synapses and for the removal of dying neurons and misfolded proteins. Disruptions to glial phagocytosis results in a range of age-dependent neurodegenerative phenotypes, primarily exemplified by vacuolization of brain tissue. Using a pre-trained deep-learning model to perform image segmentation and 3D reconstruction of vacuoles, we characterized the severity of neurodegeneration in brains lacking the phagocytic receptor Draper in glia and further demonstrated that this phenotype is attenuated by knockdown of the NF-κB transcription factor Relish in flies lacking glial Draper. Collectively, the methods and results described herein will have applications beyond the Drosophila model and the field of cell death, with important implications in understanding fertility and the underpinnings of cognitive disorders.

Female adults of many mosquito species possess distinct physiological features adapting to blood feeding for successful reproduction. The disease pathogens that are transmitted by mosquitoes have evolved to take advantages of the indispensable blood feedings to complete their transmission cycles and to survive attacks from the mosquito's innate immune system. Normal egg development and mosquito immunity are tightly controlled by tissue- and stage-specific gene expression and coordinated by many signal molecules in the mosquito. Understanding gene regulation affecting mosquito reproduction and malaria parasites infection is of paramount importance for developing novel malaria control strategies. A growing body of evidence indicates that microRNAs (miRNAs) are involved in egg maturation and immune reactions against invading pathogens in mosquitoes. However, the molecular mechanisms by which specific miRNAs selectively modulate reproduction and the survival of pathogens are largely unknown. The miRNA-induced gene-silencing pathway in mosquitoes was mostly extrapolated from the studies of flies. To explore the dynamics of miRNAs in reproduction, I used small RNAs sequencing to monitor miRNAs expression and their association with Argonaute 1 (Ago1) and Argonaute 2 (Ago2) in the malaria mosquito Anopheles gambiae (An. gambiae) during the 72-h period immediately after blood feeding. I found the abundance and Ago loading of most of the mature miRNAs were relatively stable after blood ingestion. However, miRNAs of the miR-309/286/2944 cluster were considerably upregulated after blood feeding. I confirmed that miR-309 is essential for normal egg development by depletion of endogenous miR-309 with a specific antagomir. In addition, my results showed that the Ago association of some miRNAs was not proportional to their cellular abundance implying additional regulation at miRNA integration. To investigate the functional roles of miRNAs and define context-dependent miRNA-mRNA interactions during the reproductive process, I have applied an innovative experimental approach to study miRNA-mRNA interactome. CLEAR (covalent ligation of endogenous Argonaute-bound RNAs)-CLIP can generate miRNA-mRNA chimeras from UV-irradiation stabilized Ago-miRNA-mRNA complex. My results have defined tens of thousands of miRNA-mRNA interactions in mosquitoes, including novel targets for mosquito-specific miRNAs. Verification of the predicted interactions using mRNA-seq, ribosome-profiling, and luciferase reporter assay revealed a reliable miRNA-mRNA interaction network. Based on the detected interactions, I refined the paring rules for mosquito miRNAs and illustrated the dynamic pairing between different regions of miRNAs with their targets in vivo. The miRNA-mRNA interactions were compared using this approach at multiple time points before and after blood feeding. Importantly, this study showed that the interactions were dynamic and enriched in genes that are involved in metabolisms, supporting the proposed functions of miRNAs in coordinating the gene regulation in mosquito reproduction. Plasmodium falciparum (P. falciparum) is a major human malaria parasite. To understand the functions of miRNAs in the mosquito resistance to Plasmodium infection, we analyzed the miRNA-mRNA interactions after female mosquitoes taking a P. falciparum-infected blood meal or an uninfected blood meal. Comparison of the interactions revealed enhanced miRNA-mRNA interactions after P. falciparum infection involving a group of immunity-related genes. In summary, this study has provided a systematic view and significantly advanced our understanding of the miRNA functions in mosquito reproduction and P. falciparum infection. / PHD / Female mosquito is able to transmit lots of disease to the human when it bites for blood. The blood meal provides necessary nutrient for mosquito reproduction and spread the pathogens such as malaria and Zika at the same time. Thus understanding the molecular mechanism behind this process would be greatly helpful to develop novel vector control strategy. Here, we found a distinct class of RNAs contributing to the regulation of mosquito blood meal and parasite infection. We used a novel biochemical method to decoding the special role of these kinds of RNAs in these processes. We found them regulating mosquito metabolism and immunity. This study significantly deepened our knowledge about the process of mosquito reproduction and transmitting diseases.