EasySSR: a user-friendly web application with full command-line features for large-scale batch microsatellite mining and samples comparison

Microsatellites, also known as SSRs or STRs, are polymorphic DNA regions with tandem repetitions of a nucleotide motif of size 1-6 base pairs with a broad range of applications in many fields, such as comparative genomics, molecular biology, and forensics. However, the majority of researchers do not have computational training and struggle while running command-line tools or very limited web tools for their SSR research, spending a considerable amount of time learning how to execute the software and conducting the post-processing data tabulation in other tools or manually-time that could be used directly in data analysis. We present EasySSR, a user-friendly web tool with command-line full functionality, designed for practical use in batch identifying and comparing SSRs in sequences, draft, or complete genomes, not requiring previous bioinformatic skills to run. EasySSR requires only a FASTA and an optional GENBANK file of one or more genomes to identify and compare STRs. The tool can automatically analyze and compare SSRs in whole genomes, convert GenBank to PTT files, identify perfect and imperfect SSRs and coding and non-coding regions, compare their frequencies, abundancy, motifs, flanking sequences, and iterations, producing many outputs ready for download such as PTT files, interactive charts, and Excel tables, giving the user the data ready for further analysis in minutes. EasySSR was implemented as a web application, which can be executed from any browser and is available for free at https://computationalbiology.ufpa.br/easyssr/. Tutorials, usage notes, and download links to the source code can be found at https://github.com/engbiopct/EasySSR.

Immunoexpression of stem cell markers SOX-2, NANOG AND OCT4 in ameloblastoma

Background

Ameloblastoma (AME) is characterized by a locally invasive growth pattern. In an attempt to justify the aggressiveness of neoplasms, the investigation of the role of stem cells has gained prominence. The SOX-2, NANOG and OCT4 proteins are important stem cell biomarkers.

Methodology

To verify the expression of these proteins in tissue samples of AME, dentigerous cyst (DC) and dental follicle (DF), immunohistochemistry was performed and indirect immunofluorescence were performed on the human AME (AME-hTERT) cell line.

Results

Revealed expression of SOX-2, NANOG and OCT4 in the tissue samples and AME-hTERT lineage. Greater immunostaining of the studied proteins was observed in AME compared to DC and DF (p < 0.001).

Conclusions

The presence of biomarkers indicates a probable role of stem cells in the genesis and progression of AME.

In silico approach to identify microsatellite candidate biomarkers to differentiate the biovar of Corynebacterium pseudotuberculosis genomes

Corynebacterium pseudotuberculosis is the causative bacterial agent of the zoonotic disease known as caseous lymphadenitis, and it presents several mechanisms of response to host defenses, including the presence of virulence factors (VFs). The genomes of these bacteria have several polymorphic markers known as microsatellites, or simple sequence repeats (SSRs), that can be used to characterize the genome, to study possible polymorphisms existing among strains, and to verify the effects of such polymorphic markers in coding regions and regions associated with VFs. In this study, several SSRs were identified within coding regions throughout the 54 genomes of this species, revealing possible polymorphisms associated with coding regions that could be used as strain-specific or serotype-specific identifiers of C. pseudotuberculosis. The similarities associated with SSRs amongst the different serum variants of C. pseudotuberculosis, biovars equi and ovis, were also evaluated, and it was possible to identify SSRs located in coding regions responsible for a VF enrolled in pathogenesis known to mediate bacterial adherence (SpaH-type pili virulence factor). Phylogenetic analyses revealed that strains sharing SSR patterns, including the possible polymorphisms identified in the same position of gene-coding regions, were displayed by strains with a common ancestor, corroborating with the Genome Tree Report of the NCBI. Statistical analysis showed that the microsatellite groups belonging to equi and ovis biovars have a significance of 0.006 (p-value) in similarity, thus indicating them as good biomarker candidates for C. pseudotuberculosis.

Gene expression in cell lines from human ameloblastoma immortalized using hTERT and HPV16-E6/E7

OBJECTIVE

The aim of this study was to evaluate and compare alterations in gene expression using two distinct immortalization methods (hTERT and HPV16-E6/E7) in ameloblastoma cell lines.

MATERIALS AND METHODS

A primary cell culture derived from human ameloblastoma (AME-1) was established and immortalized by two different methods using a transfection processes to hTERT and HPV-E6/E7. The RNA-seq was used to verify which immortalization method had less influence on gene expression. It was performed in four steps: extraction and collection of mRNA, PCR amplification, comparison with the human reference genome, and analysis of differential expression. The genes with differentiated expression were identified and mapped.

RESULTS

RNA-seq revealed genetic alterations in ameloblastoma cell lines after the immortalization process, including increased expression of tumor genes like MYC, E2F1, BRAF, HRAS, and HTERT, and a decrease in tumor suppressor genes like P53, P21, and Rb.

CONCLUSIONS

It is possible to affirm that cell immortalization is not an inert method regarding gene regulation mechanisms and the hTERT method (AME-TERT) presented fewer changes in gene expression levels.

Characterization of pharmacogenetic markers related to Acute Lymphoblastic Leukemia toxicity in Amazonian native Americans population

Acute Lymphoblastic Leukemia (ALL) is the most common cancer in children. Differences are found among ethnic groups in the results of the treatment of pediatric ALL. In general, children with a high level of native American ancestry tend to respond less positively to ALL treatments, which may be related to specific genomic variants found in native American groups. Despite the evidence, few data are available on the distribution of the pharmacogenomic variants relevant to the treatment of ALL in traditional Amerindian populations, such the those of the Amazon region. Given this, the present study investigated 27 molecular markers related to the treatment of ALL in Amerindians from Brazilian Amazonia and compared the frequencies with those recorded previously on five continents, that are available in the 1,000 Genomes database. The variation in the genotype frequencies among populations was evaluated using Fisher's exact test. The False Discovery Rate method was used to correct the results of the multiple analyses. Significant differences were found in the frequencies of the majority of markers between the Amerindian populations and those of other regions around the world. These findings highlight the unique genetic profile of the indigenous population of Brazilian Amazonia, which may reflect a distinct therapeutic profile for the treatment of ALL in these populations.

Bacterial diversity in deep-sea sediments under influence of asphalt seep at the São Paulo Plateau

Here we investigated the diversity of bacterial communities from deep-sea surface sediments under influence of asphalt seeps at the Sao Paulo Plateau using next-generation sequencing method. Sampling was performed at North São Paulo Plateau using the human occupied vehicle Shinkai 6500 and her support vessel Yokosuka. The microbial diversity was studied at two surficial sediment layers (0-1 and 1-4 cm) of five samples collected in cores in water depths ranging from 2456 to 2728 m. Bacterial communities were studied through sequencing of 16S rRNA gene on the Ion Torrent platform and clustered in operational taxonomic units. We observed high diversity of bacterial sediment communities as previously described by other studies. When we considered community composition, the most abundant classes were Alphaproteobacteria (27.7%), Acidimicrobiia (20%), Gammaproteobacteria (11.3%) and Deltaproteobacteria (6.6%). Most abundant OTUs at family level were from two uncultured bacteria from Actinomarinales (5.95%) and Kiloniellaceae (3.17%). The unexpected high abundance of Alphaproteobacteria and Acidimicrobiia in our deep-sea microbial communities may be related to the presence of asphalt seep at North São Paulo Plateau, since these bacterial classes contain bacteria that possess the capability of metabolizing hydrocarbon compounds.

PhageWeb - Web Interface for Rapid Identification and Characterization of Prophages in Bacterial Genomes

This study developed a computational tool with a graphical interface and a web-service that allows the identification of phage regions through homology search and gene clustering. It uses G+C content variation evaluation and tRNA prediction sites as evidence to reinforce the presence of prophages in indeterminate regions. Also, it performs the functional characterization of the prophages regions through data integration of biological databases. The performance of PhageWeb was compared to other available tools (PHASTER, Prophinder, and PhiSpy) using Sensitivity (Sn) and Positive Predictive Value (PPV) tests. As a reference for the tests, more than 80 manually annotated genomes were used. In the PhageWeb analysis, the Sn index was 86.1% and the PPV was approximately 87%, while the second best tool presented Sn and PPV values of 83.3 and 86.5%, respectively. These numbers allowed us to observe a greater precision in the regions identified by PhageWeb while compared to other prediction tools submitted to the same tests. Additionally, PhageWeb was much faster than the other computational alternatives, decreasing the processing time to approximately one-ninth of the time required by the second best software. PhageWeb is freely available at http://computationalbiology.ufpa.br/phageweb.

Genome of Leptospira borgpetersenii strain 4E, a highly virulent isolate obtained from Mus musculus in southern Brazil

A previous study by our group reported the isolation and characterisation of Leptospira borgpetersenii serogroup Ballum strain 4E. This strain is of particular interest because it is highly virulent in the hamster model. In this study, we performed whole-genome shotgun genome sequencing of the strain using the SOLiD sequencing platform. By assembling and analysing the new genome, we were able to identify novel features that have been previously overlooked in genome annotations of other strains belonging to the same species.

Draft genome sequence of Psychrobacter sp. ENNN9_III, a strain isolated from water in a polluted temperate estuarine system ( Ria de Aveiro , Portugal)

The genus Psychrobacter includes Gram-negative coccobacilli that are non-pigmented, oxidase-positive, non-motile, psychrophilic or psychrotolerant, and halotolerant. Psychrobacter strain ENNN9_III was isolated from water in a polluted temperate estuarine system, contaminated with hydrocarbons and heavy metals. The genome has a G + C content of 42.7%, 2618 open reading frames (ORFs), three copies of the rRNAs operon, and 29 tRNA genes.

Twenty-five sequences related to the degradation of aromatic compounds were predicted, as well as numerous genes related to resistance to metals or metal(loid)s. The genome sequence of Psychrobacter strain ENNN9_III provides the groundwork for further elucidation of the mechanisms of metal resistance and aromatic compounds degradation. Future studies are needed to confirm the usefulness of this strain for bioremediation proposes.

Brazilian Microbiome Project: revealing the unexplored microbial diversity--challenges and prospects

The Brazilian Microbiome Project (BMP) aims to assemble a Brazilian Metagenomic Consortium/Database. At present, many metagenomic projects underway in Brazil are widely known. Our goal in this initiative is to co-ordinate and standardize these together with new projects to come. It is estimated that Brazil hosts approximately 20 % of the entire world's macroorganism biological diversity. It is 1 of the 17 countries that share nearly 70 % of the world's catalogued animal and plant species, and is recognized as one of the most megadiverse countries. At the end of 2012, Brazil has joined GBIF (Global Biodiversity Information Facility), as associated member, to improve the access to the Brazilian biodiversity data in a free and open way. This was an important step toward increasing international collaboration and clearly shows the commitment of the Brazilian government in directing national policies toward sustainable development. Despite its importance, the Brazilian microbial diversity is still considered to be largely unknown, and it is clear that to maintain ecosystem dynamics and to sustainably manage land use, it is crucial to understand the biological and functional diversity of the system. This is the first attempt to collect and collate information about Brazilian microbial genetic and functional diversity in a systematic and holistic manner. The success of the BMP depends on a massive collaborative effort of both the Brazilian and international scientific communities, and therefore, we invite all colleagues to participate in this project.

The genome of Anopheles darlingi, the main neotropical malaria vector

Anopheles darlingi is the principal neotropical malaria vector, responsible for more than a million cases of malaria per year on the American continent. Anopheles darlingi diverged from the African and Asian malaria vectors ∼100 million years ago (mya) and successfully adapted to the New World environment. Here we present an annotated reference A. darlingi genome, sequenced from a wild population of males and females collected in the Brazilian Amazon. A total of 10 481 predicted protein-coding genes were annotated, 72% of which have their closest counterpart in Anopheles gambiae and 21% have highest similarity with other mosquito species. In spite of a long period of divergent evolution, conserved gene synteny was observed between A. darlingi and A. gambiae. More than 10 million single nucleotide polymorphisms and short indels with potential use as genetic markers were identified. Transposable elements correspond to 2.3% of the A. darlingi genome. Genes associated with hematophagy, immunity and insecticide resistance, directly involved in vector–human and vector–parasite interactions, were identified and discussed. This study represents the first effort to sequence the genome of a neotropical malaria vector, and opens a new window through which we can contemplate the evolutionary history of anopheline mosquitoes. It also provides valuable information that may lead to novel strategies to reduce malaria transmission on the South American continent. The A. darlingi genome is accessible at www.labinfo.lncc.br/index.php/anopheles-darlingi.

Whole genome sequencing of environmental Vibrio cholerae O1 from 10 nanograms of DNA using short reads

Multiple Displacement Amplification (MDA) of DNA using φ29 (phi29) DNA polymerase amplifies DNA several billion-fold, which has proved to be potentially very useful for evaluating genome information in a culture-independent manner. Whole genome sequencing using DNA from a single prokaryotic genome copy amplified by MDA has not yet been achieved due to the formation of chimeras and skewed amplification of genomic regions during the MDA step, which then precludes genome assembly. We have hereby addressed the issue by using 10 ng of genomic Vibrio cholerae DNA extracted within an agarose plug to ensure circularity as a starting point for MDA and then sequencing the amplified yield using the SOLiD platform. We successfully managed to assemble the entire genome of V. cholerae strain LMA3984-4 (environmental O1 strain isolated in urban Amazonia) using a hybrid de novo assembly strategy. Using our method, only 178 out of 16,713 (1%) of contigs were not able to be inserted into either chromosome scaffold, and out of these 178, only 3 appeared to be chimeras. The other contigs seem to be the result of template-independent non-specific amplification during MDA, yielding spurious reads. Extraction of genomic DNA within an agarose plug in order to ensure circularity of the extracted genome might be key to minimizing amplification bias by MDA for WGS.

Characterization of 23S rRNA domain V mutations in gastric biopsy patients from the eastern Amazon

Resistance of Helicobacter pylori to clarithromycin is characterised by simple point mutations in the 23S ribosomal RNA (rRNA) gene and is responsible for the majority of cases of failure to eradicate this bacterium. In this paper, we characterised the variability of the 23S rRNA gene in biopsies of patients with gastric pathologies in the eastern Amazon (Northern Region of Brazil) using PCR and sequencing. A total of 49 sequences of H. pylori strains were analysed and of those, 75.6% presented nucleotide substitutions: A2142G (3.3%), T2182C (12.9%), G2224A (6.45%), T2215C (61.3%), A2192G (3.3%), G2204C (6.4%) and T2221C (6.4%). Of the mutations identified, four are known mutations related to cases of resistance and 16.1% are not yet described, revealing a high prevalence of mutations in the H. pylori 23S rRNA gene among the strains circulating in the in the eastern Amazon. The high prevalence in individuals with gastric pathologies in the Northern Region of Brazil demonstrates the need for characterising the profile of these strains to provide correct therapy for patients, considering that mutations in this gene are normally associated with resistance to the primary medication used in controlling H. pylori infection.

An alternative genotyping method using dye-labeled universal primer to reduce unspecific amplifications

We proposed a modification the procedure of genotyping based in labeled universal primer and tailed primer. In the standard protocol, three primers are used in the same PCR reaction, a forward primer with tail added at the 5' end of the identical sequence to labeled universal primer with dye-fluorescent and a reverse primer. Unfortunately, the choice of a labeled primer characterized by a large number of complementary sequences in target genomes (which is more probable in larger genomes) result in unspecific amplifications (false positive) can cause absence or decrease amplification of the locus of interest and also false interpretation of the analysis. However, identification of possible homologies between the primer chosen for labelling and the genome is rarely possible from the available DNA data bases. In our approach, cycling is interrupted for the addition of the labeled primer only during the final cycles, thus minimizing unspecific amplification and competition between primers, resulting in the more fidelity amplification of the target regions.