Neofunctionalization driven by positive selection led to the retention of the loqs2 gene encoding an Aedes specific dsRNA binding protein.

BACKGROUND: Mosquito borne viruses, such as dengue, Zika, yellow fever and Chikungunya, cause millions of infections every year. These viruses are mostly transmitted by two urban-adapted mosquito species, Aedes aegypti and Aedes albopictus. Although mechanistic understanding remains largely unknown, Aedes mosquitoes may have unique adaptations that lower the impact of viral infection. Recently, we reported the identification of an Aedes specific double-stranded RNA binding protein (dsRBP), named Loqs2, that is involved in the control of infection by dengue and Zika viruses in mosquitoes. Preliminary analyses suggested that the loqs2 gene is a paralog of loquacious (loqs) and r2d2, two co-factors of the RNA interference (RNAi) pathway, a major antiviral mechanism in insects. RESULTS: Here we analyzed the origin and evolution of loqs2. Our data suggest that loqs2 originated from two independent duplications of the first double-stranded RNA binding domain of loqs that occurred before the origin of the Aedes Stegomyia subgenus, around 31 million years ago. We show that the loqs2 gene is evolving under relaxed purifying selection at a faster pace than loqs, with evidence of neofunctionalization driven by positive selection. Accordingly, we observed that Loqs2 is localized mainly in the nucleus, different from R2D2 and both isoforms of Loqs that are cytoplasmic. In contrast to r2d2 and loqs, loqs2 expression is stage- and tissue-specific, restricted mostly to reproductive tissues in adult Ae. aegypti and Ae. albopictus. Transgenic mosquitoes engineered to express loqs2 ubiquitously undergo developmental arrest at larval stages that correlates with massive dysregulation of gene expression without major effects on microRNAs or other endogenous small RNAs, classically associated with RNA interference. CONCLUSIONS: Our results uncover the peculiar origin and neofunctionalization of loqs2 driven by positive selection. This study shows an example of unique adaptations in Aedes mosquitoes that could ultimately help explain their effectiveness as virus vectors.

STING regulates metabolic reprogramming in macrophages via HIF-1α during Brucella infection.

Macrophages metabolic reprogramming in response to microbial insults is a major determinant of pathogen growth or containment. Here, we reveal a distinct mechanism by which stimulator of interferon genes (STING), a cytosolic sensor that regulates innate immune responses, contributes to an inflammatory M1-like macrophage profile upon Brucella abortus infection. This metabolic reprogramming is induced by STING-dependent stabilization of hypoxia-inducible factor-1 alpha (HIF-1α), a global regulator of cellular metabolism and innate immune cell functions. HIF-1α stabilization reduces oxidative phosphorylation and increases glycolysis during infection with B. abortus and, likewise, enhances nitric oxide production, inflammasome activation and IL-1ß release in infected macrophages. Furthermore, the induction of this inflammatory profile participates in the control of bacterial replication since absence of HIF-1α renders mice more susceptible to B. abortus infection. Mechanistically, activation of STING by B. abortus infection drives the production of mitochondrial reactive oxygen species (mROS) that ultimately influences HIF-1α stabilization. Moreover, STING increases the intracellular succinate concentration in infected macrophages, and succinate pretreatment induces HIF-1α stabilization and IL-1ß release independently of its cognate receptor GPR91. Collectively, these data demonstrate a pivotal mechanism in the immunometabolic regulation of macrophages during B. abortus infection that is orchestrated by STING via HIF-1α pathway and highlight the metabolic reprogramming of macrophages as a potential treatment strategy for bacterial infections.

Genomic insights into the c-di-GMP signaling and biofilm development in the saprophytic spirochete Leptospira biflexa.

C-di-GMP is a bacterial second messenger with central role in biofilm formation. Spirochete bacteria from Leptospira genus present a wide diversity, with species of medical importance and environmental species, named as saprophytic. Leptospira form biofilms in the rat's reservoir kidneys and in the environment. Here, we performed genomic analyses to identify enzymatic and effector c-di-GMP proteins in the saprophytic biofilm-forming species Leptospira biflexa serovar Patoc. We identified 40 proteins through local alignments. Amongst them, 16 proteins are potentially functional diguanylate cyclases, phosphodiesterases, or hybrid proteins. We also identified nine effectors, including PilZ proteins. Enrichment analyses suggested that c-di-GMP interacts with cAMP signaling system, CsrA system, and flagella assembly regulation during biofilm development of L. biflexa. Finally, we identified eight proteins in the pathogen Leptospira interrogans serovar Copenhageni that share high similarity with L. biflexa c-di-GMP-related proteins. This work revealed proteins related to c-di-GMP turnover and cellular response in Leptospira and their potential roles during biofilm development.

Phylogenomics and gene selection in Aspergillus welwitschiae: Possible implications in the pathogenicity in Agave sisalana.

Aspergillus welwitschiae causes bole rot disease in sisal (Agave sisalana and related species) which affects the production of natural fibers in Brazil, the main worldwide producer of sisal fibers. This fungus is a saprotroph with a broad host range. Previous research established A. welwitschiae as the only causative agent of bole rot in the field, but little is known about the evolution of this species and its strains. In this work, we performed a comparative genomics analysis of 40 Aspergillus strains. We show the conflicting molecular identity of this species, with one sisal-infecting strain sharing its last common ancestor with Aspergillus niger, having diverged only 833 thousand years ago. Furthermore, our analysis of positive selection reveals sites under selection in genes coding for siderophore transporters, Sodium­calcium exchangers, and Phosphatidylethanolamine-binding proteins (PEBPs). Herein, we discuss the possible impacts of these gene functions on the pathogenicity in sisal.

State of the Art of the Molecular Biology of the Interaction between Cocoa and Witches' Broom Disease: A Systematic Review.

Significant scientific advances to elucidate the Moniliophthora perniciosa pathosystem have been achieved in recent years, but the molecular biology of this pathogen-host interaction is still a field with many unanswered questions. In order to present insights at the molecular level, we present the first systematic review on the theme. All told, 1118 studies were extracted from public databases. Of these, 109 were eligible for the review, based on the inclusion and exclusion criteria. The results indicated that understanding the transition from the biotrophic-necrotrophic phase of the fungus is crucial for control of the disease. Proteins with strong biotechnological potential or that can be targets for pathosystem intervention were identified, but studies regarding possible applications are still limited. The studies identified revealed important genes in the M. perniciosa-host interaction and efficient molecular markers in the search for genetic variability and sources of resistance, with Theobroma cacao being the most common host. An arsenal of effectors already identified and not explored in the pathosystem were highlighted. This systematic review contributes to the understanding of the pathosystem at the molecular level, offering new insights and proposing different paths for the development of new strategies to control witches' broom disease.

BASIDIN as a New Protein Effector of the Phytopathogen Causing Witche's Broom Disease in Cocoa.

The fungus Moniliophthora perniciosa secretes protein effectors that manipulate the physiology of the host plant, but few effectors of this fungus have had their functions confirmed. We performed functional characterization of a promising candidate effector of M. perniciosa. The inoculation of rBASIDIN at 4 µmol L-1 in the mesophyll of leaflets of Solanum lycopersicum caused symptoms of shriveling within 6 h without the presence of necrosis. However, when sprayed on the plant at a concentration of 11 µmol L-1, it caused wilting symptoms only 2 h after application, followed by necrosis and cell death at 48 h. rBASIDIN applied to Theobroma cacao leaves at the same concentration caused milder symptoms. rBASIDIN caused hydrogen peroxide production in leaf tissue, damaging the leaf membrane and negatively affecting the photosynthetic rate of Solanum lycopersicum plants. Phylogenetic analysis indicated that BASIDIN has orthologs in other phytopathogenic basidiomycetes. Analysis of the transcripts revealed that BASIDIN and its orthologs are expressed in different fungal species, suggesting that this protein is differentially regulated in these basidiomycetes. Therefore, the results of applying BASIDIN allow the inference that it is an effector of the fungus M. perniciosa, with a strong potential to interfere in the defense system of the host plant.

The pan-genome of the emerging multidrug-resistant pathogen Corynebacterium striatum.

Corynebacterium striatum, a common constituent of the human skin microbiome, is now considered an emerging multidrug-resistant pathogen of immunocompromised and chronically ill patients. However, little is known about the molecular mechanisms in the transition from colonization to the multidrug-resistant (MDR) invasive phenotype in clinical isolates. This study performed a comprehensive pan-genomic analysis of C. striatum, including isolates from "normal skin microbiome" and from MDR infections, to gain insights into genetic factors contributing to pathogenicity and multidrug resistance in this species. For this, three novel genome sequences were obtained from clinical isolates of C. striatum of patients from Brazil, and other 24 complete or draft C. striatum genomes were retrieved from GenBank, including the ATCC6940 isolate from the Human Microbiome Project. Analysis of C. striatum strains demonstrated the presence of an open pan-genome (α = 0.852803) containing 3816 gene families, including 15 antimicrobial resistance (AMR) genes and 32 putative virulence factors. The core and accessory genomes included 1297 and 1307 genes, respectively. The identified AMR genes are primarily associated with resistance to aminoglycosides and tetracyclines. Of these, 66.6% are present in genomic islands, and four AMR genes, including aac(6')-ib7, are located in a class 1-integron. In conclusion, our data indicated that C. striatum possesses genomic characteristics favorable to the invasive phenotype, with high genomic plasticity, a robust genetic arsenal for iron acquisition, and important virulence determinants and AMR genes present in mobile genetic elements.

In vitro studies provide insight into effects of Dicer-2 helicase mutations in Drosophila melanogaster.

In vitro, Drosophila melanogaster Dicer-2 (Dcr-2) uses its helicase domain to initiate processing of dsRNA with blunt (BLT) termini, and its Platform•PAZ domain to initiate processing of dsRNA with 3' overhangs (ovrs). To understand the relationship of these in vitro observations to roles of Dcr-2 in vivo, we compared in vitro effects of two helicase mutations to their impact on production of endogenous and viral siRNAs in flies. Consistent with the importance of the helicase domain in processing BLT dsRNA, both point mutations eliminated processing of BLT, but not 3'ovr, dsRNA in vitro. However, the mutations had different effects in vivo. A point mutation in the Walker A motif of the Hel1 subdomain, G31R, largely eliminated production of siRNAs in vivo, while F225G, located in the Hel2 subdomain, showed reduced levels of endogenous siRNAs, but did not significantly affect virus-derived siRNAs. In vitro assays monitoring dsRNA cleavage, dsRNA binding, ATP hydrolysis, and binding of the accessory factor Loquacious-PD provided insight into the different effects of the mutations on processing of different sources of dsRNA in flies. Our in vitro studies suggest effects of the mutations in vivo relate to their effects on ATPase activity, dsRNA binding, and interactions with Loquacious-PD. Our studies emphasize the importance of future studies to characterize dsRNA termini as they exist in Drosophila and other animals.

A single unidirectional piRNA cluster similar to the flamenco locus is the major source of EVE-derived transcription and small RNAs in Aedes aegypti mosquitoes.

Endogenous viral elements (EVEs) are found in many eukaryotic genomes. Despite considerable knowledge about genomic elements such as transposons (TEs) and retroviruses, we still lack information about nonretroviral EVEs. Aedes aegypti mosquitoes have a highly repetitive genome that is covered with EVEs. Here, we identified 129 nonretroviral EVEs in the AaegL5 version of the A. aegypti genome. These EVEs were significantly associated with TEs and preferentially located in repeat-rich clusters within intergenic regions. Genome-wide transcriptome analysis showed that most EVEs generated transcripts although only around 1.4% were sense RNAs. The majority of EVE transcription was antisense and correlated with the generation of EVE-derived small RNAs. A single genomic cluster of EVEs located in a 143 kb repetitive region in chromosome 2 contributed with 42% of antisense transcription and 45% of small RNAs derived from viral elements. This region was enriched for TE-EVE hybrids organized in the same coding strand. These generated a single long antisense transcript that correlated with the generation of phased primary PIWI-interacting RNAs (piRNAs). The putative promoter of this region had a conserved binding site for the transcription factor Cubitus interruptus, a key regulator of the flamenco locus in Drosophila melanogaster Here, we have identified a single unidirectional piRNA cluster in the A. aegypti genome that is the major source of EVE transcription fueling the generation of antisense small RNAs in mosquitoes. We propose that this region is a flamenco-like locus in A. aegypti due to its relatedness to the major unidirectional piRNA cluster in Drosophila melanogaster.

Genome-wide identification of miRNAs and target regulatory network in the invasive ectoparasitic mite Varroa destructor.

Varroa destructor is an ectoparasite mite that attacks bees leading to colony disorders worldwide. microRNAs (miRNAs) are key molecules used by eukaryotes to post-transcriptional control of gene expression. Nevertheless, still lack information aboutV. destructor miRNAs and its regulatory networks. Here, we used an integrative strategy to characterize the miRNAs in the V. destructor mite. We identified 310 precursors that give rise to 500 mature miRNAs, which 257 are likely mite-specific elements. miRNAs showed canonical length ranging between 18 and 25 nucleotides and 5' uracil preference. Top 10 elements concentrated over 80% of total miRNA expression, with bantam alone representing ~50%. We also detected non-templated bases in precursor-derived small RNAs, indicative of miRNA post-transcriptional regulatory mechanisms. Finally, we note that conserved miRNAs control similar processes in different organisms, suggesting a conservative role. Altogether, our findings contribute to the better understanding of the mite biology that can assist future studies on varroosis control.

Identification of Glycycometus malaysiensis (for the first time in Brazil), Blomia tropicalis and Dermatophagoides pteronyssinus through multiplex PCR.

Blomia tropicalis and Dermatophagoides pteronyssinus play an important role in triggering allergy. Glycycometus malaysiensis causes IgE reaction in sensitive people, but is rarely reported in domestic dust, because it is morphologically similar to B. tropicalis making the identification of these species difficult. The identification of mites is mostly based on morphology, a time-consuming and ambiguous approach. Herein, we describe a multiplex polymerase chain reaction (mPCR) assay based on ribosomal DNA capable to identify mixed cultures of B. tropicalis, D. pteronyssinus and G. malaysiensis, and/or to identify these species from environmental dust. For this, the internal transcribed spacer 2 (ITS2) regions, flanked by partial sequences of the 5.8S and 28S genes, were PCR-amplified, cloned and sequenced. The sequences obtained were aligned with co-specific sequences available in the GenBank database for primer design and phylogenetic studies. Three pairs of primers were chosen to compose the mPCR assay, which was used to verify the frequency of different mites in house dust samples (n = 20) from homes of Salvador, Brazil. Blomia tropicalis was the most frequent, found in 95% of the samples, followed by G. malaysiensis (70%) and D. pteronyssinus (60%). Besides reporting for the first time the occurrence of G. malaysiensis in Brazil, our results confirm the good resolution of the ITS2 region for mite identification. Furthermore, the mPCR assay proved to be a fast and reliable tool for identifying these mites in mixed cultures and could be applied in future epidemiological studies, and for quality control of mite extract production for general use.

The insect reservoir of biodiversity for viruses and for antiviral mechanisms.

Insects are the most diverse group of animals. They can be infected by an extraordinary diversity of viruses. Among them, arthropod-borne viruses (arboviruses) can be transmitted to humans. High-throughput sequencing of small RNAs from insects provides insight into their virome, which may help understand the dynamics of vector borne infectious diseases. Furthermore, investigating the mechanisms that restrict viral infections in insects points to genetic innovations that may inspire novel antiviral strategies.

Culex quinquefasciatus mosquitoes do not support replication of Zika virus.

The rapid spread of Zika virus (ZIKV) in the Americas raised many questions about the role of Culex quinquefasciatus mosquitoes in transmission, in addition to the key role played by the vector Aedes aegypti. Here we analysed the competence of Cx. quinquefasciatus (with or without Wolbachia endosymbionts) for a ZIKV isolate. We also examined the induction of RNA interference pathways after viral challenge and the production of small virus-derived RNAs. We did not observe any infection nor such small virus-derived RNAs, regardless of the presence or absence of Wolbachia. Thus, Cx. quinquefasciatus does not support ZIKV replication and Wolbachia is not involved in producing this phenotype. In short, these mosquitoes are very unlikely to play a role in transmission of ZIKV.

Virome analyses of Hevea brasiliensis using small RNA deep sequencing and PCR techniques reveal the presence of a potential new virus.

BACKGROUND: Hevea brasiliensis is an important commercial crop due to the high quality of the latex it produces; however, little is known about viral infections in this plant. The only virus described to infect H. brasiliensis until now is a Carlavirus, which was described more than 30 years ago. Virus-derived small interfering RNA (vsiRNAs) are the product of the plant's antiviral defense triggered by dsRNA viral intermediates generated, during the replication cycle. These vsiRNAs are complementar to viral genomes and have been widely used to identify and characterize viruses in plants. METHODS: In the present study, we investigated the virome of leaf and sapwood samples from native H. brasiliensis trees collected in two geographic areas in the Brazilian Amazon. Small RNA (sRNA) deep sequencing and bioinformatic tools were used to assembly, identify and characterize viral contigs. Subsequently, PCR amplification techniques were performed to experimentally verify the presence of the viral sequences. Finally, the phylogenetic relationship of the putative new virus with related viral genomes was analyzed. RESULTS: Our strategy allowed the identification of 32 contigs with high similarity to viral reference genomes, from which 23 exhibited homology to viruses of the Tymoviridae family. The reads showed a predominant size distribution at 21 nt derived from both strands, which was consistent with the vsiRNAs profile. The presence and genome position of the viral contigs were experimentally confirmed using droplet digital PCR amplifications. A 1913 aa long fragment was obtained and used to infer the phylogenetic relationship of the putative new virus, which indicated that it is taxonomically related to the Grapevine fleck virus, genus Maculavirus. The putative new virus was named Hevea brasiliensis virus (HBrV) in reference to its host. CONCLUSION: The methodological strategy applied here proved to be efficient in detecting and confirming the presence of new viral sequences on a 'very difficult to manage' sample. This is the second time that viral sequences, that could be ascribed as a putative novel virus, associated to the rubber tree has been identified.

Sequence-independent characterization of viruses based on the pattern of viral small RNAs produced by the host.

Virus surveillance in vector insects is potentially of great benefit to public health. Large-scale sequencing of small and long RNAs has previously been used to detect viruses, but without any formal comparison of different strategies. Furthermore, the identification of viral sequences largely depends on similarity searches against reference databases. Here, we developed a sequence-independent strategy based on virus-derived small RNAs produced by the host response, such as the RNA interference pathway. In insects, we compared sequences of small and long RNAs, demonstrating that viral sequences are enriched in the small RNA fraction. We also noted that the small RNA size profile is a unique signature for each virus and can be used to identify novel viral sequences without known relatives in reference databases. Using this strategy, we characterized six novel viruses in the viromes of laboratory fruit flies and wild populations of two insect vectors: mosquitoes and sandflies. We also show that the small RNA profile could be used to infer viral tropism for ovaries among other aspects of virus biology. Additionally, our results suggest that virus detection utilizing small RNAs can also be applied to vertebrates, although not as efficiently as to plants and insects.

Functional specialization of the small interfering RNA pathway in response to virus infection.

In Drosophila, post-transcriptional gene silencing occurs when exogenous or endogenous double stranded RNA (dsRNA) is processed into small interfering RNAs (siRNAs) by Dicer-2 (Dcr-2) in association with a dsRNA-binding protein (dsRBP) cofactor called Loquacious (Loqs-PD). siRNAs are then loaded onto Argonaute-2 (Ago2) by the action of Dcr-2 with another dsRBP cofactor called R2D2. Loaded Ago2 executes the destruction of target RNAs that have sequence complementarity to siRNAs. Although Dcr-2, R2D2, and Ago2 are essential for innate antiviral defense, the mechanism of virus-derived siRNA (vsiRNA) biogenesis and viral target inhibition remains unclear. Here, we characterize the response mechanism mediated by siRNAs against two different RNA viruses that infect Drosophila. In both cases, we show that vsiRNAs are generated by Dcr-2 processing of dsRNA formed during viral genome replication and, to a lesser extent, viral transcription. These vsiRNAs seem to preferentially target viral polyadenylated RNA to inhibit viral replication. Loqs-PD is completely dispensable for silencing of the viruses, in contrast to its role in silencing endogenous targets. Biogenesis of vsiRNAs is independent of both Loqs-PD and R2D2. R2D2, however, is required for sorting and loading of vsiRNAs onto Ago2 and inhibition of viral RNA expression. Direct injection of viral RNA into Drosophila results in replication that is also independent of Loqs-PD. This suggests that triggering of the antiviral pathway is not related to viral mode of entry but recognition of intrinsic features of virus RNA. Our results indicate the existence of a vsiRNA pathway that is separate from the endogenous siRNA pathway and is specifically triggered by virus RNA. We speculate that this unique framework might be necessary for a prompt and efficient antiviral response.

SchistoDB: an updated genome resource for the three key schistosomes of humans.

The new release of SchistoDB (http://SchistoDB.net) provides a rich resource of genomic data for key blood flukes (genus Schistosoma) which cause disease in hundreds of millions of people worldwide. SchistoDB integrates whole-genome sequence and annotation of three species of the genus and provides enhanced bioinformatics analyses and data-mining tools. A simple, yet comprehensive web interface provided through the Strategies Web Development Kit is available for the mining and visualization of the data. Genomic scale data can be queried based on BLAST searches, annotation keywords and gene ID searches, gene ontology terms, sequence motifs, protein characteristics and phylogenetic relationships. Search strategies can be saved within a user's profile for future retrieval and may also be shared with other researchers using a unique web address.

Samba virus: a novel mimivirus from a giant rain forest, the Brazilian Amazon.

BACKGROUND: The identification of novel giant viruses from the nucleocytoplasmic large DNA viruses group and their virophages has increased in the last decade and has helped to shed light on viral evolution. This study describe the discovery, isolation and characterization of Samba virus (SMBV), a novel giant virus belonging to the Mimivirus genus, which was isolated from the Negro River in the Brazilian Amazon. We also report the isolation of an SMBV-associated virophage named Rio Negro (RNV), which is the first Mimivirus virophage to be isolated in the Americas. METHODS/RESULTS: Based on a phylogenetic analysis, SMBV belongs to group A of the putative Megavirales order, possibly a new virus related to Acanthamoeba polyphaga mimivirus (APMV). SMBV is the largest virus isolated in Brazil, with an average particle diameter about 574 nm. The SMBV genome contains 938 ORFs, of which nine are ORFans. The 1,213.6 kb SMBV genome is one of the largest genome of any group A Mimivirus described to date. Electron microscopy showed RNV particle accumulation near SMBV and APMV factories resulting in the production of defective SMBV and APMV particles and decreasing the infectivity of these two viruses by several logs. CONCLUSION: This discovery expands our knowledge of Mimiviridae evolution and ecology.

Circulation of bee-infecting viruses in Brazil: a call for action.

Bees are fundamental for maintaining pollination-dependent plant populations, both economically and ecologically. In Brazil, they constitute 66.3% of pollinators, contributing to an annual market value estimated at R$ 43 billion for pollination services. Unfortunately, worldwide bee populations are declining due to parasites and pathogens, more specifically viruses, alongside climate change, habitat loss, and pesticides. In this scenario, extensive research concerning bee diversity, virus diversity and surveillance, is necessary to aid the conservation of native managed pollinators and potential wild alternatives besides mitigating the emergence and spread of viral pathogens. A decrease in pollination can be a point of economic vulnerability in a country like Brazil because of its main dependence on food exports. Here we conducted a study aiming to obtain an overview of circulating viruses in bees within Brazilian territory highlighting the need for further studies to have a more realistic view of bee-infecting viruses in Brazil.

Rhabdoviral Endogenous Sequences Identified in the Leishmaniasis Vector Lutzomyia longipalpis Are Widespread in Sandflies from South America.

Sandflies are known vectors of leishmaniasis. In the Old World, sandflies are also vectors of viruses while little is known about the capacity of New World insects to transmit viruses to humans. Here, we relate the identification of RNA sequences with homology to rhabdovirus nucleocapsids (NcPs) genes, initially in the Lutzomyia longipalpis LL5 cell lineage, named NcP1.1 and NcP2. The Rhabdoviridae family never retrotranscribes its RNA genome to DNA. The sequences here described were identified in cDNA and DNA from LL-5 cells and in adult insects indicating that they are transcribed endogenous viral elements (EVEs). The presence of NcP1.1 and NcP2 in the L. longipalpis genome was confirmed in silico. In addition to showing the genomic location of NcP1.1 and NcP2, we identified another rhabdoviral insertion named NcP1.2. Analysis of small RNA molecules derived from these sequences showed that NcP1.1 and NcP1.2 present a profile consistent with elements targeted by primary piRNAs, while NcP2 was restricted to the degradation profile. The presence of NcP1.1 and NcP2 was investigated in sandfly populations from South America and the Old World. These EVEs are shared by different sandfly populations in South America while none of the Old World species studied presented the insertions.