A chromosomal-scale reference genome of the New World Screwworm, Cochliomyia hominivorax.

The New World Screwworm, Cochliomyia hominivorax (Calliphoridae), is the most important myiasis-causing species in America. Screwworm myiasis is a zoonosis that can cause severe lesions in livestock, domesticated and wild animals, and occasionally in people. Beyond the sanitary problems associated with this species, these infestations negatively impact economic sectors, such as the cattle industry. Here, we present a chromosome-scale assembly of C. hominivorax's genome, organized in 6 chromosome-length and 515 unplaced scaffolds spanning 534 Mb. There was a clear correspondence between the D. melanogaster linkage groups A-E and the chromosomal-scale scaffolds. Chromosome quotient (CQ) analysis identified a single scaffold from the X chromosome that contains most of the orthologs of genes that are on the D. melanogaster fourth chromosome (linkage group F or dot chromosome). CQ analysis also identified potential X and Y unplaced scaffolds and genes. Y-linkage for selected regions was confirmed by PCR with male and female DNA. Some of the long chromosome-scale scaffolds include Y-linked sequences, suggesting misassembly of these regions. These resources will provide a basis for future studies aiming at understanding the biology and evolution of this devastating obligate parasite.

Molecular and Developmental Signatures of Genital Size Macro-Evolution in Bugs.

Our understanding of the genetic architecture of phenotypic traits has experienced drastic growth over the last years. Nevertheless, the majority of studies associating genotypes and phenotypes have been conducted at the ontogenetic level. Thus, we still have an elusive knowledge of how these genetic-developmental architectures evolve themselves and how their evolution is mirrored in the phenotypic change across evolutionary time. We tackle this gap by reconstructing the evolution of male genital size, one of the most complex traits in insects, together with its underlying genetic architecture. Using the order Hemiptera as a model, spanning over 350 million years of evolution, we estimate the correlation between genitalia and three features: development rate, body size, and rates of DNA substitution in 68 genes associated with genital development. We demonstrate that genital size macro-evolution has been largely dependent on body size and weakly influenced by development rate and phylogenetic history. We further revealed significant correlations between mutation rates and genital size for 19 genes. Interestingly, these genes have diverse functions and participate in distinct signaling pathways, suggesting that genital size is a complex trait whose fast evolution has been enabled by molecular changes associated with diverse morphogenetic processes. Our data further demonstrate that the majority of DNA evolution correlated with the genitalia has been shaped by negative selection or neutral evolution. Thus, in terms of sequence evolution, changes in genital size are predominantly facilitated by relaxation of constraints rather than positive selection, possibly due to the high pleiotropic nature of the morphogenetic genes.

The steroid-hormone ecdysone coordinates parallel pupariation neuromotor and morphogenetic subprograms via epidermis-to-neuron Dilp8-Lgr3 signal induction.

Innate behaviors consist of a succession of genetically-hardwired motor and physiological subprograms that can be coupled to drastic morphogenetic changes. How these integrative responses are orchestrated is not completely understood. Here, we provide insight into these mechanisms by studying pupariation, a multi-step innate behavior of Drosophila larvae that is critical for survival during metamorphosis. We find that the steroid-hormone ecdysone triggers parallel pupariation neuromotor and morphogenetic subprograms, which include the induction of the relaxin-peptide hormone, Dilp8, in the epidermis. Dilp8 acts on six Lgr3-positive thoracic interneurons to couple both subprograms in time and to instruct neuromotor subprogram switching during behavior. Our work reveals that interorgan feedback gates progression between subunits of an innate behavior and points to an ancestral neuromodulatory function of relaxin signaling.

Genomic analyses of a livestock pest, the New World screwworm, find potential targets for genetic control programs.

The New World Screwworm fly, Cochliomyia hominivorax, is a major pest of livestock in South America and Caribbean. However, few genomic resources have been available for this species. A genome of 534 Mb was assembled from long read PacBio DNA sequencing of DNA from a highly inbred strain. Analysis of molecular evolution identified 40 genes that are likely under positive selection. Developmental RNA-seq analysis identified specific genes associated with each stage. We identify and analyze the expression of genes that are likely important for host-seeking behavior (chemosensory), development of larvae in open wounds in warm-blooded animals (heat shock protein, immune response) and for building transgenic strains for genetic control programs including gene drive (sex determination, germline). This study will underpin future experiments aimed at understanding the parasitic lifestyle of the screwworm fly and greatly facilitate future development of strains for efficient systems for genetic control of screwworm.

De novo construction of a transcriptome for the stink bug crop pest Chinavia impicticornis during late development.

Chinavia impicticornis is a neotropical stink bug of economic importance for various crops. Little is known about the development of the species, or the genetic mechanisms that may favor the establishment of populations in cultivated plants. Here, we conduct the first large-scale molecular study of C. impicticornis. Using tissues derived from the genitalia and the rest of the body for two immature stages of both males and females, we generated RNA-seq data, then assembled and functionally annotated a transcriptome. The de novo-assembled transcriptome contained around 400,000 contigs, with an average length of 688 bp. After pruning duplicated sequences and conducting a functional annotation, the final annotated transcriptome comprised 39,478 transcripts, of which 12,665 were assigned to Gene Ontology (GO) terms. These novel datasets will be invaluable for the discovery of molecular processes related to morphogenesis and immature biology. We hope to contribute to the growing body of research on stink bug evolution and development, as well as to the development of biorational pest management solutions.

Worldwide Phylogenetic Group Patterns of Escherichia coli from Commensal Human and Wastewater Treatment Plant Isolates.

Escherichia coli is an important microorganism in the gastrointestinal tract of warm-blooded animals. Commensal populations of E. coli consist of stable genetic isolates, which means that each individual has only one phylogenetic group (phylogroup). We evaluated the frequency of human commensal E. coli phylogroups from 116 people and observed that the majority of isolates belonged to group A. We also evaluated the frequency of phylogroups in wastewater samples and found a strong positive correlation between the phylogroup distribution in wastewater and human hosts. In order to find out if some factors, such as geographical location, and climate could influence the worldwide phylogroup distribution, we performed a meta-analysis of 39 different studies and 24 countries, including different climates, living areas, and feeding habits. Unexpectedly, our results showed no substructuring patterns of phylogroups; indicating there was no correlation between phylogroup distribution and geographic location, climate, living area, feeding habits, or date of collection.

Novel variants in GNAI3 associated with auriculocondylar syndrome strengthen a common dominant negative effect.

Auriculocondylar syndrome is a rare craniofacial disorder comprising core features of micrognathia, condyle dysplasia and question mark ear. Causative variants have been identified in PLCB4, GNAI3 and EDN1, which are predicted to function within the EDN1-EDNRA pathway during early pharyngeal arch patterning. To date, two GNAI3 variants in three families have been reported. Here we report three novel GNAI3 variants, one segregating with affected members in a family previously linked to 1p21.1-q23.3 and two de novo variants in simplex cases. Two variants occur in known functional motifs, the G1 and G4 boxes, and the third variant is one amino acid outside of the G1 box. Structural modeling shows that all five altered GNAI3 residues identified to date cluster in a region involved in GDP/GTP binding. We hypothesize that all GNAI3 variants lead to dominant negative effects.

Bacterial diversity assessment in soil of an active Brazilian copper mine using high-throughput sequencing of 16S rDNA amplicons.

Mining activities pose severe environmental risks worldwide, generating extreme pH conditions and high concentrations of heavy metals, which can have major impacts on the survival of organisms. In this work, pyrosequencing of the V3 region of the 16S rDNA was used to analyze the bacterial communities in soil samples from a Brazilian copper mine. For the analysis, soil samples were collected from the slopes (geotechnical structures) and the surrounding drainage of the Sossego mine (comprising the Sossego and Sequeirinho deposits). The results revealed complex bacterial diversity, and there was no influence of deposit geographic location on the composition of the communities. However, the environment type played an important role in bacterial community divergence; the composition and frequency of OTUs in the slope samples were different from those of the surrounding drainage samples, and Acidobacteria, Chloroflexi, Firmicutes, and Gammaproteobacteria were responsible for the observed difference. Chemical analysis indicated that both types of sample presented a high metal content, while the amounts of organic matter and water were higher in the surrounding drainage samples. Non-metric multidimensional scaling (N-MDS) analysis identified organic matter and water as important distinguishing factors between the bacterial communities from the two types of mine environment. Although habitat-specific OTUs were found in both environments, they were more abundant in the surrounding drainage samples (around 50 %), and contributed to the higher bacterial diversity found in this habitat. The slope samples were dominated by a smaller number of phyla, especially Firmicutes. The bacterial communities from the slope and surrounding drainage samples were different in structure and composition, and the organic matter and water present in these environments contributed to the observed differences.

Bacterial communities and species-specific associations with the mucus of Brazilian coral species.

We investigated the existence of species-specific associations between Brazilian coral species and bacteria. Pyrosequencing of the V3 region of the 16S rDNA was used to analyze the taxonomic composition of bacterial communities associated with the mucus of four coral species (Madracis decactis, Mussismilia hispida, Palythoa caribaeorum, and Tubastraea coccinea) in two seasons (winter and summer), which were compared with the surrounding water and sediment. The microbial communities found in samples of mucus, water, and sediment differed according to the composition and relative frequency of OTUs. The coral mucus community seemed to be more stable and resistant to seasonal variations, compared to the water and sediment communities. There was no influence of geographic location on the composition of the communities. The sediment community was extremely diverse and might act as a "seed bank" for the entire environment. Species-specific OTUs were found in P. caribaeorum, T. coccinea, and M. hispida.

Deep sequencing of New World screw-worm transcripts to discover genes involved in insecticide resistance.

BACKGROUND: The New World screw-worm (NWS), Cochliomyia hominivorax, is one of the most important myiasis-causing flies, causing severe losses to the livestock industry. In its current geographical distribution, this species has been controlled by the application of insecticides, mainly organophosphate (OP) compounds, but a number of lineages have been identified that are resistant to such chemicals. Despite its economic importance, only limited genetic information is available for the NWS. Here, as a part of an effort to characterize the C. hominivorax genome and identify putative genes involved in insecticide resistance, we sampled its transcriptome by deep sequencing of polyadenylated transcripts using the 454 sequencing technology. RESULTS: Deep sequencing on the 454 platform of three normalized libraries (larval, adult male and adult female) generated a total of 548,940 reads. Eighteen candidate genes coding for three metabolic detoxification enzyme families, cytochrome P450 monooxygenases, glutathione S-transferases and carboxyl/cholinesterases were selected and gene expression levels were measured using quantitative real-time polymerase chain reaction (qRT-PCR). Of the investigated candidates, only one gene was expressed differently between control and resistant larvae with, at least, a 10-fold down-regulation in the resistant larvae. The presence of mutations in the acetylcholinesterase (target site) and carboxylesterase E3 genes was investigated and all of the resistant flies presented E3 mutations previously associated with insecticide resistance. CONCLUSIONS: Here, we provided the largest database of NWS expressed sequence tags that is an important resource, not only for further studies on the molecular basis of the OP resistance in NWS fly, but also for functional and comparative studies among Calliphoridae flies. Among our candidates, only one gene was found differentially expressed in resistant individuals, and its role on insecticide resistance should be further investigated. Furthermore, the absence of mutations in the OP target site and the high frequency of mutant carboxylesterase E3 indicate that metabolic resistance mechanisms have evolved predominantly in this species.

A transcriptomic analysis of gene expression in the venom gland of the snake Bothrops alternatus (urutu).

BACKGROUND: The genus Bothrops is widespread throughout Central and South America and is the principal cause of snakebite in these regions. Transcriptomic and proteomic studies have examined the venom composition of several species in this genus, but many others remain to be studied. In this work, we used a transcriptomic approach to examine the venom gland genes of Bothrops alternatus, a clinically important species found in southeastern and southern Brazil, Uruguay, northern Argentina and eastern Paraguay. RESULTS: A cDNA library of 5,350 expressed sequence tags (ESTs) was produced and assembled into 838 contigs and 4512 singletons. BLAST searches of relevant databases showed 30% hits and 70% no-hits, with toxin-related transcripts accounting for 23% and 78% of the total transcripts and hits, respectively. Gene ontology analysis identified non-toxin genes related to general metabolism, transcription and translation, processing and sorting, (polypeptide) degradation, structural functions and cell regulation. The major groups of toxin transcripts identified were metalloproteinases (81%), bradykinin-potentiating peptides/C-type natriuretic peptides (8.8%), phospholipases A2 (5.6%), serine proteinases (1.9%) and C-type lectins (1.5%). Metalloproteinases were almost exclusively type PIII proteins, with few type PII and no type PI proteins. Phospholipases A2 were essentially acidic; no basic PLA2 were detected. Minor toxin transcripts were related to L-amino acid oxidase, cysteine-rich secretory proteins, dipeptidylpeptidase IV, hyaluronidase, three-finger toxins and ohanin. Two non-toxic proteins, thioredoxin and double-specificity phosphatase Dusp6, showed high sequence identity to similar proteins from other snakes. In addition to the above features, single-nucleotide polymorphisms, microsatellites, transposable elements and inverted repeats that could contribute to toxin diversity were observed. CONCLUSIONS: Bothrops alternatus venom gland contains the major toxin classes described for other Bothrops venoms based on trancriptomic and proteomic studies. The predominance of type PIII metalloproteinases agrees with the well-known hemorrhagic activity of this venom, whereas the lower content of serine proteases and C-type lectins could contribute to less marked coagulopathy following envenoming by this species. The lack of basic PLA2 agrees with the lower myotoxicity of this venom compared to other Bothrops species with these toxins. Together, these results contribute to our understanding of the physiopathology of envenoming by this species.