Genomic targets for high-resolution inference of kinship, ancestry and disease susceptibility in orang-utans (genus: Pongo).

BACKGROUND: Orang-utans comprise three critically endangered species endemic to the islands of Borneo and Sumatra. Though whole-genome sequencing has recently accelerated our understanding of their evolutionary history, the costs of implementing routine genome screening and diagnostics remain prohibitive. Capitalizing on a tri-fold locus discovery approach, combining data from published whole-genome sequences, novel whole-exome sequencing, and microarray-derived genotype data, we aimed to develop a highly informative gene-focused panel of targets that can be used to address a broad range of research questions. RESULTS: We identified and present genomic co-ordinates for 175,186 SNPs and 2315 Y-chromosomal targets, plus 185 genes either known or presumed to be pathogenic in cardiovascular (N = 109) or respiratory (N = 43) diseases in humans - the primary and secondary causes of captive orang-utan mortality - or a majority of other human diseases (N = 33). As proof of concept, we designed and synthesized 'SeqCap' hybrid capture probes for these targets, demonstrating cost-effective target enrichment and reduced-representation sequencing. CONCLUSIONS: Our targets are of broad utility in studies of orang-utan ancestry, admixture and disease susceptibility and aetiology, and thus are of value in addressing questions key to the survival of these species. To facilitate comparative analyses, these targets could now be standardized for future orang-utan population genomic studies. The targets are broadly compatible with commercial target enrichment platforms and can be utilized as published here to synthesize applicable probes.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing.

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

Tracking Cancer Evolution Reveals Constrained Routes to Metastases: TRACERx Renal.

Clear-cell renal cell carcinoma (ccRCC) exhibits a broad range of metastatic phenotypes that have not been systematically studied to date. Here, we analyzed 575 primary and 335 metastatic biopsies across 100 patients with metastatic ccRCC, including two cases sampledat post-mortem. Metastatic competence was afforded by chromosome complexity, and we identify 9p loss as a highly selected event driving metastasis and ccRCC-related mortality (p = 0.0014). Distinct patterns of metastatic dissemination were observed, including rapid progression to multiple tissue sites seeded by primary tumors of monoclonal structure. By contrast, we observed attenuated progression in cases characterized by high primary tumor heterogeneity, with metastatic competence acquired gradually and initial progression to solitary metastasis. Finally, we observed early divergence of primitive ancestral clones and protracted latency of up to two decades as a feature of pancreatic metastases.

The Arabidopsis thaliana mobilome and its impact at the species level.

Transposable elements (TEs) are powerful motors of genome evolution yet a comprehensive assessment of recent transposition activity at the species level is lacking for most organisms. Here, using genome sequencing data for 211 Arabidopsis thaliana accessions taken from across the globe, we identify thousands of recent transposition events involving half of the 326 TE families annotated in this plant species. We further show that the composition and activity of the 'mobilome' vary extensively between accessions in relation to climate and genetic factors. Moreover, TEs insert equally throughout the genome and are rapidly purged by natural selection from gene-rich regions because they frequently affect genes, in multiple ways. Remarkably, loci controlling adaptive responses to the environment are the most frequent transposition targets observed. These findings demonstrate the pervasive, species-wide impact that a rich mobilome can have and the importance of transposition as a recurrent generator of large-effect alleles.

Dual RNA-seq of parasite and host reveals gene expression dynamics during filarial worm-mosquito interactions.

BACKGROUND: Parasite biology, by its very nature, cannot be understood without integrating it with that of the host, nor can the host response be adequately explained without considering the activity of the parasite. However, due to experimental limitations, molecular studies of parasite-host systems have been predominantly one-sided investigations focusing on either of the partners involved. Here, we conducted a dual RNA-seq time course analysis of filarial worm parasite and host mosquito to better understand the parasite processes underlying development in and interaction with the host tissue, from the establishment of infection to the development of infective-stage larva. METHODOLOGY/PRINCIPAL FINDINGS: Using the Brugia malayi-Aedes aegypti system, we report parasite gene transcription dynamics, which exhibited a highly ordered developmental program consisting of a series of cyclical and state-transitioning temporal patterns. In addition, we contextualized these parasite data in relation to the concurrent dynamics of the host transcriptome. Comparative analyses using uninfected tissues and different host strains revealed the influence of parasite development on host gene transcription as well as the influence of the host environment on parasite gene transcription. We also critically evaluated the life-cycle transcriptome of B. malayi by comparing developmental stages in the mosquito relative to those in the mammalian host, providing insight into gene expression changes underpinning the mosquito-borne parasitic lifestyle of this heteroxenous parasite. CONCLUSIONS/SIGNIFICANCE: The data presented herein provide the research community with information to design wet lab experiments and select candidates for future study to more fully dissect the whole set of molecular interactions of both organisms in this mosquito-filarial worm symbiotic relationship. Furthermore, characterization of the transcriptional program over the complete life cycle of the parasite, including stages within the mosquito, could help devise novel targets for control strategies.

Tissue-enriched expression profiles in Aedes aegypti identify hemocyte-specific transcriptome responses to infection.

Hemocytes are integral components of mosquito immune mechanisms such as phagocytosis, melanization, and production of antimicrobial peptides. However, our understanding of hemocyte-specific molecular processes and their contribution to shaping the host immune response remains limited. To better understand the immunophysiological features distinctive of hemocytes, we conducted genome-wide analysis of hemocyte-enriched transcripts, and examined how tissue-enriched expression patterns change with the immune status of the host. Our microarray data indicate that the hemocyte-enriched trascriptome is dynamic and context-dependent. Analysis of transcripts enriched after bacterial challenge in circulating hemocytes with respect to carcass added a dimension to evaluating infection-responsive genes and immune-related gene families. We resolved patterns of transcriptional change unique to hemocytes from those that are likely shared by other immune responsive tissues, and identified clusters of genes preferentially induced in hemocytes, likely reflecting their involvement in cell type specific functions. In addition, the study revealed conserved hemocyte-enriched molecular repertoires, which might be implicated in core hemocyte function by cross-species meta-analysis of microarray expression data from Anopheles gambiae and Drosophila melanogaster.

A deep sequencing approach to comparatively analyze the transcriptome of lifecycle stages of the filarial worm, Brugia malayi.

BACKGROUND: Developing intervention strategies for the control of parasitic nematodes continues to be a significant challenge. Genomic and post-genomic approaches play an increasingly important role for providing fundamental molecular information about these parasites, thus enhancing basic as well as translational research. Here we report a comprehensive genome-wide survey of the developmental transcriptome of the human filarial parasite Brugia malayi. METHODOLOGY/PRINCIPAL FINDINGS: Using deep sequencing, we profiled the transcriptome of eggs and embryos, immature (≤3 days of age) and mature microfilariae (MF), third- and fourth-stage larvae (L3 and L4), and adult male and female worms. Comparative analysis across these stages provided a detailed overview of the molecular repertoires that define and differentiate distinct lifecycle stages of the parasite. Genome-wide assessment of the overall transcriptional variability indicated that the cuticle collagen family and those implicated in molting exhibit noticeably dynamic stage-dependent patterns. Of particular interest was the identification of genes displaying sex-biased or germline-enriched profiles due to their potential involvement in reproductive processes. The study also revealed discrete transcriptional changes during larval development, namely those accompanying the maturation of MF and the L3 to L4 transition that are vital in establishing successful infection in mosquito vectors and vertebrate hosts, respectively. CONCLUSIONS/SIGNIFICANCE: Characterization of the transcriptional program of the parasite's lifecycle is an important step toward understanding the developmental processes required for the infectious cycle. We find that the transcriptional program has a number of stage-specific pathways activated during worm development. In addition to advancing our understanding of transcriptome dynamics, these data will aid in the study of genome structure and organization by facilitating the identification of novel transcribed elements and splice variants.

Genome sequence of the plant-pathogenic bacterium Dickeya dadantii 3937.

Dickeya dadantii is a plant-pathogenic enterobacterium responsible for the soft rot disease of many plants of economic importance. We present here the sequence of strain 3937, a strain widely used as a model system for research on the molecular biology and pathogenicity of this group of bacteria.

Pathogenomics of Culex quinquefasciatus and meta-analysis of infection responses to diverse pathogens.

The mosquito Culex quinquefasciatus poses a substantial threat to human and veterinary health as a primary vector of West Nile virus (WNV), the filarial worm Wuchereria bancrofti, and an avian malaria parasite. Comparative phylogenomics revealed an expanded canonical C. quinquefasciatus immune gene repertoire compared with those of Aedes aegypti and Anopheles gambiae. Transcriptomic analysis of C. quinquefasciatus genes responsive to WNV, W. bancrofti, and non-native bacteria facilitated an unprecedented meta-analysis of 25 vector-pathogen interactions involving arboviruses, filarial worms, bacteria, and malaria parasites, revealing common and distinct responses to these pathogen types in three mosquito genera. Our findings provide support for the hypothesis that mosquito-borne pathogens have evolved to evade innate immune responses in three vector mosquito species of major medical importance.

Transcriptome changes in Culex quinquefasciatus (Diptera: Culicidae) salivary glands during West Nile virus infection.

Persistent West Nile virus (WNV) infection in the mosquito Culex quinquefasciatus Say (Diptera: Culicidae) is associated with pathological changes in the salivary glands, including apoptotic cell death and a corresponding reduction in virus transmission over time. The vector host response to WNV infection and the molecular basis of WNV pathogenesis in Cx. quinquefasciatus was investigated using oligonucleotide microarrays designed to detect differences in the salivary gland transcriptome between WNV-infected mosquitoes and uninfected controls. Transcripts with increased abundance in infected salivary glands included those related to immunity, transcription, protein transport and degradation, amino acid and nucleotide metabolism, signal transduction, and cellular detoxification. Microarray-based analysis detected a decrease in transcript levels of a Culex inhibitor of apoptosis gene (IAP-1) and a decrease in abundance of 11 transcripts encoding salivary gland proteins. Transcript levels for an endonuclease, a proline-rich mucin, and several D7 protein family members also decreased. Transcripts with the greatest change in abundance during infection had either no similarity to sequences found in GenBank, VectorBase, and FlyBase, or were similar to sequences with uncharacterized protein products. These transcripts represent exciting targets for future analysis. Results from this study suggest that WNV infection influences transcriptional changes in an invertebrate host target tissue that may confer an advantage to the replicating virus, induce a host defense response, and alter the composition of vector saliva. The ramifications of these changes are discussed in terms of mosquito vector competence and WNV pathogenesis.

Use of microarray hybridization to identify Brugia genes involved in mosquito infectivity.

Brugia malayi and Brugia pahangi microfilariae (mf) require a maturation period of at least 5 days in the mammalian host to successfully infect laboratory mosquitoes. This maturation process coincides with changes in the surface composition of mf that likely are associated with changes in gene expression. To test this hypothesis, we verified the differential infectivity of immature (< or =3 day) and mature (>30 day) Brugia mf for black-eyed Liverpool strain of Aedes aegypti and then assessed transcriptome changes associated with microfilarial maturation by competitively hybridizing microfilarial cDNAs to the B. malayi oligonucleotide microarray. We identified transcripts differentially abundant in immature (94 in B. pahangi and 29 in B. malayi) and mature (64 in B. pahangi and 14 in B. malayi) mf. In each case, >40% of Brugia transcripts shared no similarity to known genes or were similar to genes with unknown function; the remaining transcripts were categorized by putative function based on sequence similarity to known genes/proteins. Microfilarial maturation was not associated with demonstrable changes in the abundance of transmembrane or secreted proteins; however, immature mf expressed more transcripts associated with immune modulation, neurotransmission, transcription, and cellular cytoskeleton elements, while mature mf displayed increased transcripts potentially encoding hypodermal/muscle and surface molecules, e.g., cuticular collagens and sheath components. The results of the homologous B. malayi microarray hybridization were validated by quantitative reverse transcriptase polymerase chain reaction. These findings preliminarily lend support to the underlying hypothesis that changes in microfilarial gene expression drive maturation-associated changes that influence the parasite to develop in compatible vectors.

Mosquito infection responses to developing filarial worms.

Human lymphatic filariasis is a mosquito-vectored disease caused by the nematode parasites Wuchereria bancrofti, Brugia malayi and Brugia timori. These are relatively large roundworms that can cause considerable damage in compatible mosquito vectors. In order to assess how mosquitoes respond to infection in compatible mosquito-filarial worm associations, microarray analysis was used to evaluate transcriptome changes in Aedes aegypti at various times during B. malayi development. Changes in transcript abundance in response to the different stages of B. malayi infection were diverse. At the early stages of midgut and thoracic muscle cell penetration, a greater number of genes were repressed compared to those that were induced (20 vs. 8). The non-feeding, intracellular first-stage larvae elicited few differences, with 4 transcripts showing an increased and 9 a decreased abundance relative to controls. Several cecropin transcripts increased in abundance after parasites molted to second-stage larvae. However, the greatest number of transcripts changed in abundance after larvae molted to third-stage larvae and migrated to the head and proboscis (120 induced, 38 repressed), including a large number of putative, immunity-related genes (approximately 13% of genes with predicted functions). To test whether the innate immune system of mosquitoes was capable of modulating permissiveness to the parasite, we activated the Toll and Imd pathway controlled rel family transcription factors Rel1 and Rel2 (by RNA interference knockdown of the pathway's negative regulators Cactus and Caspar) during the early stages of infection with B. malayi. The activation of either of these immune signaling pathways, or knockdown of the Toll pathway, did not affect B. malayi in Ae. aegypti. The possibility of LF parasites evading mosquito immune responses during successful development is discussed.

Construction and characterization of an expressed sequenced tag library for the mosquito vector Armigeres subalbatus.

BACKGROUND: The mosquito, Armigeres subalbatus, mounts a distinctively robust innate immune response when infected with the nematode Brugia malayi, a causative agent of lymphatic filariasis. In order to mine the transcriptome for new insight into the cascade of events that takes place in response to infection in this mosquito, 6 cDNA libraries were generated from tissues of adult female mosquitoes subjected to immune-response activation treatments that lead to well-characterized responses, and from aging, naïve mosquitoes. Expressed sequence tags (ESTs) from each library were produced, annotated, and subjected to comparative analyses. RESULTS: Six libraries were constructed and used to generate 44,940 expressed sequence tags, of which 38,079 passed quality filters to be included in the annotation project and subsequent analyses. All of these sequences were collapsed into clusters resulting in 8,020 unique sequence clusters or singletons. EST clusters were annotated and curated manually within ASAP (A Systematic Annotation Package for Community Analysis of Genomes) web portal according to BLAST results from comparisons to Genbank, and the Anopheles gambiae and Drosophila melanogaster genome projects. CONCLUSION: The resulting dataset is the first of its kind for this mosquito vector and provides a basis for future studies of mosquito vectors regarding the cascade of events that occurs in response to infection, and thereby providing insight into vector competence and innate immunity.

Mosquito transcriptome changes and filarial worm resistance in Armigeres subalbatus.

BACKGROUND: Armigeres subalbatus is a natural vector of the filarial worm Brugia pahangi, but it rapidly and proficiently kills Brugia malayi microfilariae by melanotic encapsulation. Because B. malayi and B. pahangi are morphologically and biologically similar, the Armigeres-Brugia system serves as a valuable model for studying the resistance mechanisms in mosquito vectors. We have initiated transcriptome profiling studies in Ar. subalbatus to identify molecular components involved in B. malayi refractoriness. RESULTS: These initial studies assessed the transcriptional response of Ar. subalbatus to B. malayi at 1, 3, 6, 12, 24, 48, and 72 hrs after an infective blood feed. In this investigation, we initiated the first holistic study conducted on the anti-filarial worm immune response in order to effectively explore the functional roles of immune-response genes following a natural exposure to the parasite. Studies assessing the transcriptional response revealed the involvement of unknown and conserved unknowns, cytoskeletal and structural components, and stress and immune responsive factors. The data show that the anti-filarial worm immune response by Ar. subalbatus to be a highly complex, tissue-specific process involving varied effector responses working in concert with blood cell-mediated melanization. CONCLUSION: This initial study provides a foundation and direction for future studies, which will more fully dissect the nature of the anti-filarial worm immune response in this mosquito-parasite system. The study also argues for continued studies with RNA generated from both hemocytes and whole bodies to fully expound the nature of the anti-filarial worm immune response.

Evolutionary dynamics of immune-related genes and pathways in disease-vector mosquitoes.

Mosquitoes are vectors of parasitic and viral diseases of immense importance for public health. The acquisition of the genome sequence of the yellow fever and Dengue vector, Aedes aegypti (Aa), has enabled a comparative phylogenomic analysis of the insect immune repertoire: in Aa, the malaria vector Anopheles gambiae (Ag), and the fruit fly Drosophila melanogaster (Dm). Analysis of immune signaling pathways and response modules reveals both conservative and rapidly evolving features associated with different functional gene categories and particular aspects of immune reactions. These dynamics reflect in part continuous readjustment between accommodation and rejection of pathogens and suggest how innate immunity may have evolved.

Comparative genomics of Salmonella enterica serovar Typhi strains Ty2 and CT18.

We present the 4.8-Mb complete genome sequence of Salmonella enterica serovar Typhi strain Ty2, a human-specific pathogen causing typhoid fever. A comparison with the genome sequence of recently isolated S. enterica serovar Typhi strain CT18 showed that 29 of the 4,646 predicted genes in Ty2 are unique to this strain, while 84 genes are unique to CT18. Both genomes contain more than 200 pseudogenes; 9 of these genes in CT18 are intact in Ty2, while 11 intact CT18 genes are pseudogenes in Ty2. A half-genome interreplichore inversion in Ty2 relative to CT18 was confirmed. The two strains exhibit differences in prophages, insertion sequences, and island structures. While CT18 carries two plasmids, one conferring multiple drug resistance, Ty2 has no plasmids and is sensitive to antibiotics.

Genome sequence of Yersinia pestis KIM.

We present the complete genome sequence of Yersinia pestis KIM, the etiologic agent of bubonic and pneumonic plague. The strain KIM, biovar Mediaevalis, is associated with the second pandemic, including the Black Death. The 4.6-Mb genome encodes 4,198 open reading frames (ORFs). The origin, terminus, and most genes encoding DNA replication proteins are similar to those of Escherichia coli K-12. The KIM genome sequence was compared with that of Y. pestis CO92, biovar Orientalis, revealing homologous sequences but a remarkable amount of genome rearrangement for strains so closely related. The differences appear to result from multiple inversions of genome segments at insertion sequences, in a manner consistent with present knowledge of replication and recombination. There are few differences attributable to horizontal transfer. The KIM and E. coli K-12 genome proteins were also compared, exposing surprising amounts of locally colinear "backbone," or synteny, that is not discernible at the nucleotide level. Nearly 54% of KIM ORFs are significantly similar to K-12 proteins, with conserved housekeeping functions. However, a number of E. coli pathways and transport systems and at least one global regulator were not found, reflecting differences in lifestyle between them. In KIM-specific islands, new genes encode candidate pathogenicity proteins, including iron transport systems, putative adhesins, toxins, and fimbriae.