Immune response stability to the SARS-CoV-2 mRNA vaccine booster is influenced by differential splicing of HLA genes

Many molecular mechanisms that lead to the host antibody response to COVID-19 vaccines remain largely unknown. In this study, we used serum antibody detection combined with whole blood RNA-based transcriptome analysis to investigate variability in vaccine response in healthy recipients of a booster (third) dose schedule of the mRNA BNT162b2 vaccine against COVID-19. The cohort was divided into two groups: (1) low-stable individuals, with antibody concentration anti-SARS-CoV IgG S1 below 0.4 percentile at 180 days after boosting vaccination; and (2) high-stable individuals, with antibody values greater than 0.6 percentile of the range in the same period (median 9525 [185-80,000] AU/mL). Differential gene expression, expressed single nucleotide variants and insertions/deletions, differential splicing events, and allelic imbalance were explored to broaden our understanding of the immune response sustenance. Our analysis revealed a differential expression of genes with immunological functions in individuals with low antibody titers, compared to those with higher antibody titers, underscoring the fundamental importance of the innate immune response for boosting immunity. Our findings also provide new insights into the determinants of the immune response variability to the SARS-CoV-2 mRNA vaccine booster, highlighting the significance of differential splicing regulatory mechanisms, mainly concerning HLA alleles, in delineating vaccine immunogenicity.

Finding microbial composition and biological processes as predictive signature to access the ongoing status of mangrove preservation

Mangroves are complex land-sea transition ecosystems whose microbiota are essential for their nutrient recycling and conservation. Brazil is the third-largest estuarine area in the world and "Baía de Todos os Santos" (BTS) is one of the largest bays of the country, with wide anthropogenic exploration. Using a metagenomic approach, we investigated composition and functional adaptability as signatures of the microbiome of pristine and anthropized areas of BTS, including those under petroleum refinery influence. The taxonomic analysis showed dominance of sulfate-reducing Desulfobacteraceae, Rhodobacteraceae, and Flavobacteriaceae. Taxa were significantly diverse between pristine and disturbed areas. Disturbed mangroves showed a notary increase in abundance of halophilic, sulfur-related, and hydrocarbon-degrading genera and a decrease in diatoms compared to pristine area. The metabolic profile of BTS mangroves was correlated with the differentially abundant microbiota. Two ecological scenarios were observed: one marked by functions of central metabolism associated with biomass degradation and another by mechanisms of microbial adaptability to pollution conditions and environmental degradation. Part of the microbiome was distinct and not abundant in Brazilian estuarine soils. The microbiome signature observed in each BTS mangrove reflects how human actions impact the diversity of these ecosystems and also emphasize their role in attempting to restore disturbed mangroves. The microbiome may act as a potential biological indicator of the preservation status of these soils, despite the limitation of soil property conditions. Additionally, our data pointed to metagenomics as an additional tool for environmental assessment and reinforced the need for protective measures for the mangroves under study.

Genetic screening in a Brazilian cohort with inborn errors of immunity

BACKGROUND

Inherited genetic defects in immune system-related genes can result in Inborn Errors of Immunity (IEI), also known as Primary Immunodeficiencies (PID). Diagnosis of IEI disorders is challenging due to overlapping clinical manifestations. Accurate identification of disease-causing germline variants is crucial for appropriate treatment, prognosis, and genetic counseling. However, genetic sequencing is challenging in low-income countries like Brazil. This study aimed to perform genetic screening on patients treated within Brazil's public Unified Health System to identify candidate genetic variants associated with the patient's phenotype.

METHODS

Thirteen singleton unrelated patients from three hospitals in Rio de Janeiro were enrolled in this study. Genomic DNA was extracted from the peripheral blood lymphocytes of each patient, and whole exome sequencing (WES) analyses were conducted using Illumina NextSeq. Germline genetic variants in IEI-related genes were prioritized using a computational framework considering their molecular consequence in coding regions; minor allele frequency ≤ 0.01; pathogenicity classification based on American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines gathered from the VarSome clinical database; and IEI-related phenotype using the Franklin tool. The genes classification into IEI categories follows internationally recognized guidelines informed by the International Union of Immunological Societies Expert Committee. Additional methods for confirmation of the variant included Sanger sequencing, phasing analysis, and splice site prediction.

RESULTS

A total of 16 disease-causing variants in nine genes, encompassing six different IEI categories, were identified. X-Linked Agammaglobulinemia, caused by BTK variations, emerged as the most prevalent IEI disorder in the cohort. However, pathogenic and likely pathogenic variants were also reported in other known IEI-related genes, namely CD40LG, CARD11, WAS, CYBB, C6, and LRBA. Interestingly, two patients with suspected IEI exhibited pathogenic variants in non-IEI-related genes, ABCA12 and SLC25A13, potentially explaining their phenotypes.

CONCLUSIONS

Genetic screening through WES enabled the detection of potentially harmful variants associated with IEI disorders. These findings contribute to a better understanding of patients' clinical manifestations by elucidating the genetic basis underlying their phenotypes.

Assessing whole-exome sequencing data from undiagnosed Brazilian patients to improve the diagnostic yield of inborn errors of immunity

OBJECTIVES

Inborn error of immunity (IEI) comprises a broad group of inherited immunological disorders that usually display an overlap in many clinical manifestations challenging their diagnosis. The identification of disease-causing variants from whole-exome sequencing (WES) data comprises the gold-standard approach to ascertain IEI diagnosis. The efforts to increase the availability of clinically relevant genomic data for these disorders constitute an important improvement in the study of rare genetic disorders. This work aims to make available WES data of Brazilian patients' suspicion of IEI without a genetic diagnosis. We foresee a broad use of this dataset by the scientific community in order to provide a more accurate diagnosis of IEI disorders.

DATA DESCRIPTION

Twenty singleton unrelated patients treated at four different hospitals in the state of Rio de Janeiro, Brazil were enrolled in our study. Half of the patients were male with mean ages of 9 ± 3, while females were 12 ± 10 years old. The WES was performed in the Illumina NextSeq platform with at least 90% of sequenced bases with a minimum of 30 reads depth. Each sample had an average of 20,274 variants, comprising 116 classified as rare pathogenic or likely pathogenic according to American College of Medical Genetics and Genomics and the Association (ACMG) guidelines. The genotype-phenotype association was impaired by the lack of detailed clinical and laboratory information, besides the unavailability of molecular and functional studies which, comprise the limitations of this study. Overall, the access to clinical exome sequencing data is limited, challenging exploratory analyses and the understanding of genetic mechanisms underlying disorders. Therefore, by making these data available, we aim to increase the number of WES data from Brazilian samples despite contributing to the study of monogenic IEI-disorders.

Host genetic susceptibility underlying SARS-CoV-2-associated Multisystem Inflammatory Syndrome in Brazilian Children

BACKGROUND

Multisystem Inflammatory Syndrome in Children (MIS-C) is a life-threatening complication of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, which manifests as a hyper inflammatory process with multiorgan involvement in predominantly healthy children in the weeks following mild or asymptomatic coronavirus disease 2019 (COVID-19). However, host monogenic predisposing factors to MIS-C remain elusive.

METHODS

Herein, we used whole exome sequencing (WES) on 16 MIS-C Brazilian patients to identify single nucleotide/InDels variants as predisposition factors associated with MIS-C.

RESULTS

We identified ten very rare variants in eight genes (FREM1, MPO, POLG, C6, C9, ABCA4, ABCC6, and BSCL2) as the most promising candidates to be related to a higher risk of MIS-C development. These variants may propitiate a less effective immune response to infection or trigger the inflammatory response or yet a delayed hyperimmune response to SARS-CoV-2. Protein-Protein Interactions (PPIs) among the products of the mutated genes revealed an integrated network, enriched for immune and inflammatory response mechanisms with some of the direct partners representing gene products previously associated with MIS-C and Kawasaki disease (KD). In addition, the PPIs direct partners are also enriched for COVID-19-related gene sets. HLA alleles prediction from WES data allowed the identification of at least one risk allele in 100% of the MIS-C patients.

CONCLUSIONS

This study is the first to explore host MIS-C-associated variants in a Latin American admixed population. Besides expanding the spectrum of MIS-C-associated variants, our findings highlight the relevance of using WES for characterising the genetic interindividual variability associated with COVID-19 complications and ratify the presence of overlapping/convergent mechanisms among MIS-C, KD and COVID-19, crucial for future therapeutic management.

Buffy Coat Transcriptomic Analysis Reveals Alterations in Host Cell Protein Synthesis and Cell Cycle in Severe COVID-19 Patients

Transcriptome studies have reported the dysregulation of cell cycle-related genes and the global inhibition of host mRNA translation in COVID-19 cases. However, the key genes and cellular mechanisms that are most affected by the severe outcome of this disease remain unclear. For this work, the RNA-seq approach was used to study the differential expression in buffy coat cells of two groups of people infected with SARS-CoV-2: (a) Mild, with mild symptoms; and (b) SARS (Severe Acute Respiratory Syndrome), who were admitted to the intensive care unit with the severe COVID-19 outcome. Transcriptomic analysis revealed 1009 up-regulated and 501 down-regulated genes in the SARS group, with 10% of both being composed of long non-coding RNA. Ribosome and cell cycle pathways were enriched among down-regulated genes. The most connected proteins among the differentially expressed genes involved transport dysregulation, proteasome degradation, interferon response, cytokinesis failure, and host translation inhibition. Furthermore, interactome analysis showed Fibrillarin to be one of the key genes affected by SARS-CoV-2. This protein interacts directly with the N protein and long non-coding RNAs affecting transcription, translation, and ribosomal processes. This work reveals a group of dysregulated processes, including translation and cell cycle, as key pathways altered in severe COVID-19 outcomes.

Exploring the Bacteriome and Resistome of Humans and Food-Producing Animals in Brazil

The epidemiology of antimicrobial resistance (AMR) is complex, with multiple interfaces (human-animal-environment). In this context, One Health surveillance is essential for understanding the distribution of microorganisms and antimicrobial resistance genes (ARGs). This report describes a multicentric study undertaken to evaluate the bacterial communities and resistomes of food-producing animals (cattle, poultry, and swine) and healthy humans sampled simultaneously from five Brazilian regions. Metagenomic analysis showed that a total of 21,029 unique species were identified in 107 rectal swabs collected from distinct hosts, the highest numbers of which belonged to the domain Bacteria, mainly Ruminiclostridium spp. and Bacteroides spp., and the order Enterobacterales. We detected 405 ARGs for 12 distinct antimicrobial classes. Genes encoding antibiotic-modifying enzymes were the most frequent, followed by genes related to target alteration and efflux systems. Interestingly, carbapenemase-encoding genes such as blaAIM-1, blaCAM-1, blaGIM-2, and blaHMB-1 were identified in distinct hosts. Our results revealed that, in general, the bacterial communities from humans were present in isolated clusters, except for the Northeastern region, where an overlap of the bacterial species from humans and food-producing animals was observed. Additionally, a large resistome was observed among all analyzed hosts, with emphasis on the presence of carbapenemase-encoding genes not previously reported in Latin America. IMPORTANCE Humans and food production animals have been reported to be important reservoirs of antimicrobial resistance (AMR) genes (ARGs). The frequency of these multidrug-resistant (MDR) bacteria tends to be higher in low- and middle-income countries (LMICs), due mainly to a lack of public health policies. Although studies on AMR in humans or animals have been carried out in Brazil, this is the first multicenter study that simultaneously collected rectal swabs from humans and food-producing animals for metagenomics. Our results indicate high microbial diversity among all analyzed hosts, and several ARGs for different antimicrobial classes were also found. As far as we know, we have detected for the first time ARGs encoding carbapenemases, such as blaAIM-1, blaCAM-1, blaGIM-2, and blaHMB-1, in Latin America. Thus, our results support the importance of metagenomics as a tool to track the colonization of food-producing animals and humans by antimicrobial-resistant bacteria. In addition, a network surveillance system called GUARANI, created for this study, is ready to be expanded and to collect additional data.

Seroprevalence, Prevalence, and Genomic Surveillance: Monitoring the Initial Phases of the SARS-CoV-2 Pandemic in Betim, Brazil

The COVID-19 pandemic has created an unprecedented need for epidemiological monitoring using diverse strategies. We conducted a project combining prevalence, seroprevalence, and genomic surveillance approaches to describe the initial pandemic stages in Betim City, Brazil. We collected 3239 subjects in a population-based age-, sex- and neighborhood-stratified, household, prospective; cross-sectional study divided into three surveys 21 days apart sampling the same geographical area. In the first survey, overall prevalence (participants positive in serological or molecular tests) reached 0.46% (90% CI 0.12–0.80%), followed by 2.69% (90% CI 1.88–3.49%) in the second survey and 6.67% (90% CI 5.42–7.92%) in the third. The underreporting reached 11, 19.6, and 20.4 times in each survey. We observed increased odds to test positive in females compared to males (OR 1.88 95% CI 1.25–2.82), while the single best predictor for positivity was ageusia/anosmia (OR 8.12, 95% CI 4.72–13.98). Thirty-five SARS-CoV-2 genomes were sequenced, of which 18 were classified as lineage B.1.1.28, while 17 were B.1.1.33. Multiple independent viral introductions were observed. Integration of multiple epidemiological strategies was able to adequately describe COVID-19 dispersion in the city. Presented results have helped local government authorities to guide pandemic management.

Exome-Wide Search for Genes Associated With Central Nervous System Inflammatory Demyelinating Diseases Following CHIKV Infection: The Tip of the Iceberg

Chikungunya virus (CHIKV) is a re-emergent arbovirus that causes a disease characterized primarily by fever, rash and severe persistent polyarthralgia, although <1% of cases develop severe neurological manifestations such as inflammatory demyelinating diseases (IDD) of the central nervous system (CNS) like acute disseminated encephalomyelitis (ADEM) and extensive transverse myelitis. Genetic factors associated with host response and disease severity are still poorly understood. In this study, we performed whole-exome sequencing (WES) to identify HLA alleles, genes and cellular pathways associated with CNS IDD clinical phenotype outcomes following CHIKV infection. The cohort includes 345 patients of which 160 were confirmed for CHIKV. Six cases presented neurological manifestation mimetizing CNS IDD. WES data analysis was performed for 12 patients, including the CNS IDD cases and 6 CHIKV patients without any neurological manifestation. We identified 29 candidate genes harboring rare, pathogenic, or probably pathogenic variants in all exomes analyzed. HLA alleles were also determined and patients who developed CNS IDD shared a common signature with diseases such as Multiple sclerosis (MS) and Neuromyelitis Optica Spectrum Disorders (NMOSD). When these genes were included in Gene Ontology analyses, pathways associated with CNS IDD syndromes were retrieved, suggesting that CHIKV-induced CNS outcomesmay share a genetic background with other neurological disorders. To our knowledge, this study was the first genome-wide investigation of genetic risk factors for CNS phenotypes in CHIKV infection. Our data suggest that HLA-DRB1 alleles associated with demyelinating diseases may also confer risk of CNS IDD outcomes in patients with CHIKV infection.

Common Dysregulation of Innate Immunity Pathways in Human Primary Astrocytes Infected With Chikungunya, Mayaro, Oropouche, and Zika Viruses

Arboviruses pose a major threat throughout the world and represent a great burden in tropical countries of South America. Although generally associated with moderate febrile illness, in more severe cases they can lead to neurological outcomes, such as encephalitis, Guillain-Barré syndrome, and Congenital Syndromes. In this context astrocytes play a central role in production of inflammatory cytokines, regulation of extracellular matrix, and control of glutamate driven neurotoxicity in the central nervous system. Here, we presented a comprehensive genome-wide transcriptome analysis of human primary astrocytes infected with Chikungunya, Mayaro, Oropouche, or Zika viruses. Analyses of differentially expressed genes (DEGs), pathway enrichment, and interactomes have shown that Alphaviruses up-regulated genes related to elastic fiber formation and N-glycosylation of glycoproteins, with down-regulation of cell cycle and DNA stability and chromosome maintenance genes. In contrast, Oropouche virus up-regulated cell cycle and DNA maintenance and condensation pathways while down-regulated extracellular matrix, collagen metabolism, glutamate and ion transporters pathways. Zika virus infection only up-regulated eukaryotic translation machinery while down-regulated interferon pathways. Reactome and integration analysis revealed a common signature in down-regulation of innate immune response, antiviral response, and inflammatory cytokines associated to interferon pathway for all arboviruses tested. Validation of interferon stimulated genes by reverse transcriptase quantitative polymerase chain reaction (RT-qPCR) corroborated our transcriptome findings. Altogether, our results showed a co-evolution in the mechanisms involved in the escape of arboviruses to antiviral immune response mediated by the interferon (IFN) pathway.

Transcriptomic Analysis of the Early Strobilar Development of Echinococcus granulosus

Echinococcus granulosus has a complex life cycle involving two mammalian hosts. The transition from one host to another is accompanied by changes in gene expression, and the transcriptional events that underlie this transition have not yet been fully characterized. In this study, RNA-seq was used to compare the transcription profiles of samples from E. granulosus protoscoleces induced in vitro to strobilar development at three time points. We identified 818 differentially expressed genes, which were divided into eight expression clusters formed over the entire 24 h period. An enrichment of gene transcripts with molecular functions of signal transduction, enzymes, and protein modifications was observed upon induction and developmental progression. This transcriptomic study provides insights for understanding the complex life cycle of E. granulosus and contributes for searching for the key genes correlating with the strobilar development, which can be used to identify potential candidates for the development of anthelmintic drugs.

Genome-wide DNA methylation analysis of Metarhizium anisopliae during tick mimicked infection condition

BACKGROUND

The Metarhizium genus harbors important entomopathogenic fungi. These species have been widely explored as biological control agents, and strategies to improve the fungal virulence are under investigation. Thus, the interaction between Metarhizium species and susceptible hosts have been explored employing different methods in order to characterize putative virulence determinants. However, the impact of epigenetic modulation on the infection cycle of Metarhizium is still an open topic. Among the different epigenetic modifications, DNA methylation of cytosine bases is an important mechanism to control gene expression in several organisms. To better understand if DNA methylation can govern Metarhizium-host interactions, the genome-wide DNA methylation profile of Metarhizium anisopliae was explored in two conditions: tick mimicked infection and a saprophytic-like control.

RESULTS

Using a genome wide DNA methylation profile based on bisulfite sequencing (BS-Seq), approximately 0.60% of the total cytosines were methylated in saprophytic-like condition, which was lower than the DNA methylation level (0.89%) in tick mimicked infection condition. A total of 670 mRNA genes were found to be putatively methylated, with 390 mRNA genes uniquely methylated in the tick mimicked infection condition. GO terms linked to response to stimuli, cell wall morphogenesis, cytoskeleton morphogenesis and secondary metabolism biosynthesis were over-represented in the tick mimicked infection condition, suggesting that energy metabolism is directed towards the regulation of genes associated with infection. However, recognized virulence determinants known to be expressed at distinct infection steps, such as the destruxin backbone gene and the collagen-like protein gene Mcl1, were found methylated, suggesting that a dynamic pattern of methylation could be found during the infectious process. These results were further endorsed employing RT-qPCR from cultures treated or not with the DNA methyltransferase inhibitor 5-Azacytidine.

CONCLUSIONS

The set of genes here analyzed focused on secondary metabolites associated genes, known to be involved in several processes, including virulence. The BS-Seq pipeline and RT-qPCR analysis employing 5-Azacytidine led to identification of methylated virulence genes in M. anisopliae. The results provided evidences that DNA methylation in M. anisopliae comprises another layer of gene expression regulation, suggesting a main role of DNA methylation regulating putative virulence determinants during M. anisopliae infection cycle.

A comparative in silico linear B-cell epitope prediction and characterization for South American and African Trypanosoma vivax strains

Trypanosoma vivax is a parasite widespread across Africa and South America. Immunological methods using recombinant antigens have been developed aiming at specific and sensitive detection of infections caused by T. vivax. Here, we sequenced for the first time the transcriptome of a virulent T. vivax strain (Lins), isolated from an outbreak of severe disease in South America (Brazil) and performed a computational integrated analysis of genome, transcriptome and in silico predictions to identify and characterize putative linear B-cell epitopes from African and South American T. vivax. A total of 2278, 3936 and 4062 linear B-cell epitopes were respectively characterized for the transcriptomes of T. vivax LIEM-176 (Venezuela), T. vivax IL1392 (Nigeria) and T. vivax Lins (Brazil) and 4684 for the genome of T. vivax Y486 (Nigeria). The results presented are a valuable theoretical source that may pave the way for highly sensitive and specific diagnostic tools.

Metagenomic analysis of sediments under seaports influence in the Equatorial Atlantic Ocean

Maritime ports are anthropogenic interventions capable of causing serious alterations in coastal ecosystems. In this study, we examined the benthic microbial diversity and community structure under the influence of two maritime ports, Mucuripe (MUC) and Pecém (PEC), at Equatorial Atlantic Ocean in Northeast Brazil. Those seaports differ in architecture, time of functioning, cargo handling and contamination. The microbiomes from MUC and PEC were also compared in silico to 11 other globally distributed marine microbiomes. The comparative analysis of operational taxonomic units (OTUs) retrieved by PCR-DGGE showed that MUC presents greater richness and β diversity of Bacteria and Archaea than PEC. In line with these results, metagenomic analysis showed that MUC and PEC benthic microbial communities share the main common bacterial phyla found in coastal environments, although can be distinguish by greater abundance of Cyanobacteria in MUC and Deltaproteobacteria in PEC. Both ports differed in Archaea composition, being PEC port sediments dominated by Thaumarchaeota. The microbiomes showed little divergence in their potential metabolic pathways, although shifts on the microbial taxonomic signatures involved in nitrogen and sulphur metabolic pathways were observed. The comparative analysis of different benthic marine metagenomes from Brazil, Australia and Mexico grouped them by the geographic location rather than by the type of ecosystem, although at phylum level seaport sediments share a core microbiome constituted by Proteobacteria, Cyanobacteria, Actinobacteria, Tenericuteres, Firmicutes, Bacteriodetes and Euryarchaeota. Our results suggest that multiple physical and chemical factors acting on sediments as a result of at least 60years of port operation play a role in shaping the benthic microbial communities at taxonomic level, but not at functional level.

High expression of FOXP3 in primary melanoma is associated with tumour progression

BACKGROUND

The antitumour immune response plays an important role in the prognosis of melanoma. High numbers of circulating regulatory T cells have been associated with rapid disease progression.

OBJECTIVES

To assess the influence of forkhead box protein (FOXP)3, CD1a and langerin expression on the prognosis of primary melanoma.

METHODS

We analysed 185 primary melanomas by immunohistochemical staining for expression of the regulatory T-cell marker FOXP3 and the dendritic cell markers langerin and CD1a, and correlated marker expression with clinical outcome.

RESULTS

Disease-free survival and overall survival were significantly longer in patients expressing low levels of FOXP3 in the primary melanoma, whereas they were associated with high expression of CD1a. The negative prognostic value of FOXP3 expression was independent of the Breslow tumour thickness. Langerin expression did not correlate with the clinical outcome.

CONCLUSIONS

High expression of FOXP3 in the primary melanoma may be used as an additional independent prognostic marker for early tumour progression in patients with melanoma.

Unravelling the transcriptome profile of the Swine respiratory tract mycoplasmas

The swine respiratory ciliary epithelium is mainly colonized by Mycoplasma hyopneumoniae, Mycoplasma flocculare and Mycoplasma hyorhinis. While colonization by M. flocculare is virtually asymptomatic, M. hyopneumoniae and M. hyorhinis infections may cause respiratory disease. Information regarding transcript structure and gene abundance provides valuable insight into gene function and regulation, which has not yet been analyzed on a genome-wide scale in these Mycoplasma species. In this study, we report the construction of transcriptome maps for M. hyopneumoniae, M. flocculare and M. hyorhinis, which represent data for conducting comparative studies on the transcriptional repertory. For each species, three cDNA libraries were generated, yielding averages of 415,265, 695,313 and 93,578 reads for M. hyopneumoniae, M. flocculare and M. hyorhinis, respectively, with an average read length of 274 bp. The reads mapping showed that 92%, 98% and 96% of the predicted genes were transcribed in the M. hyopneumoniae, M. flocculare and M. hyorhinis genomes, respectively. Moreover, we showed that the majority of the genes are co-expressed, confirming the previously predicted transcription units. Finally, our data defined the RNA populations in detail, with the map transcript boundaries and transcription unit structures on a genome-wide scale.

Comparative genome analysis of entomopathogenic fungi reveals a complex set of secreted proteins

Background

Metarhizium anisopliae is an entomopathogenic fungus used in the biological control of some agricultural insect pests, and efforts are underway to use this fungus in the control of insect-borne human diseases. A large repertoire of proteins must be secreted by M. anisopliae to cope with the various available nutrients as this fungus switches through different lifestyles, i.e., from a saprophytic, to an infectious, to a plant endophytic stage. To further evaluate the predicted secretome of M. anisopliae, we employed genomic and transcriptomic analyses, coupled with phylogenomic analysis, focusing on the identification and characterization of secreted proteins.

Results

We determined the M. anisopliae E6 genome sequence and compared this sequence to other entomopathogenic fungi genomes. A robust pipeline was generated to evaluate the predicted secretomes of M. anisopliae and 15 other filamentous fungi, leading to the identification of a core of secreted proteins. Transcriptomic analysis using the tick Rhipicephalus microplus cuticle as an infection model during two periods of infection (48 and 144 h) allowed the identification of several differentially expressed genes. This analysis concluded that a large proportion of the predicted secretome coding genes contained altered transcript levels in the conditions analyzed in this study. In addition, some specific secreted proteins from Metarhizium have an evolutionary history similar to orthologs found in Beauveria/Cordyceps. This similarity suggests that a set of secreted proteins has evolved to participate in entomopathogenicity.

Conclusions

The data presented represents an important step to the characterization of the role of secreted proteins in the virulence and pathogenicity of M. anisopliae.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-822) contains supplementary material, which is available to authorized users.

De novo transcriptome analysis of Hevea brasiliensis tissues by RNA-seq and screening for molecular markers

BACKGROUND

The rubber tree, Hevea brasiliensis, is a species native to the Brazilian Amazon region and it supplies almost all the world's natural rubber, a strategic raw material for a variety of products. One of the major challenges for developing rubber tree plantations is adapting the plant to biotic and abiotic stress. Transcriptome analysis is one of the main approaches for identifying the complete set of active genes in a cell or tissue for a specific developmental stage or physiological condition.

RESULTS

Here, we report on the sequencing, assembling, annotation and screening for molecular markers from a pool of H. brasiliensis tissues. A total of 17,166 contigs were successfully annotated. Then, 2,191 Single Nucleotide Variation (SNV) and 1.397 Simple Sequence Repeat (SSR) loci were discriminated from the sequences. From 306 putative, mainly non-synonymous SNVs located in CDS sequences, 191 were checked for their ability to characterize 23 Hevea genotypes by an allele-specific amplification technology. For 172 (90%), the nucleotide variation at the predicted genomic location was confirmed, thus validating the different steps from sequencing to the in silico detection of the SNVs.

CONCLUSIONS

This is the first study of the H. brasiliensis transcriptome, covering a wide range of tissues and organs, leading to the production of the first developed SNP markers. This process could be amplified to a larger set of in silico detected SNVs in expressed genes in order to increase the marker density in available and future genetic maps. The results obtained in this study will contribute to the H. brasiliensis genetic breeding program focused on improving of disease resistance and latex yield.

Prerace venous blood acid-base values in standardbred horses

The administration of sodium bicarbonate (NaHCO3) or 'milk shakes' to Standardbred horses before racing is widespread. This study analysed nonrace day (CTL) and prerace venous acid-base values from Standardbred horses racing in Pennsylvania (PA) and New Jersey (NJ). Mean +/- s.d. CTL bicarbonate (HCO3) and base excess (BE) values, for a group of horses stabled during the 1993 racing season at Pocono Downs, Pennsylvania, were 28.6 +/- 1.9 and 2.6 +/- 1.7 mmol/l, respectively. In the same population of horses, mean +/- s.d. values for prerace HCO3- and BE values were 33.1 +/- 2.8 and 7.0 +/- 2.3 mmol/l, respectively, for horses administered frusemide (F) 4 h before race time and 31.5 +/- 2.4 and 5.5 +/- 2.0 mmol/l for the horses not administered frusemide (NF). There were differences (P < 0.05) in pHv, PvCO2, HCO3- and BE values between the CTL and prerace samples. The venous acid-base values of the CTL horses were normally distributed. The prerace acid-base values measured during 1993 were not normally distributed, indicating changes due either to the administration of alkalising substances or other manipulations of the horses on race day. Changes in all acid-base values were observed during the subsequent 1994 racing season with further changes observed when horses were placed in a secured stable 8 h before the race. The criteria used in some racing jurisdictions for disqualifying a horse were the elevations in HCO3-, Na+ and pH. The correlation coefficients (r2) for the least squares regression for HCO3- vs. pH and HCO3- vs. Na+ in prerace venous blood samples were 0.31 and 0.21, respectively, indicating a poor relationship between the 3 acid-base values. To discourage the administration of NaHCO3 to horses before a race, the use of a BE value of > or = 10 mmol/l for NF and > or = 12 mmol/l for F was adopted as a single index for the disqualification of horses from a race. The results of this study indicate that the use of a single index and, in this case BE, has curtailed the prerace administration of NaHCO3 to horses.