Donor plasma VEGF-A as a biomarker for myocardial injury and primary graft dysfunction after heart transplantation.

BACKGROUND: Vascular endothelial growth factor (VEGF)-A is an angiogenic and proinflammatory cytokine with profound effects on microvascular permeability and vasodilation. Several processes may induce VEGF-A expression in brain-dead organ donors. However, it remains unclear whether donor VEGF-A is linked to adverse outcomes after heart transplantation. METHODS: We examined plasma VEGF-A levels from 83 heart transplant donors as well as the clinical data of these donors and their respective recipients operated between 2010 and 2016. The donor plasma was analyzed using Luminex-based Multiplex and confirmed with a single-target ELISA. Based on donor VEGF-A plasma levels, the recipients were divided into 3 equal-sized groups (low VEGF <500 ng/liter, n = 28; moderate VEGF 500-3000 ng/liter, n = 28; and high VEGF >3000 ng/liter, n = 27). Biochemical and clinical parameters of myocardial injury as well as heart transplant and kidney function were followed-up for one year, while rejection episodes, development of cardiac allograft vasculopathy, and mortality were monitored for 5 years. RESULTS: Baseline parameters were comparable between the donor groups, except for age, where median ages of 40, 45, and 50 were observed for low, moderate, and high donor plasma VEGF levels groups, respectively, and therefore donor age was included as a confounding factor. High donor plasma VEGF-A levels were associated with pronounced myocardial injury (TnT and TnI), a higher inotrope score, and a higher incidence of primary graft dysfunction in the recipient after heart transplantation. Furthermore, recipients with allografts from donors with high plasma VEGF-A levels had a longer length of stay in the intensive care unit and the hospital, and an increased likelihood for prolonged renal replacement therapy. CONCLUSIONS: Our findings suggest that elevated donor plasma VEGF-A levels were associated with adverse outcomes in heart transplant recipients, particularly in terms of myocardial injury, primary graft dysfunction, and long-term renal complications. Donor VEGF-A may serve as a potential biomarker for predicting these adverse outcomes and identifying extended donor criteria.

Transcriptomic Landscape of Circulating Extracellular Vesicles in Heart Transplant Ischemia-Reperfusion.

Ischemia-reperfusion injury (IRI) is an inevitable event during heart transplantation, which is known to exacerbate damage to the allograft. However, the precise mechanisms underlying IRI remain incompletely understood. Here, we profiled the whole transcriptome of plasma extracellular vesicles (EVs) by RNA sequencing from 41 heart transplant recipients immediately before and at 12 h after transplant reperfusion. We found that the expression of 1317 protein-coding genes in plasma EVs was changed at 12 h after reperfusion. Upregulated genes of plasma EVs were related to metabolism and immune activation, while downregulated genes were related to cell survival and extracellular matrix organization. In addition, we performed correlation analyses between EV transcriptome and intensity of graft IRI (i.e., cardiomyocyte injury), as well as EV transcriptome and primary graft dysfunction, as well as any biopsy-proven acute rejection after heart transplantation. We ultimately revealed that at 12 h after reperfusion, 4 plasma EV genes (ITPKA, DDIT4L, CD19, and CYP4A11) correlated with both cardiomyocyte injury and primary graft dysfunction, suggesting that EVs are sensitive indicators of reperfusion injury reflecting lipid metabolism-induced stress and imbalance in calcium homeostasis. In conclusion, we show that profiling plasma EV gene expression may enlighten the mechanisms of heart transplant IRI.

Severe primary graft dysfunction of the heart transplant is associated with increased plasma and intragraft proinflammatory cytokine expression.

INTRODUCTION: Heart transplant results have constantly improved but primary left ventricle graft dysfunction (LV-PGD) remains a devastating complication early after transplantation. Donor and recipient systemic inflammatory response may be involved in immune activation of the transplant, and LV-PGD development. Here, we investigated donor and recipient plasma and intragraft cytokine profiles preoperatively and during LV-PGD and searched for predictive markers for LV-PGD. METHODS: Donor and recipient plasma samples (n = 74) and myocardial biopsies of heart transplants (n = 64) were analyzed. Plasma and intragraft cytokine levels were determined by multiplexed and next-generation sequencing platforms, respectively. The development of LV-PGD during the first 24 hours, and graft function and mortality up to 1 year after transplantation, were examined. RESULTS: Severe LV-PGD, but not mild or moderate LV-PGD, was significantly associated with early mortality, plasma high-sensitivity troponin elevation, and an increase in intragraft and plasma proinflammatory cytokines during reperfusion. Preoperative donor and recipient plasma cytokine levels failed to predict LV-PGD. Cytokine network analysis identified interleukins -6, -8, -10, and -18 as key players during reperfusion. Prolonged cold and total ischemia time, and increased need for red blood cell transfusions during operation were identified as clinical risk factors for severe LV-PGD. CONCLUSIONS: Severe LV-PGD was associated with a poor clinical outcome. Donor and recipient plasma cytokine profile failed to predict LV-PGD, but severe LV-PGD was associated with an increase in post-reperfusion intragraft and recipient plasma proinflammatory cytokines. Identified key cytokines may be potential therapeutic targets to improve early and long-term outcomes after heart transplantation.

Inhibition of Vascular Endothelial Growth Factor Receptors 1 and 2 Attenuates Natural Killer Cell and Innate Immune Responses in an Experimental Model for Obliterative Bronchiolitis.

Obliterative bronchiolitis (OB) after lung transplantation is a nonreversible, life-threatening complication. Herein, the role of vascular endothelial growth factor receptor (Vegfr)-1 and -2 was investigated in the development of obliterative airway disease (OAD), an experimental model for OB. The nonimmunosuppressed recipients underwent transplantation with fully major histocompatibility complex mismatched heterotopic tracheal allografts and received Vegfr1 and -2-specific monoclonal antibodies either alone or in combination, or rat IgG as a control. The treatment with Vegfr1- or -2-blocking antibody significantly decreased intragraft mRNA expression of natural killer cell activation markers early after transplantation. This was followed by reduced infiltration of Cd11b+ cells and Cd4+ T cells as well as down-regulated mRNA expression of proinflammatory chemokines and profibrotic growth factors. However, blocking of both Vegfr1 and -2 was necessary to reduce luminal occlusion. Furthermore, concomitant inhibition of the calcineurin activation pathway almost totally abolished the development of OAD. This study proposes that blocking of Vegf receptors blunted natural killer cell and innate immune responses early after transplantation and attenuated the development of OAD. The results of this study suggest that further studies on the role of Vegfr1 and -2 blocking in development of obliterative airway lesions might be rewarding.

Plasma proteome of brain-dead organ donors predicts heart transplant outcome.

BACKGROUND: The pathophysiological changes related to brain death may affect the quality of the transplanted organs and expose the recipients to risks. We probed systemic changes reflected in donor plasma proteome and investigated their relationship to heart transplant outcomes. METHODS: Plasma samples from brain-dead multi-organ donors were analyzed by label-free protein quantification using high-definition mass spectrometry. Unsupervised and supervised statistical models were used to determine proteome differences between brain-dead donors and healthy controls. Proteome variation and the corresponding biological pathways were analyzed and correlated with transplant outcomes. RESULTS: Statistical models revealed that donors had a unique but heterogeneous plasma proteome with 237 of 463 proteins being changed compared to controls. Pathway analysis showed that coagulation, gluconeogenesis, and glycolysis pathways were upregulated in donors, while complement, LXR/RXR activation, and production of nitric oxide and reactive oxygen species in macrophages pathways were downregulated. In point-biserial correlation analysis, lysine-specific demethylase 3A was moderately correlated with any grade and severe PGD. In univariate and multivariate Cox regression analyses myosin Va and proteasome activator complex subunit 2 were significantly associated with the development of acute rejections with hemodynamic compromise within 30 days. Finally, we found that elevated levels of lysine-specific demethylase 3A and moesin were identified as predictors for graft-related 1-year mortality in univariate analysis. CONCLUSIONS: We show that brain death significantly changed plasma proteome signature Donor plasma protein changes related to endothelial cell and cardiomyocyte function, inflammation, and vascular growth and arteriogenesis could predict transplant outcome suggesting a role in donor evaluation.

Genomics of asthma, allergy and chronic rhinosinusitis: novel concepts and relevance in airway mucosa.

Genome wide association studies (GWASs) have revealed several airway disease-associated risk loci. Their role in the onset of asthma, allergic rhinitis (AR) or chronic rhinosinusitis (CRS), however, is not yet fully understood. The aim of this review is to evaluate the airway relevance of loci and genes identified in GWAS studies. GWASs were searched from databases, and a list of loci associating significantly (p < 10-8) with asthma, AR and CRS was created. This yielded a total of 267 significantly asthma/AR-associated loci from 31 GWASs. No significant CRS -associated loci were found in this search. A total of 170 protein coding genes were connected to these loci. Of these, 76/170 (44%) showed bronchial epithelial protein expression in stained microscopic figures of Human Protein Atlas (HPA), and 61/170 (36%) had a literature report of having airway epithelial function. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analyses were performed, and 19 functional protein categories were found as significantly (p < 0.05) enriched among these genes. These were related to cytokine production, cell activation and adaptive immune response, and all were strongly connected in network analysis. We also identified 15 protein pathways that were significantly (p < 0.05) enriched in these genes, related to T-helper cell differentiation, virus infection, JAK-STAT signaling pathway, and asthma. A third of GWAS-level risk loci genes of asthma or AR seemed to have airway epithelial functions according to our database and literature searches. In addition, many of the risk loci genes were immunity related. Some risk loci genes also related to metabolism, neuro-musculoskeletal or other functions. Functions overlapped and formed a strong network in our pathway analyses and are worth future studies of biomarker and therapeutics.

Instability of natural selection at candidate barrier loci underlying speciation in wood ants.

Speciation underlies the generation of novel biodiversity. Yet, there is much to learn about how natural selection shapes genomes during speciation. Selection is assumed to act against gene flow at barrier loci, promoting reproductive isolation. However, evidence for gene flow and selection is often indirect and we know very little about the temporal stability of barrier loci. Here we utilize haplodiploidy to identify candidate male barrier loci in hybrids between two wood ant species. As ant males are haploid, they are expected to reveal recessive barrier loci, which can be masked in diploid females if heterozygous. We then test for barrier stability in a sample collected 10 years later and use survival analysis to provide a direct measure of natural selection acting on candidate male barrier loci. We find multiple candidate male barrier loci scattered throughout the genome. Surprisingly, a proportion of them are not stable after 10 years, natural selection apparently switching from acting against to favouring introgression in the later sample. Instability of the barrier effect and natural selection for introgressed alleles could be due to environment-dependent selection, emphasizing the need to consider temporal variation in the strength of natural selection and the stability of the barrier effect at putative barrier loci in future speciation work.

Maternal smoking during pregnancy affects adult onset of asthma in offspring: a follow up from birth to age 46†years.

RATIONALE: Environmental tobacco smoke (ETS) exposure increases asthma risk in children. There is limited knowledge of prenatal ETS for adult-onset asthma. OBJECTIVES: To determine the association between prenatal ETS and adult onset asthma. MEASUREMENTS AND MAIN RESULTS: The questionnaire and clinical data of 5200 people, free of physician-diagnosed asthma by 31 years of age, who were included in the Northern Finland Birth Cohort 1966 Study was used. The association of maternal smoking during the last 3 months of pregnancy with onset of physician-diagnosed asthma and with lung function in adult offspring was studied using adjusted multivariate regression analyses. The cumulative incidence of physician-diagnosed asthma between the ages of 31 and 46 years was 5.1% among men and 8.8% among women. Gestational smoke exposure was associated with adult-onset asthma among offspring (adjusted OR 1.54, 95% CI 1.04-2.29), namely among offspring who reported either past non-diagnosed asthma (OR 9.63, 95% CI 2.28-40.67) or past cough with wheeze (3.21, 95% CI 1.71-6.05). A significant association was detected between gestational smoke exposure and the offspring's forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) ratio at 31 years of age. In offspring with the haplotype rs11702779-AA of RUNX1, gestational smoke exposure was associated with adult-onset asthma (5.53, 95% CI 2.11-14.52, adjusted p-value for interaction 0.10). CONCLUSION: Maternal smoking during pregnancy is associated with the cumulative incidence of asthma in offspring between the ages of 31 and 46 years. The association was accentuated in offspring who at age 31, reported having past respiratory problems and/or who had haplotype rs11702779-AA. A reduction in FEV1/FVC ratio was also observed at age 31 years in offspring with gestational smoke exposure. These results could reflect the early vulnerability of offspring's airways to ETS and its putative long-term effects.

Evaluating responses to temperature during pre-metamorphosis and carry-over effects at post-metamorphosis in the wood tiger moth (Arctia plantaginis).

Insect metamorphosis is one of the most recognized processes delimiting transitions between phenotypes. It has been traditionally postulated as an adaptive process decoupling traits between life stages, allowing evolutionary independence of pre- and post-metamorphic phenotypes. However, the degree of autonomy between these life stages varies depending on the species and has not been studied in detail over multiple traits simultaneously. Here, we reared full-sib larvae of the warningly coloured wood tiger moth (Arctia plantaginis) in different temperatures and examined their responses for phenotypic (melanization change, number of moults), gene expression (RNA-seq and qPCR of candidate genes for melanization and flight performance) and life-histories traits (pupal weight, and larval and pupal ages). In the emerging adults, we examined their phenotypes (melanization and size) and compared them at three condition proxies: heat absorption (ability to engage flight), flight metabolism (ability to sustain flight) and overall flight performance. We found that some larval responses, as evidenced by gene expression and change in melanization, did not have an effect on the adult (i.e. size and wing melanization), whereas other adult traits such as heat absorption, body melanization and flight performance were found to be impacted by rearing temperature. Adults reared at high temperature showed higher resting metabolic rate, lower body melanization, faster heating rate, lower body temperature at take-off and inferior flight performance than cold-reared adults. Thus, our results did not unambiguously support the environment-matching hypothesis. Our results illustrate the importance of assessing multiple traits across life stages as these may only be partly decoupled by metamorphosis. This article is part of the theme issue 'The evolution of complete metamorphosis'.

Donor Simvastatin Treatment in Heart Transplantation.

BACKGROUND: Ischemia-reperfusion injury may compromise the short-term and long-term prognosis after heart transplantation. Experimental studies show that simvastatin administered to the organ donor is vasculoprotective and inhibits cardiac allograft ischemia-reperfusion injury. METHODS: Eighty-four multiorgan donors were randomly assigned to receive 80 mg of simvastatin (42 donors) via nasogastric tube after declaration of brain death and upon acceptance as a cardiac donor, or to receive no simvastatin (42 donors). The primary efficacy end point was postoperative plasma troponin T and I levels during the first 24 hours after heart transplantation. Secondary end points included postoperative hemodynamics, inflammation, allograft function, rejections and rejection treatments, and mortality. Results: Organ donor simvastatin treatment significantly reduced the heart recipient plasma levels of troponin T by 34% (14 900 ± 12 100 ng/L to 9800 ± 7900 ng/L, P=0.047), and troponin I by 40% (171 000 ± 151 000 ng/L to 103 000 ± 109 000 ng/L, P=0.023) at 6 hours after reperfusion, the levels of NT-proBNP (N-terminal pro-B-type natriuretic peptide) by 36% (32 800 ± 24 300 ng/L to 20 900 ± 15 900 ng/L; P=0.011) at 1 week, and the number of rejection treatments with hemodynamic compromise by 53% within the first 30 days (P=0.046). Donor simvastatin treatment did not affect donor lipid levels but was associated with a specific transplant myocardial biopsy gene expression profile, and a decrease in recipient postoperative plasma levels of CXCL10 (C-X-C motif chemokine 10), interleukin-1α, placental growth factor, and platelet-derived growth factor-BB. Postoperative hemodynamics, biopsy-proven acute rejections, and mortality were similar. No adverse effects were seen in recipients receiving noncardiac solid organ transplants from simvastatin-treated donors. CONCLUSIONS: Donor simvastatin treatment reduces biomarkers of myocardial injury after heart transplantation, and-also considering its documented general safety profile-may be used as a novel, safe, and inexpensive adjunct therapy in multiorgan donation. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov. Unique identifier: NCT01160978.

The first draft genomes of the ant Formica exsecta, and its Wolbachia endosymbiont reveal extensive gene transfer from endosymbiont to host.

BACKGROUND: Adapting to changes in the environment is the foundation of species survival, and is usually thought to be a gradual process. However, transposable elements (TEs), epigenetic modifications, and/or genetic material acquired from other organisms by means of horizontal gene transfer (HGTs), can also lead to novel adaptive traits. Social insects form dense societies, which attract and maintain extra- and intracellular accessory inhabitants, which may facilitate gene transfer between species. The wood ant Formica exsecta (Formicidae; Hymenoptera), is a common ant species throughout the Palearctic region. The species is a well-established model for studies of ecological characteristics and evolutionary conflict. RESULTS: In this study, we sequenced and assembled draft genomes for F. exsecta and its endosymbiont Wolbachia. The F. exsecta draft genome is 277.7 Mb long; we identify 13,767 protein coding genes, for which we provide gene ontology and protein domain annotations. This is also the first report of a Wolbachia genome from ants, and provides insights into the phylogenetic position of this endosymbiont. We also identified multiple horizontal gene transfer events (HGTs) from Wolbachia to F. exsecta. Some of these HGTs have also occurred in parallel in multiple other insect genomes, highlighting the extent of HGTs in eukaryotes. CONCLUSION: We present the first draft genome of ant F. exsecta, and its endosymbiont Wolbachia (wFex), and show considerable rates of gene transfer from the symbiont to the host. We expect that especially the F. exsecta genome will be valuable resource in further exploration of the molecular basis of the evolution of social organization.

Genome organization and molecular characterization of the three Formica exsecta viruses-FeV1, FeV2 and FeV4.

We present the genome organization and molecular characterization of the three Formica exsecta viruses, along with ORF predictions, and functional annotation of genes. The Formica exsecta virus-4 (FeV4; GenBank ID: MF287670) is a newly discovered negative-sense single-stranded RNA virus representing the first identified member of order Mononegavirales in ants, whereas the Formica exsecta virus-1 (FeV1; GenBank ID: KF500001), and the Formica exsecta virus-2 (FeV2; GenBank ID: KF500002) are positive single-stranded RNA viruses initially identified (but not characterized) in our earlier study. The new virus FeV4 was found by re-analyzing data from a study published earlier. The Formica exsecta virus-4 genome is 9,866 bp in size, with an overall G + C content of 44.92%, and containing five predicted open reading frames (ORFs). Our bioinformatics analysis indicates that gaps are absent and the ORFs are complete, which based on our comparative genomics analysis suggests that the genomes are complete. Following the characterization, we validate virus infection for FeV1, FeV2 and FeV4 for the first time in field-collected worker ants. Some colonies were infected by multiple viruses, and the viruses were observed to infect all castes, and multiple life stages of workers and queens. Finally, highly similar viruses were expressed in adult workers and queens of six other Formica species: F. fusca, F. pressilabris, F. pratensis, F. aquilonia, F. truncorum and F. cinerea. This research indicates that viruses can be shared between ant species, but further studies on viral transmission are needed to understand viral infection pathways.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta.

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

De novo transcriptome assembly and its annotation for the aposematic wood tiger moth (Parasemia plantaginis).

In this paper we report the public availability of transcriptome resources for the aposematic wood tiger moth (Parasemia plantaginis). A comprehensive assembly methods, quality statistics, and annotation are provided. This reference transcriptome may serve as a useful resource for investigating functional gene activity in aposematic Lepidopteran species. All data is freely available at the European Nucleotide Archive (http://www.ebi.ac.uk/ena) under study accession number: PRJEB14172.

Ancient Duplications Have Led to Functional Divergence of Vitellogenin-Like Genes Potentially Involved in Inflammation and Oxidative Stress in Honey Bees.

Protection against inflammation and oxidative stress is key in slowing down aging processes. The honey bee (Apis mellifera) shows flexible aging patterns linked to the social role of individual bees. One molecular factor associated with honey bee aging regulation is vitellogenin, a lipoglycophosphoprotein with anti-inflammatory and antioxidant properties. Recently, we identified three genes in Hymenopteran genomes arisen from ancient insect vitellogenin duplications, named vg-like-A, -B, and -C. The function of these vitellogenin homologs is unclear. We hypothesize that some of them might share gene- and protein-level similarities and a longevity-supporting role with vitellogenin. Here, we show how the structure and modifications of the vg-like genes and proteins have diverged from vitellogenin. Furthermore, all three vg-like genes show signs of positive selection, but the spatial location of the selected protein sites differ from those found in vitellogenin. We show that all these genes are expressed in both long-lived winter worker bees and in summer nurse bees with intermediate life expectancy, yet only vg-like-A shows elevated expression in winter bees as found in vitellogenin. Finally, we show that vg-like-A responds more strongly than vitellogenin to inflammatory and oxidative conditions in summer nurse bees, and that also vg-like-B responds to oxidative stress. We associate vg-like-A and, to lesser extent, vg-like-B to the antiaging roles of vitellogenin, but that vg-like-C probably is involved in some other function. Our analysis indicates that an ancient duplication event facilitated the adaptive and functional divergence of vitellogenin and its paralogs in the honey bee.

Not only for egg yolk--functional and evolutionary insights from expression, selection, and structural analyses of Formica ant vitellogenins.

Vitellogenin (Vg), a storage protein, has been extensively studied for its egg-yolk precursor role, and it has been suggested to be fundamentally involved in caste differences in social insects. More than one Vg copy has been reported in several oviparous species, including ants. However, the number and function of different Vgs, their phylogenetic relatedness, and their role in reproductive queens and nonreproductive workers have been studied in few species only. We studied caste-biased expression of Vgs in seven Formica ant species. Only one copy of conventional Vg was identified in Formica species, and three Vg homologs, derived from ancient duplications, which represent yet undiscovered Vg-like genes. We show that each of these Vg-like genes is present in all studied Hymenoptera and some of them in other insects as well. We show that after each major duplication event, at least one of the Vg-like genes has experienced a period of positive selection. This, combined with the observation that the Vg-like genes have acquired or lost specific protein domains suggests sub- or neofunctionalization between Vg and the duplicated genes. In contrast to earlier studies, Vg was not consistently queen biased in its expression, and the caste bias of the three Vg-like genes was highly variable among species. Furthermore, a truncated and Hymenoptera-specific Vg-like gene, Vg-like-C, was consistently worker biased. Multispecies comparisons are essential for Vg expression studies, and for gene expression studies in general, as we show that expression and also, putative functions cannot be generalized even among closely related species.

A metatranscriptomic approach to the identification of microbiota associated with the ant Formica exsecta.

Social insects live in cooperative colonies, often in high densities and with closely related individuals, and interact using social contact behaviours. Compared to solitary insects, social insects have evolved multi-level immunity that includes immune responses common to holometabolous insects, and social immunity, which is exclusive to social taxa. This suggests that social insects may be subject to high pathogen pressure, yet relatively little is known about the range of symbiotic and pathogenic microbial communities that associate with social insects. In this study we examined transcriptome data generated from the ant Formica exsecta for sequences identifying as microbes (or other organisms potentially of non-ant origin). Sequences showing homology to two viruses and several other potentially or obligate intracellular organisms, such as Wolbachia, Arsenophonus, Entomoplasmatales and Microsporidia, were present in the transcriptome data. These homologous sequence matches correspond to genera/species that have previously been associated with a variety of insects, including social insects. There were also sequences with identity to several other microbes such as common moulds and soil bacteria. We conclude that this sequence data provides a starting point for a deeper understanding of the biological interactions between a species of ant and the micro- and macrobiotic communities that it potentially encounters.

Identification of optimum sequencing depth especially for de novo genome assembly of small genomes using next generation sequencing data.

Next Generation Sequencing (NGS) is a disruptive technology that has found widespread acceptance in the life sciences research community. The high throughput and low cost of sequencing has encouraged researchers to undertake ambitious genomic projects, especially in de novo genome sequencing. Currently, NGS systems generate sequence data as short reads and de novo genome assembly using these short reads is computationally very intensive. Due to lower cost of sequencing and higher throughput, NGS systems now provide the ability to sequence genomes at high depth. However, currently no report is available highlighting the impact of high sequence depth on genome assembly using real data sets and multiple assembly algorithms. Recently, some studies have evaluated the impact of sequence coverage, error rate and average read length on genome assembly using multiple assembly algorithms, however, these evaluations were performed using simulated datasets. One limitation of using simulated datasets is that variables such as error rates, read length and coverage which are known to impact genome assembly are carefully controlled. Hence, this study was undertaken to identify the minimum depth of sequencing required for de novo assembly for different sized genomes using graph based assembly algorithms and real datasets. Illumina reads for E.coli (4.6 MB) S.kudriavzevii (11.18 MB) and C.elegans (100 MB) were assembled using SOAPdenovo, Velvet, ABySS, Meraculous and IDBA-UD. Our analysis shows that 50X is the optimum read depth for assembling these genomes using all assemblers except Meraculous which requires 100X read depth. Moreover, our analysis shows that de novo assembly from 50X read data requires only 6-40 GB RAM depending on the genome size and assembly algorithm used. We believe that this information can be extremely valuable for researchers in designing experiments and multiplexing which will enable optimum utilization of sequencing as well as analysis resources.